A real-space-grid QM/MM study on the ionic/radical association reaction in aqueous phase: HCHO+OH→HCHO–OH

A critical comparison of widely used solvation models is reported. It is illustrated by a study of the triazene molecule in liquid water. We consider the following approaches: (1) a continuum model based on multicentric multipole expansions of the charge distribution, (2) the averaged solvent electrostatic potential from molecular dynamics (ASEP/MD) method, and (3) molecular dynamics simulations using a combined quantum mechanics/molecular mechanics potential (QM/MM/MD). We find that the solvation induces appreciable changes in the geometry and charge distribution of triazene. These changes are only qualitatively reproduced by the dielectric continuum model, which clearly underestimates induced dipole moments and solute-solvent interaction energy. We also show that the use of effective point charges placed on solute nuclei during the classical simulations may cause significant errors in the description of the solvent structure. The addition of charges representing nitrogen atom lone pairs is compulsory to reproduce the QM/MM/MD simulation results. Moreover, our results validate the use of the mean field approximation in the study of solvent effects. A major conclusion of this study is that the ASEP/MD method constitutes a reliable alternative to the much more computationally demanding QM/MM/MD methods.

Section: Radial Distribution Functionssupporting

confidence: 80%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Comparison of three effective Hamiltonian models of increasing complexity: Triazene in water as a test case

Galván

Martı́n

Aguilar

et al. 2006

“…Similar observations have been made in literature. 33,34,73,74 Nevertheless, we stress that it is not clear at present to what extent this conclusion applies to different types of systems, e.g., enzyme reactions where local fluctuations of the MM environment may deviate significantly from the Gaussian distribution. 50…”

Section: Free Energy Profilesmentioning

confidence: 82%

Variational and perturbative formulations of quantum mechanical/molecular mechanical free energy with mean-field embedding and its analytical gradients

Yamamoto

2008

Conventional quantum chemical solvation theories are based on the mean-field embedding approximation. That is, the electronic wavefunction is calculated in the presence of the mean field of the environment. In this paper a direct quantum mechanical/molecular mechanical (QM/MM) analog of such a mean-field theory is formulated based on variational and perturbative frameworks. In the variational framework, an appropriate QM/MM free energy functional is defined and is minimized in terms of the trial wavefunction that best approximates the true QM wavefunction in a statistically averaged sense. Analytical free energy gradient is obtained, which takes the form of the gradient of effective QM energy calculated in the averaged MM potential. In the perturbative framework, the above variational procedure is shown to be equivalent to the first-order expansion of the QM energy (in the exact free energy expression) about the self-consistent reference field. This helps understand the relation between the variational procedure and the exact QM/MM free energy as well as existing QM/MM theories. Based on this, several ways are discussed for evaluating non-mean-field effects (i.e., statistical fluctuations of the QM wavefunction) that are neglected in the mean-field calculation. As an illustration, the method is applied to an S(N)2 Menshutkin reaction in water, NH(3)+CH(3)Cl-->NH(3)CH(3) (+)+Cl(-), for which free energy profiles are obtained at the Hartree-Fock, MP2, B3LYP, and BHHLYP levels by integrating the free energy gradient. Non-mean-field effects are evaluated to be <0.5 kcal/mol using a Gaussian fluctuation model for the environment, which suggests that those effects are rather small for the present reaction in water.

“…[31][32][33][34][35][36][37] As a notable feature of the code, the one-electron wave functions in the QM subsystem are represented by the real-space grids that are uni-FIG. [31][32][33][34][35][36][37] As a notable feature of the code, the one-electron wave functions in the QM subsystem are represented by the real-space grids that are uni-FIG.…”

Section: B Qm/mm Simulationsmentioning

confidence: 99%

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

Takahashi

Ohno

Kishi

et al. 2008

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.