Browsing by Author "Ojeda-May, Pedro"

Now showing 1 - 4 of 4
  • Thumbnail Image
    Item
    Assessing the accuracy of the isotropic periodic sum method through Madelung energy computation
    (2014-04) Ojeda-May, Pedro; Pu, Jingzhi
    We tested the isotropic periodic sum (IPS) method for computing Madelung energies of ionic crystals. The performance of the method, both in its nonpolar (IPSn) and polar (IPSp) forms, was compared with that of the zero-charge and Wolf potentials [D. Wolf, P. Keblinski, S. R. Phillpot, and J. Eggebrecht, J. Chem. Phys.110, 8254 (1999)]. The results show that the IPSn and IPSp methods converge the Madelung energy to its reference value with an average deviation of ∼10−4 and ∼10−7 energy units, respectively, for a cutoff range of 18–24a (a/2 being the nearest-neighbor ion separation). However, minor oscillations were detected for the IPS methods when deviations of the computed Madelung energies were plotted on a logarithmic scale as a function of the cutoff distance. To remove such oscillations, we introduced a modified IPSn potential in which both the local-region and long-range electrostatic terms are damped, in analogy to the Wolf potential. With the damped-IPSn potential, a smoother convergence was achieved. In addition, we observed a better agreement between the damped-IPSn and IPSp methods, which suggests that damping the IPSn potential is in effect similar to adding a screening potential in IPSp.
  • Thumbnail Image
    Item
    Mapping Free Energy Pathways for ATP Hydrolysis in the E. coli ABC Transporter HlyB by the String Method
    (MDPI, 2018-10-16) Zhou, Yan; Ojeda-May, Pedro; Nagaraju, Mulpuri; Kim, Bryant; Pu, Jingzhi; Chemistry and Chemical Biology, School of Science
    HlyB functions as an adenosine triphosphate (ATP)-binding cassette (ABC) transporter that enables bacteria to secrete toxins at the expense of ATP hydrolysis. Our previous work, based on potential energy profiles from combined quantum mechanical and molecular mechanical (QM/MM) calculations, has suggested that the highly conserved H-loop His residue H662 in the nucleotide binding domain (NBD) of E. coli HlyB may catalyze the hydrolysis of ATP through proton relay. To further test this hypothesis when entropic contributions are taken into account, we obtained QM/MM minimum free energy paths (MFEPs) for the HlyB reaction, making use of the string method in collective variables. The free energy profiles along the MFEPs confirm the direct participation of H662 in catalysis. The MFEP simulations of HlyB also reveal an intimate coupling between the chemical steps and a local protein conformational change involving the signature-loop residue S607, which may serve a catalytic role similar to an Arg-finger motif in many ATPases and GTPases in stabilizing the phosphoryl-transfer transition state.
  • Thumbnail Image
    Item
    Reaction Path-Force Matching in Collective Variables: Determining Ab Initio QM/MM Free Energy Profiles by Fitting Mean Force
    (American Chemical Society, 2021) Kim, Bryant; Snyder, Ryan; Nagaraju, Mulpuri; Zhou, Yan; Ojeda-May, Pedro; Keeton, Seth; Hege, Mellisa; Shao, Yihan; Pu, Jingzhi; Chemistry and Chemical Biology, School of Science
    First-principles determination of free energy profiles for condensed-phase chemical reactions is hampered by the daunting costs associated with configurational sampling on ab initio quantum mechanical/molecular mechanical (AI/MM) potential energy surfaces. Here, we report a new method that enables efficient AI/MM free energy simulations through mean force fitting. In this method, a free energy path in collective variables (CVs) is first determined on an efficient reactive aiding potential. Based on the configurations sampled along the free energy path, correcting forces to reproduce the AI/MM forces on the CVs are determined through force matching. The AI/MM free energy profile is then predicted from simulations on the aiding potential in conjunction with the correcting forces. Such cycles of correction-prediction are repeated until convergence is established. As the instantaneous forces on the CVs sampled in equilibrium ensembles along the free energy path are fitted, this procedure faithfully restores the target free energy profile by reproducing the free energy mean forces. Due to its close connection with the reaction path-force matching (RP-FM) framework recently introduced by us, we designate the new method as RP-FM in collective variables (RP-FM-CV). We demonstrate the effectiveness of this method on a type-II solution-phase SN2 reaction, NH3 + CH3Cl (the Menshutkin reaction), simulated with an explicit water solvent. To obtain the AI/MM free energy profiles, we employed the semiempirical AM1/MM Hamiltonian as the base level for determining the string minimum free energy pathway, along which the free energy mean forces are fitted to various target AI/MM levels using the Hartree-Fock (HF) theory, density functional theory (DFT), and the second-order Møller-Plesset perturbation (MP2) theory as the AI method. The forces on the bond-breaking and bond-forming CVs at both the base and target levels are obtained by force transformation from Cartesian to redundant internal coordinates under the Wilson B-matrix formalism, where the linearized FM is facilitated by the use of spline functions. For the Menshutkin reaction tested, our FM treatment greatly reduces the deviations on the CV forces, originally in the range of 12-33 to ∼2 kcal/mol/Å. Comparisons with the experimental and benchmark AI/MM results, tests of the new method under a variety of simulation protocols, and analyses of the solute-solvent radial distribution functions suggest that RP-FM-CV can be used as an efficient, accurate, and robust method for simulating solution-phase chemical reactions.
  • Thumbnail Image
    Item
    Treating electrostatics with Wolf summation in combined quantum mechanical and molecular mechanical simulations
    (AIP, 2015-11) Ojeda-May, Pedro; Pu, Jingzhi; Department of Chemistry & Chemical Biology, IU School of Science
    The Wolf summation approach [D. Wolf et al., J. Chem. Phys. 110, 8254 (1999)], in the damped shifted force (DSF) formalism [C. J. Fennell and J. D. Gezelter, J. Chem. Phys. 124, 234104 (2006)], is extended for treating electrostatics in combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulations. In this development, we split the QM/MM electrostatic potential energy function into the conventional Coulomb r −1 term and a term that contains the DSF contribution. The former is handled by the standard machinery of cutoff-based QM/MM simulations whereas the latter is incorporated into the QM/MM interaction Hamiltonian as a Fock matrix correction. We tested the resulting QM/MM-DSF method for two solution-phase reactions, i.e., the association of ammonium and chloride ions and a symmetric SN2 reaction in which a methyl group is exchanged between two chloride ions. The performance of the QM/MM-DSF method was assessed by comparing the potential of mean force (PMF) profiles with those from the QM/MM-Ewald and QM/MM-isotropic periodic sum (IPS) methods, both of which include long-range electrostatics explicitly. For ion association, the QM/MM-DSF method successfully eliminates the artificial free energy drift observed in the QM/MM-Cutoff simulations, in a remarkable agreement with the two long-range-containing methods. For the SN2 reaction, the free energy of activation obtained by the QM/MM-DSF method agrees well with both the QM/MM-Ewald and QM/MM-IPS results. The latter, however, requires a greater cutoff distance than QM/MM-DSF for a proper convergence of the PMF. Avoiding time-consuming lattice summation, the QM/MM-DSF method yields a 55% reduction in computational cost compared with the QM/MM-Ewald method. These results suggest that, in addition to QM/MM-IPS, the QM/MM-DSF method may serve as another efficient and accurate alternative to QM/MM-Ewald for treating electrostatics in condensed-phase simulations of chemical reactions.