Construction of a catalytically active iron superoxide dismutase by rational protein design

Metal-binding studies for a de novo designed calcium-binding protein

To understand the key determinants in calcium-binding affinity, a calcium-binding site with pentagonal bipyramid geometry was designed into a non-calcium-binding protein, domain 1 of CD2. This metal-binding protein has five mutations with a net charge in the coordination sphere of -5 and is termed DEEEE. Fluorescence resonance energy transfer was used to determine the metal-binding affinity of DEEEE to the calcium analog terbium. The addition of protein concentration to Tb(III) solution results in a large enhancement of Tb(III) fluorescence due to energy transfer between terbium ions and aromatic residues in CD2-D1. In addition, both calcium and lanthanum compete with terbium for the same desired metal binding pocket. Our designed protein exhibits a stronger affinity for Tb(III), with a K(d) of 21 microM, than natural calcium-binding proteins with a similar Greek key scaffold.

Enzyme-like proteins by computational design

We report the development and initial experimental validation of a computational design procedure aimed at generating enzyme-like protein catalysts called "protozymes." Our design approach utilizes a "compute and build" strategy that is based on the physical/chemical principles governing protein stability and catalytic mechanism. By using the catalytically inert 108-residue Escherichia coli thioredoxin as a scaffold, the histidine-mediated nucleophilic hydrolysis of p-nitrophenyl acetate as a model reaction, and the ORBIT protein design software to compute sequences, an active site scan identified two promising catalytic positions and surrounding active-site mutations required for substrate binding. Experimentally, both candidate protozymes demonstrated catalytic activity significantly above background. One of the proteins, PZD2, displayed "burst" phase kinetics at high substrate concentrations, consistent with the formation of a stable enzyme intermediate. The kinetic parameters of PZD2 are comparable to early catalytic Abs. But, unlike catalytic Ab design, our design procedure is independent of fold, suggesting a possible mechanism for examining the relationships between protein fold and the evolvability of protein function.

Proton and metal ion-dependent assembly of a model diiron protein

DF1 is a small, idealized model for carboxylate-bridged diiron proteins. This protein was designed to form a dimeric four-helix bundle with a dimetal ion-binding site near the center of the structure, and its crystal structure has confirmed that it adopts the intended conformation. However, the protein showed limited solubility in aqueous buffer, and access to its active site was blocked by two hydrophobic side chains. The sequence of DF1 has now been modified to provide a very soluble protein (DF2) that binds metal ions in a rapid and reversible manner. Furthermore, the DF2 protein shows significant ferroxidase activity, suggesting that its dimetal center is accessible to oxygen. The affinity of DF2 for various first-row divalent cations deviates from the Irving-Willliams series, suggesting that its structure imparts significant geometric preferences on the metal ion-binding site. Furthermore, in the absence of metal ions, the protein folds into a dimer with concomitant binding of two protons. The uptake of two protons is expected if the structure of the apo-protein is similar to that of the crystal structure of dizinc DF1. Thus, this result suggests that the active site of DF2 is retained in the absence of metal ions.

Design, synthesis, and characterization of a novel hemoprotein

Here we describe a synthetic protein (6H7H) designed to bind four heme groups via bis-histidine axial ligation. The hemes are designed to bind perpendicular to another in an orientation that mimics the relative geometry of the two heme a groups in the active site of cytochrome c oxidase. Our newly developed protein-design program, called CORE, was implemented in the design of this novel hemoprotein. Heme titration studies resolved four distinct K(D) values (K(D1) = 80 nM, K(D2) = 18 nM, K(D3) > or = 3 mM, K(D4) < or = 570 nM, with K(D3) x K(D4) = 1700); positive cooperativity in binding between the first and second heme, as well as substantial positive cooperativity between the third and forth heme, was observed. Chemical and thermal denaturation studies reveal a stable protein with native-like properties. Visible circular dichroism spectroscopy of holo-6H7H indicates excitonic coupling between heme groups. Further electrochemical and spectroscopic characterization of the holo-protein support a structure that is consistent with the predefined target structure.

Tailoring new enzyme functions by rational redesign

Site-directed mutagenesis is still a very efficient strategy to elaborate improved enzymes. Recently, advances have been made in developing rational strategies aimed at reshaping enzyme specificities and mechanisms, and at engineering biocatalysts through molecular assembling. These knowledge-based studies greatly benefit from the most recent computational analyses of enzyme structures and functions. The combination of rational and combinatorial methods opens up new vistas in the design of stable and efficient enzymes.

De novo design and structural characterization of proteins and metalloproteins

De novo protein design has recently emerged as an attractive approach for studying the structure and function of proteins. This approach critically tests our understanding of the principles of protein folding; only in de novo design must one truly confront the issue of how to specify a protein's fold and function. If we truly understand proteins, it should be possible to design receptors, enzymes, and ion channels from scratch. Further, as this understanding evolves and is further refined, it should be possible to design proteins and biomimetic polymers with properties unprecedented in nature.

Structure of a thioredoxin-like [2Fe-2S] ferredoxin from Aquifex aeolicus

The 2.3 A resolution crystal structure of a [2Fe-2S] cluster containing ferredoxin from Aquifex aeolicus reveals a thioredoxin-like fold that is novel among iron-sulfur proteins. The [2Fe-2S] cluster is located near the surface of the protein, at a site corresponding to that of the active-site disulfide bridge in thioredoxin. The four cysteine ligands are located near the ends of two surface loops. Two of these ligands can be substituted by non-native cysteine residues introduced throughout a stretch of the polypeptide chain that forms a protruding loop extending away from the cluster. The presence of homologs of this ferredoxin as components of more complex anaerobic and aerobic electron transfer systems indicates that this is a versatile fold for biological redox processes.

The aglycone specificity-determining sites are different in 2, 4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA)-glucosidase (Maize beta -glucosidase) and dhurrinase (Sorghum beta -glucosidase)

The maize beta-glucosidase isozyme Glu1 hydrolyzes a broad spectrum of substrates in addition to its natural substrate DIMBOAGlc (2-O-beta-d-glucopyranosyl-4-hydroxy-7-methoxy-1,4-benzoxazin-3-on e), whereas the sorghum beta-glucosidase isozyme Dhr1 hydrolyzes exclusively its natural substrate dhurrin (p-hydroxy-(S)-mandelonitrile-beta-d-glucose). To study the mechanism of substrate specificity further, eight chimeric beta-glucosidases were constructed by replacing peptide sequences within the C-terminal region of Glu1 with the homologous peptide sequences of Dhr1 or vice versa, where the two enzymes differ by 4 to 22 amino acid substitutions, depending on the length of the swapped regions. Five Glu1/Dhr1 chimeras hydrolyzed substrates that are hydrolyzed by both parental enzymes, including dhurrin, which is not hydrolyzed by Glu1. In contrast, three Dhr1/Glu1 chimeras hydrolyzed only dhurrin but with lower catalytic efficiency than Dhr1. Additional domain-swapping within the C-terminal domain of Glu1 showed that replacing the peptide (466)FAGFTERY(473) of Glu1 with the homologous peptide (462)SSGYTERF(469) of Dhr1 or replacing the peptide (481)NNNCTRYMKE(490) in Glu1 with the homologous peptide (477)ENGCERTMKR(486) of Dhr1 was sufficient to confer to Glu1 the ability to hydrolyze dhurrin. Data from various reciprocal chimeras, sequence comparisons, and homology modeling suggest that the Dhr1-specific Ser-462-Ser-463 and Phe-469 play a key role in dhurrin hydrolysis. Similar data suggest that DIMBOAGlc hydrolysis determinants are not located within the extreme 47-amino acid-long C-terminal domain of Glu1.

Retrostructural analysis of metalloproteins: application to the design of a minimal model for diiron proteins

De novo protein design provides an attractive approach for the construction of models to probe the features required for function of complex metalloproteins. The metal-binding sites of many metalloproteins lie between multiple elements of secondary structure, inviting a retrostructural approach to constructing minimal models of their active sites. The backbone geometries comprising the metal-binding sites of zinc fingers, diiron proteins, and rubredoxins may be described to within approximately 1 A rms deviation by using a simple geometric model with only six adjustable parameters. These geometric models provide excellent starting points for the design of metalloproteins, as illustrated in the construction of Due Ferro 1 (DF1), a minimal model for the Glu-Xxx-Xxx-His class of dinuclear metalloproteins. This protein was synthesized and structurally characterized as the di-Zn(II) complex by x-ray crystallography, by using data that extend to 2.5 A. This four-helix bundle protein is comprised of two noncovalently associated helix-loop-helix motifs. The dinuclear center is formed by two bridging Glu and two chelating Glu side chains, as well as two monodentate His ligands. The primary ligands are mostly buried in the protein interior, and their geometries are stabilized by a network of hydrogen bonds to second-shell ligands. In particular, a Tyr residue forms a hydrogen bond to a chelating Glu ligand, similar to a motif found in the diiron-containing R2 subunit of Escherichia coli ribonucleotide reductase and the ferritins. DF1 also binds cobalt and iron ions and should provide an attractive model for a variety of diiron proteins that use oxygen for processes including iron storage, radical formation, and hydrocarbon oxidation.

Rational design of nascent metalloenzymes

Understanding the early genesis of new enzymatic functions is one of the challenges in protein design, mechanistic enzymology, and molecular evolution. We have experimentally mimicked starting points in this process by introducing primitive iron and oxygen binding sites at various locations in thioredoxin, a small protein lacking metal centers, by using computational design. These rudimentary active sites show emerging enzymatic activities that select to varying degrees between different oxygen chemistries. Even within these nascent enzymes, mechanisms by which different reactions are controlled can be discerned. These involve both stabilizing and destabilizing interactions imposed on the metal center by the surrounding protein matrix.

Use of a non-rigid region in T4 lysozyme to design an adaptable metal-binding site

It is not easy to find candidate sites within a given protein where the geometry of the polypeptide chain matches that of metal-binding sites in known protein structures. By choosing a location in T4 lysozyme that is inherently flexible, it was possible to engineer a two-histidine site that binds different divalent cations. Crystallographic analysis shows that the geometry of binding of zinc is distorted tetrahedral while that of cobalt and nickel is octahedral. Insofar as spectroscopic data can be measured, they indicate that similar modes of coordination are retained in solution. The two substitutions, Thr21 --> His and Thr142 --> His, lie, respectively, on the surface of the N- and C-terminal domains on opposite sides of the active site cleft. The design takes advantage of hinge-bending motion which allows the binding site to adapt to the most favorable ligand geometry for the metal. Introduction of the two histidines increases the melting temperature of the protein by 2.0 degrees C at pH 7.4. Metal binding further increases the melting temperature, but only by a small amount (up to 1.5 degrees C). A third substitution, Gln141 --> His, which could act as a third ligand in principle, does not do so, demonstrating the difficulty in mimicking naturally occurring metal-binding sites.

From genes to protein structure and function: novel applications of computational approaches in the genomic era

The genome-sequencing projects are providing a detailed 'parts list' of life. A key to comprehending this list is understanding the function of each gene and each protein at various levels. Sequence-based methods for function prediction are inadequate because of the multifunctional nature of proteins. However, just knowing the structure of the protein is also insufficient for prediction of multiple functional sites. Structural descriptors for protein functional sites are crucial for unlocking the secrets in both the sequence and structural-genomics projects.

De novo protein design. I. In search of stability and specificity

We have developed a fully automated protein design strategy that works on the entire sequence of the protein and uses a full atom representation. At each step of the procedure, an all-atom model of the protein is built using the template protein structure and the current designed sequence. The energy of the model is used to drive a Monte Carlo optimization in sequence space: random moves are either accepted or rejected based on the Metropolis criterion. We rely on the physical forces that stabilize native protein structures to choose the optimum sequence. Our energy function includes van der Waals interactions, electrostatics and an environment free energy. Successful protein design should be specific and generate a sequence compatible with the template fold and incompatible with competing folds. We impose specificity by maintaining the amino acid composition constant, based on the random energy model. The specificity of the optimized sequence is tested by fold recognition techniques. Successful sequence designs for the B1 domain of protein G, for the lambda repressor and for sperm whale myoglobin are presented. We show that each additional term of the energy function improves the performance of our design procedure: the van der Waals term ensures correct packing, the electrostatics term increases the specificity for the correct native fold, and the environment solvation term ensures a correct pattern of buried hydrophobic and exposed hydrophilic residues. For the globin family, we show that we can design a protein sequence that is stable in the myoglobin fold, yet incompatible with the very similar hemoglobin fold.

Effect of hydrophobic core packing on sidechain dynamics

The effect of hydrophobic core packing on sidechain dynamics was analyzed by comparing the dynamics of wild-type (WT) ubiquitin to those of a variant which has seven core mutations. This variant, 1D7, was designed to resemble WT by having a well-packed core of similar volume, and we find that its overall level of dynamics is only subtly different from WT. However, the mutations caused a redistribution in the positions of core residues that are dynamic. This correlates with the tendency of these residues to populate unfavorable rotamers, suggesting that strain from poor sidechain conformations may promote increased flexibility as a mechanism to relieve unfavorable steric interactions. The results demonstrate that even when core volume is conserved, different packing arrangements in mutants can alter dynamic behavior.

Histidine placement in de novo-designed heme proteins

The effects of histidine residue placement in a de novo-designed four-alpha-helix bundle are investigated by placement of histidine residues at coiled coil heptad a positions in two distinct heptads and at each position within a single heptad repeat of our prototype heme protein maquette, [H10H24]2 [[Ac-CGGGELWKL x HEELLKK x FEELLKL x HEERLKK x L-CONH2]2]2 composed of a generic (alpha-SS-alpha)2 peptide architecture. The heme to peptide stoichiometry of variants of [H10H24]2 with either or both histidines on each helix replaced with noncoordinating alanine residues ([H10A24]2, [A10H24]2, and [A10A24]2) demonstrates the obligate requirement of histidine for biologically significant heme affinity. Variants of [A10A24]2, [[Ac-CGGGELWKL x AEELLKK x FEELLKL x AEERLKK x L-CONH2]2]2, containing a single histidine per helix in positions 9 to 15 were evaluated to verify the design based on molecular modeling. The bis-histidine site formed between heptad positions a at 10 and 10' bound ferric hemes with the highest affinity, Kd1 and Kd2 values of 1.5 and 800 nM, respectively. Placement of histidine at position 11 (heptad position b) resulted in a protein that bound a single heme with moderate affinity, Kd1 of 9.5 microM, whereas the other peptides had no measurable apparent affinity for ferric heme with Kd1 values >200 microM. The bis-histidine ligation of heme to [H10A24]2 and [H11A24]2 was confirmed by electron paramagnetic resonance spectroscopy. The protein design rules derived from this study, together with the narrow tolerances revealed, are applicable for improving future heme protein designs, for analyzing the results of randomized heme protein combinatorial libraries, as well as for implementation in automated protein design.

Hydrophobic core packing and protein design

Over the past few years, we have witnessed exciting advances in protein design. Several groups have reported success in the design of hydrophobic cores, and the principles developed in these studies have been recently applied to the full sequence design of a small protein motif and the design of a catalytically active metal center. These successes suggest that designing large, functional proteins in computero is more feasible than ever before.

The de novo design of a rubredoxin-like Fe site

A redox center similar to that of rubredoxin was designed into the 56 amino acid immunoglobulin binding B1 domain of Streptococcals protein G. The redox center in rubredoxin contains an iron ion tetrahedrally coordinated by four cysteine residues, [Fe(S-Cys)4](-1),(-2). The design criteria for the target site included taking backbone movements into account, tetrahedral metal-binding, and maintaining the structure and stability of the wild-type protein. The optical absorption spectrum of the Co(II) complex of the metal-binding variant is characteristic of tetrahedral chelation by four cysteine residues. Circular dichroism and nuclear magnetic resonance measurements reveal that the metal-free and Cd(II)-bound forms of the variant are folded correctly and are stable. The Fe(III) complex of the metal-binding mutant reproduces the optical and the electron paramagnetic resonance spectra of oxidized rubredoxin. This demonstrates that the engineered protein chelates Fe(III) in a tetrahedral array, and the resulting center is similar to that of oxidized rubredoxin.

Class-directed structure determination: foundation for a protein structure initiative

The recent sequencing of many complete genomes, combined with the development of methods that allow rapid structure determination for many proteins, has changed the way in which protein structure determinations can be approached. One-by-one determinations of individual protein structures will soon be augmented by class-directed structure analyses in which a group of proteins is targeted and structures of representative members are determined and used to represent the entire group. Such a shift in approach would be the foundation for a broad protein structure initiative targeting classes of proteins important for biotechnology and for a fundamental understanding of protein function.

From synthetic coiled coils to functional proteins: automated design of a receptor for the calmodulin-binding domain of calcineurin

A series of synthetic receptors capable of binding to the calmodulin-binding domain of calcineurin (CN393-414) was designed, synthesized and characterized. The design was accomplished by docking CN393-414 against a two-helix receptor, using an idealized three-stranded coiled coil as a starting geometry. The sequence of the receptor was chosen using a side-chain re-packing program, which employed a genetic algorithm to select potential binders from a total of 7.5x10(6) possible sequences. A total of 25 receptors were prepared, representing 13 sequences predicted by the algorithm as well as 12 related sequences that were not predicted. The receptors were characterized by CD spectroscopy, analytical ultracentrifugation, and binding assays. The receptors predicted by the algorithm bound CN393-414 with apparent dissociation constants ranging from 0.2 microM to >50 microM. Many of the receptors that were not predicted by the algorithm also bound to CN393-414. Methods to circumvent this problem and to improve the automated design of functional proteins are discussed.

The development of new biotechnologies using metalloprotein design

A variety of methodologies are under development to alter the behavior of existing metal centers or create entirely new sites within a protein framework in order to exploit the intrinsic chemical versatility of metals using the exquisite level of control that a protein matrix can exert to modulate their reactivity. Even at this relatively early stage, engineering of metal centers has led to the development of a number of emerging technologies with a wide variety of applications, including affinity purification of proteins, engineering of metal-mediated protein stability, control of protein activity, imaging and therapy, biosensors, and new catalysts.

Computer search algorithms in protein modification and design

The computer-aided design of protein sequences requires efficient search algorithms to handle the enormous combinatorial complexity involved. A variety of different algorithms have now been applied with some success. The choice of algorithm can influence the representation of the problem in several important ways--the discreteness of the configuration, the types of energy terms that can be used and the ability to find the global minimum energy configuration. The use of dead end elimination to design the complete sequence for a small protein motif and the use of genetic and mean-field algorithms to design hydrophobic cores for proteins represent the major themes of the past year.

The role of protonation and metal chelation preferences in defining the properties of mercury-binding coiled coils

To define the delicate interplay between metal chelation, protein folding and function in metalloproteins, a family of de novo-designed peptides was synthesized that self-assemble in aqueous solution to form two and three-stranded alpha-helical coiled coils. Each peptide contains a single Cys residue at an a or d position of the heptad repeat. Peptide association thus produces a Cys-rich coordination environment that has been used to bind Hg(II) ions. These peptides display a pH-dependent association, with trimers observed above the pKa of Glu side-chains and dimers below this value. Finite-difference Poisson-Boltzmann calculations suggest that the dimeric state decreases the unfavorable electrostatic interactions between positively charged Lys side-chains (relative to the trimer). The Cys-containing peptides bind Hg(II) in a position-dependent fashion. Cys at a positions form three-coordinate Hg complexes at high pH where the trimeric aggregation state predominates, and two-coordinate complexes at lower pH. A d position Cys, however, is only able to generate the two-coordinate complex, illustrating the difference in coordination geometry between the two positions in the coiled coil. The binding of Hg(II) was also shown to substantially increase the stability of the helical aggregates.

Hybrid enzymes: manipulating enzyme design

Hybrid enzymes are engineered to contain elements of two or more enzymes. Hybrid-enzyme approaches, by taking advantage of the vast array of enzymatic properties that nature has evolved, as well as the strategies that nature has used to evolve them, are becoming an increasingly important avenue for obtaining novel enzymes with desired activities and properties.

The dark side of dioxygen biochemistry

The cellular biochemistry of dioxygen is Janus-faced. The good side includes numerous enzyme-catalyzed reactions of dioxygen that occur in respiration and normal metabolism, while the dark side encompasses deleterious reactions of species derived from dioxygen that lead to damage of cellular components. These reactive oxygen species have historically been perceived almost exclusively as agents of the dark side, but it has recently become clear that they play beneficial roles as well.

Effect of four helix bundle topology on heme binding and redox properties

We have designed two alternative four helix bundle protein scaffold topologies for maquette construction to examine the effect of helix orientation on the heme binding and redox properties of our prototype heme protein maquette, (alpha-SS-alpha)2, previously described as H10H24 [Robertson, D. E., Farid, R. S., Moser, C. C., Mulholland, S. E., Pidikiti, R., Lear, J. D., Wand, A. J., DeGrado, W. F., and Dutton, P. L. (1994) Nature 368, 425]. Conversion of the disulfide-bridged di-alpha-helical monomer of (alpha-SS-alpha)2 into a single polypeptide chain results in topological reorientation of the helix dipoles and side chains within a 62 amino acid helix-loop-helix monomer, (alpha-l-alpha), which self-associates to form (alpha-l-alpha)2. Addition of an N-terminal cysteine residue to (alpha-l-alpha) with subsequent oxidation yields a 126 amino acid single molecule four helix bundle, (alpha-l-alpha-SS-alpha-l-alpha). Gel permeation chromatography demonstrated that (alpha-SS-alpha)2 and (alpha'-SS-alpha')2, a uniquely structured variant of the prototype, as well as (alpha-l-alpha)2 and (alpha'-l-alpha')2 assemble into distinct four helix bundles as designed, whereas (alpha-l-alpha-SS-alpha-l-alpha) elutes as a monomeric four alpha-helix bundle. Circular dichroism (CD) spectroscopy proves that these peptides are highly alpha-helical, and incorporation of four hemes has little effect on the helical content of the secondary structure. Four heme dissociation constants were evaluated by UV-visible spectroscopy and ranged from the 15 nM to 25 microM range for each of the peptides. The presence of Cotton effects in the visible CD illustrated that the hemes reside within the protein architecture. The equilibrium redox midpoint potentials (Em8) of the four bound hemes in each peptide are between -100 and -280 mV, as determined by redox potentiometry. The heme affinity and spectroelectrochemical properties of the hemes bound to (alpha-l-alpha)2 and (alpha-l-alpha-SS-alpha-l-alpha) are similar to those of the prototype, (alpha-SS-alpha)2, and to bis-histidine ligated b-type cytochromes, regardless of the global architectural changes imposed by these topological rearrangements. The hydrophobic cores of these peptides support local electrostatic fields which result in nativelike heme chromophore properties (spectroscopy, elevated reduction potentials, heme-heme charge interaction, and reactivity with exogenous diatomics) illustrating the utility of these non-native peptides in the study of metalloproteins.

The construction of metal centers in proteins by rational design

Metalloprotein properties result from the interplay between coordination requirements of the metal center, protein stability, and modulation of the metal center by the surrounding protein matrix. Simple metal centers, which exercise control over the protein by affecting stability or enzyme activity, have been created by rational design. Complex centers, which require control by the protein matrix, have also been constructed.

Engineering enzyme specificity

Protein engineering is the application of knowledge to design and alter protein function and structure. Although powerful methods, from specific to random, have been developed for the redesign of protein architecture, their successful application is dependent on the information known about the protein. This database of information is providing a foundation for establishing rules that govern enzyme-substrate interactions.

The dynamic properties of the M121H azurin metal site as studied by NMR of the paramagnetic Cu(II) and Co(II) metalloderivatives

The M121H azurin mutant in solution presents various species in equilibrium that can be detected and studied by 1H NMR of the Cu(II) and Co(II) paramagnetic metalloderivatives. In both cases up to three species are observed in slow exchange, the proportions of which are different for the two metalloderivatives. Above pH 5 the major species displays a tetrahedral coordination in which the His121 can be observed as a coordinated residue. Its metal site corresponds to a new type of site that is defined as a type 1.5 site. The second and third species resemble the wild type (type 1) azurin and, above pH 4.5, they are present only at a low concentration. At low pH a protonation process increases the proportion of both type 1 species at the expense of the type 1.5 species. This process, characterized by a pKa = 4.3, is assigned to the protonation of His121. At high pH the NMR spectrum of the Co(II)-M121H azurin experiences an additional transition, which is not observed in the case of the Cu(II) protein. The dynamic properties of the M121H metal site appear to be related to changes in the coordination geometry and the strength of the axial interaction between the Ndelta1 (His121) and the metal.

Synthesis of novel proteins

De novo and rational protein design are progressing towards the chemical synthesis of proteins with pre-selected structure and function. The data illustrate diverse experimental and computational approaches which test our comprehension of protein structure, hydrophobic core packing and global stability, especially of coiled-coil proteins. The incorporation of biological cofactors, including hemes, as well as active sites, such as that of iron superoxide dismutase, into designed proteins provides an exciting next step towards the synthesis of proteins with enzymatic function.

The rational design and construction of a cuboidal iron-sulfur protein

Rational protein design is an emerging approach for testing general theories of protein chemistry through the creation of new structures and functions. Here we present the first successful introduction by rational design of a [Fe4S4] cuboidal cluster into the hydrophobic core of Escherichia coli thioredoxin, a protein normally devoid of metal centers. Cuboidal [Fe4S4] is one of the stable forms of self-assembled iron-sulfur clusters that are thought to represent some of the earliest evolved biological redox centers. [Fe4S4] clusters have been recruited for use in a variety of proteins whose functions are central to many of the major biochemical processes ranging from simple soluble electron-transfer agents, to membrane-bound components of electron-transfer chains, to electron reservoirs in complex metalloenzymes such as nitrogenase. By situating an [Fe4S4] cluster into a protein environment not previously adapted by evolution we can explore the factors by which their activity is modulated by the protein matrix.