MDRefine: a Python package for refining Molecular Dynamics trajectories with experimental data

Abstract

Molecular dynamics (MD) simulations play a crucial role in resolving the underlying conformational dynamics of molecular systems. However, their capability to correctly reproduce and predict dynamics in agreement with experiments is limited by the accuracy of the force-field model. This capability can be improved by refining the structural ensembles or the force-field parameters. Furthermore, discrepancies with experimental data can be due to imprecise forward models, namely, functions mapping simulated structures to experimental observables. Here, we introduce MDRefine, a Python package aimed at implementing the refinement of the ensemble, the force-field and/or the forward model by comparing MD-generated trajectories with experimental data. The software consists of several tools that can be employed separately from each other or combined together in different ways, providing a seamless interpolation between these three different types of refinement. We use some benchmark cases to show that the combined approach is superior to separately applied refinements. Source code, documentation and examples are freely available at https://pypi.org/project/MDRefine and https://github.com/bussilab/MDRefine.

39 pages, 20 figures

References in corpus (16)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Building Water Models, A Different Approach
• Bayesian ensemble refinement by replica simulations and reweighting
• Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments
• Combining simulations and solution experiments as a paradigm for RNA force field refinement
• Shape Transitions and Chiral Symmetry Breaking in the Energy Landscape of the Mitotic Chromosome
• Towards empirical force fields that match experimental observables
• Empirical corrections to the Amber RNA force field with Target Metadynamics
• Optimizing Millions of Hyperparameters by Implicit Differentiation
• Integrating experimental data with molecular simulations to investigate RNA structural dynamics
• chemtrain: Learning Deep Potential Models via Automatic Differentiation and Statistical Physics
• Refinement of molecular dynamics ensembles using experimental data and flexible forward models
• Molecular simulations matching denaturation experiments for N6-Methyladenosine
• Boosting ensemble refinement with transferable force field corrections: synergistic optimization for molecular simulations
• Molecular simulations to investigate the impact of N6-methylation in RNA recognition: Improving accuracy and precision of binding free energy prediction