Molecular dynamics simulation studies on structural and conformational changes in tyrosine-67 nitrated cytochrome c

Protein tyrosine nitration is well-established post-translational modification occurring in a number of diseases, viz. neurodegenerative, cardiovascular diseases, ageing, etc. Tyrosine-67 (Tyr-67) nitration of cytochrome c (cyt c) was observed under oxidative stress affecting its structure and electron transfer properties. Hence, in this study, molecular dynamics (MD) simulations were carried out at room temperature to investigate the structural and conformational changes in the nitrated cyt c's. MD results revealed that the bond between FE (Heme-105) and S (Met-80) considerably weakened, radius of gyration, backbone and Cα root-mean-square deviations decreased and hydrogen bonding increased in the nitrated cyt c's relative to wild type (WT) cyt c. Ramachandran plot analysis revealed that N- and C-terminal helices also affected by nitration at CE2 carbon atom. Furthermore, essential dynamics analysis showed that amplitude of concerted motion decreased in the nitrated cyt c's, perhaps due to the increase in the hydrogen bonding interaction. Taken together, the structural and conformational changes in the active site Tyr-67 nitrated cyt c may have implications in the loss of electron/proton transfer and gain of apoptotic properties.     

Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

The relation between the structure of the first solvation shell and the IR spectra of aqueous solutions

The spectroscopic signatures of solvated anions and cations, in the O-H stretch region of water, are studied using the POLIR potential. Shifts in the spectra are shown to correlate very well with the distribution of a particular hydrogen bond angle for the waters in the first solvation shell. The results indicate that the spectral shifts might be predicted from MD simulations in a computationally convenient fashion, avoiding an explicit calculation of the spectra, as first suggested by Sharp et al. (J Chem Phys 114(4):1791–1796, 2001).     

Molecular dynamics analysis of the engrailed homeodomain-DNA recognition

Molecular dynamics (MD) simulations were performed for investigating the role of Gln50 in the engrailed homeodomain-DNA recognition. Employing the crystal structure of free engrailed homeodomain and homeodomain-DNA complex as a starting structure, we carried out MD simulations of: (i) the complex between engrailed homeodomain and a 20 base-pair DNA containing TAATTA core sequence; (ii) the free engrailed homeodomain. The simulations show that homeodomain flexibility does not depend on its ligation state. The engrailed homeodomain shows similar flexibility, and the recognition helix-3 shows very similar characteristic of high rigidity and limited conformational space in two complexation states. At the same time, DNA structure has also no obvious conformational fluctuations. These results preclude the possibility of the side chain of Gln50 forming direct hydrogen bonds to the core DNA bases. MD simulations confirm a few well-conserved sites for water-mediated hydrogen bonds from protein to DNA are occupied by water molecules, and Gln50 interacts with corresponding core DNA bases through water-mediated hydrogen bonds. So Gln50 plays a relatively modest role in determining the affinity and specificity of the engrailed homeodomain. In addition, the electrostatic interaction between homeodomain and phosphate backbone of the DNA is a main factor for N- and C-terminal arm becoming ordered upon DNA binding.     

A further investigation of the cytochrome<Emphasis Type="Italic"> b</Emphasis><Subscript>5</Subscript>–cytochrome<Emphasis Type="Italic"> c</Emphasis> complex

The interaction of reduced rabbit cytochrome b₅ with reduced yeast iso-1 cytochrome c has been studied through the analysis of ¹H–¹⁵N HSQC spectra, of ¹⁵N longitudinal (R₁) and transverse (R₂) relaxation rates, and of the solvent exchange rates of protein backbone amides. For the first time, the adduct has been investigated also from the cytochrome c side. The analysis of the NMR data was integrated with docking calculations. The result is that cytochrome b₅ has two negative patches capable of interacting with a single positive surface area of cytochrome c. At low protein concentrations and in equimolar mixture, two different 1:1 adducts are formed. At high concentration and/or with excess cytochrome c, a 2:1 adduct is formed. All the species are in fast exchange on the scale of differences in chemical shift. By comparison with literature data, it appears that the structure of one 1:1 adduct changes with the origin or primary sequence of cytochrome b₅.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.Abbreviations HSQC heteronuclear single quantum correlation spectroscopy - MD molecular dynamics     

Coupling between the BLUF and EAL domains in the blue light-regulated phosphodiesterase BlrP1

The first biochemical and structural characterization of the full-length active photoreceptor BlrP1 from Klebsiella pneumoniae was recently reported by Barends et al. [Nature 459:1015–1018, (2009)]. The light-regulated catalytic function of its C-terminal c-di-guanosine monophosphate phosphodiesterase, the EAL (Glu-Ala-Leu) domain, is activated by the N-terminal sensor of blue light using the flavin adenine dinucleotide (BLUF) domain. We performed molecular dynamics simulations on the dimeric BlrP1 protein in order to examine the coupling regions that are presumably involved in transmitting light-induced structural changes which occur in the BLUF domain to the EAL domain. According to the results of simulations and an analysis of the hydrogen bonding between the respective polypeptide chains, the region containing the site on the α3α4 loop of BLUF is responsible for communication between the photosensing and catalytic domains in the dimeric BlrP1 protein.     

The asparagine-stabilized beta-turn of apamin: contribution to structural stability from dynamics simulation and amide hydrogen exchange analysis

Molecular dynamics simulations of bee venom apamin, and an analogue having an Asn to Ala substitution at residue 2 (apamin-N2A), were analyzed to explore the contribution of hydrogen bonds involving Asn2 to local (beta-turn residues N2, C3, K4, A5) and global stability. The wild-type peptide retained a stable conformation during 2.4 ns of simulation at 67 degrees C, with high beta-turn stability characterized by backbone-side chain hydrogen bonds involving beta-turn residues K4 and A5, with the N2 side chain amide carbonyl. The loss of stabilizing interactions involving the N2 side chain resulted in the loss of the beta-turn conformation in the apamin N2A simulations (27 or 67 degrees C). This loss of beta-turn stability propagates throughout the peptide structure, with destabilization of the C-terminal helix connected to the N-terminal region by two disulfide bonds. Backbone stability in a synthetic peptide analogue (apamin-N2A) was characterized by NMR and amide hydrogen exchange measurements. Consistent with the simulations, loss of hydrogen bonds involving the N2 side chain resulted in destabilization of both the N-terminal beta-turn and the C-terminal helix. Amide exchange protection factors in the C-terminal helix were reduced by 9-11-fold in apamin N2A as compared with apamin, corresponding to free energy (deltaDeltaG(uf)) of around 1.5 kcal M(-1) at 20 degrees C. This is equivalent to the contribution of hydrogen bond interactions involving the N2 side chain to the stability of the beta-turn. Together with additional measures of exchange protection factors, the three main contributions to backbone stability in apamin that account for virtually the full thermodynamic stability of the peptide have been quantitated.     

Molecular modeling studies of the DCCD-treated cytochrome bc1 complex: predicted conformational changes and inhibition of proton translocation

Dicyclohexylcarbodiimide (DCCD) binds covalently to an acidic amino acid located in the cd loop connecting membrane-spanning helices C and D of cytochrome b resulting in an inhibition of proton translocation in the cytochrome bc ₁ complex with minimal effects on the steady state rate of electron transfer. Single turnover studies performed with the yeast cytochrome bc ₁ complex indicated that the initial phase of cytochrome b reduction was inhibited 25–45% in the DCCD-treated cytochrome bc ₁ complex, while the rate of cytochrome c ₁ reduction was unaffected. Simulations by molecular modeling predict that binding of DCCD to glutamate 163 located in the cd2 loop of cytochrome b of chicken liver mitochondria results in major conformational changes in the protein. The conformation of the cd loop and the end of helix C appeared twisted with a concomitant rearrangement of the amino acid residues of both cd1 and cd2 loops. The predicted rearrangement of the amino acid residues of the cd loop results in disruptions of the hydrogen bonds predicted to form between amino acid residues of the cd and ef loops. Simultaneously, two new hydrogen bonds are predicted to form between glutamate 272 and two residues, aspartate 253 and tyrosine 272. Formation of these new hydrogen bonds would restrict the rotation and protonation of glutamate 272, which may be necessary for the release of the second electrogenic proton obtained during ubiquinol oxidation in the bc₁ complex.     

Submolecular unfolding units of <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis> cytochrome <Emphasis Type="Italic">c</Emphasis>-551

Hydrogen exchange rates for backbone amide protons of oxidized Pseudomonas aeruginosa cytochrome c-551 (P. aeruginosa cytochrome c) have been measured in the presence of low concentrations of the denaturant guanidine hydrochloride. Analysis of the data has allowed identification of submolecular unfolding units known as foldons. The highest-energy foldon bears similarity to the proposed folding intermediate for P. aeruginosa cytochrome c. Parallels are seen to the foldons of the structurally homologous horse cytochrome c, although the heme axial methionine-bearing loop has greater local stability in P. aeruginosa cytochrome c, in accord with previous folding studies. Regions of low local stability are observed to correspond with regions that interact with redox partners, providing a link between foldon properties and function. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Kinetic and crystallographic studies of a redesigned manganese-binding site in cytochrome c peroxidase

Manganese peroxidase (MnP) from the white rot fungus Phanerochaete chrysosporium contains a manganese-binding site that plays a critical role in its function. Previously, a Mn^II-binding site was designed into cytochrome c peroxidase (CcP) based on sequence homology (Yeung et al. in Chem. Biol. 4:215–222, 1997; Gengenbach et al. in Biochemistry 38:11425–11432, 1999). Here, we report a redesign of this site based on X-ray structural comparison of MnP and CcP. The variant, CcP(D37E, V45E, H181E), displays 2.5-fold higher catalytic efficiency (k _cat/K _M) than the variant in the original design, mostly due to a stronger K _M of 1.9 mM (vs. 4.1 mM). High-resolution X-ray crystal structures of a metal-free form and a form with Co^II at the designed Mn^II site were also obtained. The metal ion in the engineered metal-binding site overlays well with Mn^II bound in MnP, suggesting that this variant is the closest structural model of the Mn^II-binding site in MnP for which a crystal structure exists. A major difference arises in the distances of the ligands to the metal; the metal–ligand interactions in the CcP variant are much weaker than the corresponding interactions in MnP, probably owing to partial occupancy of metal ion at the designed site, difference in the identity of metal ions (Co^II rather than Mn^II) and other interactions in the second coordination sphere. These results indicate that the metal ion, the ligands, and the environment around the metal-binding site play important roles in tuning the structure and function of metalloenzymes. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Role of disulfide bonds in modulating internal motions of proteins to tune their function: molecular dynamics simulation of scorpion toxin Lqh III

A series of 1-ns MD simulations were performed on the scorpion toxin Lqh III in native and disulfide bond broken states. The removal of disulfide bonds has caused hydrogen bond network alteration in the five-residue turn, the long loop, the alpha-helix, the loop connecting strands II and III, and the C-terminal region. In addition and more importantly, it has influenced the amplitude of the fluctuations of five-residue turn, loops, and C-terminal region with a minor effect on the fluctuations of the cysteines in the broken bond sites. These findings suggest that disulfide bonds are not the most important factors in rigidifying their own locations, while they have more important effects at a global scale. Furthermore, our results reveal that disulfide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the binding site residues. Therefore, we can conclude that disulfide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin molecules. Although this conclusion cannot be generalized to all peptides and proteins, it demonstrates the importance of more investigations on this aspect of disulfide bond efficacy.     

Protein folding: could hydrophobic collapse be coupled with hydrogen-bond formation?

A judicious examination of an exhaustive PDB sample of soluble globular proteins of moderate size (N<102) reveals a commensurable relationship between hydrophobic surface burial and number of backbone hydrogen bonds. An analysis of 50,000 conformations along the longest all-atom MD trajectory allows us to infer that not only the hydrophobic collapse is concurrent with the formation of backbone amide-carbonyl hydrogen bonds, they are also dynamically coupled processes. In statistical terms, hydrophobic clustering of the side chains is inevitably conducive to backbone burial and the latter process becomes thermodynamically too costly and kinetically unfeasible without amide-carbonyl hydrogen-bond formation. Furthermore, the desolvation of most hydrogen bonds is exhaustive along the pathway, implying that such bonds guide the collapse process.     

Domain-Swapped Dimer of Pseudomonas aeruginosa Cytochrome c 551: Structural Insights into Domain Swapping of Cytochrome c Family Proteins

Cytochrome c (cyt c) family proteins, such as horse cyt c, Pseudomonas aeruginosa cytochrome c ₅₅₁ (PA cyt c ₅₅₁), and Hydrogenobacter thermophilus cytochrome c ₅₅₂ (HT cyt c ₅₅₂), have been used as model proteins to study the relationship between the protein structure and folding process. We have shown in the past that horse cyt c forms oligomers by domain swapping its C-terminal helix, perturbing the Met–heme coordination significantly compared to the monomer. HT cyt c ₅₅₂ forms dimers by domain swapping the region containing the N-terminal α-helix and heme, where the heme axial His and Met ligands belong to different protomers. Herein, we show that PA cyt c ₅₅₁ also forms domain-swapped dimers by swapping the region containing the N-terminal α-helix and heme. The secondary structures of the M61A mutant of PA cyt c ₅₅₁ were perturbed slightly and its oligomer formation ability decreased compared to that of the wild-type protein, showing that the stability of the protein secondary structures is important for domain swapping. The hinge loop of domain swapping for cyt c family proteins corresponded to the unstable region specified by hydrogen exchange NMR measurements for the monomer, although the swapping region differed among proteins. These results show that the unstable loop region has a tendency to become a hinge loop in domain-swapped proteins.     

Tadpole-like Conformations of Huntingtin Exon 1 Are Characterized by Conformational Heterogeneity that Persists regardless of Polyglutamine Length

Soluble huntingtin exon 1 (Httex1) with expanded polyglutamine (polyQ) engenders neurotoxicity in Huntington's disease. To uncover the physical basis of this toxicity, we performed structural studies of soluble Httex1 for wild-type and mutant polyQ lengths. Nuclear magnetic resonance experiments show evidence for conformational rigidity across the polyQ region. In contrast, hydrogen–deuterium exchange shows absence of backbone amide protection, suggesting negligible persistence of hydrogen bonds. The seemingly conflicting results are explained by all-atom simulations, which show that Httex1 adopts tadpole-like structures with a globular head encompassing the N-terminal amphipathic and polyQ regions and the tail encompassing the C-terminal proline-rich region. The surface area of the globular domain increases monotonically with polyQ length. This stimulates sharp increases in gain-of-function interactions in cells for expanded polyQ, and one of these interactions is with the stress-granule protein Fus. Our results highlight plausible connections between Httex1 structure and routes to neurotoxicity.     

A comparison of the dynamic behavior of monomeric and dimeric insulin shows structural rearrangements in the active monomer

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.     

Crystal Structure of the Electron Carrier Domain of the Reaction Center Cytochrome cz Subunit from Green Photosynthetic Bacterium Chlorobium tepidum

In green sulfur photosynthetic bacteria, the cytochrome c_z (cyt c_z) subunit in the reaction center complex mediates electron transfer mainly from menaquinol/cytochrome c oxidoreductase to the special pair (P840) of the reaction center. The cyt c_z subunit consists of an N-terminal transmembrane domain and a C-terminal soluble domain that binds a single heme group. The periplasmic soluble domain has been proposed to be highly mobile and to fluctuate between oxidoreductase and P840 during photosynthetic electron transfer. We have determined the crystal structure of the oxidized form of the C-terminal functional domain of the cyt c_z subunit (C-cyt c_z) from thermophilic green sulfur bacterium Chlorobium tepidum at 1.3-Å resolution. The overall fold of C-cyt c_z consists of four α-helices and is similar to that of class I cytochrome c proteins despite the low similarity in their amino acid sequences. The N-terminal structure of C-cyt c_z supports the swinging mechanism previously proposed in relation with electron transfer, and the surface properties provide useful information on possible interaction sites with its electron transfer partners. Several characteristic features are observed for the heme environment: These include orientation of the axial ligands with respect to the heme plane, surface-exposed area of the heme, positions of water molecules, and hydrogen-bond network involving heme propionate groups. These structural features are essential for elucidating the mechanism for regulating the redox state of cyt c_z.     

Temperature coefficients of amide proton NMR resonance frequencies in trifluoroethanol: A monitor of intramolecular hydrogen bonds in helical peptides?

Summary 2D ¹H NMR spectroscopy of two -helical peptides which differ in their amphipathicity has been used to investigate the relationships between amide-proton chemical shifts, amide-proton exchange rates, temperature, and trifluoroethanol (TFE) concentration. In 50% TFE, in which the peptides are maximally helical, the amide-proton chemical shift and temperature coefficient patterns are very similar to each other in each peptide. Temperature coefficients from –10 to –6 ppb/K, usually indicative of the lack of intramolecular hydrogen bonds, were observed even for hydrophobic amino acids in the center of the -helices. However, slow hydrogen isotope exchange for residues from 4 to 16 in both 18-mer helices indicates intact intramolecular hydrogen bonds over most of the length of these peptides. Based on these anomalous observations, we suggest that the pattern of amide-proton shifts in -helices in H₂O/TFE solvents is dominated by bifurcated intermolecular hydrogen-bond formation between the backbone carbonyl groups and TFE. The amide-proton chemical shift changes with increasing temperature may be interpreted by a disruption of intermolecular hydrogen bonds between carbonyl groups and the TFE in TFE/water rather than by the length of intramolecular hydrogen bonds in -helices. Supplementary Material is available upon request, comprising seven pages with listings of experimental details and the NMR shift data for the two peptides.