Mechanistic aspects of biological redox reactions involving NADH 2: A combined semiempirical andab initio study of hydride-ion transfer between the NADH analogue, 1-methyl-dihydronicotinamide, and folate and dihydrofolate analogue substrates of dihydrofolate reductase

Evolution of a transition state: role of Lys100 in the active site of isocitrate dehydrogenase

An active site lysine essential to catalysis in isocitrate dehydrogenase (IDH) is absent from related enzymes. As all family members catalyze the same oxidative β-decarboxylation at the (2R)-malate core common to their substrates, it seems odd that an amino acid essential to one is not found in all. Ordinarily, hydride transfer to a nicotinamide C4 neutralizes the positive charge at N1 directly. In IDH, the negatively charged C4-carboxylate of isocitrate stabilizes the ground state positive charge on the adjacent nicotinamide N1, opposing hydride transfer. The critical lysine is poised to stabilize-and perhaps even protonate-an oxyanion formed on the nicotinamide 3-carboxamide, thereby enabling the hydride to be transferred while the positive charge at N1 is maintained. IDH might catalyze the same overall reaction as other family members, but dehydrogenation proceeds through a distinct, though related, transition state. Partial activation of lysine mutants by K(+) and NH4 (+) represents a throwback to the primordial state of the first promiscuous substrate family member.

Hybrid Quantum and Classical Simulations of the Dihydrofolate Reductase Catalyzed Hydride Transfer Reaction on an Accurate Semi-Empirical Potential Energy Surface

Dihydrofolate reductase (DHFR) catalyzes the reduction of 7,8-dihydrofolate by nicotinamide adenine dinucleotide phosphate hydride (NADPH) to form 5,6,7,8-tetrahydrofolate and oxidized nicotinamide. DHFR is a small, flexible, monomeric protein with no metals or SS bonds and serves as one of the enzymes commonly used to examine basic aspects in enzymology. In the current work, we present extensive benchmark calculations for several model reactions in the gas phase that are relevant to the DHFR catalyzed hydride transfer. To this end, we employ G4MP2 and CBS-QB3 ab initio calculations as well as numerous density functional theory methods. Using these results, we develop two specific reaction parameter (SRP) Hamiltonians based on the semiempirical AM1 method. The first generation SRP Hamiltonian does not account for dispersion, while the second generation SRP accounts for dispersion implicitly via the AM1 core-repulsion functions. These SRP semiempirical Hamiltonians are subsequently used in hybrid quantum mechanics/molecular mechanics simulations of the DHFR catalyzed reaction. Finally, kinetic isotope effects are computed using a mass-perturbation-based path-integral approach.

A theoretical study of the catalytic mechanism of formate dehydrogenase

A theoretical study of the hydride transfer between formate anion and nicotinamide adenine dinucleotide (NAD(+)) catalyzed by the enzyme formate dehydrogenase (FDH) has been carried out by a combination of two hybrid quantum mechanics/molecular mechanics techniques: statistical simulation methods and internal energy minimizations. Free energy profiles, obtained for the reaction in the enzyme active site and in solution, allow obtaining a comparative analysis of the behavior of both condensed media. Moreover, calculations of the reaction in aqueous media can be used to probe the dramatic differences between reactants state in the enzyme active site and in solution. The results suggest that the enzyme compresses the substrate and the cofactor into a conformation close to the transition structure by means of favorable interactions with the amino acid residues of the active site, thus facilitating the relative orientation of donor and acceptor atoms to favor the hydride transfer. Moreover, a permanent field created by the protein reduces the work required to reach the transition state (TS) with a concomitant polarization of the cofactor that would favor the hydride transfer. In contrast, in water the TS is destabilized with respect to the reactant species because the polarity of the solute diminishes as the reaction proceeds, and consequently the reaction field, which is created as a response to the change in the solute polarity, is also decreased. Therefore protein structure is responsible of both effects; substrate preorganization and TS stabilization thus diminishing the activation barrier. Because of the electrostatic features of the catalyzed reaction, both media preferentially stabilize the ground-state, thus explaining the small rate constant enhancement of this enzyme, but FDH does so to a much lower extent than aqueous solution. Finally, a good agreement between experimental and theoretical kinetic isotope effects is found, thus giving some credit to our results.

Freezing a Single Distal Motion in Dihydrofolate Reductase

Constraining a single motion between distal residues separated by approximately 28 A in hybrid quantum/classical molecular dynamics simulations is found to increase the free energy barrier for hydride transfer in dihydrofolate reductase by approximately 3 kcal/mol. Our analysis indicates that a single distal constraint alters equilibrium motions throughout the enzyme on a wide range of time scales. This alteration of the conformational sampling of the entire system is sufficient to significantly increase the free energy barrier and decrease the rate of hydride transfer. Despite the changes in conformational sampling introduced by the constraint, the system assumes a similar transition state conformation with a donor-acceptor distance of approximately 2.72 A to enable the hydride transfer reaction. The modified thermal sampling leads to a substantial increase in the average donor-acceptor distance for the reactant state, however, thereby decreasing the probability of sampling the transition state conformations with the shorter distances required for hydride transfer. These simulations indicate that fast thermal fluctuations of the enzyme, substrate, and cofactor lead to conformational sampling of configurations that facilitate hydride transfer. The fast thermal motions are in equilibrium as the reaction progresses along the collective reaction coordinate, and the overall average equilibrium conformational changes occur on the slower time scale measured experimentally. Recent single molecule experiments suggest that at least some of these thermally averaged equilibrium conformational changes occur on the millisecond time scale of the hydride transfer reaction. Thus, introducing a constraint that modifies the conformational sampling of an enzyme could significantly impact its catalytic activity.

Quantum-classical simulation methods for hydrogen transfer in enzymes: a case study of dihydrofolate reductase

A variety of theoretical approaches have been used to investigate hydrogen transfer in enzymatic reactions. The free energy barriers for hydrogen transfer in enzymes have been calculated using classical molecular dynamics simulations in conjunction with quantum mechanical/molecular mechanical and empirical valence bond potentials. Nuclear quantum effects have been included with several different approaches. Applications of these approaches to hydride transfer in dihydrofolate reductase are consistent with experimental measurements and provide significant insight into the protein conformational changes that facilitate the hydride transfer reaction.

Theoretical evaluation of the hydrogen kinetic isotope effect on the first step of the methylmalonyl-CoA mutase reaction

We have calculated hydrogen kinetic isotope effects (KIEs) for the first step of the methylmalonyl-CoA mutase reaction, including multidimensional tunneling correction at the zero curvature (ZCT) level, and compared them with the experimental values. Both alternative mechanisms of this step, concerted and stepwise, can be accommodated. It turned out to be essential to include Arg207 hydrogen-bonded to the reactant in the mechanism predicting simultaneous breaking of the Co-C bond of AdoCbl and hydrogen atom transfer. The consequence of the stepwise mechanism is a much larger facilitation of the homolytic dissociation of the carbon-cobalt bond by the enzyme than currently appreciated; our results suggest lowering of the activation energy by about 23 kcal mol(-1). We have also shown that large hydrogen KIEs of tunneling origin do not necessarily break the Swain-Schaad equation. Furthermore, when this equation does not hold, the exponent may be smaller in the presence of tunneling than it is at the semi-classical limit, indicating that nonclassical behavior may be a more common phenomenon than expected.

Chlorine kinetic isotope effects on the haloalkane dehalogenase reaction

We have found chlorine kinetic isotope effects on the dehalogenation catalyzed by haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 to be 1.0045 +/- 0.0004 for 1,2-dichloroethane and 1.0066 +/- 0.0004 for 1-chlorobutane. The latter isotope effect approaches the intrinsic chlorine kinetic isotope effect for the dehalogenation step. The intrinsic isotope effect has been modeled using semiempirical and DFT theory levels using the ONIOM QM/QM scheme. Our results indicate that the dehalogenation step is reversible; the overall irreversibility of the enzyme-catalyzed reaction is brought about by a step following the dehalogenation.

Orbital steering in the catalytic power of enzymes: small structural changes with large catalytic consequences

Small structural perturbations in the enzyme isocitrate dehydrogenase (IDH) were made in order to evaluate the contribution of precise substrate alignment to the catalytic power of an enzyme. The reaction trajectory of IDH was modified (i) after the adenine moiety of nicotinamide adenine dinucleotide phosphate was changed to hypoxanthine (the 6-amino was changed to 6-hydroxyl), and (ii) by replacing Mg2+, which has six coordinating ligands, with Ca2+, which has eight coordinating ligands. Both changes make large (10(-3) to 10(-5)) changes in the reaction velocity but only small changes in the orientation of the substrates (both distance and angle) as revealed by cryocrystallographic trapping of active IDH complexes. The results provide evidence that orbital overlap produced by optimal orientation of reacting orbitals plays a major quantitative role in the catalytic power of enzymes.

Mechanism of aldehyde oxidation catalyzed by horse liver alcohol dehydrogenase

The mechanism of oxidation of benzaldehyde to benzoic acid catalyzed by horse liver alcohol dehydrogenase (HLADH) has been investigated using the HLADH structure at 2.1 A resolution with NAD+ and pentafluorobenzyl alcohol in the active site [Ramaswamy et al. (1994) Biochemistry 33,5230-5237]. Constructs for molecular dynamics (MD) investigations with HLADH were obtained by a best-fit superimposition of benzaldehyde or its hydrate on the pentafluorobenzyl alcohol bound to the active site Zn(II)ion. Equilibrium bond lengths, angles, and dihedral parameters for Zn(II) bonding residues His67, Cys46, and Cys174 were obtained from small-molecule X-ray crystal structures and an ab initio-derived parameterization of zinc in HLADH [Ryde, U. (1995) Proteins: Struct., Funct., Genet. 21,40-56]. Dynamic simulations in CHARMM were carried out on the following three constructs to 100 ps: (MD1) enzyme with NAD+, benzaldehyde, and zinc-ligated HO-in the active site; (MD2) enzyme with NAD+ and hydrated benzaldehyde monoanion bound to zinc via the pro-R oxygen, with a proton residing on the pro-S oxygen; and (MD3) enzyme with NAD+ and hydrated benzaldehyde monoanion bound to zinc via the pro-S oxygen, with a proton residing on the pro-R oxygen. Analyses were done of 800 sample conformations taken in the last 40 ps of dynamics. Structures from MD1 and MD3 were used to define the initial spatial arrangements of reactive functionalities for semiempirical PM3 calculations. Using PM3, model systems were calculated of ground states and some transition states for aldehyde hydration, hydride transfer, and subsequent proton shuttling. With benzaldehyde and zinc-bound hydroxide ion in the active site, the oxygen of Zn(II)-OH resided at a distance of 2.8-5.5 A from the aldehyde carbonyl carbon during the dynamics simulation. This may be compared to the PM3 transition state for attack of the Zn(II)-OH oxygen on the benzaldehyde carbonyl carbon, which has an O...C distance of 1.877 A. HLADH catalysis of the aldehyde hydration would require very little motion aside from that in the ground state. Two simulations of benzaldehyde hydrate ligated to zinc (MD2 and MD3) both showed close approach of the aldehyde hydrate hydrogen to NAD+C4, varying from 2.3 to 3.3 A, seemingly favorable for the hydride transfer reaction. The MD2 configuration does not allow proton shuttling. On the other hand, when the pro-S oxygen is ligated to zinc (MD3), the proton on the pro-R oxygen averages 2.09 A from the hydroxyl oxygen of Ser48 such that initiation of shuttling of protons via Ser48 to the ribose 2'-hydroxyl oxygen to the 3'-hydroxyl oxygen to His51 nitrogen is sterically favorable. PM3 calculations suggest that this proton shuttle represents a stepwise reaction which occurs subsequent to hydride transfer. The PM3 transition state for hydride transfer based on the MD3 configuration has the transferring hydride 1.476 A from C4 of NAD+ and 1.433 A from the aldehyde alpha-carbon.