Beta-hairpin Folding Research Articles

Understanding how Beta-Hairpins Fold using Molecular Dynamics Simulations in Multiple Force Fields

Understanding how protein substructures fold is a crucial building block for studying large, complex protein systems. The folding mechanism of beta-hairpins has been a controversial topic over the past two decades, and many contradictory mechanisms have been proposed from both experimental and computational studies. Early computational studies of beta-hairpins were performed because of computational limits on larger structures; now, increased computational power enables us to completely dissect and understand protein substructures on the atomic level with statistically robust methods. Here, we analyze an aggregated CLN025 folding simulation dataset with simulations from several different force fields to study the beta-hairpin formation mechanism in a force-field agnostic way. We use Markov state models (MSMs) to determine that the extended state first undergoes hydrophobic collapse, and from this collapsed structure the hairpin turn is formed, which is shown to be the rate-limiting step. These mechanistic conclusions as well as the timescales of the turn formation and hydrophobic collapse events are consistent with experimental results for the order of events in CLN025 folding, which allows us to see a highly detailed process that matches up with coarser experimental data. We are also able to monitor the order of hydrogen bond formation in atomic detail. This work highlights the possibility of designing MSMs to correspond to existing experimental data, to incorporate data from multiple independent simulations in different force fields, and to coarse-grain models such that they present interpretable descriptions of beta-hairpin folding that can be verified by experiment. Finally, this work lends insight into the larger protein folding problem by demonstrating a natural means for experiment and modern computational approaches to come together to better understand biophysical phenomena: namely by the coupling of advanced theory and modern force fields with experimental validation.

Modeling the mechanism of CLN025 beta-hairpin formation.

Beta-hairpins are substructures found in proteins that can lend insight into more complex systems. Furthermore, the folding of beta-hairpins is a valuable test case for benchmarking experimental and theoretical methods. Here, we simulate the folding of CLN025, a miniprotein with a beta-hairpin structure, at its experimental melting temperature using a range of state-of-the-art protein force fields. We construct Markov state models in order to examine the thermodynamics, kinetics, mechanism, and rate-determining step of folding. Mechanistically, we find the folding process is rate-limited by the formation of the turn region hydrogen bonds, which occurs following the downhill hydrophobic collapse of the extended denatured protein. These results are presented in the context of established and contradictory theories of the beta-hairpin folding process. Furthermore, our analysis suggests that the AMBER-FB15 force field, at this temperature, best describes the characteristics of the full experimental CLN025 conformational ensemble, while the AMBER ff99SB-ILDN and CHARMM22* force fields display a tendency to overstabilize the native state.

Flow-Induced Beta-Hairpin Folding of the Glycoprotein Ib-alpha Beta-Switch

Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces beta-hairpin folding of the glycoprotein Ibalpha beta-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-micros molecular dynamics simulations, we investigate the kinetics of flow-induced beta-hairpin folding. Simulations sampling different flow velocities reveal that under flow, beta-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded beta-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the beta-switch has two stable conformations imprinted on it, namely, loop and hairpin--with flow inducing a transition between the two.

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 beta-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 beta-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 beta-hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed.

The Unfolded Ensemble and Folding Mechanism of the C-Terminal GB1 β-Hairpin

In this work, we shed new light on a much-studied case of beta-hairpin folding by means of advanced molecular dynamics simulations. A fully atomistic description of the protein and the solvent molecule is used, together with metadynamics, to accelerate the sampling and estimate free-energy landscapes. This is achieved using the path collective variables approach, which provides an adaptive description of the mechanism under study. We discover that the folding mechanism is a multiscale process where the turn region conformation leads to two different energy pathways that are connected by elongated structures. The former displays a stable 2:4 native-like structure in which an optimal hydrophobic packing and hydrogen bond pattern leads to 8 kcal/mol of stabilization. The latter shows a less-structured 3:5 beta-sheet, where hydrogen bonds and hydrophobic packing provide only 2.5 kcal/mol of stability. This perspective is fully consistent with experimental evidence that shows this to be a prototypical two-state folder, while it redefines the nature of the unfolded state.

β-Hairpin Folding by a Model Amyloid Peptide in Solution and at an Interface

The development of specific agents against amyloidoses requires an understanding of the conformational distribution of fibrillogenic peptides at a microscopic level. Here, I present molecular dynamics simulations of the model amyloid peptide LSFD with sequence LSFDNSGAITIG-NH2 in explicit water and at a water/vapor interface for a total time scale of approximately 1.8 micros. An extended structure was used as initial peptide configuration. At approximately 290 K, solvated LSFD was kinetically trapped in diverse misfolded beta-sheet/coil conformations. At 350 K, in contrast, the same type II' beta-hairpin in equilibrium with less ordered but also U-shaped conformations was observed for the core residues DNSGAITI in solution and at the interface in multiple independent simulations. The most stable structural unit of the beta-hairpin was the two residue turn (GA). The core residues exhibited a well-defined folded state in which the beta-hairpin was stabilized by a hydrogen bond between the side chain of Asn-385 and the main chain carbonyl group of Gly-387. My results suggest that beta-sheet conformations indicated from previous Fourier-transform infrared spectroscopy measurements immediately after preparation of the peptide solution may not arise from protofilaments as speculated by others but are a property of LSFD monomers. In addition, combined with previous results from X-ray scattering, my findings suggest that interfacial aggregation of LSFD implies a transition from U-shaped to extended peptide conformations. This work including the first simulations of reversible beta-hairpin folding at an interface is an essential step toward a microscopic understanding of interfacial peptide folding and self-assembly. Knowledge of the main conformation of the peptide core may facilitate the design of possible inhibitors of LSFD aggregation as a test ground for future computational therapeutic strategies against amyloid diseases.

Context‐Dependence of the Contribution of Disulfide Bonds to β‐Hairpin Stability

Incorporation of disulfide bonds to stabilize protein and peptide structures is not always a successful strategy. To advance current knowledge on the contribution of disulfide bonds to beta-hairpin stability, a previously reported beta-hairpin-forming peptide was taken as a template to design a series of Cys-containing peptides. The conformational behavior of these peptides in their oxidized, disulfide-cyclized peptides, and reduced, linear peptides, was investigated on the basis of NMR parameters: NOEs, and 1H and 13C chemical shifts. We found that the effect of disulfide bonds on beta-hairpin stability depends on its location within the beta-hairpin structure, being very small or even destabilizing when connecting two hydrogen-bonded facing residues. When the disulfide bond is linking non-hydrogen-bonded facing residues, we estimated that its contribution to the free-energy change of beta-hairpin folding is approximately -1.0 kcal mol(-1). This value is larger than those reported for most beta-hairpin-stabilizing cross-strand side-chain-side-chain interactions, except for some aromatic-aromatic interactions, in particular the Trp-Trp one, and the cation-pi interaction between Trp and the non-natural methylated Arg/Lys. As disulfide bonds are frequently used to stabilize peptide conformations, our conclusions can be useful for peptide, peptidomimetic, and protein design, and may even extend to other chemical cross-links.

Absence of reptation in the high-temperature folding of the trpzip2 β-hairpin peptide

We have carried out extensive all atom explicit solvent simulations of the high-temperature folding and unfolding of the trpzip2 beta-hairpin peptide and examined the resulting trajectories for evidence of folding via a reptation mechanism. Over 300 microcanonical simulations of 10 ns each were initiated from a Boltzmann ensemble of conformations at 425 K. Though we observed numerous folding and unfolding events, no evidence of reptation was found. The diffusional dynamics of the peptide are orders of magnitude faster than any observed reptation-like motion. Our data suggest that the dominant mechanisms for beta-hairpin folding under these conditions are hydrophobic collapse and turn formation, and that rearrangements occur via significant expansion of the polypeptide chain.

Understanding the Mechanism of β-Hairpin Folding via φ-Value Analysis

The folding kinetics of a 16-residue beta-hairpin (trpzip4) and five mutants were studied by a laser-induced temperature-jump infrared method. Our results indicate that mutations which affect the strength of the hydrophobic cluster lead to a decrease in the thermal stability of the beta-hairpin, as a result of increased unfolding rates. For example, the W45Y mutant has a phi-value of approximately zero, implying a folding transition state in which the native contacts involving Trp45 are not yet formed. On the other hand, mutations in the turn or loop region mostly affect the folding rate. In particular, replacing Asp46 with Ala leads to a decrease in the folding rate by roughly 9 times. Accordingly, the phi-value for D46A is determined to be approximately 0.77, suggesting that this residue plays a key role in stabilizing the folding transition state. This is most likely due to the fact that the main chain and side chain of Asp46 form a characteristic hydrogen bond network with other residues in the turn region. Taken together, these results support the folding mechanism we proposed before, which suggests that the turn formation is the rate-limiting step in beta-hairpin folding and, consequently, a stronger turn-promoting sequence increases the stability of a beta-hairpin primarily by increasing its folding rate, whereas a stronger hydrophobic cluster increases the stability of a beta-hairpin primarily by decreasing its unfolding rate. In addition, we have examined the compactness of the thermally denatured and urea-denatured states of another 16-residue beta-hairpin, using the method of fluorescence resonance energy transfer. Our results show that the thermally denatured state of this beta-hairpin is significantly more compact than the urea-denatured state, suggesting that the very first step in beta-hairpin folding, when initiated from an extended conformation, probably corresponds to a process of hydrophobic collapse.

Is Poisson-Boltzmann theory insufficient for protein folding simulations?

The Poisson-Boltzmann theory has been widely used in the studies of energetics and conformations of biological macromolecules. Recently, introduction of the efficient generalized Born approximation has greatly extended its applicability to areas such as protein folding simulations where highly efficient computation is crucial. However, limitations have been found in the folding simulations of a well-studied beta hairpin with several generalized Born implementations and different force fields. These studies have raised the question whether the underlining Poisson-Boltzmann theory, on which the generalized Born model is calibrated, is adequate in the treatment of polar interactions for the challenging protein folding simulations. To address the question whether the Poisson-Boltzmann theory in the current formalism might be insufficient, we directly tested our efficient numerical Poisson-Boltzmann implementation in the beta-hairpin folding simulation. Good agreement between simulation and experiment was found for the beta-hairpin equilibrium structures when the numerical Poisson-Boltzmann solvent and a recently improved generalized Born solvent were used. In addition simulated thermodynamic properties also agree well with experiment in both solvents. Finally, an overall agreement on the beta-hairpin folding mechanism was found between the current and previous studies. Thus, our simulations indicate that previously observed limitations are most likely due to imperfect calibration in previous generalized Born models but not due to the limitation of the Poisson-Boltzmann theory.

Folding mechanism of β‐hairpins studied by replica exchange molecular simulations

The folding process of trpzip2 beta-hairpin is studied by the replica exchange molecular dynamics (REMD) and normal MD simulations, aiming to understand the folding mechanism of this unique small, stable, and fast folder, as well as to reveal the general principles in the folding of beta-hairpins. According to our simulations, the TS ensemble is mainly characterized by a largely formed turn and the interaction between the inner pair of hydrophobic core residues. The folding is a zipping up of hydrogen bonds. However, the nascent turn has to be stabilized by the partially formed hydrophobic core to cross the TS. Thus our folding picture is in essence a blend of hydrogen bond-centric and hydrophobic core-centric mechanism. Our simulations provide a direct evidence for a very recent experiment (Du et al., Proc Natl Acad Sci USA 2004;101:15915-15920), which suggests that the turn formation is the rate-limiting step for beta-hairpin folding and the unfolding is mainly determined by the hydrophobic interactions. Besides, the relationship between hydrogen bond stabilities and their relative importance in folding are investigated. It is found that the hydrogen bonds with higher stabilities need not play more important roles in the folding process, and vice versa.

Calorimetric Evidence for a Two-State Unfolding of the β-Hairpin Peptide Trpzip4

beta-Sheets are a common secondary structural element found in proteins. The difficulty in studying beta-sheet folding and stability is that their formation is often dependent on the tertiary structural environment within the protein. However, the discovery of water-soluble beta-hairpin peptides has allowed them to be used as model systems because they represent the smallest units of beta-sheet structure independent of tertiary structural context. Trpzip4 has been used as a model beta-hairpin peptide to study beta-hairpin folding and stability because it is highly soluble in aqueous solutions, maintains its monomeric state, and shows reversible cooperative thermal unfolding. The previously determined thermodynamic parameters for trpzip4 thermal unfolding vary depending on the spectroscopic probe used, which questions the assumption that trpzip4 unfolds in a two-state manner. Here we provide direct calorimetric evidence that the unfolding of trpzip4 follows a two-state unfolding mode. Furthermore, the thermal unfolding of trpzip4 monitored using near- and far-UV-CD yielded thermodynamic parameters similar to those determined calorimetrically, providing additional evidence for a two-state unfolding mode.

β‐Hairpin formation in aqueous solution and in the presence of trifluoroethanol: A 1H and 13C nuclear magnetic resonance conformational study of designed peptides

In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE.

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.

@-Tides as Reporters for Molecular Associations

The 1,6-dihydro-3(2H)-pyridinone unit is an amino acid surrogate that favors the extended beta-strand conformation when incorporated in an oligopeptide ("@-tide") strand. We now report that the circular dichroism (CD) signature of the vinylogous amide in the @-unit is sensitive to conformation in organic and aqueous solvents and, therefore, is useful as a quantitative measure of @-tide association and folding processes that involve this moiety. Moreover, this method can be employed in the micromolar concentration range, which is not readily accessible using other techniques. Measurements of @-tide dimerization and beta-hairpin folding equilibria not only demonstrate the utility and generality of this approach but also provide a way to quantify amino acid side chain-side chain interactions relevant to beta-sheet stability.

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial.

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.

Expected and unexpected results from combined beta-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.

Infrared study of the stability and folding kinetics of a 15-residue beta-hairpin.

The thermal stability and folding kinetics of a 15-residue beta-hairpin (SESYINPDGTWTVTE) have been studied by using infrared (IR) spectroscopy coupled with laser-induced temperature-jump (T-jump) technique for rapid folding-unfolding initiation. An alternative method based on analyzing IR difference spectra was also introduced to obtain thermodynamic properties of beta-sheets, which complements the commonly used circular dichroism (CD) and fluorescence techniques. Equilibrium IR measurements indicate that the thermal unfolding of this beta-hairpin is fairly broad. However, it can be described by a two-state transition with a thermal melting temperature of approximately 29 degrees C. Time-resolved IR measurements following a T-jump, probed at 1634 cm(-1), indicate that the folding of this beta-hairpin follows first-order kinetics and is amazingly fast. At 300 K, the folding time is approximately 0.8 micros, which is only 2-3 times slower than that of alpha-helix formation. Additionally, the energetic barrier for folding is small (approximately 2 kcal mol(-1)). These results, in conjunction with results from other studies, support a view that the details of native contacts play a dominant role in the kinetics of beta-hairpin folding.

Transition-path sampling of beta-hairpin folding.

We examine the dynamical folding pathways of the C-terminal beta-hairpin of protein G-B1 in explicit solvent at room temperature by means of a transition-path sampling algorithm. In agreement with previous free-energy calculations, the resulting path ensembles reveal a folding mechanism in which the hydrophobic residues collapse first followed by backbone hydrogen-bond formation, starting with the hydrogen bonds inside the hydrophobic core. In addition, the path ensembles contain information on the folding kinetics, including solvent motion. Using the recently developed transition interface sampling technique, we calculate the rate constant for unfolding of the protein fragment and find it to be in reasonable agreement with experiments. The results support the validation of using all-atom force fields to study protein folding.

Thermodynamic analysis of beta-hairpin-forming peptides from the thermal dependence of (1)H NMR chemical shifts.

The temperature dependence of the (1)H chemical shifts of six designed peptides previously shown to adopt beta-hairpin structures in aqueous solution has been analyzed in terms of two-state (beta-hairpin left arrow over right arrow coil) equilibrium. The stability of the beta-hairpins formed by these peptides, as derived from their T(m) (midpoint transition temperature) values, parallels in general their ability to adopt those structures as deduced from independent NMR parameters: NOEs, Deltadelta(C)(alpha)(H), Deltadelta(C)(alpha), and Deltadelta(C)(beta) values. The observed T(m) values are dependent on the particular position within the beta-hairpin that is probed, indicating that their folding to a beta-hairpin conformation deviates from a "true" two-state transition. To obtain individual T(m) values for each hairpin region in each peptide, a simplified model of a successive uncoupled two-state equilibrium covering the entire process has been applied. The distribution of T(m) values obtained for the different beta-hairpin regions (turn, strands, backbone, side chains) in the six analyzed peptides reveals a similar pattern. A model for beta-hairpin folding is proposed on the basis of this pattern and the reasonable assumption that regions showing higher T(m) values are the last ones to unfold and, presumably, the first to form. With this assumption, the analysis suggests that turn formation is the first event in beta-hairpin folding. This is consistent with previous results on the essential role of the turn sequence in beta-hairpin folding.