Direct molecular dynamics observation of protein folding transition state ensemble.
Ontology highlight
ABSTRACT:
SUBMITTER: Ding F
PROVIDER: S-EPMC1302427 | biostudies-other | 2002 Dec
REPOSITORIES: biostudies-other
Publications
Direct molecular dynamics observation of protein folding transition state ensemble.
Ding Feng F Dokholyan Nikolay V NV Buldyrev Sergey V SV Stanley H Eugene HE Shakhnovich Eugene I EI
Biophysical journal 20021201 6
The concept of the protein transition state ensemble (TSE), a collection of the conformations that have 50% probability to convert rapidly to the folded state and 50% chance to rapidly unfold, constitutes the basis of the modern interpretation of protein engineering experiments. It has been conjectured that conformations constituting the TSE in many proteins are the expanded and distorted forms of the native state built around a specific folding nucleus. This view has been supported by a number ...[more]
Similar Datasets
Project description:Single-molecule fluorescence resonance energy transfer (smFRET) experiments are extremely useful in studying protein folding but are generally limited to time scales of greater than approximately 100 micros and distances greater than approximately 2 nm. We used single-molecule fluorescence quenching by photoinduced electron transfer, detecting short-range events, in combination with fluorescence correlation spectroscopy (PET-FCS) to investigate folding dynamics of the small binding domain BBL with nanosecond time resolution. The kinetics of folding appeared as a 10-micros decay in the autocorrelation function, resulting from stochastic fluctuations between denatured and native conformations of individual molecules. The observed rate constants were probe independent and in excellent agreement with values derived from conventional temperature-jump (T-jump) measurements. A submicrosecond relaxation was detected in PET-FCS data that reported on the kinetics of intrachain contact formation within the thermally denatured state. We engineered a mutant of BBL that was denatured under the reaction conditions that favored folding of the parent wild type ("D(phys)"). D(phys) had the same kinetic signature as the thermally denatured state and revealed segmental diffusion with a time constant of intrachain contact formation of 500 ns. This time constant was more than 10 times faster than folding and in the range estimated to be the "speed limit" of folding. D(phys) exhibited significant deviations from a random coil. The solvent viscosity and temperature dependence of intrachain diffusion showed that chain motions were slaved by the presence of intramolecular interactions. PET-FCS in combination with protein engineering is a powerful approach to study the early events and mechanism of ultrafast protein folding.
| S-EPMC2773960 | biostudies-literature
Project description:Characterizing the conformations of protein in the transition state ensemble (TSE) is important for studying protein folding. A promising approach pioneered by Vendruscolo et al. [Nature (London) 409, 641 (2001)] to study TSE is to generate conformations that satisfy all constraints imposed by the experimentally measured φ values that provide information about the native likeness of the transition states. Faísca et al. [J. Chem. Phys. 129, 095108 (2008)] generated conformations of TSE based on the criterion that, starting from a TS conformation, the probabilities of folding and unfolding are about equal through Markov Chain Monte Carlo (MCMC) simulations. In this study, we use the technique of constrained sequential Monte Carlo method [Lin et al., J. Chem. Phys. 129, 094101 (2008); Zhang et al. Proteins 66, 61 (2007)] to generate TSE conformations of acylphosphatase of 98 residues that satisfy the φ-value constraints, as well as the criterion that each conformation has a folding probability of 0.5 by Monte Carlo simulations. We adopt a two stage process and first generate 5000 contact maps satisfying the φ-value constraints. Each contact map is then used to generate 1000 properly weighted conformations. After clustering similar conformations, we obtain a set of properly weighted samples of 4185 candidate clusters. Representative conformation of each of these cluster is then selected and 50 runs of Markov chain Monte Carlo (MCMC) simulation are carried using a regrowth move set. We then select a subset of 1501 conformations that have equal probabilities to fold and to unfold as the set of TSE. These 1501 samples characterize well the distribution of transition state ensemble conformations of acylphosphatase. Compared with previous studies, our approach can access much wider conformational space and can objectively generate conformations that satisfy the φ-value constraints and the criterion of 0.5 folding probability without bias. In contrast to previous studies, our results show that transition state conformations are very diverse and are far from nativelike when measured in cartesian root-mean-square deviation (cRMSD): the average cRMSD between TSE conformations and the native structure is 9.4 Å for this short protein, instead of 6 Å reported in previous studies. In addition, we found that the average fraction of native contacts in the TSE is 0.37, with enrichment in native-like β-sheets and a shortage of long range contacts, suggesting such contacts form at a later stage of folding. We further calculate the first passage time of folding of TSE conformations through calculation of physical time associated with the regrowth moves in MCMC simulation through mapping such moves to a Markovian state model, whose transition time was obtained by Langevin dynamics simulations. Our results indicate that despite the large structural diversity of the TSE, they are characterized by similar folding time. Our approach is general and can be used to study TSE in other macromolecules.
| S-EPMC3071304 | biostudies-other
Project description:Transcription activator-like effector (TALE) proteins are a class of programmable DNA-binding proteins for which the fundamental mechanisms governing the search process are not fully understood. Here we use single-molecule techniques to directly observe TALE search dynamics along DNA templates. We find that TALE proteins are capable of rapid diffusion along DNA using a combination of sliding and hopping behaviour, which suggests that the TALE search process is governed in part by facilitated diffusion. We also observe that TALE proteins exhibit two distinct modes of action during the search process-a search state and a recognition state-facilitated by different subdomains in monomeric TALE proteins. Using TALE truncation mutants, we further demonstrate that the N-terminal region of TALEs is required for the initial non-specific binding and subsequent rapid search along DNA, whereas the central repeat domain is required for transitioning into the site-specific recognition state.
| S-EPMC4458887 | biostudies-literature
Project description:Interleukin-33 (IL-33) is currently the focus of multiple investigations into targeting pernicious inflammatory disorders. This mediator of inflammation plays a prevalent role in chronic disorders such as asthma, rheumatoid arthritis, and progressive heart disease. In order to better understand the possible link between the folding free energy landscape and functional regions in IL-33, a combined experimental and theoretical approach was applied. IL-33 is a pseudo- symmetrical protein composed of three distinct structural elements that complicate the folding mechanism due to competition for nucleation on the dominant folding route. Trefoil 1 constitutes the majority of the binding interface with the receptor whereas Trefoils 2 and 3 provide the stable scaffold to anchor Trefoil 1. We identified that IL-33 folds with a three-state mechanism, leading to a rollover in the refolding arm of its chevron plots in strongly native conditions. In addition, there is a second slower refolding phase that exhibits the same rollover suggesting similar limitations in folding along parallel routes. Characterization of the intermediate state and the rate limiting steps required for folding suggests that the rollover is attributable to a moving transition state, shifting from a post- to pre-intermediate transition state as you move from strongly native conditions to the midpoint of the transition. On a structural level, we found that initially, all independent Trefoil units fold equally well until a QCA of 0.35 when Trefoil 1 will backtrack in order to allow Trefoils 2 and 3 to fold in the intermediate state, creating a stable scaffold for Trefoil 1 to fold onto during the final folding transition. The formation of this intermediate state and subsequent moving transition state is a result of balancing the difficulty in folding the functionally important Trefoil 1 onto the remainder of the protein. Taken together our results indicate that the functional element of the protein is geometrically frustrated, requiring the more stable elements to fold first, acting as a scaffold for docking of the functional element to allow productive folding to the native state.
| S-EPMC4667907 | biostudies-literature
Project description:Statistical theories of protein folding have long predicted plausible mechanisms for reducing the vast conformational space through distinct ensembles of structures. However, these predictions have remained untested by bulk techniques, because the conformational diversity of folding molecules has been experimentally unapproachable. Owing to recent advances in single molecule force-clamp spectroscopy, we are now able to probe the structure and dynamics of the small protein ubiquitin by measuring its length and mechanical stability during each stage of folding. Here, we discover that upon hydrophobic collapse, the protein rapidly selects a subset of minimum energy structures that are mechanically weak and essential precursors of the native fold. From this much reduced ensemble, the native state is acquired through a barrier-limited transition. Our results support the validity of statistical mechanics models in describing the folding of a small protein on biological timescales.
| S-EPMC2705594 | biostudies-literature
Project description:The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-alpha, all-beta, and mixed alpha/beta and alpha +beta) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an alpha-helix maintained more structure in the TS ensemble than did residues starting in beta-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.
| S-EPMC2784564 | biostudies-literature
Project description:Small globular proteins and peptides commonly exhibit two-state folding kinetics in which the rate limiting step of folding is the surmounting of a single free energy barrier at the transition state (TS) separating the folded and the unfolded states. An intriguing question is whether the polypeptide chain reaches, and leaves, the TS by completely random fluctuations, or whether there is a directed, stepwise process. Here, the folding TS of a 15-residue β-hairpin peptide, Peptide 1, is characterized using independent 2.5 μs-long unbiased atomistic molecular dynamics (MD) simulations (a total of 15 μs). The trajectories were started from fully unfolded structures. Multiple (spontaneous) folding events to the NMR-derived conformation are observed, allowing both structural and dynamical characterization of the folding TS. A common loop-like topology is observed in all the TS structures with native end-to-end and turn contacts, while the central segments of the strands are not in contact. Non-native sidechain contacts are present in the TS between the only tryptophan (W11) and the turn region (P7-G9). Prior to the TS the turn is found to be already locked by the W11 sidechain, while the ends are apart. Once the ends have also come into contact, the TS is reached. Finally, along the reactive folding paths the cooperative loss of the W11 non-native contacts and the formation of the central inter-strand native contacts lead to the peptide rapidly proceeding from the TS to the native state. The present results indicate a directed stepwise process to folding the peptide.
| S-EPMC3169518 | biostudies-literature
Project description:The distributed computing (DC) paradigm in conjunction with the folding@home (FH) client server has been used to study the folding kinetics of small peptides and proteins, giving excellent agreement with experimentally measured folding rates, although pathways sampled in these simulations are not always consistent with the folding mechanism. In this study, we use a coarse-grain model of protein L, whose two-state kinetics have been characterized in detail by using long-time equilibrium simulations, to rigorously test a FH protocol using approximately 10,000 short-time, uncoupled folding simulations starting from an extended state of the protein. We show that the FH results give non-Poisson distributions and early folding events that are unphysical, whereas longer folding events experience a correct barrier to folding but are not representative of the equilibrium folding ensemble. Using short-time, uncoupled folding simulations started from an equilibrated denatured state ensemble (DSE), we also do not get agreement with the equilibrium two-state kinetics because of overrepresented folding events arising from higher energy subpopulations in the DSE. The DC approach using uncoupled short trajectories can make contact with traditionally measured experimental rates and folding mechanism when starting from an equilibrated DSE, when the simulation time is long enough to sample the lowest energy states of the unfolded basin and the simulated free-energy surface is correct. However, the DC paradigm, together with faster time-resolved and single-molecule experiments, can also reveal the breakdown in the two-state approximation due to observation of folding events from higher energy subpopulations in the DSE.
| S-EPMC1283817 | biostudies-literature
Project description:Both folded and unfolded conformations should be observed for a protein at its melting temperature (T(m)), where DeltaG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental T(m), the protein rapidly loses its low-temperature native structure; it then unfolds before refolding to a stable, native-like conformation. The initial unfolding follows the unfolding pathway described previously for higher-temperature simulations: the hydrophobic core is disrupted, the beta-sheet pulls apart and the alpha-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a non-native hydrophobic cluster. Refolding simply reverses this path, the side-chain interactions shift, the helix refolds, and the native packing and hydrogen bonds are recovered. The end result of this refolding is not the initial crystal structure; it contains the proper topology and the majority of the native contacts, but the structure is expanded and the contacts are long. We believe this to be the native state at elevated temperature, and the change in volume and contact lengths is consistent with experimental studies of other native proteins at elevated temperature and the chemical denaturant equivalent of T(m).
| S-EPMC1885941 | biostudies-literature
Project description:Biology relies on the precise self-assembly of its molecular components. Generic principles of protein folding have emerged from extensive studies on small, water-soluble proteins, but it is unclear how these ideas are translated into more complex situations. In particular, the one-third of cellular proteins that reside in biological membranes will not fold like water-soluble proteins because membrane proteins need to expose, not hide, their hydrophobic surfaces. Here, we apply the powerful protein engineering method of Phi-value analysis to investigate the folding transition state of the alpha-helical membrane protein, bacteriorhodopsin, from a partially unfolded state. Our results imply that much of helix B of the seven-transmembrane helical protein is structured in the transition state with single-point alanine mutations in helix B giving Phi values >0.8. However, residues Y43 and T46 give lower Phi values of 0.3 and 0.5, respectively, suggesting a possible reduction in native structure in this region of the helix. Destabilizing mutations also increase the activation energy of folding, which is accompanied by an apparent movement of the transition state toward the partially unfolded state. This apparent transition state movement is most likely due to destabilization of the structured, unfolded state. These results contrast with the Hammond effect seen for several water-soluble proteins in which destabilizing mutations cause the transition state to move toward, and become closer in energy to, the folded state. We thus introduce a classic folding analysis method to membrane proteins, providing critical insight into the folding transition state.
| S-EPMC2630066 | biostudies-literature