Combining simulations and solution experiments as a paradigm for RNA force field refinement

Abstract

Recent computational efforts have shown that the current potential energy models used in molecular dynamics are not accurate enough to describe the conformational ensemble of RNA oligomers and suggest that molecular dynamics should be complemented with experimental data. We here propose a scheme based on the maximum entropy principle to combine simulations with bulk experiments. In the proposed scheme the noise arising from both the measurements and the forward models used to back calculate the experimental observables is explicitly taken into account. The method is tested on RNA nucleosides and is then used to construct chemically consistent corrections to the Amber RNA force field that allow a large set of experimental data on nucleosides and dinucleosides to be correctly reproduced. The transferability of these corrections is assessed against independent data on tetranucleotides and displays a previously unreported agreement with experiments. This procedure can be applied to enforce multiple experimental data on multiple systems in a self-consistent framework thus suggesting a new paradigm for force field refinement.

Reprinted (adapted) with permission from J Chem Theory Comput doi: 10.1021/acs.jctc.6b00944. Copyright (2016) American Chemical Society

References in corpus (15)

• Canonical sampling through velocity-rescaling
• PLUMED 2: New feathers for an old bird
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• Bayesian ensemble refinement by replica simulations and reweighting
• Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies
• Enhanced Conformational Sampling using Replica Exchange with Collective-Variable Tempering
• Free Energy Landscape of GAGA and UUCG RNA Tetraloops
• Empirical corrections to the Amber RNA force field with Target Metadynamics
• Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments
• Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch
• RNA Folding Pathways in Stop Motion
• The inefficiency of re-weighted sampling and the curse of system size in high order path integration
• RNA unwinding from reweighted pulling simulations
• Accurate multiple time step in biased molecular simulations

Cited by in corpus (19)

• Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems
• Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments
• Analyzing and biasing simulations with PLUMED
• Exploring RNA structure and dynamics through enhanced sampling simulations
• Towards empirical force fields that match experimental observables
• Girsanov reweighting for path ensembles and Markov state models
• Predicting the kinetics of RNA oligonucleotides using Markov state models
• Assessing the accuracy of direct-coupling analysis for RNA contact prediction
• Recent Advances in Maximum Entropy Biasing Techniques for Molecular Dynamics
• Estimation of the infinitesimal generator by square-root approximation
• 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
• Inferring Spatial Source of Disease Outbreaks using Maximum Entropy
• Frontiers in integrative structural biology: modeling disordered proteins and utilizing in situ data
• Simulation-Based Inference with Approximately Correct Parameters via Maximum Entropy
• Boosting ensemble refinement with transferable force field corrections: synergistic optimization for molecular simulations
• MDRefine: a Python package for refining Molecular Dynamics trajectories with experimental data
• A Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water
• Automatic learning of hydrogen-bond fixes in an AMBER RNA force field