1
|
Koehl P, Navaza R, Tekpinar M, Delarue M. MinActionPath2: path generation between different conformations of large macromolecular assemblies by action minimization. Nucleic Acids Res 2024; 52:W256-W263. [PMID: 38783081 PMCID: PMC11223808 DOI: 10.1093/nar/gkae421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Revised: 04/25/2024] [Accepted: 05/07/2024] [Indexed: 05/25/2024] Open
Abstract
Recent progress in solving macromolecular structures and assemblies by cryogenic electron microscopy techniques enables sampling of their conformations in different states that are relevant to their biological function. Knowing the transition path between these conformations would provide new avenues for drug discovery. While the experimental study of transition paths is intrinsically difficult, in-silico methods can be used to generate an initial guess for those paths. The Elastic Network Model (ENM), along with a coarse-grained representation (CG) of the structures are among the most popular models to explore such possible paths. Here we propose an update to our software platform MinActionPath that generates non-linear transition paths based on ENM and CG models, using action minimization to solve the equations of motion. The new website enables the study of large structures such as ribosomes or entire virus envelopes. It provides direct visualization of the trajectories along with quantitative analyses of their behaviors at http://dynstr.pasteur.fr/servers/minactionpath/minactionpath2_submission.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Science and Genome Centre, University of California, Davis, CA 95616, USA
| | - Rafael Navaza
- Plateforme de Cristallographie, C2RT, Institut Pasteur, Université Paris Cité, UMR 3528 du CNRS, 75015 Paris, France
| | - Mustafa Tekpinar
- Unité Architecture et Dynamique des Macromolécules Biologiques, Institut Pasteur, Université Paris Cité, UMR 3528 du CNRS, 75015 Paris, France
| | - Marc Delarue
- Unité Architecture et Dynamique des Macromolécules Biologiques, Institut Pasteur, Université Paris Cité, UMR 3528 du CNRS, 75015 Paris, France
| |
Collapse
|
2
|
Koehl P, Orland H. Sampling constrained stochastic trajectories using Brownian bridges. J Chem Phys 2022; 157:054105. [DOI: 10.1063/5.0102295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a new method to sample conditioned trajectories of a system evolving under Langevin dynamics, based on Brownian bridges. <p>The trajectories are conditioned to end at a certain point (or in a certain region) in space.</p> <p>The bridge equations can be recast exactly in the form of a non linear stochastic integro-differential equation.</p> <p>This equation can be very well approximated when the trajectories are closely bundled together in space, i.e. at low temperature, or for transition paths. The approximate equation can be solved iteratively, using a fixed point method.</p> <p>We discuss how to choose the initial trajectories and show some examples of the performance of this method on some simple problems.</p> <p>The method allows to generate conditioned trajectories with a high accuracy.
Collapse
Affiliation(s)
- Patrice Koehl
- Computer Science and Genome Center, University of California Davis, United States of America
| | - Henri Orland
- Institut de Physique Theorique, CEA, Saclay, France
| |
Collapse
|
3
|
Das A, Kuznets-Speck B, Limmer DT. Direct Evaluation of Rare Events in Active Matter from Variational Path Sampling. PHYSICAL REVIEW LETTERS 2022; 128:028005. [PMID: 35089729 DOI: 10.1103/physrevlett.128.028005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/18/2021] [Accepted: 12/09/2021] [Indexed: 06/14/2023]
Abstract
Active matter represents a broad class of systems that evolve far from equilibrium due to the local injection of energy. Like their passive analogs, transformations between distinct metastable states in active matter proceed through rare fluctuations; however, their detailed balance violating dynamics renders these events difficult to study. Here, we present a simulation method for evaluating the rate and mechanism of rare events in generic nonequilibrium systems and apply it to study the conformational changes of a passive solute in an active fluid. The method employs a variational optimization of a control force that renders the rare event a typical one, supplying an exact estimate of its rate as a ratio of path partition functions. Using this method we find that increasing activity in the active bath can enhance the rate of conformational switching of the passive solute in a manner consistent with recent bounds from stochastic thermodynamics.
Collapse
Affiliation(s)
- Avishek Das
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - David T Limmer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawerence Berkeley National Laboratory, Berkeley, California 94720, USA
- Material Sciences Division, Lawerence Berkeley National Laboratory, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, University of California, Berkeley, California 94720, USA
| |
Collapse
|
4
|
Das A, Rose DC, Garrahan JP, Limmer DT. Reinforcement learning of rare diffusive dynamics. J Chem Phys 2021; 155:134105. [PMID: 34624994 DOI: 10.1063/5.0057323] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a method to probe rare molecular dynamics trajectories directly using reinforcement learning. We consider trajectories that are conditioned to transition between regions of configuration space in finite time, such as those relevant in the study of reactive events, and trajectories exhibiting rare fluctuations of time-integrated quantities in the long time limit, such as those relevant in the calculation of large deviation functions. In both cases, reinforcement learning techniques are used to optimize an added force that minimizes the Kullback-Leibler divergence between the conditioned trajectory ensemble and a driven one. Under the optimized added force, the system evolves the rare fluctuation as a typical one, affording a variational estimate of its likelihood in the original trajectory ensemble. Low variance gradients employing value functions are proposed to increase the convergence of the optimal force. The method we develop employing these gradients leads to efficient and accurate estimates of both the optimal force and the likelihood of the rare event for a variety of model systems.
Collapse
Affiliation(s)
- Avishek Das
- Department of Chemistry, University of California, Berkeley, California 94609, USA
| | - Dominic C Rose
- School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Juan P Garrahan
- School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - David T Limmer
- Department of Chemistry, University of California, Berkeley, California 94609, USA
| |
Collapse
|
5
|
Olla P. Description of a stochastic system by a nonadapted stochastic process. Phys Rev E 2021; 104:014139. [PMID: 34412222 DOI: 10.1103/physreve.104.014139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 07/08/2021] [Indexed: 11/07/2022]
Abstract
An approach for the description of stochastic systems is derived. Some of the variables in the system are studied forward in time, others backward in time. The approach is based on a perturbation expansion in the strength of the coupling between forward and backward variables, and it is well suited for situations in which initial and final conditions are imposed on different components of the system, and the coupling between those components is weak. The form of the stochastic equations in our approach is determined by requiring that they generate the same statistics obtained in a forward description of the dynamics. Numerical tests are carried out on a few simple two-degrees-of-freedom systems. The merit and the difficulties of the approach are discussed and compared to more traditional strategies based on transition path sampling and simple shooting algorithms.
Collapse
Affiliation(s)
- Piero Olla
- ISAC-CNR and INFN, Sez. Cagliari, I-09042 Monserrato, Italy
| |
Collapse
|
6
|
Bartolucci G, Orioli S, Faccioli P. Transition path theory from biased simulations. J Chem Phys 2018; 149:072336. [PMID: 30134709 DOI: 10.1063/1.5027253] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Transition Path Theory (TPT) provides a rigorous framework to investigate the dynamics of rare thermally activated transitions. In this theory, a central role is played by the forward committor function q+(x), which provides the ideal reaction coordinate. Furthermore, the reactive dynamics and kinetics are fully characterized in terms of two time-independent scalar and vector distributions. In this work, we develop a scheme which enables all these ingredients of TPT to be efficiently computed using the short non-equilibrium trajectories generated by means of a specific combination of enhanced path sampling techniques. In particular, first we further extend the recently introduced self-consistent path sampling algorithm in order to compute the committor q+(x). Next, we show how this result can be exploited in order to define efficient algorithms which enable us to directly sample the transition path ensemble.
Collapse
Affiliation(s)
- G Bartolucci
- Physics Department of Trento University, Via Sommarive 14, 37123 Povo (Trento), Italy
| | - S Orioli
- Physics Department of Trento University, Via Sommarive 14, 37123 Povo (Trento), Italy
| | - P Faccioli
- Physics Department of Trento University, Via Sommarive 14, 37123 Povo (Trento), Italy
| |
Collapse
|
7
|
Nguyen MK, Jaillet L, Redon S. Generating conformational transition paths with low potential-energy barriers for proteins. J Comput Aided Mol Des 2018; 32:853-867. [PMID: 30069648 DOI: 10.1007/s10822-018-0137-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2018] [Accepted: 07/19/2018] [Indexed: 10/28/2022]
Abstract
The knowledge of conformational transition paths in proteins can be useful for understanding protein mechanisms. Recently, we have introduced the As-Rigid-As-Possible (ARAP) interpolation method, for generating interpolation paths between two protein conformations. The method was shown to preserve well the rigidity of the initial conformation along the path. However, because the method is totally geometry-based, the generated paths may be inconsistent because the atom interactions are ignored. Therefore, in this article, we would like to introduce a new method to generate conformational transition paths with low potential-energy barriers for proteins. The method is composed of three processing stages. First, ARAP interpolation is used for generating an initial path. Then, the path conformations are enhanced by a clash remover. Finally, Nudged Elastic Band, a path-optimization method, is used to produce a low-energy path. Large energy reductions are found in the paths obtained from the method than in those obtained from the ARAP interpolation method alone. The results also show that ARAP interpolation is a good candidate for generating an initial path because it leads to lower potential-energy paths than two other common methods for path interpolation.
Collapse
Affiliation(s)
- Minh Khoa Nguyen
- Univ. Grenoble Alpes, Inria, CNRS, Grenoble INP (Institute of Engineering Univ. Grenoble Alpes), LJK, 38000, Grenoble, France
| | - Léonard Jaillet
- Univ. Grenoble Alpes, Inria, CNRS, Grenoble INP (Institute of Engineering Univ. Grenoble Alpes), LJK, 38000, Grenoble, France.
| | - Stéphane Redon
- Univ. Grenoble Alpes, Inria, CNRS, Grenoble INP (Institute of Engineering Univ. Grenoble Alpes), LJK, 38000, Grenoble, France
| |
Collapse
|