Next Article in Journal
Altered Development of Mesencephalic Dopaminergic Neurons in SIDS: New Insights into Understanding Sudden Infant Death Pathogenesis
Next Article in Special Issue
Identification of NSP3 (SH2D3C) as a Prognostic Biomarker of Tumor Progression and Immune Evasion for Lung Cancer and Evaluation of Organosulfur Compounds from Allium sativum L. as Therapeutic Candidates
Previous Article in Journal
CCR6 Deficiency Increases Infarct Size after Murine Acute Myocardial Infarction
Previous Article in Special Issue
Impact of Protein Corona on the Biological Identity of Nanomedicine: Understanding the Fate of Nanomaterials in the Biological Milieu
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Analysis of Punicalagin Interaction with PDIA1 and PDIA3 by Biochemical and Computational Approaches

1
Department of Biochemical Sciences “Alessandro Rossi Fanelli”, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
2
Rome Center for Molecular Design, Department of Drug Chemistry and Technology, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
3
Enrico ed Enrica Sovena Foundation, 00199 Rome, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Biomedicines 2021, 9(11), 1533; https://doi.org/10.3390/biomedicines9111533
Submission received: 30 September 2021 / Revised: 19 October 2021 / Accepted: 21 October 2021 / Published: 25 October 2021

Abstract

:
In a previous work, it was shown that punicalagin, an active ingredient of pomegranate, is able to bind to PDIA3 and inhibit its disulfide reductase activity. Here we provide evidence that punicalagin can also bind to PDIA1, the main expressed form of protein disulfide isomerase (PDI). In this comparative study, the affinity and the effect of punicalagin binding on each protein were evaluated, and a computational approach was used to identify putative binding sites. Punicalagin binds to either PDIA1 or PDIA3 with a similar affinity, but the inhibition efficacy on protein reductase activity is higher for PDIA3. Additionally, punicalagin differently affects the thermal denaturation profile of both proteins. Molecular docking and molecular dynamics simulations led to propose a punicalagin binding mode on PDIA1 and PDIA3, identifying the binding sites at the redox domains a’ in two different pockets, suggesting different effects of punicalagin on proteins’ structure. This study provides insights to develop punicalagin-based ligands, to set up a rational design for PDIA3 selective inhibitors, and to dissect the molecular determinant to modulate the protein activity.

Graphical Abstract

1. Introduction

Protein disulfide isomerase A3 (PDIA3) is a multifunctional protein acting in multiple cellular compartments [1]. PDIA3 catalyzes proper disulfide bonds and thiol groups formation on target proteins through its thiol-oxidoreductase and disulfide isomerase properties [2].
PDIA3 is a 505 amino acids protein folded in four globular domains, named, respectively, a, b, b’, and a’, which were experimentally observed to be organized in a U-shaped scaffold (PDB ID: 3F8U) [3]. Even though each domain is characterized by a thioredoxin-like structure, only the first and fourth domains, a and a’, were reported as catalytically active, while the central domains provide binding sites for PDIA3 interactors [4,5].
Recently, increased knowledge on the PDIA3 system is being gradually added to shed light on its role in different cellular functions and human pathologies. There are several milestones in the functions of PDIA3. In the endoplasmic reticulum (ER), PDIA3 modulates the synthesis of newly synthesized N-glycosylated proteins interacting with lectin-binding proteins such as calreticulin and calnexin, as described [6]. In ER, PDIA3 regulates the redox status of major histocompatibility complex (MHC) class I during antigen presentation [7]. Moreover, PDIA3 has been proven as a modulator of the 1,25-dihydroxycholecalciferol (1α,25-(OH)2 or D3 vitamin) non-genomic response in cell membranes [8], or in the nucleus as an accessory protein to regulate the transcription factors redox state, i.e., STAT3 [9]. Recently, PDIA3 activity has been correlated with influenza A virus (IAV) replication mechanism. PDIA3 is upregulated in IAV-infected mice playing a key role in the folding of IAV-hemagglutinin [10]. In human hepatocellular carcinoma, PDIA3 downregulation inhibits cell proliferation and, through STAT3 signaling, induces apoptosis in agreement with the observation that PDIA3 knockdown reduces phosphorylated STAT3 and downstream STAT3-related protein levels [11]. PDIA3 has also been associated with human ovarian cancer chemoresistance. Co-treatment of PDIA3-siRNA and paclitaxel inhibits the STAT3 signaling pathway [12]. PDIA3 downregulation inhibits proliferation and increases apoptosis in acute myeloid leukemia [13], glioma [14], and colorectal cancer cells [15].
Despite the increasing number of studies, there is not yet available any chemical tool (inhibitor or activator) to further elucidate the PDIA3 role in the above-mentioned molecular mechanisms. Therefore, the identification of a selective PDIA3 ligand represents an important issue.
In a previous study [16], punicalagin was reported to bind to PDIA3 and to inhibit the reductase activity. Punicalagin (Figure 1) is a bioactive natural compound extracted from pomegranate fruit (Punica granatum) [17], in Terminalia catappa [18], Terminalia myriocarpa [19], and in Combretum mole [20]. It is an ellagitannin belonging to the polyphenol family, and it is soluble in water, unlike many other polyphenols. Punicalagin holds a wide range of pharmacological effects such as anti-inflammatory, hepatoprotective, antioxidant, and anti-cancer effects [21,22]. To investigate the punicalagin specificity in PDIA3 binding among PDI family members, a comparative study between PDIA3 and PDIA1 was performed. PDIA1 is the archetype of the PDI family, is the most abundant PDI in the endoplasmic reticulum [23], and shares a considerable similarity in structure and enzymatic functions with PDIA3 with respect to other PDI family members [24].
Biochemical studies have been carried out in order to compare PDIA3 and PDIA1 interactions and inhibitory effects of punicalagin. Subsequently, molecular dynamics and molecular docking studies have been performed as complementary techniques. Moreover, deep analyses on PDIAs-punicalagin interactions provide insights to set up a rational design for PDIA3 selective inhibitors and to dissect the molecular determinant to modulate the protein activity.

2. Materials and Methods

2.1. Chemicals

Reagents used in this study were purchased from Sigma-Aldrich (Milan, Italy) unless otherwise stated. EDTA (ethylenediamine tetra-acetic acid) 0.5 M solution pH 8.0 was purchased from IBI Scientific (Milan, Italy), Tris(hydroxymethyl)aminomethane for buffer solutions from Merck Millipore (Milan, Italy).

2.2. Recombinant Proteins: Production and Purification

Mature human recombinant PDIA3 protein was expressed and purified according to the procedure described by Trnkova et al. [25]. PDIA3 concentration was spectrophotometrically calculated by means of the extinction coefficient (ε280 in reduced form, 44,810 M−1cm−1).
The pOLR130 plasmid encoding mature human PDIA1 with an N-terminal His6x-tag [26] was generously provided by Dr. Lloyd Ruddock (University of Oulu, Finland). E. coli cells (strain BL21) were transformed with a vector containing the human PDIA1 sequence. After IPTG induction, cells were harvested and lysed. PDIA1 was purified by ammonium sulfate fractionation (30% to 75% saturation), nickel chromatography (Protino Ni-NTA column, Macherey-Nagel, Düren, Germany), followed by anion-exchange chromatography step (Macro-Prep Q column, BioRad, Milan, Italy). SDS-PAGE and Coomassie blue staining were performed to evaluate the protein purification, and concentration was spectrophotometrically estimated using the extinction coefficient (ε280 in reduced form, 45,567 M−1cm−1).

2.3. Measurements of Tryptophan Fluorescence Quenching

Protein-ligand interactions were evaluated according to the fluorescence quenching titration method using a SPEX-FluoroMax spectrofluorometer (Horiba Scientific, Piscataway, NJ, USA). Emission spectra were recorded from 300 to 400 nm with an excitation wavelength set at 290 nm. The emission slit width was fixed to provide a 2 nm bandpass for all experiments. Scan speed was set at 120 nm min−1. A total of 50 µM PDIA1 and PDIA3 stock solutions were reduced by adding DTT 5 mM. Aliquots of freshly reduced proteins were diluted (0.1 µM final concentration) in phosphate-buffered saline (PBS) pH 7.4 containing 0.2 mM EDTA and 0.1 mM DTT in a 10 mm path length quartz fluorescence cuvette under continuous stirring. The titrations were performed by stepwise additions, at 5 min time intervals, of punicalagin solution from 0 to 4 µM (punicalagin 0.2 mM in PBS pH 7.4 freshly prepared from a 5 mM stock solution in water). Blank spectra for each experiment (no protein added) were performed in parallel for background determination. The emission spectrum of the protein was recorded three times, and the average reading values were then used for quenching analyses as already described [27].

2.4. Isothermal Titration Calorimetry

Punicalagin interaction with PDIA3 and PDIA1 was investigated using isothermal titration calorimetry (ITC) with a MicroCal ITC200 (Malvern Instruments, Malvern, UK). Aliquots of freshly reduced protein were diluted (25 µM final concentration) in 20 mM Tris-HCl pH 8.0, 150 mM NaCl and 0.5 mM DTT. Protein solution was loaded in the sample cell (0.2 mL), and the syringe was filled with punicalagin 250 µM dissolved in the same buffer. A total of 20 ligand aliquots (the first aliquot of 0.4 µL in 0.8 s and the other 19 aliquots of 2 µL in 4 s) were gradually added to the sample cell at intervals of 200 s each other to allow baseline recovery. The syringe stirring speed was set at 800 rpm. Protein and punicalagin solutions were thoroughly degassed on the MicroCal’s Thermovac machine before loading to avoid air bubble formation during the titration. Data analyses were performed by Origin software version 8.0 (OriginLab, Northampton, MA, USA). The heat of dilution for the ligand was determined in parallel, and this value was subtracted from the integrated data before curve fitting. The integration of the area under each heat burst curve was determined. All injection heats as a function of the molar ratio were plotted and analyzed using the one set of sites model to provide the dissociation constants.

2.5. Differential Scanning Calorimetry

Protein stability in solution was assayed by differential scanning calorimetry (DSC) with a Microcal VP-DSC (Microcal, Northampton, MA, USA) testing both proteins. A 5 µM protein solution (from a 50 µM stock solution freshly reduced by adding DTT 5 mM) was prepared in 20 mM Tris-HCl pH 8.0, 150 mM NaCl, containing 0.2 mM EDTA, degassed at 25 °C and filled in the measurement cell. The reference cell was filled with protein-free buffer. A heating rate of 60 °C h−1 was set, and the measurement cell worked under an excess pressure of 0.2 kPa to avoid air bubble formation during the heating. The heat absorption during the denaturation process was measured by subtracting the instrumental baseline from the sample trace and normalizing for protein concentration. To test the effect of the ligand on protein stability, analysis was also performed on both proteins in the presence of punicalagin (25 µM final concentration). The melting points (Tm, temperature at which excess heat capacity reaches a maximum) were directly observable from the thermogram. Nevertheless, data deconvolution using the Origin software provided by MicroCal was necessary to obtain Tm values.

2.6. Measurement of PDIAs Disulfide Reductase Activity

PDIAs disulfide reductase activity was monitored by using di-eosin glutathione disulfide (di-E-GSSG) as a fluorogenic substrate of PDIAs. Di-E-GSSG was synthesized by the reaction of eosin isothiocyanate with oxidized glutathione (GSSG) according to the method of Raturi and Mutus [28] with some modifications [29]. Protein activity was evaluated by monitoring the emission fluorescence increase (λem = 545 nm and λex = 525 nm), and the effect of punicalagin (concentration ranging from 0.2 to 50 µM) was tested. Inhibition constants were extrapolated by GraphPad Prism 8.0 software (GraphPad Software, San Diego, CA, USA), plotting the obtained data as logarithm dose-response curves.

3. Results

3.1. Biochemical Studies

3.1.1. Assessment of PDIA-Punicalagin Interactions by Intrinsic Fluorescence Spectroscopy

Quenching of intrinsic fluorescence of both PDIAs upon ligand binding was evaluated. Each protein possesses an intrinsic fluorescence given by aromatic amino acids content, and the tryptophan residues are the dominant source [30]. PDIA3 contains three tryptophan residues, two of them, W56 and W405, are located next to the active sites, respectively, in the a (C57-C60) and a’ (C406-C409) domains, whereas the third tryptophan, W279, is on the b’ domain and is only partially exposed to the solvent [29]. Instead, PDIA1 contains five tryptophan residues, W52, W128 sited on a domain, W364, W396, W407 on the a’ domain. W52 and W396 residues, homologous to W56 and W405 PDIA3 residues, are located next to the redox sites (C53-C56 and C397-C400) and are mostly exposed to the solvent. Considering PDIA1 3D structure, W407 residue appears buried, while W128 and W364 residues are partially and fully exposed, respectively. Upon ligand binding, the fluorescence intensity of both proteins could be quenched, suggesting that an interaction protein-ligand nearby tryptophan residues is taking place (Figure 2A).
Quenching of intrinsic tryptophan fluorescence at 338 nm has been analyzed by Stern–Volmer equation (Fo/F = 1 + KSV[L]), where F and F0 are the fluorescence intensity, respectively with or without ligand (L), and KSV is the Stern–Volmer constant value. Based on obtained results, punicalagin is able to quench the fluorescence of both proteins with a Stern–Volmer constant always greater than 1.0 × 104 M−1 (Figure 2B and Table 1). Considering that the Stern–Volmer constant is equal to the quenching constant multiplied by the average lifetime of the fluorophore (KSV = Kqτ), where the tau value for tryptophan is in the order of 1.0 × 10−8 s [31,32], Kq values are always greater than 1.0 × 1010 M1s−1, the maximum quenching rate for diffusion collision, suggesting that a static interaction occurs [33]. In order to further characterize the binding process, data were analyzed using the equation described by Bi et al. [34] and the reiterative calculation process described by Sun et al. [35], providing an estimation of the Kd values (Table 1). Punicalagin binds both PDIAs in the micro-molar order of magnitude, showing a slightly greater affinity toward PDIA3, which tryptophan residues appear more quenchable upon punicalagin binding. In addition, in both proteins, the intrinsic fluorescence analyzed in reducing conditions appear less affected by punicalagin, suggesting either a minor affinity or less exposure of tryptophan residues for the protein in the reduced conformation.

3.1.2. Analysis of PDIA-Punicalagin Interactions by Isothermal Titration Calorimetry

To a deeper characterization of PDIA-punicalagin interaction, ITC was used as it depicts the standard way to study protein-ligand interaction providing detailed thermodynamic parameters (entropy change ΔS, enthalpy change ΔH, binding constant) in solution. The enthalpy change is generally defined as the energy change resulting from all non-covalent protein-ligand interactions. Instead, the entropy change is given by the solvent entropy change upon ligand binding, as well as conformational entropy change and loss of rotational and translational degrees of freedom upon complex formation [36].
Considering the difficulty of assessing a fully oxidized state, proteins were analyzed in the presence of reducing agents after being freshly reduced. Obtained data confirmed the protein-ligand complex formation in both proteins, with an estimated Kd of 1.0 µM and 1.2 µM for PDIA1 and PDIA3, respectively (Table 2). Kd constant values were quite similar for the two proteins and slightly lower compared to the dissociation constants measured by tryptophan fluorescence quenching analysis. This may suggest that, although both proteins seem to bind punicalagin with a comparable affinity, the binding could differently involve tryptophan residues and influence protein intrinsic fluorescence. Both protein-punicalagin interactions are characterized by negative enthalpy change and positive entropy change (Table 2). Both parameters positively contribute to reducing the Gibbs free energy. Nevertheless, TΔS is much greater than the ΔH factor, meaning that the entropic factor represents the driving force of both interactions. The negative ΔH suggests that non-covalent interactions at the binding surface occur, whereas positive TΔS indicates that the rearrangement of water molecules, upon ligand binding, plays a key role by increasing the solvent degrees of freedom. However, the thermodynamic parameters obtained from ITC data seem to be different for PDIA1 and PDIA3, suggesting a slightly different interaction with punicalagin.

3.1.3. Thermally Induced Transitions of PDIAs upon Punicalagin Binding

Any conformational change induced by punicalagin binding may affect the stability of PDIA3 and PDIA1 proteins, and this was therefore evaluated by means of differential scanning calorimetry (DSC). First, the thermally induced denaturation process of both PDIAs in the reduced form was assessed. The thermal profile of both proteins displayed two endothermic peaks (Figure 3), which may represent different denaturation steps in which, respectively, partial and complete protein denaturation occurs. Deconvolution of data suggests that proteins underwent a two-step denaturation process (non-two-state transition). Despite the PDIA3 and PDIA1 sharing a suitable structure similarity, PDIA3 showed lower stability with the transition temperature (Tm) of the second step at 53.9 °C, with respect to 61.0 °C of PDIA1. Instead, the temperature difference between the first two peaks is less marked, with a Tm = 44.9 °C for PDIA3 with respect to 46.0 °C of PDIA1 (Figure 3). As reported, switching from an open to a closed conformation occurs in PDIA1, with a greater closeness between a and a’ domains in the reduced form [37]. It can be hypothesized that the first thermal transition could be the result of a partial protein unfolding leading to an open conformation similar to that adopted in the oxidized form.
Tryptophan fluorescence quenching results did not show a shift in the intrinsic fluorescence emission upon punicalagin binding. Therefore, it could be deduced that no noticeable conformational change occurs directly upon punicalagin binding. However, punicalagin binding influenced the denaturation profile of both proteins, negatively affecting PDIA3 stability. In fact, the Tm of the two endothermic peaks of PDIA3 decreases from 44.9 to 40.4 °C and from 53.9 to 51.6 °C, respectively. In the case of PDIA1, Tm values of the second peak were quite similar, 61 and 60 °C, while Tm of the first peak decreased from 46 to 42 °C. In both cases, the first thermal transition was more influenced by punicalagin binding, but only in PDIA1 did the presence of punicalagin greatly reduce the peak area associated with the first transition, corresponding to the enthalpy variation (ΔH), and broadened the second endothermic peak (Figure 3). This suggested that the binding of punicalagin can affect the stability of the two proteins differently.

3.1.4. Punicalagin Effect on PDIA3 and PDIA1 Reductase Activity

Punicalagin effects have also been studied on the enzymatic activity of PDIA3 and PDIA1. It is well known that both proteins catalyze oxidoreductase and isomerase reactions involving thiol groups and disulfide bonds [2,38]. PDIs’ reductase activity was evaluated using di-E-GSSG as fluorogenic substrate and was calculated as the initial velocity in fluorescence increased due to the reduction in the fluorogenic substrate. The previous investigation identified punicalagin as an inhibitor of PDIA3 reductase activity [16], as a non-competitive inhibitor with a Ki of 0.39 µM, and this was herein confirmed. For comparison purposes, the punicalagin inhibitory effect was extended to PDIA1. Results were analyzed as logarithmic dose response, and the half-maximal inhibitory concentrations (IC50) were evaluated (Figure 4). Obtained data show that punicalagin has a better inhibitory effect toward PDIA3 (IC50 1.5 ± 0.2 µM) compared to PDIA1 (IC50 6.1 ± 0.6 µM).

3.2. Computational Studies

To rationalize the herein presented experimental data and provide the putative binding mode for punicalagin on PDIA3 and PDIA1, a computational chemistry procedure was conducted combining molecular dynamics (MD) simulations with molecular docking simulations as summarized in the scheme depicted in Figure 5:
  • Modeling of the PDIAs crystal structures in oxidized and reduced states;
  • MD simulations of the modeled PDIAs;
  • Analysis of the MD trajectories and conformations’ sampling for the subsequent ensemble molecular docking (cross-docking) [39] simulations;
  • α and β-punicalagin MD-based conformational analysis for the molecular docking into a and a’ domains of PDIAs sampled conformations;
  • Analysis of the molecular docking results by means of statistical techniques and docking score ranking allowed binding poses selection for final rescoring with MM/GBSA [40].

3.2.1. Molecular Dynamics Simulations

As no experimental PDIA3-ligand complex is yet available, MD simulations were run to sample PDIA3 conformational flexibility to overcome the lack of ligand-protein induced-fit sampling frames to be used in molecular docking simulations to obtain insight on the putative punicalagin binding mode [41]. Starting from the PDIA3 crystal structure (PDB ID:3F8U) [3], the molecular systems were built and then subjected to MD simulation following the procedures reported in the “Computational Methods” section (Supplementary Material Sections S1.1 and S1.2). The PDIA3 and PDIA3-Tap molecular systems were modeled in both the reduced (Red) and oxidized (Ox) forms for either apoproteins or complexed with tapasin (PDIA3Ox, PDIA3Ox-Tap, PDIA3Red, and PDIA3Red-Tap). A total of 400 ns MD simulation trajectories were firstly analyzed, collecting the protein backbone root mean square deviation (RMSD) to investigate protein stability. Tapasin-containing complexes (PDIA3Ox-Tap and PDIA3Red-Tap) showed a stable RMSD profile with values ranging from 3 to 5 Å, while the free PDIA3Ox and PDIA3Red systems returned an increasing trend up to 12–14 Å (Figure S1A). This difference was also consistent with the RMSD probability density function plot (Figure S1B), likely due to the domains’ high mobility. In fact, during the PDIA3s’ MD simulations, the proteins were free to move, while in the complexed systems, tapasin induced some structural constraints on both a and a’ domains resulting in lower RMSD values. As punicalagin is also a PDIA1 inhibitor, the calculations were replicated on the PDIA1 system. To this, the PDIA1 starting crystal structure (PDB ID:4EL1) [37] was prepared as described in the “Computational Methods” section (see Supplementary Material Sections S1.1 and S1.2), solvated, and subjected to MD simulation for either oxidized (PDIA1Ox) or reduced state (PDIA1Red). Although experimentally available (PDB ID:4EKZ) [37], for consistency with the PDIA3 system, the PDIA1 reduced state was modeled starting from the oxidized PDIA1 crystal. The RMSD between the reduced PDIA1 crystal structure (PDB ID: 4EKZ) and the modeled one after the initial MD equilibration was 4.03 Å, while between the oxidized (PDB ID: 4EL1) and the reduced (PDB ID: 4EKZ) PDIA1 crystal structures was 6.87 Å. The RMSD fluctuations range (12–13 Å) and its trend observed from the MD trajectory analysis overlapped those observed for PDIA3Ox and PDIA3Red (Figure S1). Notably, PDIA1Red reached the equilibrium after 10 ns at an RMSD value of ~11 Å. Upon deeper analysis, the distance between a and a’ domains along the PDIA1 simulations was collected (Figure S4) and revealed PDIA1Red switching from an open to a closed conformation during the first 10 ns (Figure S5). The latter agreed with experimental data showing a greater closeness between PDIA1 a and a’ domains in the reduced form [37]. Further analysis on PDIA3 and PDIA1 trajectories are reported in Supplementary Material Section S2.1.

3.2.2. Molecular Docking Simulations

As introduced above, PDIAs trajectories were analyzed, and a series of snapshots (60 for PDIA3 and 30 for PDIA1) were sampled (see Supplementary Material Sections S1.3 and S2.2) to run molecular docking simulations. As α and β-punicalagin (Scheme S1) are characterized by a cyclized, highly constrained chemical structure and considering the smina limitations (the cycles are treated rigidly), molecular docking simulations were conducted using a rigid body docking procedure. To fulfill the lack of conformational flexibility, MD simulations of α and β-punicalagin were carried out to sample punicalagin conformations (details in Supplementary Material Sections S1.4 and S2.3).
The selected α and β-punicalagin conformations (40 for each) were used as keys and docked onto the a and a’ domains (locks) of PDIA3 sampled conformations. Among the resulting 24,000 docked poses for each domain, conformations were selected using either the highest Vinardo score (10 conformations) or a statistical approach applying kernel density estimation (KDE) on poses’ heavy atoms cartesian coordinates first two principal components (10 conformations). For PDIA3 a domain (Figure S11A), the α and β-punicalagin conformations selected using the two approaches displayed partial overlap, returning as the mean distance between poses’ center of mass (MDCOM) a value of 7.4 Å. For the PDIA3 a’ domain (Figure 6A), the conformations sampled using the two approaches showed a very suitable overlap, being all the predicted binding modes in the same pocket occupying the same volume (MDCOM 2.9 Å), indicating a convergent way to select the docked α and β-punicalagin conformations. The 20 selected poses were used for free energy calculations for final rescoring and selecting the most likely binding conformation.
The same docking procedure was applied for PDIA1 a and a’ domains of 30 selected MD snapshots, returning 12,000 binding poses for each PDIA1 catalytic domain. Again, the binding pose selection for free energy calculations was accomplished by means of the smina’s Vinardo score ranking and KDE algorithm. Selected docking results on PDIA1 a domain (Figure S14A) showed two distinct conformations clusters matching the two selection approaches (MDCOM 8.7 Å). Selected docking results on the a’ domain (Figure 6B) were mainly grouped into a cluster occupying a domain hydrophobic pocket (MDCOM 7.4 Å), the same identified as punicalagin binding site on PDIA3 a’ domain. However, few poses were located outside this cluster, at the a’/b’ domains interface (Figure 6B and Figure S15).

3.2.3. Punicalagin Binding Mode Selection by Free Energy Calculations

The selected α and β-punicalagin binding poses from molecular docking simulations into a and a’ domains were merged with the respective locks obtaining a number of complexes to calculate the binding free energy (ΔGbind) by means of the MM/GBSA method.
In agreement with experimental data [16], lower estimated ΔGbind values were obtained for punicalagin complexed with PDIA3Ox a’ domain. The top-ranked calculated ΔGbind for β-punicalagin on PDIA3 a domain was −39.2 kcal mol−1 while on PDIA3 a’ domain was −49.9 kcal mol−1. In both domains, β-punicalagin showed a higher affinity with respect to the α epimer. The binding conformation for β-punicalagin was directly selected for the a’ domain as the top-ranked one (Figure 7A), showing a lower ΔGbind of about 8 kcal mol−1 than the α epimer. This result lets us uniquely define a reliable binding mode for β-punicalagin on PDIA3 a’ domain where the ligand is settled in a hydrophobic pocket constituted by Val378, Val380, Val381, and Val382. Main polar interactions are with Val381, Lys433, Val378, Gly376, Asp435, Ser373 backbone atoms and with Asn439 side chain (Figure 7A).
Binding free energy calculations indicated PDIA1 a’ as the punicalagin preferred binding domain and pointed β-punicalagin as the favored epimer for both catalytic domains. Differently for PDIA3, β-punicalagin affinity was calculated to be higher for PDIA1Red. The top-ranked calculated ΔGbind for β-punicalagin on PDIA1 a domain was −42.7 kcal mol−1 while for the a’ domain was −56.8 kcal mol−1. The top-ranked binding mode on the PDIA1 a’ domain showed a ΔGbind difference of about 11 kcal mol-1 with the second one. These results let to hypothesize a reliable binding mode for β-punicalagin on PDIA1Red a’ domain where the ligand is inserted in a hydrophobic pocket at the a’/b’ interface, characteristic of the closed conformation, constituted by Pro246, Phe304, Phe305, Phe240, Phe349, Leu258, Ile318, Ile301, Arg300 and Trp396, which interacts with ellagic acid moiety trough a π–π stacking. Main polar contacts are established with Arg300, Asp297 (Figure 7B and Figure S17).

3.2.4. Punicalagin Binding Mode on PDIA3 Refinement by MD Simulation

As the predicted β-punicalagin binding mode on PDIA3 was not found in close contact with Trp405 quenching residue [16], an investigation on how β-punicalagin could modulate PDIA3 disulfide reductase activity was conducted by means of a 100 ns MD simulation. Trajectory analysis then focused on Trp405 solvent accessible surface area (SASA) collection. The β-punicalagin/PDIA3Ox complex (PDIA3Ox-Pun) showed a different Trp405 SASA trend with respect to the other trajectories (Figure 8A), showing a Trp405 burial in the second half of the simulation. A visual inspection indicated that β-punicalagin was not responsible for the PDIA3Ox-Pun SASA fluctuation trend. Nevertheless, a and a’ domains were approaching each to the other to make contacts adopting a “closed” or “collapsed” conformation (Figure 8C,D) never observed during simulations of either PDIA3Ox or PDIA3Red. Analysis of the distance fluctuations between the a and a’ domains clearly showed the under-closing conformation of PDIA3Ox-Pun, reaching distance values less than 5 Å in the second half of the simulation (Figure 8B). RMSD analysis of PDIA3Ox-Pun (Figure S13A) returned a similar trend to PDIA3Ox and PDIA3Red while the root mean square fluctuation (RMSF) (Figure S13B) indicate reduced flexibility of PDIA3Ox-Pun complex, with lower RMSF values than PDIA3Ox and PDIA3Red. These results could account in part for the quenching experiments and let to speculate on an inhibition mechanism based on the modulation of PDIA3 flexibility regulating the availability of the catalytic sites.

3.2.5. PDIA3/PDIA1 β-Punicalagin Binding Mode Comparison

Sequence identity between full-length PDIA3 and PDIA1 is 36%, while identity between PDIA3 and PDIA1 a domains is 50% and between PDIA3 and PDIA1 a’ domains is 45%, being the a’ the less conserved one. Sequence analysis was also used to visually inspect PDIs binding pocket conservation (Figure 9). Results on PDIA3 showed punicalagin settling on a conserved hydrophobic pocket approaching the 358–376 loop (Figure 9A). The latter section, corresponding to the x linker in PDIA1 [37], is also where the nonconserved interacting residues are located. Results on PDIA1 showed punicalagin binding into a hydrophobic pocket peculiar of the PDIA1Red closed conformation (Figure 9B). Although a few residues are conserved (Glu242PDIA1/Glu249PDIA3, Thr428PDIA1/Thr437PDIA3 and Trp396PDIA1/Trp405PDIA3), this pocket is not present in PDIA3 and might be responsible for the punicalagin/PDIA1 main interactions (i.e., π–π stacking with the ellagic acid moiety). The proposed β-punicalagin binding modes on PDIA3 and PDIA1 are in two different binding pockets suggesting different effects of this molecule on proteins’ structure, in agreement with presented experimental data, as well as a different inhibition mechanism. In the case of PDIA3, β-punicalagin could modulate protein flexibility and consequently the equilibrium between the “open” and the “closed” conformations (as indicated by PDIA3Ox-Pun simulation), while PDIA1 could be targeted by β-punicalagin in the “closed” state exerting a stabilizing effect through binding in the above-described pocket.
Indeed, PDIA1 flexibility was an already discussed theme in scientific literature. PDIA1 has been subjected to redox-dependent structural changes and is characterized by high interdomain mobility [37,42,43]. The b’xa’ domains are the most dynamic, and this is crucial to properly interact with substrates while exerting the chaperone function [44]. Reduced PDIA1 is more stable than the oxidized one and more resilient to thermal stress leading to hypothesize a rearrangement in a “closed” conformation in which the x linker plays a key role in the a’ domain positioning [45]. These reports agree with the herein reported PDIA1Red MD simulation (see Supplementary Material Section S2.1), in which PDIA1 adopted a “closed” conformation during the first 10 ns. Such conformation could not be able to properly interact with substrates, and β-punicalagin might stabilize it, avoiding PDIA1 opening. The binding mode proposed for PDIA3 presented β-punicalagin interacting with a hydrophobic pocket approaching the flexible loop connecting a’ and b’ domains, corresponding to PDIA1 x linker. This could explain the a/a’ approach observed in PDIA3Ox-Pun MD simulation and then the inhibition mechanism.

4. Discussion

In a previous work, it was shown that punicalagin, found in forms α and β in pomegranates as well as other natural sources, was able to bind to PDIA3 and inhibited its disulfide reductase activity. Herein we provide evidence that punicalagin can also bind to PDIA1, the main expressed form of PDI. In this comparative study, the affinity and the effect of punicalagin binding on both proteins were evaluated, and computational studies were used to identify the putative PDIA1 and PDIA3 binding sites. Although the findings indicated punicalagin capable of binding both proteins with a similar affinity, further biochemical evidence revealed several differences due to the effect of binding. Punicalagin was able to better quench the fluorescence of PDIA3 and showed a higher inhibitory effect on protein reductase activity for PDIA3 than that on PDIA1. Additionally, punicalagin differently affects the thermal denaturation profile of both proteins. The decrease in Tm for PDIA3Ox-Pun complex suggested a preferential binding for the unfolded protein conformation, while PDIA1Ox-Pun complex was characterized by a marked reduction in the peak area associated with the first thermal transition (likely the “closed-open” conformational switch that the protein in the reduced form undergoes).
Computational investigations identified PDIA1 and PDIA3 a’ redox domains as the most favorable punicalagin binding sites, which sequence and structure mostly differ between the two proteins. Furthermore, punicalagin binding sites were suggested to be in two different pockets on PDIA3 and PDA1, in agreement with the quenching and thermal denaturation profiles.
The proposed PDIA3-punicalagin binding was predicted to settle in a hydrophobic pocket on the a’ domain, and the ligand showed a preferential binding for the oxidized form. MD simulations indicate reduced flexibility of the PDIA3Ox-Pun complex, modulating the availability of the catalytic sites. Although a direct interaction with tryptophan residues was not evident, these results could account in part for the higher intrinsic fluorescence quenching observed for PDIA3-punicalagin interaction and let to speculate on an inhibition mechanism based on the modulation of PDIA3 redox domain flexibility.
The binding mode proposed for PDIA1 showed punicalagin in a hydrophobic pocket at the a’/b’ interface, experimentally observed in the closed conformation adopted by the protein in the reduced form (PDIA1Red). This agreed with the variation in the PDIA1 thermal profile induced by punicalagin. In fact, the area associated with the first endothermic peak, probably related to the closed-open transition, is greatly reduced while the second peak appears broadened. In this case, the inhibitory mechanism could be mainly based on the stabilization of the protein in the reduced form.
ITC thermodynamic data agreed with the proposed punicalagin binding modes. The PDIA1-punicalagin complex was characterized by a higher enthalpy contribution, likely resulting from several interactions stabilizing the binding (i.e., π–π stacking with the ellagic acid moiety). On the other hand, binding of punicalagin to PDIA3, as seen in PDIA3Ox-Pun MD simulation, could modify the protein flexibility and, therefore, solvent accessibility, thus explaining the higher entropic contribution in this interaction.
PDI family consists of a large number of versatile proteins that can catalyze the oxidation/reduction/exchange of disulfide bonds within substrates. PDIs are involved in a wide range of cellular functions and are differently associated with several pathological conditions, including neurodegeneration and cancer. Thus, the search for specific and selective inhibitors/modulators of their activity is important to allow an aimed therapy. Here is presented evidence that punicalagin, a naturally occurring molecule, is able to differently bind and affect the structure and activity of PDIA1 and PDIA3, two of the most abundant proteins of the family. Molecular docking and simulation analyses provided information on the different punicalagin binding modes on PDIA1 and PDIA3. These findings will be useful to develop punicalagin-based ligands and to set up a rational design for PDIs selective inhibitors.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biomedicines9111533/s1. Figure S1: Backbone RMSD over simulation time and KDE distribution plot of RMSD values. Figure S2: PDIA3 RMSD and RMSF domain-wise plots. Figure S3: PDIA1 RMSD and RMSF domain-wise plots. Figure S4: Distance between PDIA1Ox and PDIA1Red a and a’ domains over time. Figure S5: PDIA1Red in open and closed conformation. Figure S6: PCA representations of PDIA3 trajectories. Figure S7: PDIA3 conformations selected from MD simulations. Figure S8: PCA representations of PDIA1 trajectories. Figure S9: PDIA1 conformations selected from MD simulations. Figure S10: α and β-punicalagin in all-axial and all-equatorial conformations. Figure S11: Selected α and β-punicalagin binding modes for a and a’ PDIA3 domains. Figure S12: Minimized binding modes of β-punicalagin on PDIA3 a domain. Figure S13: RMSD and RMSF plots comparing all PDIA3 trajectories. Figure S14: Selected α and β-punicalagin binding modes for a and a’ PDIA1 domains. Figure S15: Selected α and β-punicalagin binding modes for a’ PDIA1 domains displayed on the whole PDIA1 structure. Figure S16: Top-ranked binding mode of β-punicalagin on PDIA1 a domain. Figure S17: Lowest energy predicted binding mode for β-punicalagin on PDIA1 a’ domain. Scheme S1: α and β-punicalagin 2D depictions. Figures S2, S3, S6, S7, S8, S9, S10, S12 and S16 are cited in the supplementary materials. References [3,16,37,40,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73] are cited in the supplementary materials.

Author Contributions

Conceptualization, R.R. and F.A.; methodology, G.P., L.A., L.C., M.S., M.M. and E.R.; software, G.P., L.A., L.C., M.S. and R.R.; validation, G.P., L.A., L.C., R.R. and F.A.; formal analysis, G.P., L.A., L.C., R.R. and F.A.; investigation, G.P., L.A., L.C., M.S., M.M. and E.R.; resources, R.R. and F.A.; data curation, G.P. and L.A.; writing—original draft preparation, G.P., L.A., L.C., R.R., M.S., M.M., E.R. and F.A.; writing—review and editing, G.P., L.A., L.C., R.R., M.M., E.R. and F.A.; visualization, G.P. and L.A.; supervision, R.R. and F.A.; project administration, R.R. and F.A.; funding acquisition, F.A. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by Sapienza Ateneo grants (F.A.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and the Supplementary Materials file.

Acknowledgments

The authors thank Lloyd Ruddock (University of Oulu, Finland) for providing pOLR130 plasmid encoding mature human PDIA1.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

PDIprotein disulfide isomerase
IAVinfluenza A virus
ITCisothermal titration calorimetry
DSCdifferential scanning calorimetry
MDmolecular dynamics
RMSDroot mean square deviation
KDEkernel density estimation
MDCOMmean distance between poses’ center of mass
SASAsolvent accessible surface area
RMSFroot mean square fluctuation

References

  1. Hettinghouse, A.; Liu, R.; Liu, C. Multifunctional Molecule ERp57: From Cancer to Neurodegenerative Diseases. Pharmacol. Ther. 2018, 181, 34–48. [Google Scholar] [CrossRef]
  2. Frickel, E.M.; Frei, P.; Bouvier, M.; Stafford, W.F.; Helenius, A.; Glockshuber, R.; Ellgaard, L. ERp57 Is a Multifunctional Thiol-Disulfide Oxidoreductase. J. Biol. Chem. 2004, 279, 18277–18287. [Google Scholar] [CrossRef] [Green Version]
  3. Dong, G.; Wearsch, P.A.; Peaper, D.R.; Cresswell, P.; Reinisch, K.M. Insights into MHC Class I Peptide Loading from the Structure of the Tapasin-ERp57 Thiol Oxidoreductase Heterodimer. Immunity 2009, 30, 21–32. [Google Scholar] [CrossRef] [Green Version]
  4. Silvennoinen, L.; Myllyharju, J.; Ruoppolo, M.; Orrù, S.; Caterino, M.; Kivirikko, K.I.; Koivunen, P. Identification and Characterization of Structural Domains of Human ERp57: Association with Calreticulin Requires Several Domains. J. Biol. Chem. 2004, 279, 13607–13615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Turano, C.; Gaucci, E.; Grillo, C.; Chichiarelli, S. ERp57/GRP58: A Protein with Multiple Functions. Cell. Mol. Biol. Lett. 2011, 16, 539–563. [Google Scholar] [CrossRef] [PubMed]
  6. Jessop, C.E.; Chakravarthi, S.; Garbi, N.; Hämmerling, G.J.; Lovell, S.; Bulleid, N.J. ERp57 Is Essential for Efficient Folding of Glycoproteins Sharing Common Structural Domains. EMBO J. 2007, 26, 28–40. [Google Scholar] [CrossRef] [PubMed]
  7. Santos, S.G.; Campbell, E.C.; Lynch, S.; Wong, V.; Antoniou, A.N.; Powis, S.J. Major Histocompatibility Complex Class I-ERp57-Tapasin Interactions within the Peptide-Loading Complex. J. Biol. Chem. 2007, 282, 17587–17593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Nemere, I.; Garbi, N.; Hammerling, G.; Hintze, K.J. Role of the 1,25D 3-MARRS Receptor in the 1,25(OH) 2D 3-Stimulated Uptake of Calcium and Phosphate in Intestinal Cells. Steroids 2012, 77, 897–902. [Google Scholar] [CrossRef] [PubMed]
  9. Chichiarelli, S.; Gaucci, E.; Ferraro, A.; Grillo, C.; Altieri, F.; Cocchiola, R.; Arcangeli, V.; Turano, C.; Eufemi, M. Role of ERp57 in the Signaling and Transcriptional Activity of STAT3 in a Melanoma Cell Line. Arch. Biochem. Biophys. 2010, 494, 178–183. [Google Scholar] [CrossRef] [PubMed]
  10. Chamberlain, N.; Korwin-Mihavics, B.R.; Nakada, E.M.; Bruno, S.R.; Heppner, D.E.; Chapman, D.G.; Hoffman, S.M.; van der Vliet, A.; Suratt, B.T.; Dienz, O.; et al. Lung Epithelial Protein Disulfide Isomerase A3 (PDIA3) Plays an Important Role in Influenza Infection, Inflammation, and Airway Mechanics. Redox Biol. 2019, 101129. [Google Scholar] [CrossRef]
  11. Kondo, R.; Ishino, K.; Wada, R.; Takata, H.; Peng, W.X.; Kudo, M.; Kure, S.; Kaneya, Y.; Taniai, N.; Yoshida, H.; et al. Downregulation of Protein Disulfide-Isomerase A3 Expression Inhibits Cell Proliferation and Induces Apoptosis through STAT3 Signaling in Hepatocellular Carcinoma. Int. J. Oncol. 2019, 54, 1409–1421. [Google Scholar] [CrossRef] [Green Version]
  12. Li, S.; Zhao, X.; Chang, S.; Li, Y.; Guo, M.; Guan, Y. ERp57-Small Interfering RNA Silencing Can Enhance the Sensitivity of Drug-Resistant Human Ovarian Cancer Cells to Paclitaxel. Int. J. Oncol. 2019, 54, 249–260. [Google Scholar] [CrossRef] [Green Version]
  13. Ye, Q.; Fu, P.; Dou, J.; Wang, N. Downregulation of PDIA3 Inhibits Proliferation and Invasion of Human Acute Myeloid Leukemia Cells. OncoTargets Ther. 2018, 11, 2925–2935. [Google Scholar] [CrossRef] [Green Version]
  14. Zou, H.; Wen, C.; Peng, Z.; Shao, Y.Y.; Hu, L.; Li, S.; Li, C.; Zhou, H.H. P4HB and PDIA3 Are Associated with Tumor Progression and Therapeutic Outcome of Diffuse Gliomas. Oncol. Rep. 2018, 39, 501–510. [Google Scholar] [CrossRef]
  15. Yang, Z.; Liu, J.; Shi, Q.; Chao, Y.; Di, Y.; Sun, J.; Zhang, J.; Huang, L.; Guo, H.; He, C. Expression of Protein Disulfide Isomerase A3 Precursor in Colorectal Cancer. OncoTargets Ther. 2018, 11, 4159–4166. [Google Scholar] [CrossRef] [Green Version]
  16. Giamogante, F.; Marrocco, I.; Cervoni, L.; Eufemi, M.; Chichiarelli, S.; Altieri, F. Punicalagin, an Active Pomegranate Component, Is a New Inhibitor of PDIA3 Reductase Activity. Biochimie 2018, 147, 122–129. [Google Scholar] [CrossRef]
  17. Khatib, M.; Innocenti, M.; Giuliani, C.; Al-Tamimi, A.; Romani, A.; Mulinacci, N. Mesocarp and Exocarp of Laffan and Wonderful Pomegranate Varieties: By-Products as a Source of Ellagitannins. Int. J. Food Nutr. Sci. 2017, 4, 60–66. [Google Scholar] [CrossRef] [Green Version]
  18. Chen, P.S.; Li, J.H.; Liu, T.Y.; Lin, T.C. Folk Medicine Terminalia Catappa and Its Major Tannin Component, Punicalagin, Are Effective against Bleomycin-Induced Genotoxicity in Chinese Hamster Ovary Cells. Cancer Lett. 2000, 152, 115–122. [Google Scholar] [CrossRef]
  19. Marzouk, M.S.A.; El-Toumy, S.A.A.; Moharram, F.A.; Shalaby, N.M.M.; Ahmed, A.A.E. Pharmacologically Active Ellagitannins from Terminalia Myriocarpa. Planta Med. 2002, 68, 523–527. [Google Scholar] [CrossRef] [PubMed]
  20. Asres, K.; Bucar, F.; Knauder, E.; Yardley, V.; Kendrick, H.; Croft, S.L. In Vitro Antiprotozoal Activity of Extract and Compounds from the Stem Bark of Combretum Molle. Phyther. Res. 2001, 15, 613–617. [Google Scholar] [CrossRef] [PubMed]
  21. Ismail, T.; Sestili, P.; Akhtar, S. Pomegranate Peel and Fruit Extracts: A Review of Potential Anti-Inflammatory and Anti-Infective Effects. J. Ethnopharmacol. 2012, 143, 397–405. [Google Scholar] [CrossRef] [PubMed]
  22. Turrini, E.; Ferruzzi, L.; Fimognari, C. Potential Effects of Pomegranate Polyphenols in Cancer Prevention and Therapy. Oxid. Med. Cell. Longev. 2015, 2015, 938475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wang, L.; Wang, X.; Wang, C.C. Protein Disulfide-Isomerase, a Folding Catalyst and a Redox-Regulated Chaperone. Free Radic. Biol. Med. 2015, 83, 305–313. [Google Scholar] [CrossRef] [PubMed]
  24. Kozlov, G.; Määttänen, P.; Thomas, D.Y.; Gehring, K. A Structural Overview of the PDI Family of Proteins. FEBS J. 2010, 277, 3924–3936. [Google Scholar] [CrossRef] [PubMed]
  25. Coppari, S.; Altieri, F.; Ferraro, A.; Chichiarelli, S.; Eufemi, M.; Turano, C. Nuclear Localization and DNA Interaction of Protein Disulfide Isomerase ERp57 in Mammalian Cells. J. Cell. Biochem. 2002, 85, 325–333. [Google Scholar] [CrossRef]
  26. Nguyen, V.D.; Saaranen, M.J.; Karala, A.R.; Lappi, A.K.; Wang, L.; Raykhel, I.B.; Alanen, H.I.; Salo, K.E.H.; Wang, C.C.; Ruddock, L.W. Two Endoplasmic Reticulum PDI Peroxidases Increase the Efficiency of the Use of Peroxide during Disulfide Bond Formation. J. Mol. Biol. 2011, 406, 503–515. [Google Scholar] [CrossRef]
  27. Giamogante, F.; Marrocco, I.; Romaniello, D.; Eufemi, M.; Chichiarelli, S.; Altieri, F. Comparative Analysis of the Interaction between Different Flavonoids and PDIA3. Oxid. Med. Cell. Longev. 2016, 2016, 4518281. [Google Scholar] [CrossRef]
  28. Raturi, A.; Mutus, B. Characterization of Redox State and Reductase Activity of Protein Disulfide Isomerase under Different Redox Environments Using a Sensitive Fluorescent Assay. Free Radic. Biol. Med. 2007, 43, 62–70. [Google Scholar] [CrossRef]
  29. Trnková, L.; Ricci, D.; Grillo, C.; Colotti, G.; Altieri, F. Green Tea Catechins Can Bind and Modify ERp57/PDIA3 Activity. Biochim. Biophys. Acta Gen. Subj. 2013, 1830, 2671–2682. [Google Scholar] [CrossRef]
  30. Ghisaidoobe, A.B.T.; Chung, S.J. Intrinsic Tryptophan Fluorescence in the Detection and Analysis of Proteins: A Focus on Förster Resonance Energy Transfer Techniques. Int. J. Mol. Sci. 2014, 15, 22518–22538. [Google Scholar] [CrossRef]
  31. Lakowicz, J.R.; Weber, G. Quenching of Fluorescence by Oxygen. Probe for Structural Fluctuations in Macromolecules. Biochemistry 1973, 12, 4161–4170. [Google Scholar] [CrossRef]
  32. Lakowicz, J.R. Principles of Fluorescence Spectroscopy; Quenching of Fluorescence; Springer: New York, NY, USA, 2006; pp. 277–330. [Google Scholar] [CrossRef]
  33. Ware, W.R. Oxygen Quenching of Fluorescence in Solution: An Experimental Study of the Diffusion Process. J. Phys. Chem. 1962, 66, 455–458. [Google Scholar] [CrossRef]
  34. Bi, S.; Song, D.; Tian, Y.; Zhou, X.; Liu, Z.; Zhang, H. Molecular Spectroscopic Study on the Interaction of Tetracyclines with Serum Albumins. Spectrochim. Acta Part. A Mol. Biomol. Spectrosc. 2005, 61, 629–636. [Google Scholar] [CrossRef]
  35. Sun, Y.; Zhang, H.; Sun, Y.; Zhang, Y.; Liu, H.; Cheng, J.; Bi, S.; Zhang, H. Study of Interaction between Protein and Main Active Components in Citrus aurantium L. by Optical Spectroscopy. J. Lumin. 2010, 130, 270–279. [Google Scholar] [CrossRef]
  36. Du, X.; Li, Y.; Xia, Y.L.; Ai, S.M.; Liang, J.; Sang, P.; Ji, X.L.; Liu, S.Q. Insights into Protein–Ligand Interactions: Mechanisms, Models, and Methods. Int. J. Mol. Sci. 2016, 17, 144. [Google Scholar] [CrossRef]
  37. Wang, C.; Li, W.; Ren, J.; Fang, J.; Ke, H.; Gong, W.; Feng, W.; Wang, C.C. Structural Insights into the Redox-Regulated Dynamic Conformations of Human Protein Disulfide Isomerase. Antioxid. Redox Signal. 2013, 19, 36–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Khan, H.A.; Mutus, B. Protein Disulfide Isomerase a Multifunctional Protein with Multiple Physiological Roles. Front. Chem. 2014, 2, 70. [Google Scholar] [CrossRef] [Green Version]
  39. Ragno, R.; Frasca, S.; Manetti, F.; Brizzi, A.; Massa, S. HIV-Reverse Transcriptase Inhibition: Inclusion of Ligand-Induced Fit by Cross-Docking Studies. J. Med. Chem. 2005, 48, 200–212. [Google Scholar] [CrossRef] [PubMed]
  40. Kollman, P.A.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W.; et al. Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models. Acc. Chem. Res. 2000, 33, 889–897. [Google Scholar] [CrossRef]
  41. Bastos, T.M.; Botelho, M.; Soares, P.; Franco, C.H.; Alc, L.; Antonini, L.; Sabatino, M.; Mautone, N.; Freitas-junior, L.H.; Moraes, C.B.; et al. Identification of Inhibitors to Trypanosoma Cruzi Sirtuins Based on Compounds Developed to Human Enzymes. Int J. Mol. Sci. 2020, 21, 3659. [Google Scholar] [CrossRef] [PubMed]
  42. Freedman, R.B.; Desmond, J.L.; Byrne, L.J.; Heal, J.W.; Howard, M.J.; Sanghera, N.; Walker, K.L.; Wallis, A.K.; Wells, S.A.; Williamson, R.A.; et al. ‘Something in the Way She Moves’: The Functional Significance of Flexibility in the Multiple Roles of Protein Disulfide Isomerase (PDI). Biochim. Biophys. Acta Proteins Proteom. 2017, 1865, 1383–1394. [Google Scholar] [CrossRef] [Green Version]
  43. Römer, R.A.; Wells, S.A.; Emilio Jimenez-Roldan, J.; Bhattacharyya, M.; Vishweshwara, S.; Freedman, R.B. The Flexibility and Dynamics of Protein Disulfide Isomerase. Proteins Struct. Funct. Bioinform. 2016, 84, 1776–1785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wang, C.; Chen, S.; Wang, X.; Wang, L.; Wallis, A.K.; Freedman, R.B.; Wang, C.C. Plasticity of Human Protein Disulfide Isomerase: Evidence for Mobility around the x-Linker Region and Its Functional Significance. J. Biol. Chem. 2010, 285, 26788–26797. [Google Scholar] [CrossRef] [Green Version]
  45. Guyette, J.; Evangelista, B.; Tatulian, S.A.; Teter, K. Stability and Conformational Resilience of Protein Disulfide Isomerase. Biochemistry 2019, 58, 3572–3584. [Google Scholar] [CrossRef] [PubMed]
  46. Webb, B.; Sali, A. Comparative protein structure modeling using MODELLER. Curr. Protoc. Bioinform. 2016, 54, 5.6.1–5.6.37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. The PyMOL Molecular Graphics System, Version 2.0; Schrödinger, LLC: New York, NY, USA. Available online: https://pymol.org/2/ (accessed on 30 September 2021).
  48. Olsson, M.H.M.; SØndergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical p K a Predictions. J. Chem. Theory Comput. 2011, 7, 525–537. [Google Scholar] [CrossRef]
  49. Darby, N.J.; Creighton, T.E. Functional Properties of the Individual Thioredoxin-like Domains of Protein Disulfide Isomerase. Biochemistry 1995, 34, 11725–11735. [Google Scholar] [CrossRef]
  50. Eastman, P.; Swails, J.; Chodera, J.D.; McGibbon, R.T.; Zhao, Y.; Beauchamp, K.A.; Wang, L.P.; Simmonett, A.C.; Harrigan, M.P.; Stern, C.D.; et al. OpenMM 7: Rapid Development of High Performance Algorithms for Molecular Dynamics. PLoS Comput. Biol. 2017, 13, 1–17. [Google Scholar] [CrossRef]
  51. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. Ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from Ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  52. Izadi, S.; Anandakrishnan, R.; Onufriev, A.V. Building Water Models: A Different Approach. J. Phys. Chem. Lett. 2014, 5, 3863–3871. [Google Scholar] [CrossRef] [Green Version]
  53. Case, D.A.; Ben-Shalom, I.Y.; Brozell, S.R.; Cerutti, D.S.; Cheatham, T.E., III; Cruzeiro, V.W.D.; Darden, T.A.; Duke, R.E.; Ghoreishi, D.; Simmerling, C.; et al. Amber 2019; University of California: San Francisco, CA, USA, 2019. [Google Scholar]
  54. Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N·log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 18, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  55. Chow, K.H.; Ferguson, D.M. Isothermal-Isobaric Molecular Dynamics Simulations with Monte Carlo Volume Sampling. Comput. Phys. Commun. 1995, 91, 283–289. [Google Scholar] [CrossRef]
  56. Lzaguirre, J.A.; Catarello, D.P.; Wozniak, J.M.; Skeel, R.D. Langevin Stabilization of Molecular Dynamics. J. Chem. Phys. 2001, 144, 2090. [Google Scholar] [CrossRef]
  57. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general Amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  58. Jakalian, A.; Bush, B.L.; Jack, D.B.; Bayly, C.I. Fast, Efficient Generation of High-Quality Atomic Charges. AM1-BCC Model: I. Method. J. Comput. Chem. 2000, 21, 132–146. [Google Scholar] [CrossRef]
  59. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Automatic atom type and bond type perception in molecular mechanical calculations. J. Mol. Graph. Model. 2006, 25, 247–260. [Google Scholar] [CrossRef] [PubMed]
  60. McGibbon, R.T.; Beauchamp, K.A.; Harrigan, M.P.; Klein, C.; Swails, J.M.; Hernández, C.X.; Schwantes, C.R.; Wang, L.P.; Lane, T.J.; Pande, V.S. MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics Trajectories. Biophys. J. 2015, 109, 1528–1532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-Learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  62. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
  63. Harris, C.R.; Millman, K.J.; vad der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array programming with NumPy. Nature 2020, 585, 357–362. [Google Scholar] [CrossRef] [PubMed]
  64. Lloyd, S.P. Least Squares Quantization in PCM. IEEE Trans. Inf. Theory 1982, 18, 129–137. [Google Scholar] [CrossRef]
  65. Scott, D.W. Multivariate Density Estimation: Theory, Practice, and Visualization, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015. [Google Scholar] [CrossRef]
  66. O’Boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An Open Chemical Toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef] [Green Version]
  67. Rocha, G.B.; Freire, R.O.; Simas, A.M.; Stewart, J.J.P. RM1: A Reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I. J. Comput. Chem. 2006, 27, 1101–1111. [Google Scholar] [CrossRef]
  68. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.; et al. General Atomic and Molecular Electronic Structure System. J. Comput. Chem. 1993, 14, 1347–1363. [Google Scholar] [CrossRef]
  69. Quiroga, R.; Villarreal, M.A. Vinardo: A Scoring Function Based on Autodock Vina Improves Scoring, Docking, and Virtual Screening. PLoS ONE 2016, 11, e0155183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Hou, T.; Wang, J.; Li, Y.; Wang, W. Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations. J. Chem. Inf. Model. 2011, 51, 69–82. [Google Scholar] [CrossRef]
  71. Miller, B.R.; McGee, T.D.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef] [PubMed]
  72. Cock, P.J.A.; Antao, T.; Chang, J.T.; Chapman, B.A.; Cox, C.J.; Dalke, A.; Friedberg, I.; Hamelryck, T.; Kauff, F.; Wilczynski, B.; et al. Biopython: Freely Available Python Tools for Computational Molecular Biology and Bioinformatics. Bioinformatics 2009, 25, 1422–1423. [Google Scholar] [CrossRef]
  73. Rousseeuw, P.J. Silhouettes: A Graphical Aid to the Interpretation and Validation of Cluster Analysis. J. Comput. Appl. Math. 1987, 20, 53–65. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Molecular structure of punicalagin.
Figure 1. Molecular structure of punicalagin.
Biomedicines 09 01533 g001
Figure 2. Protein fluorescence quenching analysis of PDIA1 and PDIA3 in the presence of punicalagin: (a) Fluorescence quenching spectra of reduced PDIA1 and PDIA3 alone (solid line) and after stepwise addition of punicalagin (dotted line) (pH 7.4, 25 °C, and λex = 290 nm, (PDI) = 0.1 × 10−6 M, (punicalagin) = 4 × 10−6 M final concentration); (b) Stern–Volmer plot of quenching data of reduced PDIs in the presence of increasing concentrations of punicalagin. Stern–Volmer constants (KSV) calculated for the different systems tested. Data represent the mean of at least three independent experiments, and error bars indicate SEM.
Figure 2. Protein fluorescence quenching analysis of PDIA1 and PDIA3 in the presence of punicalagin: (a) Fluorescence quenching spectra of reduced PDIA1 and PDIA3 alone (solid line) and after stepwise addition of punicalagin (dotted line) (pH 7.4, 25 °C, and λex = 290 nm, (PDI) = 0.1 × 10−6 M, (punicalagin) = 4 × 10−6 M final concentration); (b) Stern–Volmer plot of quenching data of reduced PDIs in the presence of increasing concentrations of punicalagin. Stern–Volmer constants (KSV) calculated for the different systems tested. Data represent the mean of at least three independent experiments, and error bars indicate SEM.
Biomedicines 09 01533 g002
Figure 3. Thermal denaturation profiles of reduced PDIA1 and PDIA3 obtained by DSC. Proteins, 5 µM concentration, were analyzed in the absence and presence of 25 µM punicalagin. The transition temperature (Tm) of each endothermic peak is reported. The dotted lines show the best fit obtained by deconvolution of each thermogram in non-two-state transitions.
Figure 3. Thermal denaturation profiles of reduced PDIA1 and PDIA3 obtained by DSC. Proteins, 5 µM concentration, were analyzed in the absence and presence of 25 µM punicalagin. The transition temperature (Tm) of each endothermic peak is reported. The dotted lines show the best fit obtained by deconvolution of each thermogram in non-two-state transitions.
Biomedicines 09 01533 g003
Figure 4. Dose-response curve of punicalagin showing the inhibition of PDIA1 and PDIA3 reductase activity using di-E-GSSG as substrate. Each point represents the average of at least three independent measurements. PDIA1, IC50 = 6.1 µM, 95% confidence from 5.5 to 6.7 µM; PDIA3, IC50 = 1.5 µM, 95% confidence from 1.3 to 1.6 µM).
Figure 4. Dose-response curve of punicalagin showing the inhibition of PDIA1 and PDIA3 reductase activity using di-E-GSSG as substrate. Each point represents the average of at least three independent measurements. PDIA1, IC50 = 6.1 µM, 95% confidence from 5.5 to 6.7 µM; PDIA3, IC50 = 1.5 µM, 95% confidence from 1.3 to 1.6 µM).
Biomedicines 09 01533 g004
Figure 5. Scheme of the computational chemistry procedure conducted to provide a putative binding mode for punicalagin on PDIA3 and PDIA1. See text for details.
Figure 5. Scheme of the computational chemistry procedure conducted to provide a putative binding mode for punicalagin on PDIA3 and PDIA1. See text for details.
Biomedicines 09 01533 g005
Figure 6. Selected α and β-punicalagin docking results for PDIA 3 (A) and PDIA1 (B) a’ domain: conformation selected by Vinardo score are depicted in red, while those selected by KDE are blue.
Figure 6. Selected α and β-punicalagin docking results for PDIA 3 (A) and PDIA1 (B) a’ domain: conformation selected by Vinardo score are depicted in red, while those selected by KDE are blue.
Biomedicines 09 01533 g006
Figure 7. Top-ranked β-punicalagin binding modes on (A) PDIA3 a’ domain (−39.2 kcal mol−1) and (B) PDIA1 a’ domain (−35.2 kcal mol−1).
Figure 7. Top-ranked β-punicalagin binding modes on (A) PDIA3 a’ domain (−39.2 kcal mol−1) and (B) PDIA1 a’ domain (−35.2 kcal mol−1).
Biomedicines 09 01533 g007
Figure 8. (A) Solvent accessible surface area (SASA) time series plot and (B) closest distance between a and a’ domains over time plot calculated along MD simulations. (C) PDIA3Ox in “open” and “closed” (D) conformation. PDIA3 a domain is depicted in red, b is yellow, b’ is orange, a’ is blue, WCGHC patterns are green lines while punicalagin is depicted as cyan sticks.
Figure 8. (A) Solvent accessible surface area (SASA) time series plot and (B) closest distance between a and a’ domains over time plot calculated along MD simulations. (C) PDIA3Ox in “open” and “closed” (D) conformation. PDIA3 a domain is depicted in red, b is yellow, b’ is orange, a’ is blue, WCGHC patterns are green lines while punicalagin is depicted as cyan sticks.
Biomedicines 09 01533 g008
Figure 9. Comparative analysis of predicted β-punicalagin binding modes on PDIA3 and PDIA1. Residues are colored by conservation: white indicates conserved regions, blue nonconserved regions, and dark gray not aligned regions. (A) Predicted β-punicalagin binding mode on PDIA3 a’ domain. Nonconserved residue names are reported in black for PDIA3 and red for PDIA1. PDIA1 nonconserved residues are depicted as lines and colored from blue to red according to their score in the alignment matrix. (B) Predicted β-punicalagin binding mode on PDIA1 a’/b’ domains interface. Nonconserved residue names are reported in black for PDIA1 and in red for PDIA3. (C) Predicted β-punicalagin binding modes for PDIA3 (wheat) and PDIA1 (green) reported on the same structure (PDIA1) for comparison. Only a’ and b’ domains are visible.
Figure 9. Comparative analysis of predicted β-punicalagin binding modes on PDIA3 and PDIA1. Residues are colored by conservation: white indicates conserved regions, blue nonconserved regions, and dark gray not aligned regions. (A) Predicted β-punicalagin binding mode on PDIA3 a’ domain. Nonconserved residue names are reported in black for PDIA3 and red for PDIA1. PDIA1 nonconserved residues are depicted as lines and colored from blue to red according to their score in the alignment matrix. (B) Predicted β-punicalagin binding mode on PDIA1 a’/b’ domains interface. Nonconserved residue names are reported in black for PDIA1 and in red for PDIA3. (C) Predicted β-punicalagin binding modes for PDIA3 (wheat) and PDIA1 (green) reported on the same structure (PDIA1) for comparison. Only a’ and b’ domains are visible.
Biomedicines 09 01533 g009
Table 1. Stern–Volmer constants (KSV) and dissociation constant (Kd) obtained by fluorescence quenching analysis of PDIA1 and PDIA3 in the presence of punicalagin. Data were calculated from fluorescence quenching analysis using both proteins (0.1 × 10−6 M) in reducing and non-reducing conditions (pH 7.4, 25 °C) and increasing concentration of punicalagin (up to 4 × 10−6 M). KSV are reported as mean and SEM of at least three independent experiments.
Table 1. Stern–Volmer constants (KSV) and dissociation constant (Kd) obtained by fluorescence quenching analysis of PDIA1 and PDIA3 in the presence of punicalagin. Data were calculated from fluorescence quenching analysis using both proteins (0.1 × 10−6 M) in reducing and non-reducing conditions (pH 7.4, 25 °C) and increasing concentration of punicalagin (up to 4 × 10−6 M). KSV are reported as mean and SEM of at least three independent experiments.
KSV (M−1 × 103)Kd (M)
Reduced form
PDIA197.9 ± 2.111.9 × 10−6
PDIA3157.1 ± 1.910.0 × 10−6
Not reduced form
PDIA1183.2 ± 2.74.9 × 10−6
PDIA3240.4 ± 3.63.9 × 10−6
Table 2. Thermodynamic data of the interaction between punicalagin and PDIAs. Both proteins (25 10−6 M) were analyzed in reduced form. Affinity constant (Kd), molar enthalpy (ΔH), and entropy (TΔS) of the reaction at 25 °C were calculated processing data obtained from the isothermal titration.
Table 2. Thermodynamic data of the interaction between punicalagin and PDIAs. Both proteins (25 10−6 M) were analyzed in reduced form. Affinity constant (Kd), molar enthalpy (ΔH), and entropy (TΔS) of the reaction at 25 °C were calculated processing data obtained from the isothermal titration.
Kd (10−6 M)ΔH (kcal/mol)TΔS (kcal/mol)
PDIA11.0 ± 0.2−2.7 ± 0.45.6 ± 0.5
PDIA31.2 ± 0.3−1.1 ± 0.26.8 ± 0.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paglia, G.; Antonini, L.; Cervoni, L.; Ragno, R.; Sabatino, M.; Minacori, M.; Rubini, E.; Altieri, F. A Comparative Analysis of Punicalagin Interaction with PDIA1 and PDIA3 by Biochemical and Computational Approaches. Biomedicines 2021, 9, 1533. https://doi.org/10.3390/biomedicines9111533

AMA Style

Paglia G, Antonini L, Cervoni L, Ragno R, Sabatino M, Minacori M, Rubini E, Altieri F. A Comparative Analysis of Punicalagin Interaction with PDIA1 and PDIA3 by Biochemical and Computational Approaches. Biomedicines. 2021; 9(11):1533. https://doi.org/10.3390/biomedicines9111533

Chicago/Turabian Style

Paglia, Giuliano, Lorenzo Antonini, Laura Cervoni, Rino Ragno, Manuela Sabatino, Marco Minacori, Elisabetta Rubini, and Fabio Altieri. 2021. "A Comparative Analysis of Punicalagin Interaction with PDIA1 and PDIA3 by Biochemical and Computational Approaches" Biomedicines 9, no. 11: 1533. https://doi.org/10.3390/biomedicines9111533

APA Style

Paglia, G., Antonini, L., Cervoni, L., Ragno, R., Sabatino, M., Minacori, M., Rubini, E., & Altieri, F. (2021). A Comparative Analysis of Punicalagin Interaction with PDIA1 and PDIA3 by Biochemical and Computational Approaches. Biomedicines, 9(11), 1533. https://doi.org/10.3390/biomedicines9111533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop