Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems

Abstract

In modern day simulations of many-body systems much of the computational complexity is shifted to the identification of slowly changing molecular order parameters called collective variables (CV) or reaction coordinates. A vast array of enhanced sampling methods are based on the identification and biasing of these low-dimensional order parameters, whose fluctuations are important in driving rare events of interest. Here describe a new algorithm for finding optimal low-dimensional collective variables for use in enhanced sampling biasing methods like umbrella sampling, metadynamics and related methods, when limited prior static and dynamic information is known about the system, and a much larger set of candidate CVs is specified. The algorithm involves estimating the best combination of these candidate CVs, as quantified by a maximum path entropy estimate of the spectral gap for dynamics viewed as a function of that CV. Through multiple practical examples, we show how this post-processing procedure can lead to optimization of CV and several orders of magnitude improvement in the convergence of the free energy calculated through metadynamics, essentially giving the ability to extract useful information even from unsuccessful metadynamics runs.

7 pages, 4 figures; corrected missing figure number and added a reference

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• Metadynamics with adaptive Gaussians
• The role of water and steric constraints in the kinetics of cavity-ligand unbinding
• Inferring transition rates on networks with incomplete knowledge

Cited by in corpus (55)

• Data-Driven Collective Variables for Enhanced Sampling
• Unsupervised machine learning in atomistic simulations, between predictions and understanding
• A variational conformational dynamics approach to the selection of collective variables in metadynamics
• State Predictive Information Bottleneck
• Enhancing entropy and enthalpy fluctuations to drive crystallization in atomistic simulations
• A unified framework for machine learning collective variables for enhanced sampling simulations: $\text{mlcolvar}$
• Maximum Caliber: a general variational principle for dynamical systems
• Special Topic: Markov Models of Molecular Kinetics
• From enhanced sampling to reaction profiles
• Recognizing Local and Global Structural Motifs at the Atomic Scale
• Deep Learning Collective Variables from Transition Path Ensemble
• Gaussian Mixture Based Enhanced Sampling
• Unfolding Hidden Barriers by Active Enhanced Sampling
• Finding Efficient Collective Variables: The Case of Crystallization
• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Predicting reaction coordinates in energy landscapes with diffusion anisotropy
• Multiscale reweighted stochastic embedding (MRSE): Deep learning of collective variables for enhanced sampling
• Making the best of a bad situation: a multiscale approach to free energy calculation
• Identification and Analysis of Transition and Metastable Markov States
• Reaction Coordinates for Conformational Transitions using Linear Discriminant Analysis on Positions
• How wet should be the reaction coordinate for ligand unbinding?
• Enhanced sampling of transition states
• Reaction coordinates and rate constants for liquid droplet nucleation: quantifying the interplay between driving force and memory
• A Deep Autoencoder Framework for Discovery of Metastable Ensembles in Biomacromolecules
• Adding alchemical variables to metadynamics to enhance sampling in free energy calculations
• Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence
• Nucleating a Different Coordination in a Crystal under Pressure: A Study of the B1-B2 Transition in NaCl by Metadynamics
• Collective Variables for Conformational Polymorphism in Molecular Crystals
• Assessing models of force-dependent unbinding rates via infrequent metadynamics
• Spectral Map: Embedding Slow Kinetics in Collective Variables
• Reweighted Manifold Learning of Collective Variables from Enhanced Sampling Simulations
• Multidimensional minimum-work control of a 2D Ising model
• Kramers turnover: from energy diffusion to spatial diffusion using metadynamics
• Kinetics and Free Energy of Ligand Dissociation Using Weighted Ensemble Milestoning
• An implementation of the maximum-caliber principle by replica-averaged time-resolved restrained simulations
• Manifold Learning in Atomistic Simulations: A Conceptual Review
• Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?
• Optimal reaction coordinates and kinetic rates from the projected dynamics of transition paths
• Adaptive Ensemble Biomolecular Simulations at Scale
• Free Energy Landscapes, Diffusion Coefficients, and Kinetic Rates from Transition Paths
• Deep learning the slow modes for rare events sampling
• maze: Heterogeneous Ligand Unbinding along Transient Protein Tunnels
• Selecting High-Dimensional Representations of Physical Systems by Reweighted Diffusion Maps
• On the Role of Solvent in Hydrophobic Cavity-ligand Recognition Kinetics
• Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles
• Thermodynamically Optimized Machine-learned Reaction Coordinates for Hydrophobic Ligand Dissociation
• Extending and validating bubble nucleation rate predictions in a Lennard-Jones fluid with enhanced sampling methods and transition state theory
• Learning Markovian Dynamics with Spectral Maps
• Identifying the leading dynamics of ubiquitin: a comparison between the tICA and the LE4PD slow fluctuations in amino acids' position
• A Sinking Approach to Explore Arbitrary Areas in Free Energy Landscapes
• Machine learning approaches for analyzing and enhancing molecular dynamics simulations
• Revealing the Atomistic Mechanism of Rare Events in Molecular Dynamics
• Using data-reduction techniques to analyse biomolecular trajectories
• Thermal Fluctuations of Singular Bar-Joint Mechanisms
• Computing rates and understanding unbinding mechanism in host-guest systems