Theoretical study of ribonucleotide reductase mechanism-based inhibition by 2'-azido-2'-deoxyribonucleoside 5'-diphosphates

2'-azido-2'-deoxyribonucleoside 5'-diphosphates are mechanism-based inhibitors of Ribonucleotide Reductase. Considerable effort has been made to elucidate their mechanism of inhibition, which is still controversial and not fully understood. Previous studies have detected the formation of a radical intermediate when the inhibitors interact with the enzyme, and several authors have proposed possible structures for this radical. We have conducted a theoretical study of the possible reactions involved, which allowed us to identify the structure of the new radical among the several proposals. A new reactional path is also proposed that is the most kinetically favored to yield this radical and ultimately inactivate the enzyme. The energetic involved in this mechanism, both for radical formation and radical decay, as well as the calculated Hyperfine Coupling Constants for the radical intermediate, are in agreement with the correspondent experimental values. This mechanistic alternative is fully coherent with remaining experimental data.     

EPR distance measurements support a model for long-range radical initiation in E. coli ribonucleotide reductase

The class I E. coli ribonucleotide reductase, composed of homodimers of R1 and R2, catalyzes the conversion of nucleoside diphosphates to deoxynucleoside diphosphates. The reduction process involves the tyrosyl radical on R2 that generates a transient thiyl radical on R1 over a proposed distance of 35 A. A mechanism-based inhibitor, 2'-azido-2'-deoxyuridine-5'-diphosphate, that reduces the tyrosyl radical on R2 and forms a nitrogen-centered radical on R1 has provided a method to measure the diagonal distance between the two subunits. PELDOR and DQC paramagnetic resonance methods give rise to a distance of 48 A, similar to that calculated from a docking model of the R1 and R2 structures.     

PELDOR spectroscopy with DOPA-beta2 and NH2Y-alpha2s: distance measurements between residues involved in the radical propagation pathway of E. coli ribonucleotide reductase

Escherichia coli ribonucleotide reductase (RNR) catalyzes the reduction of nucleotides to 2'-deoxynucleotides. The active enzyme is a 1:1 complex of two homodimeric subunits, alpha2 and beta2. The alpha2 is the site of nucleotide reduction, and beta2 harbors a diferric tyrosyl radical (Y122*) cofactor. Turnover requires formation of a cysteinyl radical (C439*) in the active site of alpha2 at the expense of the Y122* in beta2. A docking model for the alpha2beta2 interaction and a pathway for radical transfer from beta2 to alpha2 have been proposed. This pathway contains three Ys: Y356 in beta2 and Y731/Y730 in alpha2. We have previously incorporated 3-hydroxytyrosine and 3-aminotyrosine into these residues and showed that they act as radical traps. In this study, we use these alpha2/beta2 variants and PELDOR spectroscopy to measure the distance between the Y122* in one alphabeta pair and the newly formed radical in the second alphabeta pair. The results yield distances that are similar to those predicted by the docking model for radical transfer. Further, they support a long-range radical initiation process for C439* generation and provide a structural constraint for residue Y356, which is thermally labile in all beta2 structures solved to date.     

Site-specific insertion of 3-aminotyrosine into subunit alpha2 of E. coli ribonucleotide reductase: direct evidence for involvement of Y730 and Y731 in radical propagation

E. coli ribonucleotide reductase (RNR) catalyzes the production of deoxynucleotides using complex radical chemistry. Active RNR is composed of a 1:1 complex of two subunits: alpha2 and beta2. Alpha2 binds nucleoside diphosphate substrates and deoxynucleotide/ATP allosteric effectors and is the site of nucleotide reduction. Beta2 contains the stable diiron tyrosyl radical (Y122.) cofactor that initiates deoxynucleotide formation. This process is proposed to involve reversible radical transfer over >35 A between the Y122 in beta2 and C439 in the active site of alpha2. A docking model of alpha2beta2, based on structures of the individual subunits, suggests that radical initiation involves a pathway of transient, aromatic amino acid radical intermediates, including Y730 and Y731 in alpha2. In this study the function of residues Y730 and Y731 is investigated by their site-specific replacement with 3-aminotyrosine (NH2Y). Using the in vivo suppressor tRNA/aminoacyl-tRNA synthetase method, Y730NH2Y-alpha2 and Y731NH2Y-alpha2 have been generated with high fidelity in yields of 4-6 mg/g of cell paste. These mutants have been examined by stopped flow UV-vis and EPR spectroscopies in the presence of beta2, CDP, and ATP. The results reveal formation of an NH2Y radical (NH2Y730. or NH2Y731.) in a kinetically competent fashion. Activity assays demonstrate that both NH2Y-alpha2s make deoxynucleotides. These results show that the NH2Y. can oxidize C439 suggesting a hydrogen atom transfer mechanism for the radical propagation pathway within alpha2. The observed NH2Y. may constitute the first detection of an amino acid radical intermediate in the proposed radical propagation pathway during turnover.     

Investigation of reactions postulated to occur during inhibition of ribonucleotide reductases by 2'-azido-2'-deoxynucleotides

Model 3'-azido-3'-deoxynucleosides with thiol or vicinal dithiol substituents at C2' or C5' were synthesized to study reactions postulated to occur during inhibition of ribonucleotide reductases by 2'-azido-2'-deoxynucleotides. Esterification of 5'-(tert-butyldiphenylsilyl)-3'-azido-3'-deoxyadenosine and 3'-azido-3'-deoxythymidine (AZT) with 2,3-S-isopropylidene-2,3-dimercaptopropanoic acid or N-Boc-S-trityl-L-cysteine and deprotection gave 3'-azido-3'-deoxy-2'-O-(2,3-dimercaptopropanoyl or cysteinyl)adenosine and the 3'-azido-3'-deoxy-5'-O-(2,3-dimercaptopropanoyl or cysteinyl)thymidine analogs. Density functional calculations predicted that intramolecular reactions between generated thiyl radicals and an azido group on such model compounds would be exothermic by 33.6-41.2 kcal/mol and have low energy barriers of 10.4-13.5 kcal/mol. Reduction of the azido group occurred to give 3'-amino-3'-deoxythymidine, which was postulated to occur with thiyl radicals generated by treatment of 3'-azido-3'-deoxy-5'-O-(2,3-dimercaptopropanoyl)thymidine with 2,2'-azobis-(2-methyl-2-propionamidine) dihydrochloride. Gamma radiolysis of N(2)O-saturated aqueous solutions of AZT and cysteine produced 3'-amino-3'-deoxythymidine and thymine most likely by both radical and ionic processes.     

Site-specific replacement of Y356 with 3,4-dihydroxyphenylalanine in the beta2 subunit of E. coli ribonucleotide reductase

E. coli ribonucleotide reductase (RNR), composed of the homodimeric subunits alpha2 and beta2, catalyzes the conversion of nucleotides to deoxynucleotides via complex radical chemistry. The radical initiation process involves a putative proton-coupled electron transfer (PCET) pathway over 35 A between alpha2 and beta2. Y356 in beta2 has been proposed to lie on this pathway. To test this model, intein technology has been used to make beta2 semi-synthetically in which Y356 is replaced with a DOPA-amino acid. Analysis of this mutant with alpha2 and various combinations of substrate and effector by SF UV-vis spectroscopy and EPR methods demonstrates formation of a DOPA radical concomitant with disappearance of the tyrosyl radical, which initiates the reaction. The results reveal that Y356 lies on the PCET pathway and demonstrate the first kinetically competent conformational changes prior to ET. They further show that substrate binding brings about rapid conformational changes which place the complex into its active form(s) and suggest that the RNR complex is asymmetric.     

Stereodefined synthesis of O3'-labeled uracil nucleosides. 3'-[(17)O]-2'-Azido-2'-deoxyuridine 5'-diphosphate as a probe for the mechanism of inactivation of ribonucleotide reductases

Thermolysis of a 2'-[(16)O]-O-benzoyl-[(17)O]-5'-O-(tert-butyldimethylsilyl)-O(2),3'-cyclouridine derivative gave the more stable 3'-[(17)O]-O-benzoyl-[(16)O]- 5'-O-(tert-butyldimethylsilyl)-O(2),2'-cyclouridine isomer, which was converted into 3'-[(17)O]-2'-azido-2'-deoxyuridine by deprotection and nucleophilic ring opening at C2' with lithium azide. The 5'-diphosphate was prepared by nucleophilic displacement of the 5'-O-tosyl group with tris(tetrabutylammonium) hydrogen pyrophosphate. Model reactions gave (16)O and (18)O isotopomers, and base-promoted hydrolysis of an O(2),2'-cyclonucleoside gave stereodefined access to 3'-[(18)O]-1-(beta-D-arabinofuranosyl)uracil. Inactivation of ribonucleoside diphosphate reductase with 2'-azido-2'-deoxynucleotides results in appearance of EPR signals for a nitrogen-centered radical derived from azide, and 3'-[(17)O]-2'-azido-2'-deoxyuridine 5'-diphosphate provides an isotopomer to perturb EPR spectra in a predictable manner.     

Theoretical study on the inhibition of ribonucleotide reductase by 2'-mercapto-2'-deoxyribonucleoside-5'-diphosphates

Ribonucleotide reductase (RNR) is responsible for the reduction of ribonucleotides into the correspondent 2'-deoxyribonucleotides in the only physiological process that yields the monomers of DNA. The enzyme has thus become an attractive target for chemotherapies that fight proliferation-based diseases, specifically cancer and infections by some viruses and parasites. 2'-Mercapto-2'-deoxyribonucleoside-5'-diphosphates (SHdNDP) are mechanism-based inhibitors of RNR and therefore potential chemotherapeutic agents for those indications. Previous experimental studies established the in vitro and in vivo activity of SHdNDP. In the in vitro studies, it was observed that the activity was dependent on the oxidative status of the medium, with the inactivation of RNR only occurring when molecular oxygen was available. To better understand the mechanism involved in RNR inactivation by SHdNDP, we performed theoretical calculations on the possible reactions between the inhibitors and the RNR active site. As a result, we propose the possible mechanistic pathways for the chemical events that occur in the absence and in the presence of O2. They correspond to a refinement and a complement of those proposed in the literature.     

Pulsed ELDOR spectroscopy measures the distance between the two tyrosyl dadicals in the R2 subunit of the E. coli ribonucleotide reductase

Escherichia coli ribonucleotide reductase (RNR) catalyzes the conversion of nucleoside diphosphates (NDPs) to deoxynucleoside diphosphates (dNDPs). This RNR is composed of two homodimeric subunits: R1 and R2. R1 binds the NDPs in the active site, and R2 harbors the essential di-iron tyrosyl radical (Y*) cofactor. In this paper, we used PELDOR, a method that detects weak electron-electron dipolar coupling, to make the first direct measurement of the distance between the two Y*'s on each monomer of R2. In the crystal structure of R2, the Y*'s are reduced to tyrosines, and consequently R2 is inactive. In R2, where the Y*'s assume a well-defined geometry with respect to the protein backbone, the PELDOR method allows measurement of a distance of 33.1 +/- 0.2 A that compares favorably to the distance (32.4 A) between the center of mass of the spin density distribution of each Y* on each R2 monomer from the structure. The experiments provide the first direct experimental evidence for two Y*'s in a single R2 in solution.     

Rational reprogramming of the R2 subunit of Escherichia coli ribonucleotide reductase into a self-hydroxylating monooxygenase.

The outcome of O2 activation at the diiron(II) cluster in the R2 subunit of Escherichia coli (class I) ribonucleotide reductase has been rationally altered from the normal tyrosyl radical (Y122*) production to self-hydroxylation of a phenylalanine side-chain by two amino acid substitutions that leave intact the (histidine)2-(carboxylate)4 ligand set characteristic of the diiron-carboxylate family. Iron ligand Asp (D) 84 was replaced with Glu (E), the amino acid found in the cognate position of the structurally similar diiron-carboxylate protein, methane monooxygenase hydroxylase (MMOH). We previously showed that this substitution allows accumulation of a mu-1,2-peroxodiiron(III) intermediate, which does not accumulate in the wild-type (wt) protein and is probably a structural homologue of intermediate P (H(peroxo)) in O2 activation by MMOH. In addition, the near-surface residue Trp (W) 48 was replaced with Phe (F), blocking transfer of the "extra" electron that occurs in wt R2 during formation of the formally Fe(III)Fe(IV) cluster X. Decay of the mu-1,2-peroxodiiron(III) complex in R2-W48F/D84E gives an initial brown product, which contains very little Y122* and which converts very slowly (t1/2 approximately 7 h) upon incubation at 0 degrees C to an intensely purple final product. X-ray crystallographic analysis of the purple product indicates that F208 has undergone epsilon-hydroxylation and the resulting phenol has shifted significantly to become a ligand to Fe2 of the diiron cluster. Resonance Raman (RR) spectra of the purple product generated with 16O2 or 18O2 show appropriate isotopic sensitivity in bands assigned to O-phenyl and Fe-O-phenyl vibrational modes, confirming that the oxygen of the Fe(III)-phenolate species is derived from O2. Chemical analysis, experiments involving interception of the hydroxylating intermediate with exogenous reductant, and M?ssbauer and EXAFS characterization of the brown and purple species establish that F208 hydroxylation occurs during decay of the peroxo complex and formation of the initial brown product. The slow transition to the purple Fe(III)-phenolate species is ascribed to a ligand rearrangement in which mu-O2- is lost and the F208-derived phenolate coordinates. The reprogramming to F208 monooxygenase requires both amino acid substitutions, as very little epsilon-hydroxyphenylalanine is formed and pathways leading to Y122* formation predominate in both R2-D84E and R2-W48F.     

Theoretical studies on the mechanism of inhibition of Ribonucleotide Reductase by (E)-2'-Fluoromethylene-2'-deoxycitidine-5'-diphosphate

(E)-2'-Fluoromethylene-2'-deoxycitidine-5'-diphosphate (FMCDP) is a potent time-dependent inactivator of the enzyme Ribonucleotide Reductase, which operates by destructing an essential tyrosil radical and performing a covalent addition to an active site residue. Considerable effort to elucidate the inhibition mechanism has been undertaken in recent years, and some insights have been obtained. Although a mechanistic proposal has been put forward, based on a general paradigm of inhibition of RNR by 2' substituted substrate analogues, details about the mechanism have remained elusive. Namely, the exact residue that adds to the inhibitor is still not identified, although mutagenesis experiments suggest that it should correspond to the E441 residue. In this work, we performed an extensive theoretical exploration of the potential energy surface of a model system representing the active site of RNR with FMCDP, using Density Functional Theory. This study establishes the detailed mechanism of inhibition, which is considerably different from the one proposed earlier. The proposed mechanism is fully consistent with available experimental data. Energetic results reveal unambiguously that the residue adding to the inhibitor is a cysteine, most probably C439, and exclude the possibility of the addition of E441. However, the E441 residue is shown to be essential for inhibition, catalyzing both the major and minor inhibition pathways, in agreement with mutagenesis results. It is shown also that the major mode of inactivation mimics the early stages of the natural substrate pathway.     

DFT calculations for intermediate and active states of the diiron center with a tryptophan or tyrosine radical in Escherichia coli ribonucleotide reductase

Class Ia ribonucleotide reductase subunit R2 contains a diiron active site. In this paper, active-site models for the intermediate X-Trp48(?+) and X-Tyr122(?), the active Fe(III)Fe(III)-Tyr122(?), and the met Fe(III)Fe(III) states of Escherichia coli R2 are studied, using broken-symmetry density functional theory incorporated with the conductor-like screening solvation model. Different structural isomers and different protonation states have been explored. Calculated geometric, energetic, Mo?ssbauer, hyperfine, and redox properties are compared with available experimental data. Feasible detailed structures of these intermediate and active states are proposed. Asp84 and Trp48 are most likely the main contributing residues to the result that the transient Fe(IV)Fe(IV) state is not observed in wild-type class Ia E. coli R2. Asp84 is proposed to serve as a proton-transfer conduit between the diiron cluster and Tyr122 in both the tyrosine radical activation pathway and the first steps of the catalytic proton-coupled electron-transfer pathway. Proton-coupled and simple redox potential calculations show that the kinetic control of proton transfer to Tyr122(?) plays a critical role in preventing reduction from the active Fe(III)Fe(III)-Tyr122(?) state to the met state, which is potentially the reason why Tyr122(?) in the active state can be stable over a very long period.     

Ribonucleotide activation by enzyme ribonucleotide reductase: understanding the role of the enzyme

This article focuses on the first step of the catalytic mechanism for the reduction of ribonucleotides catalyzed by the enzyme Ribonucleotide Reductase (RNR). This corresponds to the activation of the substrate. In this work a large model of the active site region involving 130 atoms was used instead of the minimal gas phase models used in previous works. The ONIOM method was employed to deal with such a large system. The results gave additional information, which previous small models could not provide, allowing a much clearer evaluation of the role of the enzyme in this step. Enzyme-substrate interaction energies, specific transition state stabilization, and substrate steric strain energies were obtained. It was concluded that the transition state is stabilized in 4.0 kcal/mol by specific enzyme-substrate interactions. However, this stabilization is cancelled by the cost in conformational energy for the enzyme to adopt the transition state geometry; the overall result is that the enzyme machinery does not lead to a rate enhancement in this step. It was also found that the substrate binds to the active site with almost no steric strain, emphasizing the complementarity and specificity of the RNR active site for nucleotide binding. The main role of the enzyme at the very beginning of the catalytic cycle was concluded to be to impose stereospecifity upon substrate activation and to protect the enzyme radical from the solvent, rather than to be an reaction rate enhancement.     

Understanding the molecular interactions of different radical scavengers with ribonucleotide reductase M2 (hRRM2) domain: opening the gates and gaining access

We employed a combination of molecular docking and dynamics to understand the interaction of three different radical scavengers (SB-HSC21, ABNM13 and trimidox) with ribonucleotide reductase M2 (hRRM2) domain. On the basis of the observed results, we can propose how these ligands interact with the enzyme, and cease the radical transfer step from the di-iron center to TYR176. All the ligands alter the electron density over TYR176, -OH group by forming an extremely stable H-bond with either -NHOH group, or with phenolic hydroxyl group of the ligands. This change in electronic density disrupts the water bridge between TYR176, -OH and the di-iron center, which stops the single electron transfer process from TYR176, -OH to iron. As a consequence the enzyme is inhibited. Another interesting observation that we are reporting is the two stage gate keeping mechanism of the RR active site tunnel. We describe these as the outer Gate-1 controlled by ARG330, and the inner Gate-2 controlled by SER263, PHE240, and PHE236. We also observed a dynamic conformational shift in these residues, the incoming ligands can go through, and interact with the underlying TYR176, -OH group. From the study we found the active-site of hRRM2 is extremely flexible and shows a significant induced fit.