Density functional theory study of Fe(IV) d-d optical transitions in active-site models of class I ribonucleotide reductase intermediate X with vertical self-consistent reaction field methods

Theoretical Insights into the Generation Mechanism of the Tyr122 Radical Catalyzed by Intermediate X in Class Ia Ribonucleotide Reductase

Ribonucleotide reductase (RNR) catalyzes the reduction of ribonucleotides to deoxyribonucleotides in all organisms. There is an ∼35 Å long-range electron-hole transfer pathway during the catalytic process of class Ia RNR, which can be described as Tyr122β ↔ [Trp48β]? ↔ Tyr356β ↔ Tyr731α ↔ Tyr730α ↔ Cys439α. The formation of the Y122• radical initiates this long-range radical transfer process. However, the generation mechanism of Y122• is not yet clear due to confusion over the intermediate X structures. Based on the two reported X structures, we examined the possible mechanisms of Y122• generation by density functional theory (DFT) calculations. Our examinations revealed that the generation of the Y122• radical from the two different core structures of X was via a similar two-step reaction, with the first step of proton transfer for the formation of the proton receptor of Y122 and the second step of a proton-coupled long-range electron transfer reaction with the proton transfer from the Y122 hydroxyl group to the terminal hydroxide ligand of Fe1III and simultaneously electron transfer from the side chain of Y122 to Fe2IV. These findings provide an insight into the formation mechanism of Y122• catalyzed by the double-iron center of the β subunit of class Ia RNR.

Dioxygen Activation by Nonheme Diiron Enzymes: Diverse Dioxygen Adducts, High-Valent Intermediates, and Related Model Complexes

A growing subset of metalloenzymes activates dioxygen with nonheme diiron active sites to effect substrate oxidations that range from the hydroxylation of methane and the desaturation of fatty acids to the deformylation of fatty aldehydes to produce alkanes and the six-electron oxidation of aminoarenes to nitroarenes in the biosynthesis of antibiotics. A common feature of their reaction mechanisms is the formation of O2 adducts that evolve into more reactive derivatives such as diiron(II,III)-superoxo, diiron(III)-peroxo, diiron(III,IV)-oxo, and diiron(IV)-oxo species, which carry out particular substrate oxidation tasks. In this review, we survey the various enzymes belonging to this unique subset and the mechanisms by which substrate oxidation is carried out. We examine the nature of the reactive intermediates, as revealed by X-ray crystallography and the application of various spectroscopic methods and their associated reactivity. We also discuss the structural and electronic properties of the model complexes that have been found to mimic salient aspects of these enzyme active sites. Much has been learned in the past 25 years, but key questions remain to be answered.

Geometrical properties of the manganese(iv)/iron(iii) cofactor of Chlamydia trachomatis ribonucleotide reductase unveiled by simulations of XAS spectra

A combination of EXAFS simulations and DFT calculations, including a novel protocol to evaluate Debye–Waller factors, provide insights into the structure of the Mn(iv)/Fe(iii) cofactor ofCtR2.

Composition and Structure of the Inorganic Core of Relaxed Intermediate X(Y122F) of Escherichia coli Ribonucleotide Reductase

Activation of the diferrous center of the β2 (R2) subunit of the class 1a Escherichia coli ribonucleotide reductases by reaction with O2 followed by one-electron reduction yields a spin-coupled, paramagnetic Fe(III)/Fe(IV) intermediate, denoted X, whose identity has been sought by multiple investigators for over a quarter of a century. To determine the composition and structure of X, the present study has applied (57)Fe, (14,15)N, (17)O, and (1)H electron nuclear double resonance (ENDOR) measurements combined with quantitative measurements of (17)O and (1)H electron paramagnetic resonance line-broadening studies to wild-type X, which is very short-lived, and to X prepared with the Y122F mutant, which has a lifetime of many seconds. Previous studies have established that over several seconds the as-formed X(Y122F) relaxes to an equilibrium structure. The present study focuses on the relaxed structure. It establishes that the inorganic core of relaxed X has the composition [(OH(-))Fe(III)-O-Fe(IV)]: there is no second inorganic oxygenic bridge, neither oxo nor hydroxo. Geometric analysis of the (14)N ENDOR data, together with recent extended X-ray absorption fine structure measurements of the Fe-Fe distance (Dassama, L. M.; et al. J. Am. Chem. Soc. 2013, 135, 16758), supports the view that X contains a "diamond-core" Fe(III)/Fe(IV) center, with the irons bridged by two ligands. One bridging ligand is the oxo bridge (OBr) derived from O2 gas. Given the absence of a second inorganic oxygenic bridge, the second bridging ligand must be protein derived, and is most plausibly assigned as a carboxyl oxygen from E238.

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.

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(•).

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, Mö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.

Density functional theory analysis of structure, energetics, and spectroscopy for the Mn-Fe active site of Chlamydia trachomatis ribonucleotide reductase in four oxidation states

Models for the Mn-Fe active site structure of ribonucleotide reductase (RNR) from pathogenic bacteria Chlamydia trachomatis (Ct) in different oxidation states have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and also with finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. The detailed structures for the reduced Mn(II)-Fe(II), the met Mn(III)-Fe(III), the oxidized Mn(IV)-Fe(III) and the superoxidized Mn(IV)-Fe(IV) states are predicted. The calculated properties, including geometries, (57)Fe Mossbauer isomer shifts and quadrupole splittings, and (57)Fe and (55)Mn electron nuclear double resonance (ENDOR) hyperfine coupling constants, are compared with the available experimental data. The Mössbauer and energetic calculations show that the (mu-oxo, mu-hydroxo) models better represent the structure of the Mn(IV)-Fe(III) state than the di-mu-oxo models. The predicted Mn(IV)-Fe(III) distances (2.95 and 2.98 A) in the (mu-oxo, mu-hydroxo) models are in agreement with the extended X-ray absorption fine structure (EXAFS) experimental value of 2.92 A (Younker et al. J. Am. Chem. Soc. 2008, 130, 15022-15027). The effect of the protein and solvent environment on the assignment of the Mn metal position is examined by comparing the relative energies of alternative mono-Mn(II) active site structures. It is proposed that if the Mn(II)-Fe(II) protein is prepared with prior addition of Mn(II) or with Mn(II) richer than Fe(II), Mn is likely positioned at metal site 2, which is further from Phe127.

Calibration of scalar relativistic density functional theory for the calculation of sulfur K-edge X-ray absorption spectra

Sulfur K-edge X-ray absorption spectroscopy has been proven to be a powerful tool for investigating the electronic structures of sulfur-containing coordination complexes. The full information content of the spectra can be developed through a combination of experiment and time-dependent density functional theory (TD-DFT). In this work, the necessary calibration is carried out for a range of contemporary functionals (BP86, PBE, OLYP, OPBE, B3LYP, PBE0, TPSSh) in a scalar relativistic (0(th) order regular approximation, ZORA) DFT framework. It is shown that with recently developed segmented all-electron scalar relativistic (SARC) basis sets one obtains results that are as good as with large, uncontracted basis sets. The errors in the calibrated transition energies are on the order of 0.1 eV. The error in calibrated intensities is slightly larger, but the calculations are still in excellent agreement with experiment. The behavior of full TD-DFT linear response versus the Tamm-Dancoff approximation has been evaluated with the result that two methods are almost indistinguishable. The inclusion of relativistic effects barely changes the results for first row transition metal complexes, however, the contributions become visible for second-row transition metals and reach a maximum (of an approximately 10% change in the calibration parameters) for third row transition metal species. The protocol developed here is approximately 10 times more efficient than the previously employed protocol, which was based on large, uncontracted basis sets. The calibration strategy followed here may be readily extended to other edges.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

DFT calculations of comparative energetics and ENDOR/Mössbauer properties for two protonation states of the iron dimer cluster of ribonucleotide reductase intermediate X

Two models (I and II) for the active site structure of class-I ribonucleotide reductase (RNR) intermediate X in subunit R2 have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and with the finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. Only one of the bridging groups between the two iron centers is different between model-I and model-II. Model-I contains two mu-oxo bridges, while model-II has one bridging oxo and one bridging hydroxo. These are large active site models including up to the fourth coordination shell H-bonding residues. Mössbauer and ENDOR hyperfine property calculations show that model-I is more likely to represent the active site structure of RNR-X. However, energetically our pK(a) calculations at first highly favored the bridging oxo and hydroxo (in model-II) structure of the diiron center rather than having the di-oxo bridge (in model-I). Since the Arg236 and the nearby Lys42, which are very close to the diiron center, are on the protein surface of RNR-R2, it is highly feasible that one or two anion groups in solution would interact with the positively charged side chains of Arg236 and Lys42. The anion group(s) can be a reductant, phosphate, sulfate, nitrate, and other negatively charged groups existing in biological environments or in the buffer of the experiment. Since sulfate ions certainly exist in the buffer of the ENDOR experiment, we have examined the effect of the sulfate (SO(4)(2-), surrounded by explicit water molecules) H-bonding to the side chain of Arg236. We find that when sulfate interacts with Arg236, the carboxylate group of Asp237 tends to be protonated, and once Asp237 is protonated, the Fe(iii)Fe(iv) center in X favors the di-oxo bridge (model-I). This would explain that the ENDOR observed RNR-X active site structure is likely to be represented by model-I rather than model-II.

Quantum cluster size and solvent polarity effects on the geometries and Mössbauer properties of the active site model for ribonucleotide reductase intermediate X: a density functional theory study

In studying the properties of metalloproteins using ab initio quantum mechanical methods, one has to focus on the calculations on the active site. The bulk protein and solvent environment is often neglected, or is treated as a continuum dielectric medium with a certain dielectric constant. The size of the quantum cluster of the active site chosen for calculations can vary by including only the first-shell ligands which are directly bound to the metal centers, or including also the second-shell residues which are adjacent to and normally have H-bonding interactions with the first-shell ligands, or by including also further hydrogen bonding residues. It is not well understood how the size of the quantum cluster and the value of the dielectric constant chosen for the calculations will influence the calculated properties. In this paper, we have studied three models (A, B, and C) of different sizes for the active site of the ribonucleotide reductase intermediate X, using density functional theory (DFT) OPBE functional with broken-symmetry methodology. Each model is studied in gas-phase and in the conductor-like screening (COSMO) solvation model with different dielectric constants ε = 4, 10, 20, and 80, respectively. All the calculated Fe-ligand geometries, Heisenberg J coupling constants, and the Mössbauer isomer shifts, quadrupole splittings, and the (57)Fe, (1)H, and (17)O hyperfine tensors are compared. We find that the calculated isomer shifts are very stable. They are virtually unchanged with respect to the size of the cluster and the dielectric constant of the environment. On the other hand, certain Fe-ligand distances are sensitive to both the size of the cluster and the value of ε. ε = 4, which is normally used for the protein environment, appears too small when studying the diiron active site geometry with only the first-shell ligands as seen by comparisons with larger models.

Identification of protonated oxygenic ligands of ribonucleotide reductase intermediate X

We previously used a combination of continuous-wave (CW) and pulsed electron-nuclear double resonance (ENDOR) protocols to identify the types of protonated oxygen (OH(x)) species and their disposition within the Fe(III)/Fe(IV) cluster of intermediate X, the direct precursor of the essential diferric-tyrosyl radical cofactor of the beta2 subunit of Escherichia coli ribonucleotide reductase (RNR). We concluded that X contains the [(H(x)O)Fe(III)OFe(IV)] fragment (T model), and does not contain a mu-hydroxo bridge. When combined with a subsequent (17)O ENDOR study of X prepared with H(2)(17)O and with (17)O(2), the results led us to suggest that this fragment is the entire inorganic core of X. This has been questioned by recent reports, but these reports do not themselves agree on the core of X. One concluded that X possesses a di-mu-oxo Fe(III)/Fe(IV) core plus a terminal (H(2)O) bound to Fe(III) [e.g., Han, W.-G.; Liu, T.; Lovell, T.; Noodleman, L. J. Am. Chem. Soc. 2005, 127, 15778-15790]. The other [Mitic, N.; Clay, M. D.; Saleh, L.; Bollinger, J. M.; Solomon, E. I. J. Am. Chem. Soc. 2007, 129, 9049-9065] concluded that X contains only a single oxo bridge and postulated the presence of an additional hydroxo bridge plus a terminal hydroxyl bound to Fe(III). In this report we take advantage of improvements in 35 GHz pulsed ENDOR performance to reexamine the protonation state of oxygenic ligands of the inorganic core of X by directly probing the exchangeable proton(s) with (2)H pulsed ENDOR spectroscopy. These (2)H ENDOR measurements confirm that X contains an Fe(III)-bound terminal aqua ligand (H(x)O), but the spectra contain none of the features that would be required for the proton of a bridging hydroxyl. Thus, we confirm that X contains a terminal aqua (most likely hydroxo) ligand to Fe(III) in addition to one or two mu-oxo bridges but does not contain a mu-hydroxo bridge. The (2)H ENDOR measurements further demonstrate that this conclusion is applicable to both wild type and Y122F-beta2 mutant, and in fact we detect no difference between the properties of protons on the terminal oxygens in the two variants; likewise, (14)N ENDOR measurements of histidyl ligands bound to Fe show no difference between the two variants.

Formation and function of the Manganese(IV)/Iron(III) cofactor in Chlamydia trachomatis ribonucleotide reductase

The beta(2) subunit of a class Ia or Ib ribonucleotide reductase (RNR) is activated when its carboxylate-bridged Fe(2)(II/II) cluster reacts with O(2) to oxidize a nearby tyrosine (Y) residue to a stable radical (Y(*)). During turnover, the Y(*) in beta(2) is thought to reversibly oxidize a cysteine (C) in the alpha(2) subunit to a thiyl radical (C(*)) by a long-distance ( approximately 35 A) proton-coupled electron-transfer (PCET) step. The C(*) in alpha(2) then initiates reduction of the 2' position of the ribonucleoside 5'-diphosphate substrate by abstracting the hydrogen atom from C3'. The class I RNR from Chlamydia trachomatis (Ct) is the prototype of a newly recognized subclass (Ic), which is characterized by the presence of a phenylalanine (F) residue at the site of beta(2) where the essential radical-harboring Y is normally found. We recently demonstrated that Ct RNR employs a heterobinuclear Mn(IV)/Fe(III) cluster for radical initiation. In essence, the Mn(IV) ion of the cluster functionally replaces the Y(*) of the conventional class I RNR. The Ct beta(2) protein also autoactivates by reaction of its reduced (Mn(II)/Fe(II)) metal cluster with O(2). In this reaction, an unprecedented Mn(IV)/Fe(IV) intermediate accumulates almost stoichiometrically and decays by one-electron reduction of the Fe(IV) site. This reduction is mediated by the near-surface residue, Y222, a residue with no functional counterpart in the well-studied conventional class I RNRs. In this review, we recount the discovery of the novel Mn/Fe redox cofactor in Ct RNR and summarize our current understanding of how it assembles and initiates nucleotide reduction.

Structural analysis of the Mn(IV)/Fe(III) cofactor of Chlamydia trachomatis ribonucleotide reductase by extended X-ray absorption fine structure spectroscopy and density functional theory calculations

The class Ic ribonucleotide reductase from Chlamydia trachomatis ( Ct) uses a stable Mn(IV)/Fe(III) cofactor to initiate nucleotide reduction by a free-radical mechanism. Extended X-ray absorption fine structure (EXAFS) spectroscopy and density functional theory (DFT) calculations are used to postulate a structure for this cofactor. Fe and Mn K-edge EXAFS data yield an intermetallic distance of approximately 2.92 A. The Mn data also suggest the presence of a short 1.74 A Mn-O bond. These metrics are compared to the results of DFT calculations on 12 cofactor models derived from the crystal structure of the inactive Fe 2(III/III) form of the protein. Models are differentiated by the protonation states of their bridging and terminal OH X ligands as well as the location of the Mn(IV) ion (site 1 or 2). The models that agree best with experimental observation feature a mu-1,3-carboxylate bridge (E120), terminal solvent (H 2O/OH) to site 1, one mu-O bridge, and one mu-OH bridge. The site-placement of the metal ions cannot be discerned from the available data.

Importance of the maintenance pathway in the regulation of the activity of Escherichia coli ribonucleotide reductase

Ribonucleotide reductases (RNRs) catalyze the conversion of nucleotides to deoxynucleotides in all organisms. The Escherichia coli class Ia RNR is composed of α and β subunits that form an α₂β₂ active complex. β contains the diferric tyrosyl radical (Y^•) cofactor that is essential for the reduction process that occurs on α. [Y^•] in vitro is proportional to RNR activity, and its regulation in vivo potentially represents a mechanism for controlling RNR activity. To examine this thesis, N- and C-terminal StrepII-tagged β under the control of an l-arabinose promoter were constructed. Using these constructs and with [l-arabinose] varying from 0 to 0.5 mM in the growth medium, [β] could be varied from 4 to 3300 µM. [Y^•] in vivo and on affinity-purified Strep-β in vitro was determined by EPR spectroscopy and Western analysis. In both cases, there was 0.1–0.3 Y^• radical per β. To determine if the substoichiometric Y^• level was associated with apo β or diferric β, titrations of crude cell extracts from these growths were carried out with reduced YfaE, a 2Fe2S ferredoxin involved in cofactor maintenance and assembly. Each titration, followed by addition of O₂ to assemble the cofactor and EPR analysis to quantitate Y^•, revealed that β is completely loaded with a diferric cluster even when its concentration in vivo is 244 µM. These titrations, furthermore, resulted in 1 Y^• radical per β, the highest levels reported. Whole cell Mössbauer analysis on cells induced with 0.5 mM arabinose supports high iron loading in β. These results suggest that modulation of the level of Y^• in vivo in E. coli is a mechanism of regulating RNR activity.

Spectroscopic and electronic structure studies of intermediate X in ribonucleotide reductase R2 and two variants: a description of the FeIV-oxo bond in the FeIII-O-FeIV dimer

Spectroscopic and electronic structure studies of the class I Escherichia coli ribonucleotide reductase (RNR) intermediate X and three computationally derived model complexes are presented, compared, and evaluated to determine the electronic and geometric structure of the FeIII-FeIV active site of intermediate X. Rapid freeze-quench (RFQ) EPR, absorption, and MCD were used to trap intermediate X in R2 wild-type (WT) and two variants, W48A and Y122F/Y356F. RFQ-EPR spin quantitation was used to determine the relative contributions of intermediate X and radicals present, while RFQ-MCD was used to specifically probe the FeIII/FeIV active site, which displayed three FeIV d-d transitions between 16,700 and 22,600 cm(-1), two FeIV d-d spin-flip transitions between 23,500 and 24,300 cm(-1), and five oxo to FeIV and FeIII charge transfer (CT) transitions between 25,000 and 32,000 cm(-1). The FeIV d-d transitions were perturbed in the two variants, confirming that all three d-d transitions derive from the d-pi manifold. Furthermore, the FeIV d-pi splittings in the WT are too large to correlate with a bis-mu-oxo structure. The assignment of the FeIV d-d transitions in WT intermediate X best correlates with a bridged mu-oxo/mu-hydroxo [FeIII(mu-O)(mu-OH)FeIV] structure. The mu-oxo/mu-hydroxo core structure provides an important sigma/pi superexchange pathway, which is not present in the bis-mu-oxo structure, to promote facile electron transfer from Y122 to the remote FeIV through the bent oxo bridge, thereby generating the tyrosyl radical for catalysis.