Enzymology takes a quantum leap forward

Utilizing Molecular Dynamics Simulations, Machine Learning, Cryo-EM, and NMR Spectroscopy to Predict and Validate Protein Dynamics

Protein dynamics play a crucial role in biological function, encompassing motions ranging from atomic vibrations to large-scale conformational changes. Recent advancements in experimental techniques, computational methods, and artificial intelligence have revolutionized our understanding of protein dynamics. Nuclear magnetic resonance spectroscopy provides atomic-resolution insights, while molecular dynamics simulations offer detailed trajectories of protein motions. Computational methods applied to X-ray crystallography and cryo-electron microscopy (cryo-EM) have enabled the exploration of protein dynamics, capturing conformational ensembles that were previously unattainable. The integration of machine learning, exemplified by AlphaFold2, has accelerated structure prediction and dynamics analysis. These approaches have revealed the importance of protein dynamics in allosteric regulation, enzyme catalysis, and intrinsically disordered proteins. The shift towards ensemble representations of protein structures and the application of single-molecule techniques have further enhanced our ability to capture the dynamic nature of proteins. Understanding protein dynamics is essential for elucidating biological mechanisms, designing drugs, and developing novel biocatalysts, marking a significant paradigm shift in structural biology and drug discovery.

Interplay of structural preorganization and conformational sampling in UDP-glucuronic acid 4-epimerase catalysis

Understanding enzyme catalysis as connected to protein motions is a major challenge. Here, based on temperature kinetic studies combined with isotope effect measurements, we obtain energetic description of C-H activation in NAD-dependent UDP-glucuronic acid C4 epimerase. Approach from the ensemble-averaged ground state (GS) to the transition state-like reactive conformation (TSRC) involves, alongside uptake of heat (Δ H ‡  = 54 kJ mol-1), significant loss in entropy ( - T Δ S ‡  = 20 kJ mol-1; 298 K) and negative activation heat capacity (Δ C p ‡  = -0.64 kJ mol-1 K-1). Thermodynamic changes suggest the requirement for restricting configurational freedom at the GS to populate the TSRC. Enzyme variants affecting the electrostatic GS preorganization reveal active-site interactions important for precise TSRC sampling and H-transfer. Collectively, our study captures thermodynamic effects associated with TSRC sampling and establishes rigid positioning for C-H activation in an enzyme active site that requires conformational flexibility in fulfillment of its natural epimerase function.

Converting genetic information to non-equilibrium cellular thermodynamics in enzyme-catalyzed reactions

Living systems use genomic information to maintain a stable highly ordered state far from thermodynamic equilibrium but the specific mechanisms and general principles governing the interface of genetics and thermodynamics has not been extensively investigated. Genetic information is quantified in unitless bits termed "Shannon entropy", which does not directly relate to thermodynamic entropy or energy. Thus, it is unclear how the Shannon entropy of genetic information is converted into thermodynamic work necessary to maintain the non-equilibrium state of living systems. Here we investigate the interface of genetic information and cellular thermodynamics in enzymatic acceleration of a chemical reaction S + E → E S → E + P , where S and E are substrate and enzyme, E S is the enzyme substrate complex and P product. The rate of any intracellular chemical reaction is determined by probability functions at macroscopic (Boltzmann distribution of the reactant kinetic energies governed by temperature) or microscopic (overlap of reactant quantum wave functions) scales - described, respectively, by the Arrhenius and Knudsen equations. That is, the reaction rate, in the absence of a catalyst, is governed by temperature which determines the kinetic energy of the interacting molecules. Genetic information can act upon a when the encoded string of amino acids folds into a 3-deminsional structure that permits a lock/key spatial matching with the reactants. By optimally superposing the reactants' wave functions, the information in the enzyme increases the reaction rate by up to15 orders of magnitude under isothermal conditions. In turn, the accelerated reaction rate alters the intracellular thermodynamics environment as the products are at lower Gibbs free energy which permits thermodynamic work W m a x = - Δ G . Mathematically and biologically, the critical event that allows genetic information to produce thermodynamic work is the folding of the amino acid string specified by the gene into a 3-dimensional shape determined by its lowest energy state. Biologically, this allows the amino acid string to bind substrate and place them in an optimal spatial orientation. These key-lock are mathematically characterized by Kullback-Leibler Divergence and the interactions with the reaction channel now represent Fisher Information (the second derivative Kullback-Leibler divergence), which can take on the units of the process to which it is applied. Interestingly, Shannon is typically derived by "coarse graining" Shannon information. Thus, living system, by acting at a quantum level, "fine grain" Shannon information.

The Constrained Disorder Principle Accounts for the Variability That Characterizes Breathing: A Method for Treating Chronic Respiratory Diseases and Improving Mechanical Ventilation

Variability characterizes breathing, cellular respiration, and the underlying quantum effects. Variability serves as a mechanism for coping with changing environments; however, this hypothesis does not explain why many of the variable phenomena of respiration manifest randomness. According to the constrained disorder principle (CDP), living organisms are defined by their inherent disorder bounded by variable boundaries. The present paper describes the mechanisms of breathing and cellular respiration, focusing on their inherent variability. It defines how the CDP accounts for the variability and randomness in breathing and respiration. It also provides a scheme for the potential role of respiration variability in the energy balance in biological systems. The paper describes the option of using CDP-based artificial intelligence platforms to augment the respiratory process's efficiency, correct malfunctions, and treat disorders associated with the respiratory system.

Cardiac Contractility Modulation Therapy in Patients with Amyloid Cardiomyopathy and Heart Failure, Case Report, Review of the Biophysics of CCM Function, and AMY-CCM Registry Presentation

Cardiac amyloidosis may result in an aggressive form of heart failure (HF). Cardiac contractility modulation (CCM) has been shown to be a concrete therapeutic option in patients with symptomatic HF, but there is no evidence of its application in patients with cardiac amyloidosis. We present the case of TTR amyloidosis, where CCM therapy proved to be effective. The patient had a history of multiple HF hospitalizations due to an established diagnosis of wild type TTR-Amyloidosis with significant cardiac involvement. Since he was highly symptomatic, except during continuous dobutamine and diuretic infusion, it was opted to pursue CCM therapy device implantation. At follow up, a significant improvement in clinical status was reported with an increase of EF, functional status (6 min walk test improved from zero meters at baseline, to 270 m at 1 month and to 460 m at 12 months), and a reduction in pulmonary pressures. One year after device implantation, no other HF hospital admission was needed. CCM therapy may be effective in this difficult clinical setting. The AMY-CCM Registry, which has just begun, will evaluate the efficacy of CCM in patients with HF and diagnosed TTR amyloidosis to bring new evidence on its potential impact as a therapeutic option.

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car-Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers.

Quantum tunnelling in the context of SARS-CoV-2 infection

The SARS-CoV-2 pandemic has added new urgency to the study of viral mechanisms of infection. But while vaccines offer a measure of protection against this specific outbreak, a new era of pandemics has been predicted. In addition to this, COVID-19 has drawn attention to post-viral syndromes and the healthcare burden they entail. It seems integral that knowledge of viral mechanisms is increased through as wide a research field as possible. To this end we propose that quantum biology might offer essential new insights into the problem, especially with regards to the important first step of virus-host invasion. Research in quantum biology often centres around energy or charge transfer. While this is predominantly in the context of photosynthesis there has also been some suggestion that cellular receptors such as olfactory or neural receptors might employ vibration assisted electron tunnelling to augment the lock-and-key mechanism. Quantum tunnelling has also been observed in enzyme function. Enzymes are implicated in the invasion of host cells by the SARS-CoV-2 virus. Receptors such as olfactory receptors also appear to be disrupted by COVID-19. Building on these observations we investigate the evidence that quantum tunnelling might be important in the context of infection with SARS-CoV-2. We illustrate this with a simple model relating the vibronic mode of, for example, a viral spike protein to the likelihood of charge transfer in an idealised receptor. Our results show a distinct parameter regime in which the vibronic mode of the spike protein enhances electron transfer. With this in mind, novel therapeutics to prevent SARS-CoV-2 transmission could potentially be identified by their vibrational spectra.

De novo protein folding on computers. Benefits and challenges

There has been recent success in prediction of the three-dimensional folded native structures of proteins, most famously by the AlphaFold Algorithm running on Google's/Alphabet's DeepMind computer. However, this largely involves machine learning of protein structures and is not a de novo protein structure prediction method for predicting three-dimensional structures from amino acid residue sequences. A de novo approach would be based almost entirely on general principles of energy and entropy that govern protein folding energetics, and importantly do so without the use of the amino acid sequences and structural features of other proteins. Most consider that problem as still unsolved even though it has occupied leading scientists for decades. Many consider that it remains one of the major outstanding issues in modern science. There is crucial continuing help from experimental findings on protein unfolding and refolding in the laboratory, but only to a limited extent because many researchers consider that the speed by which real proteins folds themselves, often from milliseconds to minutes, is itself still not fully understood. This is unfortunate, because a practical solution to the problem would probably have a major effect on personalized medicine, the pharmaceutical industry, biotechnology, and nanotechnology, including for example "smaller" tasks such as better modeling of flexible "unfolded" regions of the SARS-COV-2 spike glycoprotein when interacting with its cell receptor, antibodies, and therapeutic agents. Some important ideas from earlier studies are given before moving on to lessons from periodic and aperiodic crystals, and a possible role for quantum phenomena. The conclusion is that better computation of entropy should be the priority, though that is presented guardedly.

Perspective: Defining and quantifying the role of dynamics in enzyme catalysis

Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme's conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.

Evidence that a 'dynamic knockout' in Escherichia coli dihydrofolate reductase does not affect the chemical step of catalysis

The question of whether protein motions play a role in the chemical step of enzymatic catalysis has generated much controversy in recent years. Debate has recently reignited over possible dynamic contributions to catalysis in dihydrofolate reductase, following conflicting conclusions from studies of the N23PP/S148A variant of the Escherichia coli enzyme. By investigating the temperature dependence of kinetic isotope effects, we present evidence that the reduction in the hydride transfer rate constants in this variant is not a direct result of impairment of conformational fluctuations. Instead, the conformational state of the enzyme immediately before hydride transfer, which determines the electrostatic environment of the active site, affects the rate constant for the reaction. Although protein motions are clearly important for binding and release of substrates and products, there appears to be no detectable dynamic coupling of protein motions to the hydride transfer step itself.

At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis?

Enzymes play a key role in almost all biological processes, accelerating a variety of metabolic reactions as well as controlling energy transduction, the transcription, and translation of genetic information, and signaling. They possess the remarkable capacity to accelerate reactions by many orders of magnitude compared to their uncatalyzed counterparts, making feasible crucial processes that would otherwise not occur on biologically relevant timescales. Thus, there is broad interest in understanding the catalytic power of enzymes on a molecular level. Several proposals have been put forward to try to explain this phenomenon, and one that has rapidly gained momentum in recent years is the idea that enzyme dynamics somehow contributes to catalysis. This review examines the dynamical proposal in a critical way, considering basically all reasonable definitions, including (but not limited to) such proposed effects as "coupling between conformational and chemical motions," "landscape searches" and "entropy funnels." It is shown that none of these proposed effects have been experimentally demonstrated to contribute to catalysis, nor are they supported by consistent theoretical studies. On the other hand, it is clarified that careful simulation studies have excluded most (if not all) dynamical proposals. This review places significant emphasis on clarifying the role of logical definitions of different catalytic proposals, and on the need for a clear formulation in terms of the assumed potential surface and reaction coordinate. Finally, it is pointed out that electrostatic preorganization actually accounts for the observed catalytic effects of enzymes, through the corresponding changes in the activation free energies.

Analysis of Classical and Quantum Paths for Deprotonation of Methylamine by Methylamine Dehydrogenase

The hydrogen-transfer reaction catalysed by methylamine dehydrogenase (MADH) with methylamine (MA) as substrate is a good model system for studies of proton tunnelling in enzyme reactions--an area of great current interest--for which atomistic simulations will be vital. Here, we present a detailed analysis of the key deprotonation step of the MADH/MA reaction and compare the results with experimental observations. Moreover, we compare this reaction with the related aromatic amine dehydrogenase (AADH) reaction with tryptamine, recently studied by us, and identify possible causes for the differences observed in the measured kinetic isotope effects (KIEs) of the two systems. We have used combined quantum mechanics/molecular mechanics (QM/MM) techniques in molecular dynamics simulations and variational transition state theory with multidimensional tunnelling calculations averaged over an ensemble of paths. The results reveal important mechanistic complexity. We calculate activation barriers and KIEs for the two possible proton transfers identified-to either of the carboxylate oxygen atoms of the catalytic base (Asp428beta)-and analyse the contributions of quantum effects. The activation barriers and tunnelling contributions for the two possible proton transfers are similar and lead to a phenomenological activation free energy of 16.5+/-0.9 kcal mol(-1) for transfer to either oxygen (PM3-CHARMM calculations applying PM3-SRP specific reaction parameters), in good agreement with the experimental value of 14.4 kcal mol(-1). In contrast, for the AADH system, transfer to the equivalent OD1 was found to be preferred. The structures of the enzyme complexes during reaction are analysed in detail. The hydrogen bond of Thr474beta(MADH)/Thr172beta(AADH) to the catalytic carboxylate group and the nonconserved active site residue Tyr471beta(MADH)/Phe169beta(AADH) are identified as important factors in determining the preferred oxygen acceptor. The protein environment has a significant effect on the reaction energetics and hence on tunnelling contributions and KIEs. These environmental effects, and the related clearly different preferences for the two carboxylate oxygen atoms (with different KIEs) in MADH/MA and AADH/tryptamine, are possible causes of the differences observed in the KIEs between these two important enzyme reactions.

Origin of the Temperature Dependence of Isotope Effects in Enzymatic Reactions: The Case of Dihydrofolate Reductase

The origin of the temperature dependence of kinetic isotope effects (KIEs) in enzyme reactions is a problem of general interest and a major challenge for computational chemistry. The present work simulates the nuclear quantum mechanical (NQM) effects and the corresponding KIE in dihydrofolate reductase (DHFR) and two of its mutants by using the empirical valence bond (EVB) and the quantum classical path (QCP) centroid path integral approach. Our simulations reproduce the overall observed trend while using a fully microscopic rather than a phenomenological picture and provide an interesting insight. It appears that the KIE increases when the distance between the donor and acceptor increases, in a somewhat counter intuitive way. The temperature dependence of the KIE appears to reflect mainly the temperature dependence of the distance between the donor and acceptor. This trend is also obtained from a simplified vibronic treatment, but as demonstrated here, the vibronic treatment is not valid at short and medium distances, where it is essential to use the path integral or other approaches capable of moving seamlessly from the adiabatic to the diabatic limits. It is pointed out that although the NQM effects do not contribute to catalysis in DHFR, the observed temperature dependence can be used to refine the potential of mean force for the donor and acceptor distance and its change due to distanced mutations.

Variation of atomic charges on proton transfer in strong hydrogen bonds: The case of anionic and neutral imidazole–acetate complexes

The variation of atomic charges upon proton transfer in hydrogen bonding complexes of 4-methylimidazole, in both neutral and protonated cationic forms, and acetate anion, is investigated. These complexes model the histidine (neutral and protonated)-aspartate pair present in active sites of proteases where strong N--H...O hydrogen bonds are formed. Three procedures (Merz-Kollman scheme, Natural Population Analysis, and Atoms in Molecules Method) are used to compute atomic charges and explore their variation upon H-transfer in the gas phase and in the presence of two continuum media with dielectric constants 5 (protein interiors) and 78.39 (water). The effect of electron correlation was also studied by comparing Hartree-Fock and MP2 results for both complexes in the gas phase. Greater net charge interchanged upon H-transfer is observed in the anionic complex with respect to the neutral complex. Raising the polarity of the medium increases the amount of net charge transfer in both complexes, although the neutral system exhibits a larger sensitivity to the presence of solvent. Charge transfer associated to N--H...O and N...H--O bonds reveal the ionic contribution to the interaction depending on the number of charged subunits but the presence of solvent affects little this quantity. The lack of electron correlation overestimates all the charges as well as their variations and so uncorrelated calculations should be avoided.

Importance of quantum tunneling in vacancy-hydrogen complexes in diamond

Our ab initio calculations of the hyperfine parameters for negatively charged vacancy-hydrogen and nitrogen-vacancy-hydrogen complexes in diamond compare static defect models and models which account for the quantum tunneling behavior of hydrogen. The static models give rise to hyperfine splittings that are inconsistent with the experimental electron paramagnetic resonance data. In contrast, the hyperfine parameters for the quantum dynamical models are in agreement with the experimental observations. We show that the quantum motion of the proton is crucial to the prediction of symmetry and hyperfine constants for two simple defect centers in diamond. Static a priori methods fail for these systems.

A new conceptual framework for enzyme catalysis. Hydrogen tunnelling coupled to enzyme dynamics in flavoprotein and quinoprotein enzymes

Recent years have witnessed high levels of activity in identifying enzyme systems that catalyse H-transfer by quantum tunneling. Rather than being restricted to a small number of specific enzymes as perceived initially, it has now become an accepted mechanism for H-transfer in a growing number of enzymes. Furthermore, H-tunneling is driven by the thermally induced dynamics of the enzyme. In some of those enzymes that break stable C-H bonds the reaction proceeds purely by quantum tunneling, without the need to partially ascend the barrier. Enzymes studied that fall into this category include the flavoprotein and quinoprotein amine dehydrogenases, which have proved to be excellent model systems. These enzymes have enabled us to study the relationship between barrier shape and reaction kinetics. This has involved studies with "slow" and "fast" substrates and enzymes impaired by mutagenesis. A number of key questions now remain, including the nature of the coupling between protein dynamics and quantum tunneling. The wide-ranging implications of quantum tunneling introduce a paradigm shift in the conceptual framework for enzyme catalysis, inhibition and design.

Deuterium isotope effects during carbon-hydrogen bond cleavage by trimethylamine dehydrogenase. Implications for mechanism and vibrationally assisted hydrogen tunneling in wild-type and mutant enzymes

His-172 and Tyr-169 are components of a triad in the active site of trimethylamine dehydrogenase (TMADH) comprising Asp-267, His-172, and Tyr-169. Stopped-flow kinetic studies with trimethylamine as substrate have indicated that mutation of His-172 to Gln reduces the limiting rate constant for flavin reduction approximately 10-fold (Basran, J., Sutcliffe, M. J., Hille, R., and Scrutton, N. S. (1999) Biochem. J. 341, 307-314). A kinetic isotope effect (KIE = k(H)/k(D)) accompanies flavin reduction by H172Q TMADH, the magnitude of which varies significantly with solution pH. With trimethylamine, flavin reduction by H172Q TMADH is controlled by a single macroscopic ionization (pK(a) = 6.8 +/- 0.1). This ionization is perturbed (pK(a) = 7.4 +/- 0.1) in reactions with perdeuterated trimethylamine and is responsible for the apparent variation in the KIE with solution pH. At pH 9.5, where the functional group controlling flavin reduction is fully ionized, the KIE is independent of temperature in the range 277-297 K, consistent with vibrationally assisted hydrogen tunneling during breakage of the substrate C-H bond. Y169F TMADH is approximately 4-fold more compromised than H172Q TMADH for hydrogen transfer, which occurs non-classically. Studies with Y169F TMADH suggest partial thermal excitation of substrate prior to hydrogen tunneling by a vibrationally assisted mechanism. Our studies illustrate the varied effects of compromising mutations on tunneling regimes in enzyme molecules.

Redox-dependent conformational selection in a Cys4Fe2S2 ferredoxin

Putidaredoxin (Pdx), a Cys4Fe2S2 ferredoxin from Pseudomonas putida, exhibits redox-dependent binding to its physiological redox partner, cytochrome P450(cam) (CYP101), with the reduced form of Pdx (Pdx(r)) binding with greater affinity to oxidized camphor-bound CYP101 than the oxidized form, Pdx(o). It has been previously shown that Pdx(o) is more dynamic than Pdx(r) on all accessible time scales, and it has been proposed that Pdx(r) samples only a fraction of the conformational substates populated by Pdx(o) on a time average. It is postulated that the ensemble subset populated by Pdx(r) is the same subset that binds CYP101, providing a mechanism for coupling the Pdx oxidation state to binding affinity for CYP101. Evidence from a variety of sources, including redox-dependent shifts of 15N and 13C resonances, indicates that the metal cluster binding loop of Pdx is the primary determinant of redox-dependent conformational selection. Patterns of paramagnetic effects suggest that the metal cluster binding loop contracts around the metal cluster upon reduction, possibly due to the strengthening of hydrogen bonds between the sulfur atoms of the metal cluster and the surrounding polypeptide NH and OH groups. Effects of this perturbation are then transmitted mechanically to other affected regions of the protein. A specific mutation has been introduced into the metal binding loop of Pdx, G40N, that slows conformational exchange sufficiently that the ensemble of conformational substates in Pdx(o) are directly observable as severe broadenings or splittings in affected NMR resonances. Many of the residues most affected by the mutation also show significant exchange contributions to 15N T(2) relaxation in wild-type Pdx(o). As predicted, G40N Pdx(r) shows a collapse of many of these multiplets and broadened lines to form much sharper resonances that are essentially identical to those observed in wild-type Pdx(r), indicating that Pdx(r) occupies fewer conformational substates than does Pdx(o). This is the first direct observation of such redox-dependent ensembles at slow exchange on the chemical shift time scale. These results confirm that conformational selection within the Fe2S2 cluster binding loop is the primary source of redox-dependent changes in protein dynamics in Pdx.

Toward a quantum-mechanical description of metal-assisted phosphoryl transfer in pyrophosphatase

The wealth of kinetic and structural information makes inorganic pyrophosphatases (PPases) a good model system to study the details of enzymatic phosphoryl transfer. The enzyme accelerates metal-complexed phosphoryl transfer 10(10)-fold: but how? Our structures of the yeast PPase product complex at 1.15 A and fluoride-inhibited complex at 1.9 A visualize the active site in three different states: substrate-bound, immediate product bound, and relaxed product bound. These span the steps around chemical catalysis and provide strong evidence that a water molecule (O(nu)) directly attacks PPi with a pK(a) vastly lowered by coordination to two metal ions and D117. They also suggest that a low-barrier hydrogen bond (LBHB) forms between D117 and O(nu), in part because of steric crowding by W100 and N116. Direct visualization of the double bonds on the phosphates appears possible. The flexible side chains at the top of the active site absorb the motion involved in the reaction, which may help accelerate catalysis. Relaxation of the product allows a new nucleophile to be generated and creates symmetry in the elementary catalytic steps on the enzyme. We are thus moving closer to understanding phosphoryl transfer in PPases at the quantum mechanical level. Ultra-high resolution structures can thus tease out overlapping complexes and so are as relevant to discussion of enzyme mechanism as structures produced by time-resolved crystallography.

Importance of barrier shape in enzyme-catalyzed reactions. Vibrationally assisted hydrogen tunneling in tryptophan tryptophylquinone-dependent amine dehydrogenases

C-H bond breakage by tryptophan tryptophylquinone (TTQ)-dependent methylamine dehydrogenase (MADH) occurs by vibrationally assisted tunneling (Basran, J., Sutcliffe, M. J., and Scrutton, N. S. (1999) Biochemistry 38, 3218--3222). We show here a similar mechanism in TTQ-dependent aromatic amine dehydrogenase (AADH). The rate of TTQ reduction by dopamine in AADH has a large, temperature independent kinetic isotope effect (KIE = 12.9 +/- 0.2), which is highly suggestive of vibrationally assisted tunneling. H-transfer is compromised with benzylamine as substrate and the KIE is deflated (4.8 +/- 0.2). The KIE is temperature-independent, but reaction rates are strongly dependent on temperature. With tryptamine as substrate reaction rates can be determined only at low temperature as C-H bond cleavage is rapid, and an exceptionally large KIE (54.7 +/- 1.0) is observed. Studies with deuterated tryptamine suggest vibrationally assisted tunneling is the mechanism of deuterium and, by inference, hydrogen transfer. Bond cleavage by MADH using a slow substrate (ethanolamine) occurs with an inflated KIE (14.7 +/- 0.2 at 25 degrees C). The KIE is temperature-dependent, consistent with differential tunneling of protium and deuterium. Our observations illustrate the different modes of H-transfer in MADH and AADH with fast and slow substrates and highlight the importance of barrier shape in determining reaction rate.