Property Map Collective Variable as a Useful Tool for a Force Field Correction.

Area of Science:

• Computational Chemistry
• Molecular Dynamics Simulations
• Biomolecular Modeling

Background:

• The accuracy of biomolecular simulations critically relies on the precision of molecular mechanics force fields.
• Parametrizing force fields for diverse small organic molecules presents a significant challenge.
• Existing force fields often require corrections to accurately represent molecular behavior.

Purpose of the Study:

• To develop an efficient method for estimating density functional theory (DFT)-derived force field corrections.
• To enable accurate biomolecular simulations for drug-like molecules by improving force field parameters.
• To integrate DFT-based corrections into molecular dynamics (MD) simulations without substantial computational overhead.

Main Methods:

• Proposed an approach to estimate DFT-derived force field corrections using the property map collective variable (PMCV) formula.
• Approximated the force field correction as a weighted average of corrections calculated for a limited set of reference structures.
• Validated the method by converting between seven AMBER force fields and applying it to the anticancer drug Imatinib.

Main Results:

• Demonstrated the ability to effectively convert one AMBER force field to mimic the behavior of another.
• Successfully calculated force field corrections for the anticancer drug Imatinib, showcasing practical applicability.
• The proposed method allows for force field adjustment for general drug-like molecules with minimal impact on simulation speed.

Conclusions:

• The developed method provides an efficient way to estimate DFT-derived force field corrections for molecular dynamics.
• This approach is suitable for refining force fields of drug-like molecules, enhancing simulation accuracy.
• A computational pipeline for generating these corrections is publicly available.