1
|
Harris J, Chipot C, Roux B. Statistical Mechanical Theories of Membrane Permeability. J Phys Chem B 2024; 128:9183-9196. [PMID: 39283709 DOI: 10.1021/acs.jpcb.4c05020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2024]
Abstract
A popular theoretical framework to compute the permeability coefficient of a molecule is provided by the classic Smoluchowski-Kramers treatment of the steady-state diffusive flux across a free-energy barrier. Within this framework, commonly termed "inhomogeneous solubility-diffusion" (ISD), the permeability, P, is expressed in closed form in terms of the potential of mean force and position-dependent diffusivity of the molecule of interest along the membrane normal. In principle, both quantities can be calculated from all-atom MD simulations. Although several methods exist for calculating the position-dependent diffusivity, each of these is at best an estimate. In addition, the ISD model does not account for memory effects along the chosen reaction coordinate. For these reasons, it is important to seek alternative theoretical formulations to determine the permeability coefficient that are able to account for the factors ignored by the ISD approximation. Using Green-Kubo linear response theory, we establish the familiar constitutive relation between the flux density across the membrane and the difference in the concentration of a permeant molecule, j = PΔC. On this basis, we derive a time-correlation function expression for the nonequilibrium flux across a membrane that is reminiscent of the transmission coefficient in the reactive flux formalism treatment of transition rates. An analysis based on the transition path theory framework is exploited to derive alternative expressions for the permeability coefficient. The different strategies are illustrated with stochastic simulations based on the generalized Langevin equation in addition to unbiased molecular dynamics simulations of water permeation of a lipid bilayer.
Collapse
Affiliation(s)
- Jonathan Harris
- Department of Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Université de Lorraine, Unité Mixte de Recherche n7019, B.P. 70239, 54506 cedex Vandœuvre-lès-Nancy, France
- Theoretical and Computational Biophysics Group, Beckman Institute, and Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Biochemistry and Molecular Biology, Department of Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Department of Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
2
|
Kang P, Trizio E, Parrinello M. Computing the committor with the committor to study the transition state ensemble. NATURE COMPUTATIONAL SCIENCE 2024; 4:451-460. [PMID: 38839932 DOI: 10.1038/s43588-024-00645-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2024] [Accepted: 05/14/2024] [Indexed: 06/07/2024]
Abstract
The study of the kinetic bottlenecks that hinder the rare transitions between long-lived metastable states is a major challenge in atomistic simulations. Here we propose a method to explore the transition state ensemble, which is the distribution of configurations that the system passes through as it translocates from one metastable basin to another. We base our method on the committor function and the variational principle that it obeys. We find its minimum through a self-consistent procedure that starts from information limited to the initial and final states. Right from the start, our procedure allows the sampling of very many transition state configurations. With the help of the variational principle, we perform a detailed analysis of the transition state ensemble, ranking quantitatively the degrees of freedom mostly involved in the transition and enabling a systematic approach for the interpretation of simulation results and the construction of efficient physics-informed collective variables.
Collapse
Affiliation(s)
- Peilin Kang
- Atomistic Simulations, Italian Institute of Technology, Genova, Italy
| | - Enrico Trizio
- Atomistic Simulations, Italian Institute of Technology, Genova, Italy
- Department of Materials Science, Università di Milano-Bicocca, Milano, Italy
| | | |
Collapse
|
3
|
Harris J, Chipot C, Roux B. How is Membrane Permeation of Small Ionizable Molecules Affected by Protonation Kinetics? J Phys Chem B 2024; 128:795-811. [PMID: 38227958 DOI: 10.1021/acs.jpcb.3c06765] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2024]
Abstract
According to the pH-partition hypothesis, the aqueous solution adjacent to a membrane is a mixture of the ionization states of the permeating molecule at fixed Henderson-Hasselbalch concentrations, such that each state passes through the membrane in parallel with its own specific permeability. An alternative view, based on the assumption that the rate of switching ionization states is instantaneous, represents the permeation of ionizable molecules via an effective Boltzmann-weighted average potential (BWAP). Such an assumption is used in constant-pH molecular dynamics simulations. The inhomogeneous solubility-diffusion framework can be used to compute the pH-dependent membrane permeability for each of these two limiting treatments. With biased WTM-eABF molecular dynamics simulations, we computed the potential of mean force and diffusivity of each ionization state of two weakly basic small molecules: nicotine, an addictive drug, and varenicline, a therapeutic for treating nicotine addiction. At pH = 7, the BWAP effective permeability is greater than that determined by pH-partitioning by a factor of 2.5 for nicotine and 5 for varenicline. To assess the importance of ionization kinetics, we present a Smoluchowski master equation that includes explicitly the protonation and deprotonation processes coupled with the diffusive motion across the membrane. At pH = 7, the increase in permeability due to the explicit ionization kinetics is negligible for both nicotine and varenicline. This finding is reaffirmed by combined Brownian dynamics and Markov state model simulations for estimating the permeability of nicotine while allowing changes in its ionization state. We conclude that for these molecules the pH-partition hypothesis correctly captures the physics of the permeation process. The small free energy barriers for the permeation of nicotine and varenicline in their deprotonated neutral forms play a crucial role in establishing the validity of the pH-partitioning mechanism. Essentially, BWAP fails because ionization kinetics are too slow on the time scale of membrane crossing to affect the permeation of small ionizable molecules such as nicotine and varenicline. For the singly protonated state of nicotine, the computational results agree well with experimental measurements (P1 = 1.29 × 10-7 cm/s), but the agreement for neutral (P0 = 6.12 cm/s) and doubly protonated nicotine (P2 = 3.70 × 10-13 cm/s) is slightly worse, likely due to factors associated with the aqueous boundary layer (neutral form) or leaks through paracellular pathways (doubly protonated form).
Collapse
Affiliation(s)
- Jonathan Harris
- Department of Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n◦7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
- Theoretical and Computational Biophysics Group, Beckman Institute, and Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Biochemistry and Molecular Biology, Department of Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Department of Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
4
|
Ji X, Wang R, Wang H, Liu W. On committor functions in milestoning. J Chem Phys 2023; 159:244115. [PMID: 38153148 DOI: 10.1063/5.0180513] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Accepted: 12/07/2023] [Indexed: 12/29/2023] Open
Abstract
As an optimal one-dimensional reaction coordinate, the committor function not only describes the probability of a trajectory initiated at a phase space point first reaching the product state before reaching the reactant state but also preserves the kinetics when utilized to run a reduced dynamics model. However, calculating the committor function in high-dimensional systems poses significant challenges. In this paper, within the framework of milestoning, exact expressions for committor functions at two levels of coarse graining are given, including committor functions of phase space point to point (CFPP) and milestone to milestone (CFMM). When combined with transition kernels obtained from trajectory analysis, these expressions can be utilized to accurately and efficiently compute the committor functions. Furthermore, based on the calculated committor functions, an adaptive algorithm is developed to gradually refine the transition state region. Finally, two model examples are employed to assess the accuracy of these different formulations of committor functions.
Collapse
Affiliation(s)
- Xiaojun Ji
- Research Center for Mathematics and Interdisciplinary Sciences, Shandong University, Qingdao, Shandong 266237, People's Republic of China
- Frontiers Science Center for Nonlinear Expectations (Ministry of Education), Shandong University, Qingdao, Shandong 266237, People's Republic of China
| | - Ru Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| | - Hao Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| |
Collapse
|
5
|
Chen H, Roux B, Chipot C. Discovering Reaction Pathways, Slow Variables, and Committor Probabilities with Machine Learning. J Chem Theory Comput 2023; 19:4414-4426. [PMID: 37224455 PMCID: PMC11372462 DOI: 10.1021/acs.jctc.3c00028] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
A significant challenge faced by atomistic simulations is the difficulty, and often impossibility, to sample the transitions between metastable states of the free-energy landscape associated with slow molecular processes. Importance-sampling schemes represent an appealing option to accelerate the underlying dynamics by smoothing out the relevant free-energy barriers, but require the definition of suitable reaction-coordinate (RC) models expressed in terms of compact low-dimensional sets of collective variables (CVs). While most computational studies of slow molecular processes have traditionally relied on educated guesses based on human intuition to reduce the dimensionality of the problem at hand, a variety of machine-learning (ML) algorithms have recently emerged as powerful alternatives to discover meaningful CVs capable of capturing the dynamics of the slowest degrees of freedom. Considering a simple paradigmatic situation in which the long-time dynamics is dominated by the transition between two known metastable states, we compare two variational data-driven ML methods based on Siamese neural networks aimed at discovering a meaningful RC model─the slowest decorrelating CV of the molecular process, and the committor probability to first reach one of the two metastable states. One method is the state-free reversible variational approach for Markov processes networks (VAMPnets), or SRVs─the other, inspired by the transition path theory framework, is the variational committor-based neural networks, or VCNs. The relationship and the ability of these methodologies to discover the relevant descriptors of the slow molecular process of interest are illustrated with a series of simple model systems. We also show that both strategies are amenable to importance-sampling schemes through an appropriate reweighting algorithm that approximates the kinetic properties of the transition.
Collapse
Affiliation(s)
- Haochuan Chen
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy cedex, France
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, 60637, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy cedex, France
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, 60637, United States
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, and Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
6
|
Jung H, Covino R, Arjun A, Leitold C, Dellago C, Bolhuis PG, Hummer G. Machine-guided path sampling to discover mechanisms of molecular self-organization. NATURE COMPUTATIONAL SCIENCE 2023; 3:334-345. [PMID: 38177937 PMCID: PMC10766509 DOI: 10.1038/s43588-023-00428-z] [Citation(s) in RCA: 33] [Impact Index Per Article: 33.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 03/10/2023] [Indexed: 01/06/2024]
Abstract
Molecular self-organization driven by concerted many-body interactions produces the ordered structures that define both inanimate and living matter. Here we present an autonomous path sampling algorithm that integrates deep learning and transition path theory to discover the mechanism of molecular self-organization phenomena. The algorithm uses the outcome of newly initiated trajectories to construct, validate and-if needed-update quantitative mechanistic models. Closing the learning cycle, the models guide the sampling to enhance the sampling of rare assembly events. Symbolic regression condenses the learned mechanism into a human-interpretable form in terms of relevant physical observables. Applied to ion association in solution, gas-hydrate crystal formation, polymer folding and membrane-protein assembly, we capture the many-body solvent motions governing the assembly process, identify the variables of classical nucleation theory, uncover the folding mechanism at different levels of resolution and reveal competing assembly pathways. The mechanistic descriptions are transferable across thermodynamic states and chemical space.
Collapse
Affiliation(s)
- Hendrik Jung
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Frankfurt am Main, Germany
| | - Roberto Covino
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
| | - A Arjun
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | | | | | - Peter G Bolhuis
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - Gerhard Hummer
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Frankfurt am Main, Germany.
- Institute of Biophysics, Goethe University Frankfurt, Frankfurt am Main, Germany.
| |
Collapse
|
7
|
Song K, Makarov DE, Vouga E. The effect of time resolution on the observed first passage times in diffusive dynamics. J Chem Phys 2023; 158:111101. [PMID: 36948823 DOI: 10.1063/5.0142166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023] Open
Abstract
Single-molecule and single-particle tracking experiments are typically unable to resolve fine details of thermal motion at short timescales where trajectories are continuous. We show that, when a diffusive trajectory xt is sampled at finite time intervals δt, the resulting error in measuring the first passage time to a given domain can exceed the time resolution of the measurement by more than an order of magnitude. Such surprisingly large errors originate from the fact that the trajectory may enter and exit the domain while being unobserved, thereby lengthening the apparent first passage time by an amount that is larger than δt. Such systematic errors are particularly important in single-molecule studies of barrier crossing dynamics. We show that the correct first passage times, as well as other properties of the trajectories such as splitting probabilities, can be recovered via a stochastic algorithm that reintroduces unobserved first passage events probabilistically.
Collapse
Affiliation(s)
- Kevin Song
- Department of Computer Science, University of Texas at Austin, Austin, Texas 78712, USA
| | - Dmitrii E Makarov
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Etienne Vouga
- Department of Computer Science, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
8
|
Hasyim MR, Batton CH, Mandadapu KK. Supervised learning and the finite-temperature string method for computing committor functions and reaction rates. J Chem Phys 2022; 157:184111. [DOI: 10.1063/5.0102423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [ Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. In this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.
Collapse
Affiliation(s)
- Muhammad R. Hasyim
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, USA
| | - Clay H. Batton
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, USA
| | - Kranthi K. Mandadapu
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
9
|
Abstract
The treatment of slow and rare transitions in the simulation of complex systems poses a great computational challenge. A powerful approach to tackle this challenge is the string method, which represents the transition path as a one-dimensional curve in a multidimensional space of collective variables. Commonly used strategies for pathway optimization include aligning the tangent of the string to the local mean force or to the mean drift determined from swarms of short trajectories. Here, a novel strategy is proposed, allowing the string to be optimized based on a variational principle involving the unidirectional reactive flux expressed in terms of the time-correlation function of the committor. The method is illustrated with model systems and then probed with the alanine dipeptide and a coarse-grained model of the barstar-barnase protein complex. Successive iterations variationally refine the string toward an optimal transition pathway following the gradient of the committor between two metastable states.
Collapse
Affiliation(s)
- Ziwei He
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago60637, Illinois, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche No. 7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex54506, France
| | - Benoît Roux
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago60637, Illinois, United States
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago60637, IllinoisUnited States
| |
Collapse
|
10
|
Berezhkovskii AM, Szabo A. Relations among Unidirectional Fluxes at Equilibrium, Committors, and First Passage and Transition Path Times. J Phys Chem B 2022; 126:6624-6628. [PMID: 36037104 DOI: 10.1021/acs.jpcb.2c03757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
For multidimensional diffusive dynamics, we algebraically derive remarkable analytical expressions that relate the mean first passage and transition path times between two dividing surfaces with the number of unidirectional transitions per unit time (fluxes) at equilibrium between the two surfaces and the committor (the probability of reaching one surface before the other). In one dimension, such relationships can be easily derived because analytical expressions for all the above-mentioned quantities can be found. This is not possible in higher dimensions, and at first sight, the problem seems much harder. We circumvent the difficulty that the equations determining the mean first passage and transition path times cannot be solved analytically by multiplying these equations by the committor, integrating both sides and finally using the divergence theorem. A byproduct of our derivation is an analytical expression for the starting point distribution over which individual first passage and transition path times must be averaged. It turns out that this distribution is not the Boltzmann one, but it has a simple physical interpretation.
Collapse
Affiliation(s)
- Alexander M Berezhkovskii
- Section of Molecular Transport, Eunice Kennedy Shriver National Institute of Child Health and Human development, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Attila Szabo
- Laboratory of Chemical Physics, National institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| |
Collapse
|
11
|
Krivov SV. Additive eigenvectors as optimal reaction coordinates, conditioned trajectories, and time-reversible description of stochastic processes. J Chem Phys 2022; 157:014108. [DOI: 10.1063/5.0088061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A fundamental way to analyze complex multidimensional stochastic dynamics is to describe it as diffusion on a free energy landscape—free energy as a function of reaction coordinates (RCs). For such a description to be quantitatively accurate, the RC should be chosen in an optimal way. The committor function is a primary example of an optimal RC for the description of equilibrium reaction dynamics between two states. Here, additive eigenvectors (addevs) are considered as optimal RCs to address the limitations of the committor. An addev master equation for a Markov chain is derived. A stationary solution of the equation describes a sub-ensemble of trajectories conditioned on having the same optimal RC for the forward and time-reversed dynamics in the sub-ensemble. A collection of such sub-ensembles of trajectories, called stochastic eigenmodes, can be used to describe/approximate the stochastic dynamics. A non-stationary solution describes the evolution of the probability distribution. However, in contrast to the standard master equation, it provides a time-reversible description of stochastic dynamics. It can be integrated forward and backward in time. The developed framework is illustrated on two model systems—unidirectional random walk and diffusion.
Collapse
Affiliation(s)
- Sergei V. Krivov
- University of Leeds, Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
12
|
Wu S, Li H, Ma A. A Rigorous Method for Identifying a One-Dimensional Reaction Coordinate in Complex Molecules. J Chem Theory Comput 2022; 18:2836-2844. [PMID: 35427129 DOI: 10.1021/acs.jctc.2c00132] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Understanding the mechanism of functional protein dynamics is critical to understanding protein functions. Reaction coordinate (RC) is a central topic in protein dynamics, and the grail is to find the one-dimensional RC (1D-RC) that can fully determine the value of a committor (i.e., the reaction probability in configuration space) for any protein configuration. We present a new method that, for the first time, uses a fundamental mechanical operator, the generalized work functional, to identify the rigorous 1D-RC in complex molecules. For a prototypical biomolecular isomerization reaction, the 1D-RC identified by the current method can determine the committor with an accuracy far exceeding what was achieved by previous methods. This method only requires modest computational cost and can be readily applied to large molecules. Most importantly, the generalized work functional is the physical determinant of the collectivity in functional protein dynamics and provides a tentative roadmap that connects the structure of a protein to its function.
Collapse
Affiliation(s)
- Shanshan Wu
- Richard Loan and Hill Department of Biomedical Engineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, United States
| | - Huiyu Li
- Richard Loan and Hill Department of Biomedical Engineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, United States
| | - Ao Ma
- Richard Loan and Hill Department of Biomedical Engineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, United States
| |
Collapse
|
13
|
Louwerse MD, Sivak DA. Information Thermodynamics of the Transition-Path Ensemble. PHYSICAL REVIEW LETTERS 2022; 128:170602. [PMID: 35570424 DOI: 10.1103/physrevlett.128.170602] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Accepted: 03/10/2022] [Indexed: 06/15/2023]
Abstract
The reaction coordinate describing a transition between reactant and product is a fundamental concept in the theory of chemical reactions. Within transition-path theory, a quantitative definition of the reaction coordinate is found in the committor, which is the probability that a trajectory initiated from a given microstate first reaches the product before the reactant. Here we develop an information-theoretic origin for the committor and show how selecting transition paths from a long ergodic equilibrium trajectory induces entropy production which exactly equals the information that system dynamics provide about the reactivity of trajectories. This equality of entropy production and dynamical information generation also holds at the level of arbitrary individual coordinates, providing parallel measures of the coordinate's relevance to the reaction, each of which is maximized by the committor.
Collapse
Affiliation(s)
- Miranda D Louwerse
- Department of Chemistry, Simon Fraser University, Burnaby, British Columbia V5A1S6, Canada
| | - David A Sivak
- Department of Physics, Simon Fraser University, Burnaby, British Columbia V5A1S6, Canada
| |
Collapse
|
14
|
Abstract
The kinetics of a dynamical system dominated by two metastable states is examined from the perspective of the activated-dynamics reactive flux formalism, Markov state eigenvalue spectral decomposition, and committor-based transition path theory. Analysis shows that the different theoretical formulations are consistent, clarifying the significance of the inherent microscopic lag-times that are implicated, and that the most meaningful one-dimensional reaction coordinate in the region of the transition state is along the gradient of the committor in the multidimensional subspace of collective variables. It is shown that the familiar reactive flux activated dynamics formalism provides an effective route to calculate the transition rate in the case of a narrow sharp barrier but much less so in the case of a broad flat barrier. In this case, the standard reactive flux correlation function decays very slowly to the plateau value that corresponds to the transmission coefficient. Treating the committor function as a reaction coordinate does not alleviate all issues caused by the slow relaxation of the reactive flux correlation function. A more efficient activated dynamics simulation algorithm may be achieved from a modified reactive flux weighted by the committor. Simulation results on simple systems are used to illustrate the various conceptual points.
Collapse
Affiliation(s)
- Benoît Roux
- Department of Biochemistry and Molecular Biology, Department of Chemistry, The University of Chicago, 5735 S Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
15
|
Wu S, Ma A. Mechanism for the rare fluctuation that powers protein conformational change. J Chem Phys 2022; 156:054119. [PMID: 35135246 PMCID: PMC8824576 DOI: 10.1063/5.0077444] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 01/17/2022] [Indexed: 11/14/2022] Open
Abstract
Most functional processes of biomolecules are rare events. Key to a rare event is the rare fluctuation that enables the energy activation process that precedes and powers crossing of the activation barrier. However, the physical nature of this rare fluctuation and how it enables energy activation and subsequently barrier crossing are unknown. We developed a novel metric, the reaction capacity pC, that rigorously defines the beginning and parameterizes the progress of energy activation. This enabled us to identify the rare fluctuation as a special phase-space condition that is necessary and sufficient for initiating systematic energy flow from the non-reaction coordinates into the reaction coordinates. The energy activation of a prototype biomolecular isomerization reaction is dominated by kinetic energy transferring into and accumulating in the reaction coordinates, administered by inertial forces alone. This mechanism for energy activation is fundamentally different from the mechanism suggested by Kramers theory.
Collapse
Affiliation(s)
- Shanshan Wu
- Richard Loan and Hill Department of Biomedical Engineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, USA
| | - Ao Ma
- Richard Loan and Hill Department of Biomedical Engineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, USA
| |
Collapse
|
16
|
Li W. Optimizing reaction coordinate by flux maximization in the transition path ensemble. J Chem Phys 2022; 156:054117. [DOI: 10.1063/5.0079390] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wenjin Li
- Institute for Advanced Study, Shenzhen University, Shenzhen, China
| |
Collapse
|
17
|
Berezhkovskii AM, Makarov DE. On distributions of barrier crossing times as observed in single-molecule studies of biomolecules. BIOPHYSICAL REPORTS 2021; 1:100029. [PMID: 36425456 PMCID: PMC9680812 DOI: 10.1016/j.bpr.2021.100029] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/15/2021] [Accepted: 10/19/2021] [Indexed: 06/16/2023]
Abstract
Single-molecule experiments that monitor time evolution of molecular observables in real time have expanded beyond measuring transition rates toward measuring distributions of times of various molecular events. Of particular interest is the first-passage time for making a transition from one molecular configuration ( a ) to another ( b ) and conditional first-passage times such as the transition path time, which is the first-passage time from a to b conditional upon not leaving the transition region intervening between a and b . Another experimentally accessible (but not yet studied experimentally) observable is the conditional exit time, i.e., the time to leave the transition region through a specified boundary. The distributions of such times contain a wealth of mechanistic information about the transitions in question. Here, we use the first and the second (and, if desired, higher) moments of these distributions to characterize their relative width for the model in which the experimental observable undergoes Brownian motion in a potential of mean force. We show that although the distributions of transition path times are always narrower than exponential (in that the ratio of the standard deviation to the distribution's mean is always less than 1), distributions of first-passage times and of conditional exit times can be either narrow or broad, in some cases displaying long power-law tails. The conditional exit time studied here provides a generalization of the transition path time that also allows one to characterize the temporal scales of failed barrier crossing attempts.
Collapse
Affiliation(s)
- Alexander M. Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Office of Intramural Research, Center for Information Technology, National Institutes of Health, Bethesda, Maryland
| | - Dmitrii E. Makarov
- Department of Chemistry and Biochemistry and Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas
| |
Collapse
|
18
|
Roux B. String Method with Swarms-of-Trajectories, Mean Drifts, Lag Time, and Committor. J Phys Chem A 2021; 125:7558-7571. [PMID: 34406010 PMCID: PMC8419867 DOI: 10.1021/acs.jpca.1c04110] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2021] [Revised: 07/26/2021] [Indexed: 11/29/2022]
Abstract
The kinetics of a dynamical system comprising two metastable states is formulated in terms of a finite-time propagator in phase space (position and velocity) adapted to the underdamped Langevin equation. Dimensionality reduction to a subspace of collective variables yields familiar expressions for the propagator, committor, and steady-state flux. A quadratic expression for the steady-state flux between the two metastable states can serve as a robust variational principle to determine an optimal approximate committor expressed in terms of a set of collective variables. The theoretical formulation is exploited to clarify the foundation of the string method with swarms-of-trajectories, which relies on the mean drift of short trajectories to determine the optimal transition pathway. It is argued that the conditions for Markovity within a subspace of collective variables may not be satisfied with an arbitrary short time-step and that proper kinetic behaviors appear only when considering the effective propagator for longer lag times. The effective propagator with finite lag time is amenable to an eigenvalue-eigenvector spectral analysis, as elaborated previously in the context of position-based Markov models. The time-correlation functions calculated by swarms-of-trajectories along the string pathway constitutes a natural extension of these developments. The present formulation provides a powerful theoretical framework to characterize the optimal pathway between two metastable states of a system.
Collapse
Affiliation(s)
- Benoît Roux
- Department
of Biochemistry and Molecular Biology, The
University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
19
|
Blow KE, Quigley D, Sosso GC. The seven deadly sins: When computing crystal nucleation rates, the devil is in the details. J Chem Phys 2021; 155:040901. [PMID: 34340373 DOI: 10.1063/5.0055248] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The formation of crystals has proven to be one of the most challenging phase transformations to quantitatively model-let alone to actually understand-be it by means of the latest experimental technique or the full arsenal of enhanced sampling approaches at our disposal. One of the most crucial quantities involved with the crystallization process is the nucleation rate, a single elusive number that is supposed to quantify the average probability for a nucleus of critical size to occur within a certain volume and time span. A substantial amount of effort has been devoted to attempt a connection between the crystal nucleation rates computed by means of atomistic simulations and their experimentally measured counterparts. Sadly, this endeavor almost invariably fails to some extent, with the venerable classical nucleation theory typically blamed as the main culprit. Here, we review some of the recent advances in the field, focusing on a number of perhaps more subtle details that are sometimes overlooked when computing nucleation rates. We believe it is important for the community to be aware of the full impact of aspects, such as finite size effects and slow dynamics, that often introduce inconspicuous and yet non-negligible sources of uncertainty into our simulations. In fact, it is key to obtain robust and reproducible trends to be leveraged so as to shed new light on the kinetics of a process, that of crystal nucleation, which is involved into countless practical applications, from the formulation of pharmaceutical drugs to the manufacturing of nano-electronic devices.
Collapse
Affiliation(s)
- Katarina E Blow
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David Quigley
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
20
|
Berezhkovskii AM, Gopich IV, Szabo A. Diffusive barrier crossing rates from variationally determined eigenvalues. J Chem Phys 2021; 155:034104. [PMID: 34293906 PMCID: PMC8411888 DOI: 10.1063/5.0058066] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 06/30/2021] [Indexed: 11/14/2022] Open
Abstract
Kramers' procedure for calculating the rate of activated processes involves partitioning space into reactant, barrier, and product regions by introducing two dividing surfaces. Then, a nonequilibrium steady state is established by injecting particles on one surface and removing them when they reach the other. The rate is obtained as the ratio of the steady-state flux between the surfaces and the population of the initial well. An alternative procedure that seems less artificial is to estimate the first non-zero eigenvalue of the operator that describes the dynamics and then equate its magnitude to the sum of the forward and backward rate constants. Here, we establish the relationship between these approaches for diffusive dynamics, starting with the variational principle for the eigenvalue of interest and then using a trial function involving two adjustable surfaces. We show how Kramers' flux-over-population expression for the rate constant can be obtained from our variationally determined eigenvalue in the special case where the reactant and product regions are separated by a high barrier. This work exploits the modern theory of activated rate processes where the committor (the probability of reaching one dividing surface before the other) plays a central role. Surprisingly, our upper bound for the eigenvalue can be expressed solely in terms of mean first-passage times and the mean transition-path time between the two dividing surfaces.
Collapse
Affiliation(s)
- Alexander M. Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Office of Intramural Research, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Irina V. Gopich
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 208192, USA
| | - Attila Szabo
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 208192, USA
| |
Collapse
|
21
|
Gray TH, Yong EH. An effective one-dimensional approach to calculating mean first passage time in multi-dimensional potentials. J Chem Phys 2021; 154:084103. [PMID: 33639738 DOI: 10.1063/5.0040071] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Thermally activated escape processes in multi-dimensional potentials are of interest to a variety of fields, so being able to calculate the rate of escape-or the mean first-passage time (MFPT)-is important. Unlike in one dimension, there is no general, exact formula for the MFPT. However, Langer's formula, a multi-dimensional generalization of Kramers's one-dimensional formula, provides an approximate result when the barrier to escape is large. Kramers's and Langer's formulas are related to one another by the potential of mean force (PMF): when calculated along a particular direction (the unstable mode at the saddle point) and substituted into Kramers's formula, the result is Langer's formula. We build on this result by using the PMF in the exact, one-dimensional expression for the MFPT. Our model offers better agreement with Brownian dynamics simulations than Langer's formula, although discrepancies arise when the potential becomes less confining along the direction of escape. When the energy barrier is small our model offers significant improvements upon Langer's theory. Finally, the optimal direction along which to evaluate the PMF no longer corresponds to the unstable mode at the saddle point.
Collapse
Affiliation(s)
- Thomas H Gray
- Department of Chemical Engineering and Biotechnology, West Cambridge Site, University of Cambridge, Philippa Fawcett Drive, CB3 0AS Cambridge, United Kingdom
| | - Ee Hou Yong
- Division of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| |
Collapse
|
22
|
Torrillo PA, Bogetti AT, Chong LT. A Minimal, Adaptive Binning Scheme for Weighted Ensemble Simulations. J Phys Chem A 2021; 125:1642-1649. [PMID: 33577732 PMCID: PMC8091492 DOI: 10.1021/acs.jpca.0c10724] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A promising approach for simulating rare events with rigorous kinetics is the weighted ensemble path sampling strategy. One challenge of this strategy is the division of configurational space into bins for sampling. Here we present a minimal adaptive binning (MAB) scheme for the automated, adaptive placement of bins along a progress coordinate within the framework of the weighted ensemble strategy. Results reveal that the MAB binning scheme, despite its simplicity, is more efficient than a manual, fixed binning scheme in generating transitions over large free energy barriers, generating a diversity of pathways, estimating rate constants, and sampling conformations. The scheme is general and extensible to any rare-events sampling strategy that employs progress coordinates.
Collapse
Affiliation(s)
- Paul A Torrillo
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Anthony T Bogetti
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Lillian T Chong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
23
|
Sigg D, Voelz VA, Carnevale V. Microcanonical coarse-graining of the kinetic Ising model. J Chem Phys 2020; 152:084104. [PMID: 32113343 PMCID: PMC7042020 DOI: 10.1063/1.5139228] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Accepted: 02/06/2020] [Indexed: 11/15/2022] Open
Abstract
We propose a scheme for coarse-graining the dynamics of the 2-D kinetic Ising model onto the microcanonical ensemble. At subcritical temperatures, 2-D and higher-dimensional Ising lattices possess two basins of attraction separated by a free energy barrier. Projecting onto the microcanonical ensemble has the advantage that the dependence of the crossing rate constant on environmental conditions can be obtained from a single Monte Carlo trajectory. Using various numerical methods, we computed the forward rate constants of coarse-grained representations of the Ising model and compared them with the true value obtained from brute force simulation. While coarse-graining preserves detailed balance, the computed rate constants for barrier heights between 5 kT and 9 kT were consistently 50% larger than the true value. Markovianity testing revealed loss of dynamical memory, which we propose accounts for coarse-graining error. Committor analysis did not support the alternative hypothesis that microcanonical projection is incompatible with an optimal reaction coordinate. The correct crossing rate constant was obtained by spectrally decomposing the diffusion coefficient near the free energy barrier and selecting the slowest (reactive) component. The spectral method also yielded the correct rate constant in the 3-D Ising lattice, where coarse-graining error was 6% and memory effects were diminished. We conclude that microcanonical coarse-graining supplemented by spectral analysis of short-term barrier fluctuations provides a comprehensive kinetic description of barrier crossing in a non-inertial continuous-time jump process.
Collapse
Affiliation(s)
| | - Vincent A. Voelz
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Vincenzo Carnevale
- Institute for Computational Molecular Science, College of Science and Technology, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
24
|
Chiuchiù D, Ferrare J, Pigolotti S. Assembly of heteropolymers via a network of reaction coordinates. Phys Rev E 2019; 100:062502. [PMID: 31962425 DOI: 10.1103/physreve.100.062502] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Indexed: 06/10/2023]
Abstract
In biochemistry, heteropolymers encoding biological information are assembled out of equilibrium by sequentially incorporating available monomers found in the environment. Current models of polymerization treat monomer incorporation as a sequence of discrete chemical reactions between intermediate metastable states. In this paper, we use ideas from reaction rate theory and describe nonequilibrium assembly of a heteropolymer via a continuous reaction coordinate. Our approach allows for estimating the copy error and incorporation speed from the Gibbs free energy landscape of the process. We apply our theory to several examples from a simple reaction characterized by a free energy barrier to more complex cases incorporating error correction mechanisms, such as kinetic proofreading.
Collapse
Affiliation(s)
- Davide Chiuchiù
- Biological Complexity Unit, Okinawa Institute for Science and Technology, 1919-1 Tancha, Onna, Kunigami-gun, Okinawa 904-0412, Japan
| | - James Ferrare
- Biological Complexity Unit, Okinawa Institute for Science and Technology, 1919-1 Tancha, Onna, Kunigami-gun, Okinawa 904-0412, Japan
- Tulane University, 6823 St. Charles Avenue, New Orleans, Lousiana 70118, USA
| | - Simone Pigolotti
- Biological Complexity Unit, Okinawa Institute for Science and Technology, 1919-1 Tancha, Onna, Kunigami-gun, Okinawa 904-0412, Japan
| |
Collapse
|
25
|
Yappert R, Kamat K, Peters B. The overdamped transmission coefficient: Recovering the true mean first passage time from an inaccurate reaction coordinate. J Chem Phys 2019; 151:184108. [DOI: 10.1063/1.5117237] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Ryan Yappert
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Kartik Kamat
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
26
|
Copperman J, Aristoff D, Makarov DE, Simpson G, Zuckerman DM. Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture. J Chem Phys 2019; 151:174108. [PMID: 31703496 PMCID: PMC7043855 DOI: 10.1063/1.5120511] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Accepted: 10/14/2019] [Indexed: 01/04/2023] Open
Abstract
Probability currents are fundamental in characterizing the kinetics of nonequilibrium processes. Notably, the steady-state current Jss for a source-sink system can provide the exact mean-first-passage time (MFPT) for the transition from the source to sink. Because transient nonequilibrium behavior is quantified in some modern path sampling approaches, such as the "weighted ensemble" strategy, there is strong motivation to determine bounds on Jss-and hence on the MFPT-as the system evolves in time. Here, we show that Jss is bounded from above and below by the maximum and minimum, respectively, of the current as a function of the spatial coordinate at any time t for one-dimensional systems undergoing overdamped Langevin (i.e., Smoluchowski) dynamics and for higher-dimensional Smoluchowski systems satisfying certain assumptions when projected onto a single dimension. These bounds become tighter with time, making them of potential practical utility in a scheme for estimating Jss and the long time scale kinetics of complex systems. Conceptually, the bounds result from the fact that extrema of the transient currents relax toward the steady-state current.
Collapse
Affiliation(s)
- Jeremy Copperman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239, USA
| | - David Aristoff
- Department of Mathematics, Colorado State University, Fort Collins, Colorado 80523, USA
| | - Dmitrii E Makarov
- Department of Chemistry and Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas 78712, USA
| | - Gideon Simpson
- Department of Mathematics, Drexel University, Philadelphia, Pennsylvania 19104, USA
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239, USA
| |
Collapse
|
27
|
Berezhkovskii AM, Szabo A. Committors, first-passage times, fluxes, Markov states, milestones, and all that. J Chem Phys 2019; 150:054106. [PMID: 30736684 PMCID: PMC6910584 DOI: 10.1063/1.5079742] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Accepted: 01/05/2019] [Indexed: 11/14/2022] Open
Abstract
Milestoning on a one-dimensional potential starts by choosing a set of points, called milestones, and initiating short trajectories from each milestone, which are terminated when they reach an adjacent milestone for the first time. From the average duration of these trajectories and the probabilities of where they terminate, a rate matrix can be constructed and then used to calculate the mean first-passage time (MFPT) between any two milestones. All these MFPT's turn out to be exact. Here we adopt a point of view from which this remarkable result is not unexpected. In addition, we clarify the nature of the "states" whose interconversion is described by the rate matrix constructed using information obtained from short trajectories and provide a microscopic expression for the "equilibrium population" of these states in terms of equilibrium averages of the committors.
Collapse
Affiliation(s)
- Alexander M. Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Office of Intramural Research, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Attila Szabo
- Laboratory of Chemical Physics, National institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 208192, USA
| |
Collapse
|
28
|
Abstract
A formalism is developed to describe how diffusion alters the kinetics of coupled reversible association-dissociation reactions in the presence of conformational changes that can modify the reactivity. The major difficulty in constructing a general theory is that, even to the lowest order, diffusion can change the structure of the rate equations of chemical kinetics by introducing new reaction channels (i.e., modifies the kinetic scheme). Therefore, the right formalism must be found that allows the influence of diffusion to be described in a concise and elegant way for networks of arbitrary complexity. Our key result is a set of non-Markovian rate equations involving stoichiometric matrices and net reaction rates (fluxes), in which these rates are coupled by a time-dependent pair association flux matrix, whose elements have a simple physical interpretation. Specifically, each element is the probability density that an isolated pair of reactants irreversibly associates at time t via one reaction channel on the condition that it started out with the dissociation products of another (or the same) channel. In the Markovian limit, the coupling of the chemical rates is described by committors (or splitting/capture probabilities). The committor is the probability that an isolated pair of reactants formed by dissociation at one site will irreversibly associate at another site rather than diffuse apart. We illustrate the use of our formalism by considering three reversible reaction schemes: (1) binding to a single site, (2) binding to two inequivalent sites, and (3) binding to a site whose reactivity fluctuates. In the first example, we recover the results published earlier, while in the second one we show that a new reaction channel appears, which directly connects the two bound states. The third example is particularly interesting because all species become coupled and an exchange-type bimolecular reaction appears. In the Markovian limit, some of the diffusion-modified rate constants that describe new transitions become negative, indicating that memory effects cannot be ignored.
Collapse
Affiliation(s)
- Irina V. Gopich
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Attila Szabo
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, USA
| |
Collapse
|
29
|
Das A, Makarov DE. Dynamics of Disordered Proteins under Confinement: Memory Effects and Internal Friction. J Phys Chem B 2018; 122:9049-9060. [PMID: 30092636 DOI: 10.1021/acs.jpcb.8b06112] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Many proteins are disordered under physiological conditions. How efficiently they can search for their cellular targets and how fast they can fold upon target binding is determined by their intrinsic dynamics, which have thus attracted much recent attention. Experiments and molecular simulations show that the inherent reconfiguration timescale for unfolded proteins has a solvent friction component and an internal friction component, and the microscopic origin of the latter, along with its proper mathematical description, has been a topic of considerable debate. Internal friction varies across different proteins of comparable length and increases with decreasing denaturant concentration, showing that it depends on how compact the protein is. Here we report on a systematic atomistic simulation study of how confinement, which induces a more compact unfolded state, affects dynamics and friction in disordered peptides. We find that the average reconfiguration timescales increase exponentially as the peptide's spatial dimensions are reduced; at the same time, confinement broadens the spectrum of relaxation timescales exhibited by the peptide. There are two important implications of this broadening: First, it limits applicability of the common Rouse and Zimm models with internal friction, as those models attempt to capture internal friction effects by introducing a single internal friction timescale. Second, the long-tailed distribution of relaxation times leads to anomalous diffusion effects in the dynamics of intramolecular distances. Analysis and interpretation of experimental signals from various measurements that probe intramolecular protein dynamics (such as single-molecule fluorescence correlation spectroscopy and single-molecule force spectroscopy) rely on the assumption of diffusive dynamics along the distances being probed; hence, our results suggest the need for more general models allowing for anomalous diffusion effects.
Collapse
Affiliation(s)
- Atanu Das
- Department of Chemistry , University of Texas at Austin , Austin , Texas 78712 , United States
| | - Dmitrii E Makarov
- Department of Chemistry , University of Texas at Austin , Austin , Texas 78712 , United States.,Institute for Computational Engineering and Sciences , University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
30
|
Satija R, Das A, Makarov DE. Transition path times reveal memory effects and anomalous diffusion in the dynamics of protein folding. J Chem Phys 2018; 147:152707. [PMID: 29055292 DOI: 10.1063/1.4993228] [Citation(s) in RCA: 49] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Recent single-molecule experiments probed transition paths of biomolecular folding and, in particular, measured the time biomolecules spend while crossing their free energy barriers. A surprising finding from these studies is that the transition barriers crossed by transition paths, as inferred from experimentally observed transition path times, are often lower than the independently determined free energy barriers. Here we explore memory effects leading to anomalous diffusion as a possible origin of this discrepancy. Our analysis of several molecular dynamics trajectories shows that the dynamics of common reaction coordinates used to describe protein folding is subdiffusive, at least at sufficiently short times. We capture this effect using a one-dimensional fractional Brownian motion (FBM) model, in which the system undergoes a subdiffusive process in the presence of a potential of mean force, and show that this model yields much broader distributions of transition path times with stretched exponential long-time tails. Without any adjustable parameters, these distributions agree well with the transition path times computed directly from protein trajectories. We further discuss how the FBM model can be tested experimentally.
Collapse
Affiliation(s)
- Rohit Satija
- Department of Chemistry and Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas 78712, USA
| | - Atanu Das
- Department of Chemistry and Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas 78712, USA
| | - Dmitrii E Makarov
- Department of Chemistry and Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
31
|
Berezhkovskii AM, Makarov DE. Communication: Coordinate-dependent diffusivity from single molecule trajectories. J Chem Phys 2018; 147:201102. [PMID: 29195291 DOI: 10.1063/1.5006456] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Single-molecule observations of biomolecular folding are commonly interpreted using the model of one-dimensional diffusion along a reaction coordinate, with a coordinate-independent diffusion coefficient. Recent analysis, however, suggests that more general models are required to account for single-molecule measurements performed with high temporal resolution. Here, we consider one such generalization: a model where the diffusion coefficient can be an arbitrary function of the reaction coordinate. Assuming Brownian dynamics along this coordinate, we derive an exact expression for the coordinate-dependent diffusivity in terms of the splitting probability within an arbitrarily chosen interval and the mean transition path time between the interval boundaries. This formula can be used to estimate the effective diffusion coefficient along a reaction coordinate directly from single-molecule trajectories.
Collapse
Affiliation(s)
- Alexander M Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Office of Intramural Research, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Dmitrii E Makarov
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
32
|
Krivov SV. Protein Folding Free Energy Landscape along the Committor - the Optimal Folding Coordinate. J Chem Theory Comput 2018; 14:3418-3427. [PMID: 29791148 DOI: 10.1021/acs.jctc.8b00101] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Recent advances in simulation and experiment have led to dramatic increases in the quantity and complexity of produced data, which makes the development of automated analysis tools very important. A powerful approach to analyze dynamics contained in such data sets is to describe/approximate it by diffusion on a free energy landscape - free energy as a function of reaction coordinates (RC). For the description to be quantitatively accurate, RCs should be chosen in an optimal way. Recent theoretical results show that such an optimal RC exists; however, determining it for practical systems is a very difficult unsolved problem. Here we describe a solution to this problem. We describe an adaptive nonparametric approach to accurately determine the optimal RC (the committor) for an equilibrium trajectory of a realistic system. In contrast to alternative approaches, which require a functional form with many parameters to approximate an RC and thus extensive expertise with the system, the suggested approach is nonparametric and can approximate any RC with high accuracy without system specific information. To avoid overfitting for a realistically sampled system, the approach performs RC optimization in an adaptive manner by focusing optimization on less optimized spatiotemporal regions of the RC. The power of the approach is illustrated on a long equilibrium atomistic folding simulation of HP35 protein. We have determined the optimal folding RC - the committor, which was confirmed by passing a stringent committor validation test. It allowed us to determine a first quantitatively accurate protein folding free energy landscape. We have confirmed the recent theoretical results that diffusion on such a free energy profile can be used to compute exactly the equilibrium flux, the mean first passage times, and the mean transition path times between any two points on the profile. We have shown that the mean squared displacement along the optimal RC grows linear with time as for simple diffusion. The free energy profile allowed us to obtain a direct rigorous estimate of the pre-exponential factor for the folding dynamics.
Collapse
Affiliation(s)
- Sergei V Krivov
- Astbury Center for Structural Molecular Biology , University of Leeds , Leeds LS2 9JT , United Kingdom
| |
Collapse
|
33
|
Tiwary P, Berne BJ. Predicting reaction coordinates in energy landscapes with diffusion anisotropy. J Chem Phys 2017; 147:152701. [PMID: 29055314 PMCID: PMC5446309 DOI: 10.1063/1.4983727] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2017] [Accepted: 05/03/2017] [Indexed: 01/21/2023] Open
Abstract
We consider a range of model potentials with metastable states undergoing molecular dynamics coupled to a thermal bath in the high friction regime and consider how the optimal reaction coordinate depends on the diffusion anisotropy. For this we use our recently proposed method "spectral gap optimization of order parameters (SGOOP)" [P. Tiwary and B. J. Berne, Proc. Natl. Acad. Sci. U. S. A. 113, 2839 (2016)]. We show how available information about dynamical observables in addition to static information can be incorporated into SGOOP, which can then be used to accurately determine the "best" reaction coordinate for arbitrary anisotropies. We compare our results with transmission coefficient calculations and published benchmarks wherever applicable or available, respectively.
Collapse
Affiliation(s)
- Pratyush Tiwary
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - B J Berne
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
34
|
McGibbon RT, Husic BE, Pande VS. Identification of simple reaction coordinates from complex dynamics. J Chem Phys 2017; 146:044109. [PMID: 28147508 PMCID: PMC5272828 DOI: 10.1063/1.4974306] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2016] [Accepted: 01/05/2017] [Indexed: 11/14/2022] Open
Abstract
Reaction coordinates are widely used throughout chemical physics to model and understand complex chemical transformations. We introduce a definition of the natural reaction coordinate, suitable for condensed phase and biomolecular systems, as a maximally predictive one-dimensional projection. We then show that this criterion is uniquely satisfied by a dominant eigenfunction of an integral operator associated with the ensemble dynamics. We present a new sparse estimator for these eigenfunctions which can search through a large candidate pool of structural order parameters and build simple, interpretable approximations that employ only a small number of these order parameters. Example applications with a small molecule's rotational dynamics and simulations of protein conformational change and folding show that this approach can filter through statistical noise to identify simple reaction coordinates from complex dynamics.
Collapse
Affiliation(s)
- Robert T McGibbon
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Brooke E Husic
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
35
|
Abstract
Despite its relevance in biology and engineering, the molecular mechanism driving cavitation in water remains unknown. Using computer simulations, we investigate the structure and dynamics of vapor bubbles emerging from metastable water at negative pressures. We find that in the early stages of cavitation, bubbles are irregularly shaped and become more spherical as they grow. Nevertheless, the free energy of bubble formation can be perfectly reproduced in the framework of classical nucleation theory (CNT) if the curvature dependence of the surface tension is taken into account. Comparison of the observed bubble dynamics to the predictions of the macroscopic Rayleigh-Plesset (RP) equation, augmented with thermal fluctuations, demonstrates that the growth of nanoscale bubbles is governed by viscous forces. Combining the dynamical prefactor determined from the RP equation with CNT based on the Kramers formalism yields an analytical expression for the cavitation rate that reproduces the simulation results very well over a wide range of pressures. Furthermore, our theoretical predictions are in excellent agreement with cavitation rates obtained from inclusion experiments. This suggests that homogeneous nucleation is observed in inclusions, whereas only heterogeneous nucleation on impurities or defects occurs in other experiments.
Collapse
|
36
|
Banushkina PV, Krivov SV. Optimal reaction coordinates. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1276] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Polina V. Banushkina
- Astbury Center for Structural Molecular Biology; Faculty of Biological Sciences, University of Leeds; Leeds UK
| | - Sergei V. Krivov
- Astbury Center for Structural Molecular Biology; Faculty of Biological Sciences, University of Leeds; Leeds UK
| |
Collapse
|
37
|
Affiliation(s)
- Baron Peters
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106;
| |
Collapse
|
38
|
Menzl G, Singraber A, Dellago C. S-shooting: a Bennett–Chandler-like method for the computation of rate constants from committor trajectories. Faraday Discuss 2016; 195:345-364. [DOI: 10.1039/c6fd00124f] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Mechanisms of rare transitions between long-lived stable states are often analyzed in terms of commitment probabilities, determined from swarms of short molecular dynamics trajectories. Here, we present a computer simulation method to determine rate constants from such short trajectories combined with free energy calculations. The method, akin to the Bennett–Chandler approach for the calculation of reaction rate constants, requires the definition of a valid reaction coordinate and can be applied to both under- and overdamped dynamics. We verify the correctness of the algorithm using a one-dimensional random walker in a double-well potential and demonstrate its applicability to complex transitions in condensed systems by calculating cavitation rates for water at negative pressures.
Collapse
Affiliation(s)
- Georg Menzl
- Faculty of Physics and Center for Computational Materials Science
- University of Vienna
- 1090 Vienna
- Austria
| | - Andreas Singraber
- Faculty of Physics and Center for Computational Materials Science
- University of Vienna
- 1090 Vienna
- Austria
| | - Christoph Dellago
- Faculty of Physics and Center for Computational Materials Science
- University of Vienna
- 1090 Vienna
- Austria
| |
Collapse
|
39
|
Banushkina PV, Krivov SV. Nonparametric variational optimization of reaction coordinates. J Chem Phys 2015; 143:184108. [DOI: 10.1063/1.4935180] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Polina V. Banushkina
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Sergei V. Krivov
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
40
|
Berezhkovskii AM, Szabo A, Greives N, Zhou HX. Multidimensional reaction rate theory with anisotropic diffusion. J Chem Phys 2015; 141:204106. [PMID: 25429932 DOI: 10.1063/1.4902243] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
An analytical expression is derived for the rate constant that describes diffusive transitions between two deep wells of a multidimensional potential. The expression, in contrast to the Kramers-Langer formula for the rate constant, is valid even when the diffusion is highly anisotropic. Our approach is based on a variational principle for the reactive flux and uses a trial function for the splitting probability or commitor. The theoretical result is validated by Brownian dynamics simulations.
Collapse
Affiliation(s)
- Alexander M Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Division of Computational Bioscience, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20819, USA
| | - Attila Szabo
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20819, USA
| | - Nicholas Greives
- Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, USA
| | - Huan-Xiang Zhou
- Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, USA
| |
Collapse
|
41
|
Paul S, Taraphder S. Determination of the Reaction Coordinate for a Key Conformational Fluctuation in Human Carbonic Anhydrase II. J Phys Chem B 2015; 119:11403-15. [DOI: 10.1021/acs.jpcb.5b03655] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, India
| |
Collapse
|
42
|
Gopich IV, Szabo A. Influence of diffusion on the kinetics of multisite phosphorylation. Protein Sci 2015; 25:244-54. [PMID: 26096178 DOI: 10.1002/pro.2722] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2015] [Accepted: 06/05/2015] [Indexed: 11/06/2022]
Abstract
When an enzyme modifies multiple sites on a substrate, the influence of the relative diffusive motion of the reactants cannot be described by simply altering the rate constants in the rate equations of chemical kinetics. We have recently shown that, even as a first approximation, new transitions between the appropriate species must also be introduced. The physical reason for this is that a kinase, after phosphorylating one site, can rebind and modify another site instead of diffusing away. The corresponding new rate constants depend on the capture or rebinding probabilities that an enzyme-substrate pair, which is formed after dissociation from one site, reacts at the other site rather than diffusing apart. Here we generalize our previous work to describe both random and sequential phosphorylation by considering inequivalent modification sites. In addition, anisotropic reactive sites (instead of uniformly reactive spheres) are explicitly treated by using localized sink and source terms in the reaction-diffusion equations for the enzyme-substrate pair distribution function. Finally, we show that our results can be rederived using a phenomenological approach based on introducing transient encounter complexes into the standard kinetic scheme and then eliminating them using the steady-state approximation.
Collapse
Affiliation(s)
- Irina V Gopich
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, 20892
| | - Attila Szabo
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, 20892
| |
Collapse
|
43
|
Lu J, Vanden-Eijnden E. Exact dynamical coarse-graining without time-scale separation. J Chem Phys 2015; 141:044109. [PMID: 25084883 DOI: 10.1063/1.4890367] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
A family of collective variables is proposed to perform exact dynamical coarse-graining even in systems without time scale separation. More precisely, it is shown that these variables are not slow in general, yet satisfy an overdamped Langevin equation that statistically preserves the sequence in which any regions in collective variable space are visited and permits to calculate exactly the mean first passage times from any such region to another. The role of the free energy and diffusion coefficient in this overdamped Langevin equation is discussed, along with the way they transform under any change of variable in collective variable space. These results apply both to systems with and without inertia, and they can be generalized to using several collective variables simultaneously. The view they offer on what makes collective variables and reaction coordinates optimal breaks from the standard notion that good collective variable must be slow variable, and it suggests new ways to interpret data from molecular dynamics simulations and experiments.
Collapse
Affiliation(s)
- Jianfeng Lu
- Departments of Mathematics, Physics, and Chemistry, Duke University, Box 90320, Durham, North Carolina 27708, USA
| | - Eric Vanden-Eijnden
- Courant Institute of Mathematical Sciences, New York University, 251 Mercer St., New York, New York 10012, USA
| |
Collapse
|
44
|
Affiliation(s)
- Baron Peters
- Department
of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Department
of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
45
|
Banushkina PV, Krivov SV. Fep1d: A script for the analysis of reaction coordinates. J Comput Chem 2015; 36:878-82. [DOI: 10.1002/jcc.23868] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2014] [Revised: 01/29/2015] [Accepted: 01/30/2015] [Indexed: 01/23/2023]
Affiliation(s)
- Polina V. Banushkina
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences; University of Leeds; Leeds LS2 9JT United Kingdom
| | - Sergei V. Krivov
- Astbury Center for Structural Molecular Biology, Faculty of Biological Sciences; University of Leeds; Leeds LS2 9JT United Kingdom
| |
Collapse
|
46
|
Noel JK, Whitford PC. How Simulations Reveal Dynamics, Disorder, and the Energy Landscapes of Biomolecular Function. Isr J Chem 2014. [DOI: 10.1002/ijch.201400018] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
47
|
Li W, Ma A. Recent developments in methods for identifying reaction coordinates. MOLECULAR SIMULATION 2014; 40:784-793. [PMID: 25197161 PMCID: PMC4152980 DOI: 10.1080/08927022.2014.907898] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
In the study of rare events in complex systems with many degrees of freedom, a key element is to identify the reaction coordinates of a given process. Over recent years, a number of methods and protocols have been developed to extract the reaction coordinates based on limited information from molecular dynamics simulations. In this review, we provide a brief survey over a number of major methods developed in the past decade, some of which are discussed in greater detail, to provide an overview of the problems that are partially solved and challenges that still remain. A particular emphasis has been placed on methods for identifying reaction coordinates that are related to the committor.
Collapse
Affiliation(s)
- Wenjin Li
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, IL 60607, USA
| | - Ao Ma
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, IL 60607, USA
| |
Collapse
|
48
|
Krivov SV. Method to describe stochastic dynamics using an optimal coordinate. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:062131. [PMID: 24483410 DOI: 10.1103/physreve.88.062131] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2013] [Indexed: 06/03/2023]
Abstract
A general method to describe the stochastic dynamics of Markov processes is suggested. The method aims to solve three related problems: the determination of an optimal coordinate for the description of stochastic dynamics; the reconstruction of time from an ensemble of stochastic trajectories; and the decomposition of stationary stochastic dynamics into eigenmodes which do not decay exponentially with time. The problems are solved by introducing additive eigenvectors which are transformed by a stochastic matrix in a simple way - every component is translated by a constant distance. Such solutions have peculiar properties. For example, an optimal coordinate for stochastic dynamics with detailed balance is a multivalued function. An optimal coordinate for a random walk on a line corresponds to the conventional eigenvector of the one-dimensional Dirac equation. The equation for the optimal coordinate in a slowly varying potential reduces to the Hamilton-Jacobi equation for the action function.
Collapse
Affiliation(s)
- Sergei V Krivov
- Astbury Center for Structural Molecular Biology, University of Leeds, Leeds, LS2 9JT, United Kingdom
| |
Collapse
|