Density Functional Theory Study of the Manganese-Containing Ribonucleotide Reductase from Chlamydia trachomatis: Why Manganese Is Needed in the Active Complex

Mechanistic Insights into the N-Hydroxylations Catalyzed by the Binuclear Iron Domain of SznF Enzyme: Key Piece in the Synthesis of Streptozotocin

SznF, a member of the emerging family of heme-oxygenase-like (HO-like) di-iron oxidases and oxygenases, employs two distinct domains to catalyze the conversion of Nω-methyl-L-arginine (L-NMA) into N-nitroso-containing product, which can subsequently be transformed into streptozotocin. Using unrestricted density functional theory (UDFT) with the hybrid functional B3LYP, we have mechanistically investigated the two sequential hydroxylations of L-NMA catalyzed by SznF's binuclear iron central domain. Mechanism B primarily involves the O-O bond dissociation, forming Fe(IV)=O, induced by the H+/e- introduction to the FeA side of μ-1,2-peroxo-Fe2(III/III), the substrate hydrogen abstraction by Fe(IV)=O, and the hydroxyl rebound to the substrate N radical. The stochastic addition of H+/e- to the FeB side (mechanism C) can transition to mechanism B, thereby preventing enzyme deactivation. Two other competing mechanisms, involving the direct O-O bond dissociation (mechanism A) and the addition of H2O as a co-substrate (mechanism D), have been ruled out.

Why is manganese so valuable to bacterial pathogens?

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

Interchangeable utilization of metals: New perspectives on the impacts of metal ions employed in ancient and extant biomolecules

Metal ions provide considerable functionality across biological systems, and their utilization within biomolecules has adapted through changes in the chemical environment to maintain the activity they facilitate. While ancient earth's atmosphere was rich in iron and manganese and low in oxygen, periods of atmospheric oxygenation significantly altered the availability of certain metal ions, resulting in ion replacement within biomolecules. This adaptation mechanism has given rise to the phenomenon of metal cofactor interchangeability, whereby contemporary proteins and nucleic acids interact with multiple metal ions interchangeably, with different coordinated metals influencing biological activity, stability, and toxic potential. The ability of extant organisms to adapt to fluctuating metal availability remains relevant in a number of crucial biomolecules, including the superoxide dismutases of the antioxidant defense systems and ribonucleotide reductases. These well-studied and ancient enzymes illustrate the potential for metal interchangeability and adaptive utilization. More recently, the ribosome has also been demonstrated to exhibit interchangeable interactions with metal ions with impacts on function, stability, and stress adaptation. Using these and other examples, here we review the biological significance of interchangeable metal ions from a new angle that combines both biochemical and evolutionary viewpoints. The geochemical pressures and chemical properties that underlie biological metal utilization are discussed in the context of their impact on modern disease states and treatments.

O–H and (CO)N–H bond weakening by coordination to Fe(ii)

Coordination of hydroxyl/amide groups to Fe(ii) diminishes BDFEs of O–H and (CO)N–H bonds down to 76.0 and 80.5 kcal mol⁻¹ respectively.

Theoretical Study of the Mechanism of the Nonheme Iron Enzyme EgtB

EgtB is a nonheme iron enzyme catalyzing the C-S bond formation between γ-glutamyl cysteine (γGC) and N-α-trimethyl histidine (TMH) in the ergothioneine biosynthesis. Density functional calculations were performed to elucidate and delineate the reaction mechanism of this enzyme. Two different mechanisms were considered, depending on whether the sulfoxidation or the S-C bond formation takes place first. The calculations suggest that the S-O bond formation occurs first between the thiolate and the ferric superoxide, followed by homolytic O-O bond cleavage, very similar to the case of cysteine dioxygenase. Subsequently, proton transfer from a second-shell residue Tyr377 to the newly generated iron-oxo moiety takes place, which is followed by proton transfer from the TMH imidazole to Tyr377, facilitated by two crystallographically observed water molecules. Next, the S-C bond is formed between γGC and TMH, followed by proton transfer from the imidazole CH moiety to Tyr377, which was calculated to be the rate-limiting step for the whole reaction, with a barrier of 17.9 kcal/mol in the quintet state. The calculated barrier for the rate-limiting step agrees quite well with experimental kinetic data. Finally, this proton is transferred back to the imidazole nitrogen to form the product. The alternative thiyl radical attack mechanism has a very high barrier, being 25.8 kcal/mol, ruling out this possibility.

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.

Protonation State of MnFe and FeFe Cofactors in a Ligand-Binding Oxidase Revealed by X-ray Absorption, Emission, and Vibrational Spectroscopy and QM/MM Calculations

Enzymes with a dimetal-carboxylate cofactor catalyze reactions among the top challenges in chemistry such as methane and dioxygen (O₂) activation. Recently described proteins bind a manganese-iron cofactor (MnFe) instead of the classical diiron cofactor (FeFe). Determination of atomic-level differences of homo- versus hetero-bimetallic cofactors is crucial to understand their diverse redox reactions. We studied a ligand-binding oxidase from the bacterium Geobacillus kaustophilus (R2lox) loaded with a FeFe or MnFe cofactor, which catalyzes O₂ reduction and an unusual tyrosine-valine ether cross-link formation, as revealed by X-ray crystallography. Advanced X-ray absorption, emission, and vibrational spectroscopy methods and quantum chemical and molecular mechanics calculations provided relative Mn/Fe contents, X-ray photoreduction kinetics, metal-ligand bond lengths, metal-metal distances, metal oxidation states, spin configurations, valence-level degeneracy, molecular orbital composition, nuclear quadrupole splitting energies, and vibrational normal modes for both cofactors. A protonation state with an axial water (H₂O) ligand at Mn or Fe in binding site 1 and a metal-bridging hydroxo group (μOH) in a hydrogen-bonded network is assigned. Our comprehensive picture of the molecular, electronic, and dynamic properties of the cofactors highlights reorientation of the unique axis along the Mn-OH₂ bond for the Mn1(III) Jahn-Teller ion but along the Fe-μOH bond for the octahedral Fe1(III). This likely corresponds to a more positive redox potential of the Mn(III)Fe(III) cofactor and higher proton affinity of its μOH group. Refined model structures for the Mn(III)Fe(III) and Fe(III)Fe(III) cofactors are presented. Implications of our findings for the site-specific metalation of R2lox and performance of the O₂ reduction and cross-link formation reactions are discussed.

Mechanism and selectivity of the dinuclear iron benzoyl-coenzyme A epoxidase BoxB

DFT calculations are used to elucidate the reaction mechanism and selectivity of BoxB catalyzed benzoyl-CoA epoxidation.

Benzoyl-CoA epoxidase is a dinuclear iron enzyme that catalyzes the epoxidation reaction of the aromatic ring of benzoyl-CoA with chemo-, regio- and stereo-selectivity. It has been suggested that this enzyme may also catalyze the deoxygenation reaction of epoxide, suggesting a unique bifunctionality among the diiron enzymes. We report a density functional theory study of this enzyme aimed at elucidating its mechanism and the various selectivities. The epoxidation is suggested to start with the binding of the O₂ molecule to the diferrous center to generate a diferric peroxide complex, followed by concerted O–O bond cleavage and epoxide formation. Two different pathways have been located, leading to (2S,3R)-epoxy and (2R,3S)-epoxy products, with barriers of 17.6 and 20.4 kcal mol^–1, respectively. The barrier difference is 2.8 kcal mol^–1, corresponding to a diastereomeric excess of about 99 : 1. Further isomerization from epoxide to phenol is found to have quite a high barrier, which cannot compete with the product release step. After product release into solution, fast epoxide–oxepin isomerization and racemization can take place easily, leading to a racemic mixture of (2S,3R) and (2R,3S) products. The deoxygenation of epoxide to regenerate benzoyl-CoA by a diferrous form of the enzyme proceeds via a stepwise mechanism. The C2–O bond cleavage happens first, coupled with one electron transfer from one iron center to the substrate, to form a radical intermediate, which is followed by the second C3–O bond cleavage. The first step is rate-limiting with a barrier of only 10.8 kcal mol^–1. Further experimental studies are encouraged to verify our results.

Electronic structural flexibility of heterobimetallic Mn/Fe cofactors: R2lox and R2c proteins

The electronic structure of the Mn/Fe cofactor identified in a new class of oxidases (R2lox) described by Andersson and Högbom [Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 5633] is reported. The R2lox protein is homologous to the small subunit of class Ic ribonucleotide reductase (R2c) but has a completely different in vivo function. Using multifrequency EPR and related pulse techniques, it is shown that the cofactor of R2lox represents an antiferromagnetically coupled Mn(III)/Fe(III) dimer linked by a μ-hydroxo/bis-μ-carboxylato bridging network. The Mn(III) ion is coordinated by a single water ligand. The R2lox cofactor is photoactive, converting into a second form (R2loxPhoto) upon visible illumination at cryogenic temperatures (77 K) that completely decays upon warming. This second, unstable form of the cofactor more closely resembles the Mn(III)/Fe(III) cofactor seen in R2c. It is shown that the two forms of the R2lox cofactor differ primarily in terms of the local site geometry and electronic state of the Mn(III) ion, as best evidenced by a reorientation of its unique (55)Mn hyperfine axis. Analysis of the metal hyperfine tensors in combination with density functional theory (DFT) calculations suggests that this change is triggered by deprotonation of the μ-hydroxo bridge. These results have important consequences for the mixed-metal R2c cofactor and the divergent chemistry R2lox and R2c perform.

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.

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.

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.

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.

Rapid X-ray photoreduction of dimetal-oxygen cofactors in ribonucleotide reductase

Prototypic dinuclear metal cofactors with varying metallation constitute a class of O2-activating catalysts in numerous enzymes such as ribonucleotide reductase. Reliable structures are required to unravel the reaction mechanisms. However, protein crystallography data may be compromised by x-ray photoreduction (XRP). We studied XPR of Fe(III)Fe(III) and Mn(III)Fe(III) sites in the R2 subunit of Chlamydia trachomatis ribonucleotide reductase using x-ray absorption spectroscopy. Rapid and biphasic x-ray photoreduction kinetics at 20 and 80 K for both cofactor types suggested sequential formation of (III,II) and (II,II) species and similar redox potentials of iron and manganese sites. Comparing with typical x-ray doses in crystallography implies that (II,II) states are reached in <1 s in such studies. First-sphere metal coordination and metal-metal distances differed after chemical reduction at room temperature and after XPR at cryogenic temperatures, as corroborated by model structures from density functional theory calculations. The inter-metal distances in the XPR-induced (II,II) states, however, are similar to R2 crystal structures. Therefore, crystal data of initially oxidized R2-type proteins mostly contain photoreduced (II,II) cofactors, which deviate from the native structures functional in O2 activation, explaining observed variable metal ligation motifs. This situation may be remedied by novel femtosecond free electron-laser protein crystallography techniques.

Radical-translocation intermediates and hurdling of pathway defects in "super-oxidized" (Mn(IV)/Fe(IV)) Chlamydia trachomatis ribonucleotide reductase

A class I ribonucleotide reductase (RNR) uses either a tyrosyl radical (Y(•)) or a Mn(IV)/Fe(III) cluster in its β subunit to oxidize a cysteine residue ∼35 Å away in its α subunit, generating a thiyl radical that abstracts hydrogen (H(•)) from the substrate. With either oxidant, the inter-subunit "hole-transfer" or "radical-translocation" (RT) process is thought to occur by a "hopping" mechanism involving multiple tyrosyl (and perhaps one tryptophanyl) radical intermediates along a specific pathway. The hopping intermediates have never been directly detected in a Mn/Fe-dependent (class Ic) RNR nor in any wild-type (wt) RNR. The Mn(IV)/Fe(III) cofactor of Chlamydia trachomatis RNR assembles via a Mn(IV)/Fe(IV) intermediate. Here we show that this cofactor-assembly intermediate can propagate a hole into the RT pathway when α is present, accumulating radicals with EPR spectra characteristic of Y(•)'s. The dependence of Y(•) accumulation on the presence of substrate suggests that RT within this "super-oxidized" enzyme form is gated by the protein, and the failure of a β variant having the subunit-interfacial pathway Y substituted by phenylalanine to support radical accumulation implies that the Y(•)(s) in the wt enzyme reside(s) within the RT pathway. Remarkably, two variant β proteins having pathway substitutions rendering them inactive in their Mn(IV)/Fe(III) states can generate the pathway Y(•)'s in their Mn(IV)/Fe(IV) states and also effect nucleotide reduction. Thus, the use of the more oxidized cofactor permits the accumulation of hopping intermediates and the "hurdling" of engineered defects in the RT pathway.

Metallation and mismetallation of iron and manganese proteins in vitro and in vivo: the class I ribonucleotide reductases as a case study

How cells ensure correct metallation of a given protein and whether a degree of promiscuity in metal binding has evolved are largely unanswered questions. In a classic case, iron- and manganese-dependent superoxide dismutases (SODs) catalyze the disproportionation of superoxide using highly similar protein scaffolds and nearly identical active sites. However, most of these enzymes are active with only one metal, although both metals can bind in vitro and in vivo. Iron(ii) and manganese(ii) bind weakly to most proteins and possess similar coordination preferences. Their distinct redox properties suggest that they are unlikely to be interchangeable in biological systems except when they function in Lewis acid catalytic roles, yet recent work suggests this is not always the case. This review summarizes the diversity of ways in which iron and manganese are substituted in similar or identical protein frameworks. As models, we discuss (1) enzymes, such as epimerases, thought to use Fe(II) as a Lewis acid under normal growth conditions but which switch to Mn(II) under oxidative stress; (2) extradiol dioxygenases, which have been found to use both Fe(II) and Mn(II), the redox role of which in catalysis remains to be elucidated; (3) SODs, which use redox chemistry and are generally metal-specific; and (4) the class I ribonucleotide reductases (RNRs), which have evolved unique biosynthetic pathways to control metallation. The primary focus is the class Ib RNRs, which can catalyze formation of a stable radical on a tyrosine residue in their β2 subunits using either a di-iron or a recently characterized dimanganese cofactor. The physiological roles of enzymes that can switch between iron and manganese cofactors are discussed, as are insights obtained from the studies of many groups regarding iron and manganese homeostasis and the divergent and convergent strategies organisms use for control of protein metallation. We propose that, in many of the systems discussed, "discrimination" between metals is not performed by the protein itself, but it is instead determined by the environment in which the protein is expressed.

Evidence that the β subunit of Chlamydia trachomatis ribonucleotide reductase is active with the manganese ion of its manganese(IV)/iron(III) cofactor in site 1

The reaction of a class I ribonucleotide reductase (RNR) begins when a cofactor in the β subunit oxidizes a cysteine residue ~35 Å away in the α subunit, generating a thiyl radical. In the class Ic enzyme from Chlamydia trachomatis (Ct), the cysteine oxidant is the Mn(IV) ion of a Mn(IV)/Fe(III) cluster, which assembles in a reaction between O(2) and the Mn(II)/Fe(II) complex of β. The heterodinuclear nature of the cofactor raises the question of which site, 1 or 2, contains the Mn(IV) ion. Because site 1 is closer to the conserved location of the cysteine-oxidizing tyrosyl radical of class Ia and Ib RNRs, we suggested that the Mn(IV) ion most likely resides in this site (i.e., (1)Mn(IV)/(2)Fe(III)), but a subsequent computational study favored its occupation of site 2 ((1)Fe(III)/(2)Mn(IV)). In this work, we have sought to resolve the location of the Mn(IV) ion in Ct RNR-β by correlating X-ray crystallographic anomalous scattering intensities with catalytic activity for samples of the protein reconstituted in vitro by two different procedures. In samples containing primarily Mn(IV)/Fe(III) clusters, Mn preferentially occupies site 1, but some anomalous scattering from site 2 is observed, implying that both (1)Mn(II)/(2)Fe(II) and (1)Fe(II)/(2)Mn(II) complexes are competent to react with O(2) to produce the corresponding oxidized states. However, with diminished Mn(II) loading in the reconstitution, there is no evidence for Mn occupancy of site 2, and the greater activity of these "low-Mn" samples on a per-Mn basis implies that the (1)Mn(IV)/(2)Fe(III)-β is at least the more active of the two oxidized forms and may be the only active form.

High-valent [MnFe] and [FeFe] cofactors in ribonucleotide reductases

Ribonucleotide reductases (RNRs) are essential for DNA synthesis in most organisms. In class-Ic RNR from Chlamydia trachomatis (Ct), a MnFe cofactor in subunit R2 forms the site required for enzyme activity, instead of an FeFe cofactor plus a redox-active tyrosine in class-Ia RNRs, for example in mouse (Mus musculus, Mm). For R2 proteins from Ct and Mm, either grown in the presence of, or reconstituted with Mn and Fe ions, structural and electronic properties of higher valence MnFe and FeFe sites were determined by X-ray absorption spectroscopy and complementary techniques, in combination with bond-valence-sum and density functional theory calculations. At least ten different cofactor species could be tentatively distinguished. In Ct R2, two different Mn(IV)Fe(III) site configurations were assigned either L(4)Mn(IV)(μO)(2)Fe(III)L(4) (metal-metal distance of ~2.75Å, L = ligand) prevailing in metal-grown R2, or L(4)Mn(IV)(μO)(μOH)Fe(III)L(4) (~2.90Å) dominating in metal-reconstituted R2. Specific spectroscopic features were attributed to an Fe(IV)Fe(III) site (~2.55Å) with a L(4)Fe(IV)(μO)(2)Fe(III)L(3) core structure. Several Mn,Fe(III)Fe(III) (~2.9-3.1Å) and Mn,Fe(III)Fe(II) species (~3.3-3.4Å) likely showed 5-coordinated Mn(III) or Fe(III). Rapid X-ray photoreduction of iron and shorter metal-metal distances in the high-valent states suggested radiation-induced modifications in most crystal structures of R2. The actual configuration of the MnFe and FeFe cofactors seems to depend on assembly sequences, bound metal type, valence state, and previous catalytic activity involving subunit R1. In Ct R2, the protonation of a bridging oxide in the Mn(IV)(μO)(μOH)Fe(III) core may be important for preventing premature site reduction and initiation of the radical chemistry in R1.

Mössbauer properties of the diferric cluster and the differential iron(II)-binding affinity of the iron sites in protein R2 of class Ia Escherichia coli ribonucleotide reductase: a DFT/electrostatics study

The R2 subunit of class-Ia ribonucleotide reductase (RNR) from Escherichia coli (E. coli) contains a diiron active site. Starting from the apo-protein and Fe(II) in solution at low Fe(II)/apoR2 ratios, mononuclear Fe(II) binding is observed indicating possible different Fe(II) binding affinities for the two alternative sites. Further, based on their Mössbauer spectroscopy and two-iron-isotope reaction experiments, Bollinger et al. (J. Am. Chem. Soc., 1997, 119, 5976-5977) proposed that the site Fe1, which bonds to Asp84, should be associated with the higher observed (57)Fe Mössbauer quadrupole splitting (2.41 mm s(-1)) and lower isomer shift (0.45 mm s(-1)) in the Fe(III)Fe(III) state, site Fe2, which is further from Tyr122, should have a greater affinity for Fe(II) binding than site Fe1, and Fe(IV) in the intermediate X state should reside at site Fe2. In this paper, using 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) methodologies, we have demonstrated that the observed large quadrupole splitting for the diferric state R2 does come from site Fe1(III) and it is mainly caused by the binding position of the carboxylate group of the Asp84 sidechain. Further, a series of active site clusters with mononuclear Fe(II) binding at either site Fe1 or Fe2 have been studied, which show that with a single dielectric medium outside the active site quantum region, there is no energetic preference for Fe(II) binding at one site over another. However, when including the explicit extended protein environment in the PB-SCRF model, the reaction field favors the Fe(II) binding at site Fe2 rather than at site Fe1 by ~9 kcal mol(-1). Therefore our calculations support the proposal of the previous Mössbauer spectroscopy and two-iron-isotope reaction experiments by Bollinger et al.

Class I ribonucleotide reductases: metallocofactor assembly and repair in vitro and in vivo

Incorporation of metallocofactors essential for the activity of many enyzmes is a major mechanism of posttranslational modification. The cellular machinery required for these processes in the case of mono- and dinuclear nonheme iron and manganese cofactors has remained largely elusive. In addition, many metallocofactors can be converted to inactive forms, and pathways for their repair have recently come to light. The class I ribonucleotide reductases (RNRs) catalyze the conversion of nucleotides to deoxynucleotides and require dinuclear metal clusters for activity: an Fe(III)Fe(III)-tyrosyl radical (Y•) cofactor (class Ia), a Mn(III)Mn(III)-Y• cofactor (class Ib), and a Mn(IV)Fe(III) cofactor (class Ic). The class Ia, Ib, and Ic RNRs are structurally homologous and contain almost identical metal coordination sites. Recent progress in our understanding of the mechanisms by which the cofactor of each of these RNRs is generated in vitro and in vivo and by which the damaged cofactors are repaired is providing insight into how nature prevents mismetallation and orchestrates active cluster formation in high yields.

Cyanobacterial alkane biosynthesis further expands the catalytic repertoire of the ferritin-like 'di-iron-carboxylate' proteins

Enzymes that activate dioxygen at carboxylate-bridged non-heme diiron clusters residing within ferritin-like, four-helix-bundle protein architectures have crucial roles in, among other processes, the global carbon cycle (e.g. soluble methane monooxygenase), fatty acid biosynthesis [plant fatty acyl-acyl carrier protein (ACP) desaturases], DNA biosynthesis [the R2 or β2 subunits of class Ia ribonucleotide reductases (RNRs)], and cellular iron trafficking (ferritins). Classic studies on class Ia RNRs showed long ago how this obligatorily oxidative di-iron/O2 chemistry can be used to activate an enzyme for even a reduction reaction, and more recent investigations of class Ib and Ic RNRs, coupled with earlier studies on dimanganese catalases, have shown that members of this protein family can also incorporate either one or two Mn ions and use them in place of iron for redox catalysis. These two strategies--oxidative activation for non-oxidative reactions and use of alternative metal ions--expand the catalytic repertoire of the family, probably to include activities that remain to be discovered. Indeed, a recent study has suggested that fatty aldehyde decarbonylases (ADs) from cyanobacteria, purported to catalyze a redox-neutral cleavage of a Cn aldehyde to the Cn-1 alkane (or alkene) and CO, also belong to this enzyme family and are most similar in structure to two other members with heterodinuclear (Mn-Fe) cofactors. Here, we first briefly review both the chemical principles underlying the O2-dependent oxidative chemistry of the 'classical' di-iron-carboxylate proteins and the two aforementioned strategies that have expanded their functional range, and then consider what metal ion(s) and what chemical mechanism(s) might be employed by the newly discovered cyanobacterial ADs.

Quantum Chemical Modeling of Enzymatic Reactions: The Case of Decarboxylation

We present a systematic study of the decarboxylation step of the enzyme aspartate decarboxylase with the purpose of assessing the quantum chemical cluster approach for modeling this important class of decarboxylase enzymes. Active site models ranging in size from 27 to 220 atoms are designed, and the barrier and reaction energy of this step are evaluated. To model the enzyme surrounding, homogeneous polarizable medium techniques are used with several dielectric constants. The main conclusion is that when the active site model reaches a certain size, the solvation effects from the surroundings saturate. Similar results have previously been obtained from systematic studies of other classes of enzymes, suggesting that they are of a quite general nature.

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.

Two distinct mechanisms of inactivation of the class Ic ribonucleotide reductase from Chlamydia trachomatis by hydroxyurea: implications for the protein gating of intersubunit electron transfer

Catalysis by a class I ribonucleotide reductase (RNR) begins when a cysteine (C) residue in the alpha(2) subunit is oxidized to a thiyl radical (C(*)) by a cofactor approximately 35 A away in the beta(2) subunit. In a class Ia or Ib RNR, a stable tyrosyl radical (Y(*)) is the C oxidant, whereas a Mn(IV)/Fe(III) cluster serves this function in the class Ic enzyme from Chlamydia trachomatis (Ct). It is thought that, in either case, a chain of Y residues spanning the two subunits mediates C oxidation by forming transient "pathway" Y(*)s in a multistep electron transfer (ET) process that is "gated" by the protein so that it occurs only in the ready holoenzyme complex. The drug hydroxyurea (HU) inactivates both Ia/b and Ic beta(2) subunits by reducing their C oxidants. Reduction of the stable cofactor Y(*) (Y122(*)) in Escherichia coli class Ia beta(2) is faster in the presence of alpha(2) and a substrate (CDP), leading to speculation that HU might intercept a transient ET pathway Y(*) under these turnover conditions. Here we show that this mechanism is one of two that are operant in HU inactivation of the Ct enzyme. HU reacts with the Mn(IV)/Fe(III) cofactor to give two distinct products: the previously described homogeneous Mn(III)/Fe(III)-beta(2) complex, which forms only under turnover conditions (in the presence of alpha(2) and the substrate), and a distinct, diamagnetic Mn/Fe cluster, which forms approximately 900-fold less rapidly as a second phase in the reaction under turnover conditions and as the sole outcome in the reaction of Mn(IV)/Fe(III)-beta(2) only. Formation of Mn(III)/Fe(III)-beta(2) also requires (i) either Y338, the subunit-interfacial ET pathway residue of beta(2), or Y222, the surface residue that relays the "extra electron" to the Mn(IV)/Fe(IV) intermediate during activation of beta(2) but is not part of the catalytic ET pathway, and (ii) W51, the cofactor-proximal residue required for efficient ET between either Y222 or Y338 and the cofactor. The combined requirements for the catalytic subunit, the substrate, and, most importantly, a functional surface-to-cofactor electron relay system imply that HU effects the Mn(IV)/Fe(III) --> Mn(III)/Fe(III) reduction by intercepting a Y(*) that forms when the ready holoenzyme complex is assembled, the ET gate is opened, and the Mn(IV) oxidizes either Y222 or Y338.