1
|
Del Razo MJ, Crommelin D, Bolhuis PG. Data-driven dynamical coarse-graining for condensed matter systems. J Chem Phys 2024; 160:024108. [PMID: 38193550 DOI: 10.1063/5.0177553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Accepted: 12/05/2023] [Indexed: 01/10/2024] Open
Abstract
Simulations of condensed matter systems often focus on the dynamics of a few distinguished components but require integrating the full system. A prime example is a molecular dynamics simulation of a (macro)molecule in a solution, where the molecule(s) and the solvent dynamics need to be integrated, rendering the simulations computationally costly and often unfeasible for physically/biologically relevant time scales. Standard coarse graining approaches can reproduce equilibrium distributions and structural features but do not properly include the dynamics. In this work, we develop a general data-driven coarse-graining methodology inspired by the Mori-Zwanzig formalism, which shows that macroscopic systems with a large number of degrees of freedom can be described by a few relevant variables and additional noise and memory terms. Our coarse-graining method consists of numerical integrators for the distinguished components, where the noise and interaction terms with other system components are substituted by a random variable sampled from a data-driven model. The model is parameterized using data from multiple short-time full-system simulations, and then, it is used to run long-time simulations. Applying our methodology to three systems-a distinguished particle under a harmonic and a bistable potential and a dimer with two metastable configurations-the resulting coarse-grained models are capable of reproducing not only the equilibrium distributions but also the dynamic behavior due to temporal correlations and memory effects. Remarkably, our method even reproduces the transition dynamics between metastable states, which is challenging to capture correctly. Our approach is not constrained to specific dynamics and can be extended to systems beyond Langevin dynamics, and, in principle, even to non-equilibrium dynamics.
Collapse
Affiliation(s)
- Mauricio J Del Razo
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090GD Amsterdam, The Netherlands
- Korteweg-de Vries Institute for Mathematics, University of Amsterdam, PO Box 94248, 1090GD Amsterdam, The Netherlands
- Dutch Institute for Emergent Phenomena, University of Amsterdam, Amsterdam, The Netherlands
| | - Daan Crommelin
- Korteweg-de Vries Institute for Mathematics, University of Amsterdam, PO Box 94248, 1090GD Amsterdam, The Netherlands
- Centrum Wiskunde & Informatica, 1098 XG Amsterdam, The Netherlands
| | - Peter G Bolhuis
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090GD Amsterdam, The Netherlands
| |
Collapse
|
2
|
Integrative modeling of the cell. Acta Biochim Biophys Sin (Shanghai) 2022; 54:1213-1221. [PMID: 36017893 PMCID: PMC9909318 DOI: 10.3724/abbs.2022115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A whole-cell model represents certain aspects of the cell structure and/or function. Due to the high complexity of the cell, an integrative modeling approach is often taken to utilize all available information including experimental data, prior knowledge and prior models. In this review, we summarize an emerging workflow of whole-cell modeling into five steps: (i) gather information; (ii) represent the modeled system into modules; (iii) translate input information into scoring function; (iv) sample the whole-cell model; (v) validate and interpret the model. In particular, we propose the integrative modeling of the cell by combining available (whole-cell) models to maximize the accuracy, precision, and completeness. In addition, we list quantitative predictions of various aspects of cell biology from existing whole-cell models. Moreover, we discuss the remaining challenges and future directions, and highlight the opportunity to establish an integrative spatiotemporal multi-scale whole-cell model based on a community approach.
Collapse
|
3
|
Floyd JE, Lukes JR. A neural network-assisted open boundary molecular dynamics simulation method. J Chem Phys 2022; 156:184114. [PMID: 35568556 DOI: 10.1063/5.0083198] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
A neural network-assisted molecular dynamics method is developed to reduce the computational cost of open boundary simulations. Particle influxes and neural network-derived forces are applied at the boundaries of an open domain consisting of explicitly modeled Lennard-Jones atoms in order to represent the effects of the unmodeled surrounding fluid. Canonical ensemble simulations with periodic boundaries are used to train the neural network and to sample boundary fluxes. The method, as implemented in the LAMMPS, yields temperature, kinetic energy, potential energy, and pressure values within 2.5% of those calculated using periodic molecular dynamics and runs two orders of magnitude faster than a comparable grand canonical molecular dynamics system.
Collapse
Affiliation(s)
- J E Floyd
- Department of Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - J R Lukes
- Department of Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
4
|
Del Razo MJ, Dibak M, Schütte C, Noé F. Multiscale molecular kinetics by coupling Markov state models and reaction-diffusion dynamics. J Chem Phys 2021; 155:124109. [PMID: 34598578 DOI: 10.1063/5.0060314] [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
A novel approach to simulate simple protein-ligand systems at large time and length scales is to couple Markov state models (MSMs) of molecular kinetics with particle-based reaction-diffusion (RD) simulations, MSM/RD. Currently, MSM/RD lacks a mathematical framework to derive coupling schemes, is limited to isotropic ligands in a single conformational state, and lacks multiparticle extensions. In this work, we address these needs by developing a general MSM/RD framework by coarse-graining molecular dynamics into hybrid switching diffusion processes. Given enough data to parameterize the model, it is capable of modeling protein-protein interactions over large time and length scales, and it can be extended to handle multiple molecules. We derive the MSM/RD framework, and we implement and verify it for two protein-protein benchmark systems and one multiparticle implementation to model the formation of pentameric ring molecules. To enable reproducibility, we have published our code in the MSM/RD software package.
Collapse
Affiliation(s)
- Mauricio J Del Razo
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - Manuel Dibak
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
| | | | - Frank Noé
- Department of Physics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
5
|
Tang W, Yu H, Zhao T, Qing L, Xu X, Zhao S. A dynamic reaction density functional theory for interfacial reaction-diffusion coupling at nanoscale. Chem Eng Sci 2021. [DOI: 10.1016/j.ces.2021.116513] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
6
|
Dibak M, Fröhner C, Noé F, Höfling F. Diffusion-influenced reaction rates in the presence of pair interactions. J Chem Phys 2019; 151:164105. [PMID: 31675872 DOI: 10.1063/1.5124728] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The kinetics of bimolecular reactions in solution depends, among other factors, on intermolecular forces such as steric repulsion or electrostatic interaction. Microscopically, a pair of molecules first has to meet by diffusion before the reaction can take place. In this work, we establish an extension of Doi's volume reaction model to molecules interacting via pair potentials, which is a key ingredient for interacting-particle-based reaction-diffusion (iPRD) simulations. As a central result, we relate model parameters and macroscopic reaction rate constants in this situation. We solve the corresponding reaction-diffusion equation in the steady state and derive semi-analytical expressions for the reaction rate constant and the local concentration profiles. Our results apply to the full spectrum from well-mixed to diffusion-limited kinetics. For limiting cases, we give explicit formulas, and we provide a computationally inexpensive numerical scheme for the general case, including the intermediate, diffusion-influenced regime. The obtained rate constants decompose uniquely into encounter and formation rates, and we discuss the effect of the potential on both subprocesses, exemplified for a soft harmonic repulsion and a Lennard-Jones potential. The analysis is complemented by extensive stochastic iPRD simulations, and we find excellent agreement with the theoretical predictions.
Collapse
Affiliation(s)
- Manuel Dibak
- Freie Universität Berlin, Fachbereich Mathematik und Informatik, Arnimallee 6, 14195 Berlin, Germany
| | - Christoph Fröhner
- Freie Universität Berlin, Fachbereich Mathematik und Informatik, Arnimallee 6, 14195 Berlin, Germany
| | - Frank Noé
- Freie Universität Berlin, Fachbereich Mathematik und Informatik, Arnimallee 6, 14195 Berlin, Germany
| | - Felix Höfling
- Freie Universität Berlin, Fachbereich Mathematik und Informatik, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
7
|
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
|