Constructing the equilibrium ensemble of folding pathways from short off-equilibrium simulations

Abstract

Characterizing the equilibrium ensemble of folding pathways, including their relative probability, is one of the major challenges in protein folding theory today. Although this information is in principle accessible via all-atom molecular dynamics simulations, it is difficult to compute in practice because protein folding is a rare event and the affordable simulation length is typically not sufficient to observe an appreciable number of folding events, unless very simplified protein models are used. Here we present an approach that allows for the reconstruction of the full ensemble of folding pathways from simulations that are much shorter than the folding time. This approach can be applied to all-atom protein simulations in explicit solvent. It does not use a predefined reaction coordinate but is based on partitioning the state space into small conformational states and constructing a Markov model between them. A theory is presented that allows for the extraction of the full ensemble of transition pathways from the unfolded to the folded configurations. The approach is applied to the folding of a PinWW domain in explicit solvent where the folding time is two orders of magnitude larger than the length of individual simulations. The results are in good agreement with kinetic experimental data and give detailed insights about the nature of the folding process which is shown to be surprisingly complex and parallel. The analysis reveals the existence of misfolded trap states outside the network of efficient folding intermediates that significantly reduce the folding speed. D iscovering the mechanism by which proteins fold into their native 3D structure remains an intriguing problem (1, 2). Essential questions are: How does an ensemble of denatured molecules find the same native structure, starting from different conformations? Are there particular sequences in which the structural elements of a protein are formed (3-5)? Are there multiple parallel routes by which protein structure formation can proceed (6, 7)?Full answers to these questions require one to characterize the ensemble of folding pathways, including their relative probabilities. In principle, this detailed information is accessible via molecular dynamics (MD) simulations which, when used in concert with experimental evidence, are becoming an increasingly accepted tool to understanding structural details that are not easily accessible via the experimental observables (8). MD simulations with atomistic models of proteins have been used to study the dynamics of small proteins with folding times in the microsecond range (9-13). However, even though MD simulations make the full spatiotemporal detail accessible to observation, the characterization of the pathway ensemble is computationally difficult: A brute-force approach would start simulations from an equilibrium of unfolded structures, say A, and simulate until they relax into a set of folded state B. The analysis would then only be comprised of those trajectory segments that leave A and relax to B without...

Optimal control of molecular dynamics using Markov state models

“…Recently the interest in MSMs has increased a lot, for it had been demonstrated that MSMs can be constructed even for very high dimensional systems [25]. They have been especially useful for modelling the interesting slow dynamics of biomolecules [21,[28][29][30][31][32] and materials [33] (there under the name "kinetic Monte Carlo"). If the system exhibits metastability and the jump process between the metastable sets are approximately Markovian, the corresponding MSM simply describes the Markov process that jumps between the sets with the aggregated statistics of the original process.…”

“…These transitions span large ranges of length scales, time scales and complexity, and include processes as important as folding [15,16], complex conformational rearrangements between native protein substates [17,18], and ligand binding [19]. MD simulations are becoming increasingly accepted as a tool to investigate both the structural and the dynamical features of these transitions at a level of detail that is beyond that accessible in laboratory experiments [20][21][22].…”

A network of molecular switches controls the activation of the two-component response regulator NtrC

“…MSM-weighted probability distributions, both one-and two-dimensional, are obtained by binning the raw data within each MSM state and weighting it by the MSM equilibrium state population. High-flux pathways between active and inactive states and MSM state committor values are also calculated from T ij using the transition matrix definition of transition path theory 57,58 . Committor values measure the probability of reaching the active state (at values near one) before returning to the inactive state (at values near zero).…”

Assessment and Optimization of Collective Variables for Protein Conformational Landscape: GB1 β-hairpin as a Case Study

“…Interestingly, these four metastable states obtained from hidden MSM correspond to the energy minima basin (ms1-4) obtained using REMD data. Furthermore, we used transition path theory (TPT) [53][54][55] as implemented in PyEMMA 56 package to identify the folding paths from unfolded state. The TPT analysis suggests that transition from unfolded state (ms1) to folded state (ms4) is more preferable through intermediate ms3 than ms2 (helical state).…”