Microscopic models for quantum mechanical calculations of chemical processes in solutions: LD/AMPAC and SCAAS/AMPAC calculations of solvation energies

Abstract

Our previously developed approaches for integrating quantum mechanical molecular orbital methods with microscopic solvent models are refined and examined. These approaches consider the nonlinear solute–solvent coupling in a self‐consistent way by incorporating the potential from the solvent dipoles in the solute Hamiltonian, while considering the polarization of the solvent by the potential from the solute charges. The solvent models used include the simplified Langevin Dipoles (LD) model and the much more expensive surface constrained All Atom Solvent (SCAAS) model, which is combined with a free energy pertubation (FEP) approach. Both methods are effectively integrated with the quantum mechanical AMPAC package and can be easily combined with other quantum mechanical programs. The advantages of the present approaches and their earlier versions over macroscopic reaction field models and supermolecular approaches are considered. A LD/MNDO study of solvated organic ions demonstrates that this model can yield reliable solvation energies, provided the quantum mechanical charges are scaled to have similar magnitudes to those obtained by high level ab initio methods. The incorporation of a field‐dependent hydrophobic term in the LD free energy makes the present approach capable of evaluating the free energy of transfer of polar molecules from non polar solvents to aqueous solutions. The reliability of the LD approach is examined not only by evaluating a rather standard set of solvation energies of organic ions and polar molecules, but also by considering the stringent test case of sterically hindered hydrophobic ions. In this case, we compare the LD/MNDO solvation energies to the more rigorous FEP/SCAAS/MNDO solvation energies. Both methods are found to give similar results even in this challenging test case. The FEP/SCAAS/AMPAC method is incorporated into the current version of the program ENZYMIX. This option allows one to study chemical reactions in enzymes and in solutions using the MNDO and AM1 approximations. A special procedure that uses the EVB method as a reference potential for SCF MO calculations should help in improving the reliability of such studies.