Daniel Roston

Elusive transition state of alcohol dehydrogenase unveiled

For several decades the hydride transfer catalyzed by alcohol dehydrogenase has been difficult to understand. Here we add to the large corpus of anomalous and paradoxical data collected for this reaction by measuring a normal (> 1) 2° kinetic isotope effect (KIE) for the reduction of benzaldehyde. Because the relevant equilibrium effect is inverse (< 1), this KIE eludes the traditional interpretation of 2° KIEs. It does, however, enable the development of a comprehensive model for the “tunneling ready state” (TRS) of the reaction that fits into the general scheme of Marcus-like models of hydrogen tunneling. The TRS is the ensemble of states along the intricate reorganization coordinate, where H tunneling between the donor and acceptor occurs (the crossing point in Marcus theory). It is comparable to the effective transition state implied by ensemble-averaged variational transition state theory. Properties of the TRS are approximated as an average of the individual properties of the donor and acceptor states. The model is consistent with experimental findings that previously appeared contradictory; specifically, it resolves the long-standing ambiguity regarding the location of the TRS (aldehyde-like vs. alcohol-like). The new picture of the TRS for this reaction identifies the principal components of the collective reaction coordinate and the average structure of the saddle point along that coordinate.

Molecular simulation of water and hydration effects in different environments: challenges and developments for DFTB based models.

We discuss the description of water and hydration effects that employs an approximate density functional theory, DFTB3, in either a full QM or QM/MM framework. The goal is to explore, with the current formulation of DFTB3, the performance of this method for treating water in different chemical environments, the magnitude and nature of changes required to improve its performance, and factors that dictate its applicability to reactions in the condensed phase in a QM/MM framework. A relatively minor change (on the scale of kBT) in the O–H repulsive potential is observed to substantially improve the structural properties of bulk water under ambient conditions; modest improvements are also seen in dynamic properties of bulk water. This simple change also improves the description of protonated water clusters, a solvated proton, and to a more limited degree, a solvated hydroxide. By comparing results from DFTB3 models that differ in the description of water, we confirm that proton transfer energetics are adequately described by the standard DFTB3/3OB model for meaningful mechanistic analyses. For QM/MM applications, a robust parametrization of QM-MM interactions requires an explicit consideration of condensed phase properties, for which an efficient sampling technique was developed recently and is reviewed here. The discussions help make clear the value and limitations of DFTB3 based simulations, as well as the developments needed to further improve the accuracy and transferability of the methodology.

Isotope Effects as Probes for Enzyme Catalyzed Hydrogen-Transfer Reactions

Kinetic Isotope effects (KIEs) have long served as a probe for the mechanisms of both enzymatic and solution reactions. Here, we discuss various models for the physical sources of KIEs, how experimentalists can use those models to interpret their data, and how the focus of traditional models has grown to a model that includes motion of the enzyme and quantum mechanical nuclear tunneling. We then present two case studies of enzymes, thymidylate synthase and alcohol dehydrogenase, and discuss how KIEs have shed light on the C-H bond cleavages those enzymes catalyze. We will show how the combination of both experimental and computational studies has changed our notion of how these enzymes exert their catalytic powers.

Kinetic Isotope Effects as a Probe of Hydrogen Transfers to and from Common Enzymatic Cofactors

Enzymes use a number of common cofactors as sources of hydrogen to drive biological processes, but the physics of the hydrogen transfers to and from these cofactors is not fully understood. Researchers study the mechanistically important contributions from quantum tunneling and enzyme dynamics and connect those processes to the catalytic power of enzymes that use these cofactors. Here we describe some progress that has been made in studying these reactions, particularly through the use of kinetic isotope effects (KIEs). We first discuss the general theoretical framework necessary to interpret experimental KIEs, and then describe practical uses for KIEs in the context of two case studies. The first example is alcohol dehydrogenase, which uses a nicotinamide cofactor to catalyze a hydride transfer, and the second example is thymidylate synthase, which uses a folate cofactor to catalyze both a hydride and a proton transfer.

Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site

A reactions transition state (TS) structure plays a critical role in determining reactivity and has important implications for the design of catalysts, drugs, and other applications. Here, we explore TS structure in the enzyme alkaline phosphatase using hybrid Quantum Mechanics/Molecular Mechanics simulations. We find that minor perturbations to the substrate have major effects on TS structure and the way the enzyme stabilizes the TS. Substrates with good leaving groups (LGs) have little cleavage of the phosphorus-LG bond at the TS, while substrates with poor LGs have substantial cleavage of that bond. The results predict nonlinear free energy relationships for a single rate-determining step, and substantial differences in kinetic isotope effects for different substrates; both trends were observed in previous experimental studies, although the original interpretations differed from the present model. Moreover, due to different degrees of phosphorus-LG bond cleavage at the TS for different substrates, the LG is stabilized by different interactions at the TS: while a poor LG is directly stabilized by an active site zinc ion, a good LG is mainly stabilized by active site water molecules. Our results demonstrate the considerable plasticity of TS structure and stabilization in enzymes. Furthermore, perturbations to reactivity that probe TS structure experimentally (i.e., substituent effects) may substantially perturb the TS they aim to probe, and thus classical experimental approaches such as free energy relations should be interpreted with care.

A Critical Test of the “Tunneling and Coupled Motion” Concept in Enzymatic Alcohol Oxidation

The physical mechanism of C-H bond activation by enzymes is the subject of intense study, and we have tested the predictions of two competing models for C-H activation in the context of alcohol dehydrogenase. The kinetic isotope effects (KIEs) in this enzyme have previously suggested a model of quantum mechanical tunneling and coupled motion of primary (1°) and secondary (2°) hydrogens. Here we measure the 2° H/T KIEs with both H and D at the 1° position and find that the 2° KIE is significantly deflated with D-transfer, consistent with the predictions of recent Marcus-like models of H-transfer. The results suggest that the fast dynamics of H-tunneling result in a 1° isotope effect on the structure of the tunneling ready state: the trajectory of D-transfer goes through a shorter donor-acceptor distance than that of H-transfer.

Simulations of remote mutants of dihydrofolate reductase reveal the nature of a network of residues coupled to hydride transfer

Recent experimental and theoretical studies have proposed that enzymes involve networks of coupled residues throughout the protein that participate in motions accompanying chemical barrier crossing. Here, we have examined portions of a proposed network in dihydrofolate reductase (DHFR) using quantum mechanics/molecular mechanics simulations. The simulations use a hybrid quantum mechanics‐molecular mechanics approach with a recently developed semiempirical AM1‐SRP Hamiltonian that provides accurate results for this reaction. The simulations reproduce experimentally determined catalytic rates for the wild type and distant mutants of E. coli DHFR, underscoring the accuracy of the simulation protocol. Additionally, the simulations provide detailed insight into how residues remote from the active site affect the catalyzed chemistry, through changes in the thermally averaged properties along the reaction coordinate. The mutations do not greatly affect the structure of the transition state near the bond activation, but we observe differences somewhat removed from the point of CH cleavage that affect the rate. The mutations have global effects on the thermally averaged structure that propagate throughout the enzyme and the current simulations highlight several interactions that appear to be particularly important.

Experimental and Theoretical Studies of Enzyme-Catalyzed Hydrogen-Transfer Reactions

The mechanisms of enzyme-catalyzed reactions are medicinally important and present a fascinating intellectual challenge. Many experimental and theoretical techniques can shed light on these mechanisms, and here, we shall focus on the utility of kinetic isotope effects (KIEs) to study enzymatic reactions that involve hydrogen transfers. We will provide a short background on the prevailing models to interpret KIEs and then present more detailed reviews of two model enzymes: alcohol dehydrogenase and thymidylate synthase. These two examples provide a context to describe the types of experiments and theoretical calculations that drive this field forward and the kind of information each can furnish. We emphasize the importance of cooperation between experimentalists and theoreticians to continue the progress toward a comprehensive theory of enzyme catalysis.

Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin

The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis. Comparison of the observed hydrolysis pathways and the corresponding free energy profiles leads to detailed models for the mechanism of ATP hydrolysis in the pre-powerstroke state and proposes factors that regulate the hydrolysis activity in different conformational states. In the pre-powerstroke state, the scissile Pγ-O3β bond breaks early in the reaction. Proton transfer from the lytic water to the γ-phosphate through active site residues is an important part of the kinetic bottleneck; several hydrolysis pathways that feature distinct proton transfer routes are found to have similar free energy barriers, suggesting a significant degree of plasticity in the hydrolysis mechanism. Comparison of hydrolysis in the pre-powerstroke state and the closed post-rigor model suggests that optimization of residues beyond the active site for electrostatic stabilization and preorganization is likely important to enzyme design.

Computational replication of the abnormal secondary kinetic isotope effects in a hydride transfer reaction in solution with a motion assisted H-tunneling model.

We recently reported abnormal secondary deuterium kinetic isotope effects (2° KIEs) for hydride transfer reactions from alcohols to carbocations in acetonitrile (Chem. Comm. 2012, 48, 11337). Experimental 2° KIE values were found to be inflated on the 9-C position in the xanthylium cation but deflated on the β-C position in 2-propanol with respect to the values predicted by the semi-classical transition-state theory. No primary (1°) isotope effect on 2° KIEs was observed. Herein, the KIEs were replicated by the Marcus-like H-tunneling model that requires a longer donor–acceptor distance (DAD) in a lighter isotope transfer process. The 2° KIEs for a range of potential tunneling-ready-states (TRSs) of different DADs were calculated and fitted to the experiments to find the TRS structure. The observed no effect of 1° isotope on 2° KIEs is explained in terms of the less sterically hindered TRS structure so that the change in DAD due to the change in 1° isotope does not significantly affect the reorganization of the 2° isotope and hence the 2° KIE. The effect of 1° isotope on 2° KIEs may be expected to be more pronounced and thus observable in reactions occurring in restrictive environments such as the crowded and relatively rigid active site of enzymes.