Reaction Mechanism of Zinc-Dependent Cytosine Deaminase from Escherichia coli: A Quantum-Chemical Study

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Escherichia coli cytosine deaminase: Structural and biotechnological aspects

Suicide gene therapy involves introducing viral or bacterial genes into tumor cells, which enables the conversion of a nontoxic prodrug into a toxic-lethal drug. The application of the bacterial cytosine deaminase (bCD)/5-fluorocytosine (5-FC) approach has been beneficial and progressive within the current field of cancer therapy because of the enhanced bystander effect. The basis of this method is the preferential deamination of 5-FC to 5-fluorouracil by cancer cells expressing cytosine deaminase (CD), which strongly inhibits DNA synthesis and RNA function, effectively targeting tumor cells. However, the poor binding affinity of toward 5-FC compared to the natural substrate cytosine and/or inappropriate thermostability limits the clinical applications of this gene therapy approach. Nowadays, many genetic engineering studies have been carried out to solve and improve the activity of this enzyme. In the current review, we intend to discuss the biotechnological aspects of Escherichia coli CD, including its structure, functions, molecular cloning, and protein engineering. We will also explore its relevance in cancer clinical trials. By examining these aspects, we hope to provide a thorough understanding of E. coli CD and its potential applications in cancer therapy.

Engineered deaminases as a key component of DNA and RNA editing tools

Over recent years, zinc-dependent deaminases have attracted increasing interest as key components of nucleic acid editing tools that can generate point mutations at specific sites in either DNA or RNA by combining a targeting module (such as a catalytically impaired CRISPR-Cas component) and an effector module (most often a deaminase). Deaminase-based molecular tools are already being utilized in a wide spectrum of therapeutic and research applications; however, their medical and biotechnological potential seems to be much greater. Recent reports indicate that the further development of nucleic acid editing systems depends largely on our ability to engineer the substrate specificity and catalytic activity of the editors themselves. In this review, we summarize the current trends and achievements in deaminase engineering. The presented data indicate that the potential of these enzymes has not yet been fully revealed or understood. Several examples show that even relatively minor changes in the structure of deaminases can give them completely new and unique properties.

Why is manganese so valuable to bacterial pathogens?

Apart from oxygenic photosynthesis, the extent of manganese utilization in bacteria varies from species to species and also appears to depend on external conditions. This observation is in striking contrast to iron, which is similar to manganese but essential for the vast majority of bacteria. To adequately explain the role of manganese in pathogens, we first present in this review that the accumulation of molecular oxygen in the Earth's atmosphere was a key event that linked manganese utilization to iron utilization and put pressure on the use of manganese in general. We devote a large part of our contribution to explanation of how molecular oxygen interferes with iron so that it enhances oxidative stress in cells, and how bacteria have learned to control the concentration of free iron in the cytosol. The functioning of iron in the presence of molecular oxygen serves as a springboard for a fundamental understanding of why manganese is so valued by bacterial pathogens. The bulk of this review addresses how manganese can replace iron in enzymes. Redox-active enzymes must cope with the higher redox potential of manganese compared to iron. Therefore, specific manganese-dependent isoenzymes have evolved that either lower the redox potential of the bound metal or use a stronger oxidant. In contrast, redox-inactive enzymes can exchange the metal directly within the individual active site, so no isoenzymes are required. It appears that in the physiological context, only redox-inactive mononuclear or dinuclear enzymes are capable of replacing iron with manganese within the same active site. In both cases, cytosolic conditions play an important role in the selection of the metal used. In conclusion, we summarize both well-characterized and less-studied mechanisms of the tug-of-war for manganese between host and pathogen.

Role of water coordination at zinc binding site and its catalytic pathway of dizinc creatininase: insights from quantum cluster approach

Creatininase is a key enzyme of creatinine-metabolizing pathway in mammals, and has a great potential for diagnostic application. It catalyzes the reversible conversion of creatinine to creatine. Here, we investigated its reaction mechanism with density functional theory in conjunction with the quantum cluster approach. Three reaction pathways in which several possible proton transfers assisted by either His178 or a water ligand to Zn1 (Wat2) or both were considered. DFT calculations reveal, depending on Wat2 coordination mode at Zn1, two competitive ring-opening pathways where His178 playing a central role as a proton shuttle or both His178 and Wat2 serving as a dual catalytic role as a base and an acid, respectively. Three elementary steps were proposed for the reaction: the first involves nucleophilic attack by a bridging hydroxide to the substrate and forms a gem-diolate intermediate, followed by a proton transfer from the gem-diolate to His178 (His178 protonation is a required step for efficient proton transfers). Finally, the second proton transfer from the protonated His178 or Wat2 to the amide of substrate leads to the ring opening. The first proton transfer is the rate-limiting step of the whole reaction, in consistent with previous experimental and computational studies. A detailed understanding of the reaction mechanism of the creatininase enzyme family will also be helpful for developing a biosensor for kidney function.

Combined Theoretical, Bioinformatic, and Biochemical Analyses of RNA Editing by Adenine Base Editors

Adenine base editors (ABEs) have been subjected to multiple rounds of mutagenesis with the goal of optimizing their function as efficient and precise genome editing agents. Despite an ever-expanding data set of ABE mutants and their corresponding DNA or RNA-editing activity, the molecular mechanisms defining these changes remain to be elucidated. In this study, we provide a systematic interpretation of the nature of these mutations using an entropy-based classification model that relies on evolutionary data from extant protein sequences. Using this model in conjunction with experimental analyses, we identify two previously reported mutations that form an epistatic pair in the RNA-editing functional landscape of ABEs. Molecular dynamics simulations reveal the atomistic details of how these two mutations affect substrate-binding and catalytic activity, via both individual and cooperative effects, hence providing insights into the mechanisms through which these two mutations are epistatically coupled.

The current toolbox for APOBEC drug discovery

Mutational processes driving genome evolution and heterogeneity contribute to immune evasion and therapy resistance in viral infections and cancer. APOBEC3 (A3) enzymes promote such mutations by catalyzing the deamination of cytosines to uracils in single-stranded DNA. Chemical inhibition of A3 enzymes may yield an antimutation therapeutic strategy to improve the durability of current drug therapies that are prone to resistance mutations. A3 small-molecule drug discovery efforts to date have been restricted to a single high-throughput biochemical activity assay; however, the arsenal of discovery assays has significantly expanded in recent years. The assays used to study A3 enzymes are reviewed here with an eye towards their potential for small-molecule discovery efforts.

Insights into the Dual Shuttle Catalytic Mechanism of Guanine Deaminase

Guanine deaminases (GD) are essential enzymes that help in regulating the nucleobase pool. Since the deamination reaction can result in the accumulation of mutagenic bases that can lead to genomic instability, these enzymes are tightly regulated and are nonpromiscuous. Here, we delineate the basis of their substrate fidelity via entailing the reaction mechanism of deamination by employing density functional theory (DFT) calculations on NE0047, a GD from Nitrosomonas europaea. The results show that, unlike pyrimidine deaminases, which require a single glutamic acid as a proton shuttle, GDs involve two amino acids, E79 and E143 (numbering in NE0047), which control its reactivity. The hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have shown that the first Zn-bound proton transfer to the N3 atom of the substrate is mediated by the E79 residue, and the second proton is transferred to the amine nitrogen of substrate via E143. Moreover, cluster models reveal that the crystallographic water molecules near the active site control the reactivity. A comparison with human GD reveals that the proposed catalytic mechanism is generic, and the knowledge generated here can be effectively applied to design selective inhibitors.

Structure guided mutagenesis reveals the substrate determinants of guanine deaminase

Guanine deaminases (GDs) are essential enzymes that regulate the overall nucleobase pool. Since the deamination of guanine to xanthine results in the production of a mutagenic base, these enzymes have evolved to be very specific in nature. Surprisingly, they accept structurally distinct triazine ammeline, an intermediate in the melamine pathway, as one of the moonlighting substrates. Here, by employing NE0047 (a GD from Nitrosomonas europaea), we delineate the nuance in the catalytic mechanism that allows these two distinct substrates to be catalyzed. A combination of enzyme kinetics, X-ray crystallographic, and calorimetric studies reveal that GDs operate via a dual proton shuttle mechanism with two glutamates, E79 and E143, crucial for deamination. Additionally, N66 appears to be central for substrate anchoring and participates in catalysis. The study highlights the importance of closure of the catalytic loop and of maintenance of the hydrophobic core by capping residues like F141 and F48 for the creation of an apt environment for activation of the zinc-assisted catalysis. This study also analyzes evolutionarily distinct GDs and asserts that GDs incorporate subtle variations in the active site architectures while keeping the most critical active site determinants conserved.

Structural Determinants for Substrate Selectivity in Guanine Deaminase Enzymes of the Amidohydrolase Superfamily

Guanine deaminase is a metabolic enzyme, found in all forms of life, which catalyzes the conversion of guanine to xanthine. Despite the availability of several crystal structures, the molecular determinants of substrate orientation and mechanism remain to be elucidated for the amidohydrolase family of guanine deaminase enzymes. Here, we report the crystal structures of Escherichia coli and Saccharomyces cerevisiae guanine deaminase enzymes (EcGuaD and Gud1, respectively), both members of the amidohydrolase superfamily. EcGuaD and Gud1 retain the overall TIM barrel tertiary structure conserved among amidohydrolase enzymes. Both proteins also possess a single zinc cation with trigonal bipyrimidal coordination geometry within their active sites. We also determined a liganded structure of Gud1 bound to the product, xanthine. Analysis of this structure, along with kinetic data of native and site-directed mutants of EcGuaD, identifies several key residues that are responsible for substrate recognition and catalysis. In addition, after a small library of compounds had been screened, two guanine derivatives, 8-azaguanine and 1-methylguanine, were identified as EcGuaD substrates. Interestingly, both EcGuaD and Gud1 also exhibit secondary ammeline deaminase activity. Overall, this work details key structural features of substrate recognition and catalysis of the amidohydrolase family of guanine deaminase enzymes in support of our long-term goal to engineer these enzymes with altered activity and substrate specificity.

Nucleobase deaminases: a potential enzyme system for new therapies

This review presents an overview of the structure, function and mechanism of CDA deaminases and their potential as enzyme systems for development of new antimicrobial therapies.

Measurement and Prediction of Chlorine Kinetic Isotope Effects in Enzymatic Systems

Approaches to determine chlorine kinetic isotope effects (Cl-KIEs) on enzymatic dehalogenations are discussed and illustrated by representative examples. Three aspects are considered. First methodology for experimental measurement of Cl-KIEs, with stress being on FAB-IRMS technique developed in our laboratory, is described. Subsequently, we concentrate our discussion on the consequences of reaction complexity in the interpretation of experimental values, a problem especially important in cases of polychlorinated reactants. The most fruitful studies of enzymatic dehalogenations by Cl-KIEs require their theoretical evaluation, hence the computational focus of the second part of this chapter.

Computational Biochemistry-Enzyme Mechanisms Explored

Understanding enzyme mechanisms is a major task to achieve in order to comprehend how living cells work. Recent advances in biomolecular research provide huge amount of data on enzyme kinetics and structure. The analysis of diverse experimental results and their combination into an overall picture is, however, often challenging. Microscopic details of the enzymatic processes are often anticipated based on several hints from macroscopic experimental data. Computational biochemistry aims at creation of a computational model of an enzyme in order to explain microscopic details of the catalytic process and reproduce or predict macroscopic experimental findings. Results of such computations are in part complementary to experimental data and provide an explanation of a biochemical process at the microscopic level. In order to evaluate the mechanism of an enzyme, a structural model is constructed which can be analyzed by several theoretical approaches. Several simulation methods can and should be combined to get a reliable picture of the process of interest. Furthermore, abstract models of biological systems can be constructed combining computational and experimental data. In this review, we discuss structural computational models of enzymatic systems. We first discuss various models to simulate enzyme catalysis. Furthermore, we review various approaches how to characterize the enzyme mechanism both qualitatively and quantitatively using different modeling approaches.

Product release mechanism and the complete enzyme catalysis cycle in yeast cytosine deaminase (yCD): A computational study

Yeast cytosine deaminase (yCD) is critical in gene-directed enzyme prodrug therapy as it catalyzes the hydrolytic deamination of cytosine. The product (uracil) release process is considered as rate-limiting in the whole enzymatic catalysis and includes the cleavage of the uracil-metal bond and the delivery of free uracil out of the reactive site. Herein extensive combined random acceleration molecular dynamics (RAMD) and molecular dynamics (MD) simulations coupled with the umbrella sampling technique have been performed to study the product transport mechanism. Five channels have been identified, and the thermodynamic and dynamic characterizations for the two most favorable channels have been determined and analyzed. The free energy barrier for the most beneficial pathway is about 13kcal/mol and mainly results from the cleavage of hydrogen bonds between the ligand uracil and surrounding residues Asn51, Glu64, and Asp155. The conjugated rings of Phe114 and Trp152 play gating and guiding roles in the product delivery via π⋯π van der Waals interactions with the product. Finally, the full cycle of the enzymatic catalysis has been determined, making the whole process computationally more precise.

Effect of nucleobase change on cytosine deamination through DNA photo-cross-linking reaction via 3-cyanovinylcarbazole nucleoside

Photo-chemical deamination of cytosine using 3-cyanovinylcarbazole nucleoside (^CNVK) mediated photo-cross-linking is a technique for site-directed mutagenesis. Using this technique in vivo requires the elimination of a high-temperature incubation step; instead, incubation should be carried out under physiological conditions. To improve the reactivity of ^CNVK mediated photo-cross-link induced deamination of cytosine under physiological conditions, an evaluation of base pairing in cytosine was carried out with respect to its deamination. Guanine was replaced with 4 different counter bases (inosine, 2-aminopurine, 5-nitroindole, and nebularine), showing distinct hydrogen bonding patterns with target cytosine, which was incorporated at the -1 position with respect to ^CNVK in the ^CNVK-modified photo-responsive oligodeoxyribonucleotides to ascertain the role of hydrogen bonding in deamination under physiological conditions. Among the counter bases, inosine showed the highest acceleration towards the photo-induced deamination reaction.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

Evaluation of the Substrate Scope of Benzoic Acid (De)carboxylases According to Chemical and Biochemical Parameters

The enzymatic carboxylation of phenolic compounds has been attracting increasing interest in recent years, owing to its regioselectivity and technical potential as a biocatalytic equivalent for the Kolbe-Schmitt reaction. Mechanistically the reaction was demonstrated to occur through electrophilic aromatic substitution/water elimination with bicarbonate as a cosubstrate. The effects of the substituents on the phenolic ring have not yet been elucidated in detail, but this would give detailed insight into the substrate-activity relationship and would provide predictability for the acceptance of future substrates. In this report we show how the kinetic and (apparent) thermodynamic behavior can be explained through the evaluation of linear free energy relationships based on electronic, steric, and geometric parameters and through the consideration of enzyme-ligand interactions. Moreover, the similarity between the benzoic acid decarboxylases and the amidohydrolases superfamily is investigated, and promiscuous hydrolytic activity of the decarboxylase in the context of the hydrolysis of an activated ester bond has been established.

Enzyme catalysis by entropy without Circe effect

Entropic effects have often been invoked to explain the extraordinary catalytic power of enzymes. In particular, the hypothesis that enzymes can use part of the substrate-binding free energy to reduce the entropic penalty associated with the subsequent chemical transformation has been very influential. The enzymatic reaction of cytidine deaminase appears to be a distinct example. Here, substrate binding is associated with a significant entropy loss that closely matches the activation entropy penalty for the uncatalyzed reaction in water, whereas the activation entropy for the rate-limiting catalytic step in the enzyme is close to zero. Herein, we report extensive computer simulations of the cytidine deaminase reaction and its temperature dependence. The energetics of the catalytic reaction is first evaluated by density functional theory calculations. These results are then used to parametrize an empirical valence bond description of the reaction, which allows efficient sampling by molecular dynamics simulations and computation of Arrhenius plots. The thermodynamic activation parameters calculated by this approach are in excellent agreement with experimental data and indeed show an activation entropy close to zero for the rate-limiting transition state. However, the origin of this effect is a change of reaction mechanism compared the uncatalyzed reaction. The enzyme operates by hydroxide ion attack, which is intrinsically associated with a favorable activation entropy. Hence, this has little to do with utilization of binding free energy to pay the entropic penalty but rather reflects how a preorganized active site can stabilize a reaction path that is not operational in solution.

Molecular recognition of thiaclopride by Aplysia californica AChBP: new insights from a computational investigation

The binding of thiaclopride (THI), a neonicotinoid insecticide, with Aplysia californica acetylcholine binding protein (Ac-AChBP), the surrogate of the extracellular domain of insects nicotinic acetylcholine receptors, has been studied with a QM/QM' hybrid methodology using the ONIOM approach (M06-2X/6-311G(d):PM6). The contributions of Ac-AChBP key residues for THI binding are accurately quantified from a structural and energetic point of view. The importance of water mediated hydrogen-bond (H-bond) interactions involving two water molecules and Tyr55 and Ser189 residues in the vicinity of the THI nitrile group, is specially highlighted. A larger stabilization energy is obtained with the THI-Ac-AChBP complex compared to imidacloprid (IMI), the forerunner of neonicotinoid insecticides. Pairwise interaction energy calculations rationalize this result with, in particular, a significantly more important contribution of the pivotal aromatic residues Trp147 and Tyr188 with THI through CH···π/CH···O and π-π stacking interactions, respectively. These trends are confirmed through a complementary non-covalent interaction (NCI) analysis of selected THI-Ac-AChBP amino acid pairs.