Toward a new approach for determination of solute's charge distribution to analyze interatomic electrostatic interactions in quantum mechanical/molecular mechanical simulations

We present an alternative approach to determine "density-dependent property"-derived charges for molecules in the condensed phase. In the case of a solution, it is essential to take into consideration the electron polarization of molecules in the active site of this system. The solute and solvent molecules in this site have to be described by a quantum mechanical technique and the others are allowed to be treated by a molecular mechanical method (QM/MM scheme). For calculations based on this scheme, using the forces and interaction energy as density-dependent property our charges from interaction energy and forces (CHIEF) approach can provide the atom-centered charges on the solute atoms. These charges reproduce well the electrostatic potentials around the solvent molecules and present properly the picture of the electron density of the QM subsystem in the solution system. Thus, the CHIEF charges can be considered as the atomic charges under the conditions of the QM/MM simulation, and then enable one to analyze electrostatic interactions between atoms in the QM and MM regions. This approach would give a view of the QM nuclei and electrons different from the conventional methods.     

Toward a consistent treatment of polarization in model QM/MM calculations

The concept of model chemistries within hybrid QM/MM calculations has been addressed through analysis of the polarization energy determined by two distinct approaches based on (i) induced charges and (ii) induced dipoles. The quantum mechanical polarization energy for four configurations of the water dimer has been determined for a range of basis sets using Morokuma energy decomposition analysis. This benchmark value has been compared to the fully classical polarization energy determined using the induced dipole approach, and the molecular mechanics polarization energy calculated using induced charges within the MM region of hybrid QM/MM calculations. From the water dimer calculations, it is concluded that the induced charge approach is consistent with medium sized basis set calculations whereas the induced dipole approach is consistent with large basis set calculations. This result is highly relevant to the concept of QM/MM model chemistries.     

A novel quantum mechanical/molecular mechanical approach to the free energy calculation for isomerization of glycine in aqueous solution

The free energy change associated with the isomerization reaction of glycine in water solution has been studied by a hybrid quantum mechanical/molecular mechanical (QM/MM) approach combined with the theory of energy representation (QM/MM-ER) recently developed. The solvation free energies for both neutral and zwitterionic form of glycine have been determined by means of the QM/MM-ER simulation. The contributions of the electronic polarization and the fluctuation of the QM solute to the solvation free energy have been investigated. It has been found that the contribution of the density fluctuation of the zwitterionic solute is estimated as -4.2 kcal/mol in the total solvation free energy of -46.1 kcal/mol, while that of the neutral form is computed as -3.0 kcal/mol in the solvation free energy of -15.6 kcal/mol. The resultant free energy change associated with the isomerization of glycine in water has been obtained as -7.8 kcal/mol, in excellent agreement with the experimental data of -7.3 or -7.7 kcal/mol, implying the accuracy of the QM/MM-ER approach. The results have also been compared with those computed by other methodologies such as the polarizable continuum model and the classical molecular simulation. The efficiency and advantage of the QM/MM-ER method has been discussed.     

Polarization and charge-transfer effects in aqueous solution via ab initio QM/MM simulations

Combined ab initio quantum mechanical and molecular mechanical (QM/MM) simulations coupled with the block-localized wave function energy decomposition (BLW-ED) method have been conducted to study the solvation of two prototypical ionic systems, acetate and methylammonium ions in aqueous solution. Calculations reveal that the electronic polarization between the targeted solutes and water is the primary many-body effect, whereas the charge-transfer term only makes a small fraction of the total solute-solvent interaction energy. In particular, the polarization effect is dominated by the solvent (water) polarization.     

Two-Photon Solvatochromism II: Experimental and Theoretical Study of Solvent Effects on the Two-Photon Absorption Spectrum of Reichardt’s Dye

In this study, we report on the influence of solvent on the two-photon absorption (2PA) spectra of Reichardt’s dye (RD). The measurement of 2PA cross-sections is performed for three solvents (chloroform, dimethyl formamide, and dimethyl sulfoxide) using the Z-scan technique. The key finding of this study is the observation that the cross-section, corresponding to the 2PA of the intramolecular charge-transfer state, diminishes substantially upon increasing the solvent polarity. To unravel the solvent dependence of the 2PA cross-section, the electronic structure of RD is determined using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach, in which polarization between the solute and solvent is taken into account by using a self-consistent scheme in the solvent polarization. The two-state approximation proves to be adequate for the studied system, and allowed the observed solvent-polarity-induced decrease of the 2PA cross-section to be related to the decrease of the transition moment and the increase in the excitation energy.     

An application of the novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: An ionic dissociation of a water molecule in the supercritical water

A novel quantum chemical approach recently developed has been applied to an ionic dissociation of a water molecule (2H(2)O-->H(3)O(+)+OH(-)) in ambient and supercritical water. The method is based on the quantum mechanical/molecular mechanical (QM/MM) simulations combined with the theory of energy representation (QM/MM-ER), where the energy distribution function of MM solvent molecules around a QM solute serves as a fundamental variable to determine the hydration free energy of the solute according to the rigorous framework of the theory of energy representation. The density dependence of the dissociation free energy in the supercritical water has been investigated for the density range from 0.1 to 0.6 g/cm(3) with the temperature fixed at a constant. It has been found that the product ionic species significantly stabilizes in the high density region as compared with the low density. Consequently, the dissociation free energy decreases monotonically as the density increases. The decomposition of the hydration free energy has revealed that the entropic term (-TDeltaS) strongly depends on the density of the solution and dominates the behavior of the dissociation free energy with respect to the variation of the density. The increase in the entropic term in the low density region can be attributed to the decrease in the translational degrees of freedom brought about by the aggregation of solvent water molecules around the ionic solute.     

Simple and exact approach to the electronic polarization effect on the solvation free energy: formulation for quantum-mechanical/molecular-mechanical system and its applications to aqueous solutions

We present a simple and exact numerical approach to compute the free energy contribution δμ in solvation due to the electron density polarization and fluctuation of a quantum-mechanical solute in the quantum-mechanical/molecular-mechanical (QM/MM) simulation combined with the theory of the energy representation (QM/MM-ER). Since the electron density fluctuation is responsible for the many-body QM-MM interactions, the standard version of the energy representation method cannot be applied directly. Instead of decomposing the QM-MM polarization energy into the pairwise additive and non-additive contributions, we take sum of the polarization energies in the QM-MM interaction and adopt it as a new energy coordinate for the method of energy representation. Then, it is demonstrated that the free energy δμ can be exactly formulated in terms of the energy distribution functions for the solution and reference systems with respect to this energy coordinate. The benchmark tests were performed to examine the numerical efficiency of the method with respect to the changes in the individual properties of the solvent and the solute. Explicitly, we computed the solvation free energy of a QM water molecule in ambient and supercritical water, and also the free-energy change associated with the isomerization reaction of glycine from neutral to zwitterionic structure in aqueous solution. In all the systems examined, it was demonstrated that the computed free energy δμ agrees with the experimental value, irrespective of the choice of the reference electron density of the QM solute. The present method was also applied to a prototype reaction of adenosine 5'-triphosphate hydrolysis where the effect of the electron density fluctuation is substantial due to the excess charge. It was demonstrated that the experimental free energy of the reaction has been accurately reproduced with the present approach.     

Predicting hydration free energies with a hybrid QM/MM approach: an evaluation of implicit and explicit solvation models in SAMPL4

The correct representation of solute-water interactions is essential for the accurate simulation of most biological phenomena. Several highly accurate quantum methods are available to deal with solvation by using both implicit and explicit solvents. So far, however, most evaluations of those methods were based on a single conformation, which neglects solute entropy. Here, we present the first test of a novel approach to determine hydration free energies that uses molecular mechanics (MM) to sample phase space and quantum mechanics (QM) to evaluate the potential energies. Free energies are determined by using re-weighting with the Non-Boltzmann Bennett (NBB) method. In this context, the method is referred to as QM-NBB. Based on snapshots from MM sampling and accounting for their correct Boltzmann weight, it is possible to obtain hydration free energies that incorporate the effect of solute entropy. We evaluate the performance of several QM implicit solvent models, as well as explicit solvent QM/MM for the blind subset of the SAMPL4 hydration free energy challenge. While classical free energy simulations with molecular dynamics give root mean square deviations (RMSD) of 2.8 and 2.3 kcal/mol, the hybrid approach yields an improved RMSD of 1.6 kcal/mol. By selecting an appropriate functional and basis set, the RMSD can be reduced to 1 kcal/mol for calculations based on a single conformation. Results for a selected set of challenging molecules imply that this RMSD can be further reduced by using NBB to reweight MM trajectories with the SMD implicit solvent model.     

Solvent effects on the nπ* transition of pyrimidine in aqueous solution

A hybrid quantum mechanical and molecular mechanical potential is used in Monte Carlo simulations to examine the solvent effects on the electronic excitation energy for the n→π* transition of pyrimidine in aqueous solution. In the present study, the pyrimidine molecule is described by the semi-empirical AM1 model, while the solvent molecules are treated classically. Two sets of calculations are performed: the first involves the use of the pairwise three-point charge TIP3P model for water, and the second computation employs a polarizable many-body potential for the solvent. The latter calculation takes into account the effect of solvent polarization following the solute electronic excitation, and makes a correction to the energies determined using pairwise potentials, which neglects such fast polarization effects and overestimates the solute-solvent interactions on the Franck-Condon excited states. Our simulation studies of pyrimidine in water indicate that the solvent charge redistribution following the solute electronic excitation makes modest corrections (about −130␣cm⁻¹) to the energy predicted by using pairwise potentials. Specific hydrogen bonding interactions between pyrimidine and water are important for the prediction of solvatochromic shifts for pyrimidine. The computed n→π* blue shift is 2275±110 cm⁻¹, which may be compared with the experimental value (2700 cm⁻¹) from isooctane to water. Received: 14 January 1997 / Accepted: 21 February 1997     

Novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: free energy calculation for the Beckmann rearrangement promoted by proton transfers in the supercritical water

The Beckmann rearrangement of acetone oxime promoted by proton transfers in the supercritical water has been investigated by means of the hybrid quantum mechanical/molecular mechanical approach combined with the theory of energy representation (QM/MM-ER) recently developed. The transition state (TS) structures have been explored by ab initio calculations for the reaction of hydrated acetone oxime on the assumption that the reaction is catalyzed by proton transfers along the hydrogen bonds connecting the solute and the solvent water molecules. Up to two water molecules have been considered as reactants that take part in the proton transfers. As a result of the density functional theory calculations with B3LYP functional and aug-cc-pVDZ basis set, it has been found that participation of two water molecules in the reaction reduces the activation free energy by -12.3 kcal/mol. Furthermore, the QM/MM-ER simulations have revealed that the TS is more stabilized than the reactant state in the supercritical water by 2.7 kcal/mol when two water molecules are involved in the reaction. Solvation free energies of the reactant and the TS have been decomposed into terms due to the electronic polarization of the solute, electron density fluctuation, and others to elucidate the origin of the stabilization of the TS as compared with the reactant. It has been revealed that the promotion of the chemical reaction due to the hydration mainly originates from the interaction between the nonpolarized solute and the solvent water molecules at the supercritical state.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

A noniterative mean-field QM/MM-type approach with a linear response approximation toward an efficient free-energy evaluation

Mean-field treatment of solvent provides an efficient technique to investigate chemical processes in solution in quantum mechanics/molecular mechanics (QM/MM) framework. In the algorithm, an iterative calculation is required to obtain the self-consistency between QM and MM regions, which is a time-consuming step. In the present study, we have proposed a noniterative approach by introducing a linear response approximation (LRA) into the solvation term in the one-electron part of Fock matrix in a hybrid approach between molecular-orbital calculations and a three-dimensional (3D) integral equation theory for molecular liquids (multicenter molecular Ornstein–Zernike self-consistent field [MC-MOZ-SCF]; Kido et al., J. Chem. Phys. 2015, 143, 014103). To save the computational time, we have also developed a fast method to generate electrostatic potential map near solute and the solvation term in Fock matrix, using Fourier transformation (FT) and real spherical harmonics expansion (RSHE). To numerically validate the LRA and FT-RSHE method, we applied the present approach to water, carbonic acid, and their ionic species in aqueous solution. Molecular properties of the solutes were evaluated by the present approach with four different types of initial wave functions and compared with those by the original (MC-MOZ-SCF). We found that an initial wave function considering solvation effects is needed to appropriately reproduce the properties by MC-MOZ-SCF. Furthermore, a benchmark test for 32 solute molecules was performed to evaluate the accuracy of the present approach for solvation free energy (SFE) and measure the speedup ratio for MC-MOZ-SCF. The error of SFE for MC-MOZ-SCF does not correlate with the SFE but increases in proportion to the electronic reorganization energy. Similar to water and carbonic acid, an initial wave function with solvation effects is also important to make the error small. From the averaged speed up ratio, the present approach is 13.5 times faster than MC-MOZ-SCF. © 2019 Wiley Periodicals, Inc.     

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

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.