- Research article
- Open Access
Dynamic changes in the secondary structure of ECE-1 and XCE account for their different substrate specificities
BMC Bioinformaticsvolume 13, Article number: 285 (2012)
X-converting enzyme (XCE) involved in nervous control of respiration, is a member of the M13 family of zinc peptidases, for which no natural substrate has been identified yet. In contrast, it’s well characterized homologue endothelin-converting enzyme-1 (ECE-1) showed broad substrate specificity and acts as endopeptidase as well as dipeptidase. To explore the structural differences between XCE and ECE-1, homology model of XCE was built using the complex structure of ECE-1 with phosphoramidon (pdb-id: 3DWB) as template. Phosphoramidon was docked into the binding site of XCE whereas phosphate oxygen of the inhibitor was used as water molecule to design the apo forms of both enzymes. Molecular dynamics simulation of both enzymes was performed to analyze the dynamic nature of their active site residues in the absence and presence of the inhibitor.
Homology model of XCE explained the role of non-conserved residues of its S2’ subsite. Molecular dynamics (MD) simulations identified the flexible transitions of F149/I150, N566/N571, W714/W719, and R145/R723 residues of ECE-1/XCE for the strong binding of the inhibitor. Secondary structure calculations using DSSP method reveals the folding of R145/R723 residue of ECE-1/XCE into β-sheet structure while unfolding of the S2’ subsite residues in aECE-1 and sustained compact folding of that of aXCE. The results evaluated are in good agreement with available experimental data, thus providing detailed molecular models which can explain the structural and specificities differences between both zinc peptidases.
Secondary structure changes of both enzymes during the simulation time revealed the importance of β-sheet structure of R145/R723 for its binding with the terminal carboxylate group of the inhibitor. Unfolding of the α-helix comprising the S2’ subsite residues in aECE-1 correlate well with its endopeptidase activity while their compact folding in aXCE may account for the inactivity of the enzyme towards large C-terminal containing substrates.
Zinc peptidases such as matrix metalloproteinases (MMPs)[1, 2] angiotensin converting enzyme (ACE)[3, 4] and neutral endopeptidase (NEP) are involved in peptide metabolism. The peptide metabolism is activated by degradation of a wide range of bioactive peptides and therefore specific inhibitors of them have therapeutic values[6, 7]. One of the important classes of M13 family (of zinc peptidases; classification according to MEROPS database) is “gluzincins” which is defined by a HExxH motif including two histidines and a glutamic acid as zinc-coordinating ligands. Zinc peptidases of neprilysin family are gluzincins that include several enzymes for instance, neutral endopeptidase (NEP), NEP2/soluble secreted endopeptidase (SEP)/neprilysin-like enzyme 1 (NL1)/membrane metalloendopeptidase-like 2 (MMEL2), endothelin-converting enzymes ECE-1 and ECE-2[14, 15], the KELL blood group protein, the phosphate-regulating neutral endopeptidase on the X chromosome (PHEX), and X-converting enzyme (XCE)/endothelin-converting enzyme-like 1 (ECEL1) /rodent homologue damaged-induced neuronal endopeptidase (DINE).
XCE (nowadays are known as ECEL1 but we used XCE in this manuscript to differentiate with ECE-1) is expressed in the nervous system, particularly in the medulla oblongata and in the spinal cord, presumably by cholinergic neurons such as motor neurons or striatum interneurons. The physiological function of XCE was first reported from the inactivation of the corresponding gene in mice, which described the enzyme as to play a vital role in the nervous control of respiration. Benoit et al. found that XCE is mostly located in the cellular endoplasmic reticulum (ER) consistent with the high mol. wt smears observed on SDS-PAGE, while only less than 10% portion of the enzyme reaches the cell surface. They also suggested the enzyme’s function in both compartments i.e., ER and the cell surface. ER is considered as a target organelle for XCE to suppress stress due to responding of the enzyme to nerve injury. The presence of XCE at the cell surface refers one of the enzyme’s functions to the regulation of the activity of extracellular peptides like other family members such as ECE-1 and NEP performed. Several candidate neuropeptides such as endothelin, galanin, calcitonin, bradykinin, met- and leu-enkephalins, and somatostatin were considered as potential substrates for XCE activity. Particularly, galanin was supposed to be the most probable one because of its presence and degradation in spinal cord/cerebrospinal fluid. Inhibition of the galanin degradation was reported to be done only by phosphoramidon, but none other inhibitors were found active against XCE. The enzyme was found active only against a synthetic tripeptide substrate, but so far, no obvious reason of the inactivity of XCE against potential neuropeptides has been reported in the literature inspite of they linked its subcellular localization with the lack of success in the identification of extracellular substrate. More work is therefore needed to identify its extracellular substrates to answer questions concerning the enzyme’s precise function.
Contrary to XCE, its well studied homologue ECE-1 can efficiently cleaves variety of substrates. ECE-1 exists as a disulfide-linked homodimer in vivo and cleaves the W21-V22 bond in big endothelin-1 (ET-1), a potent vasoconstrictor. The monomeric C412S mutant of rat ECE-1 (C428S in human) has been shown to have much lower efficiency for the cleavage of big ET-1 as compare to the wild type showing dimerization of ECE-1 which is preferred for effective conversion of big ET-1 into ET-1. Furthermore, M. V. Hoang and A. J. Turner established that ECE-1 also cleaves the unrelated bradykinins (BK) at a significant rate, in addition to its in-vivo substrate big ET-1 (endopeptidase action), thereby acting as a peptidyl dipeptidase. The lack of sequence similarity in the BK peptides and the peptidyl dipeptidase revealed broad specificity and additional physiological roles for ECE-1 possibly linked to its subcellular location. Moreover, recombinant ECE-1 was found to have minimal activity against small substrates (smaller than hexapeptides), such as Leu-enkephalin. However, large peptides such as neurotensin, substance P, bradykinin, and the oxidized insulin B chain were also observed to be hydrolyzed by the enzyme as efficiently as the big ET-1 was. In spite of various natural substrates known upto date, the detailed mechanism of the cleavage of the substrates for such a variable length and unrelated sequences is still lacking for ECE-1. Hence, the understanding of extracellular substrate activities of ECE-1 and XCE remain challenging despite of their several activity profiling experiments because of the lack of appropriate structural knowledge.
Because of the difficulties encountered in the crystallization of proteins specially of the membrane proteins, homology modeling is being served as a valuable tool since last two decades to solve the three-dimensional structures of proteins having at least one X-ray crystal structure of homologous protein[28, 29]. For neprilysin family, initially crystal structure knowledge of thermolysin, a bacterial protein, was used to model the structural features of these proteins. Crystal structure of neutral endopeptidase (pdb-id: 1DMT), the only well characterized member opened new door to acquire structural knowledge of several members of this class of the protein. The three-dimensional structure of ECE-1 was also solved with the co-crystallized metalloprotease inhibitor phosphoramidon (pdb-id: 3DWB). XCE shared much identity with ECE-1 and was also identified as ECE-like protein, therefore the crystal structure coordinates of ECE-1 was utilized to rationalize the structural facts of XCE. Here, the three-dimensional structures of XCE with and without the metalloprotease inhibitor i.e., phosphoramidon were modelled. However, a molecular-level understanding of the function of a biological macromolecule requires knowledge of both its structural and dynamical properties. The dynamical flexibility of both proteins in their inhibitor bound and unbound forms were explored by molecular dynamics (MD) simulations. Our aim is to identify the structural differences between the modeled XCE and ECE-1 by elucidating conformational flexibility of both the enzymes. The characterization of structure and dynamics of both proteins and in particular their subsites with respect to phosphoramidon would be significant to understand the structural differences responsible for their different substrate specificities and to design specific inhibitors against both zinc peptidases.
Results and discussion
In the absence of experimental structures, computational methods are used to predict 3D protein models to provide insight into the structure and function of proteins. There are several successfull examples where homology modeling has aided in the prediction of protein function[32–35]. However, the choice of template, inaccurate alignments and inefficient refinement methods are still the main sources of errors in homology modeling. We therefore critically checked the modeling of XCE at every step from template alignment to refinement before validating the final model from external sources.
Sequence conservation between both proteins
Primary sequence alignment showed that amino acid residues 1–100 in XCE could not be modeled due to lack of equivalent residues in ECE-1 structure. Therefore, R102 is considered as the first amino acid of the ectodomain of XCE model which is aligned with S101 of ECE-1 as depicted from Figure1. Final alignment comprised of 4 gaps in XCE sequence, and these gaps are not longer than 3 amino acids. Three characteristic motifs, 565VNAYY569, 607HELTH611, and 667ENIAD671 of ECE-1, hallmarks of the M13 family of zinc peptidases, are conserved and correspond to 570LNAYY574, 612HELTH616, and 672ENIAD676 motifs in XCE respectively. The presence of VNAYY motif in ECE-1 is responsible for the specificity towards big endothelins binding whereas semi-conserved LNAYY motif in XCE may have influence on the binding of substrate. Two histidines H612 and H616 in the HELTH motif and one glutamate E672 in the ENIAD motif coordinate to the catalytic zinc (Zn2+) in XCE. Corresponding to E608 of ECE-1, E613 of XCE is believed to act as a nucleophile which promotes the attack of zinc-bound water molecule on the scissile peptide. The resulting transition state formed during hydrolysis of substrate gets stabilized by H737 in XCE corresponding to H732 of ECE-1. D614 and D671 of ECE-1 form Asp-His-Zinc triad to place zinc atom appropriately for catalysis and identical role can also be predicted for the conserved D619 and D676 in XCE. Finally, the evaluation of stereochemical properties is presented by the Ramachandran plot provided as Additional file1: Figure S1.
Comparison between protein-inhibitor complexes
The docked conformation of phosphoramidon in the modelled XCE is compared in terms of its interaction with the X-ray bound conformation of ECE-1. Figure2A represents the docked conformation of phosphoramidon in complex with XCE whereas the X-ray bound conformation of phosphoramidon with ECE-1 is shown in Figure2B. The plausibility of the enzyme inhibitor complex in zinc peptidases is primarily strengthened by the backbone hydrogen bond contributors from “NAYY” motif. The carbonyl oxygen of A572 in XCE created a weak hydrogen bond with the P1’ amide nitrogen of phosphoramidon in a similar manner as formed by A567 in ECE-1. Two hydrogen bonds were found for N566 of ECE-1 with the P2’ amide nitrogen and oxygen atoms of the terminal carboxylate while its equivalent N571 of XCE maintained only one H-bond with the P2’ amide nitrogen. In XCE, however, V565 of ECE-1 was replaced with L570 amino acid, but hydrogen bonds between the P2’ indole moiety of phosphoramidon and their carbonyl oxygen atoms are conserved in both proteins. R743 of XCE was also involved in the hydrogen bond with the P1’ carbonyl oxygen of the inhibitor as observed for R738 of ECE-1. Hence, all hydrogen bonding interactions are also found in the XCE-bound docked conformation of the inhibitor.
The specificity of all metalloproteases depends on the interactions of the P1’ and P2’ moieties of their substrates/inhibitors with the corresponding S1’ and S2’ subsites of the enzyme. We used the well characterized subsite specificity information of ECE-1 to understand the structural difference in XCE present at the subsite level. The amino acid residues I587, L592, I606, T608, I609, F715, N718 and W719 constitute the S1’ subsite in XCE corresponding to the large pocket of ECE-1 comprising I582, F587, I601, V603, V604, F710, V713 and W714 residues as illustrated from Figures2C and2D, respectively. The P1’ group with large hydrophobic side-chains indicates the specificity of ECE-1 due to the large hydrophobic S1’ subsite. The presence of non-conserved hydrophilic T608 in XCE corresponding to V603 of ECE-1 did not change the overall hydrophobic environment of the S1’ recognition pocket because of the formation of strong hydrogen bond with G578 located outside the binding pocket. Replacement of V604 with I609 in XCE provided more hydrophobic environment than that of ECE-1, thus resulting in a favorable interaction with the P1’ hydrophobic side-chain of the inhibitor. N718 in XCE occupied more space which increases the depth and flexibility of the S1’ subsite in XCE as compared to ECE-1 with small V713.
The S2’ subsites of XCE and ECE-1 contain all non-conserved amino acids except one arginine as shown in Figures2E and2F, respectively. The side-chain of R145, S150, and W153 in ECE-1 occupied more space and volume than their counterpart T146, A151, and G154 residues in XCE. Schulz et al. suggested that the sidechains of the S2’ subsite residues in ECE-1 were disordered and did not establish any hydrogen bond due to which, indole moiety of phosphoramidon was not appropriately anchored into the S2’ pocket of ECE-1. On the other hand, all residues of the S2’ subsite in XCE are found very close, thus leading to very narrow binding pocket responsible for strong binding of phosphoramidon at this region. R145 of ECE-1 formed water mediated hydrogen bonds with the terminal carboxylate of phosphoramidon, which are lacking in XCE due to the non-conserved T146 residue. Table1 summarizes the effect of non-conserved subsite residues of both enzymes and their possible consequences on the inhibitor binding. The binding pattern of phosphoramidon and its relative stability into the active site of both proteins are further explored with the help of MD simulations.
Molecular dynamics simulations
Molecular dynamics simulations permit characterization of biomolecular processes such as the conformational transitions associated with the functions of protein at molecular and atomic level. Nowadays, this technique is being extensively used to study the conformational flexibility of biological macromolecules because of the advances in computer hardware, software, and MD algorithms. We calculated fluctuations in protein structure for all four systems aXCE, cXCE, aECE-1, cECE-1 by means of nanosecond timescale MD simulations which are proved to be a complementary technique for experiments. Modern computer architectures provide means to perform MD simulations of varying time length that can range from nanoseconds to microsecond-timescale or even at millisecond timescale. However, long timescale MD simulations based on all-atom force fields are still not so convincing to be employed. Since most of the force fields were parameterized previously and are therefore not extensively applied in this time regime. As an illustration of the problems related to unrealistic and irreversible structural changes, the DNA simulation of the timescale beyond 20–30 ns time could be exemplified. Such long timescale MD simulations of biological macromolecules utilizing enhanced sampling methods and/or coarse-graining methods were frequently reported. On the other hand short timescale MD simulations such as of few nanoseconds have serious concerns about dealing the process of relaxation and the proper sampling of the equilibrium dynamics. These two significant factors must be kept in mind when presenting results based on the short MD run. Therefore, to ensure about sufficient amount of sampling which is a hard task, the use of replicate simulations could be helpful to address this issue. In the absence of replicate simulations the situation becomes worse to test proper sampling as well as the reproducibility of the results obtained from the simulations could be challenged. Could be an alternative of the replicate simulations, the use of better analysis of the equilibration process would be significant. Plotting root mean square deviation (RMSD) from the time-averaged structure (rmsdtime−averaged) calculated from the simulation trajectories of all systems would also adequately answer other relevant questions about sampling time as opposed to rmsdinit.
Root mean square deviation (RMSD) of the protein backbone
An initial evaluation involves analysis of the protein structure via root mean square deviations (RMSD) of backbone atoms with respect to the equilibrated structure (rmsdinit) to perform sampling of trajectories and to evaluate the stability of proteins throughout the simulation time as illustrated in Figure3A. The rmsdinit for apo XCE (aXCE) fluctuated between 1.0 and 1.8 Å with an average value of 1.45 Å whereas rmsdinit of the inhibitor bound XCE (cXCE) appeared from the same starting point as in aXCE but after 1 ns it exhibited a large deviation up to 2.5 Å, which indicates changes in the protein structure due to protein-inhibitor complex formation. The template protein ECE-1 in their apo (aECE-1) and complexed form (cECE-1) showed an increased rmsdinit value due to more flexibility compared to its homologue XCE in their corresponding forms (aXCE and cXCE). The complex formation in ECE (cECE-1) resulted in rmsdinit difference of ~0.48 Å from the average rmsdinit value calculated for its apo form (~1.72 Å). The difference in rmsdinit values for aXCE and cXCE resembles to ECE-1 proteins which indicates the similar conformational changes in both proteins from the binding of phosphoramidon which will be described in latter sections. Figure3B illustrates the root mean square deviation with respect to the X-ray structure of all heavy atoms of ECE-1 which exhibited a fluctuation between ~2.5 and ~3.0 Å with an average value of 2.6 Å. The identical approach was used to analyze the structural deviation in terms of rmsd with respect to the initial homology model of XCE (Figure3B), which varied between ~2.0 and ~3.0 Å, thus showing a significant structural deviation from the homologue structure. Besides these two rmsd plots such as rmsdinit and rmsd with respect to X-ray structure which are considered poor metrics to observe equilibration process, the rmsdtime−averaged plots of all proteins forming V-shape could therefore be a fair substitute to observe the under-sampled simulations. The smallest rmsd values in the plots shown in Figure3C are concentrated around the center of the simulations of 10 ns showing that the equilibration is in progress. U-shaped rmsdtime−averaged plots were obtained but not V-shaped that show that the equilibration is in progress along with some contaminations of relaxation processes. Nevertheless, the behavior shown by the simulations having both equilibration and a small amount relaxation processes can still be regarded as providing biologically significant information. A question can also arisen that why the short MD run could provide a conformation close to that of the native state of the protein as determined via X-ray structure (pdb-id : 3DWB). This issue was addressed by Ramachandran plots (Additional file2: Figure S2, Additional file3: Figure S3, Additional file4: Figure S4, Additional file5: Figure S5) developed for all the equilibrated structures that are in excellent agreement with experimental data. This validates the simulation approach which was used to yield the nearly native conformation of the protein as all simulations adopted the near crystal structure distributions with some variation in the secondary structure elements. The approach is therefore, applicable to analyze other properties of the given size of protein.
B-factor and secondary structure calculations
The fluctuations in the local structure of XCE and ECE-1 protein in their apo and complexed forms were further investigated employing the B-factors (Bi) calculation for αC atom of each residue. Figure4A and4B depict calculated B-factors of (apo and complexed) XCE and ECE-1, respectively reflecting significant fluctuations in the loop regions. The calculated average B-factor of cECE-1 (24.3 Å2) was lower than average crystallographic B-factors for αC atom (49.3 Å2) resulting from a higher degree of hydration in the simulation compared to the crystal structure. Considering the information obtained from comparison between theoretical and experimental B-factors of cECE-1 protein, the similar approach was used to calculate average B-factors of each XCE amino acids. The average cXCE B-factors was calculated as 24.18 Å2 which is comparable to the average B-factor of cECE-1 indicating a bit lower dynamic flexibility of cXCE compared to cECE-1. Secondary structure calculations were also performed to correlate fluctuations of residues with the corresponding changes in the secondary structure elements for aXCE and cXCE, as shown in Figures5A and5B, respectively. Amino acid residues 137–147 of aXCE located in the region L1 comprised mostly of loops along with a turn which underwent antiparallel β sheet conformation due to complexation with the inhibitor (cXCE). In aXCE, residues 145–146 were folded into β strand after a period of 4 ns whereas β strand of 143–145 residues in cXCE was degraded into loops during the simulation observed via visual inspection. Furthermore, amino acid residues 151–154 of the S2’ subsite passed from α-helix to 3–10 helix, which in turn affects the region L1 of XCE because it is found at the immediate vicinity of the S2’ subsite. Loop area L1 in cXCE was more fluctuated than that of aXCE, which indicates conformational modifications of this region after the inhibitor binding (Figure4A). The region L2 (166–173) of both XCE comprised of a loop (166–170) and an α-helix (171–173), respectively. Secondary structure plots of XCE proteins (Figures5A and5B) suggest that the unfolding of the α-helix (171–173) results in an increased fluctuation of the region L2 in aXCE which is reduced upon binding with the inhibitor. Moreover, regions L3 (residues 190–194) and L4 (residues 221–224) are situated at far bottom of the active site and exhibited high level of fluctuations in the secondary structure from α-helix to 3-helix/turn and from 3–10 helix to loop/turn, respectively (Figures5A and5B). At the front of the binding site, most significant regions L5 (residues 330–350), L6 (residues 368–383) and L7 (residues 410–445) are comprised of loops only and therefore, are expected to show high flexibility. These regions exhibited more fluctuations in cXCE compared to aXCE indicating the direct influence of the inhibitor binding (Figure4A). Major deviations in the secondary structure of the residues 330–338 of L5, 378–383 of L6, and 424–426 of L7 are observed as they are folded into α-helix in aXCE whereas they showed random transitions from 3–10 helix/turns to unstructured loops in cXCE (Figures5A and5B). Figure4B illustrates B-factors of each amino acid residue of apo and complexed ECE-1 correlating the extent of structural modifications produced by the inhibitor. During simulation, the region L1 (residues 143–155) including residues of the S2’ subsite experiences changes in its secondary structure elements for both forms of ECE-1 as shown in Figures5C and5D. Amino acid residues 143–144 are present as turn in aECE-1 whereas as parallel β-sheet in cECE-1. R145 belonging to the S2’ subsite was folded into antiparallel β-sheet at initial stage of the simulation in both forms of ECE-1. It is suggested that the folding of R145 is necessary for its strong interaction with the terminal carboxylate. The conformation adopted by amino acid residues 150–155 in aECE-1 was a compact α-helix but during simulation this α-helix was transformed into a turn due to its drifting from the active site (visual inspection), however, it is found as a compact α-helix throughout the simulation in cECE-1 (Figures5C and5D). Amino acid residues 274–275 close to the active site constitute a loop (L2) which showed large fluctuations in both forms of ECE-1. In cECE-1, L3 (residues 296–297) and L4 (residues 421–428) get stabilized after the inhibitor binding attributed by their low fluctuations in B-factors while L5 (residues 646–647) and L6 (residues 681–691) showed higher fluctuations compared to those of aECE-1 (Figure4B). Secondary structure calculations revealed conformational transformations into turns/β sheet/α-helix in the L4 region only whereas L3, L5, and L6 regions remain unstructured as turns/loops in both forms of ECE-1 (Figures5C and5D). Thus, inhibitor binding accounts for high fluctuations of the loop regions and affects overall secondary structure elements surrounding the active site of both XCE and ECE-1 proteins.
Conformational changes in the active site amino acids of XCE
The conformational changes induced by the inhibitor on the active site of zinc peptidases have not been understood yet in detail due to the nonavailability of their 3D apo structures. The separate sampling of the bound and unbound states is the key to describe conformational differences upon binding. Figure6A shows the superimposed aXCE structures before and averaged over 10 ns simulation whereas similar coordinates for cXCE are displayed in Figure6B. Significant conformational changes in the side-chains of I150, L570, N571 and R723 are observed in both forms as a result of simulation. During simulation, the side-chain of I150 acquired cis conformation in cXCE, being close to the P1’ leucyl moiety of the inhibitor, whereas its trans conformation is stable in aXCE. N571 rotated its side-chain towards the P2’ carbonyl oxygen in cXCE during simulation. The inside bending of N571 is found initially in both forms of XCE while its outside bending is only observed during the simulation in cXCE which results in the formation of hydrogen bond (Figure6B). Inside and outside transition of N571 drastically affects the hydrogen bonding interaction of L570 with the P2’ indole moiety of the inhibitor. All residues of the S1’ and S2’ subsites are stable in their initial position, however, W719 approached more close to the P1’ leucyl moiety of the inhibitor during simulation. The positively charged R723 in cXCE acquired extended conformation and attracted towards the negatively charged carboxylate of the inhibitor as the simulation proceed. DSSP calculations also revealed the folding of R723 into antiparallel β-sheet during simulation.
Conformational changes in the active site amino acids of ECE-1
Figure7A depicts the active site residues of the superimposed aECE-1 structures before and averaged over 10 ns simulation while similar coordinates for cECE-1 are shown in Figure7B. In cECE-1, the side-chain of R145 approached very close to the carboxylate group of the inhibitor during simulation, however, it was found far apart before simulation. DSSP calculations also illustrated the folding of R145 into antiparallel β-sheet during simulation. Side-chain nitrogen atom of N566 rotated towards the P1’ amidic NH group of the inhibitor in cECE-1 during simulation causes loss of its interaction which was found in the crystal structure (Figure7B). Before simulation, F149, F587, and W714 are perpendicularly inclined towards the P1’ leucyl moiety of the inhibitor and known to form an aromatic wall which can block its deeper approach into the S1’ pocket. During simulation, W714 folded into 3–10 helix and approached more closer to the P1’ group of the inhibitor in cECE-1 as compared to the crystal structure. All residues of the S1’ subsite are stable in their initial position but I582 changes trans conformation in cECE-1. In aECE-1, remarkable conformational drift in the S2’ residues (149–153) from the active site is observed (Figure7A), however, in cECE-1 this drifting is not found due to the presence of the inhibitor which can hold them towards itself. It is not identified yet that in order to fit into the interior of the enzyme whether big endothelins, the in-vivo substrates of ECE-1, forces an opening of the enzyme wall or requires a folded conformation of their C-terminal part. Hence, the conformational drift in the S2’ residues identified in our simulation study correlates well with the activity of aECE-1 towards the large C-terminals containing natural substrates.
Stability of enzyme-inhibitor interactions
Both enzymes contain tetrahedral zinc ion (Zn2+) coordinated with H607/H612, H611/H616, E667/E672 and the PO1 atom of phosphoramidon in ECE-1/XCE, respectively. Table2 summarizes the average coordinate bond distances whereas average hydrogen bond and hydrophobic interactions evaluated over stable period of 2–8 ns for the enzyme-inhibitor complexes are given in Table3. The stability of the hydrogen bonds between the inhibitor and protein are probed by their percent occupancy. Backbone hydrogen bonding residues of XCE are L570, N571, and A572, which are equivalent to V565, N566, and A567 of ECE-1. The P1’ amidic nitrogen, the P2’ amidic nitrogen and the carbonyl oxygen of phosphoramidon strongly bound to the carbonyl oxygen of A572, the carbonyl oxygen and amidic nitrogen of N571 in XCE, also attributed by their high percent occupancy, whereas the sustainability of these interactions from A567 and N566 of ECE-1 is low (Table3). During simulation, non-conserved residues L570 of XCE and V565 of ECE-1 got changes in their backbone dihedral angles which led to the loss of their interaction with the NH group of the indole moiety of phosphoramidon. Thus, it was analyzed that not even a single backbone residue of ECE-1 maintain its vital hydrogen bonding, which is also corroborated by the inhibitory potency of phosphoramidon in μM range against ECE-1. An hydrogen bond from R743 of XCE with the P1’ carbonyl oxygen showed very low occupancy compared to the corresponding bond from R738 of ECE-1 which comparatively attained higher stability (Table3). I582 and V604 of ECE-1 remained close to the P1’ group during simulation, however, equivalent residues I587 and I609 of XCE were found far apart before simulation and only I609 established hydrophobic interaction during simulation (Table3). F149 of ECE-1 formed van der Waal contact and π-π stacking to the P2’ indole moiety while I150 of XCE did not interact with the P2’ group but prefered hydrophobic interaction with the P1’ moiety due to its cis conformation (Table3). Before simulation, R145 of ECE-1 was involved in water mediated hydrogen bonds with the terminal carboxylate but after 1500 ps these interactions are converted into strong salt bridge. Johnson et al. determined carboxydipeptidase activity of ECE-1 during the cleavage of bradykinin and substance P revealing that these substrates with COOH-terminal carboxylates are highly preferred over their cognate amides and esters. Interestingly, formation of salt bridge from R145 coorelate with the carboxydipeptidase activity of ECE-1, however, Johnson et al. did not find any significant decrease in activity when measured against R133L mutant enzyme. The responsibility of the terminal carboxylate binding in XCE is taken by R723 (equivalent to R718 of ECE-1) and in doing so, N571 and R723 of XCE behave like the terminal carboxylate binding residues N542 and R102 of neutral endopeptidase (NEP). During simulation, the binding mode of phosphoramidon in XCE is found similar to its binding in NEP and it can be speculated that phosphoramidon binds firmly with XCE just like its strong binding to NEP. The different binding conformation of phosphoramidon in XCE and ECE-1 strengthens our hypothesis that although XCE is more close to ECE-1 but still minor differences in the behavior of their equivalent residues of the active site and its surrounding affect their tendency towards natural substrates.
The great interest in members of the M13 family of zinc peptidases as putative therapeutic targets and the discovery of new members of this family lead us to investigate the three-dimensional structure and dynamic flexibility of XCE (a novel member of the M13 family) and enable us to compare it with already well established member ECE-1. Homology modeling of XCE and docking of phosphoramidon into the modeled structure provided the primary information about its structural similarities and differences with ECE-1. Simulation of apo and phosphoramidon bound forms of XCE and ECE-1 revealed that the inhibitor induced distinct conformational changes surrounding the active site of both enzymes. The stability of I150, I609 and W719 of XCE and I582, V604 and W714 of ECE-1 suggest that these residues behave as the central determinant in the hydrophobic lock of their S1’ pockets. W714/W719 maintains 3–10 helical structure whereas R145/R723 and F149/T146 of ECE-1/XCE acquire antiparallel β-sheet structure for their strong binding to the inhibitor. Secondary structure changes of the S2’ residues (150–155) from α-helix to turn in aECE-1 explain the accomodation of large C-terminals of natural substrates into this region (endopeptidase activity) while their compactness in helical form clarify the difference in specificity of aXCE at this region. The carboxydipeptidase activity of ECE-1 is also justified by the folding of R145 into antiparallel β-sheet and similar role can be speculated for its structural counterpart R723 of XCE. Hence, larger variations and greater fluctuations in the S2’ subsite of both enzymes would be instrumental in determining their different substrate specificities. To the best of our knowledge, this is the first theoretical study of ECE-1 and XCE in context of their comparison. The structural features of both proteins and their dynamical behavior with phosphoramidon obtained from this study provide valuable information that can be used in the rational design of specific inhibitors for ECE-like proteins. However, molecular mechanism of substrate specificity in XCE needs to be ascertained carefully with the help of detailed experiments in the light of our preliminary theoretical study.
The human XCE sequence (accession number O95672-1) obtained from SWISS-PROT was searched for a suitable template using Blast at NCBI database. As expected, the human ECE-1 phosphoramidon complex was identified as the most befitting template with E-value of 2 × 10−153, thus exhibiting 39% identity and 60% similarity with the query sequence. The sequence of the ectodomain of human ECE-1 (90–770) was taken directly from Protein Data Bank and aligned to XCE sequence to compare the main features of both proteins by utilizing the program Jalview. Homology modeling was carried out by utilizing the homology module implemented in the molecular operating environment (MOE) suite of programs. MOE align facility was used to give the best alignment as input for modeling. Initially, 10 models were built and evaluated within the MOE package by a residue packing quality function, which depends on the number of buried nonpolar side-chain groups and on hydrogen bonding. Energy minimizations were carried out to remove geometrical strain using Amber99 potential with RMS gradient of 0.01. The zinc ion and water molecule bound to zinc were positioned in a similar manner as observed in the ECE-1 structure by superimposing the homology model and template structure. The best intermediate model with the lowest energy was chosen for the quality assessment using Procheck.
Molecular dynamics simulations
Starting structures for MD simulations
The apo form of both proteins denoted as aECE-1 and aXCE were constructed by coordinating the zinc ion to water molecule ‘Wat762’ and ‘Wat776’, respectively by using the zinc bound oxygen coordinates of phosphoramidon present in the pdb-id 3DWB (the complex of ECE-1 with the phosphoramidon inhibitor that was termed as cECE-1). The bioactive conformation of phosphoramidon was manually transferred into the modeled XCE to form the complexed form of XCE named as cXCE. Energy minimizations were carried out on the phosphoramidon and its interacting protein amino acids in cXCE to remove geometrical strain using Amber99 potential with the restraints applied on rest of the protein. To establish potentials for the complexed forms, electrostatic potential (ESP) for all atoms of phosphoramidon were calculated at HF level of theory using 6-31G* basis sets with Gaussian03. The ESP charges were then fitted to obtain restrained electrostatic potential (RESP) charges with antechamber module of AMBER 10. It had been reported that the interaction energies between zinc and its coordinates were markedly underestimated by molecular mechanics calculations if the zinc ion was modeled octahedral traditionally. To minimize the underestimation that would hamper the stability and interactions of zinc coordinates, the Cationic Dummy Atom (CaDA) approach was used. The approach contains four identical dummy atoms tetrahedrically attached to the zinc ion and transfers all the atomic charge of the zinc divalent cation evenly to the four dummy atoms. The atoms are dummy in that sense they interact with other atoms electrostatically but not sterically, thus mimicking 4sp3 vacant orbitals of zinc atom accommodating its lone-pair electrons. The force field parameters of the tetrahedron-shaped zinc were used from the published protocol in molecular dynamics (MD) simulations[51, 52]. In case of ECE-1, the first-shell zinc residues H607, H611, E667 and Wat762 while that of XCE, H612, H616, E672 and Wat776 were deprotonated as histidinate, glutamate, and hydroxide, respectively. In the second hydration shell of zinc ion, amino acids E608, D614, D671, H732 and E613, D619, D676, H737 were kept protonated in ECE-1 and XCE, respectively[53–56]. These structures were then solvated with cubic box of TIP3P water model extended to 8 Å beyond the surface of protein and neutralized with sodium ions. After preliminary modifications of all four systems, i.e. aXCE, cXCE, aECE-1, and cECE-1, each system was subjected to 10 ns molecular dynamics (MD) simulation using SANDER module of AMBER installed on a cluster computing facility at University of Karachi consisting 10 nodes.
Molecular dynamics simulation protocol
Force field parameters such as Duan et al. and the generalized amber force field (gaff) were utilized for protein and phosphoramidon, respectively. System topologies and coordinate files were generated with the xleap module of AMBER. Prior to MD simulations, stepwise energy minimizations were performed, which consist of 1000 steps on phosphoramidon with a positional constraint of 10 kcal mol−1 applied to rest of the complex. 50 steps energy minimizations on the tetrahedron shaped zinc divalent cation and another 50 steps energy minimizations on the zinc alongwith its coordinating residues and phosphoramidon were performed with a positional constraint applied to rest of the system followed by 2000 steps of unrestrained energy minimizations. Each system was gradually heated in the NPT ensemble from 0 to 300 K over 50 ps. All simulations used a dielectric constant of 1.0, a periodic boundary condition at a constant temperature of 300 K, a constant pressure of 1 atm with isotropic molecule based scaling and the SHAKE algorithm applied to constrain all bonds involving the hydrogen atom. The temperature was rescaled with the Berendsen coupling algorithm and the particle mesh ewald method[61, 62] was used to calculate long range electrostatic interactions using a grid size of 80 × 96 × 90 (grid spacing of 1.0 Angstrom) combined with a fourth-order B-spline interpolation to compute the potential and forces in between grid points. Each system was equilibrated by gradually reducing the constraints and finally 10 ns simulation was performed with 2.0 fs time step. The resulting trajectories were analyzed using the PTRAJ included in AmberTools 1.4 package. Hydrogen bonds were defined to have a maximum hydrogen-acceptor distance of 3.5 Å and a minimum donor-hydrogen-acceptor angle of 120°. Hydrogen bonds possessing greater than 40% occupancy are classified as strong, 20–40% are weak, and less than 20% are termed as unstable over 10 ns simulation. Chimera and VMD were used to visualize the structures and trajectories obtained from MD simulation. Secondary structure analyses over 10 ns simulation were carried out employing the defined secondary structure of proteins (DSSP) method.
Visse R, Nagase H: Matrix metalloproteinases and tissue inhibitors of metalloproteinases: structure, function, and biochemistry. Circ Res 2003, 92: 827–839. 10.1161/01.RES.0000070112.80711.3D
Takaishi H, Kimura T, Dalal S, Okada Y, D’Armiento J: Joint diseases and matrix metalloproteinases: a role for MMP-13. Curr Pharm Biotechnol 2008, 9: 47–54. 10.2174/138920108783497659
Matchar DB, McCrory DC, Orlando LA, Patel MR, Patel UD, Patwardhan MB, Powers B, Samsa GP, Gray RN: Systematic review: comparative effectiveness of angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers for treating essential hypertension. Ann Intern Med 2008, 148: 16–29.
Doulton TWR, He FJ, MacGregor GA: Systematic review of combined angiotensin-converting enzyme inhibition and angiotensin receptor blockade in hypertension. Hypertens 2005, 45: 880–886. 10.1161/01.HYP.0000161880.59963.da
Turner AJ, Isaac RE, Coates D: The neprilysin (NEP) family of zinc metalloendopeptidases: genomics and function. BioEssays 2001, 23: 261–269. 10.1002/1521-1878(200103)23:3<261::AID-BIES1036>3.0.CO;2-K
Molinaro G, Rouleau JL, Adam A: Vasopeptidase inhibitors: a new class of dual zinc metallopeptidase inhibitors for cardiorenal therapeutics. Curr Opin Pharmacol 2002, 2: 131–141. 10.1016/S1471-4892(02)00138-8
Whitworth JA: Emerging drugs in the management of hypertension. Expert Opin Emerg Drugs 2003, 8: 377–388. 10.1517/1472822.214.171.1247
Rawlings ND, Morton FR, Barrett AJ: MEROPS: the peptidase database. Nucleic Acid Res 2006, 34: D270-D272. 10.1093/nar/gkj089
Roques BP, Noble F, Dauge V, Fournie-Zaluski MC, Beaumont A: Neutral endopeptidase 24.11: structure, inhibition, and experimental and clinical pharmacology. Pharmacol Rev 1993, 45: 87–146.
Ouimet T, Facchinetti P, Rose C, Bonhomme MC, Gros C, Schwartz JC: Neprilysin II: A putative novel metalloprotease and its isoforms in CNS and testis. Biochem Biophys Res Commun 2000, 271: 565–570. 10.1006/bbrc.2000.2664
Ikeda K, Emoto N, Raharjo SB, Nurhantari Y, Saiki K, Yokoyama M, Matsuo M: Molecular identification and characterization of novel membrane-bound metalloprotease, the soluble secreted form of which hydrolyzes a variety of vasoactive peptides. J Biol Chem 1999, 274: 32469–32477. 10.1074/jbc.274.45.32469
Ghaddar G, Ruchon AF, Carpentier M, Marcinkiewicz M, Seidah NG, Crine P, Desgroseillers L, Boileau G: Molecular cloning and biochemical characterization of a new mouse testis soluble-zinc-metallopeptidase of the neprilysin family. Biochem J 2000, 347: 419–429. 10.1042/0264-6021:3470419
Bonvouloir N, Lemieux N, Crine P, Boileau G, DesGroseillers L: Molecular cloning, tissue distribution, and chromosomal localization of MMEL2, a gene coding for a novel human member of the neutral endopeptidase-24.11 family. DNA Cell Biol 2001, 20: 493–498. 10.1089/104454901316976127
Xu D, Emoto N, Giaid A, Slaughter C, Kaw S, DeWit D, Yanagisawa M: ECE-1: a membrane-bound metalloprotease that catalyzes the proteolytic activation of big endothelin-1. Cell 1994, 78: 473–485. 10.1016/0092-8674(94)90425-1
Emoto N, Yanagisawa M: Endothelin-converting enzyme-2 is a membrane bound, phosphoramidon-sensitive metalloprotease with acidic pH optimum. J Biol Chem 1995, 70: 15262–15268.
Lee S, Zambas ED, Marsh WL, Redman CM: Molecular cloning and primary structure of Kell blood group protein. Proc Natl Acad Sci USA 1991, 88: 6353–6357. 10.1073/pnas.88.14.6353
Francis F, Hennig S, Korn B, Reinhardt R, Jong PD, Poustka A, Lehrach H, Rowe PSN, Goulding JN, Summerfield T: A gene (PEX) with homologies to endopeptidases is mutated in patients with X–linked hypophosphatemic rickets. Nat Genet 1995, 11: 130–136. 10.1038/ng1095-130
Valdenaire O, Richards JG, Faull RLM, Schweizer A: XCE, a new member of the endothelin-converting enzyme and neutral endopeptidase family, is preferentially expressed in the CNS. Mol Brain Res 1999, 64: 211–221. 10.1016/S0169-328X(98)00321-0
Valdenaire O, Rohrbacher E, Langeveld A, Schweizer A, Meijers C: Organization and chromosomal localization of the human ECEL1 (XCE) gene encoding a zinc metallopeptidase involved in the nervous control of respiration. Biochem J 2000, 346: 611–616. 10.1042/0264-6021:3460611
Kiryu-Seo S, Sasaki M, Yokohama H, Nakagomi S, Hirayama T, Aoki S, Wada K, Kiyama H: Damage-induced neuronal endopeptidase (DINE) is a unique metallopeptidase expressed in response to neuronal damage and activates superoxide scavengers. Proc Natl Acad Sci USA 2000, 97: 4345–4350. 10.1073/pnas.070509897
Schweizer A, Valdenaire O, Koster A, Lang Y, Schmitt G, Lenz B, Bluethmann H, Rohrer J: Neonatal lethality in mice deficient in XCE, a novel member of the endothelin-converting enzyme and neutral endopeptidase family. J Biol Chem 1999, 274: 20450–20456. 10.1074/jbc.274.29.20450
Benoit A, Vargas MA, Desgroseillers L, Boileau G: Endothelin-converting enzyme-like 1 (ECEL1) is present both in the plasma membrane and in the endoplasmic reticulum. Biochem J 2004, 380: 8812–8818.
Kiryu-Seo S, Kiyama H: DINE (damage-induced neuronal endopeptidase). Protein Pept Lett 2004, 11: 451–460. 10.2174/0929866043406526
Bedecs K, Langel U, Bartfai T: Metabolism of galanin and galanin (1–16) in isolated cerebrospinal fluid and spinal cord membranes from rat. Neuropeptides 1995, 29: 137–143. 10.1016/0143-4179(95)90015-2
Shimada K, Takahashi M, Turner AJ, Tanzawa K: Rat endothelin-converting enzyme-1 forms a dimer through Cys412 with a similar catalytic mechanism and a distinct substrate binding mechanism compared with neutral endopeptidase. Biochem J 1996, 315: 863–867.
Hoang MV, Turner AJ: Novel activity of endothelin-converting enzyme: hydrolysis of bradykinin. Biochem J 1997, 327: 23–26.
Johnson GD, Stevenson T, Ahn K: Hydrolysis of peptide hormones by endothelin-converting enzyme-1. J Biol Chem 1999, 274: 4053–4058. 10.1074/jbc.274.7.4053
Wiltgen M, Tilz PG: Homology modelling: a review about the method on hand of the diabetic antigen GAD 65 structure prediction. WMW Wiener Medizinische Wochenschrift 2009, 159: 112–125. 10.1007/s10354-009-0662-z
Cavasotto CN, Phatak SS: Homology modeling in drug discovery: current trends and applications. Drug Discov Today 2009, 14: 676–683. 10.1016/j.drudis.2009.04.006
Oefner C, D’Arcy A, Hennig M, Winkler FK, Dale GE: Structure of human neutral endopeptidase (neprilysin) complexed with phosphoramidon. J Mol Biol 2000, 296: 341–349. 10.1006/jmbi.1999.3492
Schulz H, Dale GE, Karimi-Nejad Y, Oefner C: Structure of human endothelin-converting enzyme I complexed with phosphoramidon. J Mol Biol 2009, 385: 178–187. 10.1016/j.jmb.2008.10.052
Song L, Kalyanaraman C, Fedorov AA, Fedorov E, Glasner M, Brown S, Imker H, Babbitt P, Almo S, Jacobson M: Prediction and assignment of function for a divergent N-succinyl amino acid racemase. Nat Chem Biol 2007, 3: 486–491. 10.1038/nchembio.2007.11
Guimaraes AJ, Hamilton AJ, Guedes HLM, Nosanchuk DJD, Zancope-Oliveira R: Biological function and molecular mapping of M antigen in yeast phase of Histoplasma capsulatum. PLoS One 2008, 3: e34491-e344917.
Sun W, Gerth C, Maeda A, Lodowski DT, Kraak LVD, Saperstein DA, Heon E, Palczewski K: Novel RDH12 mutations associated with Leber congenital amaurosis and cone-rod dystrophy: biochemical and clinical evaluations. Vision Res 2007, 47: 2055–2066. 10.1016/j.visres.2007.04.005
Proell M, Riedl SJ, Fritz JH, Rojas AM, Schwarzenbacher R: The Nod-like receptor (NLR) family: a tale of similarities and differences. PLoS One 2008, 3: e21191-e211911.
Larsson P, Wallner B, Lindahl E, Elofsson A: Using multiple templates to improve quality of homology models in automated homology modeling. Protein Sci 2008, 17: 990–1002. 10.1110/ps.073344908
Okada K, Arai Y, Hata M, Matsuyama K, Yano M: Big endothelin-1 structure important for specific processing by endothelin-converting enzyme of bovine endothelial cells. Eur J Biochem 1993, 218: 493–498. 10.1111/j.1432-1033.1993.tb18401.x
Klepeis JL, Lindorff-Larsen K, Dror RO, Shaw DE: Long-timescale molecular dynamics simulations of protein structure and function. Curr Opin Struct Biol 2009, 19: 120–127. 10.1016/j.sbi.2009.03.004
Prez A, Marchn I, Svozil D, Sponer J, Cheatham TE 3rd, Laughton CA, Orozco M: Refinement of the AMBER force field for nucleic acids: improving the description of conformers. Biophys J 2007, 92: 3817–3829. 10.1529/biophysj.106.097782
Sansom CE, Hoang MV, Turner AJ: Molecular modelling and site-directed mutagenesis of the active site of endothelin-converting enzyme. Prot Eng 1998, 11: 1235–1241. 10.1093/protein/11.12.1235
Johnson GD, Swenson HR, Ramage R, Ahn K: Mapping the active site of endothelin-converting enzyme-1 through subsite specificity and mutagenesis studies: a comparison with neprilysin. Arch Biochem Biophys 2002, 398: 240–248. 10.1006/abbi.2001.2708
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acid Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389
Clamp M, Cuff J, Searle SM, Barton GJ: The jalview java alignment editor. Bioinformatics 2004, 20: 426–427. 10.1093/bioinformatics/btg430
Molecular Operating Environment (MOE), 2011.10. Canada: ChemicalComputing Group. Inc; 2011. 1010 Sherbooke St. West, Suite #910 Montreal, QC, H3A 2R7 1010 Sherbooke St. West, Suite #910 Montreal, QC, H3A 2R7
Wang J, Cieplak P, Kollman PA: How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J Comput Chem 2000, 21: 1049–1074. 10.1002/1096-987X(200009)21:12<1049::AID-JCC3>3.0.CO;2-F
Laskowski RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr 1993, 26: 283–291. 10.1107/S0021889892009944
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA: Gaussian 09, revision A. 02. Wallingford: Gaussian. Inc; 2009.
Cornell WD, Cieplak P, Bayly CI, Kollmann PA: Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation. J Am Chem Soc 1993, 115: 9620–9631. 10.1021/ja00074a030
Case DA III, Cheatham TE III, Darden T, Gohlke H, Luo R, Merz KM Jr, Onufriev A, Simmerling C, Wang B, Woods RJ: The Amber biomolecular simulation programs. J Comput Chem 2005, 26: 1668–1688. 10.1002/jcc.20290
Aqvist J, Warshel A: Free energy relationships in metalloenzyme-catalyzed reactions. Calculations of the effects of metal ion substitutions in staphylococcal nuclease. J Am Chem Soc 1990, 112: 2860–2868. 10.1021/ja00164a003
Pang YP: Novel zinc protein molecular dynamics simulations: Steps toward antiangiogenesis for cancer treatment. J Mol Mod 1999, 5: 196–202. 10.1007/s008940050119
Pang YP, Xu K, Yazal JEG, Prendergast F: Successful molecular dynamics simulation of the zinc-bound farnesyltransferase using the cationic dummy atom approach. Protein Sci 2000, 9: 1857–1865.
Pang YP: Successful molecular dynamics simulation of two zinc complexes bridged by a hydroxide in phosphotriesterase using the cationic dummy atom method. Protein: Struc Func Bioinfo 2001, 45: 183–189. 10.1002/prot.1138
Yazal JE, Pang YP: Ab initio calculations of proton dissociation energies of zinc ligands: hypothesis of imidazolate as zinc ligand in proteins. J Phys Chem B 1999, 103: 8773–8779. 10.1021/jp991787m
Yazal JE, Roe RR, Pang YP: Zinc’s affect on proton transfer between imidazole and acetate predicted by ab initio calculations. J Phys Chem B 2000, 104: 6662–6667. 10.1021/jp994283x
Yazal JE, Pang YP: Comparison of DFT, Møller-Plesset, and coupled cluster calculations of the proton dissociation energies of imidazole and N-methylacetamide in the presence of zinc (II). J Mol Struct (THEOCHEM) 2001, 545: 271–274. 10.1016/S0166-1280(01)00405-5
Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T: A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J Comput Chem 2003, 24: 1999–2012. 10.1002/jcc.10349
Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA: Development and testing of a general amber force field. J Comput Chem 2004, 25: 1157–1174. 10.1002/jcc.20035
Miyamoto S, Kollman PA: SETTLE: an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J Comput Chem 1992, 13: 952–962. 10.1002/jcc.540130805
Berendsen HJC, Postma JPM, Gunsteren WFV, DiNola A, Haak JR: Molecular dynamics with coupling to an external bath. J Chem Phys 1984, 81: 3684–3690. 10.1063/1.448118
Darden T, York D, Pedersen L: Particle mesh Ewald: An N log (N) method for Ewald sums in large systems. J Chem Phys 1993, 98: 10089–10092. 10.1063/1.464397
Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG: A smooth particle mesh Ewald method. J Chem Phys 1995, 103: 8577–8593. 10.1063/1.470117
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimeraa visualization system for exploratory research and analysis. J Comput Chem 2004, 25: 1605–1612. 10.1002/jcc.20084
Humphrey W, Dalke A, Schulten K: VMD: visual molecular dynamics. J Mol Graph 1996, 14: 33–38. 10.1016/0263-7855(96)00018-5
Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 1983, 22: 2577–2637. 10.1002/bip.360221211
Financial support from Higher Education Commission, Pakistan for providing funds (Project # 20-1509/R&D/10/3218) is highly acknowledged. We also greatly acknowledge the technical support provided by Prof. Bernd M. Rode (University of Innsbruck) during this research. The authors are also grateful to the AMBER supporting team for providing us AMBER software and technical support.
The authors declare that they have no competing interests.
SI and STM designed and performed the experiments, and analyzed the data. SI wrote the paper and STM and Z-ul-H proof read the paper. All research has been done under the supervision of Z-ul-H. All authors read and approved the final manuscript.
Sadaf Iqbal and Syed Tarique Moin contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.