Ab initio molecular dynamics-based assignment of the protonation state of pepstatin A/HIV-1 protease cleavage site

A recent 13C NMR experiment (Smith et al. Nature Struct. Biol. 1996, 3, 946-950) on the Asp 25-Asp25' dyad in pepstatin A/HIV-1 protease measured two separate resonance lines, which were interpreted as being a singly protonated dyad. We address this issue by performing ab initio molecular dynamics calculations on models for this site accompanied by calculations of 13C NMR chemical shifts and isotopic shifts. We find that already on the picosecond time-scale the model proposed by Smith et al. is not stable and evolves toward a different monoprotonated form whose NMR pattern differs from the experimental one. We suggest, instead, a different protonation state in which both aspartic groups are protonated. Despite the symmetric protonation state, the calculated 13C NMR properties are in good agreement with the experiment. We rationalize this result using a simple valence bond model, which explains the chemical inequality of the two C sites. The model calculations, together with our calculations on the complex, allow also the rationalization of 13C NMR properties on other HIV-1 PR/inhibitor complexes. Both putative binding of the substrate to the free enzyme, which has the dyad singly protonated (Piana, S.; Carloni, P. Proteins: Struct., Funct., Genet. 2000, 39, 26-36), and pepstatin A binding to the diprotonated form are consistent with the inverse solvent isotope effect on the onset of inhibition of pepsin by pepstatin and the kinetic iso-mechanism proposed for aspartic proteases (Cho, T.-K.; Rebholz, K.; Northrop, D.B. Biochemistry 1994, 33, 9637-9642).     

Hydride transfer catalyzed by xylose isomerase: mechanism and quantum effects

We have applied molecular dynamics umbrella-sampling simulation and ensemble-averaged variational transition state theory with multidimensional tunneling (EA-VTST/MT) to calculate the reaction rate of xylose-to- xylulose isomerization catalyzed by xylose isomerase in the presence of two Mg2+ ions. The calculations include determination of the free energy of activation profile and ensemble averaging in the transmission coefficient. The potential energy function is approximated by a combined QM/MM/SVB method involving PM3 for the quantum mechanical (QM) subsystem, CHARMM22 and TIP3P for the molecular mechanical (MM) environment, and a simple valence bond (SVB) local function of two bond distances for the hydride transfer reaction. The simulation confirms the essential features of a mechanism postulated on the basis of kinetics and X-ray data by Whitlow et al. (Whitlow, M.; Howard, A. J.; Finzel, B. C.; Poulos, T. L.; Winborne, E.; Gilliland, G. L. Proteins 1991, 9, 153) and Ringe, Petsko, and coworkers (Labie, A.; Allen, K.-N.; Petsko, G. A.; Ringe, D. Biochemistry 1994, 33, 5469). This mechanism involves a rate-determining 1,2-hydride shift with prior and post proton transfers. Inclusion of quantum mechanical vibrational energy is important for computing the free energy of activation, and quantum mechanical tunneling effects are essential for computing kinetic isotope effects (KIEs). It is found that 85% of the reaction proceeds by tunneling and 15% by overbarrier events. The computed KIE for the ratio of hydride to deuteride transfer is in good agreement with the experimental results. The molecular dynamics simulations reveal that proton and hydride transfer reactions are assisted by breathing motions of the mobile Mg2+ ion in the active site, providing evidence for concerted motion of Mg2+ during the hydride transfer step.     

Deep tunneling dominates the biologically important hydride transfer reaction from NADH to FMN in morphinone reductase

The temperature dependence of the primary kinetic isotope effect (KIE), combined temperature-pressure studies of the primary KIE, and studies of the alpha-secondary KIE previously led us to infer that hydride transfer from nicotinamide adenine dinucleotide to flavin mononucleotide in morphinone reductase proceeds via environmentally coupled hydride tunneling. We present here a computational analysis of this hydride transfer reaction using QM/MM molecular dynamics simulations and variational transition-state theory calculations. Our calculated primary and secondary KIEs are in good agreement with the corresponding experimental values. Although the experimentally observed KIE lies below the semiclassical limit, our calculations suggest that approximately 99% of the reaction proceeds via tunneling: this is the first "deep tunneling" reaction observed for hydride transfer. We also show that the dominant tunneling mechanism is controlled by the isotope at the primary rather than the secondary position: with protium in the primary position, large-curvature tunneling dominates, whereas with deuterium in this position, small-curvature tunneling dominates. Also, our study is consistent with tunneling being preceded by reorganization: in the reactant, the rings of the nicotinamide and isoalloxazine moieties are stacked roughly parallel to each other, and as the system moves toward a "tunneling-ready" configuration, the nicotinamide ring rotates to become almost perpendicular to the isoalloxazine ring.     

Hydride transfer reaction catalyzed by hyperthermophilic dihydrofolate reductase is dominated by quantum mechanical tunneling and is promoted by both inter- and intramonomeric correlated motions

Simulations of hydride and deuteride transfer catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) are presented. TmDHFR was modeled with its active homodimeric quaternary structure, where each monomer has three subdomains. The potential energy function was a combined quantum mechanical and molecular mechanical potential (69 atoms were treated quantum mechanically, and 35 287, by molecular mechanics). The calculations of the rate constants by ensemble-averaged variational transition state theory with multidimensional tunneling predicted that hydride and deuteride transfer at 278 K proceeded with 81 and 80% by tunneling. These percentages decreased to 50 and 49% at 338 K. The kinetic isotope effect was dominated by contributions of bound vibrations and decreased from 3.0 to 2.2 over the temperature range. The calculated rates for hydride and deuteride transfer catalyzed by the hypothetical monomer were smaller by approximately 2 orders of magnitude. At 298 K tunneling contributed 73 and 66% to hydride and deuteride transfer in the monomer. The decreased catalytic efficiency of the monomer was therefore not the result of a decrease of the tunneling contribution but an increase in the quasi-classical activation free energy. The catalytic effect was associated in the dimer with correlated motions between domains as well as within and between subunits. The intrasubunit correlated motions were decreased in the monomer when compared to both native dimeric TmDHFR and monomeric E. coli enzyme. TmDHFR and its E. coli homologue involve similar patterns of correlated interactions that affect the free energy barrier of hydride transfer despite only 27% sequence identity and different quaternary structures.     

NMR study of C-kinetic isotope effects at C natural abundance to characterize oxidations and an enzyme-catalyzed reduction

¹³C-kinetic isotope effects (KIEs) of four cinnamyl alcohol oxidations and a xylose reductase-catalyzed cinnamyl aldehyde reduction have been determined by ¹³C NMR using competition reactions with reactants at natural ¹³C-abundance. Differences in KIEs among oxidations indicate dissimilarities between the respective hydrogen transfers. Their mechanistic implications are discussed. A low primary KIE of the enzymatic reduction is consistent with a kinetically complex mechanism in which steps other than the chemical step of hydride transfer from NADH are slow.     

Tunneling proton transfer in biological systems. Role of temperature and pressure

The possibilities of hydrogen atom tunneling transfer in biological liquids are discussed. Basic mechanisms of temperature and pressure effects on the tunneling rate constant are considered: the reorganization of reagents and the medium due to the transfer of H atoms and changes in the value and shape of the chemical reaction potential barrier upon intermolecular and soft intramolecular vibrations. Expressions are derived for the tunneling transition rate constant and kinetic isotopic effect as functions of temperature and pressure. It is found that the temperature dependence of the isotope effect is mainly affected by the second mechanism only. The theory is compared with the literature??s experimental data on the temperature dependence of the isotope effect. It is shown that experiments are described well by the theory at sensible values of the fitting parameters.     

Ground-state proton-transfer dynamics governed by configurational optimization

The ground-state proton transfer (GSPT) of 7-hydroxyquinoline along a hydrogen-bonded alcohol chain has been investigated in n-alkanes using time-resolved transient-absorption spectroscopy with variation of alcohols, media, isotopes, and temperatures. As a 7-hydroxyquinoline molecule associates with two alcohol molecules via hydrogen bonding to form a cyclic complex in a nonpolar aprotic medium, the intrinsic GSPT dynamics of the cyclic complex in a n-alkane has been observed directly without being interfered with by solvent association to form the cyclic complex. GSPT occurs concertedly without accumulating any reaction intermediate and yet asymmetrically with a rate-determining tunneling process. Both the rate constant and the kinetic isotope effect of GSPT increase rapidly with the proton-donating ability of the alcohol but decrease greatly with the molecular size of the alcohol. The reorganization of the hydrogen-bond bridge to form an optimal precursor configuration for efficient proton tunneling takes place prior to intrinsic GSPT, and configurational optimization becomes more important as the molecular size of the alcohol increases. Consequently, the larger contribution of configurational optimization to GSPT leads to the weaker asymmetric character of GSPT.     

Proton tunneling in aromatic amine dehydrogenase is driven by a short-range sub-picosecond promoting vibration: consistency of simulation and theory with experiment

Hydrogen transfer, an essential component of most biological reactions, is a quantum problem. However, the proposed role of compressive motion in promoting enzymatic H-transfer is contentious. Using molecular dynamics simulations and density functional theory (DFT) calculations, we show that, during proton tunneling in the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase (AADH), a sub-picosecond promoting vibration is inherent to the iminoquinone intermediate. We show by numerical modeling that this short-range vibration, with a frequency of approximately 165 cm-1, is consistent with "gating" motion in the hydrogen tunneling model of Kuznetsov and Ulstrup (Kuznetsov, A. M.; Ulstrup, J. Can. J. Chem. 1999, 77, 1085) in an enzymatic reaction with an observed protium/deuterium kinetic isotope effect that is not measurably temperature-dependent.     

The effect of pressure on hydrogen transfer reactions with quinones

The effect of pressure on the oxidation of hydroarenes 3-9 with 2,3-dichloro-5,6-dicyano-1,4-quinone (DDQ; 1 a) or o-chloranil (10), leading to the corresponding arenes, has been investigated. The activation volumes were determined from the pressure dependence of the rate constants of these reactions monitored by on-line UV/Vis spectroscopic measurements in an optical high-pressure cell (up to 3500 bar). The finding that they are highly negative and only moderately dependent on the solvent polarity (DeltaV( not equal ) = -13 to -25 in MTBE and -15 to -29 cm(3) mol(-1) in MeCN/AcOEt, 1:1) rules out the formation of ionic species in the rate-determining step and is good evidence for a hydrogen atom transfer mechanism leading to a pair of radicals in the rate-determining step, as was also suggested by kinetic measurements, studies of kinetic isotope effects, and spin-trapping experiments. The strong pressure dependence of the kinetic deuterium isotope effect for the reaction of 9,10-dihydroanthracene 5/5-9,9,10,10-D(4) with DDQ (1 a) can be attributed to a tunneling component in the hydrogen transfer. In the case of formal 1,3-dienes and enes possessing two vicinal C--H bonds, which have to be cleaved during the dehydrogenation, a pericyclic hydrogen transfer has to considered as one mechanistic alternative. The comparison of the kinetic deuterium isotope effects determined for the oxidation of tetralin 9/9-1,1,4,4-D(4)/9-2,2,3,3-D(4)/9-D(12) either with DDQ (1 a) or with thymoquinone 1 c indicates that the reaction with DDQ (1 a) proceeds in a stepwise manner through hydrogen atom transfer, analogously to the oxidations of 1,4-dihydroarenes, whereas the reaction with thymoquinone 1 c is concerted, following the course of a pericyclic hydrogen transfer. The difference in the mechanistic courses of these two reactions may be explained by the effect of the CN and Cl substituents in 1 a, which stabilize a radical intermediate better than the alkyl groups in 1 c. The mechanistic conclusions are substantiated by DFT calculations.     

Mutagenesis of morphinone reductase induces multiple reactive configurations and identifies potential ambiguity in kinetic analysis of enzyme tunneling mechanisms

We have identified multiple reactive configurations (MRCs) of an enzyme-coenzyme complex that have measurably different kinetic properties. In the complex formed between morphinone reductase (MR) and the NADH analogue 1,4,5,6-tetrahydro-NADH (NADH4) the nicotinamide moiety is restrained close to the FMN isoalloxazine ring by hydrogen bonds from Asn-189 and His-186 as determined from the X-ray crystal structure. Molecular dynamic simulations indicate that removal of one of these hydrogen bonds in the N189A MR mutant allows the nicotinamide moiety to occupy a region of configurational space not accessible in wild-type enzyme. Using stopped-flow spectroscopy, we show that reduction of the FMN cofactor by NADH in N189A MR is multiphasic, identifying at least four different reactive configurations of the MR-NADH complex. This contrasts with wild-type MR in which hydride transfer occurs by environmentally coupled tunneling in a single kinetic phase [Pudney et al. J. Am. Chem. Soc. 2006, 128, 14053-14058]. Values for primary and alpha-secondary kinetic isotope effects, and their temperature dependence, for three of the kinetic phases in the N189A MR are consistent with hydride transfer by tunneling. Our analysis enables derivation of mechanistic information concerning different reactive configurations of the same enzyme-coenzyme complex using ensemble stopped-flow methods. Implications for the interpretation from kinetic data of tunneling mechanisms in enzymes are discussed.     

Hydrogen tunneling in peptidylglycine alpha-hydroxylating monooxygenase

The temperature dependence of the primary and secondary intrinsic isotope effects for the C-H bond cleavage catalyzed by peptidylglycine alpha-hydroxylating monooxygenase has been determined. Analysis of the magnitude and Arrhenius behavior of the intrinsic isotope effects provides strong evidence for the use of tunneling as a primary catalytic strategy for this enzyme. Modeling of the isotope effect data allows for a comparison to the hydrogen transfer catalyzed by soybean lipoxygenase in terms of environmental reorganization energy and frequency of the protein vibration that controls the hydrogen transfer.     

Effective thermal oxidation of isopropanol by an NAD model

The reaction of 10-methylacridinium cation (MA⁺) with isopropanol in the parent alcohol medium under dark, oxygen-free, and refluxing conditions gave hydride transfer product 10-methyl-9,10-dihydroacridine (MAH). The kinetics of the alcoholic oxidation reaction, including the kinetic isotope effect and the kinetic temperature effect, were determined. Hydride transfer is involved in the rate-determining step.     

Intrinsic isotope effects on benzylic hydroxylation by the aromatic amino acid hydroxylases: evidence for hydrogen tunneling, coupled motion, and similar reactivities

Deuterium kinetic isotope effects for hydroxylation of the methyl group of 4-methylphenylalanine have been used as a probe of the relative reactivities of the hydroxylating intermediates in the aromatic amino acid hydroxylases phenylalanine, tyrosine, and tryptophan hydroxylase. When there are three deuterium atoms in the methyl group, all three enzymes exhibit an intrinsic isotope effect of about 13. The temperature dependence of the isotope effect is consistent with moderate tunneling, with the extent of tunneling identical for all three enzymes. In the case of phenylalanine hydroxylase, the presence of the regulatory domain has no effect on the values. The intrinsic primary and secondary isotope effects were determined using 4-methylphenylalanine containing one or two deuterium atoms in the methyl group. With one deuterium atom, the intrinsic primary and secondary effects have average values of 10 and 1.1, respectively. With two deuterium atoms, the primary effects decrease to 7.4 and the secondary effect increases to 1.3, consistent with coupled motion of the primary and secondary hydrogens. The results with all three enzymes are consistent with a hydrogen abstraction mechanism. The similarities of the isotope effects and extent of tunneling establish that the reactivities of the hydroxylating intermediates in the three enzymes are essentially identical.     

Mechanism investigation of ketone hydrogenation catalyzed by ruthenium bifunctional catalysts: insights from a DFT study

In this paper, the mechanism of ketone hydrogenation catalyzed by five Ru bifunctional catalysts with different structural frameworks was studied in detail using density functional theory (DFT). This mechanism contains hydrogen transfer, dehydrogenation of alcohol, and dihydrogen activation fundamental reactions. The involvement of alcohol is also discussed and found with different activities in hydrogen transfer, dehydrogenation and dihydrogen activation steps in five systems. Our calculated results indicate that the weak Ru-H bond, stronger basicity of hydride and stronger X-H acidity will decrease the barrier of the HT step, and that the polar micro-environment of dihydrogen coordinating with Ru catalysts and short hydrogen transfer distance would be able to facilitate the heterolytic splitting of dihydrogen in the dihydrogen activation step.     

Kinetic isotope effect evidence for a concerted hydrogen transfer mechanism in transfer hydrogenations catalyzed by [p-(Me2CH)C6H4Me]Ru- (NHCHPhCHPhNSO2C6H4-p-CH3)

The isotope effects in the reaction of [p-(Me2CH)C6H4Me]Ru(NHCHPhCHPhNSO2C6H4-p-CH3) (1) with isopropyl alcohol were 1.79 for transfer of hydrogen from OH to N and 2.86 for transfer from CH to Ru. The isotope effect for transfer of deuterium from doubly labeled material (kCHOH/kCDOD = 4.88) was within experimental error of the product of the two individual isotope effects. These isotope effects provide convincing evidence for a mechanism involving concurrent transfer of hydrogen from oxygen to nitrogen and from carbon to ruthenium.     

Tunneling transfer of atomic particles in chemical and biological reactions: The role of intermolecular vibrations and media reorganization

Various types of chemical and biological tunneling reactions in a condensed phase are discussed. The analytical expressions for the rate constants in different temperature ranges are given. Experimental data on such low-temperature processes as hydrogen transfer from a molecule to a radical between two molecules and intramolecular transformations are considered. Data on the kinetic isotope effect upon the transfer of atomic particles in the solid phase and biological liquids are presented. The effect pressure has on different tunneling reactions is also considered; where possible, experimental results are compared with theory.     

Small temperature dependence of the kinetic isotope effect for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase

The H/D primary kinetic isotope effect (KIE) for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase (ecDHFR) is calculated as a function of temperature employing ensemble-averaged variational transition-state theory with multidimensional tunneling. The calculated KIEs display only a small temperature dependence over the temperature range of 5 to 45 degrees C. We identify two key features that contribute to canceling most of the temperature dependence of the KIE that would be expected on the basis of simpler models. Related issues such as the isotope effects on Arrhenius preexponential factors, large differences between free energies of activation and Arrhenius activation energy, and fluctuations of effective barriers are also discussed.     

Modulation of ligand-field parameters by heme ruffling in cytochromes c revealed by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectra of variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht c-552) and Pseudomonas aeruginosa cytochrome c(551) (Pa c-551) are analyzed to determine the effect of heme ruffling on ligand-field parameters. Mutations introduced at positions 13 and 22 in Ht c-552 were previously demonstrated to influence hydrogen bonding in the proximal heme pocket and to tune reduction potential (E(m)) over a range of 80 mV [Michel, L. V.; Ye, T.; Bowman, S. E. J.; Levin, B. D.; Hahn, M. A.; Russell, B. S.; Elliott, S. J.; Bren, K. L. Biochemistry 2007, 46, 11753-11760]. These mutations are shown here to also increase heme ruffling as E(m) decreases. The primary effect on electronic structure of increasing heme ruffling is found to be a decrease in the axial ligand-field term Δ/λ, which is proposed to arise from an increase in the energy of the d(xy) orbital. Mutations at position 7, previously demonstrated to influence heme ruffling in Pa c-551 and Ht c-552, are utilized to test this correlation between molecular and electronic structure. In conclusion, the structure of the proximal heme pocket of cytochromes c is shown to play a role in determining heme conformation and electronic structure.