Investigation of the dominant hydration structures among the ionic species in aqueous solution: Novel quantum mechanics/molecular mechanics simulations combined with the theory of energy representation

Omi

Morita

et al. 2012

Self Cite

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.

Section: A Qm/mm-er Approachmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

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

Omi

Morita

et al. 2012

Self Cite

“…(9) were examined by various applications both in the QM/MM and classical simulations. [31][32][33][34][35][36][37] We refer the readers to the previous works 19 for theoretical details about formulation of the free energy μ.…”

Section: A the Qm/mm-er Methodsmentioning

confidence: 99%

Computation of the free energy due to electron density fluctuation of a solute in solution: A QM/MM method with perturbation approach combined with a theory of solutions

Suzuoka

Morita

2014

Self Cite

We developed a perturbation approach to compute solvation free energy Δμ within the framework of QM (quantum mechanical)/MM (molecular mechanical) method combined with a theory of energy representation (QM/MM-ER). The energy shift η of the whole system due to the electronic polarization of the solute is evaluated using the second-order perturbation theory (PT2), where the electric field formed by surrounding solvent molecules is treated as the perturbation to the electronic Hamiltonian of the isolated solute. The point of our approach is that the energy shift η, thus obtained, is to be adopted for a novel energy coordinate of the distribution functions which serve as fundamental variables in the free energy functional developed in our previous work. The most time-consuming part in the QM/MM-ER simulation can be, thus, avoided without serious loss of accuracy. For our benchmark set of molecules, it is demonstrated that the PT2 approach coupled with QM/MM-ER gives hydration free energies in excellent agreements with those given by the conventional method utilizing the Kohn-Sham SCF procedure except for a few molecules in the benchmark set. A variant of the approach is also proposed to deal with such difficulties associated with the problematic systems. The present approach is also advantageous to parallel implementations. We examined the parallel efficiency of our PT2 code on multi-core processors and found that the speedup increases almost linearly with respect to the number of cores. Thus, it was demonstrated that QM/MM-ER coupled with PT2 deserves practical applications to systems of interest.

“…In the conventional theory of solutions, 10 the spatial distribution of the solvent molecules around the solute serves as a fundamental variable to evaluate the solvation free energy. In previous papers, [19][20][21][22][23][24][25] we performed QM/ MM-ER simulations for various chemical events in aqueous solutions to test the accuracy and efficiency of the method and found that our approach is quite adequate to compute free energy change for reactions in solution. 17,18 The advantage to combine the theory of energy representation with the QM/MM method is obvious.…”

Section: Introductionmentioning

confidence: 97%

Computation of the free energy change associated with one-electron reduction of coenzyme immersed in water: A novel approach within the framework of the quantum mechanical/molecular mechanical method combined with the theory of energy representation

Ohno

Kishi

et al. 2008

Self Cite

The isoalloxazine ring (flavin ring) is a part of the coenzyme flavin adenine dinucleotide and acts as an active site in the oxidation of a substrate. We have computed the free energy change Deltamicro(red) associated with one-electron reduction of the flavin ring immersed in water by utilizing the quantum mechanical/molecular mechanical method combined with the theory of energy representation (QM/MM-ER method) recently developed. As a novel treatment in implementing the QM/MM-ER method, we have identified the excess charge to be attached on the flavin ring as a solute while the remaining molecules, i.e., flavin ring and surrounding water molecules, are treated as solvent species. Then, the reduction free energy can be decomposed into the contribution Deltamicro(red)(QM) due to the oxidant described quantum chemically and the free energy Deltamicro(red)(MM) due to the water molecules represented by a classical model. By the sum of these contributions, the total reduction free energy Deltamicro(red) has been given as -80.1 kcal/mol. To examine the accuracy and efficiency of this approach, we have also conducted the Deltamicro(red) calculation using the conventional scheme that Deltamicro(red) is constructed from the solvation free energies of the flavin rings at the oxidized and reduced states. The conventional scheme has been implemented with the QM/MM-ER method and the calculated Deltamicro(red) has been estimated as -81.0 kcal/mol, showing excellent agreement with the value given by the new approach. The present approach is efficient, in particular, to compute free energy change for the reaction occurring in a protein since it enables ones to circumvent the numerical problem brought about by subtracting the huge solvation free energies of the proteins in two states before and after the reduction.