Arginine Kinase Activates Arginine for Phosphorylation by Pyramidalization and Polarization.

Arginine phosphorylation plays numerous roles throughout biology. Arginine kinase (AK) catalyzes the delivery of an anionic phosphoryl group (PO3-) from ATP to a planar, trigonal nitrogen in a guanidinium cation. Density functional theory (DFT) calculations have yielded a model of the transition state (TS) for the AK-catalyzed reaction. They reveal a network of over 50 hydrogen bonds that delivers unprecedented pyramidalization and out-of-plane polarization of the arginine guanidinium nitrogen (Nη2) and aligns the electron density on Nη2 with the scissile P-O bond, leading to in-line phosphoryl transfer via an associative mechanism. In the reverse reaction, the hydrogen-bonding network enforces the conformational distortion of a bound phosphoarginine substrate to increase the basicity of Nη2. This enables Nη2 protonation, which triggers PO3- migration to generate ATP. This polarization-pyramidalization of nitrogen in the arginine side chain is likely a general phenomenon that is exploited by many classes of enzymes mediating the post-translational modification of arginine.

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer.

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

A non-canonical nucleophile unlocks a new mechanistic pathway in a designed enzyme.

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 (BHMeHis1.8) that is more than an order of magnitude more active than our earlier MBHase. Interestingly, although MeHis23 to His substitution in BHMeHis1.8 reduces activity by 4-fold, the resulting His containing variant is still a potent MBH biocatalyst. However, analysis of the BHMeHis1.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.

Mechanistic implications of the ternary complex structural models for the photoenzyme protochlorophyllide oxidoreductase.

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.

Decoding Catalysis by Terpene Synthases.

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 Rev2017, 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.

Tuning of B12 photochemistry in the CarH photoreceptor to avoid radical photoproducts.

Time-resolved infrared spectroscopy reveals the flow of electron density through coenzyme B12 in the light-activated, bacterial transcriptional regulator, CarH. The protein stabilises a series of charge transfer states that result in a photoresponse that avoids reactive, and potentially damaging, radical photoproducts.

Direct Comparison between Förster Resonance Energy Transfer and Light-Induced Triplet-Triplet Electron Resonance Spectroscopy.

To carry out reliable and comprehensive structural investigations, the exploitation of different complementary techniques is required. Here, we report that dual triplet-spin/fluorescent labels enable the first parallel distance measurements by electron spin resonance (ESR) and Förster resonance energy transfer (FRET) on exactly the same molecules with orthogonal chromophores, allowing for direct comparison. An improved light-induced triplet-triplet electron resonance method with 2-color excitation is used, improving the signal-to-noise ratio of the data and yielding a distance distribution that provides greater insight than the single distance resulting from FRET.

Redox driven B12-ligand switch drives CarH photoresponse.

CarH is a coenzyme B12-dependent photoreceptor involved in regulating carotenoid biosynthesis. How light-triggered cleavage of the B12 Co-C bond culminates in CarH tetramer dissociation to initiate transcription remains unclear. Here, a series of crystal structures of the CarH B12-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 B12-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 B12 photo- and redox-chemistry to drive large-scale conformational changes through stepwise Co-ligation changes.

Heterometallic lanthanide complexes with site-specific binding that enable simultaneous visible and NIR-emission.

Macrocyclic lanthanide complexes have become widely developed due to their distinctive luminescence characteristics and wide range of applications in biological imaging. However, systems with sufficient brightness and metal selectivity can be difficult to produce on a molecular scale. Presented herein is the stepwise introduction of differing lanthanide ions in a bis-DO3A/DTPA scaffold to afford three trinuclear bimetallic [Ln2Ln'] lanthanide complexes with site-specific, controlled binding [(Yb2Tb), (Eu2Tb), (Yb2Eu)]. The complexes display simultaneous emission from all LnIII centers across the visible (TbIII, EuIII) and near infra-red (YbIII) spectrum when excited via phenyl ligand sensitization at a wide range of temperatures and are consequently of interest for exploiting imaging in the near infra-red II biological window. Analysis of lifetime data over a range of excitation regimes reveals intermetallic communication between TbIII and EuIII centers and further develops the understanding of multimetallic lanthanide complexes.

Efficient NADPH-dependent dehalogenation afforded by a self-sufficient reductive dehalogenase.

Reductive dehalogenases are corrinoid and iron-sulfur cluster-containing enzymes that catalyze the reductive removal of a halogen atom. The oxygen-sensitive and membrane-associated nature of the respiratory reductive dehalogenases has hindered their detailed kinetic study. In contrast, the evolutionarily related catabolic reductive dehalogenases are oxygen tolerant, with those that are naturally fused to a reductase domain with similarity to phthalate dioxygenase presenting attractive targets for further study. We present efficient heterologous expression of a self-sufficient catabolic reductive dehalogenase from Jhaorihella thermophila in Escherichia coli. Combining the use of maltose-binding protein as a solubility-enhancing tag with the btuCEDFB cobalamin uptake system affords up to 40% cobalamin occupancy and a full complement of iron-sulfur clusters. The enzyme is able to efficiently perform NADPH-dependent dehalogenation of brominated and iodinated phenolic compounds, including the flame retardant tetrabromobisphenol, under both anaerobic and aerobic conditions. NADPH consumption is tightly coupled to product formation. Surprisingly, corresponding chlorinated compounds only act as competitive inhibitors. Electron paramagnetic resonance spectroscopy reveals loss of the Co(II) signal observed in the resting state of the enzyme under steady-state conditions, suggesting accumulation of Co(I)/(III) species prior to the rate-limiting step. In vivo reductive debromination activity is readily observed, and when the enzyme is expressed in E. coli strain W, supports growth on 3-bromo-4-hydroxyphenylacetic as a sole carbon source. This demonstrates the potential for catabolic reductive dehalogenases for future application in bioremediation.

How Is Substrate Halogenation Triggered by the Vanadium Haloperoxidase from Curvularia inaequalis?

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 H2O2. 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.

Reactivity Differences of Trigonal Pyramidal Nonheme Iron(IV)-Oxo and Iron(III)-Oxo Complexes: Experiment and Theory.

High-valent metal-oxo species play critical roles in enzymatic catalysis yet their properties are still poorly understood. In this work we report a combined experimental and computational study into biomimetic iron(IV)-oxo and iron(III)-oxo complexes with tight second-coordination sphere environments that restrict substrate access. The work shows that the second-coordination sphere slows the hydrogen atom abstraction step from toluene dramatically and the kinetics is zeroth order in substrate. However, the iron(II)-hydroxo that is formed has a low reduction potential and hence cannot do OH rebound favorably. The tolyl radical in solution then reacts further with alternative reaction partners. By contrast, the iron(IV)-oxo species reacts predominantly through OH rebound to form alcohol products. Our studies show that the oxidation state of the metal influences reactivities and selectivities with substrate dramatically and that enzymes will likely need an iron(IV) center to catalyze C-H hydroxylation reactions.

Photoinduced Electron Transfer from a 1,4,5,6-Tetrahydro Nicotinamide Adenine Dinucleotide (Phosphate) Analogue to Oxidized Flavin in an Ene-Reductase Flavoenzyme.

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)H4 to FMN, generating a NAD(P)H4 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.

How a 10-epi-Cubebol Synthase Avoids Premature Reaction Quenching to Form a Tricyclic Product at High Purity.

Terpenes are the largest class of natural products and are attractive targets in the fuel, fragrance, pharmaceutical, and flavor industries. Harvesting terpenes from natural sources is environmentally intensive and often gives low yields and purities, requiring further downstream processing. Engineered terpene synthases (TSs) offer a solution to these problems, but the low sequence identity and high promiscuity among TSs are major challenges for targeted engineering. Rational design of TSs requires identification of key structural and chemical motifs that steer product outcomes. Producing the sesquiterpenoid 10-epi-cubebol from farnesyl pyrophosphate (FPP) requires many steps and some of Nature's most difficult chemistry. 10-epi-Cubebol synthase from Sorangium cellulosum (ScCubS) guides a highly reactive carbocationic substrate through this pathway, preventing early quenching and ensuring correct stereochemistry at every stage. The cyclizations carried out by ScCubS potentially represent significant evolutionary expansions in the chemical space accessible by TSs. Here, we present the high-resolution crystal structure of ScCubS in complex with both a trinuclear magnesium cluster and pyrophosphate. Computational modeling, experiment, and bioinformatic analysis identified residues important in steering the reaction chemistry. We show that S206 is crucial in 10-epi-cubebol synthesis by enlisting the nearby F211 to shape the active site contour and prevent the formation of early escape cadalane products. We also show that N327 and F104 control the distribution between several early-stage cations and whether the final product is derived from the germacrane, cadalane, or cubebane hydrocarbon scaffold. Using these insights, we reengineered ScCubS so that its main product was germacradien-4-ol, which derives from the germacrane, rather than the cubebane, scaffold. Our work emphasizes that mechanistic understanding of cation stabilization in TSs can be used to guide catalytic outcomes.

A Vitamin B2 -Photocatalysed Approach to Methionine Analogues.

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 B2 ) 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.

Combined Pulsed Electron Double Resonance EPR and Molecular Dynamics Investigations of Calmodulin Suggest Effects of Crowding Agents on Protein Structures.

Calmodulin (CaM) is a highly dynamic Ca2+-binding protein that exhibits large conformational changes upon binding Ca2+ and target proteins. Although it is accepted that CaM exists in an equilibrium of conformational states in the absence of target protein, the physiological relevance of an elongated helical linker region in the Ca2+-replete form has been highly debated. In this study, we use PELDOR (pulsed electron-electron double resonance) EPR measurements of a doubly spin-labeled CaM variant to assess the conformational states of CaM in the apo-, Ca2+-bound, and Ca2+ plus target peptide-bound states. Our findings are consistent with a three-state conformational model of CaM, showing a semi-open apo-state, a highly extended Ca2+-replete state, and a compact target protein-bound state. Molecular dynamics simulations suggest that the presence of glycerol, and potentially other molecular crowding agents, has a profound effect on the relative stability of the different conformational states. Differing experimental conditions may explain the discrepancies in the literature regarding the observed conformational state(s) of CaM, and our PELDOR measurements show good evidence for an extended conformation of Ca2+-replete CaM similar to the one observed in early X-ray crystal structures.

How Photoactivation Triggers Protochlorophyllide Reduction: Computational Evidence of a Stepwise Hydride Transfer during Chlorophyll Biosynthesis.

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.

Chelator-Based Parameterization of the 12-6-4 Lennard-Jones Molecular Mechanics Potential for More Realistic Metal Ion-Protein Interactions.

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 Cij coefficients. Here, we demonstrate a method of generating Cij parameters that allows parametrization of specific metal ion-ligating groups in order to tune binding energies computed by thermodynamic integration. The new Cij 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.

Molecular Determinants of Carbocation Cyclisation in Bacterial Monoterpene Synthases.

Monoterpene synthases are often promiscuous enzymes, yielding product mixtures rather than pure compounds due to the nature of the branched reaction mechanism involving reactive carbocations. Two previously identified bacterial monoterpene synthases, a linalool synthase (bLinS) and a cineole synthase (bCinS), produce nearly pure linalool and cineole from geranyl diphosphate, respectively. We used a combined experimental and computational approach to identify critical residues involved in bacterial monoterpenoid synthesis. Phe77 is essential for bCinS activity, guiding the linear carbocation intermediate towards the formation of the cyclic α-terpinyl intermediate; removal of the aromatic ring results in variants that produce acyclic products only. Computational chemistry confirmed the importance of Phe77 in carbocation stabilisation. Phe74, Phe78 and Phe179 are involved in maintaining the active site shape in bCinS without a specific role for the aromatic ring. Phe295 in bLinS, and the equivalent Ala301 in bCinS, are essential for linalool and cineole formation, respectively. Where Phe295 places steric constraints on the carbocation intermediates, Ala301 is essential for bCinS initial cyclisation and activity. Our multidisciplinary approach gives unique insights into how carefully placed amino acid residues in the active site can direct carbocations down specific paths, by placing steric constraints or offering stabilisation via cation-π interactions.

Engineering an efficient and enantioselective enzyme for the Morita-Baylis-Hillman reaction.

The combination of computational design and directed evolution could offer a general strategy to create enzymes with new functions. So far, this approach has delivered enzymes for a handful of model reactions. Here we show that new catalytic mechanisms can be engineered into proteins to accelerate more challenging chemical transformations. Evolutionary optimization of a primitive design afforded an efficient and enantioselective enzyme (BH32.14) for the Morita-Baylis-Hillman (MBH) reaction. BH32.14 is suitable for preparative-scale transformations, accepts a broad range of aldehyde and enone coupling partners and is able to promote selective monofunctionalizations of dialdehydes. Crystallographic, biochemical and computational studies reveal that BH32.14 operates via a sophisticated catalytic mechanism comprising a His23 nucleophile paired with a judiciously positioned Arg124. This catalytic arginine shuttles between conformational states to stabilize multiple oxyanion intermediates and serves as a genetically encoded surrogate of privileged bidentate hydrogen-bonding catalysts (for example, thioureas). This study demonstrates that elaborate catalytic devices can be built from scratch to promote demanding multi-step processes not observed in nature.