We present an atomic-level description of the reaction chemistry of an enzyme-catalyzed reaction dominated by proton tunneling. By solving structures of reaction intermediates at near-atomic resolution, we have identified the reaction pathway for tryptamine oxidation by aromatic amine dehydrogenase. Combining experiment and computer simulation, we show proton transfer occurs predominantly to oxygen O2 of Asp ¹²⁸ β in a reaction dominated by tunneling over ∼0.6 angstroms. The role of long-range coupled motions in promoting tunneling is controversial. We show that, in this enzyme system, tunneling is promoted by a short-range motion modulating proton-acceptor distance and no long-range coupled motion is required.

Congenital hypogonadotropic hypogonadism (CHH) is a condition characterized by absent puberty and infertility due to gonadotropin releasing hormone (GnRH) deficiency, which is often associated with anosmia (Kallmann syndrome, KS). We identified loss-of-function heterozygous mutations in anti-Müllerian hormone (AMH) and its receptor, AMHR2, in 3% of CHH probands using whole-exome sequencing. We showed that during embryonic development, AMH is expressed in migratory GnRH neurons in both mouse and human fetuses and unconvered a novel function of AMH as a pro-motility factor for GnRH neurons. Pathohistological analysis of Amhr2-deficient mice showed abnormal development of the peripheral olfactory system and defective embryonic migration of the neuroendocrine GnRH cells to the basal forebrain, which results in reduced fertility in adults. Our findings highlight a novel role for AMH in the development and function of GnRH neurons and indicate that AMH signaling insufficiency contributes to the pathogenesis of CHH in humans.

It is now widely accepted that enzyme-catalysed C–H bond breakage occurs by quantum mechanical tunnelling. This paradigm shift in the conceptual framework for these reactions away from semi-classical transition state theory (TST, i.e. including zero-point energy, but with no tunnelling correction) has been driven over the recent years by experimental studies of the temperature dependence of kinetic isotope effects (KIEs) for these reactions in a range of enzymes, including the tryptophan tryptophylquinone-dependent enzymes such as methylamine dehydrogenase and aromatic amine dehydrogenase, and the flavoenzymes such as morphinone reductase and pentaerythritol tetranitrate reductase, which produced observations that are also inconsistent with the simple Bell-correction model of tunnelling. However, these data—especially, the strong temperature dependence of reaction rates and the variable temperature dependence of KIEs—are consistent with other tunnelling models (termed full tunnelling models), in which protein and/or substrate fluctuations generate a configuration compatible with tunnelling. These models accommodate substrate/protein (environment) fluctuations required to attain a configuration with degenerate nuclear quantum states and, when necessary, motion required to increase the probability of tunnelling in these states. Furthermore, tunnelling mechanisms in enzymes are supported by atomistic computational studies performed within the framework of modern TST, which incorporates quantum nuclear effects.

The active site of an enzyme is surrounded by a fluctuating environment of protein and solvent conformational states, and a realistic calculation of chemical reaction rates and kinetic isotope effects of enzyme-catalyzed reactions must take account of this environmental diversity. Ensemble-averaged variational transition state theory with multidimensional tunneling (EA-VTST/MT) was developed as a way to carry out such calculations. This theory incorporates ensemble averaging, quantized vibrational energies, energy, tunneling, and recrossing of transition state dividing surfaces in a systematic way. It has been applied successfully to a number of hydrogen-, proton-, and hydride-transfer reactions. The theory also exposes the set of effects that should be considered in reliable rate constants calculations. We first review the basic theory and the steps in the calculation. A key role is played by the generalized free energy of activation profile, which is obtained by quantizing the classical potential of mean force as a function of a reaction coordinate because the one-way flux through the transition state dividing surface can be written in terms of the generalized free energy of activation. A recrossing transmission coefficient accounts for the difference between the one-way flux through the chosen transition state dividing surface and the net flux, and a tunneling transmission coefficient converts classical motion along the reaction coordinate to quantum mechanical motion. The tunneling calculation is multidimensional, accounting for the change in vibrational frequencies along the tunneling path and shortening of the tunneling path with respect to the minimum energy path (MEP), as promoted by reaction-path curvature. The generalized free energy of activation and the transmission coefficients both involve averaging over an ensemble of reaction paths and conformations, and this includes the coupling of protein motions to the rearrangement of chemical bonds in a statistical mechanically correct way. The standard deviations of the transmissions coefficients provide information on the diversity of the distribution of reaction paths, barriers, and protein conformations along the members of an ensemble of reaction paths passing through the transition state. We first illustrate the theory by discussing the application to both wild-type and mutant Escherichia coli dihydrofolate reductase and hyperthermophilic Thermotoga maritima dihydrofolate reductase (DHFR); DHFR is of special interest because the protein conformational changes have been widely studied. Then we present shorter discussions of several other applications of EA-VTST/MT to transfer of protons, hydrogen atoms, and hydride ions and their deuterated analogs. Systems discussed include hydride transfer in alcohol dehydrogenase, xylose isomerase, and thymidylate synthase, proton transfer in methylamine dehydrogenase, hydrogen atom transfer in methylmalonyl-CoA mutase, and nucleophilic substitution in haloalkane dehalogenase and two-dimensional potentials of mean force for potentially coupled proton and hydride transfer in the β-oxidation of butyryl-coenzyme A catalyzed by short-chain acyl-CoA dehydrogenase and in the pyruvate to lactate transformation catalyzed by lactate dehydrogenase.

ALOX15 (12/15-lipoxygenase) orthologs have been implicated in maturational degradation of intracellular organelles and in the biosynthesis of antiinflammatory and proresolving eicosanoids. Here we hypothesized that lower mammals (mice, rats, pigs) express 12-lipoxygenating ALOX15 orthologs. In contrast, 15-lipoxygenating isoforms are found in higher primates (orangutans, men), and these results suggest an evolution of ALOX15 specificity. To test this hypothesis we first cloned and characterized ALOX15 orthologs of selected Catarrhini representing different stages of late primate evolution and found that higher primates (men, chimpanzees) express 15-lipoxygenating orthologs. In contrast, lower primates (baboons, rhesus monkeys) express 12-lipoxygenating enzymes. Gibbons, which are flanked in evolution by rhesus monkeys (12-lipoxygenating ALOX15) and orangutans (15-lipoxygenating ALOX15), express an ALOX15 ortholog with pronounced dual specificity. To explore the driving force for this evolutionary alterations, we quantified the lipoxin synthase activity of 12-lipoxygenating (rhesus monkey, mouse, rat, pig, humIle418Ala) and 15-lipoxygenating (man, chimpanzee, orangutan, rabbit, ratPhe353Ala) ALOX15 variants and found that, when normalized to their arachidonic acid oxygenase activities, the lipoxin synthase activities of 15-lipoxygenating ALOX15 variants were more than fivefold higher (P < 0.01). Comparative molecular dynamics simulations and quantum mechanics/molecular mechanics calculations indicated that, for the 15-lipoxygenating rabbit ALOX15, the energy barrier for C13-hydrogen abstraction (15-lipoxygenation) was 17 kJ/mol lower than for arachidonic acid 12-lipoxygenation. In contrast, for the 12-lipoxygenating Ile418Ala mutant, the energy barrier for 15-lipoxygenation was 10 kJ/mol higher than for 12-lipoxygenation. Taken together, our data suggest an evolution of ALOX15 specificity, which is aimed at optimizing the biosynthetic capacity for antiinflammatory and proresolving lipoxins.