-
Force field optimization by imposing kinetic constraints with path reweighting
Authors:
P. G. Bolhuis,
Z. F. Brotzakis,
B. G. Keller
Abstract:
Empirical force fields employed in molecular dynamics simulations of complex systems can be optimised to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the rates of interconversion between metastable states in such systems, is hardly ever incorporated in a force field, due to a lack of an efficient approach. Here, we introduce…
▽ More
Empirical force fields employed in molecular dynamics simulations of complex systems can be optimised to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the rates of interconversion between metastable states in such systems, is hardly ever incorporated in a force field, due to a lack of an efficient approach. Here, we introduce such a framework, based on the relationship between dynamical observables such as rate constants, and the underlying force field parameters, using the statistical mechanics of trajectories. Given a prior ensemble of molecular trajectories produced with imperfect force field parameters, the approach allows the optimal adaption of these parameters, such that the imposed constraint of equal predicted and experimental rate constant is obeyed. To do so, the method combines the continuum path ensemble Maximum Caliber approach with path reweighting methods for stochastic dynamics. When multiple solutions are found, the method selects automatically the combination that corresponds to the smallest perturbation of the entire path ensemble, as required by the Maximum Entropy principle. To show the validity of the approach we illustrate the method on simple test systems undergoing rare event dynamics. Next to simple 2D potentials we explore particle models representing molecular isomerisation reactions as well as protein-ligand unbinding. Besides optimal interaction parameters the methodology gives physical insight into what parts of the model are most sensitive to the kinetics. We discuss the generality and broad implications of the methodology.
△ Less
Submitted 10 July, 2022;
originally announced July 2022.
-
A method of incorporating rate constants as kinetic constraints in molecular dynamics simulations
Authors:
Z. Faidon Brotzakis,
Michele Vendruscolo,
Peter. G. Bolhuis
Abstract:
From the point of view of statistical mechanics, a full characterisation of a molecular system requires the experimental determination of its possible states, their populations and the respective interconversion rates. Well-established methods can incorporate in molecular dynamics simulations experimental information about states using structural restraints, and about populations using thermodynam…
▽ More
From the point of view of statistical mechanics, a full characterisation of a molecular system requires the experimental determination of its possible states, their populations and the respective interconversion rates. Well-established methods can incorporate in molecular dynamics simulations experimental information about states using structural restraints, and about populations using thermodynamic restraints. However, it is still unclear how to include experimental knowledge of interconversion rates. Here we introduce a method of imposing known rate constants as constraints in molecular dynamics simulations, which is based on a combination of the maximum entropy and maximum caliber principles. Starting from an existing ensemble of trajectories, obtained from either molecular dynamics or enhanced trajectory sampling, this method provides a minimally perturbed path distribution consistent with the kinetic constraints, as well as a modified free energy and committor landscape. We illustrate the application of the method to simple toy systems, as well as to all atom molecular simulations of peptide association and folding. We find that by combining experimental rate coefficient data and molecular dynamics simulations we are able to determine new transition states, reaction mechanisms and free energies. For instance, in the case of chignolin protein folding we find that imposing a slower folding rate shifts the transition state to more native like conformations, while it increases the stability of the unfolded region.
△ Less
Submitted 1 June, 2020;
originally announced June 2020.
-
Approximating Free Energy and Committor Landscapes in Standard Transition Path Sampling using Virtual Interface Exchange
Authors:
Z. Faidon Brotzakis,
Peter G. Bolhuis
Abstract:
Transition path sampling (TPS) is a powerful technique for investigating rare transitions, especially when the mechanism is unknown and one does not have access to the reaction coordinate. Straightforward application of TPS does not directly provide the free energy landscape nor the kinetics, which motivated the development of path sampling extensions, such as transition interface sampling (TIS),…
▽ More
Transition path sampling (TPS) is a powerful technique for investigating rare transitions, especially when the mechanism is unknown and one does not have access to the reaction coordinate. Straightforward application of TPS does not directly provide the free energy landscape nor the kinetics, which motivated the development of path sampling extensions, such as transition interface sampling (TIS), and the reweighted paths ensemble (RPE), that are able to simultaneously access both kinetics and thermodynamics. However, performing TIS is more involved than TPS, and still requires (some) insight in the reaction to define interfaces. While packages that can efficiently compute path ensembles for TIS are now available, it would be useful to directly compute the free energy from a single TPS simulation. To achieve this, we developed an approximate method, denoted Virtual Interface Exchange, that makes use of the rejected pathways in a form of waste recycling. The method yields an approximate reweighted path ensemble that allows an immediate view of the free energy landscape from a single TPS, as well as enables a full committor analysis.
△ Less
Submitted 9 July, 2019;
originally announced July 2019.
-
Augmented Harmonic Linear Discriminant Analysis
Authors:
Z. Faidon Brotzakis,
Dan Mendels,
Michele Parrinello
Abstract:
Many processes of scientific and technological interest are characterized by time scales that render their simulation impossible if one uses present day simulation capabilities. To overcome this challenge a variety of enhanced simulation methods has been developed. A much-used class of methods relies on the use of collective variables. The efficiency of these methods relies critically on an educat…
▽ More
Many processes of scientific and technological interest are characterized by time scales that render their simulation impossible if one uses present day simulation capabilities. To overcome this challenge a variety of enhanced simulation methods has been developed. A much-used class of methods relies on the use of collective variables. The efficiency of these methods relies critically on an educated guess of the collective variables. For this reason much effort has been devoted to the construction and improvement of collective variables. Among the many methods proposed, harmonic linear discriminant analysis has proven effective. This method builds the collective coordinates solely from the knowledge of the fluctuations in the different metastable state. In this Letter we propose to improve upon the harmonic linear discriminant analysis by adding to the construction of the collective coordinates an extra bit of information, namely that of the transition state ensemble. Configurations belonging to the transition state ensemble are harnessed by the use of the spring shooting transition path sampling algorithm. We show on a challenging example that these coordinates thus augmented not only perform better in the calculation of the static properties, but also accelerate considerably the calculation of reaction rates.
△ Less
Submitted 23 February, 2019;
originally announced February 2019.