Ring Contraction Catalyzed by the Metal-Dependent Radical SAM Enzyme: 7-Carboxy-7-deazaguanine Synthase from B. multivorans. Theoretical Insights into the Reaction Mechanism and the Influence of Metal Ions

Cysteine Radical and Glutamate Collaboratively Enable C-H Bond Activation and C-N Bond Cleavage in a Glycyl Radical Enzyme HplG

Hydroxyprolines are abundant in nature and widely utilized by many living organisms. Isomerization of trans-4-hydroxy-d-proline (t4D-HP) to generate 2-amino-4-ketopentanoate has been found to need a glycyl radical enzyme HplG, which catalyzes the cleavage of the C-N bond, while dehydration of trans-4-hydroxy-l-proline involves a homologous enzyme of HplG. Herein, molecular dynamics simulations and quantum mechanics/molecular mechanics (QM/MM) calculations are employed to understand the reaction mechanism of HplG. Two possible reaction pathways of HplG have been explored to decipher the origin of its chemoselectivity. The QM/MM calculations reveal that the isomerization proceeds via an initial hydrogen shift from the Cγ site of t4D-HP to a catalytic cysteine radical, followed by cleavage of the Cδ-N bond in t4D-HP to form a radical intermediate that captures a hydrogen atom from the cysteine. Activation of the Cδ-H bond in t4D-HP to bring about dehydration of t4D-HP possesses an extremely high energy barrier, thus rendering the dehydration pathway implausible in HplG. On the basis of the current calculations, conserved residue Glu429 plays a pivotal role in the isomerization pathway: the hydrogen bonding between it and t4D-HP weakens the hydroxyalkyl Cγ-Hγ bond, and it acts as a proton acceptor to trigger the cleavage of the C-N bond in t4D-HP. Our current QM/MM calculations rationalize the origin of the experimentally observed chemoselectivity of HplG and propose an H-bond-assisted bond activation strategy in radical-containing enzymes. These findings have general implications on radical-mediated enzymatic catalysis and expand our understanding of how nature wisely and selectively activates the C-H bond to modulate catalytic selectivity.

Queuosine biosynthetic enzyme, QueE moonlights as a cell division regulator

In many organisms, stress responses to adverse environments can trigger secondary functions of certain proteins by altering protein levels, localization, activity, or interaction partners. Escherichia coli cells respond to the presence of specific cationic antimicrobial peptides by strongly activating the PhoQ/PhoP two-component signaling system, which regulates genes important for growth under this stress. As part of this pathway, a biosynthetic enzyme called QueE, which catalyzes a step in the formation of queuosine (Q) tRNA modification is upregulated. When cellular QueE levels are high, it co-localizes with the central cell division protein FtsZ at the septal site, blocking division and resulting in filamentous growth. Here we show that QueE affects cell size in a dose-dependent manner. Using alanine scanning mutagenesis of amino acids in the catalytic active site, we pinpoint particular residues in QueE that contribute distinctly to each of its functions - Q biosynthesis or regulation of cell division, establishing QueE as a moonlighting protein. We further show that QueE orthologs from enterobacteria like Salmonella typhimurium and Klebsiella pneumoniae also cause filamentation in these organisms, but the more distant counterparts from Pseudomonas aeruginosa and Bacillus subtilis lack this ability. By comparative analysis of E. coli QueE with distant orthologs, we elucidate a unique region in this protein that is responsible for QueEs secondary function as a cell division regulator. A dual-function protein like QueE is an exception to the conventional model of one gene, one enzyme, one function, which has divergent roles across a range of fundamental cellular processes including RNA modification and translation to cell division and stress response.

Queuosine biosynthetic enzyme, QueE moonlights as a cell division regulator

In many organisms, stress responses to adverse environments can trigger secondary functions of certain proteins by altering protein levels, localization, activity, or interaction partners. Escherichia coli cells respond to the presence of specific cationic antimicrobial peptides by strongly activating the PhoQ/PhoP two-component signaling system, which regulates genes important for growth under this stress. As part of this pathway, a biosynthetic enzyme called QueE, which catalyzes a step in the formation of queuosine (Q) tRNA modification is upregulated. When cellular QueE levels are high, it co-localizes with the central cell division protein FtsZ at the septal site, blocking division and resulting in filamentous growth. Here we show that QueE affects cell size in a dose-dependent manner. Using alanine scanning mutagenesis of amino acids in the catalytic active site, we pinpoint residues in QueE that contribute distinctly to each of its functions-Q biosynthesis or regulation of cell division, establishing QueE as a moonlighting protein. We further show that QueE orthologs from enterobacteria like Salmonella typhimurium and Klebsiella pneumoniae also cause filamentation in these organisms, but the more distant counterparts from Pseudomonas aeruginosa and Bacillus subtilis lack this ability. By comparative analysis of E. coli QueE with distant orthologs, we elucidate a unique region in this protein that is responsible for QueE's secondary function as a cell division regulator. A dual-function protein like QueE is an exception to the conventional model of "one gene, one enzyme, one function", which has divergent roles across a range of fundamental cellular processes including RNA modification and translation to cell division and stress response.

If It Is Hard, It Is Worth Doing: Engineering Radical Enzymes from Anaerobes

With a pressing need for sustainable chemistries, radical enzymes from anaerobes offer a shortcut for many chemical transformations and deliver highly sought-after functionalizations such as late-stage C-H functionalization, C-C bond formation, and carbon-skeleton rearrangements, among others. The challenges in handling these oxygen-sensitive enzymes are reflected in their limited industrial exploitation, despite what they may deliver. With an influx of structures and mechanistic understanding, the scope for designed radical enzymes to deliver wanted processes becomes ever closer. Combined with new advances in computational methods and workflows for these complex systems, the outlook for an increased use of radical enzymes in future processes is exciting.

Activation of Glycyl Radical Enzymes─Multiscale Modeling Insights into Catalysis and Radical Control in a Pyruvate Formate-Lyase-Activating Enzyme

Pyruvate formate-lyase (PFL) is a glycyl radical enzyme (GRE) playing a pivotal role in the metabolism of strict and facultative anaerobes. Its activation is carried out by a PFL-activating enzyme, a member of the radical S-adenosylmethionine (rSAM) superfamily of metalloenzymes, which introduces a glycyl radical into the Gly radical domain of PFL. The activation mechanism is still not fully understood and is structurally based on a complex with a short model peptide of PFL. Here, we present extensive molecular dynamics simulations in combination with quantum mechanics/molecular mechanics (QM/MM)-based kinetic and thermodynamic reaction evaluations of a more complete activation model comprising the 49 amino acid long C-terminus region of PFL. We reveal the benefits and pitfalls of the current activation model, providing evidence that the bound peptide conformation does not resemble the bound protein-protein complex conformation with PFL, with implications for the activation process. Substitution of the central glycine with (S)- and (R)-alanine showed excellent binding of (R)-alanine over unstable binding of (S)-alanine. Radical stabilization calculations indicate that a higher radical stability of the glycyl radical might not be the sole origin of the evolutionary development of GREs. QM/MM-derived radical formation kinetics further demonstrate feasible activation barriers for both peptide and C-terminus activation, demonstrating why the crystalized model peptide system is an excellent inhibitory system for natural activation. This new evidence supports the theory that GREs converged on glycyl radical formation due to the better conformational accessibility of the glycine radical loop, rather than the highest radical stability of the formed peptide radicals.

Computational Approaches: An Underutilized Tool in the Quest to Elucidate Radical SAM Dynamics

Enzymes are biological catalysts whose dynamics enable their reactivity. Visualizing conformational changes, in particular, is technically challenging, and little is known about these crucial atomic motions. This is especially problematic for understanding the functional diversity associated with the radical S-adenosyl-L-methionine (SAM) superfamily whose members share a common radical mechanism but ultimately catalyze a broad range of challenging reactions. Computational chemistry approaches provide a readily accessible alternative to exploring the time-resolved behavior of these enzymes that is not limited by experimental logistics. Here, we review the application of molecular docking, molecular dynamics, and density functional theory, as well as hybrid quantum mechanics/molecular mechanics methods to the study of these enzymes, with a focus on understanding the mechanistic dynamics associated with turnover.

Radical Stabilization Energies for Enzyme Engineering: Tackling the Substrate Scope of the Radical Enzyme QueE

Experimental assessment of catalytic reaction mechanisms and profiles of radical enzymes can be severely challenging due to the reactive nature of the intermediates and sensitivity of cofactors such as iron-sulfur clusters. Here, we present an enzyme-directed computational methodology for the assessment of thermodynamic reaction profiles and screening for radical stabilization energies (RSEs) for the assessment of catalytic turnovers in radical enzymes. We have applied this new screening method to the radical S-adenosylmethione enzyme 7-carboxy-7-deazaguanine synthase (QueE), following a detailed molecular dynamics (MD) analysis that clarifies the role of both specific enzyme residues and bound Mg²⁺, Ca²⁺, or Na⁺. The MD simulations provided the basis for a statistical approach to sample different conformational outcomes. RSE calculation at the M06-2X/6-31+G* level of theory provided the most computationally cost-effective assessment of enzyme-based energies, facilitated by an initial triage using semiempirical methods. The impact of intermolecular interactions on RSE was clearly established, and application to the assessment of potential alternative substrates (focusing on radical clock type rearrangements) proposes a selection of carbon-substituted analogues that would react to afford cyclopropylcarbinyl radical intermediates as candidates for catalytic turnover by QueE.

Metabolic functions of the human gut microbiota: the role of metalloenzymes

Covering: up to the end of 2017 The human body is composed of an equal number of human and microbial cells. While the microbial community inhabiting the human gastrointestinal tract plays an essential role in host health, these organisms have also been connected to various diseases. Yet, the gut microbial functions that modulate host biology are not well established. In this review, we describe metabolic functions of the human gut microbiota that involve metalloenzymes. These activities enable gut microbial colonization, mediate interactions with the host, and impact human health and disease. We highlight cases in which enzyme characterization has advanced our understanding of the gut microbiota and examples that illustrate the diverse ways in which metalloenzymes facilitate both essential and unique functions of this community. Finally, we analyze Human Microbiome Project sequencing datasets to assess the distribution of a prominent family of metalloenzymes in human-associated microbial communities, guiding future enzyme characterization efforts.

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Catalytic mechanism of the PrhA (V150L/A232S) double mutant involved in the fungal meroterpenoid biosynthetic pathway: a QM/MM study

QM/MM calculations reveal the mechanism of a nonheme Fe(ii)/α-ketoglutarate-dependent oxygenase involved in the fungal meroterpenoid biosynthetic pathway.

Mechanical insights into the oxidative cleavage of resveratrol catalyzed by dioxygenase NOV1 from Novosphingobium aromaticivorans: confirmation of dioxygenase mechanism by QM/MM calculations

QM/MM calculations confirm that the oxidative cleavage of resveratrol catalyzed by dioxygenase NOV1 follows the dioxygenase mechanism.

A QM/MM study of the catalytic mechanism of SAM methyltransferase RlmN from Escherichia coli

RlmN is a radical S-adenosylmethionine (SAM) enzyme that catalyzes the C2 methylation of adenosine 2503 (A2503) in 23S rRNA and adenosine 37 (A37) in several Escherichia coli transfer RNAs (tRNA). The catalytic reaction of RlmN is distinctly different from that of typical SAM-dependent methyltransferases that employs an S_N 2 mechanism, but follows a ping-pong mechanism which involves the intermediate methylation of a conserved cysteine residue. Recently, the x-ray structure of a key intermediate in the RlmN reaction has been reported, allowing us to perform combined quantum mechanics and molecular mechanics (QM/MM) calculations to delineate the reaction details of RlmN at atomic level. Starting from the Cross-Linked RlmN C118A-tRNA complex, the possible mechanisms for both the formation and the resolution of the cross-linked species (IM2) have been illuminated. On the basis of our calculations, IM2 is formed by the attack of the C355-based methylene radical on the sp² -hybridized C2 of the adenosine ring, corresponding to energy barrier of 14.4 kcal/mol, and the resolution of IM2 is confirmed to follow a radical fragmentation mechanism. The cleavage of C'-S' bond of mC355-A37 cross-link is in concert with the deprotonation of C2 by C118 residue, which is the rate-limiting step with an energy barrier of 17.4 kcal/mol. Moreover, the cleavage of C'-S' bond of IM2 can occur independently, that is, it does not require the loss of an electron of IM2 and the formation of disulfide bond between C355 and C118 as precondition. These findings would deepen the understanding of the catalysis of RlmN.

Radical Reaction Control in the AdoMet Radical Enzyme CDG Synthase (QueE): Consolidate, Destabilize, Accelerate

Controlling radical intermediates and thus catalysing and directing complex radical reactions is a central feature of S-adensosylmethionine (SAM)-dependent radical enzymes. We report ab initio and DFT calculations highlighting the specific influence of ion complexation, including Mg2+ , identified as a key catalytic component on radical stability and reaction control in 7-carboxy-7-deazaguanine synthase (QueE). Radical stabilisation energies (RSEs) of key intermediates and radical clock-like model systems of the enzyme-catalysed rearrangement of 6-carboxytetrahydropterin (CPH4), reveals a directing role of Mg2+ in destabilising both the substrate-derived radical and corresponding side reactions, with the effect that the experimentally-observed rearrangement becomes dominant over possible alternatives. Importantly, this is achieved with minimal disruption of the thermodynamics of the substrate itself, affording a novel mechanism for an enzyme to both maintain binding potential and accelerate the rearrangement step. Other mono and divalent ions were probed with only dicationic species achieving the necessary radical conformation to facilitate the reaction.

Tryptophan lyase (NosL): mechanistic insights into amine dehydrogenation and carboxyl fragment migration by QM/MM calculations

QM/MM calculations suggest two feasible pathways for the breaking of the C–C bond of the substrate. The breaking of the Cα–Cβ bond leads to the final product, whereas the cleavage of the Cα–C bond will terminate in the EPR-trapped radical intermediate.

Cleavage mechanism of the aliphatic C–C bond catalyzed by 2,4′-dihydroxyacetophenone dioxygenase from Alcaligenes sp. 4HAP: a QM/MM study

Calculations suggest that the reactant complex may firstly undergo a triplet–quintet crossing to initiate the reaction and then the subsequent chemistry occurs on the multiple-states surfaces. The key C–C bond cleavage is accompanied by an insertion reaction of oxygen radical.

Exploring the substrate specificity and catalytic mechanism of imidazolonepropionase (HutI) from Bacillus subtilis

QM/MM calculations reveal that imidazolonepropionase (HutI) specifically catalyzes the hydrolytic cleavage of (S)-enantiomer of substrate to yieldl-formiminoglutamic acid.

Insights into the catalytic mechanism of N-acetylglucosaminidase glycoside hydrolase from Bacillus subtilis: a QM/MM study

QM/MM calculations on NagZs fromBacillus subtilisfurther confirm NagZs to be glycoside phosphorylases rather than glycoside hydrolases.

Quantum mechanics and molecular mechanics study of the reaction mechanism of quorum quenching enzyme: N-acyl homoserine lactonase with C6-HSL

N-Acyl-homoserine lactonase fromOchrobactrumsp. strain (AidH) is a novel AHL (N-acyl-homoserine lactone)-lactonase that hydrolyzes the ester bond of the homoserine lactone ring of AHLs.

Radical-mediated ring contraction in the biosynthesis of 7-deazapurines

Pyrrolopyrimidine containing natural products are widely distributed in Nature. The biosynthesis of the 7-deazapurine moiety that is common to all pyrrolopyrimidines entails multiple steps, one of which is a complex radical-mediated ring contraction reaction catalyzed by CDG synthase. Herein we review the biosynthetic pathways of deazapurines, focusing on the biochemical and structural insights into CDG synthase.