A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations

Simulation of the enzyme reaction mechanism of malate dehydrogenase

A hybrid numerical method, which employs molecular mechanics to describe the bulk of the solvent-protein matrix and a semiempirical quantum-mechanical treatment for atoms near the reactive site, was utilized to simulate the minimum energy surface and reaction pathway for the interconversion of malate and oxaloacetate catalyzed by the enzyme malate dehydrogenase (MDH). A reaction mechanism for proton and hydride transfers associated with MDH and cofactor nicotinamide adenine dinucleotide (NAD) is deduced from the topology of the calculated energy surface. The proposed mechanism consists of (1) a sequential reaction with proton transfer preceding hydride transfer (malate to oxaloacetate direction), (2) the existence of two transition states with energy barriers of approximately 7 and 15 kcal/mol for the proton and hydride transfers, respectively, and (3) reactant (malate) and product (oxaloacetate) states that are nearly isoenergetic. Simulation analysis of the calculated energy profile shows that solvent effects due to the protein matrix dramatically alter the intrinsic reactivity of the functional groups involved in the MDH reaction, resulting in energetics similar to that found in aqueous solution. An energy decomposition analysis indicates that specific MDH residues (Arg-81, Arg-87, Asn-119, Asp-150, and Arg-153) in the vicinity of the substrate make significant energetic contributions to the stabilization of proton transfer and destabilization of hydride transfer. This suggests that these amino acids play an important role in the catalytic properties of MDH.

Structure and dynamics of a proton wire: a theoretical study of H+ translocation along the single-file water chain in the gramicidin A channel

The rapid translocation of H+ along a chain of hydrogen-bonded water molecules, or proton wire, is thought to be an important mechanism for proton permeation through transmembrane channels. Computer simulations are used to study the properties of the proton wire formed by the single-file waters in the gramicidin A channel. The model includes the polypeptidic dimer, with 22 water molecules and one excess proton. The dissociation of the water molecules is taken into account by the "polarization model" of Stillinger and co-workers. The importance of quantum effects due to the light mass of the hydrogen nuclei is examined with the use of discretized Feynman path integral molecular dynamics simulations. Results show that the presence of an excess proton in the pore orients the single-file water molecules and affects the geometry of water-water hydrogen bonding interactions. Rather than a well-defined hydronium ion OH3+ in the single-file region, the protonated species is characterized by a strong hydrogen bond resembling that of O2H5+. The quantum dispersion of protons has a small but significant effect on the equilibrium structure of the hydrogen-bonded water chain. During classical trajectories, proton transfer between consecutive water molecules is a very fast spontaneous process that takes place in the subpicosecond time scale. The translocation along extended regions of the chain takes place neither via a totally concerted mechanism in which the donor-acceptor pattern would flip over the entire chain in a single step, nor via a succession of incoherent hops between well-defined intermediates. Rather, proton transfer in the wire is a semicollective process that results from the subtle interplay of rapid hydrogen-bond length fluctuations along the water chain. These rapid structural fluctuations of the protonated single file of waters around an average position and the slow movements of the average position of the excess proton along the channel axis occur on two very different time scales. Ultimately, it is the slow reorganization of hydrogen bonds between single-file water molecules and channel backbone carbonyl groups that, by affecting the connectivity and the dynamics of the single-file water chain, also limits the translocation of the proton across the pore.

Model assembly study of the ligand binding by p-hydroxybenzoate hydroxylase: correlation between the calculated binding energies and the experimental dissociation constants

The energies of binding of seven ligands by p-hydroxybenzoate hydroxylase (PHBH) were calculated theoretically. Direct enzyme-ligand interaction energies were calculated using the ab initio quantum mechanical model assembly of the active site at the 3-21G level. Solvation energies of the ligands needed in the evaluation of the binding energies were calculated with the semiempirical AM1-SM2 method and the long-range electrostatic interaction energies between the ligands and the protein matrix classically using the static charge distributions of the ligands and the protein. Energies for proton-transfer between the ligands' OH or SH substituent at position 4 and the active-site tyrosine within the ab initio model assemblies were calculated and compared to the corresponding pKas in aqueous solution. Excluding 3,4-dihydroxybenzoate, the natural product of PHBH, a linear relationship between the calculated binding energies and the experimental binding free energies was found with a correlation coefficient of 0.90. Contributions of the direct enzyme-ligand interaction energies, solvation energies and the long-range electrostatic interaction energies to the calculated binding energies were analyzed. The proton-transfer energies of the ligands with substituents ortho to the ionized OH were found to be perturbed less in the model calculations than the energies of their meta isomers as deduced from the corresponding pKas.

The hydration and solvent polarization effects of nucleotide bases

A combined Monte Carlo quantum mechanical and molecular mechanical (QM/MM) simulation method is used to determine the free energy of hydration and the solvent polarization effect for the nucleotide bases. In the present AM1/TIP3P model, the solute molecule is characterized by valence electrons and effective nucleus cores with Hartree-Fock molecular orbital theory incorporating a solute-solvent interaction Hamiltonian. It is found that polarization energy contributes up to 37%-61% of the total solute-solvent interaction for the systems considered. The computed free energies of hydration are compared with previous theoretical results.