Abstract:Interest in atomically detailed simulations has grown significantly with recent advances in computational hardware and Markov state modeling (MSM) methods, yet outstanding questions remain that hinder their widespread adoption. Namely, how do alternative sampling strategies explore conformational space and how might this influence predictions generated from the data? Here, we seek to answer these questions for four commonly used sampling methods: (1) a single long simulation, (2) many short simulations run in parallel, (3) adaptive sampling, and (4) our recently developed goal-oriented sampling algorithm, FAST. We first develop a theoretical framework for analytically calculating the probability of discovering select states on simple landscapes, where we uncover the drastic effects of varying the number and length of simulations. We then use kinetic Monte Carlo simulations on a variety of physically inspired landscapes to characterize the probability of discovering particular states and transition pathways for each of the four methods. Consistently, we find that FAST simulations discover each target state with the highest probability, while traversing realistic pathways. Furthermore, we uncover the potential pathology that short parallel simulations sometimes predict an incorrect transition pathway by crossing large energy barriers that long simulations would typically circumnavigate. We refer to this pathology as "pathway tunneling". To protect against this phenomenon when using adaptive-sampling and FAST simulations, we introduce the FAST-string method. This method enhances sampling along the highest-flux transition paths to refine an MSMs transition probabilities and discriminate between competing pathways. Additionally, we compare the performance of a variety of MSM estimators in describing accurate thermodynamics and kinetics. For adaptive sampling, we recommend simply normalizing the transition counts out of each state after adding small pseudocounts to avoid creating sources or sinks. Lastly, we evaluate whether our insights from simple landscapes hold for all atom molecular dynamics simulations of the folding of the lambda-repressor protein. Remarkably, we find that FAST-contacts predicts the same folding pathway as a set of long simulations but with orders of magnitude less simulation time.

Molecular search with conformational change: One-dimensional discrete-state stochastic model

On the Stochastic Dynamics of Molecular Conformation

Learning Effective Molecular Models from Experimental Observables.

How conformational dynamics influences the protein search for targets on DNA

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

Large-Scale Biomolecular Conformational Transitions Explored by a Combined Elastic Network Model and Enhanced Sampling Molecular Dynamics.

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Computer Chemical Simulation— A New Method of Molecular Conformational Analysis

A Combined Systematic-Stochastic Algorithm for the Conformational Search in Flexible Acyclic Molecules

A Conformational Search Method for Protein Systems Using Genetic Crossover and Metropolis Criterion

Molecular Conformation Generation via Shifting Scores

Choice of adaptive sampling strategy impacts state discovery, transition probabilities, and the apparent mechanism of conformational changes.

First-principles molecular structure search with a genetic algorithm

Elucidating molecular mechanisms of functional conformational changes of proteins via Markov state models

A Variable Neighbourhood Descent Heuristic for Conformational Search Using a Quantum Annealer

Investigation of Markovian and Non-Markovian Search Processes of Monomers of a Rouse Chain Confined in a Spherical Cavity

Exploring Conformational Landscapes Along Anharmonic Low-Frequency Vibrations

Multiple Simulated Annealing-Molecular Dynamics (MSA-MD) for Conformational Space Search of Peptide and Miniprotein

Markov Chain Monte Carlo in the Analysis of Single-Molecule Experimental Data

Exploring conformational search protocols for ligand-based virtual screening and 3-D QSAR modeling