Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease.

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long‐range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51–56 [IG(51–56)] and 50–57 [IG(50–57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the β‐turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283–323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the β‐hairpin structure of peptides corresponding to the C‐terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long‐range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the β‐hairpin structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Structural properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1‐40) (Aβ(1‐40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was ?93.56 ± 6.341 kJ mol^?1. Prevalence of random coil conformations was found for both chains with the exceptions of increased β‐sheet content from residues 16‐21 and 29‐32 of chain A and residues 15‐21 and 30‐33 of chain B with β‐turn/β‐bend conformations in both chains from residues 1‐16, 21‐29 of chain A, 1‐16, and 21‐29 of chain B. There is a mixed β‐turn/β‐sheet region from residues 33‐38 of both chains. Analysis of intra‐ and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β‐sheet configuration with contacts between residues 17‐21 of chain A with residues 17‐21 and 31‐36 of chain B as well as residues 31‐36 of chain A with residues 17‐21 and 31‐36 of chain B. Based on dihedral principal component analysis, the antiparallel β‐sheet‐loop‐β‐sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1‐40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C‐terminal hydrophobic cores with the expected β‐sheet‐loop‐β‐sheet conformational motif. Proteins 2017; 85:1024–1045. © 2017 Wiley Periodicals, Inc.     

Exploring the free energy landscape of a model β‐hairpin peptide and its isoform

Secondary structural transitions from α‐helix to β‐sheet conformations are observed in several misfolding diseases including Alzheimer's and Parkinson's. Determining factors contributing favorably to the formation of each of these secondary structures is therefore essential to better understand these disease states. β‐hairpin peptides form basic components of anti‐parallel β‐sheets and are suitable model systems for characterizing the fundamental forces stabilizing β‐sheets in fibrillar structures. In this study, we explore the free energy landscape of the model β‐hairpin peptide GB1 and its E2 isoform that preferentially adopts α‐helical conformations at ambient conditions. Umbrella sampling simulations using all‐atom models and explicit solvent are performed over a large range of end‐to‐end distances. Our results show the strong preference of GB1 and the E2 isoform for β‐hairpin and α‐helical conformations, respectively, consistent with previous studies. We show that the unfolded states of GB1 are largely populated by misfolded β‐hairpin structures which differ from each other in the position of the β‐turn. We discuss the energetic factors contributing favorably to the formation of α‐helix and β‐hairpin conformations in these peptides and highlight the energetic role of hydrogen bonds and non‐bonded interactions. Proteins 2014; 82:2394–2402. © 2014 Wiley Periodicals, Inc.     

VCD spectroscopic properties of the β‐hairpin forming miniprotein CLN025 in various solvents

Electronic and vibrational circular dichroism are often used to determine the secondary structure of proteins, because each secondary structure has a unique spectrum. Little is known about the vibrational circular dichroic spectroscopic features of the β‐hairpin. In this study, the VCD spectral features of a decapeptide, YYDPETGTWY (CLN025), which forms a stable β‐hairpin that is stabilized by intramolecular weakly polar interactions and hydrogen bonds were determined. Molecular dynamics simulations and ECD spectropolarimetry were used to confirm that CLN025 adopts a β‐hairpin in water, TFE, MeOH, and DMSO and to examine differences in the secondary structure, hydrogen bonds, and weakly polar interactions. CLN025 was synthesized by microwave‐assisted solid phase peptide synthesis with N^α‐Fmoc protected amino acids. The VCD spectra displayed a (?,+,?) pattern with bands at 1640 to 1656 cm^?1, 1667 to 1687 cm^?1, and 1679 to 1686 cm^?1 formed by the overlap of a lower frequency negative couplet and a higher frequency positive couplet. A maximum IR absorbance was observed at 1647 to 1663 cm^?1 with component bands at 1630 cm^?1, 1646 cm^?1, 1658 cm^?1, and 1675 to 1680 cm^?1 that are indicative of the β‐sheet, random meander, either random meander or loop and turn, respectively. These results are similar to the results of others, who examined the VCD spectra of β‐hairpins formed by ^DPro‐Xxx turns and indicated that observed pattern is typical of β‐hairpins. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 442–450, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Structural and sequence features of two residue turns in beta‐hairpins

Beta‐turns in beta‐hairpins have been implicated as important sites in protein folding. In particular, two residue β‐turns, the most abundant connecting elements in beta‐hairpins, have been a major target for engineering protein stability and folding. In this study, we attempted to investigate and update the structural and sequence properties of two residue turns in beta‐hairpins with a large data set. For this, 3977 beta‐turns were extracted from 2394 nonhomologous protein chains and analyzed. First, the distribution, dihedral angles and twists of two residue turn types were determined, and compared with previous data. The trend of turn type occurrence and most structural features of the turn types were similar to previous results, but for the first time Type II turns in beta‐hairpins were identified. Second, sequence motifs for the turn types were devised based on amino acid positional potentials of two‐residue turns, and their distributions were examined. From this study, we could identify code‐like sequence motifs for the two residue beta‐turn types. Finally, structural and sequence properties of beta‐strands in the beta‐hairpins were analyzed, which revealed that the beta‐strands showed no specific sequence and structural patterns for turn types. The analytical results in this study are expected to be a reference in the engineering or design of beta‐hairpin turn structures and sequences. Proteins 2014; 82:1721–1733. © 2014 Wiley Periodicals, Inc.     

Propensities of peptides containing the Asn‐Gly segment to form β‐turn and β‐hairpin structures

The propensities of peptides that contain the Asn‐Gly segment to form β‐turn and β‐hairpin structures were explored using the density functional methods and the implicit solvation model in CH₂Cl₂ and water. The populations of preferred β‐turn structures varied depending on the sequence and solvent polarity. In solution, β‐hairpin structures with βI′ turn motifs were most preferred for the heptapeptides containing the Asn‐Gly segment regardless of the sequence of the strands. These preferences in solution are consistent with the corresponding X‐ray structures. The sequence, H‐bond strengths, solvent polarity, and conformational flexibility appeared to interact to determine the preferred β‐hairpin structure of each heptapeptide, although the β‐turn segments played a role in promoting the formation of β‐hairpin structures and the β‐hairpin propensity varied. In the heptapeptides containing the Asn‐Gly segment, the β‐hairpin formation was enthalpically favored and entropically disfavored at 25°C in water. The calculated results for β‐turns and β‐hairpins containing the Asn‐Gly segment imply that these structural preferences may be useful for the design of bioactive macrocyclic peptides containing β‐hairpin mimics and the design of binding epitopes for protein–protein and protein–nucleic acid recognitions. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 653–664, 2016.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. III. Dynamics of long‐range hydrophobic interactions

A 20‐residue peptide, IG(42–61), derived from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptoccocus was studied using circular dichroism, nuclear magnetic resonance (NMR) spectroscopy at various temperatures and by differential scanning calorimetry (DSC). Unlike other related peptides studied so far, this peptide displays two heat capacity peaks in DSC measurements (at a scanning rate of 1.5 deg/min at a peptide concentration of 0.07 mM), which suggests a three‐state folding/unfolding process. The results from DSC and NMR measurements suggest the formation of a dynamic network of hydrophobic interactions stabilizing the structure, which resembles a β‐hairpin shape over a wide range of temperatures (283–313 K). Our results show that IG (42–61) possesses a well‐organized three‐dimensional structure stabilized by long‐range hydrophobic interactions (Tyr50 ··· Phe57 and Trp48 ··· Val59) at T = 283 K and (Trp48 ··· Val59) at 305 and 313 K. The mechanism of β‐hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Folding of a DNA hairpin loop structure in explicit solvent using replica-exchange molecular dynamics simulations

Hairpin loop structures are common motifs in folded nucleic acids. The 5'-GCGCAGC sequence in DNA forms a characteristic and stable trinucleotide hairpin loop flanked by a two basepair stem helix. To better understand the structure formation of this hairpin loop motif in atomic detail, we employed replica-exchange molecular dynamics (RexMD) simulations starting from a single-stranded DNA conformation. In two independent 36 ns RexMD simulations, conformations in very close agreement with the experimental hairpin structure were sampled as dominant conformations (lowest free energy state) during the final phase of the RexMDs ( approximately 35% at the lowest temperature replica). Simultaneous compaction and accumulation of folded structures were observed. Comparison of the GCA trinucleotides from early stages of the simulations with the folded topology indicated a variety of central loop conformations, but arrangements close to experiment that are sampled before the fully folded structure also appeared. Most of these intermediates included a stacking of the C(2) and G(3) bases, which was further stabilized by hydrogen bonding to the A(5) base and a strongly bound water molecule bridging the C(2) and A(5) in the DNA minor groove. The simulations suggest a folding mechanism where these intermediates can rapidly proceed toward the fully folded hairpin and emphasize the importance of loop and stem nucleotide interactions for hairpin folding. In one simulation, a loop motif with G(3) in syn conformation (dihedral flip at N-glycosidic bond) accumulated, resulting in a misfolded hairpin. Such conformations may correspond to long-lived trapped states that have been postulated to account for the folding kinetics of nucleic acid hairpins that are slower than expected for a semiflexible polymer of the same size.     

Interplaying factors for the formation of photoswitchable β‐hairpins: the advantage of a flexible switch

A series of peptidomimetics intended to promote the β‐hairpin motif have been studied. Structural variations include a turn region with and without a photoswitchable chromophore, and strands with amino acid side chains supporting various degrees of interstrand interactions for hairpin stabilisation. The propensity of the compounds to form β‐hairpins was evaluated experimentally by NMR spectroscopy, translational self‐diffusion studies and CD spectroscopy. In the presence of hairpin stabilising interstrand interactions, the structurally flexible stilbene chromophore appeared to be well compatible with the imposed secondary structure. Copyright © 2008 European Peptide Society and John Wiley & Sons, Ltd.     

β‐Hairpin stabilization in a 28‐residue peptide derived from the β‐subunit sequence of human chorionic gonadotropin hormone

The β‐subunit of the human chorionic gonadotropin (hCG) hormone, which is believed to be related to certain types of cancer, contains three hairpin‐like fragments. To investigate the role of β‐hairpin formation in the early stages of the hCGβ folding, a 28‐residue peptide with the sequence RDVRFESIRLPGSPRGVNPVVSYAVALS, corresponding to the H3‐β hairpin fragment (residues 60–87) of the hCGβ subunit, was studied under various conditions using three optical spectroscopic methods: Fourier transform ir spectroscopy, electronic CD, and vibrational CD. Environmental conditions are critical factors for formation of secondary structure in this peptide. TFE : H₂O mixed solvents induced helical formation. Formation of β‐structure in this peptide, which may be related to the native β‐hairpin formation in the intact hormone, was found to be induced only under conditions such as high concentration, high temperature, and the presence of nonmicellar sodium dodecyl sulfate concentrations. These findings support a protein folding mechanism for the hCGβ subunit in which an initial hydrophobic collapse, which increases intermolecular interactions in hCGβ, is needed to induce the H3‐β hairpin formation. © 1999 John Wiley & Sons, Inc. Biopoly 50: 413–423, 1999     

Optimization of a β‐sheet‐cap for long loop closure

Protein loops make up a large portion of the secondary structure in nature. But very little is known concerning loop closure dynamics and the effects of loop composition on fold stability. We have designed a small system with stable β‐sheet structures, including features that allow us to probe these questions. Using paired Trp residues that form aromatic clusters on folding, we are able to stabilize two β‐strands connected by varying loop lengths and composition (an example sequence: R W ITVTI – loop – KKIRV W E). Using NMR and CD, both fold stability and folding dynamics can be investigated for these systems. With the 16 residue loop peptide (sequence: R W ITVTI‐(GGGGKK)₂GGGG‐KKIRV W E) remaining folded (ΔG_U = 1.6 kJ/mol at 295K). To increase stability and extend the series to longer loops, we added an additional Trp/Trp pair in the loop flanking position. With this addition to the strands, the 16 residue loop (sequence: R W ITVRI W ‐(GGGGKK)₂GGGG‐ W KTIRV W E) supports a remarkably stable β‐sheet (ΔG_U = 6.3 kJ/mol at 295 K, T_m = ～55°C). Given the abundance of loops in binding motifs and between secondary structures, these constructs can be powerful tools for peptide chemists to study loop effects; with the Trp/Trp pair providing spectroscopic probes for assessing both stability and dynamics by NMR.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of β‐hairpin structure

We previously studied a 16‐amino acid‐residue fragment of the C‐terminal β‐hairpin of the B3 domain (residues 46–61), [IG(46–61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46–61) by systematically shortening the peptide by one residue at a time from both the C‐ and the N‐terminus. To determine the structure and stability of two resulting 12‐ and 14‐amino acid‐residue peptides, IG(48–59) and IG(47–60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48–59) possesses organized three‐dimensional structure stabilized by hydrophobic interactions (Tyr50–Phe57 and Trp48–Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (T_m = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47–60) determined by DSC is T_m = 330 K and its structure is similar to that of the native β‐hairpin at all (lower) temperatures examined (283–313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a β‐hairpin structural element. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Refolding simulations of an isolated fragment of barnase into a native-like β hairpin: Evidence for compactness and hydrogen bonding as concurrent stabilizing factors

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a β hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptide conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the β hairpin. All this lends reason to believe that the 85–102 barnase fragment is a strong initiation site for folding. Proteins 29:212–227, 1997. © 1997 Wiley-Liss, Inc.     

Effects of force fields on the conformational and dynamic properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1‐40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER‐ff99sb‐ILDN, AMBER‐ff99sb*‐ILDN, AMBER‐ff99sb‐NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER‐ff99sb‐ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α‐helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER‐ff99sb‐NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER‐ff99sb‐NMR force field, the others tended to under estimate the expected amount of β‐sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER‐ff99sb‐NMR, reproduce a theoretically expected β‐sheet‐turn‐β‐sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C‐terminal hydrophobic cores from residues 17‐21 and 30‐36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.     

Conformational dynamics of a short antigenic peptide in its free and antibody bound forms gives insight into the role of β‐turns in peptide immunogenicity

Earlier immunological experiments with a synthetic 36‐residue peptide (75‐110) from Influenza hemagglutinin have been shown to elicit anti‐peptide antibodies (Ab) which could cross‐react with the parent protein. In this article, we have studied the conformational features of a short antigenic (Ag) peptide (⁹⁸YPYDVPDYASLRS¹¹⁰) from Influenza hemagglutinin in its free and antibody (Ab) bound forms with molecular dynamics simulations using GROMACS package and OPLS‐AA/L all‐atom force field at two different temperatures (293 K and 310 K). Multiple simulations for the free Ag peptide show sampling of ordered conformations and suggest different conformational preferences of the peptide at the two temperatures. The free Ag samples a conformation crucial for Ab binding (β‐turn formed by “DYAS” sequence) with greater preference at 310 K while, it samples a native‐like conformation with relatively greater propensity at 293 K. The sequence “DYAS” samples β‐turn conformation with greater propensity at 310 K as part of the hemagglutinin protein also. The bound Ag too samples the β‐turn involving “DYAS” sequence and in addition it also samples a β‐turn formed by the sequence “YPYD” at its N‐terminus, which seems to be induced upon binding to the Ab. Further, the bound Ag displays conformational flexibility at both 293 K and 310 K, particularly at terminal residues. The implications of these results for peptide immunogenicity and Ag–Ab recognition are discussed. Proteins 2015; 83:1352–1367. © 2015 Wiley Periodicals, Inc.     

RNA simulations: probing hairpin unfolding and the dynamics of a GNRA tetraloop

Simulations of an RNA hairpin containing a GNRA tetraloop were conducted to allow the characterization of its secondary structure formation and dynamics. Ten 10 ns trajectories of the folded hairpin 5'-GGGC[GCAA]GCCU-3' were generated using stochastic dynamics and the GB/SA implicit solvent model at 300 K. Overall, we find the stem to be a very stable subunit of this molecule, whereas multiple loop conformations and transitions between them were observed. These trajectories strongly suggest that extension of the C6 base away from the loop occurs cooperatively with an N-type-->S-type sugar pucker conversion in that residue and that similar pucker transitions are necessary to stabilize other looped-out bases. In addition, a short-lived conformer with an extended fourth loop residue (A8) lacking this stabilizing 2'-endo pucker mode was observed. Results of thermal perturbation at 400 K support this model of loop dynamics. Unfolding trajectories were produced using this same methodology at temperatures of 500 to 700 K. The observed unfolding events display three-state behavior kinetically (including native, globular, and unfolded populations) and, based on these observations, we propose a folding mechanism that consists of three distinct events: (i) collapse of the random unfolded structure and sampling of the globular state; (ii) passage into the folded region of configurational space as stem base-pairs form and gain helicity; and (iii) attainment of proper loop geometry and organization of loop pairing and stacking interactions. These results are considered in the context of current experimental knowledge of this and similar nucleic acid hairpins.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

The ensemble folding kinetics of the FBP28 WW domain revealed by an all-atom Monte Carlo simulation in a knowledge-based potential

In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous P_fold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc.