Search

The photoenzyme protochlorophyllide oxidoreductase (POR) is an important enzyme for understanding biological H-transfer mechanisms. It uses light to catalyse the reduction of protochlorophyllide to chlorophyllide, a key step in chlorophyll biosynthesis. Although a wealth of spectroscopic data have provided crucial mechanistic insight, a structural rationale for POR photocatalysis has proved challenging and remains hotly debated. Recent structural models of the ternary enzyme-substrate complex, derived from crystal and electron microscopy data, show differences in the orientation of the protochlorophyllide substrate and the architecture of the POR active site, with significant implications for the catalytic mechanism. Here, we use a combination of computational and experimental approaches to investigate the compatibility of each structural model with the hypothesised reaction mechanisms and propose an alternative structural model for the cyanobacterial POR ternary complex. We show that a strictly conserved tyrosine, previously proposed to act as the proton donor in POR photocatalysis, is unlikely to be involved in this step of the reaction but is crucial for Pchlide binding. Instead, an active site cysteine is important for both hydride and proton transfer reactions in POR and is proposed to act as the proton donor, either directly or through a water-mediated network. Moreover, a conserved glutamine is important for Pchlide binding and ensuring efficient photochemistry by tuning its electronic properties, likely by interacting with the central Mg atom of the substrate. This optimal 'binding pose' for the POR ternary enzyme-substrate complex illustrates how light energy can be harnessed to facilitate enzyme catalysis by this unique enzyme., (© 2023 The Authors. The FEBS Journal published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Photoreceptor proteins utilise chromophores to sense light and trigger a biological response. The discovery that adenosylcobalamin (or coenzyme B 12 ) can act as a light-sensing chromophore heralded a new field of B 12 -photobiology. Although microbial genome analysis indicates that photoactive B 12 -binding domains form part of more complex protein architectures, regulating a range of molecular-cellular functions in response to light, experimental evidence is lacking. Here we identify and characterise a sub-family of multi-centre photoreceptors, termed photocobilins, that use B 12 and biliverdin (BV) to sense light across the visible spectrum. Crystal structures reveal close juxtaposition of the B 12 and BV chromophores, an arrangement that facilitates optical coupling. Light-triggered conversion of the B 12 affects quaternary structure, in turn leading to light-activation of associated enzyme domains. The apparent widespread nature of photocobilins implies involvement in light regulation of a wider array of biochemical processes, and thus expands the scope for B 12 photobiology. Their characterisation provides inspiration for the design of broad-spectrum optogenetic tools and next generation bio-photocatalysts., (© 2024. The Author(s).)

Directed evolution of computationally designed enzymes has provided new insights into the emergence of sophisticated catalytic sites in proteins. In this regard, we have recently shown that a histidine nucleophile and a flexible arginine can work in synergy to accelerate the Morita-Baylis-Hillman (MBH) reaction with unrivalled efficiency. Here, we show that replacing the catalytic histidine with a non-canonical N δ -methylhistidine (MeHis23) nucleophile leads to a substantially altered evolutionary outcome in which the catalytic Arg124 has been abandoned. Instead, Glu26 has emerged, which mediates a rate-limiting proton transfer step to deliver an enzyme (BH MeHis 1.8) that is more than an order of magnitude more active than our earlier MBHase. Interestingly, although MeHis23 to His substitution in BH MeHis 1.8 reduces activity by 4-fold, the resulting His containing variant is still a potent MBH biocatalyst. However, analysis of the BH MeHis 1.8 evolutionary trajectory reveals that the MeHis nucleophile was crucial in the early stages of engineering to unlock the new mechanistic pathway. This study demonstrates how even subtle perturbations to key catalytic elements of designed enzymes can lead to vastly different evolutionary outcomes, resulting in new mechanistic solutions to complex chemical transformations., (© 2024. The Author(s).)

The review by Christianson, published in 2017 on the twentieth anniversary of the emergence of the field, summarizes the foundational discoveries and key advances in terpene synthase/cyclase (TS) biocatalysis (Christianson, D. W. Chem Rev 2017 , 117 (17), 11570-11648. DOI: 10.1021/acs.chemrev.7b00287). Here, we review the TS literature published since then, bringing the field up to date and looking forward to what could be the near future of TS rational design. Many revealing discoveries have been made in recent years, building on the knowledge and fundamental principles uncovered during those initial two decades of study. We use these to explore TS reaction chemistry and see how a combined experimental and computational approach helps to decipher the complexities of TS catalysis. Revealed are a suite of catalytic motifs which control product outcome in TSs, some obvious, some more subtle. We examine each in detail, using the most recent papers and insights to illustrate how exactly this fascinating class of enzymes takes a single acyclic substrate and turns it into the many thousands of complex terpenoids found in Nature. We then explore some of the recent strategies for TS engineering, including machine learning and other data-driven approaches. From this, rational and predictive engineering of TSs, "designer terpene synthases", will begin to emerge as a realistic goal., Competing Interests: The authors declare the following competing financial interest(s): NSS is founder and shareholder of C3 Biotechnologies Ltd., (© 2023 The Authors. Published by American Chemical Society.)

CarH is a coenzyme B 12 -dependent photoreceptor involved in regulating carotenoid biosynthesis. How light-triggered cleavage of the B 12 Co-C bond culminates in CarH tetramer dissociation to initiate transcription remains unclear. Here, a series of crystal structures of the CarH B 12 -binding domain after illumination suggest formation of unforeseen intermediate states prior to tetramer dissociation. Unexpectedly, in the absence of oxygen, Co-C bond cleavage is followed by reorientation of the corrin ring and a switch from a lower to upper histidine-Co ligation, corresponding to a pentacoordinate state. Under aerobic conditions, rapid flash-cooling of crystals prior to deterioration upon illumination confirm a similar B 12 -ligand switch occurs. Removal of the upper His-ligating residue prevents monomer formation upon illumination. Combined with detailed solution spectroscopy and computational studies, these data demonstrate the CarH photoresponse integrates B 12 photo- and redox-chemistry to drive large-scale conformational changes through stepwise Co-ligation changes., (© 2023. Springer Nature Limited.)

Vanadium haloperoxidases (VHPOs) are unique enzymes in biology that catalyze a challenging halogen transfer reaction and convert a strong aromatic C-H bond into C-X (X = Cl, Br, I) with the use of a vanadium cofactor and H 2 O 2 . The VHPO catalytic cycle starts with the conversion of hydrogen peroxide and halide (X = Cl, Br, I) into hypohalide on the vanadate cofactor, and the hypohalide subsequently reacts with a substrate. However, it is unclear whether the hypohalide is released from the enzyme or otherwise trapped within the enzyme structure for the halogenation of organic substrates. A substrate-binding pocket has never been identified for the VHPO enzyme, which questions the role of the protein in the overall reaction mechanism. Probing its role in the halogenation of small molecules will enable further engineering of the enzyme and expand its substrate scope and selectivity further for use in biotechnological applications as an environmentally benign alternative to current organic chemistry synthesis. Using a combined experimental and computational approach, we elucidate the role of the vanadium haloperoxidase protein in substrate halogenation. Activity studies show that binding of the substrate to the enzyme is essential for the reaction of the hypohalide with substrate. Stopped-flow measurements demonstrate that the rate-determining step is not dependent on substrate binding but partially on hypohalide formation. Using a combination of molecular mechanics (MM) and molecular dynamics (MD) simulations, the substrate binding area in the protein is identified and even though the selected substrates (methylphenylindole and 2-phenylindole) have limited hydrogen-bonding abilities, they are found to bind relatively strongly and remain stable in a binding tunnel. A subsequent analysis of the MD snapshots characterizes two small tunnels leading from the vanadate active site to the surface that could fit small molecules such as hypohalide, halide, and hydrogen peroxide. Density functional theory studies using electric field effects show that a polarized environment in a specific direction can substantially lower barriers for halogen transfer. A further analysis of the protein structure indeed shows a large dipole orientation in the substrate-binding pocket that could enable halogen transfer through an applied local electric field. These findings highlight the importance of the enzyme in catalyzing substrate halogenation by providing an optimal environment to lower the energy barrier for this challenging aromatic halide insertion reaction., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

Recent reports have described the use of ene-reductase flavoenzymes to catalyze non-natural photochemical reactions. These studies have focused on using reduced flavoenzyme, yet oxidized flavins have superior light harvesting properties. In a binary complex of the oxidized ene-reductase pentaerythritol tetranitrate reductase with the nonreactive nicotinamide coenzyme analogs 1,4,5,6-tetrahydro NAD(P)H, visible photoexcitation of the flavin mononucleotide (FMN) leads to one-electron transfer from the NAD(P)H 4 to FMN, generating a NAD(P)H 4 cation radical and anionic FMN semiquinone. This electron transfer occurs in ∼1 ps and appears to kinetically outcompete reductive quenching from aromatic residues in the active site. Time-resolved infrared measurements show that relaxation processes appear to be largely localized on the FMN and the charge-separated state is short-lived, with relaxation, presumably via back electron transfer, occurring over ∼3-30 ps. While this demonstrates the potential for non-natural photoactivity, useful photocatalysis will likely require longer-lived excited states, which may be accessible by enzyme engineering and/or a judicious choice of substrate.

To rationally engineer the substrate scope and selectivity of flavin-dependent halogenases (FDHs), it is essential to first understand the reaction mechanism and substrate interactions in the active site. FDHs have long been known to achieve regioselectivity through an electrophilic aromatic substitution at C7 of the natural substrate Trp, but the precise role of a key active-site Lys residue remains ambiguous. Formation of hypochlorous acid (HOCl) at the cofactor-binding site is achieved by the direct reaction of molecular oxygen and a single chloride ion with reduced FAD and flavin hydroxide, respectively. HOCl is then guided 10 Å into the halogenation active site. Lys79, located in this site, has been proposed to direct HOCl toward Trp C7 through hydrogen bonding or a direct reaction with HOCl to form an -NH 2 Cl + intermediate. Here, we present the most likely mechanism for halogenation based on molecular dynamics (MD) simulations and active-site density functional theory "cluster" models of FDH PrnA in complex with its native substrate l-tryptophan, hypochlorous acid, and the FAD cofactor. MD simulations with different protonation states for key active-site residues suggest that Lys79 directs HOCl through hydrogen bonding, which is confirmed by calculations of the reaction profiles for both proposed mechanisms., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

Access to new non-canonical amino acid residues is crucial for medicinal chemistry and chemical biology. Analogues of the amino acid methionine have been far less explored-despite their use in biochemistry, pharmacology and peptide bioconjugation. This is largely due to limited synthetic access. Herein, we exploit a new disconnection to access non-natural methionines through the development of a photochemical method for the radical α-C-H functionalization of sulfides with alkenes, in water, using inexpensive and commercially-available riboflavin (vitamin B 2 ) as a photocatalyst. Our photochemical conditions allow the two-step synthesis of novel methionine analogues-by radical addition to unsaturated amino acid derivatives-and the chemoselective modification of peptide side-chains to yield non-natural methionine residues within small peptides. The mechanism of the bio-inspired flavin photocatalysis has been probed by experimental, DFT and TDDFT studies., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Angewandte Chemie published by Wiley-VCH GmbH.)

Metal ions are associated with a variety of proteins and play critical roles in a wide range of biochemical processes. There are multiple ways to study and quantify protein-metal ion interactions, including molecular dynamics simulations. Recently, the AMBER molecular mechanics forcefield was modified to include a 12-6-4 Lennard-Jones potential, which allows for a better description of nonbonded terms through the additional pairwise C ij coefficients. Here, we demonstrate a method of generating C ij parameters that allows parametrization of specific metal ion-ligating groups in order to tune binding energies computed by thermodynamic integration. The new C ij coefficients were tested on a series of chelators: ethylenediaminetetraacetic acid, nitrilotriacetic acid, egtazic acid, and the EF1 loop peptides from the proteins lanmodulin and calmodulin. The new parameters show significant improvements in computed binding energies relative to existing force fields and produce coordination numbers and ion-oxygen distances that are in good agreement with experimental values. This parametrization method should be extensible to a range of other systems and could be readily adapted to tune properties other than binding energies.

The photochemical reaction catalyzed by enzyme protochlorophyllide oxidoreductase (POR), a rare example of a photoactivated enzyme, is a crucial step during chlorophyll biosynthesis and involves the fastest known biological hydride transfer. Structures of the enzyme with bound substrate protochlorophyllide (PChlide) and coenzyme nicotinamide adenine dinucleotide phosphate (NADPH) have recently been published, opening up the possibility of using computational approaches to provide a comprehensive understanding of the excited state chemistry. Herein, we propose a complete mechanism for the photochemistry between PChlide and NADPH based on density functional theory (DFT) and time-dependent DFT calculations that is consistent with recent experimental data. In this multi-step mechanism, photoexcitation of PChlide leads to electron transfer from NADPH to PChlide, which in turn facilitates hydrogen atom transfer by weakening the breaking C-H bond. This work rationalizes how photoexcitation facilitates hydride transfer in POR and has more general implications for biological hydride transfer reactions., Competing Interests: The authors declare no competing financial interest., (© 2022 American Chemical Society.)

Organisms across the natural world respond to their environment through the action of photoreceptor proteins. The vitamin B 12 -dependent photoreceptor, CarH, is a bacterial transcriptional regulator that controls the biosynthesis of carotenoids to protect against photo-oxidative stress. The binding of B 12 to CarH monomers in the dark results in the formation of a homo-tetramer that complexes with DNA; B 12 photochemistry results in tetramer dissociation, releasing DNA for transcription. Although the details of the response of CarH to light are beginning to emerge, the biophysical mechanism of B 12 -binding in the dark and how this drives domain assembly is poorly understood. Here - using a combination of molecular dynamics simulations, native ion mobility mass spectrometry and time-resolved spectroscopy - we reveal a complex picture that varies depending on the availability of B 12 . When B 12 is in excess, its binding drives structural changes in CarH monomers that result in the formation of head-to-tail dimers. The structural changes that accompany these steps mean that they are rate-limiting. The dimers then rapidly combine to form tetramers. Strikingly, when B 12 is scarcer, as is likely in nature, tetramers with native-like structures can form without a B 12 complement to each monomer, with only one apparently required per head-to-tail dimer. We thus show how a bulky chromophore such as B 12 shapes protein/protein interactions and in turn function, and how a protein can adapt to a sub-optimal availability of resources. This nuanced picture should help guide the engineering of B 12 -dependent photoreceptors as light-activated tools for biomedical applications., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Enzymatic photocatalysis is seldom used in biology. Photocatalysis by light-dependent protochlorophyllide oxidoreductase (LPOR)-one of only a few natural light-dependent enzymes-is an exception, and is responsible for the conversion of protochlorophyllide to chlorophyllide in chlorophyll biosynthesis. Photocatalysis by LPOR not only regulates the biosynthesis of the most abundant pigment on Earth but it is also a 'master switch' in photomorphogenesis in early plant development. Following illumination, LPOR promotes chlorophyll production, plastid membranes are transformed and the photosynthetic apparatus is established. Given these remarkable, light-induced pigment and morphological changes, the LPOR-catalysed reaction has been extensively studied from catalytic, physiological and plant development perspectives, highlighting vital, and multiple, cellular roles of this intriguing enzyme. Here, we offer a perspective in which the link between LPOR photocatalysis and plant photomorphogenesis is explored. Notable breakthroughs in LPOR structural biology have uncovered the structural-mechanistic basis of photocatalysis. These studies have clarified how photon absorption by the pigment protochlorophyllide-bound in a ternary LPOR-protochlorophyllide-NADPH complex-triggers photocatalysis and a cascade of complex molecular and cellular events that lead to plant morphological changes. Photocatalysis is therefore the master switch responsible for early-stage plant development and ultimately life on Earth.

Molecular dynamics (MD) simulations are a popular method of studying protein structure and function, but are unable to reliably sample all relevant conformational space in reasonable computational timescales. A range of enhanced sampling methods are available that can improve conformational sampling, but these do not offer a complete solution. We present here a proof-of-principle method of combining MD simulation with machine learning to explore protein conformational space. An autoencoder is used to map snapshots from MD simulations onto a user-defined conformational landscape defined by principal components analysis or specific structural features, and we show that we can predict, with useful accuracy, conformations that are not present in the training data. This method offers a new approach to the prediction of new low energy/physically realistic structures of conformationally dynamic proteins and allows an alternative approach to enhanced sampling of MD simulations., (© 2021 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Protochlorophyllide oxidoreductase (POR) catalyses reduction of protochlorophyllide (Pchlide) to chlorophyllide, a light-dependent reaction of chlorophyll biosynthesis. POR is also important in plant development as it is the main constituent of prolamellar bodies in etioplast membranes. Prolamellar bodies are highly organised, paracrystalline structures comprising aggregated oligomeric structures of POR-Pchlide-NADPH complexes. How these oligomeric structures are formed and the role of Pchlide in oligomerisation remains unclear. POR crystal structures highlight two peptide regions that form a 'lid' to the active site, and undergo conformational change on binding Pchlide. Here, we show that Pchlide binding triggers formation of large oligomers of POR using size exclusion chromatography. A POR 'octamer' has been isolated and its structure investigated by cryo-electron microscopy at 7.7 Å resolution. This structure shows that oligomer formation is most likely driven by the interaction of amino acid residues in the highly conserved lid regions. Computational modelling indicates that Pchlide binding stabilises exposure of hydrophobic surfaces formed by the lid regions, which supports POR dimerisation and ultimately oligomer formation. Studies with variant PORs demonstrate that lid residues are involved in substrate binding and photocatalysis. These highly conserved lid regions therefore have a dual function. The lid residues position Pchlide optimally to enable photocatalysis. Following Pchlide binding, they also enable POR oligomerisation - a process that is reversed through subsequent photocatalysis in the early stages of chloroplast development., (© 2020 The Authors. The FEBS Journal published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

The cytochrome P450 monooxygenase P450 BM3 (BM3) is a biotechnologically important and versatile enzyme capable of producing important compounds such as the medical drugs pravastatin and artemether, and the steroid hormone testosterone. BM3 is a natural fusion enzyme comprising two major domains: a cytochrome P450 (heme-binding) catalytic domain and a NADPH-cytochrome P450 reductase (CPR) domain containing FAD and FMN cofactors in distinct domains of the CPR. A crystal structure of full-length BM3 enzyme is not available in its monomeric or catalytically active dimeric state. In this study, we provide detailed insights into the protein-protein interactions that occur between domains in the BM3 enzyme and characterize molecular interactions within the BM3 dimer by using several hybrid mass spectrometry (MS) techniques, namely native ion mobility MS (IM-MS), collision-induced unfolding (CIU), and hydrogen-deuterium exchange MS (HDX-MS). These methods enable us to probe the structure, stoichiometry, and domain interactions in the ∼240 kDa BM3 dimeric complex. We obtained high-sequence coverage (88-99%) in the HDX-MS experiments for full-length BM3 and its component domains in both the ligand-free and ligand-bound states. We identified important protein interaction sites, in addition to sites corresponding to heme-CPR domain interactions at the dimeric interface. These findings bring us closer to understanding the structure and catalytic mechanism of P450 BM3., (© 2020 Jeffreys et al.)

Monoterpenoids are industrially important natural products with applications in the flavours, fragrances, fuels and pharmaceutical industries. Most monoterpenoids are produced by plants, but recently two bacterial monoterpene synthases have been identified, including a cineole synthase (bCinS). Unlike plant cineole synthases, bCinS is capable of producing nearly pure cineole from geranyl diphosphate in a complex cyclisation cascade that is tightly controlled. Here we have used a multidisciplinary approach to show that Asn305 controls water attack on the α-terpinyl cation and subsequent cyclisation and deprotonation of the α-terpineol intermediate, key steps in the cyclisation cascade which direct product formation towards cineole. Mutation of Asn305 results in variants that no longer produce α-terpineol or cineole. Molecular dynamics simulations revealed that water coordination is disrupted in all variants tested. Quantum mechanics calculations indicate that Asn305 is most likely a (transient) proton acceptor for the final deprotonation step. Our synergistic approach gives unique insight into how a single residue, Asn305, tames the promiscuous chemistry of monoterpene synthase cyclisation cascades. It does this by tightly controlling the final steps in cineole formation catalysed by bCinS to form a single hydroxylated monoterpene product., (© 2019 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

While it is well established that thermally-activated quantum mechanical tunnelling of light particles (electrons and light atoms, typically hydrogen) plays a role in many enzyme-catalysed reactions, there are few definitive experimental signatures of atomic tunnelling and no clear methods of directly estimating the relative tunnelling contribution from typical experimental data. As most enzyme reactions involve the binding/capture of freely diffusing substrate(s), reactions are typically initiated by mixing and experimental conditions must then be compatible with liquid water (the solvent). This precludes the classic test of tunnelling: the observation of temperature-independent rate constants at cryogenic temperatures. Instead, H-tunnelling is usually inferred from kinetic isotope effects that are larger than the semiclassical limit. Often, the temperature dependence of the reaction is also measured over the experimentally accessible range (∼278-313 K for mesophilic enzymes), with resulting data analysed and interpreted using variations of Arrhenius, Eyring or Marcus theory. The apparent Arrhenius and Eyring activation parameters allow some quantitative comparison of different reactions, but do not directly provide any information about tunnelling, while the validity of parameters derived from non-adiabatic models such as Marcus theory are questionable due to the partially adiabatic nature of these reactions. Here, we use the correlation found between apparent activation enthalpy and entropy across several series of enzyme variants and tunnelling contributions determined using computational chemistry in an attempt to question and define new signatures of hydrogen tunnelling, which can be used to interpret typical experimental kinetic data measured for enzyme-catalysed reactions.

The enzyme protochlorophyllide oxidoreductase (POR) catalyses a light-dependent step in chlorophyll biosynthesis that is essential to photosynthesis and, ultimately, all life on Earth 1-3 . POR, which is one of three known light-dependent enzymes 4,5 , catalyses reduction of the photosensitizer and substrate protochlorophyllide to form the pigment chlorophyllide. Despite its biological importance, the structural basis for POR photocatalysis has remained unknown. Here we report crystal structures of cyanobacterial PORs from Thermosynechococcus elongatus and Synechocystis sp. in their free forms, and in complex with the nicotinamide coenzyme. Our structural models and simulations of the ternary protochlorophyllide-NADPH-POR complex identify multiple interactions in the POR active site that are important for protochlorophyllide binding, photosensitization and photochemical conversion to chlorophyllide. We demonstrate the importance of active-site architecture and protochlorophyllide structure in driving POR photochemistry in experiments using POR variants and protochlorophyllide analogues. These studies reveal how the POR active site facilitates light-driven reduction of protochlorophyllide by localized hydride transfer from NADPH and long-range proton transfer along structurally defined proton-transfer pathways.

Catechol- O -methyltransferase (COMT) is a model S-adenosyl-l-methionine (SAM) dependent methyl transferase, which catalyzes the methylation of catecholamine neurotransmitters such as dopamine in the primary pathway of neurotransmitter deactivation in animals. Despite extensive study, there is no consensus view of the physical basis of catalysis in COMT. Further progress requires experimental data that directly probes active site geometry, protein dynamics and electrostatics, ideally in a range of positions along the reaction coordinate. Here we establish that sinefungin, a fungal-derived inhibitor of SAM-dependent enzymes that possess transition state-like charge on the transferring group, can be used as a transition state analog of COMT when combined with a catechol. X-ray crystal structures and NMR backbone assignments of the ternary complexes of the soluble form of human COMT containing dinitrocatechol, Mg 2+ and SAM or sinefungin were determined. Comparison and further analysis with the aid of density functional theory calculations and molecular dynamics simulations provides evidence for active site "compaction", which is driven by electrostatic stabilization between the transferring methyl group and "equatorial" active site residues that are orthogonal to the donor-acceptor (pseudo reaction) coordinate. We propose that upon catecholamine binding and subsequent proton transfer to Lys 144, the enzyme becomes geometrically preorganized, with little further movement along the donor-acceptor coordinate required for methyl transfer. Catalysis is then largely facilitated through stabilization of the developing charge on the transferring methyl group via "equatorial" H-bonding and electrostatic interactions orthogonal to the donor-acceptor coordinate., Competing Interests: The authors declare no competing financial interest.

The scope for biocatalytic modification of non-native carvone derivatives for speciality intermediates has hitherto been limited. Additionally, caprolactones are important feedstocks with diverse applications in the polymer industry and new non-native terpenone-derived biocatalytic caprolactone syntheses are thus of potential value for industrial biocatalytic materials applications. Biocatalytic reduction of synthetic analogues of R-(-)-carvone with additional substituents at C3 or C6, or both C3 and C6, using three types of OYEs (OYE2, PETNR and OYE3) shows significant impact of both regio-substitution and the substrate diastereomer. Bioreduction of (-)-carvone derivatives substituted with a Me and/or OH group at C6 is highly dependent on the diastereomer of the substrate. Derivatives bearing C6 substituents larger than methyl moieties are not substrates. Computer docking studies of PETNR with both (6S)-Me and (6R)-Me substituted (-)-carvone provides a model consistent with the outcomes of bioconversion. The products of bioreduction were efficiently biotransformed by the Baeyer-Villiger monooxygenase (BVase) CHMO_Phi1 to afford novel trisubstituted lactones with complete regioselectivity to provide a new biocatalytic entry to these chiral caprolactones. This provides both new non-native polymerization feedstock chemicals, but also with enhanced efficiency and selectivity over native (+)-dihydrocarvone Baeyer-Villigerase expansion. Optimum enzymatic reactions were scaled up to 60-100 mg, demonstrating the utility for preparative biocatalytic synthesis of both new synthetic scaffold-modified dihydrocarvones and efficient biocatalytic entry to new chiral caprolactones, which are potential single-isomer chiral polymer feedstocks., (© 2019 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

UVR8 is a plant photoreceptor protein that regulates photomorphogenic and protective responses to UV light. The inactive, homodimeric state absorbs UV-B light, resulting in dissociation into monomers, which are considered to be the active state and comprise a β-propeller core domain and intrinsically disordered N- and C-terminal tails. The C terminus is required for functional binding to signaling partner COP1. To date, however, structural studies have only been conducted with the core domain where the terminal tails have been truncated. Here, we report structural investigations of full-length UVR8 using native ion mobility mass spectrometry adapted for photoactivation. We show that, while truncated UVR8 photoconverts from a single conformation of dimers to a single monomer conformation, the full-length protein exists in numerous conformational families. The full-length dimer adopts both a compact state and an extended state where the C terminus is primed for activation. In the monomer the extended C terminus destabilizes the core domain to produce highly extended yet stable conformations, which we propose are the fully active states that bind COP1. Our results reveal the conformational diversity of full-length UVR8. We also demonstrate the potential power of native mass spectrometry to probe functionally important structural dynamics of photoreceptor proteins throughout nature., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Aggregation of therapeutic proteins is a key factor in the generation of unwanted immunogenicity, and can result in reduced serum half-life, neutralization of function and adverse health effects. There is currently little information regarding how aggregates interact with B-cell receptors or cognate antibodies at the protein sequence level, or whether non-native, aggregate-induced epitopes predominate in these interactions. Using an antibody fragment (single chain antibody variable fragment; scFv) that forms aggregates readily at low temperature, anti-scFv IgG antibody responses were generated by intraperitoneal injection of BALB/c strain mice with monomer or aggregate preparations. Aggregate-specific immunosignatures were identified by oligo-peptide microarray fine epitope mapping, using overlapping 15mer peptides based on the linear sequence of scFv, printed onto glass slides. IgG antibodies from mice immunized with aggregated scFv preferentially recognized a patch of overlapping peptides. This region mapped to a β-strand located at the interface between the V H and V L domains. Molecular dynamics simulations indicated that the V L domain is less stable than the V H domain, suggesting the interface region between the two domains becomes exposed during partial unfolding of the scFv during aggregate formation. These data are consistent with the hypothesis that epitopes from partially unfolded states are revealed, or are more fully exposed, in the aggregated state, and that this can augment the IgG antibody response. This observation offers the theoretical possibility that epitopes preferentially associated with aggregates can be identified from the anti-drug antibody serum IgG response which may, in turn, lead to better methods for detection of anti-drug antibody responses, and improved design of therapeutic proteins to control immunogenicity., (Copyright © 2018 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The realization of a synthetic biology approach to microbial (1 R ,2 S ,5 R )-( - )-menthol ( 1 ) production relies on the identification of a gene encoding an isopulegone isomerase (IPGI), the only enzyme in the Mentha piperita biosynthetic pathway as yet unidentified. We demonstrate that Δ5-3-ketosteroid isomerase (KSI) from Pseudomonas putida can act as an IPGI, producing ( R )-(+)-pulegone (( R )- 2 ) from (+)- cis -isopulegone ( 3 ). Using a robotics-driven semirational design strategy, we identified a key KSI variant encoding four active site mutations, which confer a 4.3-fold increase in activity over the wild-type enzyme. This was assisted by the generation of crystal structures of four KSI variants, combined with molecular modeling of 3 binding to identify key active site residue targets. The KSI variant was demonstrated to function efficiently within cascade biocatalytic reactions with downstream Mentha enzymes pulegone reductase and (-)-menthone:(-)-menthol reductase to generate 1 from 3 . This study introduces the use of a recombinant IPGI, engineered to function efficiently within a biosynthetic pathway for the production of 1 in microorganisms., Competing Interests: The authors declare the following competing financial interest(s): A.C., H.S.T., and N.S.S. are named inventors on a submitted patent describing the work in this paper.

Aromatic amine dehydrogenase (AADH) and related enzymes are at the heart of debates on the roles of quantum tunneling and protein dynamics in catalysis. The reaction of tryptamine in AADH involves significant quantum tunneling in the rate-limiting proton transfer step, shown by large H/D primary kinetic isotope effects (KIEs), with unusual temperature dependence. We apply correlated ab initio combined quantum mechanics/molecular mechanics (QM/MM) methods, at levels up to local coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), to calculate accurate potential energy surfaces for this reaction, which are necessary for quantitative analysis of tunneling contributions and reaction dynamics. Different levels of QM/MM treatment are tested. Multiple pathways are calculated with fully flexible transition state optimization by the climbing-image nudged elastic band method at the density functional QM/MM level. The average LCCSD(T) potential energy barriers to proton transfer are 16.7 and 14.0 kcal/mol for proton transfer to the two carboxylate atoms of the catalytic base, Asp128β. The results show that two similar, but distinct pathways are energetically accessible. These two pathways have different barriers, exothermicity and curvature, and should be considered in analyses of the temperature dependence of reaction and KIEs in AADH and other enzymes. These results provide a benchmark for this prototypical enzyme reaction and will be useful for developing empirical models, and analyzing experimental data, to distinguish between different conceptual models of enzyme catalysis.

All cryptochromes are currently classified as flavoproteins. In animals their best-described role is as components of the circadian clock. This circadian function is variable, and can be either light-dependent or -independent; the molecular origin of this difference is unknown. Type I animal cryptochromes are photoreceptors that entrain an organism's clock to its environment, whereas Type II (including mammals) regulate circadian timing in a light-independent manner. Here, we reveal that, in contrast to Type I, Type II animal cryptochromes lack the structural features to securely bind the photoactive flavin cofactor. We provide a molecular basis for the distinct circadian roles of different animal cryptochromes, which also has significant implications for the putative role of Type II cryptochromes in animal photomagnetoreception.

Distinct mechanisms have been proposed for the biological dehalogenation catalyzed by cobalamin-dependent enzymes, with two recent crystallographic studies suggesting different mechanisms based on the observed interaction between the organohalide substrate and cobalamin. In one case, involving an aromatic dibromide substrate in NpRdhA, a novel Co II -Br interaction was observed using EPR, suggesting a mechanism involving a [CoXR] adduct. However, in the case of trichloroethylene in PceA, a significantly longer Co-Cl distance was observed in X-ray crystal structures, suggesting a dissociative electron transfer mechanism. Subsequent DFT models of these reactions have not reproduced these differences in binding modes. Here, we have performed molecular docking and DFT calculations to investigate and compare the interaction between different organohalides and cobalamin in both NpRdhA and PceA. In each case, despite differences in binding in the Co II state, the reaction likely proceeds via formation of a [CoXR] adduct in the Co I state that weakens the breaking carbon-halide bond, suggesting this could be a general mechanism for cobalamin-dependent dehalogenation.

Protochlorophyllide (Pchlide), an intermediate in the biosynthesis of chlorophyll, is the substrate for the light-driven enzyme protochlorophyllide oxidoreductase. Pchlide has excited-state properties that allow it to initiate photochemistry in the enzyme active site, which involves reduction of Pchlide by sequential hydride and proton transfer. The basis of this photochemical behavior has been investigated here using a combination of time-resolved spectroscopies and density functional theory calculations of a number of Pchlide analogues with modifications to various substituent groups. A keto group on ring E is essential for excited-state charge separation in the molecule, which is the driving force for the photoreactivity of the pigment. Vibrational "fingerprints" of specific regions of the Pchlide chromophore have been assigned, allowing identification of the modes that are crucial for excited-state chemistry in the enzyme. This work provides an understanding of the structural determinants of Pchlide that are important for harnessing light energy.

Photoinitiated proton-coupled electron transfer (PCET) kinetics has been investigated in a series of four modified tyrosines linked to a ruthenium photosensitizer in acetonitrile, with each tyrosine bearing an internal hydrogen bond to a covalently linked pyridine or benzimidazole base. After correcting for differences in driving force, it is found that the intrinsic PCET rate constant still varies by 2 orders of magnitude. The differences in rates, as well as the magnitude of the kinetic isotope effect (KIE = kH/kD), both generally correlate with DFT calculated proton donor-acceptor distances. An Arrhenius analysis of temperature dependent data shows that the difference in reactivity arises primarily from differences in activation energies. We use this kinetic data to evaluate a commonly employed theoretical model for proton tunneling which includes a harmonic distribution of proton donor-acceptor distances due to vibrational motions of the molecule. Applying this model to the experimental data yields the conclusion that donor-acceptor compression is more facile in the compounds with shorter PT distance; however, this is contrary to independent calculations for the same compounds. This discrepancy is likely because the assumption in the model of Morse-shaped proton potential energy surfaces is inappropriate for (strongly) hydrogen-bonded systems. These results question the general applicability of this model. The results also suggest that a correlation of rate vs proton tunneling distance for the series of compounds is complicated by a concomitant variation of other relevant parameters.

Enzyme-catalysed H-transfer reactions are ubiquitous, yet fundamental details of these reactions remain unresolved. In this perspective, we discuss the roles of nuclear quantum tunnelling and (compressive) dynamics during these reactions. Evidence for the coupling of specific substrate and/or protein vibrations to the chemical coordinate is considered and a case is made for the combination of multiple experimental and computational/theoretical approaches when studying these reactions.

The Fe(II)- and 2-oxoglutarate (2-OG)-dependent dioxygenases comprise a large and diverse enzyme superfamily the members of which have multiple physiological roles. Despite this diversity, these enzymes share a common chemical mechanism and a core structural fold, a double-stranded β-helix (DSBH), as well as conserved active site residues. The prolyl hydroxylases are members of this large superfamily. Prolyl hydroxylases are involved in collagen biosynthesis and oxygen sensing in mammalian cells. Structural-mechanistic studies with prolyl hydroxylases have broader implications for understanding mechanisms in the Fe(II)- and 2-OG-dependent dioxygenase superfamily. Here, we describe crystal structures of an N-terminally truncated viral collagen prolyl hydroxylase (vCPH). The crystal structure shows that vCPH contains the conserved DSBH motif and iron binding active site residues of 2-OG oxygenases. Molecular dynamics simulations are used to delineate structural changes in vCPH upon binding its substrate. Kinetic investigations are used to report on reaction cycle intermediates and compare them to the closest homologues of vCPH. The study highlights the utility of vCPH as a model enzyme for broader mechanistic analysis of Fe(II)- and 2-OG-dependent dioxygenases, including those of biomedical interest.

The coenzyme B12-dependent photoreceptor protein, CarH, is a bacterial transcriptional regulator that controls the biosynthesis of carotenoids in response to light. On binding of coenzyme B12 the monomeric apoprotein forms tetramers in the dark, which bind operator DNA thus blocking transcription. Under illumination the CarH tetramer dissociates, weakening its affinity for DNA and allowing transcription. The mechanism by which this occurs is unknown. Here we describe the photochemistry in CarH that ultimately triggers tetramer dissociation; it proceeds via a cob(III)alamin intermediate, which then forms a stable adduct with the protein. This pathway is without precedent and our data suggest it is independent of the radical chemistry common to both coenzyme B12 enzymology and its known photochemistry. It provides a mechanistic foundation for the emerging field of B12 photobiology and will serve to inform the development of a new class of optogenetic tool for the control of gene expression.

Coupling of fast protein dynamics to enzyme chemistry is controversial and has ignited considerable debate, especially over the past 15 years in relation to enzyme-catalyzed H-transfer. H-transfer can occur by quantum tunneling, and the temperature dependence of kinetic isotope effects (KIEs) has emerged as the "gold standard" descriptor of these reactions. The anomalous temperature dependence of KIEs is often rationalized by invoking fast motions to facilitate H-transfer, yet crucially, direct evidence for coupled motions is lacking. The fast motions hypothesis underpinning the temperature dependence of KIEs is based on inference. Here, we have perturbed vibrational motions in pentaerythritol tetranitrate reductase (PETNR) by isotopic substitution where all non-exchangeable atoms were replaced with the corresponding heavy isotope ((13)C, (15)N, and (2)H). The KIE temperature dependence is perturbed by heavy isotope labeling, demonstrating a direct link between (promoting) vibrations in the protein and the observed KIE. Further we show that temperature-independent KIEs do not necessarily rule out a role for fast dynamics coupled to reaction chemistry. We show causality between fast motions and enzyme chemistry and demonstrate how this impacts on experimental KIEs for enzyme reactions.

The rate and kinetic isotope effect (KIE) on proton transfer during the aromatic amine dehydrogenase-catalyzed reaction with phenylethylamine shows complex pressure and temperature dependences. We are able to rationalize these effects within an environmentally coupled tunneling model based on constant pressure molecular dynamics (MD) simulations. As pressure appears to act anisotropically on the enzyme, perturbation of the reaction coordinate (donor-acceptor compression) is, in this case, marginal. Therefore, while we have previously demonstrated that pressure and temperature dependences can be used to infer H-tunneling and the involvement of promoting vibrations, these effects should not be used in the absence of atomistic insight, as they can vary greatly for different enzymes. We show that a pressure-dependent KIE is not a definitive hallmark of quantum mechanical H-tunneling during an enzyme-catalyzed reaction and that pressure-independent KIEs cannot be used to exclude tunneling contributions or a role for promoting vibrations in the enzyme-catalyzed reaction. We conclude that coupling of MD calculations with experimental rate and KIE studies is required to provide atomistic understanding of pressure effects in enzyme-catalyzed reactions.

Protein domain motion is often implicated in biological electron transfer, but the general significance of motion is not clear. Motion has been implicated in the transfer of electrons from human cytochrome P450 reductase (CPR) to all microsomal cytochrome P450s (CYPs). Our hypothesis is that tight coupling of motion with enzyme chemistry can signal "ready and waiting" states for electron transfer from CPR to downstream CYPs and support vectorial electron transfer across complex redox chains. We developed a novel approach to study the time-dependence of dynamical change during catalysis that reports on the changing conformational states of CPR. FRET was linked to stopped-flow studies of electron transfer in CPR that contains donor-acceptor fluorophores on the enzyme surface. Open and closed states of CPR were correlated with key steps in the catalytic cycle which demonstrated how redox chemistry and NADPH binding drive successive opening and closing of the enzyme. Specifically, we provide evidence that reduction of the flavin moieties in CPR induces CPR opening, whereas ligand binding induces CPR closing. A dynamic reaction cycle was created in which CPR optimizes internal electron transfer between flavin cofactors by adopting closed states and signals "ready and waiting" conformations to partner CYP enzymes by adopting more open states. This complex, temporal control of enzyme motion is used to catalyze directional electron transfer from NADPH→FAD→FMN→heme, thereby facilitating all microsomal P450-catalysed reactions. Motions critical to the broader biological functions of CPR are tightly coupled to enzyme chemistry in the human NADPH-CPR-CYP redox chain. That redox chemistry alone is sufficient to drive functionally necessary, large-scale conformational change is remarkable. Rather than relying on stochastic conformational sampling, our study highlights a need for tight coupling of motion to enzyme chemistry to give vectorial electron transfer along complex redox chains., Competing Interests: The authors have declared that no competing interests exist.

The role of dynamical effects in enzyme catalysis is both complex and widely debated. Understanding how dynamics can influence the barrier to an enzyme catalyzed reaction requires the development of new methodologies and tools. In particular compressive dynamics-the focus of this study-may decrease both the height and width of a reaction barrier. By making targeted mutations in the active site of morphinone reductase we are able to alter the equilibrium of conformational states for the reactive complex in turn altering the donor-acceptor (D-A) distance for H-transfer. The sub-A changes which we induce are monitored using novel spectroscopic and kinetic "rulers". This new way of detecting variation in D-A distance allows us to analyze trends between D-A distance and the force constant of a compressive dynamical mode. We find that as the D-A distance decreases, the force constant for a compressive mode increases. Further, we demonstrate that-contrary to current dogma-compression may not always cause the magnitude of the primary kinetic isotope effect to decrease.

It is generally accepted that enzymes catalyze reactions by lowering the apparent activation energy by transition state stabilization or through destabilization of ground states. A more controversial proposal is that enzymes can also accelerate reactions through barrier compression-an idea that has emerged from studies of H-tunneling reactions in enzyme systems. The effects of barrier compression on classical (over-the-barrier) reactions, and the partitioning between tunneling and classical reaction paths, have largely been ignored. We performed theoretical and computational studies on the effects of barrier compression on the shape of potential energy surfaces/reaction barriers for model (malonaldehyde and methane/methyl radical anion) and enzymatic (aromatic amine dehydrogenase) proton transfer systems. In all cases, we find that barrier compression is associated with an approximately linear decrease in the activation energy. For partially nonadiabatic proton transfers, we show that barrier compression enhances, to similar extents, the rate of classical and proton tunneling reactions. Our analysis suggests that barrier compression-through fast promoting vibrations, or other means-could be a general mechanism for enhancing the rate of not only tunneling, but also classical, proton transfers in enzyme catalysis., (Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Proton-coupled electron transfer (PCET) was studied in two biomimetic covalently linked Ru(bpy)(3)-tyrosine complexes with the phenolic proton hydrogen-bonded to an internal carboxylate group. The phenolic group is either a salicylic acid (o-hydroxybenzoic acid, SA) or an o-hydroxyphenyl-acetic acid (PA), where the former gives a resonance-assisted hydrogen bond. Transient absorption data allowed direct determination of the rate constant for these intramolecular, bidirectional, and concerted PCET (CEP) reactions, as a function of temperature and H/D isotope. We found, unexpectedly, that the hydrogen bond in SA is in fact weaker than the hydrogen bond in the complex with PA, which forced us to reassess an earlier hypothesis that the proton coupling term for CEP with SA is increased by a stronger hydrogen bond. Consequently, the kinetic data was modeled numerically using a quantum mechanical rate expression. Sufficient experimentally determined observables were available to give robust and well-determined parameter values. This analysis, coupled with DFT/B3LYP and MP2 calculations and MD simulations, gave a detailed insight into the parameters that control the CEP reactions, and the effect of internal hydrogen bonds. We observed that a model with a static proton-tunneling distance is unable to describe the reaction correctly, requiring unrealistic values for the equilibrium proton-tunneling distances. Instead, when promoting vibrations that modulate the proton donor-acceptor distance were included, satisfactory fits to the experimental data were obtained, with parameter values that agree with DFT calculations and MD simulations. According to these results, it is in fact the weaker hydrogen bond of SA which increases the proton coupling. The inner reorganization energy of the phenolic groups is a significant factor contributing to the CEP barriers, but this is reduced by the hydrogen bonds to 0.35 and 0.50 eV for the two complexes. The promoting vibrations increase the rate of CEP by over 2 orders of magnitude, and dramatically reduce the kinetic isotope effect from ca. 40 for the static case to a modest value of 2-3.

Electron-transfer-induced linkage isomerization was investigated in a series of bis-tridentate Ru polypyridyl complexes [Ru(L-X-OH)(Y-tpy)](2+) with ambidentate ligand L-X-OH=bpy-C(R)(OH)-py (bpy=2,2'-bipyridine; py=pyridine; R=H, Me, Ph, or tBu) and spectator ligand Y-tpy (tpy=2,2':6',2''-terpyridine; Y=p-tolyl, p-PhCO(2)Me, Cl, OEt, N-pyrrolidine). The ligand-bound ambidentate motif switches reversibly between N and O coordination in the Ru(II) and Ru(III) state, respectively. The potentials of the Ru(III/II) couple differ by about 0.5 V between the isomers, and this results in a bistable electrochemical response of the molecular switches. The effects of structural modifications in form of substituents on the linking carbon atom of the ambidentate ligand and on the central pyridine moiety of the spectator ligand were investigated by electrochemical and computational methods. Differences in isomerization behavior span six orders of magnitude in rate constants and two orders of magnitude in equilibrium constants. The results can be interpreted in terms of steric and electronic substituent effects and their influence on rotational barriers, ligation geometry, and electron deficiency of the metal center.

The role of promoting vibrations in enzymic reactions involving hydrogen tunnelling is contentious. While models incorporating such promoting vibrations have successfully reproduced and explained experimental observations, it has also been argued that such vibrations are not part of the catalytic effect. In this study, we have employed combined quantum mechanical/molecular mechanical methods with molecular dynamics and potential energy surface calculations to investigate how enzyme and substrate motion affects the energy barrier to proton transfer for the rate-limiting H-transfer step in aromatic amine dehydrogenase (AADH) with tryptamine as substrate. In particular, the conformation of the iminoquinone adduct induced by AADH was found to be essential for a promoting vibration identified previously-this lowers significantly the 'effective' potential energy barrier, that is the barrier which remains to be surmounted following collective, thermally equilibrated motion attaining a quantum degenerate state of reactants and products. When the substrate adopts a conformation similar to that in the free iminoquinone, this barrier was found to increase markedly. This is consistent with AADH facilitating the H-transfer event by holding the substrate in a conformation that induces a promoting vibration.

The physical basis of the catalytic power of enzymes remains contentious despite sustained and intensive research efforts. Knowledge of enzyme catalysis is predominantly descriptive, gained from traditional protein crystallography and solution studies. Our goal is to understand catalysis by developing a complete and quantitative picture of catalytic processes, incorporating dynamic aspects and the role of quantum tunnelling. Embracing ideas that we have spearheaded from our work on quantum mechanical tunnelling effects linked to protein dynamics for H-transfer reactions, we review our recent progress in mapping macroscopic kinetic descriptors to an atomistic understanding of dynamics linked to biological H-tunnelling reactions.

Hydrogen transfer, an essential component of most biological reactions, is a quantum problem. However, the proposed role of compressive motion in promoting enzymatic H-transfer is contentious. Using molecular dynamics simulations and density functional theory (DFT) calculations, we show that, during proton tunneling in the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase (AADH), a sub-picosecond promoting vibration is inherent to the iminoquinone intermediate. We show by numerical modeling that this short-range vibration, with a frequency of approximately 165 cm-1, is consistent with "gating" motion in the hydrogen tunneling model of Kuznetsov and Ulstrup (Kuznetsov, A. M.; Ulstrup, J. Can. J. Chem. 1999, 77, 1085) in an enzymatic reaction with an observed protium/deuterium kinetic isotope effect that is not measurably temperature-dependent.

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.

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(128)beta in a reaction dominated by tunneling over approximately 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.

Using the definitions of protein folds encoded in a text string, a dynamic programming algorithm was devised to compare these and identify their largest common substructure and calculate the distance (in terms of the number of edit operations) that this lay from each structure. This provided a metric on which the folds were clustered into a 'phylogenetic' tree. This construction differs from previous automatic structure clustering algorithms as it has explicit representation of the structures at 'ancestral' branching nodes, even when these have no corresponding known structure. The resulting tree was compared with that compiled by an 'expert' in the field and while there was broad agreement, differences were found that resulted from differing degrees of emphasis being placed on the types of operations that can be used to transform structures. Some concluding speculations on the relationship of such trees to the evolutionary history and folding of the proteins are advanced.