Spectroscopic and computational studies of the de novo designed protein DF2t: correlation to the biferrous active site of ribonucleotide reductase and factors that affect O2 reactivity

Allosteric cooperation in a de novo-designed two-domain protein

We describe the de novo design of an allosterically regulated protein, which comprises two tightly coupled domains. One domain is based on the DF (Due Ferri in Italian or two-iron in English) family of de novo proteins, which have a diiron cofactor that catalyzes a phenol oxidase reaction, while the second domain is based on PS1 (Porphyrin-binding Sequence), which binds a synthetic Zn-porphyrin (ZnP). The binding of ZnP to the original PS1 protein induces changes in structure and dynamics, which we expected to influence the catalytic rate of a fused DF domain when appropriately coupled. Both DF and PS1 are four-helix bundles, but they have distinct bundle architectures. To achieve tight coupling between the domains, they were connected by four helical linkers using a computational method to discover the most designable connections capable of spanning the two architectures. The resulting protein, DFP1 (Due Ferri Porphyrin), bound the two cofactors in the expected manner. The crystal structure of fully reconstituted DFP1 was also in excellent agreement with the design, and it showed the ZnP cofactor bound over 12 Å from the dimetal center. Next, a substrate-binding cleft leading to the diiron center was introduced into DFP1. The resulting protein acts as an allosterically modulated phenol oxidase. Its Michaelis-Menten parameters were strongly affected by the binding of ZnP, resulting in a fourfold tighter Km and a 7-fold decrease in kcat These studies establish the feasibility of designing allosterically regulated catalytic proteins, entirely from scratch.

De Novo Design of Four-Helix Bundle Metalloproteins: One Scaffold, Diverse Reactivities

De novo protein design represents an attractive approach for testing and extending our understanding of metalloprotein structure and function. Here, we describe our work on the design of DF (Due Ferri or two-iron in Italian), a minimalist model for the active sites of much larger and more complex natural diiron and dimanganese proteins. In nature, diiron and dimanganese proteins protypically bind their ions in 4-Glu, 2-His environments, and they catalyze diverse reactions, ranging from hydrolysis, to O₂-dependent chemistry, to decarbonylation of aldehydes. In the design of DF, the position of each atom-including the backbone, the first-shell ligands, the second-shell hydrogen-bonded groups, and the well-packed hydrophobic core-was bespoke using precise mathematical equations and chemical principles. The first member of the DF family was designed to be of minimal size and complexity and yet to display the quintessential elements required for binding the dimetal cofactor. After thoroughly characterizing its structural, dynamic, spectroscopic, and functional properties, we added additional complexity in a rational stepwise manner to achieve increasingly sophisticated catalytic functions, ultimately demonstrating substrate-gated four-electron reduction of O₂ to water. We also briefly describe the extension of these studies to the design of proteins that bind nonbiological metal cofactors (a synthetic porphyrin and a tetranuclear cluster), and a Zn²⁺/proton antiporting membrane protein. Together these studies demonstrate a successful and generally applicable strategy for de novo metalloprotein design, which might indeed mimic the process by which primordial metalloproteins evolved. We began the design process with a highly symmetrical backbone and binding site, by using point-group symmetry to assemble the secondary structures that position the amino acid side chains required for binding. The resulting models provided a rough starting point and initial parameters for the subsequent precise design of the final protein using modern methods of computational protein design. Unless the desired site is itself symmetrical, this process requires reduction of the symmetry or lifting it altogether. Nevertheless, the initial symmetrical structure can be helpful to restrain the search space during assembly of the backbone. Finally, the methods described here should be generally applicable to the design of highly stable and robust catalysts and sensors. There is considerable potential in combining the efficiency and knowledge base associated with homogeneous metal catalysis with the programmability, biocompatibility, and versatility of proteins. While the work reported here focuses on testing and learning the principles of natural metalloproteins by designing and studying proteins one at a time, there is also considerable potential for using designed proteins that incorporate both biological and nonbiological metal ion cofactors for the evolution of novel catalysts.

Design and engineering of artificial oxygen-activating metalloenzymes

Many efforts are being made in the design and engineering of metalloenzymes with catalytic properties fulfilling the needs of practical applications. Progress in this field has recently been accelerated by advances in computational, molecular and structural biology. This review article focuses on the recent examples of oxygen-activating metalloenzymes, developed through the strategies of de novo design, miniaturization processes and protein redesign. Considerable progress in these diverse design approaches has produced many metal-containing biocatalysts able to adopt the functions of native enzymes or even novel functions beyond those found in Nature.

Structure/function correlations over binuclear non-heme iron active sites

Binuclear non-heme iron enzymes activate O2 to perform diverse chemistries. Three different structural mechanisms of O2 binding to a coupled binuclear iron site have been identified utilizing variable-temperature, variable-field magnetic circular dichroism spectroscopy (VTVH MCD). For the μ-OH-bridged Fe(II)2 site in hemerythrin, O2 binds terminally to a five-coordinate Fe(II) center as hydroperoxide with the proton deriving from the μ-OH bridge and the second electron transferring through the resulting μ-oxo superexchange pathway from the second coordinatively saturated Fe(II) center in a proton-coupled electron transfer process. For carboxylate-only-bridged Fe(II)2 sites, O2 binding as a bridged peroxide requires both Fe(II) centers to be coordinatively unsaturated and has good frontier orbital overlap with the two orthogonal O2 π* orbitals to form peroxo-bridged Fe(III)2 intermediates. Alternatively, carboxylate-only-bridged Fe(II)2 sites with only a single open coordination position on an Fe(II) enable the one-electron formation of Fe(III)-O2 (-) or Fe(III)-NO(-) species. Finally, for the peroxo-bridged Fe(III)2 intermediates, further activation is necessary for their reactivities in one-electron reduction and electrophilic aromatic substitution, and a strategy consistent with existing spectral data is discussed.

Design of dinuclear manganese cofactors for bacterial reaction centers

A compelling target for the design of electron transfer proteins with novel cofactors is to create a model for the oxygen-evolving complex, a Mn4Ca cluster, of photosystem II. A mononuclear Mn cofactor can be added to the bacterial reaction center, but the addition of multiple metal centers is constrained by the native protein architecture. Alternatively, metal centers can be incorporated into artificial proteins. Designs for the addition of dinuclear metal centers to four-helix bundles resulted in three artificial proteins with ligands for one, two, or three dinuclear metal centers able to bind Mn. The three-dimensional structure determined by X-ray crystallography of one of the Mn-proteins confirmed the design features and revealed details concerning coordination of the Mn center. Electron transfer between these artificial Mn-proteins and bacterial reaction centers was investigated using optical spectroscopy. After formation of a light-induced, charge-separated state, the experiments showed that the Mn-proteins can donate an electron to the oxidized bacteriochlorophyll dimer of modified reaction centers, with the Mn-proteins having additional metal centers being more effective at this electron transfer reaction. Modeling of the structure of the Mn-protein docked to the reaction center showed that the artificial protein likely binds on the periplasmic surface similarly to cytochrome c2, the natural secondary donor. Combining reaction centers with exogenous artificial proteins provides the opportunity to create ligands and investigate the influence of inhomogeneous protein environments on multinuclear redox-active metal centers. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Systematic Perturbations of Binuclear Non-heme Iron Sites: Structure and Dioxygen Reactivity of de Novo Due Ferri Proteins

DFsc (single-chain due ferri) proteins allow for modeling binuclear non-heme iron enzymes with a similar fold. Three 4A → 4G variants of DFsc were studied to investigate the effects of (1) increasing the size of the substrate/solvent access channel (G4DFsc), (2) including an additional His residue in the first coordination sphere along with three additional helix-stabilizing mutations [3His-G4DFsc(Mut3)], and (3) the three helix-stabilizing mutations alone [G4DFsc(Mut3)] on the biferrous structures and their O2 reactivities. Near-infrared circular dichroism and magnetic circular dichroism (MCD) spectroscopy show that the 4A → 4G mutations increase coordination of the diiron site from 4-coordinate/5-coordinate to 5-coordinate/5-coordinate, likely reflecting increased solvent accessibility. While the three helix-stabilizing mutations [G4DFsc(Mut3)] do not affect the coordination number, addition of the third active site His residue [3His-G4DFsc(Mut3)] results in a 5-coordinate/6-coordinate site. Although all 4A→ 4G variants have significantly slower pseudo-first-order rates when reacting with excess O2 than DFsc (∼2 s(-1)), G4DFsc and 3His-G4DFsc(Mut3) have rates (∼0.02 and ∼0.04 s(-1)) faster than that of G4DFsc(Mut3) (∼0.002 s(-1)). These trends in the rate of O2 reactivity correlate with exchange coupling between the Fe(II) sites and suggest that the two-electron reduction of O2 occurs through end-on binding at one Fe(II) rather than through a peroxy-bridged intermediate. UV-vis absorption and MCD spectroscopies indicate that an Fe(III)Fe(III)-OH species first forms in all three variants but converts into an Fe(III)-μ-OH-Fe(III) species only in the 2-His forms, a process inhibited by the additional active site His ligand that coordinatively saturates one of the iron centers in 3His-G4DFsc(Mut3).

Spectroscopic studies of single and double variants of M ferritin: lack of conversion of a biferrous substrate site into a cofactor site for O2 activation

Ferritin has a binuclear non-heme iron active site that functions to oxidize iron as a substrate for formation of an iron mineral core. Other enzymes of this class have tightly bound diiron cofactor sites that activate O₂ to react with substrate. Ferritin has an active site ligand set with 1-His/4-carboxylate/1-Gln rather than the 2-His/4-carboxylate set of the cofactor site. This ligand variation has been thought to make a major contribution to this biferrous substrate rather than cofactor site reactivity. However, the Q137E/D140H double variant of M ferritin, has a ligand set that is equivalent to most of the diiron cofactor sites, yet did not rapidly react with O₂ or generate the peroxy intermediate observed in the cofactor sites. Therefore, in this study, a combined spectroscopic methodology of circular dichroism (CD)/magnetic CD (MCD)/variable temperature, variable field (VTVH) MCD has been applied to evaluate the factors required for the rapid O₂ activation observed in cofactor sites. This methodology defines the coordination environment of each iron and the bridging ligation of the biferrous active sites in the double and corresponding single variants of frog M ferritin. Based on spectral changes, the D140H single variant has the new His ligand binding, and the Q137E variant has the new carboxylate forming a μ-1,3 bridge. The spectra for the Q137E/D140H double variant, which has the cofactor ligand set, however, reflects a site that is more coordinately saturated than the cofactor sites in other enzymes including ribonucleotide reductase, indicating the presence of additional water ligation. Correlation of this double variant and the cofactor sites to their O₂ reactivities indicates that electrostatic and steric changes in the active site and, in particular, the hydrophobic nature of a cofactor site associated with its second sphere protein environment, make important contributions to the activation of O₂ by the binuclear non-heme iron enzymes.

Designing functional metalloproteins: from structural to catalytic metal sites

Metalloenzymes efficiently catalyze some of the most important and difficult reactions in nature. For many years, coordination chemists have effectively used small molecule models to understand these systems. More recently, protein design has been shown to be an effective approach for mimicking metal coordination environments. Since the first designed proteins were reported, much success has been seen for incorporating metal sites into proteins and attaining the desired coordination environment but until recently, this has been with a lack of significant catalytic activity. Now there are examples of designed metalloproteins that, although not yet reaching the activity of native enzymes, are considerably closer. In this review, we highlight work leading up to the design of a small metalloprotein containing two metal sites, one for structural stability (HgS3) and the other a separate catalytic zinc site to mimic carbonic anhydrase activity (ZnN3O). The first section will describe previous studies that allowed for a high affinity thiolate site that binds heavy metals in a way that stabilizes three-stranded coiled coils. The second section will examine ways of preparing histidine rich environments that lead to metal based hydrolytic catalysts. We will also discuss other recent examples of the design of structural metal sites and functional metalloenzymes. Our work demonstrates that attaining the proper first coordination geometry of a metal site can lead to a significant fraction of catalytic activity, apparently independent of the type of secondary structure of the surrounding protein environment. We are now in a position to begin to meet the challenge of building a metalloenzyme systematically from the bottom-up by engineering and analyzing interactions directly around the metal site and beyond.

Intermediate P* from soluble methane monooxygenase contains a diferrous cluster

During a single turnover of the hydroxylase component (MMOH) of soluble methane monooxygenase from Methylosinus trichosporium OB3b, several discrete intermediates are formed. The diiron cluster of MMOH is first reduced to the Fe(II)Fe(II) state (H(red)). O₂ binds rapidly at a site away from the cluster to form the Fe(II)Fe(II) intermediate O, which converts to an Fe(III)Fe(III)-peroxo intermediate P and finally to the Fe(IV)Fe(IV) intermediate Q. Q binds and reacts with methane to yield methanol and water. The rate constants for these steps are increased by a regulatory protein, MMOB. Previously reported transient kinetic studies have suggested that an intermediate P* forms between O and P in which the g = 16 EPR signal characteristic of the reduced diiron cluster of H(red) and O is lost. This was interpreted as signaling oxidation of the cluster, but a low level of accumulation of P* prevented further characterization. In this study, three methods for directly detecting and trapping P* are applied together to allow its spectroscopic and kinetic characterization. First, the MMOB mutant His33Ala is used to specifically slow the decay of P* without affecting its formation rate, leading to its nearly quantitative accumulation. Second, spectra-kinetic data collection is used to provide a sensitive measure of the formation and decay rate constants of intermediates as well as their optical spectra. Finally, the substrate furan is included to react with Q and quench its strong chromophore. The optical spectrum of P* closely mimics those of H(red) and O, but it is distinctly different from that of P. The reaction cycle rate constants allowed prediction of the times for maximal accumulation of the intermediates. Mössbauer spectra of rapid freeze-quench samples at these times show that the intermediates are formed at almost exactly the predicted levels. The Mössbauer spectra show that the diiron cluster of P*, quite unexpectedly, is in the Fe(II)Fe(II) state. Thus, the loss of the g = 16 EPR signal results from a change in the electronic structure of the Fe(II)Fe(II) center rather than oxidation. The similarity of the optical and Mössbauer spectra of H(red), O, and P* suggests that only subtle changes occur in the electronic and physical structure of the diiron cluster as P* forms. Nevertheless, the changes that do occur are necessary for O₂ to be activated for hydrocarbon oxidation.

Structure and dynamics of a dizinc metalloprotein: effect of charge transfer and polarization

Structures and dynamics of a recently designed dizinc metalloprotein (DFsc) (J. Mol. Biol. 2003, 334, 1101) are studied by molecular dynamics simulation using a dynamically adapted polarized force field derived from fragment quantum calculation for protein in solvent. To properly describe the effect of charge transfer and polarization in the present approach, quantum chemistry calculation of the zinc-binding group is periodically performed (on-the-fly) to update the atomic charges of the zinc-binding group during the MD simulation. Comparison of the present result with those obtained from simulations under standard AMBER force field reveals that charge transfer and polarization are critical to maintaining the correct asymmetric metal coordination in the DFsc. Detailed analysis of the result also shows that dynamic fluctuation of the zinc-binding group facilitates solvent interaction with the zinc ions. In particular, the dynamic fluctuation of the zinc-zinc distance is shown to be an important feature of the catalytic function of the di-ion zinc-binding group. Our study demonstrates that the dynamically adapted polarization approach is computationally practical and can be used to study other metalloprotein systems.

Computational protein design: engineering molecular diversity, nonnatural enzymes, nonbiological cofactor complexes, and membrane proteins

Computational and theoretical methods are advancing protein design as a means to create and investigate proteins. Such efforts further our capacity to control, design and understand biomolecular structure, sequence and function. Herein, the focus is on some recent applications that involve using theoretical and computational methods to guide the design of protein sequence ensembles, new enzymes, proteins with novel cofactors, and membrane proteins.

CD and MCD spectroscopic studies of the two Dps miniferritin proteins from Bacillus anthracis: role of O2 and H2O2 substrates in reactivity of the diiron catalytic centers

DNA protection during starvation (Dps) proteins are miniferritins found in bacteria and archaea that provide protection from uncontrolled Fe(II)/O radical chemistry; thus the catalytic sites are targets for antibiotics against pathogens, such as anthrax. Ferritin protein cages synthesize ferric oxymineral from Fe(II) and O(2)/H(2)O(2), which accumulates in the large central cavity; for Dps, H(2)O(2) is the more common Fe(II) oxidant contrasting with eukaryotic maxiferritins that often prefer dioxygen. To better understand the differences in the catalytic sites of maxi- versus miniferritins, we used a combination of NIR circular dichroism (CD), magnetic circular dichroism (MCD), and variable-temperature, variable-field MCD (VTVH MCD) to study Fe(II) binding to the catalytic sites of the two Bacillus anthracis miniferritins: one in which two Fe(II) react with O(2) exclusively (Dps1) and a second in which both O(2) or H(2)O(2) can react with two Fe(II) (Dps2). Both result in the formation of iron oxybiomineral. The data show a single 5- or 6-coordinate Fe(II) in the absence of oxidant; Fe(II) binding to Dps2 is 30× more stable than Dps1; and the lower limit of K(D) for binding a second Fe(II), in the absence of oxidant, is 2-3 orders of magnitude weaker than for the binding of the single Fe(II). The data fit an equilibrium model where binding of oxidant facilitates formation of the catalytic site, in sharp contrast to eukaryotic M-ferritins where the binuclear Fe(II) centers are preformed before binding of O(2). The two different binding sequences illustrate the mechanistic range possible for catalytic sites of the family of ferritins.

Computational protein design: Advances in the design and redesign of biomolecular nanostructures

Computational protein design facilitates the continued development of methods for the design of biomolecular structure, sequence and function. Recent applications include the design of novel protein sequences and structures, proteins incorporating nonbiological components, protein assemblies, soluble variants of membrane proteins, and proteins that modulate membrane function.

Spectroscopic and metal-binding properties of DF3: an artificial protein able to accommodate different metal ions

The design, synthesis, and metal-binding properties of DF3, a new de novo designed di-iron protein model are described ("DF" represents due ferri, Italian for "two iron," "di-iron"). DF3 is the latest member of the DF family of synthetic proteins. They consist of helix-loop-helix hairpins, designed to dimerize and form an antiparallel four-helix bundle that encompasses a metal-binding site similar to those of non-heme carboxylate-bridged di-iron proteins. Unlike previous DF proteins, DF3 is highly soluble in water (up to 3 mM) and forms stable complexes with several metal ions (Zn, Co, and Mn), with the desired secondary structure and the expected stoichiometry of two ions per protein. UV-vis studies of Co(II) and Fe(III) complexes confirm a metal-binding environment similar to previous di-Co(II)- and di-Fe(III)-DF proteins, including the presence of a mu-oxo-di-Fe(III) unit. Interestingly, UV-vis, EPR, and resonance Raman studies suggest the interaction of a tyrosine adjacent to the di-Fe(III) center. The design of DF3 was aimed at increasing the accessibility of small molecules to the active site of the four-helix bundle. Indeed, binding of azide to the di-Fe(III) site demonstrates a more accessible metal site compared with previous DFs. In fact, fitting of the binding curve to the Hill equation allows us to quantify a 150% accessibility enhancement, with respect to DF2. All these results represent a significant step towards the development of a functional synthetic DF metalloprotein.

Spectroscopic definition of the biferrous and biferric sites in de novo designed four-helix bundle DFsc peptides: implications for O2 reactivity of binuclear non-heme iron enzymes

DFsc is a single chain de novo designed four-helix bundle peptide that mimics the core protein fold and primary ligand set of various binuclear non-heme iron enzymes. DFsc and the E11D, Y51L, and Y18F single amino acid variants have been studied using a combination of near-IR circular dichroism (CD), magnetic circular dichroism (MCD), variable temperature variable field MCD (VTVH MCD), and X-ray absorption (XAS) spectroscopies. The biferrous sites are all weakly antiferromagnetically coupled with mu-1,3 carboxylate bridges and one 4-coordinate and one 5-coordinate Fe, very similar to the active site of class I ribonucleotide reductase (R2) providing open coordination positions on both irons for dioxygen to bridge. From perturbations of the MCD and VTVH MCD the iron proximal to Y51 can be assigned as the 4-coordinate center, and XAS results show that Y51 is not bound to this iron in the reduced state. The two open coordination positions on one iron in the biferrous state would become occupied by dioxygen and Y51 along the O(2) reaction coordinate. Subsequent binding of Y51 functions as an internal spectral probe of the O(2) reaction and as a proton source that would promote loss of H(2)O(2). Coordination by a ligand that functions as a proton source could be a structural mechanism used by natural binuclear iron enzymes to drive their reactions past peroxo biferric level intermediates.

CD and MCD studies of the effects of component B variant binding on the biferrous active site of methane monooxygenase

The multicomponent soluble form of methane monooxygenase (sMMO) catalyzes the oxidation of methane through the activation of O 2 at a nonheme biferrous center in the hydroxylase component, MMOH. Reactivity is limited without binding of the sMMO effector protein, MMOB. Past studies show that mutations of specific MMOB surface residues cause large changes in the rates of individual steps in the MMOH reaction cycle. To define the structural and mechanistic bases for these observations, CD, MCD, and VTVH MCD spectroscopies coupled with ligand-field (LF) calculations are used to elucidate changes occurring near and at the MMOH biferrous cluster upon binding of MMOB and the MMOB variants. Perturbations to both the CD and MCD are observed upon binding wild-type MMOB and the MMOB variant that similarly increases O 2 reactivity. MMOB variants that do not greatly increase O 2 reactivity fail to cause one or both of these changes. LF calculations indicate that reorientation of the terminal glutamate on Fe2 reproduces the spectral perturbations in MCD. Although this structural change allows O 2 to bridge the diiron site and shifts the redox active orbitals for good overlap, it is not sufficient for enhanced O 2 reactivity of the enzyme. Binding of the T111Y-MMOB variant to MMOH induces the MCD, but not CD changes, and causes only a small increase in reactivity. Thus, both the geometric rearrangement at Fe2 (observed in MCD) coupled with a more global conformational change that may control O 2 access (probed by CD), induced by MMOB binding, are critical factors in the reactivity of sMMO.

Spectroscopic definition of the ferroxidase site in M ferritin: comparison of binuclear substrate vs cofactor active sites

Maxi ferritins, 24 subunit protein nanocages, are essential in humans, plants, bacteria, and other animals for the concentration and storage of iron as hydrated ferric oxide, while minimizing free radical generation or use by pathogens. Formation of the precursors to these ferric oxides is catalyzed at a nonheme biferrous substrate site, which has some parallels with the cofactor sites in other biferrous enzymes. A combination of circular dichroism (CD), magnetic circular dichroism (MCD), and variable-temperature, variable-field MCD (VTVH MCD) has been used to probe Fe(II) binding to the substrate active site in frog M ferritin. These data determined that the active site within each subunit consists of two inequivalent five-coordinate (5C) ferrous centers that are weakly antiferromagnetically coupled, consistent with a mu-1,3 carboxylate bridge. The active site ligand set is unusual and likely includes a terminal water bound to each Fe(II) center. The Fe(II) ions bind to the active sites in a concerted manner, and cooperativity among the sites in each subunit is observed, potentially providing a mechanism for the control of ferritin iron loading. Differences in geometric and electronic structure--including a weak ligand field, availability of two water ligands at the biferrous substrate site, and the single carboxylate bridge in ferritin--coincide with the divergent reaction pathways observed between this substrate site and the previously studied cofactor active sites.

Oxygen reactivity of the biferrous site in the de novo designed four helix bundle peptide DFsc: nature of the "intermediate" and reaction mechanism

The DFsc and DFscE11D de novo designed protein scaffolds support biomimetic diiron cofactor sites that react with dioxygen forming a 520 nm "intermediate" species with an apparent pseudo-first-order formation rate constant of 2.2 and 4.8 s-1, respectively. Resonance Raman spectroscopy shows that this absorption feature is due to a phenolate-to-ferric charge transfer transition arising from a single tyrosine residue coordinating terminally to one of the ferric ions in the site. Phenol coordination could provide a proton to promote rapid loss of a putative peroxo species.

Geometric and electronic structure studies of the binuclear nonheme ferrous active site of toluene-4-monooxygenase: parallels with methane monooxygenase and insight into the role of the effector proteins in O2 activation

Multicomponent monooxygenases, which carry out a variety of highly specific hydroxylation reactions, are of great interest as potential biocatalysts in a number of applications. These proteins share many similarities in structure and show a marked increase in O2 reactivity upon addition of an effector component. In this study, circular dichroism (CD), magnetic circular dichroism (MCD), and variable-temperature, variable-field (VTVH) MCD have been used to gain spectroscopic insight into the Fe(II)Fe(II) active site in the hydroxylase component of Toluene-4 monoxygenase (T4moH) and the complex of T4moH bound by its effector protein, T4moD. These results have been correlated to spectroscopic data and density functional theory (DFT) calculations on MmoH and its interaction with MmoB. Together, these data provide further insight into the geometric and electronic structure of these biferrous active sites and, in particular, the perturbation associated with component B/D binding. It is found that binding of the effector protein changes the geometry of one iron center and orientation of its redox active orbital to accommodate the binding of O2 in a bridged structure for efficient 2-electron transfer that can form a peroxo intermediate.

Solution NMR structure of a designed metalloprotein and complementary molecular dynamics refinement

We report the solution NMR structure of a designed dimetal-binding protein, di-Zn(II) DFsc, along with a secondary refinement step employing molecular dynamics techniques. Calculation of the initial NMR structural ensemble by standard methods led to distortions in the metal-ligand geometries at the active site. Unrestrained molecular dynamics using a nonbonded force field for the metal shell, followed by quantum mechanical/molecular mechanical dynamics of DFsc, were used to relax local frustrations at the dimetal site that were apparent in the initial NMR structure and provide a more realistic description of the structure. The MD model is consistent with NMR restraints, and in good agreement with the structural and functional properties expected for DF proteins. This work demonstrates that NMR structures of metalloproteins can be further refined using classical and first-principles molecular dynamics methods in the presence of explicit solvent to provide otherwise unavailable insight into the geometry of the metal center.

High-spin versus broken symmetry—Effect of DFT spin density representation on the geometries of three diiron (III) model compounds

Unrestricted density functional theory calculations have been conducted on three diiron(III) synthetic model compounds containing antiferromagnetically coupled high-spin (HS) irons for which crystallographic structures and Raman spectral data are available. Three density functionals have been employed: BPW91, PWC, and BOP. The study compares the effects on optimized geometries and harmonic vibrational frequencies of spin-paired (SP) low-spin, HS, and broken symmetry antiferromagnetically coupled singlet representations of the spin density distribution. The geometries around the diiron centers in the HS and broken symmetry (BS) representations are found to be similar, both markedly different from those arising from the SP representation. Small differences between the HS and BS results are seen in bond lengths, angles, Raman frequencies, and spin densities associated with oxo and peroxo bridges between the irons.

Computational protein design: structure, function and combinatorial diversity

Computational protein design has blossomed with the development of methods for addressing the complexities involved in specifying the structure, sequence and function of proteins. Recent applications include the design of novel functional membrane and soluble proteins, proteins incorporating non-biological components and protein combinatorial libraries.

Spin exchange effects on the physicochemical properties of tetraoxolene-bridged bimetallic complexes

The synthesis, physical, and spectroscopic properties of a series of metal complexes bridged by the redox-active chloranilate ligand are described. Compounds containing the (CAcat,cat)4- ligand, where (CAcat,cat)4- represents the fully reduced aromatic form of chloranilate, have been prepared by two different routes from H2CA and H4CA starting materials; the corresponding (CAsq,cat)3- analogue was obtained by one-electron oxidation with decamethylferrocenium tetrafluoroborate. Homo- and heterobimetallic complexes containing CrIII and GaIII with chloranilate have been prepared, yielding the following six complexes: [Ga2(tren)2(CAcat,cat)](BPh4)2 (1), [Ga2(tren)2(CAsq,cat)](BPh4)2(BF4) (2), [GaCr(tren)2(CAcat,cat)](BPh4)2 (3), [GaCr(tren)2(CAsq,cat)](BPh4)2(BF4) (4), [Cr2(tren)2(CAcat,cat)] (BPh4)2 (5), and [Cr2(tren)2(CAsq,cat)](BPh4)2(BF4) (6) (where tren is tris(2-aminoethyl)amine). Single-crystal X-ray structures have been obtained for complexes 1, 3, and 5; nearly identical C-C bond distances within the quinoidal ligand confirm the aromatic character of the bridge in each case. Complex 2 exhibits a temperature-independent magnetic moment of microeff = 1.64 +/- 0.04 microB in the solid state between 4 and 350 K, consistent with the CAsq,cat formulation of the ligand and an S = 1/2 ground state for complex 2. Complex 3 exhibits a value of microeff = 3.44 +/- 0.09 microB that is also temperature-independent, indicating an S = 3/2 ground state. Complexes 4-6 are all influenced by Heisenberg spin exchange. The temperature-independent behavior of complexes 4 and 6 indicate the presence of strong antiferromagnetic exchange between the CrIII and the (sq,cat) bridging radical yielding well-isolated ground states of S = 1 and 5/2 for 4 and 6, respectively. In contrast, complex 5 exhibits a weak intramolecular antiferromagnetic exchange interaction between the two CrIII centers (J = -2 cm-1 for H = -2Jŝ1.ŝ2) via superexchange through the diamagnetic CAcat,cat bridge. The absorption spectra of the CAsq,cat-containing complexes exhibit a number of sharp, relatively intense features in fluid solution. Group theoretical arguments coupled with a qualitative ligand-field analysis including the effects of Heisenberg spin exchange suggest that several of the observed transitions are a consequence of exchange interactions in both the ground- and excited-state manifolds of the compounds.

Computational study of iron(II) and -(III) complexes with a simple model human H ferritin ferroxidase center

Interaction of iron ions with a six-amino acid model of the ferroxidase center of human H chain ferritin has been examined in density functional theory calculations. The model, based on experimental studies of oxidation of Fe2+ at the center, consists of Glu27, Glu62, His65, Glu107, Gln141, and Ala144. Reasonable structures are obtained in a survey of types of iron complexes inferred to occur in the ferroxidase reaction. Structures of complexes of the model center with one and two Fe2+ ions, with diiron(III) bridged by peroxide and bridged by oxide-peroxide combinations, have been optimized. Calculations on diiron(III) complexes confirm that stable peroxide-bridged complexes can form and that the Fe-Fe distance in at least one is consistent with the experimental Fe-Fe distance observed in the blue peroxodiferric complex of ferritin.

Deoxyhypusine hydroxylase is a Fe(II)-dependent, HEAT-repeat enzyme. Identification of amino acid residues critical for Fe(II) binding and catalysis [corrected]

Deoxyhypusine hydroxylase (DOHH) catalyzes the final step in the post-translational synthesis of hypusine (N(epsilon)-(4-amino-2-hydroxybutyl)lysine) in eIF5A. DOHH is a HEAT-repeat protein with eight tandem helical hairpins in a symmetrical dyad. It contains two potential iron coordination sites (one on each dyad) composed of two strictly conserved His-Glu motifs. The purified human recombinant DOHH was a mixture of active holoenzyme containing 2 mol of iron/mol of DOHH and inactive metal-free apoenzyme. The two species could be distinguished by their different mobilities upon native gel electrophoresis. The DOHH apoenzyme exhibited markedly reduced levels of iron and activity. DOHH activity could be restored only by the addition of Fe2+ to the apoenzyme but not by other metals including Cd2+,Co2+,Cr2+,Cu2+,Mg2+,Mn2+,Ni2+, and Zn2+. The role of the strictly conserved His-Glu residues was evaluated by site-directed mutagenesis. Substitution of any single amino acid in the four His-Glu motifs with alanine abolished the enzyme activity. Of these eight alanine substitutions, six, including H56A, H89A, E90A, H207A, H240A, and E241A, caused a severe reduction in the iron content. Our results provide strong evidence that Fe(II) is the active-site-bound metal critical for DOHH catalysis and that the strictly conserved His-Glu motifs are essential for iron binding and catalysis. Furthermore, the iron to DOHH stoichiometry and dependence of iron binding on each of the four conserved His-Glu motifs suggest a binuclear iron mediated reaction mechanism, distinct from that of other Fe(II)-dependent protein hydroxylases, such as prolyl 4-hydroxylase or lysyl hydroxylases.

Avoiding high-valent iron intermediates: superoxide reductase and rubrerythrin

The Fenton or Fenton-type reaction between aqueous ferrous ion and hydrogen peroxide generates a highly oxidizing species, most often formulated as hydroxyl radical or ferryl ([Fe(IV)O](2+)). Intracellular Fenton-type chemistry can be lethal if not controlled. Nature has, therefore, evolved enzymes to scavenge superoxide and hydrogen peroxide, the reduced dioxygen species that initiate intracellular Fenton-type chemistry. Two such enzymes found predominantly in air-sensitive bacteria and archaea, superoxide reductase (SOR) and rubrerythrin (Rbr), functioning as a peroxidase (hydrogen peroxide reductase), contain non-heme iron. The iron coordination spheres in these enzymes contain five or six protein ligands from His and Glu residues, and, in the case of SOR, a Cys residue. SOR contains a mononuclear active site that is designed to protonate and rapidly expel peroxide generated as a product of the enzymatic reaction. The ferrous SOR reacts adventitiously but relatively slowly (several seconds to a few minutes) with exogenous hydrogen peroxide, presumably in a Fenton-type reaction. The diferrous active site of Rbr reacts more rapidly with hydrogen peroxide but can divert Fenton-type reactions towards the two-electron reduction of hydrogen peroxide to water. Proximal aromatic residues may function as radical sinks for Fenton-generated oxidants. Fenton-initiated damage to these iron active sites may become apparent only under extremely oxidizing intracellular conditions.