Polarizable Force Field for Molecular Ions Based on the Classical Drude Oscillator

Area of Science:

• Computational chemistry
• Molecular modeling
• Biophysics

Background:

• Accurate polarizable models are essential for simulating biological macromolecules.
• Developing force field parameters for molecular ions is a critical step.
• Existing models often lack sufficient accuracy for ionic moieties.

Purpose of the Study:

• To develop and optimize accurate force field parameters for molecular ions within a polarizable energy function.
• To compare these parameters with the nonpolarizable CHARMM general force field (CGenFF).
• To establish a foundation for future polarizable force fields for biological systems.

Main Methods:

• Hierarchical approach: parameter optimization for small molecules serving as building blocks.
• Model compounds: cationic (ammonium derivatives, imidazolium, guanidinium) and anionic (acetate, phenolate, alkanethiolates).
• Development within the SWM4-NDP polarizable water model, ensuring hydration free energies match atomic ion data.

Main Results:

• Optimized parameters show good agreement with quantum mechanical (QM) and experimental data.
• Significant differences in water structure around ions observed between Drude and additive force fields.
• Highlights the crucial role of polarization in aqueous solvation.

Conclusions:

• The developed parameters are foundational for future Drude polarizable force fields.
• These parameters will enhance simulations of biological macromolecules and drug-like molecules.
• Polarization effects are vital for accurately modeling molecular solvation.