Bridging the 12-6-4 Model and the Fluctuating Charge Model

Area of Science:

• Computational chemistry and biophysics
• Molecular modeling and simulation
• Biomolecular simulations

Background:

• Classical force fields are crucial for atomic-level molecular dynamics (MD) simulations of biological systems.
• The standard 12-6 Lennard-Jones model fails to accurately predict hydration free energy and ion-oxygen distances for highly charged metal ions (+2 or higher).
• Previous 12-6-4 models addressed ion-induced dipole interactions but overestimated coordination numbers for certain ions.

Purpose of the Study:

• To develop a refined molecular dynamics model that accurately simulates metal ion behavior in biological systems.
• To address the limitations of existing models, specifically the overestimation of coordination numbers for highly charged metal ions.
• To improve the physical realism of MD simulations for metal ions, particularly in the context of metalloproteins.

Main Methods:

• A modified 12-6 model was employed, adjusting atomic charges of first solvation shell water molecules to incorporate ion-induced dipole interactions.
• This approach accounts for both ion-dipole interactions and increased repulsion among water molecules in the first solvation shell.
• Simulations were validated against experimental hydration free energy and ion-oxygen distance data for various metal ions (Mg2+, Zn2+, Al3+, Fe3+, In3+).

Main Results:

• The new model successfully reproduced experimental hydration free energy and ion-oxygen distance values for Mg2+, Zn2+, Al3+, Fe3+, and In3+.
• The overestimation of coordination numbers observed with the 12-6-4 model for Fe3+ and In3+ was resolved.
• Simulation results demonstrated strong agreement with ab initio MD simulations and revealed a physically meaningful relationship between the C4 parameter and induced dipole moment.

Conclusions:

• The developed model provides a more accurate and physically grounded approach for simulating metal ions in biological systems using molecular dynamics.
• This work bridges the gap between the 12-6-4 model and fluctuating charge models, offering improved simulation capabilities.
• The findings have significant implications for enhancing the accuracy of molecular dynamics simulations of metalloproteins.