Combining simulations and solution experiments as a paradigm for RNA force field refinement
arXiv:1611.06156 · doi:10.1021/acs.jctc.6b00944
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 (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
- Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies
- Free Energy Landscape of GAGA and UUCG RNA Tetraloops
- RNA unwinding from reweighted pulling simulations