1
|
Blow KE, Tribello GA, Sosso GC, Quigley D. Interplay of multiple clusters and initial interface positioning for forward flux sampling simulations of crystal nucleation. J Chem Phys 2023; 158:2895225. [PMID: 37290068 DOI: 10.1063/5.0152343] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Accepted: 05/19/2023] [Indexed: 06/10/2023] Open
Abstract
Forward flux sampling (FFS) is a path sampling technique widely used in computer simulations of crystal nucleation from the melt. In such studies, the order parameter underpinning the progress of the FFS algorithm is often the size of the largest crystalline nucleus. In this work, we investigate the effects of two computational aspects of FFS simulations, using the prototypical Lennard-Jones liquid as our computational test bed. First, we quantify the impact of the positioning of the liquid basin and first interface in the space of the order parameter. In particular, we demonstrate that these choices are key to ensuring the consistency of the FFS results. Second, we focus on the frequently encountered scenario where the population of crystalline nuclei is such that there are multiple clusters of size comparable to the largest one. We demonstrate the contribution of clusters other than the largest cluster to the initial flux; however, we also show that they can be safely ignored for the purposes of converging a full FFS calculation. We also investigate the impact of different clusters merging, a process that appears to be facilitated by substantial spatial correlations-at least at the supercooling considered here. Importantly, all of our results have been obtained as a function of system size, thus contributing to the ongoing discussion on the impact of finite size effects on simulations of crystal nucleation. Overall, this work either provides or justifies several practical guidelines for performing FFS simulations that can also be applied to more complex and/or computationally expensive models.
Collapse
Affiliation(s)
- Katarina E Blow
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Gareth A Tribello
- Centre for Quantum Materials and Technologies, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - David Quigley
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
2
|
Reiner M, Bachmair B, Tiefenbacher MX, Mai S, González L, Marquetand P, Dellago C. Nonadiabatic Forward Flux Sampling for Excited-State Rare Events. J Chem Theory Comput 2023; 19:1657-1671. [PMID: 36856706 PMCID: PMC10061683 DOI: 10.1021/acs.jctc.2c01088] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Indexed: 03/02/2023]
Abstract
We present a rare event sampling scheme applicable to coupled electronic excited states. In particular, we extend the forward flux sampling (FFS) method for rare event sampling to a nonadiabatic version (NAFFS) that uses the trajectory surface hopping (TSH) method for nonadiabatic dynamics. NAFFS is applied to two dynamically relevant excited-state models that feature an avoided crossing and a conical intersection with tunable parameters. We investigate how nonadiabatic couplings, temperature, and reaction barriers affect transition rate constants in regimes that cannot be otherwise obtained with plain, traditional TSH. The comparison with reference brute-force TSH simulations for limiting cases of rareness shows that NAFFS can be several orders of magnitude cheaper than conventional TSH and thus represents a conceptually novel tool to extend excited-state dynamics to time scales that are able to capture rare nonadiabatic events.
Collapse
Affiliation(s)
- Madlen
Maria Reiner
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Maximilian Xaver Tiefenbacher
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Leticia González
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Christoph Dellago
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Faculty
of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
Collapse
|
3
|
Hall SW, Díaz Leines G, Sarupria S, Rogal J. Practical guide to replica exchange transition interface sampling and forward flux sampling. J Chem Phys 2022; 156:200901. [DOI: 10.1063/5.0080053] [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
Path sampling approaches have become invaluable tools to explore the mechanisms and dynamics of the so-called rare events that are characterized by transitions between metastable states separated by sizable free energy barriers. Their practical application, in particular to ever more complex molecular systems, is, however, not entirely trivial. Focusing on replica exchange transition interface sampling (RETIS) and forward flux sampling (FFS), we discuss a range of analysis tools that can be used to assess the quality and convergence of such simulations, which is crucial to obtain reliable results. The basic ideas of a step-wise evaluation are exemplified for the study of nucleation in several systems with different complexities, providing a general guide for the critical assessment of RETIS and FFS simulations.
Collapse
Affiliation(s)
- Steven W. Hall
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Sapna Sarupria
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| | - Jutta Rogal
- Department of Chemistry, New York University, New York, New York 10003, USA
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
4
|
Berkov D, Semenova EK, Gorn NL. Concept of the Energy-Dependent Temperature for Direct Langevin Dynamics Simulations of Rare Events in Systems with Arbitrary High Energy Barriers. PHYSICAL REVIEW LETTERS 2021; 127:247201. [PMID: 34951811 DOI: 10.1103/physrevlett.127.247201] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Revised: 10/22/2021] [Accepted: 11/10/2021] [Indexed: 06/14/2023]
Abstract
We suggest an algorithm which allows single-stage direct Langevin dynamics simulations of transitions over arbitrarily high energy barriers. For this purpose, we propose a concept of the energy-dependent temperature (EDT): near the energy minima this temperature is high, but it tends toward room temperature for energies approaching the barrier value. In the resulting algorithm simulation time required for the computation of the escape rate over the barrier does not increase with barrier height. Switching times computed via our EDT algorithm agree very well with those obtained with the forward flux sampling (FFS). As the simulation time required by our method does not increase with the energy barrier, we achieve a very large speed-up compared even to our strongly optimized version of FFS (and all other multistage algorithms). In addition, our approach is free from the instability occurring in all multistage "climbing" methods where a product of a large number of transition probabilities between the interfaces must be computed.
Collapse
Affiliation(s)
- Dmitry Berkov
- General Numerics Research Lab, Leutragraben 1, D-07743 Jena, Germany
| | - Elena K Semenova
- General Numerics Research Lab, Leutragraben 1, D-07743 Jena, Germany
| | - Natalia L Gorn
- General Numerics Research Lab, Leutragraben 1, D-07743 Jena, Germany
| |
Collapse
|
5
|
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
|
6
|
Albaugh A, Gingrich TR. Estimating reciprocal partition functions to enable design space sampling. J Chem Phys 2020; 153:204102. [PMID: 33261473 DOI: 10.1063/5.0025358] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Reaction rates are a complicated function of molecular interactions, which can be selected from vast chemical design spaces. Seeking the design that optimizes a rate is a particularly challenging problem since the rate calculation for any one design is itself a difficult computation. Toward this end, we demonstrate a strategy based on transition path sampling to generate an ensemble of designs and reactive trajectories with a preference for fast reaction rates. Each step of the Monte Carlo procedure requires a measure of how a design constrains molecular configurations, expressed via the reciprocal of the partition function for the design. Although the reciprocal of the partition function would be prohibitively expensive to compute, we apply Booth's method for generating unbiased estimates of a reciprocal of an integral to sample designs without bias. A generalization with multiple trajectories introduces a stronger preference for fast rates, pushing the sampled designs closer to the optimal design. We illustrate the methodology on two toy models of increasing complexity: escape of a single particle from a Lennard-Jones potential well of tunable depth and escape from a metastable tetrahedral cluster with tunable pair potentials.
Collapse
Affiliation(s)
- Alex Albaugh
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Todd R Gingrich
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| |
Collapse
|
7
|
Hussain S, Haji-Akbari A. Studying rare events using forward-flux sampling: Recent breakthroughs and future outlook. J Chem Phys 2020; 152:060901. [DOI: 10.1063/1.5127780] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Sarwar Hussain
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
8
|
Klein MC, Roberts E. Automatic error control during forward flux sampling of rare events in master equation models. J Chem Phys 2020; 152:035102. [PMID: 31968949 DOI: 10.1063/1.5129461] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Enhanced sampling methods, such as forward flux sampling (FFS), have great capacity for accelerating stochastic simulations of nonequilibrium biochemical systems involving rare events. However, the description of the tradeoffs between simulation efficiency and error in FFS remains incomplete. We present a novel and mathematically rigorous analysis of the errors in FFS that, for the first time, covers the contribution of every phase of the simulation. We derive a closed form expression for the optimally efficient count of samples to take in each FFS phase in terms of a fixed constraint on sampling error. We introduce a new method, forward flux pilot sampling (FFPilot), that is designed to take full advantage of our optimizing equation without prior information or assumptions about the phase weights and costs along the transition path. In simulations of both single and multidimensional gene regulatory networks, FFPilot is able to completely control sampling error. We then discuss how memory effects can introduce additional error when relaxation along the transition path is slow. This extra error can be traced to correlations between the FFS phases and can be controlled by monitoring the covariance between them. Finally, we show that, in sets of simulations with matched error, FFPilot is on the order of tens-to-hundreds of times faster than direct sampling and noticeably more efficient than previous FFS methods.
Collapse
Affiliation(s)
- Max C Klein
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Elijah Roberts
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
9
|
Affiliation(s)
- Hossein Rezvantalab
- Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Ronald G. Larson
- Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
10
|
Mökkönen H, Ala-Nissila T, Jónsson H. Efficient dynamical correction of the transition state theory rate estimate for a flat energy barrier. J Chem Phys 2016; 145:094901. [DOI: 10.1063/1.4962167] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Harri Mökkönen
- Department of Applied Physics, Aalto University School of Science, P.O. Box 1100, FIN-00076 Aalto, Espoo, Finland
- COMP CoE, Aalto University School of Science, P.O. Box 1100, FIN-00076 Aalto, Espoo, Finland
- Faculty of Physical Sciences and Science Institute, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Tapio Ala-Nissila
- Department of Applied Physics, Aalto University School of Science, P.O. Box 1100, FIN-00076 Aalto, Espoo, Finland
- COMP CoE, Aalto University School of Science, P.O. Box 1100, FIN-00076 Aalto, Espoo, Finland
- Department of Physics, Brown University, Providence Rhode Island 02912-1843, USA
| | - Hannes Jónsson
- Department of Applied Physics, Aalto University School of Science, P.O. Box 1100, FIN-00076 Aalto, Espoo, Finland
- Faculty of Physical Sciences and Science Institute, University of Iceland, VR-III, 107 Reykjavík, Iceland
| |
Collapse
|
11
|
Jungblut S, Dellago C. Pathways to self-organization: Crystallization via nucleation and growth. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2016; 39:77. [PMID: 27498980 DOI: 10.1140/epje/i2016-16077-6] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Accepted: 07/13/2016] [Indexed: 06/06/2023]
Abstract
Crystallization, a prototypical self-organization process during which a disordered state spontaneously transforms into a crystal characterized by a regular arrangement of its building blocks, usually proceeds by nucleation and growth. In the initial stages of the transformation, a localized nucleus of the new phase forms in the old one due to a random fluctuation. Most of these nuclei disappear after a short time, but rarely a crystalline embryo may reach a critical size after which further growth becomes thermodynamically favorable and the entire system is converted into the new phase. In this article, we will discuss several theoretical concepts and computational methods to study crystallization. More specifically, we will address the rare event problem arising in the simulation of nucleation processes and explain how to calculate nucleation rates accurately. Particular attention is directed towards discussing statistical tools to analyze crystallization trajectories and identify the transition mechanism.
Collapse
Affiliation(s)
- S Jungblut
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Wien, Austria
| | - C Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090, Wien, Austria.
| |
Collapse
|
12
|
Thapar V, Escobedo FA. Simultaneous estimation of free energies and rates using forward flux sampling and mean first passage times. J Chem Phys 2015; 143:244113. [DOI: 10.1063/1.4938248] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Vikram Thapar
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| | - Fernando A. Escobedo
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
13
|
|
14
|
Kratzer K, Arnold A, Allen RJ. Automatic, optimized interface placement in forward flux sampling simulations. J Chem Phys 2013; 138:164112. [DOI: 10.1063/1.4801866] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
|
15
|
Becker NB, Allen RJ, ten Wolde PR. Non-stationary forward flux sampling. J Chem Phys 2012; 136:174118. [DOI: 10.1063/1.4704810] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
16
|
Abstract
The drying of hydrophobic cavities is believed to play an important role in biophysical phenomena such as the folding of globular proteins, the opening and closing of ligand-gated ion channels, and ligand binding to hydrophobic pockets. We use forward flux sampling, a molecular simulation technique, to compute the rate of capillary evaporation of water confined between two hydrophobic surfaces separated by nanoscopic gaps, as a function of gap, surface size, and temperature. Over the range of conditions investigated (gaps between 9 and 14 Å and surface areas between 1 and 9 nm(2)), the free energy barrier to evaporation scales linearly with the gap between hydrophobic surfaces, suggesting that line tension makes the predominant contribution to the free energy barrier. The exponential dependence of the evaporation rate on the gap between confining surfaces causes a 10 order-of-magnitude decrease in the rate when the gap increases from 9 to 14 Å. The computed free energy barriers are of the order of 50 kT and are predominantly enthalpic. Evaporation rates per unit area are found to be two orders of magnitude faster in confinement by the larger (9 nm(2)) than by the smaller (1 nm(2)) surfaces considered here, at otherwise identical conditions. We show that this rate enhancement is a consequence of the dependence of hydrophobic hydration on the size of solvated objects. For sufficiently large surfaces, the critical nucleus for the evaporation process is a gap-spanning vapor tube.
Collapse
|
17
|
Du WN, Marino KA, Bolhuis PG. Multiple state transition interface sampling of alanine dipeptide in explicit solvent. J Chem Phys 2011; 135:145102. [DOI: 10.1063/1.3644344] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Araque JC, Panagiotopoulos AZ, Robert MA. Lattice model of oligonucleotide hybridization in solution. I. Model and thermodynamics. J Chem Phys 2011; 134:165103. [PMID: 21528982 DOI: 10.1063/1.3568145] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
A coarse-grained lattice model of DNA oligonucleotides is proposed to investigate the general mechanisms by which single-stranded oligonucleotides hybridize to their complementary strands in solution. The model, based on a high-coordination cubic lattice, is simple enough to allow the direct simulation of DNA solutions, yet capturing how the fundamental thermodynamic processes are microscopically encoded in the nucleobase sequences. Physically relevant interactions are considered explicitly, such as interchain excluded volume, anisotropic base-pairing and base-stacking, and single-stranded bending rigidity. The model is studied in detail by a specially adapted Monte Carlo simulation method, based on parallel tempering and biased trials, which is designed to overcome the entropic and enthalpic barriers associated with the sampling of hybridization events of multiple single-stranded chains in solution. This methodology addresses both the configurational complexity of bringing together two complementary strands in a favorable orientation (entropic barrier) and the energetic penalty of breaking apart multiple associated bases in a double-stranded state (enthalpic barrier). For strands with sequences restricted to nonstaggering association and homogeneous pairing and stacking energies, base-pairing is found to dominate the hybridization over the translational and conformational entropy. For strands with sequence-dependent pairing corresponding to that of DNA, the complex dependence of the model's thermal stability on concentration, sequence, and degree of complementarity is shown to be qualitatively and quantitatively consistent both with experiment and with the predictions of statistical mechanical models.
Collapse
Affiliation(s)
- Juan C Araque
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, USA
| | | | | |
Collapse
|
19
|
Borrero EE, Weinwurm M, Dellago C. Optimizing transition interface sampling simulations. J Chem Phys 2011; 134:244118. [DOI: 10.1063/1.3601919] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
20
|
Schlick T, Collepardo-Guevara R, Halvorsen LA, Jung S, Xiao X. Biomolecularmodeling and simulation: a field coming of age. Q Rev Biophys 2011; 44:191-228. [PMID: 21226976 PMCID: PMC3700731 DOI: 10.1017/s0033583510000284] [Citation(s) in RCA: 100] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We assess the progress in biomolecular modeling and simulation, focusing on structure prediction and dynamics, by presenting the field’s history, metrics for its rise in popularity, early expressed expectations, and current significant applications. The increases in computational power combined with improvements in algorithms and force fields have led to considerable success, especially in protein folding, specificity of ligand/biomolecule interactions, and interpretation of complex experimental phenomena (e.g. NMR relaxation, protein-folding kinetics and multiple conformational states) through the generation of structural hypotheses and pathway mechanisms. Although far from a general automated tool, structure prediction is notable for proteins and RNA that preceded the experiment, especially by knowledge-based approaches. Thus, despite early unrealistic expectations and the realization that computer technology alone will not quickly bridge the gap between experimental and theoretical time frames, ongoing improvements to enhance the accuracy and scope of modeling and simulation are propelling the field onto a productive trajectory to become full partner with experiment and a field on its own right.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, 100 Washington Square East, Silver Building, New York, NY 10003, USA.
| | | | | | | | | |
Collapse
|
21
|
Zwier MC, Kaus JW, Chong LT. Efficient Explicit-Solvent Molecular Dynamics Simulations of Molecular Association Kinetics: Methane/Methane, Na+/Cl−, Methane/Benzene, and K+/18-Crown-6 Ether. J Chem Theory Comput 2011; 7:1189-97. [DOI: 10.1021/ct100626x] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Matthew C. Zwier
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Joseph W. Kaus
- 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
|
22
|
Borrero EE, Dellago C. Overcoming barriers in trajectory space: mechanism and kinetics of rare events via Wang-Landau enhanced transition path sampling. J Chem Phys 2011; 133:134112. [PMID: 20942528 DOI: 10.1063/1.3496376] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Within the framework of transition path sampling (TPS), activation energies can be computed as path ensemble averages without a priori information about the reaction mechanism [C. Dellago and P. G. Bolhuis, Mol. Simul. 30, 795 (2004)]. Activation energies computed for different conditions can then be used to determine by numerical integration the rate constant for a system of interest from the rate constant known for a reference system. However, in systems with complex potential energy surfaces, multiple reaction pathways may exist making ergodic sampling of trajectory space difficult. Here, we present a combination of TPS with the Wang-Landau (WL) flat-histogram algorithm for an efficient sampling of the transition path ensemble. This method, denoted by WL-TPS, has the advantage that from one single simulation, activation energies at different temperatures can be determined even for systems with multiple reaction mechanisms. The proposed methodology for rate constant calculations does not require the knowledge of the reaction coordinate and is generally applicable to Arrhenius and non-Arrhenius processes. We illustrate the applicability of this technique by studying a two-dimensional toy system consisting of a triatomic molecule immersed in a fluid of repulsive soft disks. We also provide an expression for the calculation of activation volumes from path averages such that the pressure dependence of the rate constant can be obtained by numerical integration.
Collapse
Affiliation(s)
- Ernesto E Borrero
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | | |
Collapse
|
23
|
Velez-Vega C, Borrero EE, Escobedo FA. Kinetics and mechanism of the unfolding native-to-loop transition of Trp-cage in explicit solvent via optimized forward flux sampling simulations. J Chem Phys 2011; 133:105103. [PMID: 20849192 DOI: 10.1063/1.3474803] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The native-to-loop (N-L) unfolding transition of Trp-cage protein was studied via optimized forward flux sampling (FFS) methods with trajectories evolved using molecular dynamics. The rate constant calculated from our simulations is in good agreement with the experimental value for the native-to-unfolded transition of this protein; furthermore, the trajectories sampled a phase region consistent with that reported in previous studies for the N-L transition using transition path sampling and transition interface sampling. A new variant of FFS is proposed and implemented that allows a better control of a constant flux of partial paths. A reaction coordinate model was obtained, at no extra cost, from the transition path ensemble generated by FFS, through iterative use of the FFS-least-square estimation method [E. E. Borrero and F. A. Escobedo, J. Chem. Phys. 127, 164101 (2007)] and an adaptive staging optimization algorithm [E. E. Borrero and F. A. Escobedo, J. Chem. Phys. 129, 024115 (2008)]. Finally, we further elucidate the unfolding mechanism by correlating the unfolding progress with changes in the root mean square deviation from the α carbons of the native state, the root mean square deviation from an ideal α-helix, and other structural properties of the protein.
Collapse
Affiliation(s)
- Camilo Velez-Vega
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| | | | | |
Collapse
|
24
|
Adams DA, Ziff RM, Sander LM. Computation of nucleation at a nonequilibrium first-order phase transition using a rare-event algorithm. J Chem Phys 2010; 133:174107. [DOI: 10.1063/1.3499321] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Dickson A, Warmflash A, Dinner AR. Separating forward and backward pathways in nonequilibrium umbrella sampling. J Chem Phys 2010; 131:154104. [PMID: 20568844 DOI: 10.1063/1.3244561] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Umbrella sampling enforces uniform sampling of steady-state distributions that are functions of arbitrary numbers of order parameters. The key to applying such methods to nonequilibrium processes is the accumulation of fluxes between regions. A significant difference between microscopically reversible and irreversible systems is that, in the latter case, the transition path ensemble for a reaction can be significantly different for "forward" and "backward" trajectories. Here, we show how to separately treat forward and backward pathways in nonequilibrium umbrella sampling simulations by working in an extended space. In this extended space, the exact rate (for equilibrium or nonequilibrium processes) can be calculated "for free" as a flux in phase space. We compare the efficiency of this rate calculation with forward flux sampling for a two-dimensional potential and show that nonequilibrium umbrella sampling is more efficient when an intermediate is present. We show that this technique can also be used to describe steady-state limit cycles by examining a simulation of circadian oscillations. We obtain the path of the limit cycle in a space of 22 order parameters, as well as the oscillation period. The relation of our method to others is discussed.
Collapse
Affiliation(s)
- Alex Dickson
- James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | | | | |
Collapse
|
26
|
Adams DA, Sander LM, Ziff RM. The barrier method: A technique for calculating very long transition times. J Chem Phys 2010; 133:124103. [DOI: 10.1063/1.3485285] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
27
|
Ryu S, Cai W. Numerical tests of nucleation theories for the Ising models. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:011603. [PMID: 20866625 DOI: 10.1103/physreve.82.011603] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2010] [Revised: 06/21/2010] [Indexed: 05/29/2023]
Abstract
The classical nucleation theory (CNT) is tested systematically by computer simulations of the two-dimensional (2D) and three-dimensional (3D) Ising models with a Glauber-type spin flip dynamics. While previous studies suggested potential problems with CNT, our numerical results show that the fundamental assumption of CNT is correct. In particular, the Becker-Döring theory accurately predicts the nucleation rate if the correct droplet free energy function is provided as input. This validates the coarse graining of the system into a one dimensional Markov chain with the largest droplet size as the reaction coordinate. Furthermore, in the 2D Ising model, the droplet free energy predicted by CNT matches numerical results very well, after a logarithmic correction term from Langer's field theory and a constant correction term are added. But significant discrepancies are found between the numerical results and existing theories on the magnitude of the logarithmic correction term in the 3D Ising model. Our analysis underscores the importance of correctly accounting for the temperature dependence of surface energy when comparing numerical results and nucleation theories.
Collapse
Affiliation(s)
- Seunghwa Ryu
- Department of Physics, Stanford University, Stanford, California 94305, USA
| | | |
Collapse
|
28
|
Borrero EE, Contreras Martínez LM, DeLisa MP, Escobedo FA. Kinetics and reaction coordinates of the reassembly of protein fragments via forward flux sampling. Biophys J 2010; 98:1911-20. [PMID: 20441755 PMCID: PMC2862158 DOI: 10.1016/j.bpj.2009.12.4329] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2009] [Revised: 12/07/2009] [Accepted: 12/15/2009] [Indexed: 11/16/2022] Open
Abstract
We studied the mechanism of the reassembly and folding process of two fragments of a split lattice protein by using forward flux sampling (FFS). Our results confirmed previous thermodynamics and kinetics analyses that suggested that the disruption of the critical core (of an unsplit protein that folds by a nucleation mechanism) plays a key role in the reassembly mechanism of the split system. For several split systems derived from a parent 48-mer model, we estimated the reaction coordinates in terms of collective variables by using the FFS least-square estimation method and found that the reassembly transition is best described by a combination of the total number of native contacts, the number of interchain native contacts, and the total conformational energy of the split system. We also analyzed the transition path ensemble obtained from FFS simulations using the estimated reaction coordinates as order parameters to identify the microscopic features that differentiate the reassembly of the different split systems studied. We found that in the fastest folding split system, a balanced distribution of the original-core amino acids (of the unsplit system) between protein fragments propitiates interchain interactions at early stages of the folding process. Only this system exhibits a different reassembly mechanism from that of the unsplit protein, involving the formation of a different folding nucleus. In the slowest folding system, the concentration of the folding nucleus in one fragment causes its early prefolding, whereas the second fragment tends to remain as a detached random coil. We also show that the reassembly rate can be either increased or decreased by tuning interchain cooperativeness via the introduction of a single point mutation that either strengthens or weakens one of the native interchain contacts (prevalent in the transition state ensemble).
Collapse
Affiliation(s)
| | | | | | - Fernando A. Escobedo
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York
| |
Collapse
|
29
|
Affiliation(s)
- Alex Dickson
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637;
| | - Aaron R. Dinner
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637;
| |
Collapse
|
30
|
Zhang BW, Jasnow D, Zuckerman DM. The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures. J Chem Phys 2010; 132:054107. [PMID: 20136305 PMCID: PMC2830257 DOI: 10.1063/1.3306345] [Citation(s) in RCA: 133] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2009] [Accepted: 01/12/2010] [Indexed: 11/14/2022] Open
Abstract
The "weighted ensemble" method, introduced by Huber and Kim [Biophys. J. 70, 97 (1996)], is one of a handful of rigorous approaches to path sampling of rare events. Expanding earlier discussions, we show that the technique is statistically exact for a wide class of Markovian and non-Markovian dynamics. The derivation is based on standard path-integral (path probability) ideas, but recasts the weighted-ensemble approach as simple "resampling" in path space. Similar reasoning indicates that arbitrary nonstatic binning procedures, which merely guide the resampling process, are also valid. Numerical examples confirm the claims, including the use of bins which can adaptively find the target state in a simple model.
Collapse
Affiliation(s)
- Bin W Zhang
- Department of Computational Biology, School of Medicine, University of Pittsburgh, Pennsylvania 15260, USA
| | | | | |
Collapse
|
31
|
Allen RJ, Valeriani C, Rein Ten Wolde P. Forward flux sampling for rare event simulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:463102. [PMID: 21715864 DOI: 10.1088/0953-8984/21/46/463102] [Citation(s) in RCA: 229] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Rare events are ubiquitous in many different fields, yet they are notoriously difficult to simulate because few, if any, events are observed in a conventional simulation run. Over the past several decades, specialized simulation methods have been developed to overcome this problem. We review one recently developed class of such methods, known as forward flux sampling. Forward flux sampling uses a series of interfaces between the initial and final states to calculate rate constants and generate transition paths for rare events in equilibrium or nonequilibrium systems with stochastic dynamics. This review draws together a number of recent advances, summarizes several applications of the method and highlights challenges that remain to be overcome.
Collapse
Affiliation(s)
- Rosalind J Allen
- SUPA, School of Physics and Astronomy, University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, UK
| | | | | |
Collapse
|
32
|
Zheng L, Chen M, Yang W. Simultaneous escaping of explicit and hidden free energy barriers: application of the orthogonal space random walk strategy in generalized ensemble based conformational sampling. J Chem Phys 2009; 130:234105. [PMID: 19548709 DOI: 10.1063/1.3153841] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
To overcome the pseudoergodicity problem, conformational sampling can be accelerated via generalized ensemble methods, e.g., through the realization of random walks along prechosen collective variables, such as spatial order parameters, energy scaling parameters, or even system temperatures or pressures, etc. As usually observed, in generalized ensemble simulations, hidden barriers are likely to exist in the space perpendicular to the collective variable direction and these residual free energy barriers could greatly abolish the sampling efficiency. This sampling issue is particularly severe when the collective variable is defined in a low-dimension subset of the target system; then the "Hamiltonian lagging" problem, which reveals the fact that necessary structural relaxation falls behind the move of the collective variable, may be likely to occur. To overcome this problem in equilibrium conformational sampling, we adopted the orthogonal space random walk (OSRW) strategy, which was originally developed in the context of free energy simulation [L. Zheng, M. Chen, and W. Yang, Proc. Natl. Acad. Sci. U.S.A. 105, 20227 (2008)]. Thereby, generalized ensemble simulations can simultaneously escape both the explicit barriers along the collective variable direction and the hidden barriers that are strongly coupled with the collective variable move. As demonstrated in our model studies, the present OSRW based generalized ensemble treatments show improved sampling capability over the corresponding classical generalized ensemble treatments.
Collapse
Affiliation(s)
- Lianqing Zheng
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, USA
| | | | | |
Collapse
|
33
|
Escobedo FA, Borrero EE, Araque JC. Transition path sampling and forward flux sampling. Applications to biological systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:333101. [PMID: 21828593 DOI: 10.1088/0953-8984/21/33/333101] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The last decade has seen a rapid growth in the number of simulation methods and applications dealing with the sampling of transition pathways of rare nanoscale events. Such studies are crucial, for example, for understanding the mechanism and kinetics of conformational transitions and enzymatic events associated with the function of biomolecules. In this review, a broad account of transition path sampling approaches is provided, starting from the general concepts, progressing to the specific principles that underlie some of the most important methods, and eventually singling out the so-called forward flux sampling method for a more detailed description. This is done because forward flux sampling, despite its appealing simplicity and potential efficiency, has thus far received limited attention from practitioners. While path sampling methods have a widespread application to many types of rare transitional events, here only recent applications involving biomolecules are reviewed, including isomerization, protein folding, and enzyme catalysis.
Collapse
|
34
|
Schlick T. Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules. F1000 BIOLOGY REPORTS 2009; 1:51. [PMID: 20948633 PMCID: PMC2948272 DOI: 10.3410/b1-51] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
The rugged energy landscape of biomolecules together with shortcomings of traditional molecular dynamics (MD) simulations require specialized methods for capturing large-scale, long-time configurational changes along with chemical dynamics behavior. In this report, MD-based methods for biomolecules are surveyed, involving modification of the potential, simulation protocol, or algorithm as well as global reformulations. While many of these methods are successful at probing the thermally accessible configuration space at the expense of altered kinetics, more sophisticated approaches like transition path sampling or Markov chain models are required to obtain mechanistic information, reaction pathways, and/or reaction rates. Divide-and-conquer methods for sampling and for piecing together reaction rate information are especially suitable for readily available computer cluster networks. Successful applications to biomolecules remain a challenge.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University 251 Mercer Street, New York, NY 10012 USA.
| |
Collapse
|
35
|
Velez-Vega C, Borrero EE, Escobedo FA. Kinetics and reaction coordinate for the isomerization of alanine dipeptide by a forward flux sampling protocol. J Chem Phys 2009; 130:225101. [DOI: 10.1063/1.3147465] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
36
|
Borrero EE, Escobedo FA. Simulating the Kinetics and Thermodynamics of Transitions via Forward Flux/Umbrella Sampling. J Phys Chem B 2009; 113:6434-45. [DOI: 10.1021/jp809103k] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ernesto E. Borrero
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853
| | - Fernando A. Escobedo
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853
| |
Collapse
|