Summarizing lecture: factors influencing enzymatic H-transfers, analysis of nuclear tunnelling isotope effects and thermodynamic versus specific effects

Solvent and Temperature Probes of the Long-Range Electron-Transfer Step in Tyramine β-Monooxygenase: Demonstration of a Long-Range Proton-Coupled Electron-Transfer Mechanism

Tyramine β-monooxygenase (TβM) belongs to a family of physiologically important dinuclear copper monooxygenases that function with a solvent-exposed active site. To accomplish each enzymatic turnover, an electron transfer (ET) must occur between two solvent-separated copper centers. In wild-type TβM, this event is too fast to be rate limiting. However, we have recently shown [Osborne, R. L.; et al. Biochemistry2013, 52, 1179] that the Tyr216Ala variant of TβM leads to rate-limiting ET. In this study, we present a pH–rate profile study of Tyr216Ala, together with deuterium oxide solvent kinetic isotope effects (KIEs). A solvent KIE of 2 on k_cat is found in a region where k_cat is pH/pD independent. As a control, the variant Tyr216Trp, for which ET is not rate determining, displays a solvent KIE of unity. We conclude, therefore, that the observed solvent KIE arises from the rate-limiting ET step in the Tyr216Ala variant, and show how small solvent KIEs (ca. 2) can be fully accommodated from equilibrium effects within the Marcus equation. To gain insight into the role of the enzyme in the long-range ET step, a temperature dependence study was also pursued. The small enthalpic barrier of ET (E_a = 3.6 kcal/mol) implicates a significant entropic barrier, which is attributed to the requirement for extensive rearrangement of the inter-copper environment during PCET catalyzed by the Tyr216Ala variant. The data lead to the proposal of a distinct inter-domain pathway for PCET in the dinuclear copper monooxygenases.

Factors affecting hydrogen-tunneling contribution in hydroxylation reactions promoted by oxoiron(IV) porphyrin π-cation radical complexes

Hydrogen atom transfer with a tunneling effect (H-tunneling) has been proposed to be involved in aliphatic hydroxylation reactions catalyzed by cytochrome P450 and synthetic heme complexes as a result of the observation of large hydrogen/deuterium kinetic isotope effects (KIEs). In the present work, we investigate the factors controlling the H-tunneling contribution to the H-transfer process in hydroxylation reaction by examining the kinetics of hydroxylation reactions at the benzylic positions of xanthene and 1,2,3,4-tetrahydronaphthalene by oxoiron(IV) 5,10,15,20-tetramesitylporphyrin π-cation radical complexes ((TMP(+•))Fe(IV)O(L)) under single-turnover conditions. The Arrhenius plots for these hydroxylation reactions of H-isotopomers have upwardly concave profiles. The Arrhenius plots of D-isotopomers, clear isosbestic points, and product analysis rule out the participation of thermally dependent other reaction processes in the concave profiles. These results provide evidence for the involvement of H-tunneling in the rate-limiting H-transfer process. These profiles are simulated using an equation derived from Bell's tunneling model. The temperature dependence of the KIE values (k(H)/k(D)) determined for these reactions indicates that the KIE value increases as the reaction temperature becomes lower, the bond dissociation energy (BDE) of the C-H bond of a substrate becomes higher, and the reactivity of (TMP(+•))Fe(IV)O(L) decreases. In addition, we found correlation of the slope of the ln(k(H)/k(D)) - 1/T plot and the bond strengths of the Fe═O bond of (TMP(+•))Fe(IV)O(L) estimated from resonance Raman spectroscopy. These observations indicate that these factors modulate the extent of the H-tunneling contribution by modulating the ratio of the height and thickness of the reaction barrier.

Concerted versus stepwise mechanism in thymidylate synthase

Thymidylate synthase (TSase) catalyzes the intracellular de novo formation of thymidylate (a DNA building block) in most living organisms, making it a common target for chemotherapeutic and antibiotic drugs. Two mechanisms have been proposed for the rate-limiting hydride transfer step in TSase catalysis: a stepwise mechanism in which the hydride transfer precedes the cleavage of the covalent bond between the enzymatic cysteine and the product and a mechanism where both happen concertedly. Striking similarities between the enzyme-bound enolate intermediates formed in the initial and final step of the reaction supported the first mechanism, while QM/MM calculations favored the concerted mechanism. Here, we experimentally test these two possibilities using secondary kinetic isotope effect (KIE), mutagenesis study, and primary KIEs. The findings support the concerted mechanism and demonstrate the critical role of an active site arginine in substrate binding, activation of enzymatic nucleophile, and the hydride transfer studied here. The elucidation of this reduction/substitution sheds light on the critical catalytic step in TSase and may aid future drug or biomimetic catalyst design.

Dynamics and mechanism of repair of ultraviolet-induced (6-4) photoproduct by photolyase

One of the detrimental effects of UV radiation on DNA is the formation of the (6-4) photoproduct (6-4PP) between two adjacent pyrimidines1. This lesion interferes with replication and transcription and may result in mutation and cell death2. In many organisms a flavoenzyme called photolyase uses blue light energy to repair the 6-4PP3. The molecular mechanism of the repair reaction is poorly understood. Here, we use ultrafast spectroscopy to show that the key step in the repair photocycle is a cyclic proton transfer between the enzyme and the substrate. By femtosecond synchronization of the enzymatic dynamics with the repair function, we followed the function evolution and observed direct electron transfer from the excited flavin cofactor to the 6-4PP in 225 ps but surprisingly fast back electron transfer in 50 ps without repair. Strikingly, we found that the catalytic proton transfer between a histidine residue in the active site and the 6-4PP, induced by the initial photoinduced electron transfer from the excited flavin cofactor to 6-4PP, occurs in 425 ps and leads to 6-4PP repair in tens of nanoseconds. These key dynamics define the repair photocycle and explain the underlying molecular mechanism of the enzyme’s modest efficiency.

Contribution of flavin covalent linkage with histidine 99 to the reaction catalyzed by choline oxidase

The FAD-dependent choline oxidase has a flavin cofactor covalently attached to the protein via histidine 99 through an 8alpha-N(3)-histidyl linkage. The enzyme catalyzes the four-electron oxidation of choline to glycine betaine, forming betaine aldehyde as an enzyme-bound intermediate. The variant form of choline oxidase in which the histidine residue has been replaced with asparagine was used to investigate the contribution of the 8alpha-N(3)-histidyl linkage of FAD to the protein toward the reaction catalyzed by the enzyme. Decreases of 10-fold and 30-fold in the k(cat)/K(m) and k(cat) values were observed as compared with wild-type choline oxidase at pH 10 and 25 degrees C, with no significant effect on k(cat)/K(O) using choline as substrate. Both the k(cat)/K(m) and k(cat) values increased with increasing pH to limiting values at high pH consistent with the participation of an unprotonated group in the reductive half-reaction and the overall turnover of the enzyme. The pH independence of both (D)(k(cat)/K(m)) and (D)k(cat), with average values of 9.2 +/- 3.3 and 7.4 +/- 0.5, respectively, is consistent with absence of external forward and reverse commitments to catalysis, and the chemical step of CH bond cleavage being rate-limiting for both the reductive half-reaction and the overall enzyme turnover. The temperature dependence of the (D)k(red) values suggests disruption of the preorganization in the asparagine variant enzyme. Altogether, the data presented in this study are consistent with the FAD-histidyl covalent linkage being important for the optimal positioning of the hydride ion donor and acceptor in the tunneling reaction catalyzed by choline oxidase.

Hydrogen atom transfer reactions of a ruthenium imidazole complex: hydrogen tunneling and the applicability of the Marcus cross relation

The reaction of Ru(II)(acac)2(py-imH) (Ru(II)imH) with TEMPO(*) (2,2,6,6-tetramethylpiperidine-1-oxyl radical) in MeCN quantitatively gives Ru(III)(acac)2(py-im) (Ru(III)im) and the hydroxylamine TEMPO-H by transfer of H(*) (H(+) + e(-)) (acac = 2,4-pentanedionato, py-imH = 2-(2'-pyridyl)imidazole). Kinetic measurements of this reaction by UV-vis stopped-flow techniques indicate a bimolecular rate constant k(3H) = 1400 +/- 100 M(-1) s(-1) at 298 K. The reaction proceeds via a concerted hydrogen atom transfer (HAT) mechanism, as shown by ruling out the stepwise pathways of initial proton or electron transfer due to their very unfavorable thermochemistry (Delta G(o)). Deuterium transfer from Ru(II)(acac)2(py-imD) (Ru(II)imD) to TEMPO(*) is surprisingly much slower at k(3D) = 60 +/- 7 M(-1) s(-1), with k(3H)/k(3D) = 23 +/- 3 at 298 K. Temperature-dependent measurements of this deuterium kinetic isotope effect (KIE) show a large difference between the apparent activation energies, E(a3D) - E(a3H) = 1.9 +/- 0.8 kcal mol(-1). The large k(3H)/k(3D) and DeltaE(a) values appear to be greater than the semiclassical limits and thus suggest a tunneling mechanism. The self-exchange HAT reaction between Ru(II)imH and Ru(III)im, measured by (1)H NMR line broadening, occurs with k(4H) = (3.2 +/- 0.3) x 10(5) M(-1) s(-1) at 298 K and k(4H)/k(4D) = 1.5 +/- 0.2. Despite the small KIE, tunneling is suggested by the ratio of Arrhenius pre-exponential factors, log(A(4H)/A(4D)) = -0.5 +/- 0.3. These data provide a test of the applicability of the Marcus cross relation for H and D transfers, over a range of temperatures, for a reaction that involves substantial tunneling. The cross relation calculates rate constants for Ru(II)imH(D) + TEMPO(*) that are greater than those observed: k(3H,calc)/k(3H) = 31 +/- 4 and k(3D,calc)/k(3D) = 140 +/- 20 at 298 K. In these rate constants and in the activation parameters, there is a better agreement with the Marcus cross relation for H than for D transfer, despite the greater prevalence of tunneling for H. The cross relation does not explicitly include tunneling, so close agreement should not be expected. In light of these results, the strengths and weaknesses of applying the cross relation to HAT reactions are discussed.

Incorporation of hydrostatic pressure into models of hydrogen tunneling highlights a role for pressure-modulated promoting vibrations

Hydrostatic pressure offers an alternative to temperature as an experimental probe of hydrogen-transfer reactions. H tunneling reactions have been shown to exhibit kinetic isotope effects (KIEs) that are sensitive to pressure, and environmentally coupled H tunneling reactions, those reactions in which H transfer is coupled to atomic fluctuations (a promoting vibration) along the reaction coordinate, often have quite temperature-dependent KIEs. We present here a theoretical treatment of the combined effect of temperature and pressure on environmentally coupled H tunneling reactions. We develop a generalized expression for the KIE, which can be used as a simple fitting function for combined experimental temperature- and pressure-dependent KIE data sets. With this expression, we are able to extract information about the pressure dependence of both the apparent tunneling distance and the frequency of the promoting vibration. The KIE expression is tested on two data sets {the reduction of chloranil by leuco crystal violet [Isaacs, N. S., Javaid, K., and Rannala, E. (1998) J. Chem. Soc., Perkin Trans. 2, 709-711] and the reduction of morphinone reductase by NADH [Hay, S., Sutcliffe, M. J., and Scrutton, N. S. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 507-512]} and suggests that hydrostatic pressure is a sensitive probe of nuclear quantum mechanical effects in H-transfer reactions.

Contributions of the Protein Environment to the Midpoint Potentials of the A1 Phylloquinones and the FX Iron−Sulfur Cluster in Photosystem I

Electrostatic calculations have predicted that the partial negative charge associated with D575PsaB plays a significant role in modulating the midpoint potentials of the A1A and A1B phylloquinones in photosystem I. To test this prediction, the side chain of residue 575PsaB was changed from negatively charged (D) to neutral (A) and to positively charged (K). D566PsaB, which is located at a considerable distance from either A1A or A1B, and should affect primarily the midpoint potential of FX, was similarly changed. In the 575PsaB variants, the rate of electron transfer from A1A to FX is observed to decrease slightly according to the sequence D/A/K. In the 566PsaB variants, the opposite effect of a slight increase in the rate is observed according to the same sequence D/A/K. These results are consistent with the expectation that changing these residues will shift the midpoint potentials of nearby cofactors to more positive values and that the magnitude of this shift will depend on the distance between the cofactors and the residues being changed. Thus, the midpoint potentials of A1A and A1B should experience a larger shift than will FX in the 575PsaB variants, while FX should experience a larger shift than will either A1A or A1B in the 566PsaB variants. As a result, the driving energy for electron transfer from A1A and A1B to FX will be decreased in the former and increased in the latter. This rationalization of the changes in kinetics is compared with the results of electrostatic calculations. While the altered amino acids shift the midpoint potentials of A1A, A1B, and FX by different amounts, the difference in the shifts between A1A and FX or between A1B and FX is small so that the overall effect on the electron transfer rate between A1A and FX or between A1B and FX is predicted to be small. These conclusions are borne out by experiment.

H and other transfers in enzymes and in solution: theory and computations, a unified view. 2. Applications to experiment and computations

Equations obtained in part I for the free-energy barrier to one-step enzymatic reactions between bound reactants are discussed. The rate is expressed in terms of lambdao (protein reorganization energy), DeltaG(o) (standard free energy of reaction of the H-transfer step), bond breaking/bond forming term, w (work terms), and H-transmission property. Two alternative approximations for the coupling of the bond breaking/bond forming and protein are distinguished experimentally in favorable cases by the DeltaG(o) where the maximum deuterium kinetic isotope effect occurs. Plots of log rate versus DeltaG(o) and properties such as DeltaS* and DeltaS(o) are discussed. The weak or zero T-dependence of the kinetic isotope effect for wild-type enzymes operating under physiological conditions is interpreted in terms of vanishing (or isotopically insensitive) w plus transfer from the lowest H-state. Static and dynamic protein flexibility is discussed. While the many correlations accessible for electron transfers are not available for H-transfers in enzymes, a combination of experiment, computation, and analytical approaches can assist in evaluating the utility of the present equations and in suggesting further experiments and computations. A protein reorganization energy lambdao is obtained in the literature from the extended valence bond formalism where diabatic electronic states are used. A method is suggested for extracting it when instead a bond distance difference coordinate is used. The results may provide a bridge between the two approaches.

Enzymatic catalysis and transfers in solution. I. Theory and computations, a unified view

The transfer of hydride, proton, or H atom between substrate and cofactor in enzymes has been extensively studied for many systems, both experimentally and computationally. A simple equation for the reaction rate, an analog of an equation obtained earlier for electron transfer rates, is obtained, but now containing an approximate analytic expression for the bond rupture-bond forming feature of these H transfers. A "symmetrization," of the potential energy surfaces is again introduced [R. A. Marcus, J. Chem. Phys. 43, 679 (1965); J. Phys. Chem. 72, 891 (1968)], together with Gaussian fluctuations of the remaining coordinates of the enzyme and solution needed for reaching the transition state. Combining the two expressions for the changes in the difference of the two bond lengths of the substrate-cofactor subsystem and in the fluctuation coordinates of the protein leading to the transition state, an expression is obtained for the free energy barrier. To this end a two-dimensional reaction space (m,n) is used that contains the relative coordinates of the H in the reactants, the heavy atoms to which it is bonded, and the protein/solution reorganization coordinate, all leading to the transition state. The resulting expression may serve to characterize in terms of specific parameters (two "reorganization" terms, thermodynamics, and work terms), experimental and computational data for different enzymes, and different cofactor-substrate systems. A related characterization was used for electron transfers. To isolate these factors from nuclear tunneling, when the H-tunneling effect is large, use of deuterium and tritium transfers is of course helpful, although tunneling has frequently and understandably dominated the discussions. A functional form is suggested for the dependence of the deuterium kinetic isotope effect (KIE) on DeltaG degrees and a different form for the 13C KIE. Pressure effects on deuterium and 13C KIEs are also discussed. Although formulated for a one-step transfer of a light particle in an enzyme, the results would also apply to single-step transfers of other atoms and groups in enzymes and in solution.