Molecular dynamics simulations of beta-hairpin folding

A systematic analysis of the beta hairpin motif in the Protein Data Bank

The beta hairpin motif is a ubiquitous protein structural motif that can be found in molecules across the tree of life. This motif, which is also popular in synthetically designed proteins and peptides, is known for its stability and adaptability to broad functions. Here, we systematically probe all 49,000 unique beta hairpin substructures contained within the Protein Data Bank (PDB) to uncover key characteristics correlated with stable beta hairpin structure, including amino acid biases and enriched interstrand contacts. We find that position specific amino acid preferences, while seen throughout the beta hairpin structure, are most evident within the turn region, where they depend on subtle turn dynamics associated with turn length and secondary structure. We also establish a set of broad design principles, such as the inclusion of aspartic acid residues at a specific position and the careful consideration of desired secondary structure when selecting residues for the turn region, that can be applied to the generation of libraries encoding proteins or peptides containing beta hairpin structures.

A computational approach to get insights into multiple faces of additives in modulation of protein aggregation pathways

An enormous population worldwide is presently confronted with debilitating neurodegenerative diseases. The etiology of the disease is connected to protein aggregation and the events involved therein. Thus, a complete understanding of an inhibitor at different stages in the process is imperative for the formulation of a drug molecule. This review presents a detailed summary of the current status of different cosolvents. It further develops how the complex aggregation pathway can be simplified into three steps common to all proteins and the way computer simulations can be exploited to gain insights into the ways by which known inhibitors can affect all these stages. Computation of theoretical parameters in this regard and their correlation with experimental techniques is accentuated. In addition to providing an outline of the scope of different additives, this review showcases the way by which the problem of analyzing an effect of an additive can be addressed effectively via MD simulations.

Interplay of secondary structures and side-chain contacts in the denatured state of BBA1

The denatured state of a miniprotein BBA1 is studied under the native condition with the AMBER/Poisson-Boltzmann energy model and with the self-guided enhanced sampling technique. Forty independent trajectories are collected to sample the highly diversified denatured structures. Our simulation data show that the denatured BBA1 contains high percentage of native helix and native turn, but low percentage of native hairpin. Conditional population analysis indicates that the native helix formation and the native hairpin formation are not cooperative in the denatured state. Side-chain analysis shows that the native hydrophobic contacts are more preferred than the non-native hydrophobic contacts in the denatured BBA1. In contrast, the salt-bridge contacts are more or less nonspecific even if their populations are higher than those of hydrophobic contacts. Analysis of the trajectories shows that the native helix mostly initiates near the N terminus and propagates to the C terminus, and mostly forms from 3(10)-helix/turn to alpha helix. The same analysis shows that the native turn is important but not necessary in its formation in the denatured BBA1. In addition, the formations of the two strands in the native hairpin are rather asymmetric, demonstrating the likely influence of the protein environment. Energetic analysis shows that the native helix formation is largely driven by electrostatic interactions in denatured BBA1. Further, the native helix formation is associated with the breakup of non-native salt-bridge contacts and the accumulation of native salt-bridge contacts. However, the native hydrophobic contacts only show a small increase upon the native helix formation while the non-native hydrophobic contacts stay essentially the same, different from the evolution of hydrophobic contacts observed in an isolated helix folding.

β-Hairpins with native-like and non-native hydrogen bonding patterns could form during the refolding of staphylococcal nuclease

Refolding of staphylococcal nuclease has been studied recently by hydrogen-deuterium exchange and NMR spectroscopy. These studies infer that beta-hairpin formed by strand 2 and strand 3 connected by reverse turn forms early during the refolding of nuclease. Typically, hydrogen-deuterium exchange NMR techniques are usually carried out on a time scale of milliseconds whereas beta-hairpins are known to fold on a much shorter time scale. It follows that in the experiments, the hydrogen-deuterium exchange protection patterns could be arising from a significant population of fully formed hairpins. In order to demonstrate it is the fully formed hairpins which gives rise to the hydrogen-deuterium exchange protection patterns, we have considered molecular dynamics simulation of the peptide (21)DTVKLMYKGQPMTFR(35) from staphylococcal nuclease corresponding to the beta-hairpin region, using GROMOS96 force field under NVT conditions. Starting from unfolded conformational states, the peptide folds into hairpin conformations with native-like and non-native hydrogen bonding patterns. Subsequent to folding, equilibrium conditions prevail. The computed protection factors and atom depth values, at equilibrium, of the various amide protons agree qualitatively with experimental observations. A collection of molecules following the trajectories observed in the simulations can account for experimental observations. These simulations provide a molecular picture of the formed hairpins and their conformational features during the refolding experiments on nuclease, monitored by hydrogen-deuterium exchange.

Role of hydrophobic interactions and salt-bridges in β-hairpin folding

Beta-hairpins are the simplest form of beta-sheets which, due to the presence of long-range interactions, can be considered as tertiary structures. Molecular dynamics simulation is a powerful tool that can unravel whole pathways of protein folding/unfolding at atomic resolution. We have performed several molecular dynamics simulations, to a total of over 250 ns, of a beta-hairpin peptide in water using GROMACS. We show that hydrophobic interactions are necessary for initiating the folding of the peptide. Once formed, the peptide is stabilized by hydrogen bonds and disruption of hydrophobic interactions in the folded peptide does not denature the structure. In the absence of hydrophobic interactions, the peptide fails to fold. However, the introduction of a salt-bridge compensates for the loss of hydrophobic interactions to a certain extent.

Thermodynamic and kinetic characterization of a β-hairpin peptide in solution: An extended phase space sampling by molecular dynamics simulations in explicit water

The folding of the amyloidogenic H1 peptide MKHMAGAAAAGAVV taken from the syrian hamster prion protein is explored in explicit aqueous solution at 300 K using long time scale all-atom molecular dynamics simulations for a total simulation time of 1.1 mus. The system, initially modeled as an alpha-helix, preferentially adopts a beta-hairpin structure and several unfolding/refolding events are observed, yielding a very short average beta-hairpin folding time of approximately 200 ns. The long time scale accessed by our simulations and the reversibility of the folding allow to properly explore the configurational space of the peptide in solution. The free energy profile, as a function of the principal components (essential eigenvectors) of motion, describing the main conformational transitions, shows the characteristic features of a funneled landscape, with a downhill surface toward the beta-hairpin folded basin. However, the analysis of the peptide thermodynamic stability, reveals that the beta-hairpin in solution is rather unstable. These results are in good agreement with several experimental evidences, according to which the isolated H1 peptide adopts very rapidly in water beta-sheet structure, leading to amyloid fibril precipitates [Nguyen et al., Biochemistry 1995;34:4186-4192; Inouye et al., J Struct Biol 1998;122:247-255]. Moreover, in this article we also characterize the diffusion behavior in conformational space, investigating its relations with folding/unfolding conditions.

Complex folding pathways in a simple beta-hairpin

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin.

Beta-hairpin folding mechanism of a nine-residue peptide revealed from molecular dynamics simulations in explicit water

The beta-hairpin fold mechanism of a nine-residue peptide, which is modified from the beta-hairpin of alpha-amylase inhibitor tendamistat (residues 15-23), is studied through direct folding simulations in explicit water at native folding conditions. Three 300-nanosecond self-guided molecular dynamics (SGMD) simulations have revealed a series of beta-hairpin folding events. During these simulations, the peptide folds repeatedly into a major cluster of beta-hairpin structures, which agree well with nuclear magnetic resonance experimental observations. This major cluster is found to have the minimum conformational free energy among all sampled conformations. This peptide also folds into many other beta-hairpin structures, which represent some local free energy minimum states. In the unfolded state, the N-terminal residues of the peptide, Tyr-1, Gln-2, and Asn-3, have a confined conformational distribution. This confinement makes beta-hairpin the only energetically favored structure to fold. The unfolded state of this peptide is populated with conformations with non-native intrapeptide interactions. This peptide goes through fully hydrated conformations to eliminate non-native interactions before folding into a beta-hairpin. The folding of a beta-hairpin starts with side-chain interactions, which bring two strands together to form interstrand hydrogen bonds. The unfolding of the beta-hairpin is not simply the reverse of the folding process. Comparing unfolding simulations using MD and SGMD methods demonstrate that SGMD simulations can qualitatively reproduce the kinetics of the peptide system.

Expected and unexpected results from combined β-hairpin design elements

A model beta-hairpin dodecapeptide [EFGWVpGKWTIK] was designed by including a favorable D-ProGly Type II' beta-turn sequence and a Trp-zip interaction, while also incorporating a beta-strand unfavorable glycine residue in the N-terminal strand. This peptide is highly folded and monomeric in aqueous solution as determined by combined analysis with circular dichroism and 1H NMR spectroscopy. A peptide representing the folded conformation of the model beta-hairpin [cyclic(EFGWVpGKWTIKpG)] and a linear peptide representing the unfolded conformation [EFGWVPGKWTIK] yield unexpected relative deviations between the CD and 1H NMR spectroscopic results that are attributed to variations in the packing interactions of the aromatic side chains. Mutational analysis of the model beta-hairpin indicates that the Trp-zip interaction favors folding and stability relative to an alternate hydrophobic cluster between Trp and Tyr residues [EFGYVpGKWTIK]. The significance of select diagonal interactions in the model beta-hairpin was tested by rearranging the cross-strand hydrophobic interactions to provide a folded peptide [EWFGIpGKTYWK] displaying evidence of an unusual backbone conformation at the hydrophobic cluster. This unusual conformation does not appear to be a result of the glycine residue in the beta-strand, as replacement with a serine results in a peptide [EWFSIpGKTYWK] with a similar and seemingly characteristic CD spectrum. However, an alternate arrangement of hydrophobic residues with a Trp-zip interaction in a similar position to the parent beta-hairpin [EGFWVpGKWITK] results in a folded beta-hairpin conformation. The differences between side chain packing of these peptides precludes meaningful thermodynamic analysis and illustrates the caution necessary when interpreting beta-hairpin folding thermodynamics that are driven, at least in part, by aromatic cross strand interactions.

Possible locally driven folding pathways of TC5b, a 20-residue protein

A novel computational procedure for modeling possible locally driven folding pathways by stepwise elongations of the peptide chain was successfully applied to TC5b, a 20-residue miniprotein. Systematic exploration of the possible locally driven pathways showed that the Trp-cage structure of TC5b could be obtained by stepwise elongation starting from the noncentral local nucleation centers preexisting in the unfolded state of TC5b. The probable locally driven folding pathway starts with folding of alpha-helical fragment 4-9, followed by formation of the proper three-dimensional structure of fragment 4-12, and then 4-18. Accordingly, the Trp-cage-forming interactions emerge successively, first Trp(6)-Pro(12), then Trp(6)-Pro(18), and then Trp(6)-Tyr(3). The Trp-cage-like structures of TC5b found in this study by independent energy calculations are in excellent agreement with the NMR experimental data. The same procedure rationalizes the incomplete Trp-cage formation observed for two analogs of TC5b. Generally, the success of this novel approach is encouraging and provides some justification for the use of computational simulations of locally driven protein folding.

Real-time NMR kinetic studies provide global and residue-specific information on the non-cooperative unfolding of the beta-trefoil protein, interleukin-1beta

The interleukin-1beta (IL-1beta) structural motif is a beta-trefoil super fold created by six two-stranded beta-hairpins. Turns are thus particularly important in creating the topology and the arrangement of beta-strands in this structural motif. In contrast to the signals observed in optical studies, real-time NMR kinetic investigations of the denaturant-induced unfolding of interleukin-1beta provide direct, global, and residue-specific information on the structural nature of the unfolding reaction. Heterogeneity in the individual amino acid residue kinetics reveals a rugged unfolding landscape. The relative kinetic stability of native-like turns supports low temperature molecular dynamics predictions of turn-controlled unfolding.

Alpha-hairpin stability and folding of transmembrane segments

Molecular Dynamics (MD) simulations at low dielectric constant have been carried out for peptides matching the double spanning segments of transmembrane proteins. Different folding dynamics have been observed. The peptides folded into the stable helix-turn-helix conformation-alpha-hairpin-with antiparallel-oriented strands or unstable alpha-hairpin conformation that unfolded later into the straight helical structure. The peptide having flexible residues in the TM helices often misfolded into a tangled structure that can be avoided by restricting the flexibility of these residues. General conclusions can be drawn from the observed folding dynamics. The stability and folding of some double spanning transmembrane fragments are self-assembling. The following and/or neighboring peptide chains of the protein may support the stability of the hairpin structure of other fragments. The stability of the TM helices containing flexible residues could be maintained due to contacts with neighboring TM segments.

Native and non-native interactions along protein folding and unfolding pathways

Gō-type models, which include only native contact interactions in the energy function, are being used increasingly to describe the protein folding reaction. To investigate the validity of such models, we determine the role of native and non-native interactions along folding and unfolding pathways. For this purpose, we use a molecular mechanics energy function with an implicit solvation model (an effective energy function or potential of mean force) that can be expressed in a pairwise decomposable form. We find that for the native state and a wide range of other configurations, the contact energy is an accurate description, in part due to the cancellation of non-zero contributions from more distant residues. However, significant errors in the energy are introduced for non-native structures if the energy is calculated from the native contacts alone. Non-native contacts tend to make a significant contribution, particularly for molten globules and collapsed states along the unfolding pathways. The implication of these results for the use of Gō-type models in studies of protein folding are discussed.

Role of hydrophilic and hydrophobic contacts in folding of the second beta-hairpin fragment of protein G: molecular dynamics simulation studies of an all-atom model

Predicting the folding mechanism of the second beta-hairpin fragment of the Ig-binding domain B of streptococcal protein G is unexpectedly challenging for simplified reduced models because the models developed so far indicated a different folding mechanism from what was suggested from high-temperature unfolding and equilibrium free-energy surface analysis based on established all-atom empirical force fields in explicit or implicit solvent. This happened despite the use of empirical residue-based interactions, multibody hydrophobic interactions, and inclusions of hydrogen bonding effects in the simplified models. This article employs a recently developed all-atom (except nonpolar hydrogens) model interacting with simple square-well potentials to fold the peptide fragment by molecular dynamics simulation methods. In this study, 193 out of 200 trajectories are folded at two reduced temperatures (3.5 and 3.7) close to the transition temperature T* approximately 4.0. Each simulation takes <7 h of CPU time on a Pentium 800-MHz PC. Folding of the new all-atom model is found to be initiated by collapse before the formation of main-chain hydrogen bonds. This verifies the mechanism proposed from previous all-atom unfolding and equilibrium simulations. The new model further predicts that the collapse is initiated by two nucleation contacts (a hydrophilic contact between D46 and T49 and a hydrophobic contact between Y45 and F52), in agreement with recent NMR measurements. The results suggest that atomic packing and native contact interactions play a dominant role in folding mechanism.

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.

Effects of turn residues on beta-hairpin folding--a molecular dynamics study

Folding of beta-hairpin structures of synthetic peptides has been simulated using the molecular dynamics method with a solvent-referenced potential. Two similar sequences, Ac-MQIFVKS(D)PGKTITLKV-NH(2) and Ac-MQIFVKS(L)PGKTITLKV-NH(2), derived from the N-terminal beta-hairpin of ubiquitin, were used to study the effects of turn residues in beta-hairpin folding. The simulations were carried out for 80 ns at 297 K. With extended initial conformation, the (D)P-containing peptide folded into a stable 2:2 beta-hairpin conformation with a type II' beta-turn at (D)PG. The overall beta-hairpin ratio, calculated by the DSSP algorithm, was 32.6%. With randomly generated initial conformations, the peptide also formed the stable 2:2 beta-hairpin conformation. The interactions among the side chains in the 2:2 beta-hairpin were almost identical to those in the native protein. These interactions reduced the solvation energy upon folding and stabilized the beta-hairpin conformation. Without the solvent effect, the peptide did not fold into stable beta-hairpin structures. The solvent effect is crucial for the formation of the beta-hairpin conformation. The effect of the temperature has also been studied. The (L)P-containing peptide did not fold into a stable beta-hairpin conformation and had a much lower beta-hairpin ratio (16.6%). The( L)P-containing peptide has similar favorable side-chain interactions, but the turn formed by (L)PG does not connect well with the right-handed twist of the beta-strands. For comparison, the isolated N-terminal peptide of ubiquitin, Ac-MQIFVKTLTGKTITLEV-NH(2), was also simulated and its beta-hairpin ratio was low, indicating that the beta-hairpin in the native structure is stabilized by the interaction with the protein environment. These simulation results agreed qualitatively with the available experimental findings.