Radical initiation in the class I ribonucleotide reductase: long-range proton-coupled electron transfer?

Hydrogen bonding networks tune proton-coupled redox steps during the enzymatic six-electron conversion of nitrite to ammonia

Multielectron multiproton reactions play an important role in both biological systems and chemical reactions involved in energy storage and manipulation. A key strategy employed by nature in achieving such complex chemistry is the use of proton-coupled redox steps. Cytochrome c nitrite reductase (ccNiR) catalyzes the six-electron seven-proton reduction of nitrite to ammonia. While a catalytic mechanism for ccNiR has been proposed on the basis of studies combining computation and crystallography, there have been few studies directly addressing the nature of the proton-coupled events that are predicted to occur along the nitrite reduction pathway. Here we use protein film voltammetry to directly interrogate the proton-coupled steps that occur during nitrite reduction by ccNiR. We find that conversion of nitrite to ammonia by ccNiR adsorbed to graphite electrodes is defined by two distinct phases; one is proton-coupled, and the other is not. Mutation of key active site residues (H257, R103, and Y206) modulates these phases and specifically alters the properties of the detected proton-dependent step but does not inhibit the ability of ccNiR to conduct the full reduction of nitrite to ammonia. We conclude that the active site residues examined are responsible for tuning the protonation steps that occur during catalysis, likely through an extensive hydrogen bonding network, but are not necessarily required for the reaction to proceed. These results provide important insight into how enzymes can specifically tune proton- and electron transfer steps to achieve high turnover numbers in a physiological pH range.

Assembly of nonheme Mn/Fe active sites in heterodinuclear metalloproteins

The ferritin superfamily contains several protein groups that share a common fold and metal coordinating ligands. The different groups utilize different dinuclear cofactors to perform a diverse set of reactions. Several groups use an oxygen-activating di-iron cluster, while others use di-manganese or heterodinuclear Mn/Fe cofactors. Given the similar primary ligand preferences of Mn and Fe as well as the similarities between the binding sites, the basis for metal specificity in these systems remains enigmatic. Recent data for the heterodinuclear cluster show that the protein scaffold per se is capable of discriminating between Mn and Fe and can assemble the Mn/Fe center in the absence of any potential assembly machineries or metal chaperones. Here we review the current understanding of the assembly of the heterodinuclear cofactor in the two different protein groups in which it has been identified, ribonucleotide reductase R2c proteins and R2-like ligand-binding oxidases. Interestingly, although the two groups form the same metal cluster they appear to employ partly different mechanisms to assemble it. In addition, it seems that both the thermodynamics of metal binding and the kinetics of oxygen activation play a role in achieving metal specificity.

Photo-induced proton-coupled electron transfer reactions of acridine orange: comprehensive spectral and kinetics analysis

The triplet excited state of acridine orange ((3)*AO) undergoes a proton-coupled electron transfer (PCET) reaction with tri-tert-butylphenol ((ttb)PhOH) in acetonitrile. Each of the reaction components possesses a spectroscopic signature, providing a rare opportunity to monitor the individual proton transfer, electron transfer, and H(•)-transfer components in parallel via transient absorption spectroscopy. This enhanced optical tracking, along with excited-state thermochemical analysis, facilitates assignment of the mechanism of excited-state PCET reactivity. (3)*AO is quenched via concerted proton-electron transfer (CPET) from (ttb)PhOH to form acridine radical (AOH(•)) and (ttb)PhO(•) (kCPET = 3.7 × 10(8) M(-1) s(-1), KIE = 1.3). Subsequently, AOH(•) reduces the phenoxyl radical (kET = 5.5 × 10(9) M(-1) s(-1)), forming AOH(+) and (ttb)PhO(-), followed by proton transfer (kPT = 1.0 × 10(9) M(-1) s(-1)) to regenerate the starting reactants.

Molybdopterin biosynthesis: trapping of intermediates for the MoaA-catalyzed reaction using 2'-deoxyGTP and 2'-chloroGTP as substrate analogues

MoaA is a radical S-adenosylmethionine (AdoMet) enzyme that catalyzes a complex rearrangement of guanosine-5'-triphosphate (GTP) in the first step of molybdopterin biosynthesis. In this paper, we provide additional characterization of the MoaA reaction product, describe the use of 2′-chloroGTP to trap the GTP C3′ radical, generated by hydrogen atom transfer to the 5′-deoxyadenosyl radical, and the use of 2′-deoxyGTP to block a late step in the reaction sequence. These probes, coupled with the previously reported trapping of an intermediate in which C3′ of the ribose is linked to C8 of the purine, allow us to propose a plausible mechanism for the MoaA-catalyzed reaction.

Biological charge transfer via flickering resonance

Biological electron-transfer (ET) reactions are typically described in the framework of coherent two-state electron tunneling or multistep hopping. However, these ET reactions may involve multiple redox cofactors in van der Waals contact with each other and with vibronic broadenings on the same scale as the energy gaps among the species. In this regime, fluctuations of the molecular structures and of the medium can produce transient energy level matching among multiple electronic states. This transient degeneracy, or flickering electronic resonance among states, is found to support coherent (ballistic) charge transfer. Importantly, ET rates arising from a flickering resonance (FR) mechanism will decay exponentially with distance because the probability of energy matching multiple states is multiplicative. The distance dependence of FR transport thus mimics the exponential decay that is usually associated with electron tunneling, although FR transport involves real carrier population on the bridge and is not a tunneling phenomenon. Likely candidates for FR transport are macromolecules with ET groups in van der Waals contact: DNA, bacterial nanowires, multiheme proteins, strongly coupled porphyrin arrays, and proteins with closely packed redox-active residues. The theory developed here is used to analyze DNA charge-transfer kinetics, and we find that charge-transfer distances up to three to four bases may be accounted for with this mechanism. Thus, the observed rapid (exponential) distance dependence of DNA ET rates over distances of ≲ 15 Å does not necessarily prove a tunneling mechanism.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

Electrostatics, hydration, and proton transfer dynamics in the membrane domain of respiratory complex I

Complex I serves as the primary electron entry point into the mitochondrial and bacterial respiratory chains. It catalyzes the reduction of quinones by electron transfer from NADH, and couples this exergonic reaction to the translocation of protons against an electrochemical proton gradient. The membrane domain of the enzyme extends ∼180 Å from the site of quinone reduction to the most distant proton pathway. To elucidate possible mechanisms of the long-range proton-coupled electron transfer process, we perform large-scale atomistic molecular dynamics simulations of the membrane domain of complex I from Escherichia coli. We observe spontaneous hydration of a putative proton entry channel at the NuoN/K interface, which is sensitive to the protonation state of buried glutamic acid residues. In hybrid quantum mechanics/classical mechanics simulations, we find that the observed water wires support rapid proton transfer from the protein surface to the center of the membrane domain. To explore the functional relevance of the pseudosymmetric inverted-repeat structures of the antiporter-like subunits NuoL/M/N, we constructed a symmetry-related structure of a possible alternate-access state. In molecular dynamics simulations, we find the resulting structural changes to be metastable and reversible at the protein backbone level. However, the increased hydration induced by the conformational change persists, with water molecules establishing enhanced lateral connectivity and pathways for proton transfer between conserved ionizable residues along the center of the membrane domain. Overall, the observed water-gated transitions establish conduits for the unidirectional proton translocation processes, and provide a possible coupling mechanism for the energy transduction in complex I.

Solid-state magnetic switching triggered by proton-coupled electron-transfer assisted by long-distance proton-alkali cation transport

Acidity of water molecules coordinated to Co ions in CoFe Prussian blue analogues (PBA) has been used to reversibly activate the Co(III)Fe(II) ↔ Co(II)Fe(III) electron transfer. The study of the structure and the electronic structure shows that the process implies an original PCET reaction between a solid-state porous coordination polymer and hydroxide ions in solution. The PCET reaction spreads throughout the solid network thanks to a long-range H(+) and Rb(+) transport within the pore channels of PBA taking advantage of the hydrogen-bonding network of zeolitic water molecules acting as proton wires.

Conserved electron donor complex Dre2-Tah18 is required for ribonucleotide reductase metallocofactor assembly and DNA synthesis

Eukaryotic ribonucleotide reductases (RNRs) require a diferric-tyrosyl radical (Fe(III)2-Y•) cofactor to produce deoxynucleotides essential for DNA replication and repair. This metallocofactor is an important target of RNR-based therapeutics, although mechanisms of in vivo cofactor assembly, inactivation, and reactivation are poorly understood. Here, we demonstrate that the conserved Fe-S protein-diflavin reductase complex, Dre2-Tah18, plays a critical role in RNR cofactor biosynthesis. Depletion of Dre2 affects both RNR gene transcription and mRNA turnover through the activation of the DNA-damage checkpoint and the Aft1/Aft2-controlled iron regulon. Under conditions of comparable RNR protein levels, cells with diminishing Dre2 have significantly reduced ability to make deoxynucleotides. Furthermore, the kinetics and levels of in vivo reconstitution of the RNR cofactor are severely impaired in two conditional tah18 mutants. Together, these findings provide insight into RNR cofactor formation and reveal a shared mechanism underlying assembly of the Fe(III)2-Y• cofactor in RNR and the Fe-S clusters in cytosolic and nuclear proteins.

Structural and molecular basis of the peroxynitrite-mediated nitration and inactivation of Trypanosoma cruzi iron-superoxide dismutases (Fe-SODs) A and B: disparate susceptibilities due to the repair of Tyr35 radical by Cys83 in Fe-SODB through intramolecular electron transfer

Trypanosoma cruzi, the causative agent of Chagas disease, contains exclusively iron-dependent superoxide dismutases (Fe-SODs) located in different subcellular compartments. Peroxynitrite, a key cytotoxic and oxidizing effector biomolecule, reacted with T. cruzi mitochondrial (Fe-SODA) and cytosolic (Fe-SODB) SODs with second order rate constants of 4.6 ± 0.2 × 10(4) M(-1) s(-1) and 4.3 ± 0.4 × 10(4) M(-1) s(-1) at pH 7.4 and 37 °C, respectively. Both isoforms are dose-dependently nitrated and inactivated by peroxynitrite. Susceptibility of T. cruzi Fe-SODA toward peroxynitrite was similar to that reported previously for Escherichia coli Mn- and Fe-SODs and mammalian Mn-SOD, whereas Fe-SODB was exceptionally resistant to oxidant-mediated inactivation. We report mass spectrometry analysis indicating that peroxynitrite-mediated inactivation of T. cruzi Fe-SODs is due to the site-specific nitration of the critical and universally conserved Tyr(35). Searching for structural differences, the crystal structure of Fe-SODA was solved at 2.2 Å resolution. Structural analysis comparing both Fe-SOD isoforms reveals differences in key cysteines and tryptophan residues. Thiol alkylation of Fe-SODB cysteines made the enzyme more susceptible to peroxynitrite. In particular, Cys(83) mutation (C83S, absent in Fe-SODA) increased the Fe-SODB sensitivity toward peroxynitrite. Molecular dynamics, electron paramagnetic resonance, and immunospin trapping analysis revealed that Cys(83) present in Fe-SODB acts as an electron donor that repairs Tyr(35) radical via intramolecular electron transfer, preventing peroxynitrite-dependent nitration and consequent inactivation of Fe-SODB. Parasites exposed to exogenous or endogenous sources of peroxynitrite resulted in nitration and inactivation of Fe-SODA but not Fe-SODB, suggesting that these enzymes play distinctive biological roles during parasite infection of mammalian cells.

Clofarabine targets the large subunit (α) of human ribonucleotide reductase in live cells by assembly into persistent hexamers

Clofarabine (ClF) is a drug used in the treatment of leukemia. One of its primary targets is human ribonucleotide reductase (hRNR), a dual-subunit, (α(2))(m)(β(2))(n), regulatory enzyme indispensable in de novo dNTP synthesis. We report that, in live mammalian cells, ClF targets hRNR by converting its α-subunit into kinetically stable hexamers. We established mammalian expression platforms that enabled isolation of functional α and characterization of its altered oligomeric associations in response to ClF treatment. Size exclusion chromatography and electron microscopy documented persistence of in-cell-assembled-α(6). Our data validate hRNR as an important target of ClF, provide evidence that in vivo α's quaternary structure can be perturbed by a nonnatural ligand, and suggest small-molecule-promoted, persistent hexamerization as a strategy to modulate hRNR activity. These studies lay foundations for documentation of RNR oligomeric state within a cell.

Redox-dependent structural coupling between the α2 and β2 subunits in E. coli ribonucleotide reductase

Ribonucleotide reductase (RNR) catalyzes the production of deoxyribonucleotides in all cells. In E. coli class Ia RNR, a transient α2β2 complex forms when a ribonucleotide substrate, such as CDP, binds to the α2 subunit. A tyrosyl radical (Y122O•)-diferric cofactor in β2 initiates substrate reduction in α2 via a long-distance, proton-coupled electron transfer (PCET) process. Here, we use reaction-induced FT-IR spectroscopy to describe the α2β2 structural landscapes, which are associated with dATP and hydroxyurea (HU) inhibition. Spectra were acquired after mixing E. coli α2 and β2 with a substrate, CDP, and the allosteric effector, ATP. Isotopic chimeras, (13)Cα2β2 and α2(13)Cβ2, were used to define subunit-specific structural changes. Mixing of α2 and β2 under turnover conditions yielded amide I (C═O) and II (CN/NH) bands, derived from each subunit. The addition of the inhibitor, dATP, resulted in a decreased contribution from amide I bands, attributable to β strands and disordered structures. Significantly, HU-mediated reduction of Y122O• was associated with structural changes in α2, as well as β2. To define the spectral contributions of Y122O•/Y122OH in the quaternary complex, (2)H4 labeling of β2 tyrosines and HU editing were performed. The bands of Y122O•, Y122OH, and D84, a unidentate ligand to the diferric cluster, previously identified in isolated β2, were observed in the α2β2 complex. These spectra also provide evidence for a conformational rearrangement at an additional β2 tyrosine(s), Yx, in the α2β2/CDP/ATP complex. This study illustrates the utility of reaction-induced FT-IR spectroscopy in the study of complex enzymes.

Long-range electron tunneling

Electrons have so little mass that in less than a second they can tunnel through potential energy barriers that are several electron-volts high and several nanometers wide. Electron tunneling is a critical functional element in a broad spectrum of applications, ranging from semiconductor diodes to the photosynthetic and respiratory charge transport chains. Prior to the 1970s, chemists generally believed that reactants had to collide in order to effect a transformation. Experimental demonstrations that electrons can transfer between reactants separated by several nanometers led to a revision of the chemical reaction paradigm. Experimental investigations of electron exchange between redox partners separated by molecular bridges have elucidated many fundamental properties of these reactions, particularly the variation of rate constants with distance. Theoretical work has provided critical insights into the superexchange mechanism of electronic coupling between distant redox centers. Kinetics measurements have shown that electrons can tunnel about 2.5 nm through proteins on biologically relevant time scales. Longer-distance biological charge flow requires multiple electron tunneling steps through chains of redox cofactors. The range of phenomena that depends on long-range electron tunneling continues to expand, providing new challenges for both theory and experiment.

Photochemical reaction mechanism of UV-B-induced monomerization of UVR8 dimers as the first signaling event in UV-B-regulated gene expression in plants

The Arabidopsis thaliana UV RESISTANCE LOCUS8 (UVR8) protein has been identified to specifically mediate photomorphogenic UV-B responses by acting as a UV-B photoreceptor. The dimeric structure of the UVR8 protein dissociates into signaling-active monomers upon UV-B exposure, and the monomers rapidly interact with downstream signaling components to regulate gene expression. UVR8 monomers revert to dimers in the absence of UV-B radiation, thereby reversing transcription activation. UVR8 amino acid residues W233 and W285 have been identified to play critical roles in the UVR8 dimer for the response to UV-B irradiation. In the present work, the photoreaction mechanism for UVR8 monomerization is explored with quantum chemical cluster calculations and evaluated by molecular dynamics simulations using the wild-type UVR8 dimer and novel force field parameters developed for intermediate radicals formed in the photochemical process. Three different models are investigated, which show that the preferred mechanism for UVR8 monomerization involves electron transfer from residue W233 to W285 and onward to R338 initiated by UV-B irradiation, coupled to simultaneous proton transfer from W233 to D129 leading to the formation of protonated D129, a deprotonated W233 radical, and a neutral R338 radical. Due to the formation of the neutral R338 radical, salt bridges involving this residue are disrupted together with the concomitant interruption of several other key salt bridges R286-D96, R286-D107, R338-D44, R354-E43, and R354-E53. The resulting large decrease in protein-protein interaction energy arising from this sequence of events leads to the monomerization of the UVR8 dimer. The mechanism presented is in accord with all experimental data available to date.

Breaking bonds with electrons and protons. Models and examples

Besides its theoretical interest, the attention currently aroused by proton-coupled electron transfers (PCET reactions) has two main motives. One is a better understanding of biological processes in which PCET reactions are involved, Photosystem II as well as a myriad of other natural systems. The other is directed toward synthetic processes, many of which are related to global energy challenges. Until recently, the analyses of the mechanism and reactivity of PCET reactions have focused on outersphere transfers, those in which no bond between heavy atoms (all atoms with the exception of H) is concomitantly formed or broken. Conversely, reactions in which electron transfer triggers the breaking of a heavy-atom bond with no proton transfer have been extensively analyzed, both theoretically and experimentally. In both cases, strategies have been developed to distinguish between stepwise and concerted pathways. In each case, kinetic models have been devised, allowing the relation between activation and thermodynamic driving force to be established by means of parameters pertaining to the initial and final state. Although many natural and artificial processes include electron transfer, proton transfer, and heavy-atom bond breaking (/formation), no means were offered until recently to analyze the mechanism of such reactions, notably to establish the degree of concertedness of the three constitutive events. Likewise, no kinetic models were available to describe reactions where the three events are concerted. In this Account, we discuss the strategies to distinguish stepwise, partially concerted (when two of the three events are concerted), and totally concerted pathways in these reactions that include electron transfer, proton transfer, and heavy-atom bond breaking. These mechanism analysis methods are illustrated and validated by three examples. First we describe the electrochemical cleavage of an O-O bond in an aliphatic peroxide molecule with a pendant carboxylic acid group that can serve as proton donor for electron transfer and bond breaking. In the second example, we examine the breaking of one of the C-O bonds of CO2 within a multistep process where the reduction of CO2 into CO is catalyzed by an electrogenerated iron(0) porphyrin in the presence of various Brönsted acids. In this case, an intramolecular electron transfer triggers proton transfer and bond cleavage. In the first two examples, all three events are concerted. The third example also involves catalysis. It describes the cleavage of a cobalt-carbon bond in the reduction of chloroacetonitrile catalyzed by an electrogenerated cobalt(I) porphyrin. It illustrates the rather common case where the intermediate formed by the reaction of a transition metal complex with the substrate has to be cleaved to close the catalytic cycle. In the first two examples, all three events are concerted, whereas, in the last case, a partially concerted pathway takes place: proton transfer and bond-breaking (Co-C cleavage) are concerted after an initial electron transfer step. The all-concerted cases require a model that connects the kinetics to the thermodynamic driving force of the reaction. Starting from previous models of outersphere electron transfer, concerted proton-electron transfer, and concerted dissociative electron transfer, we describe a model for all-concerted proton-electron-bond breaking reactions. These pathways skip the high-energy intermediates that occur in stepwise pathways, but could introduce kinetic penalties. The all-concerted model allows one to assess these penalties and the way in which they can be fought by the supplement of driving force offered by concerted proton transfer.

EXAFS simulation refinement based on broken-symmetry DFT geometries for the Mn(IV)-Fe(III) center of class I RNR from Chlamydia trachomatis

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides into deoxyribonucleotides necessary for DNA biosynthesis. Unlike the conventional class Ia RNRs which use a diiron cofactor in their subunit R2, the active site of the RNR-R2 from Chlamydia trachomatis (Ct) contains a Mn/Fe cofactor. The detailed structure of the Mn/Fe core has yet to be established. In this paper we evaluate six different structural models of the Ct RNR active site in the Mn(iv)/Fe(iii) state by using Mössbauer parameter calculations and simulations of Mn/Fe extended X-ray absorption fine structure (EXAFS) spectroscopy, and we identify a structure similar to a previously proposed DFT-optimized model that shows quantitative agreement with both EXAFS and Mössbauer spectroscopic data.

How protein structure affects redox reactivity: example of Human centrin 2

Human centrin 2 is a protein very sensitive to oxidative stress. Protein reactivity is unraveled by gamma radiolysis and electrochemical techniques.

Trifluoroacetic acid catalyzed thiophenylmethylation and thioalkylmethylation of lactams and phenols via domino three-component reaction in water

A new method for thiophenylmethylation and thioalkylmethylation of lactams and phenols has been described.

Solvent effect on electron and proton transfer in the excited state of a hydrogen bonded phenol–imidazole complex

Density functional theory calculations have been carried out for the ground state (S₀) and the first excited state (S₁) of the H-bonded phenol and imidazole complex as a model system for the active site of photosystem II.

Femtosecond dynamics of short-range protein electron transfer in flavodoxin

Intraprotein electron transfer (ET) in flavoproteins is important for understanding the correlation of their redox, configuration, and reactivity at the active site. Here, we used oxidized flavodoxin as a model system and report our complete characterization of a photoinduced redox cycle from the initial charge separation in 135-340 fs to subsequent charge recombination in 0.95-1.6 ps and to the final cooling relaxation of the product(s) in 2.5-4.3 ps. With 11 mutations at the active site, we observed that these ultrafast ET dynamics, much faster than active-site relaxation, mainly depend on the reduction potentials of the electron donors with minor changes caused by mutations, reflecting a highly localized ET reaction between the stacked donor and acceptor at a van der Waals distance and leading to a gas-phase type of bimolecular ET reaction confined in the active-site nanospace. Significantly, these ultrafast ET reactions ensure our direct observation of vibrationally excited reaction product(s), suggesting that the back ET barrier is effectively reduced because of the decrease in the total free energy in the Marcus inverted region, leading to the accelerated charge recombination. Such vibrationally coupled charge recombination should be a general feature of flavoproteins with similar configurations and interactions between the cofactor flavin and neighboring aromatic residues.

Oxidative damage in MauG: implications for the control of high-valent iron species and radical propagation pathways

The di-heme enzyme MauG catalyzes the oxidative biosynthesis of a tryptophan tryptophylquinone cofactor on a precursor of the enzyme methylamine dehydrogenase (preMADH). Reaction of H2O2 with the diferric form of MauG, or reaction of O2 with diferrous MauG, forms the catalytic intermediate known as bis-Fe(IV), which acts as the key oxidant during turnover. The site of substrate oxidation is more than 40 Å from the high-spin heme iron where H2O2 initially reacts, and catalysis relies on radical hopping through an interfacial residue, Trp199 of MauG. In the absence of preMADH, the bis-Fe(IV) intermediate is remarkably stable, but repeated exposure to H2O2 results in suicide inactivation. Using mass spectrometry, we show that this process involves the oxidation of three Met residues (108, 114, and 116) near the high-spin heme through ancillary electron transfer pathways engaged in the absence of substrate. The mutation of a conserved Pro107 in the distal pocket of the high-spin heme results in a dramatic increase in the level of oxidation of these Met residues. These results illustrate structural mechanisms by which MauG controls reaction with its high-valent heme cofactor and limits uncontrolled oxidation of protein residues and loss of catalytic activity. The conservation of Met residues near the high-spin heme among MauG homologues from different organisms suggests that eventual deactivation of MauG may function in a biological context. That is, methionine oxidation may represent a protective mechanism that prevents the generation of reactive oxygen species by MauG in the absence of preMADH.

The role of water-based hydrogen atom wires in long-range electron-transfer reactions in aqueous media for the FeII-FeIII self-exchange and related systems

Extended reach: A calculated mechanism for long-range proton-coupled electron transfer (PCET, see picture) through an array of structured water molecules between Fe(III)-Fe(III) complexes accounts for the reaction enthalpy and kinetic isotope effect previously measured for this reaction. This mechanism may be general and occur for other hydroxo-metal complexes.

The conserved Lys-95 charged residue cluster is critical for the homodimerization and enzyme activity of human ribonucleotide reductase small subunit M2

Ribonucleotide reductase (RR) catalyzes the reduction of ribonucleotides to deoxyribonucleotides for DNA synthesis. Human RR small subunit M2 exists in a homodimer form. However, the importance of the dimer form to the enzyme and the related mechanism remain unclear. In this study, we tried to identify the interfacial residues that may mediate the assembly of M2 homodimer by computational alanine scanning based on the x-ray crystal structure. Co-immunoprecipitation, size exclusion chromatography, and RR activity assays showed that the K95E mutation in M2 resulted in dimer disassembly and enzyme activity inhibition. In comparison, the charge-exchanging double mutation of K95E and E98K recovered the dimerization and activity. Structural comparisons suggested that a conserved cluster of charged residues, including Lys-95, Glu-98, Glu-105, and Glu-174, at the interface may function as an ionic lock for M2 homodimer. Although the measurements of the radical and iron contents showed that the monomer (the K95E mutant) was capable of generating the diiron and tyrosyl radical cofactor, co-immunoprecipitation and competitive enzyme inhibition assays indicated that the disassembly of M2 dimer reduced its interaction with the large subunit M1. In addition, the immunofluorescent and fusion protein-fluorescent imaging analyses showed that the dissociation of M2 dimer altered its subcellular localization. Finally, the transfection of the wild-type M2 but not the K95E mutant rescued the G1/S phase cell cycle arrest and cell growth inhibition caused by the siRNA knockdown of M2. Thus, the conserved Lys-95 charged residue cluster is critical for human RR M2 homodimerization, which is indispensable to constitute an active holoenzyme and function in cells.

Reversible, long-range radical transfer in E. coli class Ia ribonucleotide reductase

Ribonucleotide reductases (RNRs) catalyze the conversionof nucleotides to 2'-deoxynucleotides and are classified on the basis of the metallo-cofactor used to conduct this chemistry. The class Ia RNRs initiate nucleotide reduction when a stable diferric-tyrosyl radical (Y•, t1/2 of 4 days at 4 °C) cofactor in the β2 subunit transiently oxidizes a cysteine to a thiyl radical (S•) in the active site of the α2 subunit. In the active α2β2 complex of the class Ia RNR from E. coli , researchers have proposed that radical hopping occurs reversibly over 35 Å along a specific pathway comprised of redox-active aromatic amino acids: Y122• ↔ [W48?] ↔ Y356 in β2 to Y731 ↔ Y730 ↔ C439 in α2. Each step necessitates a proton-coupled electron transfer (PCET). Protein conformational changes constitute the rate-limiting step in the overall catalytic scheme and kinetically mask the detailed chemistry of the PCET steps. Technology has evolved to allow the site-selective replacement of the four pathway tyrosines with unnatural tyrosine analogues. Rapid kinetic techniques combined with multifrequency electron paramagnetic resonance, pulsed electron-electron double resonance, and electron nuclear double resonance spectroscopies have facilitated the analysis of stable and transient radical intermediates in these mutants. These studies are beginning to reveal the mechanistic underpinnings of the radical transfer (RT) process. This Account summarizes recent mechanistic studies on mutant E. coli RNRs containing the following tyrosine analogues: 3,4-dihydroxyphenylalanine (DOPA) or 3-aminotyrosine (NH2Y), both thermodynamic radical traps; 3-nitrotyrosine (NO2Y), a thermodynamic barrier and probe of local environmental perturbations to the phenolic pKa; and fluorotyrosines (FnYs, n = 2 or 3), dual reporters on local pKas and reduction potentials. These studies have established the existence of a specific pathway spanning 35 Å within a globular α2β2 complex that involves one stable (position 122) and three transient (positions 356, 730, and 731) Y•s. Our results also support that RT occurs by an orthogonal PCET mechanism within β2, with Y122• reduction accompanied by proton transfer from an Fe1-bound water in the diferric cluster and Y356 oxidation coupled to an off-pathway proton transfer likely involving E350. In α2, RT likely occurs by a co-linear PCET mechanism, based on studies of light-initiated radical propagation from photopeptides that mimic the β2 subunit to the intact α2 subunit and on [(2)H]-ENDOR spectroscopic analysis of the hydrogen-bonding environment surrounding a stabilized NH2Y• formed at position 730. Additionally, studies on the thermodynamics of the RT pathway reveal that the relative reduction potentials decrease according to Y122 < Y356 < Y731 ≈ Y730 ≤ C439, and that the pathway in the forward direction is thermodynamically unfavorable. C439 oxidation is likely driven by rapid, irreversible loss of water during the nucleotide reduction process. Kinetic studies of radical intermediates reveal that RT is gated by conformational changes that occur on the order of >100 s(-1) in addition to the changes that are rate-limiting in the wild-type enzyme (∼10 s(-1)). The rate constant of one of the PCET steps is ∼10(5) s(-1), as measured in photoinitiated experiments.

Geometric and electronic structure of the Mn(IV)Fe(III) cofactor in class Ic ribonucleotide reductase: correlation to the class Ia binuclear non-heme iron enzyme

The class Ic ribonucleotide reductase (RNR) from Chlamydia trachomatis (Ct) utilizes a Mn/Fe heterobinuclear cofactor, rather than the Fe/Fe cofactor found in the β (R2) subunit of the class Ia enzymes, to react with O2. This reaction produces a stable Mn(IV)Fe(III) cofactor that initiates a radical, which transfers to the adjacent α (R1) subunit and reacts with the substrate. We have studied the Mn(IV)Fe(III) cofactor using nuclear resonance vibrational spectroscopy (NRVS) and absorption (Abs)/circular dichroism (CD)/magnetic CD (MCD)/variable temperature, variable field (VTVH) MCD spectroscopies to obtain detailed insight into its geometric/electronic structure and to correlate structure with reactivity; NRVS focuses on the Fe(III), whereas MCD reflects the spin-allowed transitions mostly on the Mn(IV). We have evaluated 18 systematically varied structures. Comparison of the simulated NRVS spectra to the experimental data shows that the cofactor has one carboxylate bridge, with Mn(IV) at the site proximal to Phe127. Abs/CD/MCD/VTVH MCD data exhibit 12 transitions that are assigned as d-d and oxo and OH(-) to metal charge-transfer (CT) transitions. Assignments are based on MCD/Abs intensity ratios, transition energies, polarizations, and derivative-shaped pseudo-A term CT transitions. Correlating these results with TD-DFT calculations defines the Mn(IV)Fe(III) cofactor as having a μ-oxo, μ-hydroxo core and a terminal hydroxo ligand on the Mn(IV). From DFT calculations, the Mn(IV) at site 1 is necessary to tune the redox potential to a value similar to that of the tyrosine radical in class Ia RNR, and the OH(-) terminal ligand on this Mn(IV) provides a high proton affinity that could gate radical translocation to the α (R1) subunit.

A 2.8 Å Fe-Fe separation in the Fe2(III/IV) intermediate, X, from Escherichia coli ribonucleotide reductase

A class Ia ribonucleotide reductase (RNR) employs a μ-oxo-Fe2(III/III)/tyrosyl radical cofactor in its β subunit to oxidize a cysteine residue ~35 Å away in its α subunit; the resultant cysteine radical initiates substrate reduction. During self-assembly of the Escherichia coli RNR-β cofactor, reaction of the protein's Fe2(II/II) complex with O2 results in accumulation of an Fe2(III/IV) cluster, termed X, which oxidizes the adjacent tyrosine (Y122) to the radical (Y122(•)) as the cluster is converted to the μ-oxo-Fe2(III/III) product. As the first high-valent non-heme-iron enzyme complex to be identified and the key activating intermediate of class Ia RNRs, X has been the focus of intensive efforts to determine its structure. Initial characterization by extended X-ray absorption fine structure (EXAFS) spectroscopy yielded a Fe-Fe separation (d(Fe-Fe)) of 2.5 Å, which was interpreted to imply the presence of three single-atom bridges (O(2-), HO(-), and/or μ-1,1-carboxylates). This short distance has been irreconcilable with computational and synthetic models, which all have d(Fe-Fe) ≥ 2.7 Å. To resolve this conundrum, we revisited the EXAFS characterization of X. Assuming that samples containing increased concentrations of the intermediate would yield EXAFS data of improved quality, we applied our recently developed method of generating O2 in situ from chlorite using the enzyme chlorite dismutase to prepare X at ~2.0 mM, more than 2.5 times the concentration realized in the previous EXAFS study. The measured d(Fe-Fe) = 2.78 Å is fully consistent with computational models containing a (μ-oxo)2-Fe2(III/IV) core. Correction of the d(Fe-Fe) brings the experimental data and computational models into full conformity and informs analysis of the mechanism by which X generates Y122(•).

Hypermutability and error catastrophe due to defects in ribonucleotide reductase

The enzyme ribonucleotide reductase (RNR) plays a critical role in the production of deoxynucleoside-5'-triphosphates (dNTPs), the building blocks for DNA synthesis and replication. The levels of the cellular dNTPs are tightly controlled, in large part through allosteric control of RNR. One important reason for controlling the dNTPs relates to their ability to affect the fidelity of DNA replication and, hence, the cellular mutation rate. We have previously isolated a set of mutants of Escherichia coli RNR that are characterized by altered dNTP pools and increased mutation rates (mutator mutants). Here, we show that one particular set of RNR mutants, carrying alterations at the enzyme's allosteric specificity site, is characterized by relatively modest dNTP pool deviations but exceptionally strong mutator phenotypes, when measured in a mutational forward assay (>1,000-fold increases). We provide evidence indicating that this high mutability is due to a saturation of the DNA mismatch repair system, leading to hypermutability and error catastrophe. The results indicate that, surprisingly, even modest deviations of the cellular dNTP pools, particularly when the pool deviations promote particular types of replication errors, can have dramatic consequences for mutation rates.

Tryptophan-accelerated electron flow across a protein-protein interface

We report a new metallolabeled blue copper protein, Re126W122Cu(I) Pseudomonas aeruginosa azurin, which has three redox sites at well-defined distances in the protein fold: Re(I)(CO)3(4,7-dimethyl-1,10-phenanthroline) covalently bound at H126, a Cu center, and an indole side chain W122 situated between the Re and Cu sites (Re-W122(indole) = 13.1 Å, dmp-W122(indole) = 10.0 Å, Re-Cu = 25.6 Å). Near-UV excitation of the Re chromophore leads to prompt Cu(I) oxidation (<50 ns), followed by slow back ET to regenerate Cu(I) and ground-state Re(I) with biexponential kinetics, 220 ns and 6 μs. From spectroscopic measurements of kinetics and relative ET yields at different concentrations, it is likely that the photoinduced ET reactions occur in protein dimers, (Re126W122Cu(I))2 and that the forward ET is accelerated by intermolecular electron hopping through the interfacial tryptophan: *Re//←W122←Cu(I), where // denotes a protein-protein interface. Solution mass spectrometry confirms a broad oligomer distribution with prevalent monomers and dimers, and the crystal structure of the Cu(II) form shows two Re126W122Cu(II) molecules oriented such that redox cofactors Re(dmp) and W122-indole on different protein molecules are located at the interface at much shorter intermolecular distances (Re-W122(indole) = 6.9 Å, dmp-W122(indole) = 3.5 Å, and Re-Cu = 14.0 Å) than within single protein folds. Whereas forward ET is accelerated by hopping through W122, BET is retarded by a space jump at the interface that lacks specific interactions or water molecules. These findings on interfacial electron hopping in (Re126W122Cu(I))2 shed new light on optimal redox-unit placements required for functional long-range charge separation in protein complexes.

Targeting ribonucleotide reductase for cancer therapy

INTRODUCTION

Ribonucleotide reductase (RR) is a unique enzyme, because it is responsible for reducing ribonucleotides to their corresponding deoxyribonucleotides, which are the building blocks required for DNA replication and repair. Dysregulated RR activity is associated with genomic instability, malignant transformation and cancer development. The use of RR inhibitors, either as a single agent or combined with other therapies, has proven to be a promising approach for treating solid tumors and hematological malignancies.

AREAS COVERED

This review covers recent publications in the area of RR, which include: i) the structure, function and regulation of RR; ii) the roles of RR in cancer development; iii) the classification, mechanisms and clinical application of RR inhibitors for cancer therapy and iv) strategies for developing novel RR inhibitors in the future.

EXPERT OPINION

Exploring the possible nonenzymatic roles of RR subunit proteins in carcinogenesis may lead to new rationales for developing novel anticancer drugs. Updated information about the structure and holoenzyme models of RR will help in identifying potential sites in the protein that could be targets for novel RR inhibitors. Determining RR activity and subunit levels in clinical samples will provide a rational platform for developing personalized cancer therapies that use RR inhibitors.

A small-molecule blocking ribonucleotide reductase holoenzyme formation inhibits cancer cell growth and overcomes drug resistance

Ribonucleotide reductase (RNR) is an attractive target for anticancer agents given its central function in DNA synthesis, growth, metastasis, and drug resistance of cancer cells. The current clinically established RNR inhibitors have the shortcomings of short half-life, drug resistance, and iron chelation. Here, we report the development of a novel class of effective RNR inhibitors addressing these issues. A novel ligand-binding pocket on the RNR small subunit (RRM2) near the C-terminal tail was proposed by computer modeling and verified by site-directed mutagenesis and nuclear magnetic resonance (NMR) techniques. A compound targeting this pocket was identified by virtual screening of the National Cancer Institute (NCI) diverse small-molecule database. By lead optimization, we developed the novel RNR inhibitor COH29 that acted as a potent inhibitor of both recombinant and cellular human RNR enzymes. COH29 overcame hydroxyurea and gemcitabine resistance in cancer cells. It effectively inhibited proliferation of most cell lines in the NCI 60 human cancer panel, most notably ovarian cancer and leukemia, but exerted little effect on normal fibroblasts or endothelial cells. In mouse xenograft models of human cancer, COH29 treatment reduced tumor growth compared with vehicle. Site-directed mutagenesis, NMR, and surface plasmon resonance biosensor studies confirmed COH29 binding to the proposed ligand-binding pocket and offered evidence for assembly blockade of the RRM1-RRM2 quaternary structure. Our findings offer preclinical validation of COH29 as a promising new class of RNR inhibitors with a new mechanism of inhibition, with broad potential for improved treatment of human cancer.

Structure of a dinuclear iron cluster-containing β-hydroxylase active in antibiotic biosynthesis

A family of dinuclear iron cluster-containing oxygenases that catalyze β-hydroxylation tailoring reactions in natural product biosynthesis by nonribosomal peptide synthetase (NRPS) systems was recently described [Makris, T. M., Chakrabarti, M., Münck, E., and Lipscomb, J. D. (2010) Proc. Natl. Acad. Sci. U.S.A. 107, 15391-15396]. Here, the 2.17 Å X-ray crystal structure of the archetypal enzyme from the family, CmlA, is reported. CmlA catalyzes β-hydroxylation of l-p-aminophenylalanine during chloramphenicol biosynthesis. The fold of the N-terminal domain of CmlA is unlike any previously reported, but the C-terminal domain has the αββα fold of the metallo-β-lactamase (MBL) superfamily. The diiron cluster bound in the C-terminal domain is coordinated by an acetate, three His residues, two Asp residues, one Glu residue, and a bridging oxo moiety. One of the Asp ligands forms an unusual monodentate bridge. No other oxygen-activating diiron enzyme utilizes this ligation or the MBL protein fold. The N-terminal domain facilitates dimerization, but using computational docking and a sequence-based structural comparison to homologues, we hypothesize that it likely serves additional roles in NRPS recognition and the regulation of O2 activation.

Uncoupling of allosteric and oligomeric regulation in a functional hybrid enzyme constructed from Escherichia coli and human ribonucleotide reductase

An N-terminal-domain (NTD) and adjacent catalytic body (CB) make up subunit-α of ribonucleotide reductase (RNR), the rate-limiting enzyme for de novo dNTP biosynthesis. A strong linkage exists between ligand binding at the NTD and oligomerization-coupled RNR inhibition, inducible by both dATP and nucleotide chemotherapeutics. These observations have distinguished the NTD as an oligomeric regulation domain dictating the assembly of inactive RNR oligomers. Inactive states of RNR differ between eukaryotes and prokaryotes (α6 in human versus α4β4 in Escherichia coli , wherein β is RNR's other subunit); however, the NTD structurally interconnects individual α2 or α2 and β2 dimeric motifs within the respective α6 or α4β4 complexes. To elucidate the influence of NTD ligand binding on RNR allosteric and oligomeric regulation, we engineered a human- E. coli hybrid enzyme (HE) where human-NTD is fused to E. coli -CB. Both the NTD and the CB of the HE bind dATP. The HE specifically partners with E. coli -β to form an active holocomplex. However, although the NTD is the sole physical tether to support α2 and/or β2 associations in the dATP-bound α6 or α4β4 fully inhibited RNR complexes, the binding of dATP to the HE NTD only partially suppresses HE activity and fully precludes formation of higher-order HE oligomers. We postulate that oligomeric regulation is the ultimate mechanism for potent RNR inhibition, requiring species-specific NTD-CB interactions. Such interdomain cooperativity in RNR oligomerization is unexpected from structural studies alone or biochemical studies of point mutants.

A novel clade of unique eukaryotic ribonucleotide reductase R2 subunits is exclusive to apicomplexan parasites

Apicomplexa are protist parasites of tremendous medical and economic importance, causing millions of deaths and billions of dollars in losses each year. Apicomplexan-related diseases may be controlled via inhibition of essential enzymes. Ribonucleotide reductase (RNR) provides the only de novo means of synthesizing deoxyribonucleotides, essential precursors for DNA replication and repair. RNR has long been the target of antibacterial and antiviral therapeutics. However, targeting this ubiquitous protein in eukaryotic pathogens may be problematic unless these proteins differ significantly from that of their respective host. The typical eukaryotic RNR enzymes belong to class Ia, and the holoenzyme consists minimally of two R1 and two R2 subunits (α₂β₂). We generated a comparative, annotated, structure-based, multiple-sequence alignment of R2 subunits, identified a clade of R2 subunits unique to Apicomplexa, and determined its phylogenetic position. Our analyses revealed that the apicomplexan-specific sequences share characteristics with both class I R2 and R2lox proteins. The putative radical-harboring residue, essential for the reduction reaction by class Ia R2-containing holoenzymes, was not conserved within this group. Phylogenetic analyses suggest that class Ia subunits are not monophyletic and consistently placed the apicomplexan-specific clade sister to the remaining class Ia eukaryote R2 subunits. Our research suggests that the novel apicomplexan R2 subunit may be a promising candidate for chemotherapeutic-induced inhibition as it differs greatly from known eukaryotic host RNRs and may be specifically targeted.

Structural basis for assembly of the Mn(IV)/Fe(III) cofactor in the class Ic ribonucleotide reductase from Chlamydia trachomatis

The class Ic ribonucleotide reductase (RNR) from Chlamydia trachomatis (Ct) employs a Mn(IV)/Fe(III) cofactor in each monomer of its β2 subunit to initiate nucleotide reduction. The cofactor forms by reaction of Mn(II)/Fe(II)-β2 with O2. Previously, in vitro cofactor assembly from apo β2 and divalent metal ions produced a mixture of two forms, with Mn at site 1 (Mn(IV)/Fe(III)) or site 2 (Fe(III)/Mn(IV)), of which the more active Mn(IV)/Fe(III) product predominates. Here we have addressed the basis for metal site selectivity by determining X-ray crystal structures of apo, Mn(II), and Mn(II)/Fe(II) complexes of Ct β2. A structure obtained anaerobically with equimolar Mn(II), Fe(II), and apoprotein reveals exclusive incorporation of Mn(II) at site 1 and Fe(II) at site 2, in contrast to the more modest site selectivity achieved previously. Site specificity is controlled thermodynamically by the apoprotein structure, as only minor adjustments of ligands occur upon metal binding. Additional structures imply that, by itself, Mn(II) binds in either site. Together, the structures are consistent with a model for in vitro cofactor assembly in which Fe(II) specificity for site 2 drives assembly of the appropriately configured heterobimetallic center, provided that Fe(II) is substoichiometric. This model suggests that use of a Mn(IV)/Fe(III) cofactor in vivo could be an adaptation to Fe(II) limitation. A 1.8 Å resolution model of the Mn(II)/Fe(II)-β2 complex reveals additional structural determinants for activation of the cofactor, including a proposed site for side-on (η(2)) addition of O2 to Fe(II) and a short (3.2 Å) Mn(II)-Fe(II) interionic distance, promoting formation of the Mn(IV)/Fe(IV) activation intermediate.

Modulation of Y356 photooxidation in E. coli class Ia ribonucleotide reductase by Y731 across the α2:β2 interface

Substrate turnover in class Ia ribonucleotide reductase (RNR) requires reversible radical transport across two subunits over 35 Å, which occurs by a multistep proton-coupled electron-transfer mechanism. Using a photooxidant-labeled β2 subunit of Escherichia coli class Ia RNR, we demonstrate photoinitiated oxidation of a tyrosine in an α2:β2 complex, which results in substrate turnover. Using site-directed mutations of the redox-active tyrosines at the subunit interface, Y356F(β) and Y731F(α), this oxidation is identified to be localized on Y356. The rate of Y356 oxidation depends on the presence of Y731 across the interface. This observation supports the proposal that unidirectional PCET across the Y356(β)-Y731(α)-Y730(α) triad is crucial to radical transport in RNR.

Distance-independent charge recombination kinetics in cytochrome c-cytochrome c peroxidase complexes: compensating changes in the electronic coupling and reorganization energies

Charge recombination rate constants vary no more than 3-fold for interprotein ET in the Zn-substituted wild type (WT) cytochrome c peroxidase (CcP):cytochrome c (Cc) complex and in complexes with four mutants of the Cc protein (i.e., F82S, F82W, F82Y, and F82I), despite large differences in the ET distance. Theoretical analysis indicates that charge recombination for all complexes involves a combination of tunneling and hopping via Trp191. For three of the five structures (WT and F82S(W)), the protein favors hopping more than that in the other two structures that have longer heme → ZnP distances (F82Y(I)). Experimentally observed biexponential ET kinetics is explained by the complex locking in alternative coupling pathways, where the acceptor hole state is either primarily localized on ZnP (slow phase) or on Trp191 (fast phase). The large conformational differences between the CcP:Cc interface for the F82Y(I) mutants compared to that the WT and F82S(W) complexes are predicted to change the reorganization energies for the CcP:Cc ET reactions because of changes in solvent exposure and interprotein ET distances. Since the recombination reaction is likely to occur in the inverted Marcus regime, an increased reorganization energy compensates the decreased role for hopping recombination (and the longer transfer distance) in the F82Y(I) mutants. Taken together, coupling pathway and reorganization energy effects for the five protein complexes explain the observed insensitivity of recombination kinetics to donor-acceptor distance and docking pose and also reveals how hopping through aromatic residues can accelerate long-range ET.

An Iron(II)(1,3-bis(2'-pyridylimino)isoindoline) Complex as a Catalyst for Substrate Oxidation with H2O2. Evidence for a Transient Peroxodiiron(III) Species

The complex [Fe(indH)(solvent)3](ClO4)2 (1) has been isolated from the reaction of equimolar amounts of 1,3-bis(2'-pyridylimino)isoindoline (indH) and Fe(ClO4)2 in acetonitrile and characterized by X-ray crystallography and several spectroscopic techniques. It is a suitable catalyst for the oxidation of thioanisoles and benzyl alcohols with H2O2 as the oxidant. Hammett correlations and kinetic isotope effect experiments support the involvement of an electrophilic metal-based oxidant. A metastable green species (2) is observed when 1 is reacted with H2O2 at -40 °C, which has been characterized to have a FeIII(μ-O)(μ-O2)FeIII core on the basis of UV-Vis, electron paramagnetic resonance, resonance Raman, and X-ray absorption spectroscopic data.