A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions.

Area of Science:

• Computational Chemistry
• Theoretical Chemistry
• Molecular Modeling

Background:

• The three-dimensional reference interaction site model integral equation theory (3D-RISM) combined with ab initio molecular orbital (MO) theory (3D-RISM-SCF) is a powerful tool for studying solvation structures.
• Calculating the electrostatic potential in 3D-RISM-SCF is computationally intensive, limiting its application to large systems.
• Accurate solvation structure determination is crucial for understanding macromolecular behavior in solution.

Purpose of the Study:

• To develop a computationally efficient procedure for calculating electrostatic potential and solvated Fock matrix in 3D-RISM-SCF.
• To reduce the computational cost associated with the ab initio molecular orbital calculations in macromolecular solvation studies.
• To enable more feasible application of 3D-RISM-SCF to complex solvated macromolecular systems.

Main Methods:

• A novel procedure was developed to optimize the calculation of electrostatic potential and solvated Fock matrix based on solute-solvent distance.
• The method avoids calculations within atomic repulsive cores and uses classical approximations in distant regions.
• Direct integration of solute molecular orbitals is employed for the intermediate region.

Main Results:

• The proposed procedure significantly reduces computational cost for electrostatic potential and solvated Fock matrix calculations.
• Electronic structure, energy gradients of Methionine-Enkephalin, and solvation structure were accurately estimated in aqueous solution.
• Results obtained using the new procedure show good agreement with traditional methods and continuum models.

Conclusions:

• The developed procedure offers a substantial improvement in computational efficiency for 3D-RISM-SCF calculations of solvated macromolecules.
• This advancement facilitates more accurate and feasible studies of solvation effects on macromolecular systems.
• The method provides a valuable alternative for researchers investigating molecular solvation and its impact on chemical and biological processes.