Automated Parametrization of Biomolecular Force Fields from Quantum Mechanics/Molecular Mechanics (QM/MM) Simulations through Force Matching
Citations Over TimeTop 11% of 2007 papers
Abstract
We introduce a novel procedure to parametrize biomolecular force fields. We perform finite-temperature quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, with the fragment or moiety that has to be parametrized being included in the QM region. By applying a force-matching algorithm, we derive a force field designed in order to reproduce the steric, electrostatic, and dynamic properties of the QM subsystem. The force field determined in this manner has an accuracy that is comparable to the one of the reference QM/MM calculation, but at a greatly reduced computational cost. This allows calculating quantities that would be prohibitive within a QM/MM approach, such as thermodynamic averages involving slow motions of a protein. The method is tested on three different systems in aqueous solution: dihydrogenphosphate, glycyl-alanine dipeptide, and a nitrosyl-dicarbonyl complex of technetium(I). Molecular dynamics simulations with the optimized force field show overall excellent performance in reproducing properties such as structures and dipole moments of the solutes as well as their solvation pattern.
Related Papers
- → Automated Parametrization of Biomolecular Force Fields from Quantum Mechanics/Molecular Mechanics (QM/MM) Simulations through Force Matching(2007)85 cited
- → Molecular Dynamics Simulations of Proteins: Can the Explicit Water Model Be Varied?(2007)71 cited
- → NMR Solvent Shifts of Adenine in Aqueous Solution from Hybrid QM/MM Molecular Dynamics Simulations(2007)35 cited
- → Benchmarking Hybrid Atomistic/Coarse-Grained Schemes for Proteins with an Atomistic Water Layer(2019)6 cited
- → Use of Crystal MD Simulations to Speed Up Evaluation of Binding Free Energies of Dimannose Deoxy Analogs With M4-P51g-Cyanovirin-N(2010)