Electrostatic embedding in large-scale first principles quantum mechanical calculations on biomolecules.

Area of Science:

• Computational Chemistry
• Biomolecular Modeling
• Quantum Mechanics

Background:

• Atomistic biomolecular simulations require high accuracy for binding free energies and enzymatic reactions.
• Current force fields face limitations in accuracy and transferability, especially for unusual ligands.
• First-principles methods like density functional theory (DFT) offer accuracy but are computationally expensive.

Purpose of the Study:

• To develop and evaluate an electrostatic embedding scheme for accurate and efficient biomolecular simulations.
• To enable the treatment of entire biomolecules and large solvent portions using DFT.
• To assess the accuracy of this method for calculating ligand-biomolecule and ligand-solvent interaction energies.

Main Methods:

• Implementation of a self-consistent electrostatic embedding scheme within a linear-scaling DFT program.
• Representation of the environment using localized charge distributions electrostatically coupled to the quantum system.
• Testing the scheme on interactions of neutral and charged ligands with biomolecules and solvent.

Main Results:

• The electrostatic embedding scheme allows for accurate calculations of interaction energies.
• Achieved accuracy comparable to full DFT calculations for the entire system.
• Demonstrated applicability to a variety of neutral and charged species.

Conclusions:

• Electrostatic embedding offers a computationally efficient approach to high-accuracy biomolecular simulations.
• This method overcomes limitations of traditional force fields for complex systems.
• Enables accurate DFT-level treatment of large biomolecular systems and their environments.