Bayesian ensemble refinement by replica simulations and reweighting

Abstract

We describe different Bayesian ensemble refinement methods, examine their interrelation, and discuss their practical application. With ensemble refinement, the properties of dynamic and partially disordered (bio)molecular structures can be characterized by integrating a wide range of experimental data, including measurements of ensemble-averaged observables. We start from a Bayesian formulation in which the posterior is a functional that ranks different configuration space distributions. By maximizing this posterior, we derive an optimal Bayesian ensemble distribution. For discrete configurations, this optimal distribution is identical to that obtained by the maximum entropy "ensemble refinement of SAXS" (EROS) formulation. Bayesian replica ensemble refinement enhances the sampling of relevant configurations by imposing restraints on averages of observables in coupled replica molecular dynamics simulations. We show that the strength of the restraint should scale linearly with the number of replicas to ensure convergence to the optimal Bayesian result in the limit of infinitely many replicas. In the "Bayesian inference of ensembles" (BioEn) method, we combine the replica and EROS approaches to accelerate the convergence. An adaptive algorithm can be used to sample directly from the optimal ensemble, without replicas. We discuss the incorporation of single-molecule measurements and dynamic observables such as relaxation parameters. The theoretical analysis of different Bayesian ensemble refinement approaches provides a basis for practical applications and a starting point for further investigations.

Paper submitted to The Journal of Chemical Physics (15 pages, 4 figures); updated references; expanded discussions of related formalisms, error treatment, and ensemble refinement with and without replicas; appendix

References in corpus (1)

• Beyond Boltzmann-Gibbs statistics: Maximum entropy hyperensembles out-of-equilibrium

Cited by in corpus (17)

• Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments
• Combining simulations and solution experiments as a paradigm for RNA force field refinement
• Towards empirical force fields that match experimental observables
• Learning to Evolve Structural Ensembles of Unfolded and Disordered Proteins Using Experimental Solution Data
• Recent Advances in Maximum Entropy Biasing Techniques for Molecular Dynamics
• An implementation of the maximum-caliber principle by replica-averaged time-resolved restrained simulations
• Refinement of molecular dynamics ensembles using experimental data and flexible forward models
• Molecular simulations minimally restrained by experimental data
• Frontiers in integrative structural biology: modeling disordered proteins and utilizing in situ data
• Boosting ensemble refinement with transferable force field corrections: synergistic optimization for molecular simulations
• Simulation-Based Inference with Approximately Correct Parameters via Maximum Entropy
• MDRefine: a Python package for refining Molecular Dynamics trajectories with experimental data
• Progress in deep Markov State Modeling: Coarse graining and experimental data restraints
• Optimal weights and priors in simultaneous fitting of multiple small-angle scattering datasets
• Automatic learning of hydrogen-bond fixes in an AMBER RNA force field
• Computational Methods to Investigate Intrinsically Disordered Proteins and their Complexes
• Determination of protein structural ensembles using cryo-electron microscopy