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.     

Detection of a hidden folding intermediate in the focal adhesion target domain: Implications for its function and folding

The focal adhesion target (FAT) domain of focal adhesion kinase has a four-helix bundle structure. Based on a hydrogen exchange-constrained computer simulation study and some indirect experimental results, it has been suggested that a partially unfolded state of the FAT domain with the N-terminal helix unfolded plays an important role in its biological function. Here, using a native-state hydrogen exchange method, we directly detected an intermediate with the N-terminal helix unfolded in a mutant (Y925E) of the FAT domain. In addition, kinetic folding studies on the FAT domain suggest that this intermediate exists on the native side of the rate-limiting transition state for folding. These results provide more direct evidence of the existence of the proposed intermediate and help to understand the folding mechanism of small single domain proteins.     

Characterisation of the conformational properties of urea-unfolded Im7: implications for the early stages of protein folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.     

Consequences of stabilizing the natively disordered f helix for the folding pathway of apomyoglobin

The F helix region of sperm whale apomyoglobin is disordered, undergoing conformational fluctuations between a folded helical conformation and one or more locally unfolded states. To examine the effects of F helix stabilization on the folding pathway of apomyoglobin, we have introduced mutations to augment intrinsic helical structure in the F helix of the kinetic folding intermediate and to increase its propensity to fold early in the pathway, using predictions based on plots of the average area buried upon folding (AABUF) derived from the primary sequence. Two mutant proteins were prepared: a double mutant, P88K/S92K (F2), and a quadruple mutant, P88K/A90L/S92K/A94L (F4). Whereas the AABUF for F2 predicts that the F helix will not fold early in the pathway, the F helix in F4 shows a significantly increased AABUF and is therefore predicted to fold early. Protection of amide protons by formation of hydrogen-bonded helical structure during the early folding events has been analyzed by pH-pulse labeling. Consistent with the AABUF prediction, many of the F helix residues for F4 are significantly protected in the kinetic intermediate but are not protected in the F2 mutant. F4 folds via a kinetically trapped burst-phase intermediate that contains stabilized secondary structure in the A, B, F, G, and H helix regions. Rapid folding of the F helix stabilizes the central core of the misfolded intermediate and inhibits translocation of the H helix back to its native position, thereby decreasing the overall folding rate.     

All-atom Monte Carlo simulation of GCAA RNA folding

We report a detailed all-atom simulation of the folding of the GCAA RNA tetraloop. The GCAA tetraloop motif is a very common and thermodynamically stable secondary structure in natural RNAs. We use our simulation methods to study the folding behavior of a 12-base GCAA tetraloop structure with a four-base helix adjacent to the tetraloop proper. We implement an all-atom Monte Carlo (MC) simulation of RNA structural dynamics using a Go potential. Molecular dynamics (MD) simulation of RNA and protein has realistic energetics and sterics, but is extremely expensive in terms of computational time. By coarsely treating non-covalent energetics, but retaining all-atom sterics and entropic effects, all-atom MC techniques are a useful method for the study of protein and now RNA. We observe a sharp folding transition for this structure, and in simulations at room temperature the state histogram shows three distinct minima: an unfolded state (U), a more narrow intermediated state (I), and a narrow folded state (F). The intermediate consists primarily of structures with the GCAA loop and some helix hydrogen bonds formed. Repeated kinetic folding simulations reveal that the number of helix base-pairs forms a simple 1D reaction coordinate for the I-->N transition.     

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 μs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 ? was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.     

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.     

Microsecond Scale Replica Exchange Molecular Dynamic Simulation of Villin Headpiece: An Insight into the Folding Landscape

Abstract

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 μs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

The folding thermodynamics and kinetics of crambin using an all-atom Monte Carlo simulation

We present a novel Monte Carlo simulation of protein folding, in which all heavy atoms are represented as interacting hard spheres. This model includes all degrees of freedom relevant to folding, all side-chain and backbone torsions, and uses a Go potential. In this study, we focus on the 46 residue alpha/beta protein crambin and two of its structural components, the helix and helix hairpin. For a wide range of temperatures, we recorded multiple folding events of these three structures from random coils to native conformations that differ by less than 1 A C(alpha) dRMS from their crystal structure coordinates. The thermodynamics and kinetic mechanism of the helix-coil transition obtained from our simulation shows excellent agreement with currently available experimental and molecular dynamics data. Based on insights obtained from folding its smaller structural components, a possible folding mechanism for crambin is proposed. We observed that the folding occurs via a cooperative, first order-like process, and that many folding pathways to the native state exist. One particular sequence of events constitutes a "fast-folding" pathway where kinetic traps are avoided. At very low temperatures, a kinetic trap arising from the incorrect packing of side-chains was observed. These results demonstrate that folding to the native state can be observed in a reasonable amount of time on desktop computers even when an all-atom representation is used, provided the energetics sufficiently stabilize the native state.     

De novo simulations of the folding thermodynamics of the GCN4 leucine zipper.

Entropy Sampling Monte Carlo (ESMC) simulations were carried out to study the thermodynamics of the folding transition in the GCN4 leucine zipper (GCN4-lz) in the context of a reduced model. Using the calculated partition functions for the monomer and dimer, and taking into account the equilibrium between the monomer and dimer, the average helix content of the GCN4-lz was computed over a range of temperatures and chain concentrations. The predicted helix contents for the native and denatured states of GCN4-lz agree with the experimental values. Similar to experimental results, our helix content versus temperature curves show a small linear decline in helix content with an increase in temperature in the native region. This is followed by a sharp transition to the denatured state. van't Hoff analysis of the helix content versus temperature curves indicates that the folding transition can be described using a two-state model. This indicates that knowledge-based potentials can be used to describe the properties of the folded and unfolded states of proteins.     

Simulation of folding of a small alpha-helical protein in atomistic detail using worldwide-distributed computing

By employing thousands of PCs and new worldwide-distributed computing techniques, we have simulated in atomistic detail the folding of a fast-folding 36-residue alpha-helical protein from the villin headpiece. The total simulated time exceeds 300 micros, orders of magnitude more than previous simulations of a molecule of this size. Starting from an extended state, we obtained an ensemble of folded structures, which is on average 1.7A and 1.9A away from the native state in C(alpha) distance-based root-mean-square deviation (dRMS) and C(beta) dRMS sense, respectively. The folding mechanism of villin is most consistent with the hydrophobic collapse view of folding: the molecule collapses non-specifically very quickly ( approximately 20ns), which greatly reduces the size of the conformational space that needs to be explored in search of the native state. The conformational search in the collapsed state appears to be rate-limited by the formation of the aromatic core: in a significant fraction of our simulations, the C-terminal phenylalanine residue packs improperly with the rest of the hydrophobic core. We suggest that the breaking of this interaction may be the rate-determining step in the course of folding. On the basis of our simulations we estimate the folding rate of villin to be approximately 5micros. By analyzing the average features of the folded ensemble obtained by simulation, we see that the mean folded structure is more similar to the native fold than any individual folded structure. This finding highlights the need for simulating ensembles of molecules and averaging the results in an experiment-like fashion if meaningful comparison between simulation and experiment is to be attempted. Moreover, our results demonstrate that (1) the computational methodology exists to simulate the multi-microsecond regime using distributed computing and (2) that potential sets used to describe interatomic interactions may be sufficiently accurate to reach the folded state, at least for small proteins. We conclude with a comparison between our results and current protein-folding theory.     

Core formation in apomyoglobin: probing the upper reaches of the folding energy landscape

An acid-destabilized form of apomyoglobin, the so-called E state, consists of a set of heterogeneous structures that are all characterized by a stable hydrophobic core composed of 30-40 residues at the intersection of the A, G, and H helices of the protein, with little other secondary structure and no other tertiary structure. Relaxation kinetics studies were carried out to characterize the dynamics of core melting and formation in this protein. The unfolding and/or refolding response is induced by a laser-induced temperature jump between the folded and unfolded forms of E, and structural changes are monitored using the infrared amide I' absorbance at 1648-1651 cm(-1) that reports on the formation of solvent-protected, native-like helix in the core and by fluorescence emission changes from apomyoglobin's Trp14, a measure of burial of the indole group of this residue. The fluorescence kinetics data are monoexponential with a relaxation time of 14 micros. However, infrared kinetics data are best fit to a biexponential function with relaxation times of 14 and 59 micros. These relaxation times are very fast, close to the limits placed on folding reactions by diffusion. The 14 micros relaxation time is weakly temperature dependent and thus represents a pathway that is energetically downhill. The appearance of this relaxation time in both the fluorescence and infrared measurements indicates that this folding event proceeds by a concomitant formation of compact secondary and tertiary structures. The 59 micros relaxation time is much more strongly temperature dependent and has no fluorescence counterpart, indicating an activated process with a large energy barrier wherein nonspecific hydrophobic interactions between helix A and the G and H helices cause some helix burial but Trp14 remains solvent exposed. These results are best fit by a multiple-pathway kinetic model when U collapses to form the various folded core structures of E. Thus, the results suggest very robust dynamics for core formation involving multiple folding pathways and provide significant insight into the primary processes of protein folding.     

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.     

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.     

Enhanced Sampling of Coarse-Grained Transmembrane-Peptide Structure Formation from Hydrogen-Bond Replica Exchange

Protein structure formation in the membrane highlights a grand challenge of sampling in computer simulations, because kinetic traps and slow dynamics make it difficult to find the native state. Exploiting increased fluctuations at higher temperatures can help overcome free-energy barriers, provided the membrane’s structure remains stable. In this work, we apply Hamiltonian replica-exchange molecular dynamics, where we only tune the backbone hydrogen-bond strength to help reduce the propensity of long-lived misfolded states. Using a recently developed coarse-grained model, we illustrate the robustness of the method by folding different WALP transmembrane helical peptides starting from stretched, unstructured conformations. We show the efficiency of the method by comparing to simulations without enhanced sampling, achieving folding in one example after significantly longer simulation times. Analysis of the bilayer structure during folding provides insight into the local membrane deformation during helix formation as a function of chain length (from 16 to 23 residues). Finally, we apply our method to fold the 50-residue-long major pVIII coat protein (fd coat) of the filamentous fd bacteriophage. Our results agree well with experimental structures and atomistic simulations based on implicit membrane models, suggesting that our explicit CG folding protocol can serve as a starting point for better-refined atomistic simulations in a multiscale framework.     

Exploring the Minimally Frustrated Energy Landscape of Unfolded ACBP

The unfolded state of globular proteins is not well described by a simple statistical coil due to residual structural features, such as secondary structure or transiently formed long-range contacts. The principle of minimal frustration predicts that the unfolded ensemble is biased toward productive regions in the conformational space determined by the native structure. Transient long-range contacts, both native-like and non-native-like, have previously been shown to be present in the unfolded state of the four-helix-bundle protein acyl co-enzyme binding protein (ACBP) as seen from both perturbations in nuclear magnetic resonance (NMR) chemical shifts and structural ensembles generated from NMR paramagnetic relaxation data. To study the nature of the contacts in detail, we used paramagnetic NMR relaxation enhancements, in combination with single-point mutations, to obtain distance constraints for the acid-unfolded ensemble of ACBP. We show that, even in the acid-unfolded state, long-range contacts are specific in nature and single-point mutations affect the free-energy landscape of the unfolded protein. Using this approach, we were able to map out concerted, interconnected, and productive long-range contacts. The correlation between the native-state stability and compactness of the denatured state provides further evidence for native-like contact formation in the denatured state. Overall, these results imply that, even in the earliest stages of folding, ACBP dynamics are governed by native-like contacts on a minimally frustrated energy landscape.     

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.