The nuclear-electronic orbital ͑NEO͒ method for the calculation of mixed nuclear-electronic wave functions is presented. Both electronic and nuclear molecular orbitals are expressed as linear combinations of Gaussian basis functions. In the NEO-HF ͑Hartree-Fock͒ method, the energy corresponding to the single-configurational mixed nuclear-electronic wave function is minimized with respect to the molecular orbitals. Multiconfigurational approaches are implemented to include significant correlation effects. In the NEO-CI ͑configuration interaction͒ method, the energy corresponding to the multiconfigurational mixed nuclear-electronic wave function is minimized with respect to the CI coefficients. In the NEO-MCSCF ͑multiconfigurational self-consistent-field͒ method, the energy is minimized with respect to the molecular orbitals as well as the CI coefficients. Analytic gradient expressions are presented for NEO-HF and NEO-MCSCF. These analytic gradients allow the variational optimization of the centers of the nuclear basis functions. They also enable the location and characterization of geometry stationary points and the generation of minimum energy paths and dynamic reaction paths. The advantages of the NEO approach are that nuclear quantum effects are incorporated during the electronic structure calculation, the Born-Oppenheimer separation of electrons and nuclei is avoided, excited vibrational-electronic states may be calculated, and its accuracy may be improved systematically. Initial applications are presented to illustrate the computational feasibility and accuracy of this approach.
The quantum dynamics of the hydride transfer reaction catalyzed by liver alcohol dehydrogenase (LADH) are studied with real-time dynamical simulations including the motion of the entire solvated enzyme. The electronic quantum effects are incorporated with an empirical valence bond potential, and the nuclear quantum effects of the transferring hydrogen are incorporated with a mixed quantum/classical molecular dynamics method in which the transferring hydrogen nucleus is represented by a three-dimensional vibrational wave function. The equilibrium transition state theory rate constants are determined from the adiabatic quantum free energy profiles, which include the free energy of the zero point motion for the transferring nucleus. The nonequilibrium dynamical effects are determined by calculating the transmission coefficients with a reactive flux scheme based on real-time molecular dynamics with quantum transitions (MDQT) surface hopping trajectories. The values of nearly unity for these transmission coefficients imply that nonequilibrium dynamical effects such as barrier recrossings are not dominant for this reaction. The calculated deuterium and tritium kinetic isotope effects for the overall rate agree with experimental results. These simulations elucidate the fundamental nature of the nuclear quantum effects and provide evidence of hydrogen tunneling in the direction along the donor-acceptor axis. An analysis of the geometrical parameters during the equilibrium and nonequilibrium simulations provides insight into the relation between specific enzyme motions and enzyme activity. The donor-acceptor distance, the catalytic zinc-substrate oxygen distance, and the coenzyme (NAD(+)/NADH) ring angles are found to strongly impact the activation free energy barrier, while the donor-acceptor distance and one of the coenzyme ring angles are found to be correlated to the degree of barrier recrossing. The distance between VAL-203 and the reactive center is found to significantly impact the activation free energy but not the degree of barrier recrossing. This result indicates that the experimentally observed effect of mutating VAL-203 on the enzyme activity is due to the alteration of the equilibrium free energy difference between the transition state and the reactant rather than nonequilibrium dynamical factors. The promoting motion of VAL-203 is characterized in terms of steric interactions involving THR-178 and the coenzyme.
A hybrid approach for simulating proton and hydride transfer reactions in enzymes is presented. The electronic quantum effects are incorporated with an empirical valence bond approach. The nuclear quantum effects of the transferring hydrogen are included with a mixed quantum/classical molecular dynamics method in which the hydrogen nucleus is described as a multidimensional vibrational wave function. The free energy profiles are obtained as functions of a collective reaction coordinate. A perturbation formula is derived to incorporate the vibrationally adiabatic nuclear quantum effects into the free energy profiles. The dynamical effects are studied with the molecular dynamics with quantum transitions (MDQT) surface hopping method, which incorporates nonadiabatic transitions among the adiabatic hydrogen vibrational states. The MDQT method is combined with a reactive flux approach to calculate the transmission coefficient and to investigate the real-time dynamics of reactive trajectories. This hybrid approach includes nuclear quantum effects such as zero point energy, hydrogen tunneling, and excited vibrational states, as well as the dynamics of the complete enzyme and solvent. The nuclear quantum effects are incorporated during the generation of the free energy profiles and dynamical trajectories rather than subsequently added as corrections. Moreover, this methodology provides detailed mechanistic information at the molecular level and allows the calculation of rates and kinetic isotope effects. An initial application of this approach to the enzyme liver alcohol dehydrogenase is also presented.
The OH radical is one of the most important oxidants in the atmosphere due to its high reactivity. The study of hydrogen-bonded complexes of OH with the water molecules is a topic of significant current interest. In this work, we present the development of a new analytical functional form for the interaction potential between the rigid OH radical and H(2)O molecules. To do this we fit a selected functional form to a set of high level ab initio data. Since there is a low-lying excited state for the H(2)O.OH complex, the impact of the excited state on the chemical behavior of the OH radical can be very important. We perform a potential energy surface scan using the CCSD(T)/aug-cc-pVTZ level of electronic structure theory for both excited and ground states. To model the physics of the unpaired electron in the OH radical, we develop a tensor polarizability generalization of the Thole-type all-atom polarizable rigid potential for the OH radical, which effectively describes the interaction of OH with H(2)O for both ground and excited states. The stationary points of (H(2)O)(n)OH clusters were identified as a benchmark of the potential.
The methodology for a vibrational analysis within the nuclear-electronic orbital ͑NEO͒ framework is presented. In the NEO approach, specified nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wave functions are calculated variationally with molecular orbital methods. Both electronic and nuclear molecular orbitals are expressed as linear combinations of Gaussian basis functions. The NEO potential energy surface depends on only the classical nuclei, and each point on this surface is optimized variationally with respect to all molecular orbitals as well as the centers of the nuclear basis functions. The NEO vibrational analysis involves the calculation, projection, and diagonalization of a numerical Hessian to obtain the harmonic vibrational frequencies corresponding to the classical nuclei. This analysis allows the characterization of stationary points on the NEO potential energy surface. It also enables the calculation of zero point energy corrections and thermodynamic properties such as enthalpy, entropy, and free energy for chemical reactions on the NEO potential energy surface. Illustrative applications of this vibrational analysis to a series of molecules and to a nucleophilic substitution reaction are presented.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.