Cryo-EM and Molecular Dynamics Simulations Reveal Hidden Conformational Dynamics Controlling Ammonia Transport in Human Asparagine Synthetase

How dynamical motions in enzymes might be linked to catalytic function is of significant general interest, although almost all relevant experimental data, to date, has been obtained for enzymes with a single active site. Recent advances in X-ray crystallography and cryogenic electron microscopy offer the promise of elucidating dynamical motions for proteins that are not amenable to study using solution-phase NMR methods. Here we use 3D variability analysis (3DVA) of an EM structure for human asparagine synthetase (ASNS) in combination with atomistic molecular dynamics (MD) simulations to detail how dynamic motions of a single side chain mediates interconversion of the open and closed forms of a catalytically relevant intramolecular tunnel, thereby regulating catalytic function. Our 3DVA results are consistent with those obtained independently from MD simulations, which further suggest that formation of a key reaction intermediate acts to stabilize the open form of the tunnel in ASNS to permit ammonia translocation and asparagine formation. This conformational selection mechanism for regulating ammonia transfer in human ASNS contrasts sharply with those employed in other glutamine-dependent amidotransferases that possess a homologous glutaminase domain. Our work illustrates the power of cryo-EM to identify localized conformational changes and hence dissect the conformational landscape of large proteins. When combined with MD simulations, 3DVA is a powerful approach to understanding how conformational dynamics regulate function in metabolic enzymes with multiple active sites.

Reactivity and mechanism in chemical and synthetic biology

Physical organic chemistry and mechanistic thinking provide a strong intellectual framework for understanding the chemical logic of evolvable informational macromolecules and metabolic transformations in living organisms. These concepts have also led to numerous successes in designing and applying tools to delineate biological function in health and disease, chemical ecology and possible alternative chemistries employed by extraterrestrial life. A symposium at the 2020 Pacifichem meeting was scheduled in December 2020 to discuss designing and exploiting expanded genetic alphabets, methods to understand the biosynthesis of natural products and re-engineering primary metabolism in bacteria. The COVID-19 pandemic led to postponement of in-person discussions, with the symposium eventually being held on 20-21 December 2021 as an online event. This issue is a written record of work presented on biosynthetic pathways and enzyme catalysis, engineering microorganisms with new metabolic capabilities, and the synthesis of non-canonical, nucleobases for medical applications and for studies of alternate chemistries for living organisms. The variety of opinion pieces, reviews and original research articles provide a starting point for innovations that clarify how complex biological systems emerge from the rules of chemical reactivity and mechanism. This article is part of the themed issue 'Reactivity and mechanism in chemical and synthetic biology'.

Experimental and computational snapshots of C-C bond formation in a C-nucleoside synthase

The biosynthetic enzyme, ForT, catalyses the formation of a C-C bond between 4-amino-1H-pyrazoledicarboxylic acid and MgPRPP to produce a C-nucleoside precursor of formycin A. The transformation catalysed by ForT is of chemical interest because it is one of only a few examples in which C-C bond formation takes place via an electrophilic substitution of a small, aromatic heterocycle. In addition, ForT is capable of discriminating between the aminopyrazoledicarboxylic acid and an analogue in which the amine is replaced by a hydroxyl group; a remarkable feat given the steric and electronic similarities of the two molecules. Here we report biophysical measurements, structural biology and quantum chemical calculations that provide a detailed molecular picture of ForT-catalysed C-C bond formation and the conformational changes that are coupled to catalysis. Our findings set the scene for employing engineered ForT variants in the biocatalytic production of novel, anti-viral C-nucleoside and C-nucleotide analogues.

Novel Dihydropyrimidinone Derivatives as Potential P-Glycoprotein Modulators

P-glycoprotein (Pgp), an ATP binding cassette (ABC) transporter, is an ATP-dependent efflux pump responsible for cancer multidrug resistance. As part of efforts to identify human Pgp (hPgp) inhibitors, we prepared a series of novel triazole-conjugated dihydropyrimidinones using a synthetic approach that is well suited for obtaining compound libraries. Several of these dihydropyrimidinone derivatives modulate human P-glycoprotein (hPgp) activity with low micromolar EC₅₀ values. Molecular docking studies suggest that these compounds bind to the M-site of the transporter.

Functional Characterization of Two PLP-Dependent Enzymes Involved in Capsular Polysaccharide Biosynthesis from Campylobacter jejuni

Campylobacter jejuni is a Gram-negative, pathogenic bacterium that causes campylobacteriosis, a form of gastroenteritis. C. jejuni is the most frequent cause of food-borne illness in the world, surpassing Salmonella and E. coli. Coating the surface of C. jejuni is a layer of sugar molecules known as the capsular polysaccharide that, in C. jejuni NCTC 11168, is composed of a repeating unit of d-glycero-l-gluco-heptose, d-glucuronic acid, d-N-acetyl-galactosamine, and d-ribose. The d-glucuronic acid moiety is further amidated with either serinol or ethanolamine. It is unknown how these modifications are synthesized and attached to the polysaccharide. Here, we report the catalytic activities of two previously uncharacterized, pyridoxal phosphate (PLP)-dependent enzymes, Cj1436 and Cj1437, from C. jejuni NCTC 11168. Using a combination of mass spectrometry and nuclear magnetic resonance, we determined that Cj1436 catalyzes the decarboxylation of l-serine phosphate to ethanolamine phosphate. Cj1437 was shown to catalyze the transamination of dihydroxyacetone phosphate to (S)-serinol phosphate in the presence of l-glutamate. The probable routes to the ultimate formation of the glucuronamide substructures in the capsular polysaccharides of C. jejuni are discussed.

Building better polymerases: Engineering the replication of expanded genetic alphabets

DNA polymerases are today used throughout scientific research, biotechnology, and medicine, in part for their ability to interact with unnatural forms of DNA created by synthetic biologists. Here especially, natural DNA polymerases often do not have the "performance specifications" needed for transformative technologies. This creates a need for science-guided rational (or semi-rational) engineering to identify variants that replicate unnatural base pairs (UBPs), unnatural backbones, tags, or other evolutionarily novel features of unnatural DNA. In this review, we provide a brief overview of the chemistry and properties of replicative DNA polymerases and their evolved variants, focusing on the Klenow fragment of Taq DNA polymerase (Klentaq). We describe comparative structural, enzymatic, and molecular dynamics studies of WT and Klentaq variants, complexed with natural or noncanonical substrates. Combining these methods provides insight into how specific amino acid substitutions distant from the active site in a Klentaq DNA polymerase variant (ZP Klentaq) contribute to its ability to replicate UBPs with improved efficiency compared with Klentaq. This approach can therefore serve to guide any future rational engineering of replicative DNA polymerases.

Improving the Treatment of Acute Lymphoblastic Leukemia

l-Asparaginase (EC 3.5.1.1) was first used as a component of combination drug therapies to treat acute lymphoblastic leukemia (ALL), a cancer of the blood and bone marrow, almost 50 years ago. Administering this enzyme to reduce asparagine levels in the blood is a cornerstone of modern clinical protocols for ALL; indeed, this remains the only successful example of a therapy targeted against a specific metabolic weakness in any form of cancer. Three problems, however, constrain the clinical use of l-asparaginase. First, a type II bacterial variant of l-asparaginase is administered to patients, the majority of whom are children, which produces an immune response thereby limiting the time over which the enzyme can be tolerated. Second, l-asparaginase is subject to proteolytic degradation in the blood. Third, toxic side effects are observed, which may be correlated with the l-glutaminase activity of the enzyme. This Perspective will outline how asparagine depletion negatively impacts the growth of leukemic blasts, discuss the structure and mechanism of l-asparaginase, and briefly describe the clinical use of chemically modified forms of clinically useful l-asparaginases, such as Asparlas, which was recently given FDA approval for use in children (babies to young adults) as part of multidrug treatments for ALL. Finally, we review ongoing efforts to engineer l-asparaginase variants with improved therapeutic properties and briefly detail emerging, alternate strategies for the treatment of forms of ALL that are resistant to asparagine depletion.

Uncovering the chemistry of C-C bond formation in C-nucleoside biosynthesis: crystal structure of a C-glycoside synthase/PRPP complex

The enzyme ForT catalyzes C-C bond formation between 5'-phosphoribosyl-1'-pyrophosphate (PRPP) and 4-amino-1H-pyrazole-3,5-dicarboxylate to make a key intermediate in the biosynthesis of formycin A 5'-phosphate by Streptomyces kaniharaensis. We report the 2.5 Å resolution structure of the ForT/PRPP complex and locate active site residues critical for PRPP recognition and catalysis.

Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet

Evolution has yielded biopolymers that are constructed from exactly four building blocks and are able to support Darwinian evolution. Synthetic biology aims to extend this alphabet, and we recently showed that 8-letter (hachimoji) DNA can support rule-based information encoding. One source of replicative error in non-natural DNA-like systems, however, is the occurrence of alternative tautomeric forms, which pair differently. Unfortunately, little is known about how structural modifications impact free-energy differences between tautomers of the non-natural nucleobases used in the hachimoji expanded genetic alphabet. Determining experimental tautomer ratios is technically difficult, and so, strategies for improving hachimoji DNA replication efficiency will benefit from accurate computational predictions of equilibrium tautomeric ratios. We now report that high-level quantum-chemical calculations in aqueous solution by the embedded cluster reference interaction site model, benchmarked against free-energy molecular simulations for solvation thermodynamics, provide useful quantitative information on the tautomer ratios of both Watson-Crick and hachimoji nucleobases. In agreement with previous computational studies, all four Watson-Crick nucleobases adopt essentially only one tautomer in water. This is not the case, however, for non-natural nucleobases and their analogues. For example, although the enols of isoguanine and a series of related purines are not populated in water, these heterocycles possess N₁-H and N₃-H keto tautomers that are similar in energy, thereby adversely impacting accurate nucleobase pairing. These robust computational strategies offer a firm basis for improving experimental measurements of tautomeric ratios, which are currently limited to studying molecules that exist only as two tautomers in solution.

Building better enzymes: Molecular basis of improved non-natural nucleobase incorporation by an evolved DNA polymerase

Obtaining semisynthetic microorganisms that exploit the information density of "hachimoji" DNA requires access to engineered DNA polymerases. A KlenTaq variant has been reported that incorporates the "hachimoji" P:Z nucleobase pair with a similar efficiency to that seen for Watson-Crick nucleobase incorporation by the wild type (WT) KlenTaq DNA polymerase. The variant polymerase differs from WT KlenTaq by only four amino acid substitutions, none of which are located within the active site. We now report molecular dynamics (MD) simulations on a series of binary complexes aimed at elucidating the contributions of the four amino acid substitutions to altered catalytic activity. These simulations suggest that WT KlenTaq is insufficiently flexible to be able to bind AEGIS DNA correctly, leading to the loss of key protein/DNA interactions needed to position the binary complex for efficient incorporation of the "hachimoji" Z nucleobase. In addition, we test literature hypotheses about the functional roles of each amino acid substitution and provide a molecular description of how individual residue changes contribute to the improved activity of the KlenTaq variant. We demonstrate that MD simulations have a clear role to play in systematically screening DNA polymerase variants capable of incorporating different types of nonnatural nucleobases thereby limiting the number that need to be characterized by experiment. It is now possible to build DNA molecules containing nonnatural nucleobase pairs in addition to A:T and G:C. Exploiting this development in synthetic biology requires engineered DNA polymerases that can replicate nonnatural nucleobase pairs. Computational studies on a DNA polymerase variant reveal how amino acid substitutions outside of the active site yield an enzyme that replicates nonnatural nucleobase pairs with high efficiency. This work will facilitate efforts to obtain bacteria possessing an expanded genetic alphabet.

PMP-diketopiperazine adducts form at the active site of a PLP dependent enzyme involved in formycin biosynthesis

ForI is a PLP-dependent enzyme from the biosynthetic pathway of the C-nucleoside antibiotic formycin. Cycloserine is thought to inhibit PLP-dependent enzymes by irreversibly forming a PMP-isoxazole. We now report that ForI forms novel PMP-diketopiperazine derivatives following incubation with both d and l cycloserine. This unexpected result suggests chemical diversity in the chemistry of cycloserine inhibition.

High-resolution crystal structure of human asparagine synthetase enables analysis of inhibitor binding and selectivity

Expression of human asparagine synthetase (ASNS) promotes metastatic progression and tumor cell invasiveness in colorectal and breast cancer, presumably by altering cellular levels of L-asparagine. Human ASNS is therefore emerging as a bona fide drug target for cancer therapy. Here we show that a slow-onset, tight binding inhibitor, which exhibits nanomolar affinity for human ASNS in vitro, exhibits excellent selectivity at 10 μM concentration in HCT-116 cell lysates with almost no off-target binding. The high-resolution (1.85 Å) crystal structure of human ASNS has enabled us to identify a cluster of negatively charged side chains in the synthetase domain that plays a key role in inhibitor binding. Comparing this structure with those of evolutionarily related AMP-forming enzymes provides insights into intermolecular interactions that give rise to the observed binding selectivity. Our findings demonstrate the feasibility of developing second generation human ASNS inhibitors as lead compounds for the discovery of drugs against metastasis.

Wen Zhu et al. report the crystal structure of human asparagine synthetase at a 1.85 Å resolution, enabling computational analysis of inhibitor binding. They also find new insights into the intermolecular interactions contributing to binding specificity of inhibitors.

Assessing the Influence of Mutation on GTPase Transition States by Using X-ray Crystallography, 19 F NMR, and DFT Approaches

We report X-ray crystallographic and ¹⁹ F NMR studies of the G-protein RhoA complexed with MgF₃^- , GDP, and RhoGAP, which has the mutation Arg85'Ala. When combined with DFT calculations, these data permit the identification of changes in transition state (TS) properties. The X-ray data show how Tyr34 maintains solvent exclusion and the core H-bond network in the active site by relocating to replace the missing Arg85' sidechain. The ¹⁹ F NMR data show deshielding effects that indicate the main function of Arg85' is electronic polarization of the transferring phosphoryl group, primarily mediated by H-bonding to O^3G and thence to P^G . DFT calculations identify electron-density redistribution and pinpoint why the TS for guanosine 5'-triphosphate (GTP) hydrolysis is higher in energy when RhoA is complexed with RhoGAP_Arg85'Ala relative to wild-type (WT) RhoGAP. This study demonstrates that ¹⁹ F NMR measurements, in combination with X-ray crystallography and DFT calculations, can reliably dissect the response of small GTPases to site-specific modifications.

Toward an Expanded Genome: Structural and Computational Characterization of an Artificially Expanded Genetic Information System

Although the fundamental properties of DNA as first proposed by Watson and Crick in 1953 provided a basic understanding of how duplex DNA was organized and might be replicated, it was not until the first crystal structures of DNA (Z-DNA in 1979, B-DNA in 1980, and A-DNA in 1982) that the true complexity of the molecule began to be appreciated. Many crystal structures of oligonucleotides have since shed light on the helical forms that "Watson-Crick" DNA can adopt, their associated groove widths, and the properties of the nucleobase pairs and their interactions in all three helical forms. Additional understanding of the properties of Watson-Crick DNA has been provided by computational studies employing a variety of theoretical methods. Together with these studies devoted to understanding Watson-Crick DNA, recent efforts to expand the genetic alphabet have founded a new field in synthetic biology. One of these efforts, the artificially expanded genetic information system (AEGIS) developed by Steven Benner and co-workers, takes advantage of orthogonal hydrogen bonding to produce DNA comprised of six nucleobase pairs, of which the most extensively studied is referred to as P:Z with P being 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) and Z being 6-amino-5-nitro-2(1H)-pyridone. P:Z forms three edge-on hydrogen bonds that differ from standard Watson-Crick pairs in the arrangement of acceptors and donor groups; P presents acceptor, acceptor, donor, and Z presents donor, donor, acceptor. Z is unique among the AEGIS nucleobases in having a nitro group present in the major groove. PZ-containing DNA has been exploited in a number of clinical applications and is being used to develop receptors and catalysts. Ultimately, the grand challenge will be to create a semisynthetic organism with an expanded genome. Furthermore, just as our understanding of the properties of natural DNA have benefited from structural and computational characterization, so too will our understanding of artificial DNA. This Account focuses on the structural and biophysical properties of AEGIS DNA containing P:Z pairs. We begin with the fundamental properties of P:Z nucleobase pairs, including their electrostatic potential and hydrogen-bonding energies, as elucidated by quantum mechanical calculations. We then examine the impact of including multiple consecutive P:Z pairs into duplex DNA providing an opportunity to investigate stacking interactions between P:Z pairs. The self-complementary 5'-CTTATPPTAZZATAAG was crystallized in B-form using the host-guest system along with analogous natural sequences including Gs or As. Use of the host-guest system to characterize B-DNA obviates a number of limitations on the structural characterization of sequences of interest; these include the ability to crystallize the desired sequences and to distinguish structural effects imparted by the lattice constraints from those inherent in the sequence itself. On the other hand, 3/6ZP, 5'-CTTATPPPZZZATAAG, was crystallized in A-form in a DNA-only lattice allowing a comparative analysis of P:Z pairs in two of the biologically relevant helical forms: A- and B-DNA. Computational studies on the 3/6ZP sequence starting in A-form provide additional evidence for a more energetically favorable stacking interaction, which we term the "slide" conformer, observed in the A-form crystal structure; this unusual stacking interaction plays a major role in altering the conformational dynamics observed for the PZ-containing duplex as compared to a GC-containing "control" duplex in long time scale molecular dynamics simulations. This combined use of structural and computational strategies paves the way for obtaining a detailed description of artificial DNA, both in how it differs from Watson-Crick DNA and in the rational discovery of proteins, such as endonucleases, transcription factors, and polymerases, which can specifically manipulate DNA containing AEGIS nucleobase pairs.

How to name atoms in phosphates, polyphosphates, their derivatives and mimics, and transition state analogues for enzyme-catalysed phosphoryl transfer reactions (IUPAC Recommendations 2016)

Procedures are proposed for the naming of individual atoms, P, O, F, N, and S in phosphate esters, amidates, thiophosphates, polyphosphates, their mimics, and analogues of transition states for enzyme-catalyzed phosphoryl transfer reactions. Their purpose is to enable scientists in very different fields, e.g. biochemistry, biophysics, chemistry, computational chemistry, crystallography, and molecular biology, to share standard protocols for the labelling of individual atoms in complex molecules. This will facilitate clear and unambiguous descriptions of structural results, as well as scientific intercommunication concerning them. At the present time, perusal of the Protein Data Bank (PDB) and other sources shows that there is a limited degree of commonality in nomenclature, but a large measure of irregularity in more complex structures. The recommendations described here adhere to established practice as closely as possible, in particular to IUPAC and IUBMB recommendations and to “best practice” in the PDB, especially to its atom labelling of amino acids, and particularly to Cahn-Ingold-Prelog rules for stereochemical nomenclature. They are designed to work in complex enzyme sites for binding phosphates but also to have utility for non-enzymatic systems. Above all, the recommendations are designed to be easy to comprehend and user-friendly.

Synthesis and anti-inflammatory activity evaluation of novel triazolyl-isatin hybrids

New isatin-triazole based hybrids have been synthesized and evaluated for their inhibitory activity of TNF-α induced expression of Intercellular Adhesion Molecule-1 (ICAM-1) on the surface of human endothelial cells. Structure-activity relationship (SAR) studies revealed that the presence of the electron-attracting bromo substituent at position-5 of the isatin moiety played an important role in enhancing the anti-inflammatory potential of the synthesized compounds. Z-1-[3-(1H-1,2,4-Triazol-1-yl)propyl]-5-bromo-3-[2-(4-methoxyphenyl)hydrazono]indolin-2-one (19) with an IC50 = 20 μM and 89% ICAM-1 inhibition with MTD at 200 μM was found to be the most potent of all the synthesized derivatives. Introduction of 1,2,4-triazole ring and electron-donating methoxy group on the phenylhydrazone moiety resulted in four-fold increase of the anti-inflammatory activity.

Substrate Binding Mode and Molecular Basis of a Specificity Switch in Oxalate Decarboxylase

Oxalate decarboxylase (OxDC) catalyzes the conversion of oxalate into formate and carbon dioxide in a remarkable reaction that requires manganese and dioxygen. Previous studies have shown that replacing an active-site loop segment Ser¹⁶¹-Glu¹⁶²-Asn¹⁶³-Ser¹⁶⁴ in the N-terminal domain of OxDC with the cognate residues Asp¹⁶¹-Ala¹⁶²-Ser-¹⁶³-Asn¹⁶⁴ of an evolutionarily related, Mn-dependent oxalate oxidase gives a chimeric variant (DASN) that exhibits significantly increased oxidase activity. The mechanistic basis for this change in activity has now been investigated using membrane inlet mass spectrometry (MIMS) and isotope effect (IE) measurements. Quantitative analysis of the reaction stoichiometry as a function of oxalate concentration, as determined by MIMS, suggests that the increased oxidase activity of the DASN OxDC variant is associated with only a small fraction of the enzyme molecules in solution. In addition, IE measurements show that C–C bond cleavage in the DASN OxDC variant proceeds via the same mechanism as in the wild-type enzyme, even though the Glu¹⁶² side chain is absent. Thus, replacement of the loop residues does not modulate the chemistry of the enzyme-bound Mn(II) ion. Taken together, these results raise the possibility that the observed oxidase activity of the DASN OxDC variant arises from an increased level of access of the solvent to the active site during catalysis, implying that the functional role of Glu¹⁶² is to control loop conformation. A 2.6 Å resolution X-ray crystal structure of a complex between oxalate and the Co(II)-substituted ΔE162 OxDC variant, in which Glu¹⁶² has been deleted from the active site loop, reveals the likely mode by which the substrate coordinates the catalytically active Mn ion prior to C–C bond cleavage. The “end-on” conformation of oxalate observed in the structure is consistent with the previously published V/K IE data and provides an empty coordination site for the dioxygen ligand that is thought to mediate the formation of Mn(III) for catalysis upon substrate binding.

(19)F NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA-Catalyzed GTP Hydrolysis

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use (19)F NMR chemical shifts in the MgF3(-) transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ-PO3(-) oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ-PO3(-) group based on the structure of a RhoA/GAP-GDP-MgF3(-) TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85' side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln-63 carboxamide by Arg85' to preclude direct hydrogen bonding from water to the target γ-PO3(-) group.

Substrate Distortion and the Catalytic Reaction Mechanism of 5-Carboxyvanillate Decarboxylase

5-Carboxyvanillate decarboxylase (LigW) catalyzes the conversion of 5-carboxyvanillate to vanillate in the biochemical pathway for the degradation of lignin. This enzyme was shown to require Mn²⁺ for catalytic activity and the kinetic constants for the decarboxylation of 5-carboxyvanillate by the enzymes from Sphingomonas paucimobilis SYK-6 (k_cat = 2.2 s^–1 and k_cat/K_m = 4.0 × 10⁴ M^–1 s^–1) and Novosphingobium aromaticivorans (k_cat = 27 s^–1 and k_cat/K_m = 1.1 × 10⁵ M^–1 s^–1) were determined. The three-dimensional structures of both enzymes were determined in the presence and absence of ligands bound in the active site. The structure of LigW from N. aromaticivorans, bound with the substrate analogue, 5-nitrovanillate (K_d = 5.0 nM), was determined to a resolution of 1.07 Å. The structure of this complex shows a remarkable enzyme-induced distortion of the nitro-substituent out of the plane of the phenyl ring by approximately 23°. A chemical reaction mechanism for the decarboxylation of 5-carboxyvanillate by LigW was proposed on the basis of the high resolution X-ray structures determined in the presence ligands bound in the active site, mutation of active site residues, and the magnitude of the product isotope effect determined in a mixture of H₂O and D₂O. In the proposed reaction mechanism the enzyme facilitates the transfer of a proton to C5 of the substrate prior to the decarboxylation step.

Formation of Hexacoordinate Mn(III) in Bacillus subtilis Oxalate Decarboxylase Requires Catalytic Turnover

Oxalate decarboxylase (OxDC) catalyzes the disproportionation of oxalic acid monoanion into CO2 and formate. The enzyme has long been hypothesized to utilize dioxygen to form mononuclear Mn(III) or Mn(IV) in the catalytic site during turnover. Recombinant OxDC, however, contains only tightly bound Mn(II), and direct spectroscopic detection of the metal in higher oxidation states under optimal catalytic conditions (pH 4.2) has not yet been reported. Using parallel mode electron paramagnetic resonance spectroscopy, we now show that substantial amounts of Mn(III) are indeed formed in OxDC, but only in the presence of oxalate and dioxygen under acidic conditions. These observations provide the first direct support for proposals in which Mn(III) removes an electron from the substrate to yield a radical intermediate in which the barrier to C-C bond cleavage is significantly decreased. Thus, OxDC joins a small list of enzymes capable of stabilizing and controlling the reactivity of the powerful oxidizing species Mn(III).

Synthesis of macromolecular systems via lipase catalyzed biocatalytic reactions: applications and future perspectives

This review highlights the application of lipases in the synthesis of pharmaceutically important small molecules and polymers for diverse applications.

Gas phase RDX decomposition pathways using coupled cluster theory

Electronic and free energy barriers for a series of gas-phase RDX decomposition mechanisms have been obtain using coupled cluster singles, doubles, and perturbative triples with complete basis set (CCSD(T)/CBS) electronic energies for MBPT(2)/cc-pVTZ structures.

Intrinsic Carbon-Carbon Bond Reactivity at the Manganese Center of Oxalate Decarboxylase from Density Functional Theory

The detailed manganese-dependent chemistry employed by oxalate decarboxylase (OxDC) to catalyze the nonoxidative decarboxylation of oxalic acid remains poorly understood. For example, enzyme activity requires the presence of dioxygen even though this compound is not a formal substrate in the reaction. We now report density functional theory (DFT) calculations upon a series of hypothetical OxDC active site model structures. Our results suggest that the function of the metal ion may be to position dioxygen and oxalate such that electrons can be shuttled directly between these species, thereby removing the need for the existence of Mn(III) as an intermediate in the mechanism. These calculations also indicate that the intrinsic, gas-phase reactivity of the Bacillus subtilis oxalate decarboxylase active center is to oxidize oxalate. Since this reactivity is not observed for OxDC, our DFT results suggest that protein environment modulates the intrinsic metallocenter reactivity, presumably by affecting the electronic distribution at the manganese center during catalysis.

Structural basis for a six nucleotide genetic alphabet

Expanded genetic systems are most likely to work with natural enzymes if the added nucleotides pair with geometries that are similar to those displayed by standard duplex DNA. Here, we present crystal structures of 16-mer duplexes showing this to be the case with two nonstandard nucleobases (Z, 6-amino-5-nitro-2(1H)-pyridone and P, 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)one) that were designed to form a Z:P pair with a standard "edge on" Watson-Crick geometry, but joined by rearranged hydrogen bond donor and acceptor groups. One duplex, with four Z:P pairs, was crystallized with a reverse transcriptase host and adopts primarily a B-form. Another contained six consecutive Z:P pairs; it crystallized without a host in an A-form. In both structures, Z:P pairs fit canonical nucleobase hydrogen-bonding parameters and known DNA helical forms. Unique features include stacking of the nitro group on Z with the adjacent nucleobase ring in the A-form duplex. In both B- and A-duplexes, major groove widths for the Z:P pairs are approximately 1 Å wider than those of comparable G:C pairs, perhaps to accommodate the large nitro group on Z. Otherwise, ZP-rich DNA had many of the same properties as CG-rich DNA, a conclusion supported by circular dichroism studies in solution. The ability of standard duplexes to accommodate multiple and consecutive Z:P pairs is consistent with the ability of natural polymerases to biosynthesize those pairs. This, in turn, implies that the GACTZP synthetic genetic system can explore the entire expanded sequence space that additional nucleotides create, a major step forward in this area of synthetic biology.

Facile C(sp(2))-C(sp(2)) bond cleavage in oxalic acid-derived radicals

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. Many questions remain about the catalytic mechanism of OxDC although it has been proposed that the reaction proceeds via substrate-based radical intermediates. Using coupled cluster theory combined with implicit solvation models we have examined the effects of radical formation on the structure and reactivity of oxalic acid-derived radicals in aqueous solution. Our results show that the calculated solution-phase free-energy barrier for C-C bond cleavage to form CO2 is decreased from 34.2 kcal/mol for oxalic acid to only 9.3 kcal/mol and a maximum of 3.5 kcal/mol for the cationic and neutral oxalic acid-derived radicals, respectively. These studies also show that the C-C σ bonding orbital of the radical cation contains only a single electron, giving rise to an elongated C-C bond distance of 1.7 Å; a similar lengthening of the C-C bond is not observed for the neutral radical. This study provides new chemical insights into the structure and stability of plausible intermediates in the catalytic mechanism of OxDC, and suggests that removal of an electron to form a radical (with or without the concomitant loss of a proton) may be a general strategy for cleaving the unreactive C-C bonds between adjacent sp(2)-hybridized carbon atoms.

Observation of superoxide production during catalysis of Bacillus subtilis oxalate decarboxylase at pH 4

This contribution describes the trapping of the hydroperoxyl radical at a pH of 4 during turnover of wild-type oxalate decarboxylase and its T165V mutant using the spin-trap BMPO. Radicals were detected and identified by a combination of EPR and mass spectrometry. Superoxide, or its conjugate acid, the hydroperoxyl radical, is expected as an intermediate in the decarboxylation and oxidation reactions of the oxalate monoanion, both of which are promoted by oxalate decarboxylase. Another intermediate, the carbon dioxide radical anion was also observed. The quantitative yields of superoxide trapping are similar in the wild type and the mutant while it is significantly different for the trapping of the carbon dioxide radical anion. This suggests that the two radicals are released from different sites of the protein.

A mechanochemical switch to control radical intermediates

B₁₂-dependent enzymes employ radical species with exceptional prowess to catalyze some of the most chemically challenging, thermodynamically unfavorable reactions. However, dealing with highly reactive intermediates is an extremely demanding task, requiring sophisticated control strategies to prevent unwanted side reactions. Using hybrid quantum mechanical/molecular mechanical simulations, we follow the full catalytic cycle of an AdoB₁₂-dependent enzyme and present the details of a mechanism that utilizes a highly effective mechanochemical switch. When the switch is "off", the 5'-deoxyadenosyl radical moiety is stabilized by releasing the internal strain of an enzyme-imposed conformation. Turning the switch "on," the enzyme environment becomes the driving force to impose a distinct conformation of the 5'-deoxyadenosyl radical to avoid deleterious radical transfer. This mechanochemical switch illustrates the elaborate way in which enzymes attain selectivity of extremely chemically challenging reactions.

Assigning the EPR fine structure parameters of the Mn(II) centers in Bacillus subtilis oxalate decarboxylase by site-directed mutagenesis and DFT/MM calculations

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. EPR-based strategies for investigating the catalytic mechanism of decarboxylation are complicated by the difficulty of assigning the signals associated with the two Mn(II) centers located in the N- and C-terminal cupin domains of the enzyme. We now report a mutational strategy that has established the assignment of EPR fine structure parameters to each of these Mn(II) centers at pH 8.5. These experimental findings are also used to assess the performance of a multistep strategy for calculating the zero-field splitting parameters of protein-bound Mn(II) ions. Despite the known sensitivity of calculated D and E values to the computational approach, we demonstrate that good estimates of these parameters can be obtained using cluster models taken from carefully optimized DFT/MM structures. Overall, our results provide new insights into the strengths and limitations of theoretical methods for understanding electronic properties of protein-bound Mn(II) ions, thereby setting the stage for future EPR studies on the electronic properties of the Mn(II) centers in OxDC and site-specific variants.

Rhodium-catalyzed direct synthesis of unprotected NH-sulfoximines from sulfoxides

A novel rhodium–catalyzed imination of sulfoxides using O-(2,4-dinitrophenyl)hydroxylamine is developed under mild conditions with good functional group tolerance.

Computational, structural, and kinetic evidence that Vibrio vulnificus FrsA is not a cofactor-independent pyruvate decarboxylase

The fermentation-respiration switch (FrsA) protein in Vibrio vulnificus was recently reported to catalyze the cofactor-independent decarboxylation of pyruvate. We now report quantum mechanical/molecular mechenical calculations that examine the energetics of C-C bond cleavage for a pyruvate molecule bound within the putative active site of FrsA. These calculations suggest that the barrier to C-C bond cleavage in the bound substrate is 28 kcal/mol, which is similar to that estimated for the uncatalyzed decarboxylation of pyruvate in water at 25 °C. In agreement with the theoretical predictions, no pyruvate decarboxylase activity was detected for recombinant FrsA protein that could be crystallized and structurally characterized. These results suggest that the functional annotation of FrsA as a cofactor-independent pyruvate decarboxylase is incorrect.

A sulfoximine-based inhibitor of human asparagine synthetase kills L-asparaginase-resistant leukemia cells

An adenylated sulfoximine transition-state analogue 1, which inhibits human asparagine synthetase (hASNS) with nanomolar potency, has been reported to suppress the proliferation of an l-asparagine amidohydrolase (ASNase)-resistant MOLT-4 leukemia cell line (MOLT-4R) when l-asparagine is depleted in the medium. We now report the synthesis and biological activity of two new sulfoximine analogues of 1 that have been studied as part of systematic efforts to identify compounds with improved cell permeability and/or metabolic stability. One of these new analogues, an amino sulfoximine 5 having no net charge at cellular pH, is a better hASNS inhibitor (K(I)(∗)=8 nM) than 1 and suppresses proliferation of MOLT-4R cells at 10-fold lower concentration (IC(50)=0.1mM). More importantly, and in contrast to the lead compound 1, the presence of sulfoximine 5 at concentrations above 0.25 mM causes the death of MOLT-4R cells even when ASNase is absent in the culture medium. The amino sulfoximine 5 exhibits different dose-response behavior when incubated with an ASNase-sensitive MOLT-4 cell line (MOLT-4S), supporting the hypothesis that sulfoximine 5 exerts its effect by inhibiting hASNS in the cell. Our work provides further evidence for the idea that hASNS represents a chemotherapeutic target for the treatment of leukemia, and perhaps other cancers, including those of the prostate.

A structural element that facilitates proton-coupled electron transfer in oxalate decarboxylase

The conformational properties of an active-site loop segment, defined by residues Ser(161)-Glu(162)-Asn(163)-Ser(164), have been shown to be important for modulating the intrinsic reactivity of Mn(II) in the active site of Bacillus subtilis oxalate decarboxylase. We now detail the functional and structural consequences of removing a conserved Arg/Thr hydrogen-bonding interaction by site-specific mutagenesis. Hence, substitution of Thr-165 by a valine residue gives an OxDC variant (T165V) that exhibits impaired catalytic activity. Heavy-atom isotope effect measurements, in combination with the X-ray crystal structure of the T165V OxDC variant, demonstrate that the conserved Arg/Thr hydrogen bond is important for correctly locating the side chain of Glu-162, which mediates a proton-coupled electron transfer (PCET) step prior to decarboxylation in the catalytically competent form of OxDC. In addition, we show that the T165V OxDC variant exhibits a lower level of oxalate consumption per dioxygen molecule, consistent with the predictions of recent spin-trapping experiments [Imaram et al. (2011) Free Radicals Biol. Med. 50, 1009-1015]. This finding implies that dioxygen might participate as a reversible electron sink in two putative PCET steps and is not merely used to generate a protein-based radical or oxidized metal center.

Membrane inlet for mass spectrometric measurement of catalysis by enzymatic decarboxylases

Membrane inlet mass spectrometry (MIMS) uses diffusion across a permeable membrane to detect in solution uncharged molecules of small molecular weight. We point out here the application of MIMS to determine catalytic properties of decarboxylases using as an example catalysis by oxalate decarboxylase (OxDC) from Bacillus subtilis. The decarboxylase activity generates carbon dioxide and formate from the nonoxidative reaction but is accompanied by a concomitant oxidase activity that consumes oxalate and oxygen and generates CO(2) and hydrogen peroxide. The application of MIMS in measuring catalysis by OxDC involves the real-time and continuous detection of oxygen and product CO(2) from the ion currents of their respective mass peaks. Steady-state catalytic constants for the decarboxylase activity obtained by measuring product CO(2) using MIMS are comparable to those acquired by the traditional endpoint assay based on the coupled reaction with formate dehydrogenase, and measuring consumption of O(2) using MIMS also estimates the oxidase activity. The use of isotope-labeled substrate ((13)C(2)-enriched oxalate) in MIMS provides a method to characterize the catalytic reaction in cell suspensions by detecting the mass peak for product (13)CO(2) (m/z 45), avoiding inaccuracies due to endogenous (12)CO(2).

Characterization of Ceriporiopsis subvermispora bicupin oxalate oxidase expressed in Pichia pastoris

Oxalate oxidase (E.C. 1.2.3.4) catalyzes the oxygen-dependent oxidation of oxalate to carbon dioxide in a reaction that is coupled with the formation of hydrogen peroxide. Although there is currently no structural information available for oxalate oxidase from Ceriporiopsis subvermispora (CsOxOx), sequence data and homology modeling indicate that it is the first manganese-containing bicupin enzyme identified that catalyzes this reaction. Interestingly, CsOxOx shares greatest sequence homology with bicupin microbial oxalate decarboxylases (OxDC). We show that CsOxOx activity directly correlates with Mn content and other metals do not appear to be able to support catalysis. EPR spectra indicate that the Mn is present as Mn(II), and are consistent with the coordination environment expected from homology modeling with known X-ray crystal structures of OxDC from Bacillus subtilis. EPR spin-trapping experiments support the existence of an oxalate-derived radical species formed during turnover. Acetate and a number of other small molecule carboxylic acids are competitive inhibitors for oxalate in the CsOxOx catalyzed reaction. The pH dependence of this reaction suggests that the dominant contribution to catalysis comes from the monoprotonated form of oxalate binding to a form of the enzyme in which an active site carboxylic acid residue must be unprotonated.

EPR spin trapping of an oxalate-derived free radical in the oxalate decarboxylase reaction

EPR spin trapping experiments on bacterial oxalate decarboxylase from Bacillus subtilis under turn-over conditions are described. The use of doubly (13)C-labeled oxalate leads to a characteristic splitting of the observed radical adducts using the spin trap N-tert-butyl-α-phenylnitrone linking them directly to the substrate. The radical was identified as the carbon dioxide radical anion which is a key intermediate in the hypothetical reaction mechanism of both decarboxylase and oxidase activities. X-ray crystallography had identified a flexible loop, SENS161-4, which acts as a lid to the putative active site. Site directed mutagenesis of the hinge amino acids, S161 and T165 was explored and showed increased radical trapping yields compared to the wild type. In particular, T165V shows approximately ten times higher radical yields while at the same time its decarboxylase activity was reduced by about a factor of ten. This mutant lacks a critical H-bond between T165 and R92 resulting in compromised control over its radical chemistry allowing the radical intermediate to leak into the surrounding solution.

Nitric oxide reversibly inhibits Bacillus subtilis oxalate decarboxylase

Membrane inlet mass spectrometry (MIMS) has been employed to assay the catalytic activity of oxalate decarboxylase (OxDC), allowing us to demonstrate that nitric oxide (NO) reversibly inhibits the enzyme under dioxygen-depleted conditions. X-band EPR measurements do not provide any direct evidence for the interaction of NO with either of the Mn(II) centers in OxDC raising the possibility that there is a separate dioxygen-binding pocket in the enzyme.

A conserved glutamate controls the commitment to acyl-adenylate formation in asparagine synthetase

Inhibitor docking studies have implicated a conserved glutamate residue (Glu-348) as a general base in the synthetase active site of the enzyme asparagine synthetase B from Escherichia coli (AS-B). We now report steady-state kinetic, isotope transfer, and positional isotope exchange experiments for a series of site-directed AS-B mutants in which Glu-348 is substituted by conservative amino acid replacements. We find that formation of the β-aspartyl-AMP intermediate, and therefore the eventual production of asparagine, is dependent on the presence of a carboxylate side chain at this position in the synthetase active site. In addition, Glu-348 may also play a role in mediating the conformational changes needed to (i) coordinate, albeit weakly, the glutaminase and synthetase activities of the enzyme and (ii) establish the structural integrity of the intramolecular tunnel along which ammonia is translocated. The importance of Glu-348 in mediating acyl-adenylate formation contrasts with the functional role of the cognate residues in β-lactam synthetase (BLS) and carbapenem synthetase (CPS) (Tyr-348 and Tyr-345, respectively), which both likely evolved from asparagine synthetase. Given the similarity of the chemistry catalyzed by AS-B, BLS, and CPS, our work highlights the difficulty of predicting the functional outcome of single site mutations on enzymes that catalyze almost identical chemical transformations.

Sampling long time scale protein motions: OSRW simulation of active site loop conformational free energies in formyl-CoA:oxalate CoA transferase

X-ray crystallographic snapshots have shown that conformational changes of a tetraglycine loop in the active site of formyl-CoA:oxalate CoA transferase (FRC) play an important role in the catalytic cycle of the enzyme. Orthogonal space random walk (OSRW) simulations have been applied to obtain quantitative computational estimates of the relative free energy of the "open" and "closed" conformations of this loop together with the energetic barrier for interconversion of these states in wild type FRC. These OSRW calculations not only show that the two conformations have similar free energies but also predict a barrier that is consistent with the observed turnover number of the enzyme. In an effort to quantitate the importance of specific residues in the tetraglycine loop, OSRW simulations have also been performed on the G258A, G259A, G260A, and G261A FRC variants both to examine the energetic effects of replacing each glycine residue and to correlate the computed energies with kinetic and structural observations. In enzymes with substantially reduced catalytic efficiency (k(cat)/K(M)), the OSRW simulations reveal the adoption of additional low energy loop conformations. In the case of the G260A FRC variant, the new conformation identified by simulation is similar to that observed in the X-ray crystal structure of the protein. These results provide further evidence for the power of the OSRW method in sampling conformational space and, hence, in providing quantitative free energy estimates for the conformations adopted by functionally important active site loops. In addition, these simulations model the motions of side chains that are correlated with changes in loop conformation thereby permitting access to long time-scale motions through the use of nanosecond simulations.

Metal dependence of oxalate decarboxylase activity

Bacillus subtilis oxalate decarboxylase (OxDC) catalyzes the conversion of oxalate into CO(2) and formate. The enzyme is composed of two cupin domains, each of which contains a Mn(II) ion. Although there is general agreement that Mn(II) in the N-terminal domain mediates OxDC-catalyzed decarboxylation, legitimate questions have been raised concerning the function (if any) of the Mn(II) bound in the C-terminal cupin domain. We have investigated this problem using a series of OxDC mutants in which Mn(II) binding is perturbed by mutagenesis of Glu-101 and Glu-280, which coordinate the metal in the N-terminal and C-terminal domains, respectively. We now demonstrate that decarboxylase activity and total manganese content are sensitive to modifications in either metal-binding glutamate residue. These findings, in combination with EPR measurements, raise the possibility that the C-terminal Mn(II) center can catalyze the decarboxylation reaction. Further support for this conclusion has been provided from a combination of in vivo and in vitro strategies for preparing wild-type OxDC in which Mn(II) is incorporated to a variety of extents. Kinetic characterization of these variants shows that OxDC activity is linearly correlated with manganese content, as might be expected if both sites can catalyze the breakdown of oxalate into formate and CO(2). These studies also represent the first unequivocal demonstration that OxDC activity is uniquely mediated by manganese.

Shining a light on post-translational modification

Post-translational modification, such as phosphorylation or glycosylation, provides a mechanism for increasing the diversity of protein structures in the cell and regulating biological activity. In addition, such modifications may result in the localization of proteins to specific cellular organelles, with incorrect targeting being associated with a number of diseases. The simplest strategy to identify the functional importance of post-translational modifications is to use mutagenesis methods to replace the residue that is post-translationally modified by one that cannot undergo the relevant chemical transformation. Merely causing "loss of function" does not, however, address questions concerning how cellular function depends on the timing of post-translational changes andor the movement of modified proteins between organelles. The recent demonstration that genetically encoded "photocaged" proteins can be employed to resolve such issues therefore represents an exciting advance in this research area, and is an elegant illustration of the power of combining the power of chemical synthesis and methods for manipulating the biological machinery of protein synthesis.

Differential substrate specificity and kinetic behavior of Escherichia coli YfdW and Oxalobacter formigenes formyl coenzyme A transferase

The yfdXWUVE operon appears to encode proteins that enhance the ability of Escherichia coli MG1655 to survive under acidic conditions. Although the molecular mechanisms underlying this phenotypic behavior remain to be elucidated, findings from structural genomic studies have shown that the structure of YfdW, the protein encoded by the yfdW gene, is homologous to that of the enzyme that mediates oxalate catabolism in the obligate anaerobe Oxalobacter formigenes, O. formigenes formyl coenzyme A transferase (FRC). We now report the first detailed examination of the steady-state kinetic behavior and substrate specificity of recombinant, wild-type YfdW. Our studies confirm that YfdW is a formyl coenzyme A (formyl-CoA) transferase, and YfdW appears to be more stringent than the corresponding enzyme (FRC) in Oxalobacter in employing formyl-CoA and oxalate as substrates. We also report the effects of replacing Trp-48 in the FRC active site with the glutamine residue that occupies an equivalent position in the E. coli protein. The results of these experiments show that Trp-48 precludes oxalate binding to a site that mediates substrate inhibition for YfdW. In addition, the replacement of Trp-48 by Gln-48 yields an FRC variant for which oxalate-dependent substrate inhibition is modified to resemble that seen for YfdW. Our findings illustrate the utility of structural homology in assigning enzyme function and raise the question of whether oxalate catabolism takes place in E. coli upon the up-regulation of the yfdXWUVE operon under acidic conditions.

Reinvestigation of the catalytic mechanism of formyl-CoA transferase, a class III CoA-transferase

Formyl-coenzyme A transferase from Oxalobacter formigenes belongs to the Class III coenzyme A transferase family and catalyzes the reversible transfer of a CoA carrier between formyl-CoA and oxalate, forming oxalyl-CoA and formate. Formyl-CoA transferase has a unique three-dimensional fold composed of two interlaced subunits locked together like rings of a chain. We here present an intermediate in the reaction, formyl-CoA transferase containing the covalent beta-aspartyl-CoA thioester, adopting different conformations in the two active sites of the dimer, which was identified through crystallographic freeze-trapping experiments with formyl-CoA and oxalyl-CoA in the absence of acceptor carboxylic acid. The formation of the enzyme-CoA thioester was also confirmed by mass spectrometric data. Further structural data include a trapped aspartyl-formyl anhydride protected by a glycine loop closing down over the active site. In a crystal structure of the beta-aspartyl-CoA thioester of an inactive mutant variant, oxalate was found bound to the open conformation of the glycine loop. Together with hydroxylamine trapping experiments and kinetic as well as mutagenesis data, the structures of these formyl-CoA transferase complexes provide new information on the Class III CoA-transferase family and prompt redefinition of the catalytic steps and the modified reaction mechanism of formyl-CoA transferase proposed here.