Abstract:Quantifying the conformational ensembles of biomolecules is fundamental to describing mechanisms of processes such as ligand binding and allosteric regulation. Accurate quantification of these ensembles remains a challenge for all but the simplest molecules. One such challenge is insufficient sampling which enhanced sampling approaches, such as metadynamics, were designed to overcome; however, the non-uniform frame weights that result from many of these approaches present an additional challenge to ensemble quantification techniques such as Markov State Modeling or structural clustering. Here, we present rigorous inclusion of non-uniform frame weights into a structural clustering method entitled shapeGMM. The shapeGMM method fits a Gaussian mixture model to particle positions, and here we advance that approach by incorporating nonuniform frame weights in the estimates of all parameters of the model. The resulting models are high dimensional probability densities for the unbiased systems from which we can compute important thermodynamic properties such as relative free energies and configurational entropy. The accuracy of this approach is demonstrated by the quantitative agreement between GMMs computed by Hamiltonian reweighting and direct simulation of a coarse-grained helix model system. Furthermore, the relative free energy computed from a high dimensional probability density of alanine dipeptide reweighted from a metadynamics simulation quantitatively reproduces the metadynamics free energy in the basins. Finally, the method identifies hidden structures along the actin globular to filamentous-like structural transition from a metadynamics simulation on a linear discriminant analysis coordinate trained on GMM states, demonstrating the broad applicability of combining our prior and new methods, and illustrating how structural clustering of biased data can lead to biophysical insight. Combined, these results demonstrate that frame-weighted shapeGMM is a powerful approach to quantify biomolecular ensembles from biased simulations.

Lattice Models for Magnetic Fluids: Correlations Between Order Parameters

On the Stochastic Dynamics of Molecular Conformation

Multiscale macromolecular simulation: role of evolving ensembles

Equations for Stochastic Macromolecular Mechanics of Single Proteins: Equilibrium Fluctuations, Transient Kinetics and Nonequilibrium Steady-State

Structuring and sampling complex conformation space: Weighted ensemble dynamics simulations

Ensemble annealing of complex physical systems

Extended Ensemble Monte Carlo

Quantifying Unbiased Conformational Ensembles from Biased Simulations Using ShapeGMM

Stochastic simulations for the time evolution of systems which obey generalized statistics: Fractional exclusion statistics and Gentile's statistics

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Global Langevin model of multidimensional biomolecular dynamics

Molecular dynamics in a grand ensemble: Bergmann–Lebowitz model and adaptive resolution simulation

Simulating rare events using a Weighted Ensemble-based string method

A molecular dynamics method for simulations in the canonical ensemble

Integrating Molecular Simulation and Experimental Data: A Bayesian/Maximum Entropy reweighting approach

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions

A stochastic algorithm for the dynamic evolution process by particle deposition

Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation

Ensembles Are Required to Handle Aleatoric and Parametric Uncertainty in Molecular Dynamics Simulation

Stochastic Optimization Based Study of Dimerization Kinetics

Atomically Detailed Simulations of Helix Formation with the Stochastic Difference Equation