Dehydration of ribonucleotides catalyzed by ribonucleotide reductase: the role of the enzyme

From quantum-derived principles underlying cysteine reactivity to combating the COVID-19 pandemic

The COVID‐19 pandemic poses a challenge in coming up with quick and effective means to counter its cause, the SARS‐CoV‐2. Here, we show how the key factors governing cysteine reactivity in proteins derived from combined quantum mechanical/continuum calculations led to a novel multi‐targeting strategy against SARS‐CoV‐2, in contrast to developing potent drugs/vaccines against a single viral target such as the spike protein. Specifically, they led to the discovery of reactive cysteines in evolutionary conserved Zn²⁺‐sites in several SARS‐CoV‐2 proteins that are crucial for viral polypeptide proteolysis as well as viral RNA synthesis, proofreading, and modification. These conserved, reactive cysteines, both free and Zn²⁺‐bound, can be targeted using the same Zn‐ejector drug (disulfiram/ebselen), which enables the use of broad‐spectrum anti‐virals that would otherwise be removed by the virus's proofreading mechanism. Our strategy of targeting multiple, conserved viral proteins that operate at different stages of the virus life cycle using a Zn‐ejector drug combined with other broad‐spectrum anti‐viral drug(s) could enhance the barrier to drug resistance and antiviral effects, as compared to each drug alone. Since these functionally important nonstructural proteins containing reactive cysteines are highly conserved among coronaviruses, our proposed strategy has the potential to tackle future coronaviruses.

This article is categorized under:

Structure and Mechanism > Reaction Mechanisms and Catalysis

Structure and Mechanism > Computational Biochemistry and Biophysics

Electronic Structure Theory > Density Functional Theory

Selective Transformations of Carbohydrates Inspired by Radical-Based Enzymatic Mechanisms

Enzymes are a longstanding source of inspiration for synthetic reaction development. However, enzymatic reactivity and selectivity are frequently untenable in a synthetic context, as the principles that govern control in an enzymatic setting often do not translate to small molecule catalysis. Recent synthetic methods have revealed the viability of using small molecule catalysts to promote highly selective radical-mediated transformations of minimally protected sugar substrates. These transformations share conceptual similarities with radical SAM enzymes found in microbial carbohydrate biosynthesis and present opportunities for synthetic chemists to access microbial and unnatural carbohydrate building blocks without the need for protecting groups or lengthy synthetic sequences. Here, we highlight strategies through which radical reaction pathways can enable the site-, regio-, and diastereoselective transformation of minimally protected carbohydrates in both synthetic and enzymatic systems.

Identifying Drug Repurposing Opportunities in Oncology

The never-ending explosion in the cost of new oncology drugs is reducing in many countries the access to the most recent, effective anticancer therapies and represents a significant obstacle to the design and realization of combinatorial trials. Already approved, anticancer and nonanticancer drugs can be considered for in silico, preclinical, and clinical repurposing approaches and offer the significant advantages of a potentially cheaper, faster, and safer validation. This review discusses recent advances and challenges in the field.

Catalytic Hydrolysis Mechanism of Cocaine by Human Carboxylesterase 1: An Orthoester Intermediate Slows Down the Reaction

Human carboxylesterase 1 (hCES1) is a major carboxylesterase in the human body and plays important roles in the metabolism of a wide variety of substances, including lipids and drugs, and therefore is attracting more and more attention from areas including lipid metabolism, pharmacokinetics, drug–drug interactions, and prodrug activation. In this work, we studied the catalytic hydrolysis mechanism of hCES1 by the quantum mechanics computation method, using cocaine as a model substrate. Our results support the four-step theory of the esterase catalytic hydrolysis mechanism, in which both the acylation stage and the deacylation stage include two transition states and a tetrahedral intermediate. The roles and cooperation of the catalytic triad, S221, H468, and E354, were also analyzed in this study. Moreover, orthoester intermediates were found in hCES1-catalyzed cocaine hydrolysis reaction, which significantly elevate the free energy barrier and slow down the reaction. Based on this finding, we propose that hCES1 substrates with β-aminocarboxylester structure might form orthoester intermediates in hCES1-catalyzed hydrolysis, and therefore prolong their in vivo half-life. Thus, this study helps to clarify the catalytic mechanism of hCES1 and elucidates important details of its catalytic process, and furthermore, provides important insights into the metabolism of hCES1 substrates and drug designing.

Mechanism and Inhibitor Exploration with Binuclear Mg Ketol-Acid Reductoisomerase: Targeting the Biosynthetic Pathway of Branched-Chain Amino Acids

Binuclear Mg ketol-acid reductoisomerase (KARI), which converts (S)-2-acetolactate into (R)-2,3-dihydroxyisovalerate, is responsible for the second step of the biosynthesis of branched-chain amino acids in plants and microorganisms and thus serves as a key inhibition target potentially without effects on mammals. Here, through the use of density functional calculations and a chemical model, the KARI-catalyzed reaction has been demonstrated to include the initial deprotonation of the substrate C2 hydroxy group, bridged by the two Mg ions, alkyl migration from the C2-alkoxide carbon atom to the C3-carbonyl carbon atom, and hydride transfer from a nicotinamide adenine dinucleotide phosphate [NAD(P)H] cofactor to C2. A dead-end mechanism with a hydride transferred to the C3 carbonyl group has been ruled out. The nucleophilicity (migratory aptitude) of the migrating carbon atom and the provision of additional negative charge to the di-Mg coordination sphere have significant effects on the steps of alkyl migration and hydride transfer, respectively. Other important mechanistic characteristics are also revealed. Inspired by the mechanism, an inhibitor (2-carboxylate-lactic acid) was designed and predicted by barrier analysis to be effective in inactivating KARI, hence probably enriching the antifungal and antibacterial library. Two types of slow substrate analogues (2-trihalomethyl acetolactic acids and 2-glutaryl lactic acid) were also found.

Effector-Binding-Directed Dimerization and Dynamic Communication between Allosteric Sites of Ribonucleotide Reductase

Proteins forming dimers or larger complexes can be strongly influenced by their effector-binding status. We investigated how the effector-binding event is coupled with interface formation via computer simulations, and we quantified the correlation of two types of contact interactions: between the effector and its binding pocket and between protein monomers. This was achieved by connecting the protein dynamics at the monomeric level with the oligomer interface information. We applied this method to ribonucleotide reductase (RNR), an essential enzyme for de novo DNA synthesis. RNR contains two important allosteric sites, the s-site (specificity site) and the a-site (activity site), which bind different effectors. We studied these different binding states with atomistic simulation and used their coarse-grained contact information to analyze the protein dynamics. The results reveal that the effector-protein dynamics at the s-site and dimer interface formation are positively coupled. We further quantify the resonance level between these two events, which can be applied to other similar systems. At the a-site, different effector-binding states (ATP vs dATP) drastically alter the protein dynamics and affect the activity of the enzyme. On the basis of these results, we propose a new mechanism of how the a-site regulates enzyme activation.

How does Mo-dependent perchlorate reductase work in the decomposition of oxyanions?

The mechanisms of Mo-dependent perchlorate reductase (PcrAB)-catalyzed decomposition of perchlorate, bromate, iodate, and nitrate were revealed by density functional calculations.

How does binuclear zinc amidohydrolase FwdA work in the initial step of methanogenesis: From formate to formyl-methanofuran

The initial step of methanogenesis is the fixation of CO₂ to formyl-methanofuran (formyl-MFR) catalyzed by formyl-MFR dehydrogenase, which can be divided into two half reactions. Herein, the second half reaction catalyzed by FwdA (formyl-methanofuran dehydrogenase subunit A), i.e., from formate to formyl-methanofuran, has been investigated using density functional theory and a chemical model based on the X-ray crystal structure. The calculations indicate that, compared with other well-known di-zinc hydrolases, the FwdA reaction employs a reverse mechanism, including the nucleophilic attack of MFR amine on formate carbon leading to a tetrahedral gem-diolate intermediate, two steps of proton transfer from amine to formate moieties assisted by the Asp385, and the CO bond dissociation to form the formyl-MFR product. The second step of proton transfer from the amine moiety to the Asp385 is rate-limiting with an overall barrier of 21.2 kcal/mol. The two zinc ions play an important role in stabilizing the transition states and intermediates, in particular the negative charge at the formate moiety originated from the nucleophilic attack of the MFR amine. The work here appends a crucial piece in the methanogenic mechanistics and advances the understanding of the global carbon cycle.

The Decarboxylation of α,β-Unsaturated Acid Catalyzed by Prenylated FMN-Dependent Ferulic Acid Decarboxylase and the Enzyme Inhibition

Ferulic acid decarboxylase (Fdc1) is able to catalyze the decarboxylation of α,β-unsaturated acids using a novel cofactor, prenylated flavin mononucleotide (PrFMN). Using density functional theory calculations, we here have investigated the Fdc1 reaction mechanism with the substrate of α-methylcinnamic acid. It is demonstrated that Fdc1 employs a 1,3-dipolar cycloaddition mechanism involving four concerted steps, where the Glu282 acts as a crucial proton donor to protonate the α carbon (C_α). The last step, the decomposition of a pyrrolidine species, is rate-limiting with an overall barrier of 18.9 kcal mol^-1. Furthermore, when α-hydroxycinnamic acid is used, the Glu282 is found to have another face to transport the hydroxyl proton to the C_β atom to promote the tautomerization from enol intermediate to ketone species leading to the inhibition of the Fdc1 enzyme. The PrFMN roles are also discussed in detail.

Non-enzymatic ribonucleotide reduction in the prebiotic context

Model studies of prebiotic chemistry have revealed compelling routes for the formation of the building blocks of proteins and RNA, but not DNA. Today, deoxynucleotides required for the construction of DNA are produced by reduction of nucleotides catalysed by ribonucleotide reductases, which are radical enzymes. This study considers potential non-enzymatic routes via intermediate radicals for the ancient formation of deoxynucleotides. In this context, several mechanisms for ribonucleotide reduction, in a putative H2 S/HS(.) environment, are characterized using computational chemistry. A bio-inspired mechanistic cycle involving a keto intermediate and HSSH production is found to be potentially viable. An alternative pathway, proceeding through an enol intermediate is found to exhibit similar energetic requirements. Non-cyclical pathways, in which HSS(.) is generated in the final step instead of HS(.) , show a markedly increased thermodynamic driving force (ca. 70 kJ mol(-1) ) and thus warrant serious consideration in the context of the prebiotic ribonucleotide reduction.

The catalytic mechanism of carboxylesterases: a computational study

The catalytic mechanism of carboxylesterases (CEs, EC 3.1.1.1) is explored by computational means. CEs hydrolyze ester, amide, and carbamate bonds found in xenobiotics and endobiotics. They can also perform transesterification, a reaction important, for instance, in cholesterol homeostasis. The catalytic mechanisms with three different substrates (ester, thioester, and amide) have been established at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory. It was found that the reactions proceed through a mechanism involving four steps instead of two as is generally proposed: (i) nucleophilic attack of serine to the substrate, forming the first tetrahedral intermediate, (ii) formation of the acyl-enzyme complex concomitant with the release of the alcohol product, (iii) nucleophilic attack of a water or alcohol molecule forming the second tetrahedral intermediate, and (iv) the release of the second product of the reaction. The results agree very well with the available experimental data and show that the hydrolytic and the transesterification reactions are competitive processes when the substrate is an ester. In all the other studied substrates (thioester or amide), the hydrolytic and transesterification process are less favorable and some of them might not even take place under in vivo conditions.

Isomerization of Δ5-androstene-3,17-dione into Δ4-androstene-3,17-dione catalyzed by human glutathione transferase A3-3: a computational study identifies a dual role for glutathione

Glutathione transferases (GSTs) are important enzymes in the metabolism of electrophilic xenobiotic and endobiotic toxic compounds. In addition, human GST A3-3 also catalyzes the double bond isomerization of Δ5-androstene-3,17-dione (Δ(5)-AD) and Δ(5)-pregnene-3,20-dione (Δ(5)-PD), which are the immediate precursors of testosterone and progesterone. In fact, GST A3-3 is the most efficient human enzyme known to exist in the catalysis of these reactions. In this work, we have used density functional theory (DFT) calculations to propose a refined mechanism for the isomerization of Δ(5)-AD catalyzed by GST A3-3. In this mechanism the glutathione (GSH) thiol and Tyr9 catalyze the proton transfer from the Δ(5)-AD C4 atom to the Δ(5)-AD C6 atom, with a rate limiting activation energy of 15.8 kcal · mol(-1). GSH has a dual function, because it is also responsible for stabilizing the negative charge that is formed in the O3 atom of the enolate intermediate. The catalytic role of Tyr9 depends on significant conformational rearrangements of its side chain. Neither of these contributions to catalysis has been observed before. Residues Phe10, Leu111, Ala 208, and Ala 216 complete the list of the important catalytic residues. The mechanism detailed here is based on the GST A3-3:GSH:Δ(4)-AD crystal structure and is consistent with all available experimental data.

Binding mode prediction and identification of new lead compounds from natural products as renin and angiotensin converting enzyme inhibitors

In this study a novel renin and ACE inhibitor was developed from natural products using computational techniques. Molecular dynamic simulations showed that the new lead compound has significant binding to the targets.

Mechanism of glutathione transferase P1-1-catalyzed activation of the prodrug canfosfamide (TLK286, TELCYTA)

Canfosfamide (TLK286, TELCYTA) is a prodrug that upon activation by glutathione transferase P1-1 (GST P1-1) yields an anticancer alkylating agent and a glutathione derivative. The rationale underlying the use of TLK286 in chemotherapy is that tumor cells overexpressing GST P1-1 will be locally exposed to the released alkylating agent with limited collateral toxicity to the surrounding normal tissues. TLK286 has demonstrated clinical effects in phase II and III clinical trials for the treatment of malignancies, such as ovarian cancer, nonsmall cell lung cancer, and breast cancer, as a single agent and in combination with other chemotherapeutic agents. In spite of these promising results, the detailed mechanism of GST P1-1 activation of the prodrug has not been elucidated. Here, we propose a mechanism for the TLK286 activation by GST P1-1 on the basis of density functional theory (DFT) and on potential of mean force (PMF) calculations. A catalytic water molecule is instrumental to the activation by forming a network of intermolecular interactions between the active-site Tyr7 hydroxyl and the sulfone and COO(-) groups of TLK286. The results obtained are consistent with the available experimental kinetic data and provide an atomistic understanding of the TLK286 activation mechanism.

The sulfur shift: an activation mechanism for periplasmic nitrate reductase and formate dehydrogenase

A structural rearrangement known as sulfur shift occurs in some Mo-containing enzymes of the DMSO reductase family. This mechanism is characterized by the displacement of a coordinating cysteine thiol (or SeCys in Fdh) from the first to the second shell of the Mo-coordination sphere metal. The hexa-coordinated Mo ion found in the as-isolated state cannot bind directly any exogenous ligand (substrate or inhibitors), while the penta-coordinated ion, attained upon sulfur shift, has a free binding site for direct coordination of the substrate. This rearrangement provides an efficient mechanism to keep a constant coordination number throughout an entire catalytic pathway. This mechanism is very similar to the carboxylate shift observed in Zn-dependent enzymes, and it has been recently detected by experimental means. In the present paper, we calculated the geometries and energies involved in the sulfur-shift mechanism using QM-methods (M06/(6-311++G(3df,2pd),SDD)//B3LYP/(6-31G(d),SDD)). The results indicated that the sulfur-shift mechanism provides an efficient way to enable the metal ion for substrate coordination.

Computational enzymatic catalysis

Computational methodologies are playing increasingly important roles in elucidating and presenting the complete and detailed mechanisms of enzymatic reactions because of their capacity to determine and characterize intermediates and transition states from both structural and energetics points of view, independent of their reduced lifetimes and without interfering with the natural reactional flux. These features are turning the field into an active and interesting area of research, involving a diverse range of studies, mostly directed at understanding the ways in which enzymes function under certain circumstances and predicting how they will behave under others. The accuracy of the computational data obtained for a given mechanistic hypothesis depends essentially on three mutually exclusive factors: the accuracy of the Hamiltonian of the reaction mechanism, consideration of the modulating aspect of the enzyme's structure in the energetics of the active center, and consideration of the enzyme's conformational fluctuations and dynamics. Although, unfortunately, it is impossible at present to optimize these crucial factors simultaneously, the success of any enzymatic mechanistic study depends on the level of equilibrium achieved among them. Different authors adopt different solutions, and this Account summarizes the most favored, with emphasis placed on our own preferences. Another crucial aspect in computational enzymatic catalysis is the model used in the calculations. Our aim is to build the simplest model that captures the essence of the catalytic power of an enzyme, allowing us to apply the highest possible theoretical level and minimize accidental errors. The choice is, however, far from obvious, ranging from simple models containing tens of atoms up to models of full enzymes plus solvent. Many factors underlie the choice of an appropriate model; here, examples are presented of very different modeling strategies that have been employed to obtain meaningful results. One particular case study, that of enzyme ribonucleotide reductase (RNR), a radical enzyme that catalyzes the reduction of ribonucleotides into deoxyribonucleotides, is one of the examples illustrating how the successive increase of the system's size does not dramatically change the thermodynamics and kinetics of the reaction. The values obtained and presented speak for themselves in that the only ones that are distinctly different are those calculated using an exceedingly small model, which omitted the amino acids that establish hydrogen bonds with the reactive unit of the substrate. This Account also describes our computational analysis of the mechanism of farnesyltransferase, a heterodimeric zinc metalloenzyme that is currently one of the most fascinating targets in cancer research. We focus on the present methodologies that we have been using, our models and understanding of the problem, and the accuracy of results and associated problems within this area of research.

Understanding Ribonucleotide Reductase Inactivation by Gemcitabine

This paper focuses on the inhibition of ribonucleotide reductase (RNR) by gemcitabine, 2',2'-difluoro-2'-deoxycytidine (dFdC), a deoxycytidine analogue that is a very active drug against solid tumors and is currently commercialized as gemzar. RNR inactivation is reductant-dependent and occurs in a very different way from that of other known substrate analogues. In the presence of reductants monomer R1 of RNR is inhibited, whereas in the absence of reductants the radical is lost and monomer R2 is inhibited. As inside the cell reductants are available, it is likely that R1 inactivation is the most favorable mechanism responsible for drug cytotoxicity. This inhibition pathway has been unknown to date, but we have conducted a theoretical study that has led us to the first proposal of a mechanism for RNR inhibition by dFdC in the presence of reductants.

Enzyme ribonucleotide reductase: unraveling an enigmatic paradigm of enzyme inhibition by furanone derivatives

Several 2'-substituted-2'-deoxyribonucleotides are potent inactivators of the enzyme ribonucleotide reductase (RNR), by destroying the essential tyrosyl radical located in subunit R2 or/and covalently alkylating the subunit R1. In the absence of external reductants, the inactivation is achieved by alkylation of subunit R1 by a methylene-3(2H)-furanone. The furanone is generated in solution through degradation of a keto-deoxyribonucleotide intermediate, produced during the inhibitory mechanism of a wide group of 2'-substituted inhibitors, and is easily detected experimentally by UV spectroscopy. Interestingly, the same keto-deoxyribonucleotide is also a proposed intermediate of the normal substrate pathway, but by some unknown reason, it does not dissociate from the active site and does not inactivate the enzyme. Therefore, if the currently accepted mechanism for substrate reduction is correct, there must be some specific reason that makes such a reactive intermediate behave differently, not dissociating from the active site during substrate reduction. In this article, we propose to validate the current substrate mechanism by showing that the keto-deoxyribonucleotide dissociates from the active site only in the case of the inhibitors, and therefore, it corresponds to a viable intermediate in the substrate mechanism. Furthermore, we answer unexplained experimental observations that concern the predomination of the normal reduction mechanism over the abnormal ketone formation in the FdNDP and the release of F(-), either in the normal or in the abnormal turnover. For that purpose, we have investigated the interaction between the enzyme and this keto-deoxyribonucleotide generated from the normal substrate and from two widely studied representative inhibitors. A model containing 140 atoms was used to represent the desired structures. The results allowed us to conclude that the solvation free energy of the 2'-substituents, its influence inside the active site, and the charge transfer mechanism from a protein side chain to solution are the thermodynamic driving forces for the intermediate dissociation and subsequent RNR inhibition. Such charge transfer cannot be accomplished by the natural substrate, preventing its dissociation. These results elucidate a paradox which has been unexplained for more than 20 years and further validates both the proposed substrate and inhibition chemical mechanisms.

The Carboxylate Shift in Zinc Enzymes: A Computational Study

Zinc is the second most abundant transition element in biology and the only metal known to be represented in enzymes from each one of the six classes established by the International Union of Biochemistry. The flexible coordination geometry, the fast ligand exchange, the lack of redox activity, and its role as Lewis acid are just some of the features that make zinc an invaluable element in biological catalysis. In this study, we have analyzed the importance in mononuclear Zn enzymes of an interesting mechanistic phenomenon known as carboxylate shift, which is characterized by a change in the coordination mode of a carboxylate group (mono to bidentate or vice versa) with both ligand entrance or exit from the metal coordination sphere. Using B3LYP calculations, we were able to unveil in detail patterns relating the intrinsic characteristics of a given Zn coordination sphere with the existence or not of a carboxylate-shift mechanism and the additional energy stabilization arising from it. In particular, a specific Zn coordination sphere containing a carboxylate ligand (Asp or Glu), a cysteine, and a histidine has been shown to have the most favorable combination of amino acid residues that ensures a fast ligand exchange.

A Combined ONIOM Quantum Chemical−Molecular Dynamics Study of Zinc−Uracil Bond Breaking in Yeast Cytosine Deaminase

A QM/MM method that combines ONIOM quantum chemistry and molecular dynamics is developed and applied to a step in the deamination of cytosine to uracil in yeast cytosine deaminase (yCD). A two-layer ONIOM calculation is used for the reaction complex, with an inner part treated at a high level for the chemical reaction (bond breaking) and a middle part treated at a lower level for relevant protein residues that are frozen in the quantum optimization. An outer layer (protein and solvent) is treated using MD. Configurations for the entire system are generated using MD and optimized with ONIOM. The method permits the use of high-level quantum calculations along with sufficient configurational sampling to approximate the potential of mean force for certain bond-breaking reactions. A previously proposed reaction mechanism for deamination (Sklenak, S.; Yao, L. S.; Cukier, R. I.; Yan, H. G. J. Am. Chem. Soc. 2004, 126, 14879) requires breaking the bond between a catalytic zinc and the O4 of uracil in order to permit product release. Using an ONIOM approach, direct bond cleavage was found to be energetically unfavorable. In the work presented here, the combined ONIOM MD method is used to show that the barrier for bond cleavage is small, approximately 3 kcal/mol, and, consequently, should not be the rate-limiting step in the reaction.

Protein–cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: Effect of hydrogen bonding on the electronic and geometric structure of the primary quinone. A density functional theory study

The effect of hydrogen bonding to the primary quinone (Q(A) and Q(*)(-)(A)) in bacterial reaction centers was studied using density functional theory (DFT) calculations. The charge neutral state Q(A) was investigated by optimizing the hydrogen atom positions of model systems extracted from 15 different X-ray structures. From this analysis, mean values of the H-bond lengths and directions were derived. It was found that the N(delta)-H of His M219 forms a shorter H-bond to Q(A) than the N-H of Ala M260. The H-bond of His M219 is linear and more twisted out of the quinone plane. The radical anion Q(*)(-)(A) in the protein environment was investigated by using a mixed quantum mechanics/molecular mechanics (QM/MM) approach. Two geometry optimizations with a different number of flexible atoms were performed. H-bond lengths were obtained and spectroscopic parameters calculated, i.e. the hyperfine and nuclear quadrupole couplings of magnetic nuclei coupled to the radical. Good agreement was found with the results provided by EPR/ENDOR spectroscopy. This implies that the calculated lengths and directions of the H-bonds to Q(*)(-)(A) are reliable values. From a comparison of the neutral and reduced state of Q(A) it was concluded that the H-bond distances are shortened by approximately 0.17 Angstroms (His M219) and approximately 0.13 Angstroms (Ala M260) upon single reduction of the quinone. It is shown that the point-dipole approximation can not be used for an estimation of H-bond lengths from measured hyperfine couplings in a system with out-of-plane H-bonding. In contrast, the evaluation of the nuclear quadrupole couplings of (2)H nuclei substituted in the hydrogen bonds yields H-bond lengths close to the values that were deduced from DFT geometry optimizations. The significance of hydrogen bonding to the quinone cofactors in biological systems is discussed.