scholarly journals Full Kinetics of CO Entry, Internal Diffusion, and Exit in Myoglobin from Transition-Path Theory Simulations

2015 ◽  
Vol 137 (8) ◽  
pp. 3041-3050 ◽  
Author(s):  
Tang-Qing Yu ◽  
Mauro Lapelosa ◽  
Eric Vanden-Eijnden ◽  
Cameron F. Abrams
2016 ◽  
Vol 113 (33) ◽  
pp. 9193-9198 ◽  
Author(s):  
Yilin Meng ◽  
Diwakar Shukla ◽  
Vijay S. Pande ◽  
Benoît Roux

Nonreceptor tyrosine kinases of the Src family are large multidomain allosteric proteins that are crucial to cellular signaling pathways. In a previous study, we generated a Markov state model (MSM) to simulate the activation of c-Src catalytic domain, used as a prototypical tyrosine kinase. The long-time kinetics of transition predicted by the MSM was in agreement with experimental observations. In the present study, we apply the framework of transition path theory (TPT) to the previously constructed MSM to characterize the main features of the activation pathway. The analysis indicates that the activating transition, in which the activation loop first opens up followed by an inward rotation of the αC-helix, takes place via a dense set of intermediate microstates distributed within a fairly broad “transition tube” in a multidimensional conformational subspace connecting the two end-point conformations. Multiple microstates with negligible equilibrium probabilities carry a large transition flux associated with the activating transition, which explains why extensive conformational sampling is necessary to accurately determine the kinetics of activation. Our results suggest that the combination of MSM with TPT provides an effective framework to represent conformational transitions in complex biomolecular systems.


2020 ◽  
Vol 30 (6) ◽  
pp. 3321-3366
Author(s):  
Luzie Helfmann ◽  
Enric Ribera Borrell ◽  
Christof Schütte ◽  
Péter Koltai

Abstract Given two distinct subsets A, B in the state space of some dynamical system, transition path theory (TPT) was successfully used to describe the statistical behavior of transitions from A to B in the ergodic limit of the stationary system. We derive generalizations of TPT that remove the requirements of stationarity and of the ergodic limit and provide this powerful tool for the analysis of other dynamical scenarios: periodically forced dynamics and time-dependent finite-time systems. This is partially motivated by studying applications such as climate, ocean, and social dynamics. On simple model examples, we show how the new tools are able to deliver quantitative understanding about the statistical behavior of such systems. We also point out explicit cases where the more general dynamical regimes show different behaviors to their stationary counterparts, linking these tools directly to bifurcations in non-deterministic systems.


Author(s):  
Weikang Wang ◽  
Jianhua Xing

ABSTRACTA problem ubiquitous in almost all scientific areas is escape from a metastable state, or relaxation from one stationary distribution to a new one1. More than a century of studies lead to celebrated theoretical and computational developments such as the transition state theory and reactive flux formulation. Modern transition path sampling and transition path theory focus on an ensemble of trajectories that connect the initial and final states in a state space2, 3. However, it is generally unfeasible to experimentally observe these trajectories in multiple dimensions and compare to theoretical results. Here we report and analyze single cell trajectories of human A549 cells undergoing TGF-β induced epithelial-to-mesenchymal transition (EMT) in a combined morphology and protein texture space obtained through time lapse imaging. From the trajectories we identify parallel reaction paths with corresponding reaction coordinates and quasi-potentials. Studying cell phenotypic transition dynamics will provide testing grounds for nonequilibrium reaction rate theories.


2016 ◽  
Vol 113 (12) ◽  
pp. 3263-3268 ◽  
Author(s):  
Robert B. Best ◽  
Gerhard Hummer

All-atom molecular dynamics simulations now allow us to create movies of proteins folding and unfolding. However, it is difficult to assess the accuracy of the folding mechanisms observed because experiments cannot yet directly resolve events occurring along the transition paths between unfolded and folded states. Protein folding ϕ-values provide residue-resolved information about folding mechanisms by comparing the effects of mutations on folding rates and stability, but determining ϕ-values by separately simulating mutant proteins would be computationally demanding and prone to large statistical errors. Here we use transition path theory to develop a method for computing ϕ-values directly from the transition path ensemble, without the need for additional simulations. This path-based approach uses the full transition path information available from equilibrium folding and unfolding trajectories, or from transition path sampling, and does not require identification of folding transition states. Applying our approach to a set of simulations of 10 small proteins by Shaw and coworkers [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517–520; Piana S, Lindorff-Larsen K, Shaw DE (2011) Biophys J 100(9):L47–L49; and Piana S, Lindorff-Larsen K, Shaw DE (2013) Proc Natl Acad Sci USA 110(15):5915–5920], we find good agreement with experiments in most cases where data are available. We can further resolve the contributions to fractional ϕ-values coming from partial contact formation versus transition path heterogeneity. Although in some cases, there is substantial heterogeneity of folding mechanism, in others, such as Ubiquitin, the mechanism is strongly conserved.


2016 ◽  
Vol 195 ◽  
pp. 443-468 ◽  
Author(s):  
Ralf Banisch ◽  
Eric Vanden-Eijnden

A computational procedure is proposed to generate directly loop-erased transition paths in the context of non-equilibrium reactions, i.e. reactions that occur in systems whose dynamics is not in detailed balance. The procedure builds on results from Transition Path Theory (TPT), and it avoids altogether the need to generate reactive trajectories, either by brute-force calculations or using importance sampling schemes such as Transition Path Sampling (TPS). This is computationally advantageous since these reactive trajectories can themselves be very long and intricate in complex reactions. The loop-erased transition paths, on the other hand, are shorter and simpler because, by construction, they are pruned of all the detours typical reactive trajectories make and contain only their productive pieces that carry the effective current of the reaction. As a result they give direct access to the reaction rate and mechanism.


Sign in / Sign up

Export Citation Format

Share Document