Accurately reproducing ab initio electrostatic potentials with multipoles and fragmentation.

Area of Science:

• Computational chemistry
• Molecular modeling
• Biophysics

Background:

• Accurate calculation of molecular electrostatic potentials (ESPs) is crucial for understanding molecular interactions.
• Traditional methods like ab initio and density functional theory (DFT) are computationally expensive for large systems.
• Point charge models offer a computationally cheaper alternative but often lack accuracy.

Purpose of the Study:

• To develop and validate an energy-based fragmentation method for calculating ESPs.
• To compare the accuracy of this fragmentation method against full ab initio/DFT calculations.
• To assess the performance of various point charge models in reproducing ESPs.

Main Methods:

• An energy-based fragmentation approach (Bettens & Lee, 2006) was employed.
• Calculations were performed on peptides and other biologically relevant molecules.
• Comparison with full ab initio and DFT ESPs was conducted.
• Various point charge models were evaluated.

Main Results:

• The fragmentation method accurately reproduced ESPs for charged and uncharged peptides and other molecules.
• The fragmentation method significantly outperformed standard point charge models.
• The approach was successfully applied to a large system (neuraminidase tetramer, ~24,000 atoms).
• Using distributed monopoles, dipoles, and quadrupoles yielded ESPs comparable to the fragmentation method.

Conclusions:

• The energy-based fragmentation method provides an accurate and efficient way to compute ESPs for large biomolecules.
• This method offers a significant improvement over traditional point charge models.
• The approach is scalable and suitable for complex biological systems.