The Folding Transition-State Ensemble of a Four-Helix Bundle Protein: Helix Propensity as a Determinant and Macromolecular Crowding as a Probe

The four-helix bundle protein Rd-apocyt b₅₆₂, a redesigned stable variant of apocytochrome b₅₆₂, exhibits two-state folding kinetics. Its transition-state ensemble has been characterized by Φ-value analysis. To elucidate the molecular basis of the transition-state ensemble, we have carried out high-temperature molecular dynamics simulations of the unfolding process. In six parallel simulations, unfolding started with the melting of helix I and the C-terminal half of helix IV, and followed by helix III, the N-terminal half of helix IV and helix II. This ordered melting of the helices is consistent with the conclusion from native-state hydrogen exchange, and can be rationalized by differences in intrinsic helix propensity. Guided by experimental Φ-values, a putative transition-state ensemble was extracted from the simulations. The residue helical probabilities of this transition-state ensemble show good correlation with the Φ-values. To further validate the putative transition-state ensemble, the effect of macromolecular crowding on the relative stability between the unfolded ensemble and the transition-state ensemble was calculated. The resulting effect of crowding on the folding kinetics agrees well with experimental observations. This study shows that molecular dynamics simulations combined with calculation of crowding effects provide an avenue for characterize the transition-state ensemble in atomic details.     

General dynamic properties of Abeta12-36 amyloid peptide involved in Alzheimer's disease from unfolding simulation

To study the folding/unfolding properties of a beta-amyloid peptide Abeta(12-36) of Alzheimer's disease, five molecular dynamics simulations of Abeta(12-36) in explicit water were done at 450 K starting from a structure that is stable in trifluoroethanol/water at room temperature with two alpha-helices. Due to high temperature, the initial helical structure unfolded during the simulation. The observed aspects of the unfolding were as follows. 1) One helix (helix 1) had a longer life than the other (helix 2), which correlates well with the theoretically computed Phi values. 2) Temporal prolongation of helix 1 was found before unfolding. 3) Hydrophobic cores formed frequently with rearrangement of amino-acid residues in the hydrophobic cores. The formation and rearrangement of the hydrophobic cores may be a general aspect of this peptide in the unfolded state, and the structural changes accompanied by the hydrophobic-core rearrangement may lead the peptide to the most stable structure. 4) Concerted motions (collective modes) appeared to unfold helix 1. The collective modes were similar with those observed in another simulation at 300 K. The analysis implies that the conformation moves according to the collective modes when the peptide is in the initial stage of protein unfolding and in the final stage of protein folding.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.     

Simulation and experiment at high temperatures: ultrafast folding of a thermophilic protein by nucleation-condensation

We used Phi-value analysis to characterise the transition state for folding of a thermophilic protein at the relatively high temperature of 325 K. PhiF values for the folding of the three-helix bundle, peripheral subunit binding domain from Bacillus stearothermophilus (E3BD) were determined by temperature-jump experiments in the absence of chemical denaturants. E3BD folded in microseconds through a highly diffuse transition state. Excellent agreement was observed between experiment and the results from eight (independent) molecular dynamics simulations of unfolding at 373 K. We used a combination of heteronuclear NMR experiments and molecular dynamics simulations to characterise the denatured ensemble, and found that it contained very little persistent, residual structure. However, those regions that adopt helical structure in the native state were found by simulation to be poised for helix formation in the denatured state. These regions also had significant structure in the transition state for folding. The overall folding pathway appears to be nucleation-condensation.     

Folding of a peptide corresponding to the alpha-helix in bovine pancreatic trypsin inhibitor

A short peptide corresponding to the alpha-helical region of BPTI shows partial folding in aqueous solution (pH 7) as judged by circular dichroism (CD). Folding is temperature and denaturant sensitive, and the peptide is monomeric. The difference CD spectrum, obtained from spectra at two temperatures, indicates that the peptide folds as an alpha-helix. Difference CD spectroscopy provides a sensitive assay for helix formation in peptides exhibiting small amounts of structure. Helix stability in this peptide shows a marked pH dependence which is consistent with stabilizing charged side-chain interactions with the helix dipole and/or salt bridge formation.     

Thermodynamics of the alpha-helix-coil transition of amphipathic peptides in a membrane environment: implications for the peptide-membrane binding equilibrium

Amphipathic alpha-helices are the membrane binding motif in many proteins. The corresponding peptides are often random coil in solution but are folded into an alpha-helix upon interaction with the membrane. The energetics of this ubiquitous folding process are still a matter of conjecture. Here, we present a new method to quantitatively analyze the thermodynamics of peptide folding at the membrane interface. We have systematically varied the helix content of a given amphipathic peptide when bound to the membrane and have correlated the thermodynamic binding parameters determined by isothermal titration calorimetry with the alpha-helix content obtained by circular dichroism spectroscopy. The peptides investigated were the antibiotic magainin 2 amide and three analogs in which two adjacent amino acid residues were substituted by their d-enantiomers. The thermodynamic parameters controlling the alpha-helix formation were found to be linearly related to the helicity of the membrane-bound peptides. Helix formation at the membrane surface is characterized by an enthalpy change of DeltaH(helix) approximately -0.7 kcal/mol per residue, an entropy change of DeltaS(helix) approximately -1.9 cal/molK residue and a free energy change of DeltaG(helix)=-0.14 kcal/mol residue. Helix formation is a strong driving force of peptide insertion into the membrane and accounts for about 50 % of the free energy of binding. An increase in temperature entails an unfolding of the membrane-bound helix. The temperature dependence can be described with the Zimm-Bragg theory and the enthalpy of unfolding agrees with that deduced from isothermal titration calorimetry.     

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

Nucleation and propagation of the collagen triple helix in single-chain and trimerized peptides: transition from third to first order kinetics

The kinetics of triple helix formation from single non-crosslinked peptide chains were studied for the collagen models (ProProGly)10 and (ProHypGly)10 in a broad concentration range and compared with those in nucleated trimers. At very low peptide concentrations the reaction order is 3 but decreases at higher concentrations. For (ProProGly)10 the third order rate constant is 800 M(-2) x s(-1) at 7 degrees C, which corresponds to a very long half time of 15 hours at 60 microM chain concentration. For (ProHypGly)10 the rate constant is about 1000-fold higher, which is consistent with the stabilizing effect of 4-hydroxyproline in collagens. The concentration dependence of the reaction order is explained by a nucleation mechanism in which a very unstable dimer is in fast equilibrium with the monomeric chains and addition of the third chain occurs in a rate-limiting step. At high concentrations nucleation is faster than propagation of helix formation and propagation becomes rate-limiting. To test this hypothesis an artificial nucleus was introduced by fusion of (ProProGly)10 with the trimeric foldon domain of T4 phage or the crosslinking domain of collagen III GlyProProGlyProCysCysGlyGlyGly. These domains were recombinantly attached to the C terminus of (GlyProPro)10 and link the three chains in a similar way to the C-terminal propeptide domain in collagen III. This results in a local intrinsic chain concentration of about 1 M. A first order reaction is observed for the folding of the triple helix in (GlyProPro)10foldon with a half time of 8.3 minutes, which approximately matches the rate of folding from single chains at 1 M peptide concentration. A high activation energy of 54 kJ/mol is found for this reaction, whereas the temperature dependence of the nucleation step is close to zero, confirming earlier findings on natural collagens that cis-trans isomerization of peptide bonds is the rate-limiting step in propagation.     

Full Reconstruction of a Vectorial Protein Folding Pathway by Atomic Force Microscopy and Molecular Dynamics Simulations

During co-translational folding, the nascent polypeptide chain is extruded sequentially from the ribosome exit tunnel and, under severe conformational constraints, is dictated by its one-dimensional geometry. How do such vectorial constraints impact the folding pathway? Here, we combine single-molecule atomic force spectroscopy and steered molecular dynamics simulations to examine protein folding in the presence of one-dimensional constraints that are similar to those imposed on the nascent polypeptide chain. The simulations exquisitely reproduced the experimental unfolding and refolding force extension relationships and led to the full reconstruction of the vectorial folding pathway of a large polypeptide, the 253-residue consensus ankyrin repeat protein, NI6C. We show that fully stretched and then relaxed NI6C starts folding by the formation of local secondary structures, followed by the nucleation of three N-terminal repeats. This rate-limiting step is then followed by the vectorial and sequential folding of the remaining repeats. However, after partial unfolding, when allowed to refold, the C-terminal repeats successively regain structures without any nucleation step by using the intact N-terminal repeats as a template. These results suggest a pathway for the co-translational folding of repeat proteins and have implications for mechanotransduction.     

Thermal unfolding simulations of a multimeric protein--transition state and unfolding pathways

The folding of an oligomeric protein poses an extra challenge to the folding problem because the protein not only has to fold correctly; it has to avoid nonproductive aggregation. We have carried out over 100 molecular dynamics simulations using an implicit solvation model at different temperatures to study the unfolding of one of the smallest known tetramers, p53 tetramerization domain (p53tet). We found that unfolding started with disruption of the native tetrameric hydrophobic core. The transition state for the tetramer to dimer transition was characterized as a diverse ensemble of different structures using Phi value analysis in quantitative agreement with experimental data. Despite the diversity, the ensemble was still native-like with common features such as partially exposed tetramer hydrophobic core and shifts in the dimer-dimer arrangements. After passing the transition state, the secondary and tertiary structures continued to unfold until the primary dimers broke free. The free dimer had little secondary structure left and the final free monomers were random-coil like. Both the transition states and the unfolding pathways from these trajectories were very diverse, in agreement with the new view of protein folding. The multiple simulations showed that the folding of p53tet is a mixture of the framework and nucleation-condensation mechanisms and the folding is coupled to the complex formation. We have also calculated the entropy and effective energy for the different states along the unfolding pathway and found that the tetramerization is stabilized by hydrophobic interactions.     

Kinetic analysis of molecular dynamics simulations reveals changes in the denatured state and switch of folding pathways upon single-point mutation of a beta-sheet miniprotein

The effects of a single-point mutation on folding thermodynamics and kinetics are usually interpreted by focusing on the native structure and the transition state. Here, the entire conformational spaces of a 20-residue three-stranded antiparallel beta-sheet peptide (double hairpin) and of its single-point mutant W10V are sampled close to the melting temperature by equilibrium folding-unfolding molecular dynamics simulations for a total of 40 micros. The folded state as well as the most populated free energy basins in the denatured state are isolated by grouping conformations according to fast relaxation at equilibrium. Such kinetic analysis provides more detailed and useful information than a simple projection of the free energy. The W10V mutant has the same native structure as the wild type peptide, and similar folding rate and stability. In the denatured state, the N-terminal hairpin is about 20% more structured in W10V than the wild type mainly because of van der Waals interactions. Notably, the W10V mutation influences also the van der Waals energy at the transition state ensemble causing a shift in the ratio of fluxes between two different transition state regions on parallel folding pathways corresponding to nucleation at either of the two beta-hairpins. Previous experimental studies have focused on the effects of denaturant-dependent or temperature-dependent changes in the structure of the denatured state. The atomistic simulations show that a single-point mutation in the central strand of a beta-sheet peptide results in remarkable changes in the topography of the denatured state ensemble. These changes modulate the relative accessibility of parallel folding pathways because of kinetic partitioning of the denatured state. Therefore, the observed dependence of the folding process on the starting ensemble raises questions on the biological significance of in vitro folding studies under strongly denaturing conditions.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

The unfolded state of the villin headpiece helical subdomain: computational studies of the role of locally stabilized structure

The 36 residue villin headpiece helical subdomain (HP36) is one of the fastest cooperatively folding proteins, folding on the microsecond timescale. HP36's simple three helix topology, fast folding and small size have made it an attractive model system for computational and experimental studies of protein folding. Recent experimental studies have explored the denatured state of HP36 using fragment analysis coupled with relatively low-resolution spectroscopic techniques. These studies have shown that there is apparently only a small tendency to form locally stabilized secondary structure. Here, we complement the experimental studies by using replica exchange molecular dynamics with explicit solvent to investigate the structural features of these peptide models of unfolded HP36. To ensure convergence, two sets of simulations for each fragment were performed with different initial structures, and simulations were continued until these generated very similar final ensembles. These simulations reveal low populations of native-like structure and early folding events that cannot be resolved by experiment. For each fragment, calculated J-coupling constants and helical propensities are in good agreement with experimental trends. HP-1, corresponding to residues 41 to 53 and including the first alpha-helix, contains the highest helical population. HP-3, corresponding to residues 62 through 75 and including the third alpha-helix, contains a small population of helical turn residing at the N terminus while HP-2, corresponding to residues 52 through 61 and including the second alpha-helix, formed little to no structure in isolation. Overall, HP-1 was the only fragment to adopt a native-like conformation, but the low population suggests that formation of significant structure only occurs after formation of specific tertiary interactions.     

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.     

Exploring the helix-coil transition via all-atom equilibrium ensemble simulations

The ensemble folding of two 21-residue alpha-helical peptides has been studied using all-atom simulations under several variants of the AMBER potential in explicit solvent using a global distributed computing network. Our extensive sampling, orders of magnitude greater than the experimental folding time, results in complete convergence to ensemble equilibrium. This allows for a quantitative assessment of these potentials, including a new variant of the AMBER-99 force field, denoted AMBER-99 phi, which shows improved agreement with experimental kinetic and thermodynamic measurements. From bulk analysis of the simulated AMBER-99 phi equilibrium, we find that the folding landscape is pseudo-two-state, with complexity arising from the broad, shallow character of the "native" and "unfolded" regions of the phase space. Each of these macrostates allows for configurational diffusion among a diverse ensemble of conformational microstates with greatly varying helical content and molecular size. Indeed, the observed structural dynamics are better represented as a conformational diffusion than as a simple exponential process, and equilibrium transition rates spanning several orders of magnitude are reported. After multiple nucleation steps, on average, helix formation proceeds via a kinetic "alignment" phase in which two or more short, low-entropy helical segments form a more ideal, single-helix structure.     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

Molecular dynamics simulations of a fibrillogenic peptide derived from apolipoprotein C-II

The pathway to amyloid fibril formation in proteins involves specific structural changes leading to the combination of misfolded intermediates into oligomeric assemblies. Recent NMR studies showed the presence of “turns” in amyloid peptides, indicating that turn formation may play an important role in the nucleation of the intramolecular folding and possible assembly of amyloid. Fully solvated all-atom molecular dynamics simulations were used to study the structure and dynamics of the apolipoprotein C-II peptide 56 to 76, associated with the formation of amyloid fibrils. The peptide populated an ensemble of turn structures, stabilized by hydrogen bonds and hydrophobic interactions enabling the formation of a strong hydrophobic core which may provide the conditions required to initiate aggregation. Two competing mechanisms discussed in the literature were observed. This has implications in understanding the mechanism of amyloid formation in not only apoC-II and its fragments, but also in other amyloidogenic peptides.