Atomic Charge Parameters for the Finite Difference Poisson-Boltzmann Method Using Electronegativity Neutralization.

Area of Science:

• Computational Chemistry
• Physical Chemistry

Background:

• Accurate calculation of solvation free energies is crucial for understanding molecular interactions.
• Existing methods like charge equilibration (QEq) require optimization for specific computational approaches.

Purpose of the Study:

• To optimize the Rappe and Goddard charge equilibration (QEq) method for efficient and accurate solvation free energy calculations.
• To validate the optimized QEq method using experimental data and compare it with existing models.

Main Methods:

• Optimization of atomic electronegativity values within the QEq framework.
• Utilizing the finite difference Poisson-Boltzmann (FDPB) method for solvation energy calculations.
• Testing against experimental solvation free energies for small molecules.

Main Results:

• Achieved excellent agreement with experimental solvation free energies (mean unsigned error of ~0.5 kcal/mol).
• Demonstrated the QEq method's ability to capture conformational and solvent induction effects on charge distribution.
• Compared QEq results with fixed CHARMm charges for alanine dipeptide in water.

Conclusions:

• The optimized QEq method provides a fast and accurate approach for calculating solvation free energies.
• This method enhances molecular dynamics studies by better reflecting charge distribution changes.
• The optimized QEq method is a valuable tool for computational chemistry research.