Fluctuations of primary ubiquitin folding intermediates in a force clamp

Folding experiments of single ubiquitin molecules under force clamp using an atomic force microscope revealed a dynamic long-lived intermediate with nanometer scale end-to-end distance fluctuations along an unexpectedly complex folding pathway. To examine the nature of this intermediate at the atomic level as well as the driving forces that give rise to the observed fluctuations, we performed molecular dynamics refolding simulations of unfolded ubiquitin under constant force. After an initial fast collapse, we find a highly dynamic, broad ensemble of conformations with partial and continuously changing secondary structure and side chain interactions. This ensemble resembles a molten-globule-like state, similar in nature to the previously described non-native state of ubiquitin in solution, but stretched by the external force. The scale of the end-to-end distance fluctuations derived from the simulations compares well with experiment. Transient formation of unspecific and metastable hydrophobic clusters along the chain are found to give rise to the observed end-to-end distance fluctuations. These distinct collapses, interpreted as folding attempts, imply an upper limit for the folding attempt frequency of approximately 10 ns. Our results suggest possible relations between force-induced unfolding and temperature or chemically induced denaturation.     

Early events in protein folding: Is there something more than hydrophobic burst?

The presence of native contacts in the denatured state of many proteins suggests that elements of the biologically active structure of these molecules are formed during the initial stage of the folding process. The rapidity with which these events take place makes it difficult to study them in vitro, but, by the same token, suitable for studies in silico. With the help of all-atom, explicit solvent, molecular dynamics simulations we have followed in time, starting from elongated structureless conformations, the early events in the folding of src-SH3 domain and of proteins G, L, and CI2. It is observed that within the first 50 ns two important events take place, essentially independent of each other: hydrophobic collapse and formation of a few selected native contacts. The same contacts are also found in simulations carried out in the presence of guanidinium chloride in order to reproduce the conditions used to characterize experimentally the denatured state and testify to the fact that these contacts are to be considered a resilient characterizing property of the denaturated state.     

Does water play a structural role in the folding of small nucleic acids?

Nucleic acid structure and dynamics are known to be closely coupled to local environmental conditions and, in particular, to the ionic character of the solvent. Here we consider what role the discrete properties of water and ions play in the collapse and folding of small nucleic acids. We study the folding of an experimentally well-characterized RNA hairpin-loop motif (sequence 5'-GGGC[GCAA]GCCU-3') via ensemble molecular dynamics simulation and, with nearly 500 micros of aggregate simulation time using an explicit representation of the ionic solvent, report successful ensemble folding simulations with a predicted folding time of 8.8(+/-2.0) micros, in agreement with experimental measurements of approximately 10 micros. Comparing our results to previous folding simulations using the GB/SA continuum solvent model shows that accounting for water-mediated interactions is necessary to accurately characterize the free energy surface and stochastic nature of folding. The formation of the secondary structure appears to be more rapid than the fastest ionic degrees of freedom, and counterions do not participate discretely in observed folding events. We find that hydrophobic collapse follows a predominantly expulsive mechanism in which a diffusion-search of early structural compaction is followed by the final formation of native structure that occurs in tandem with solvent evacuation.     

beta-hairpin stability and folding: molecular dynamics studies of the first beta-hairpin of tendamistat

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors. High temperature simulations starting from widely different initial configurations (native beta-hairpin, alpha and left-handed helical and extended conformations) begin sampling similar regions of the conformational space within tens of nanoseconds, and both native and non-native beta-hairpin conformations are recovered. Transitions between conformational clusters are accompanied by an increase in energy fluctuations, which is consistent with the increase in heat capacity measured experimentally upon protein folding. The folding events observed in the various simulations support a model for beta-hairpin formation in which the turn is formed first, followed by hydrogen bond formation closing the hairpin, and subsequent stabilization by side-chain hydrophobic interactions.     

Breaking non-native hydrophobic clusters is the rate-limiting step in the folding of an alanine-based peptide

The formation mechanism of an alanine-based peptide has been studied by all-atom molecular dynamics simulations with a recently developed all-atom point-charge force field and the Generalize Born continuum solvent model at an effective salt concentration of 0.2M. Thirty-two simulations were conducted. Each simulation was performed for 100 ns. A surprisingly complex folding process was observed. The development of the helical content can be divided into three phases with time constants of 0.06-0.08, 1.4-2.3, and 12-13 ns, respectively. Helices initiate extreme rapidly in the first phase similar to that estimated from explicit solvent simulations. Hydrophobic collapse also takes place in this phase. A folding intermediate state develops in the second phase and is unfolded to allow the peptide to reach the transition state in the third phase. The folding intermediate states are characterized by the two-turn short helices and the transition states are helix-turn-helix motifs-both of which are stabilized by hydrophobic clusters. The equilibrium helical content, calculated by both the main-chain Phi-Psi torsion angles and the main-chain hydrogen bonds, is 64-66%, which is in remarkable agreement with experiments. After corrected for the solvent viscosity effect, an extrapolated folding time of 16-20 ns is obtained that is in qualitative agreement with experiments. Contrary to the prevailing opinion, neither initiation nor growth of the helix is the rate-limiting step. Instead, the rate-limiting step for this peptide is breaking the non-native hydrophobic clusters in order to reach the transition state. The implication to the folding mechanisms of proteins is also discussed.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Kinetics of cytochrome C folding: atomically detailed simulations

The vast range of time scales (from nanoseconds to seconds) during protein folding is a challenge for experiments and computations. To make concrete predictions on folding mechanisms, atomically detailed simulations of protein folding, using potentials derived from chemical physics principles, are desired. However, due to their computational complexity, straightforward molecular dynamics simulations of protein folding are impossible today. An alternative algorithm is used that makes it possible to compute approximate atomically detailed long time trajectories (the Stochastic Difference Equation in Length). This algorithm is used to compute 26 atomically detailed folding trajectories of cytochrome c (a millisecond process). The early collapse of the protein chain (with marginal formation of secondary structure), and the earlier formation of the N and C helices (compare to the 60's helix) are consistent with the experiment. The existence of an energy barrier upon entry to the molten globule is examined as well. In addition to (favorable) comparison to experiments, we show that non-native contacts drive the formation of the molten globule. In contrast to popular folding models, the non-native contacts do not form off-pathway kinetic traps in cytochrome c.     

Specific non-native hydrophobic interactions in a hidden folding intermediate: implications for protein folding

Structures of intermediates and transition states in protein folding are usually characterized by amide hydrogen exchange and protein engineering methods and interpreted on the basis of the assumption that they have native-like conformations. We were able to stabilize and determine the high-resolution structure of a partially unfolded intermediate that exists after the rate-limiting step of a four-helix bundle protein, Rd-apocyt b(562), by multidimensional NMR methods. The intermediate has partial native-like secondary structure and backbone topology, consistent with our earlier native state hydrogen exchange results. However, non-native hydrophobic interactions exist throughout the structure. These and other results in the literature suggest that non-native hydrophobic interactions may occur generally in partially folded states. This can alter the interpretation of mutational protein engineering results in terms of native-like side chain interactions. In addition, since the intermediate exists after the rate-limiting step and Rd-apocyt b(562) folds very rapidly (k(f) approximately 10(4) s(-1)), these results suggest that non-native hydrophobic interactions, in the absence of topological misfolding, are repaired too rapidly to slow folding and cause the accumulation of folding intermediates. More generally, these results illustrate an approach for determining the high-resolution structure of folding intermediates.     

Early events in the folding of an amphipathic peptide: A multinanosecond molecular dynamics study.

Folding of the capped LQQLLQQLLQL peptide is investigated at the water-hexane interface by molecular dynamics simulations for 161.5 ns. Initially placed in the aqueous phase as a beta-strand, the peptide rapidly adsorbs to the interface, where it adopts an amphipathic conformation. The marginal presence of nonamphipathic structures throughout the complete trajectory indicates that the corresponding conformations are strongly disfavored at the interface. It is further suggestive that folding in an interfacial environment proceeds through a pathway of successive amphipathic intermediates. The energetic and entropic penalties involved in the conformational changes along this pathway markedly increase the folding time scales of LQQLLQQLLQL, explaining why the alpha-helix, the hypothesized lowest free energy structure for a sequence with a hydrophobic periodicity of 3.6, has not been reached yet. The formation of a type I beta-turn at the end of the simulation confirms the importance of such motifs as initiation sites allowing the peptide to coalesce towards a secondary structure. Proteins 1999;36:383-399.     

Free energy landscape and folding mechanism of a beta-hairpin in explicit water: a replica exchange molecular dynamics study

The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart.     

The folding pathway of ubiquitin from all-atom molecular dynamics simulations

The folding (unfolding) pathway of ubiquitin is probed using all-atom molecular dynamics simulations. We dissect the folding pathway using two techniques: first, we probe the folding pathway of ubiquitin by calculating the evolution of structural properties over time and second, we identify the rate determining transition state for folding. The structural properties that we look at are hydrophobic solvent accessible surface area (SASA) and Calpha-root-mean-square deviation (rmsd). These properties on their own tell us relatively little about the folding pathway of ubiquitin; however, when plotted against each other, they become powerful tools for dissecting ubiquitin's folding mechanism. Plots of Calpha-rmsd against SASA serve as a phase space trajectories for the folding of ubiquitin. In this study, these plots show that ubiquitin folds to the native state via the population of an intermediate state. This is shown by an initial hydrophobic collapse phase followed by a second phase of secondary structure arrangement. Analysis of the structure of the intermediate state shows that it is a collapsed species with very little secondary structure. In reconciling these observations with recent experimental data, the transition that we observe in our simulations from the unfolded state (U) to the intermediate state (I) most likely occurs in the dead-time of the stopped flow instrument. The folding pathway of ubiquitin is probed further by identification of the rate-determining transition state for folding. The method used for this is essential dynamics, which utilizes a principal component analysis (PCA) on the atomic fluctuations throughout the simulation. The five transition state structures identified in silico are in good agreement with the experimentally determined transition state. The calculation of phi-values from the structures generated in the simulations is also carried out and it shows a good correlation with the experimentally measured values. An initial analysis of the denatured state shows that it is compact with fluctuating regions of nonnative secondary structure. It is found that the compactness in the denatured state is due to the burial of some hydrophobic residues. We conclude by looking at a correlation between folding kinetics and residual structure in the denatured state. A hierarchical folding mechanism is then proposed for ubiquitin.     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.     

Dynamics of the core tryptophan during the formation of a productive molten globule intermediate of barstar

The evolution of the nanosecond dynamics of the core tryptophan, Trp53, of barstar has been monitored during the induction of collapse and structure formation in the denatured D form at pH 12, by addition of increasing concentrations of the stabilizing salt Na(2)SO(4). Time-resolved fluorescence methods have been used to monitor the dynamics of Trp53 in the intermediates that are populated during the salt-induced transition of the D form to the molten globule B form. The D form approximates a random coil and displays two rotational correlation times. A long rotational correlation time of 2.54 ns originates from segmental mobility, and a short correlation time of 0.26 ns originates from independent motion of the tryptophan side chain. Upon addition of approximately 0.1 M Na(2)SO(4), the long rotational correlation time increases to approximately 6.4 ns, as the chain collapses and the segmental motions merge to reflect the global tumbling motion of a pre-molten globule P form. The P form exists as an expanded form with approximately 30% greater volume than the native (N) state. The persistence of an approximately 50% contribution to anisotropy decay by the short rotational correlation time suggests that the core of the P form is highly molten and permits free rotation of the Trp side chain. With increasing salt concentrations, tight core packing is achieved before secondary and tertiary structure formation is complete, an observation which agrees well with earlier kinetic folding studies. Thus, the equilibrium model developed here for describing the formation of structure during folding faithfully captures snapshots of transient kinetic intermediates observed on the folding pathway of barstar. A comparison of the refolding kinetics at pH 7, when refolding is initiated from the D, P, and B forms, suggests that formation of a collapsed state with a rigid core and approximately 30% secondary and tertiary structure, which presumably defines a coarse native-like topology, constitutes the intrinsic barrier in the folding of barstar.     

Cooperative folding mechanism of a beta-hairpin peptide studied by a multicanonical replica-exchange molecular dynamics simulation

G-peptide is a 16-residue peptide of the C-terminal end of streptococcal protein G B1 domain, which is known to fold into a specific beta-hairpin within 6 micros. Here, we study molecular mechanism on the stability and folding of G-peptide by performing a multicanonical replica-exchange (MUCAREM) molecular dynamics simulation with explicit solvent. Unlike the preceding simulations of the same peptide, the simulation was started from an unfolded conformation without any experimental information on the native conformation. In the 278-ns trajectory, we observed three independent folding events. Thus MUCAREM can be estimated to accelerate the folding reaction more than 60 times than the conventional molecular dynamics simulations. The free-energy landscape of the peptide at room temperature shows that there are three essential subevents in the folding pathway to construct the native-like beta-hairpin conformation: (i) a hydrophobic collapse of the peptide occurs with the side-chain contacts between Tyr45 and Phe52, (ii) then, the native-like turn is formed accompanying with the hydrogen-bonded network around the turn region, and (iii) finally, the rest of the backbone hydrogen bonds are formed. A number of stable native hydrogen bonds are formed cooperatively during the second stage, suggesting the importance of the formation of the specific turn structure. This is also supported by the accumulation of the nonnative conformations only with the hydrophobic cluster around Tyr45 and Phe52. These simulation results are consistent with high phi-values of the turn region observed by experiment.     

Atomic-level structure characterization of an ultrafast folding mini-protein denatured state

Atomic-level analyses of non-native protein ensembles constitute an important aspect of protein folding studies to reach a more complete understanding of how proteins attain their native form exhibiting biological activity. Previously, formation of hydrophobic clusters in the 6 M urea-denatured state of an ultrafast folding mini-protein known as TC5b from both photo-CIDNP NOE transfer studies and FCS measurements was observed. Here, we elucidate the structural properties of this mini-protein denatured in 6 M urea performing (15)N NMR relaxation studies together with a thorough NOE analysis. Even though our results demonstrate that no elements of secondary structure persist in the denatured state, the heterogeneous distribution of R(2) rate constants together with observing pronounced heteronuclear NOEs along the peptide backbone reveals specific regions of urea-denatured TC5b exhibiting a high degree of structural rigidity more frequently observed for native proteins. The data are complemented with studies on two TC5b point mutants to verify the importance of hydrophobic interactions for fast folding. Our results corroborate earlier findings of a hydrophobic cluster present in urea-denatured TC5b comprising both native and non-native contacts underscoring their importance for ultra rapid folding. The data assist in finding ways of interpreting the effects of pre-existing native and/or non-native interactions on the ultrafast folding of proteins; a fact, which might have to be considered when defining the starting conditions for molecular dynamics simulation studies of protein folding.     

Aromatic and methyl NOEs highlight hydrophobic clustering in the unfolded state of an SH3 domain

The N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) exists in equilibrium between a folded (F(exch)) state and a relatively compact unfolded (U(exch)) state under nondenaturing conditions. Selectively labeled samples of the domain have been analyzed by NOESY NMR experiments to probe residual hydrophobic clustering in the U(exch) state. The labeling strategy included selective protonation of aromatic rings or delta-methyl groups on Ile and Leu residues in a highly deuterated background. Combined with long mixing times, the methods permitted observation of significant numbers of long-range interactions between hydrophobic side chains, providing evidence for multiple conformers involving non-native hydrophobic clusters around the Trp 36 indole. Comparison of these data with previously reported HN-HN NOEs yields structural insight into the diversity of structures within the U(exch) ensemble in the drkN SH3 domain. Many of the HN-HN NOEs are consistent with models containing compact residual nativelike secondary structure and greater exposure of the Trp 36 indole to solvent, similar to kinetic intermediates formed in the hierarchic condensation model of folding. However, the methyl and aromatic NOE data better fit conformations with non-native burial of the Trp indole surrounded by hydrophobic groups and more loosely formed beta-structure; these structural characteristics are more consistent with those of kinetic intermediates formed during the hydrophobic collapse mechanism of folding. This suite of NOE data provides a more complete picture of the structures that span the U(exch) state ensemble, from conformers with non-native structure but long-range contacts to those that are highly nativelike. Together, the results are also consistent with the folding funnel view involving multiple folding pathways for this molecule.     

The chicken-egg scenario of protein folding revisited

What is the first step in protein folding - hydrophobic collapse (compaction) or secondary structure formation? It is still not clear if the major driving force in protein folding is hydrogen bonding or hydrophobic interactions or both. We analyzed data on the conformational characteristics of 41 globular proteins in native and partially folded conformational states. Our analysis shows that a good correlation exists between relative decrease in hydrodynamic volume and increase in secondary structure content. No compact equilibrium intermediates lacking secondary structure, or highly ordered non-compact species, were found. This correlation provides experimental support for the hypothesis that hydrophobic collapse occurs simultaneously with formation of secondary structure in the early stages of the protein folding.     

Lack of coupling between secondary structure formation and collapse in a model polypeptide that mimics early folding intermediates, the F2 fragment of the Escherichia coli tryptophan-synthase beta chain.

The isolated, 101-residue long C-terminal (so called F2) fragment of the beta chain from Escherichia coli tryptophan synthase was shown previously to fold into an ensemble of conformations that are condensed, to contain large amounts of highly dynamic secondary structures, and to behave as a good model of structured intermediates that form at the very early stages of protein folding. Here, solvent perturbations were used to investigate the forces that are involved in stabilizing the secondary structure (monitored by far-UV CD) and the condensation of the polypeptide chain (monitored by dynamic light scattering) in isolated F2. It was observed that neither the ionic strength, nor the pH (between 7 and 10), nor salts of the Hofmeister series affected the global secondary structure contents of F2, whereas some of these salts affected the collapse slightly. Addition of trifluoroethanol resulted in a large increase in both the amount of secondary structure and the Stokes radius of F2. Conversely, F2 became more condensed upon raising the temperature from 4 to 60 degrees C, whereas in this temperature range, the secondary structure undergoes significant melting. These observations lead to the conclusion that, in isolated F2, there is no coupling between the hydrophobic collapse and the secondary structure. This finding will be discussed in terms of early events in protein folding.