@article{KieningerDonatiKeller, author = {Kieninger, Stefanie and Donati, Luca and Keller, Bettina G.}, title = {Dynamical reweighting methods for Markov models}, series = {Current Opinion in Structural Biology}, volume = {61}, journal = {Current Opinion in Structural Biology}, doi = {10.1016/j.sbi.2019.12.018}, pages = {124 -- 131}, abstract = {Conformational dynamics is essential to biomolecular processes. Markov State Models (MSMs) are widely used to elucidate dynamic properties of molecular systems from unbiased Molecular Dynamics (MD). However, the implementation of reweighting schemes for MSMs to analyze biased simulations is still at an early stage of development. Several dynamical reweighing approaches have been proposed, which can be classified as approaches based on (i) Kramers rate theory, (ii) rescaling of the probability density flux, (iii) reweighting by formulating a likelihood function, (iv) path reweighting. We present the state-of-the-art and discuss the methodological differences of these methods, their limitations and recent applications.}, language = {en} } @article{QuerDonatiKelleretal., author = {Quer, Jannes and Donati, Luca and Keller, Bettina and Weber, Marcus}, title = {An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates}, series = {SIAM Journal on Scientific Computing}, volume = {40}, journal = {SIAM Journal on Scientific Computing}, number = {2}, doi = {10.1137/17m1124772}, pages = {A653 -- A670}, abstract = {In this article we propose an adaptive importance sampling scheme for dynamical quantities of high dimensional complex systems which are metastable. The main idea of this article is to combine a method coming from Molecular Dynamics Simulation, Metadynamics, with a theorem from stochastic analysis, Girsanov's theorem. The proposed algorithm has two advantages compared to a standard estimator of dynamic quantities: firstly, it is possible to produce estimators with a lower variance and, secondly, we can speed up the sampling. One of the main problems for building importance sampling schemes for metastable systems is to find the metastable region in order to manipulate the potential accordingly. Our method circumvents this problem by using an assimilated version of the Metadynamics algorithm and thus creates a non-equilibrium dynamics which is used to sample the equilibrium quantities.}, language = {en} } @article{DonatiHartmannKeller, author = {Donati, Luca and Hartmann, Carsten and Keller, Bettina G.}, title = {Girsanov reweighting for path ensembles and Markov state models}, series = {The Journal of Chemical Physics}, volume = {146}, journal = {The Journal of Chemical Physics}, number = {24}, doi = {10.1063/1.4989474}, pages = {244112}, abstract = {The sensitivity of molecular dynamics on changes in the potential energy function plays an important role in understanding the dynamics and function of complex molecules. We present a method to obtain path ensemble averages of a perturbed dynamics from a set of paths generated by a reference dynamics. It is based on the concept of path probability measure and the Girsanov theorem, a result from stochastic analysis to estimate a change of measure of a path ensemble. Since Markov state models (MSMs) of the molecular dynamics can be formulated as a combined phase-space and path ensemble average, the method can be extended to reweight MSMs by combining it with a reweighting of the Boltzmann distribution. We demonstrate how to efficiently implement the Girsanov reweighting in a molecular dynamics simulation program by calculating parts of the reweighting factor "on the fly" during the simulation, and we benchmark the method on test systems ranging from a two-dimensional diffusion process and an artificial many-body system to alanine dipeptide and valine dipeptide in implicit and explicit water. The method can be used to study the sensitivity of molecular dynamics on external perturbations as well as to reweight trajectories generated by enhanced sampling schemes to the original dynamics.}, language = {en} } @article{GhysbrechtDonatiKeller, author = {Ghysbrecht, Simon and Donati, Luca and Keller, Bettina G.}, title = {Accuracy of reaction coordinate based rate theories for modelling chemical reactions: insights from the thermal isomerization in retinal}, series = {Submitted to The Journal of Physical Chemistry A}, journal = {Submitted to The Journal of Physical Chemistry A}, abstract = {Modern potential energy surfaces have shifted attention to molecular simulations of chemical reactions. While various methods can estimate rate constants for conformational transitions in molecular dynamics simulations, their applicability to studying chemical reactions remains uncertain due to the high and sharp energy barriers and complex reaction coordinates involved. This study focuses on the thermal cis-trans isomerization in retinal, employing molecular simulations and comparing rate constant estimates based on one-dimensional rate theories with those based on sampling transitions and grid-based models for low-dimensional collective variable spaces. Even though each individual method to estimate the rate passes its quality tests, the rate constant estimates exhibit disparities of up to four orders of magnitude. Rate constant estimates based on one-dimensional reaction coordinates prove challenging to converge, even if the reaction coordinate is optimized. However, consistent estimates of the rate constant are achieved by sampling transitions and by multi-dimensional grid-based models.}, language = {en} }