Tertiary templates for the design of diiron proteins

Rational design of photosynthetic reaction center protein maquettes

New technologies for efficient solar-to-fuel energy conversion will help facilitate a global shift from dependence on fossil fuels to renewable energy. Nature uses photosynthetic reaction centers to convert photon energy into a cascade of electron-transfer reactions that eventually produce chemical fuel. The design of new reaction centers de novo deepens our understanding of photosynthetic charge separation and may one day allow production of biofuels with higher thermodynamic efficiency than natural photosystems. Recently, we described the multi-step electron-transfer activity of a designed reaction center maquette protein (the RC maquette), which can assemble metal ions, tyrosine, a Zn tetrapyrrole, and heme into an electron-transport chain. Here, we detail our modular strategy for rational protein design and show that the intended RC maquette design agrees with crystal structures in various states of assembly. A flexible, dynamic apo-state collapses by design into a more ordered holo-state upon cofactor binding. Crystal structures illustrate the structural transitions upon binding of different cofactors. Spectroscopic assays demonstrate that the RC maquette binds various electron donors, pigments, and electron acceptors with high affinity. We close with a critique of the present RC maquette design and use electron-tunneling theory to envision a path toward a designed RC with a substantially higher thermodynamic efficiency than natural photosystems.

Catalysis and Electron Transfer in De Novo Designed Metalloproteins

One of the hallmark advances in our understanding of metalloprotein function is showcased in our ability to design new, non-native, catalytically active protein scaffolds. This review highlights progress and milestone achievements in the field of de novo metalloprotein design focused on reports from the past decade with special emphasis on de novo designs couched within common subfields of bioinorganic study: heme binding proteins, monometal- and dimetal-containing catalytic sites, and metal-containing electron transfer sites. Within each subfield, we highlight several of what we have identified as significant and important contributions to either our understanding of that subfield or de novo metalloprotein design as a discipline. These reports are placed in context both historically and scientifically. General suggestions for future directions that we feel will be important to advance our understanding or accelerate discovery are discussed.

Computational design of a sensitive, selective phase-changing sensor protein for the VX nerve agent

The VX nerve agent is one of the deadliest chemical warfare agents. Specific, sensitive, real-time detection methods for this neurotoxin have not been reported. The creation of proteins that use biological recognition to fulfill these requirements using directed evolution or library screening methods has been hampered because its toxicity makes laboratory experimentation extraordinarily expensive. A pair of VX-binding proteins were designed using a supercharged scaffold that couples a large-scale phase change from unstructured to folded upon ligand binding, enabling fully internal binding sites that present the maximum surface area possible for high affinity and specificity in target recognition. Binding site residues were chosen using a new distributed evolutionary algorithm implementation in protCAD. Both designs detect VX at parts per billion concentrations with high specificity. Computational design of fully buried molecular recognition sites, in combination with supercharged phase-changing chassis proteins, enables the ready development of a new generation of small-molecule biosensors.

De novo metalloprotein design

Natural metalloproteins perform many functions - ranging from sensing to electron transfer and catalysis - in which the position and property of each ligand and metal, is dictated by protein structure. De novo protein design aims to define an amino acid sequence that encodes a specific structure and function, providing a critical test of the hypothetical inner workings of (metallo)proteins. To date, de novo metalloproteins have used simple, symmetric tertiary structures - uncomplicated by the large size and evolutionary marks of natural proteins - to interrogate structure-function hypotheses. In this Review, we discuss de novo design applications, such as proteins that induce complex, increasingly asymmetric ligand geometries to achieve function, as well as the use of more canonical ligand geometries to achieve stability. De novo design has been used to explore how proteins fine-tune redox potentials and catalyse both oxidative and hydrolytic reactions. With an increased understanding of structure-function relationships, functional proteins including O2-dependent oxidases, fast hydrolases, and multi-proton/multi-electron reductases, have been created. In addition, proteins can now be designed using xeno-biological metals or cofactors and principles from inorganic chemistry to derive new-to-nature functions. These results and the advances in computational protein design suggest a bright future for the de novo design of diverse, functional metalloproteins.

A de novo binuclear zinc enzyme with DNA cleavage activity

Metallohydrolases are broadly used throughout biology, often to catalyze the degradation of macromolecules such as DNA and proteins. Many of these enzymes function with zinc in their active site, and an important subset of these enzymes utilize a binuclear zinc active site. Mimics of these enzymes have been developed, some of which catalyze the digestion of DNA. However, the majority of the mimics that utilize zinc are small molecules, and most are mononuclear. Herein, we report DNA cleavage activity by the de novo designed Due Ferri single-chain (DFsc) protein containing a binuclear zinc active site. This binuclear zinc-protein complex is able to digest plasmid DNA at rates up to 50 ng/h, and these cleavage rates are affected by changes to amino acid residues near the zinc-binding site. These results indicate that the DFsc scaffold is a good model system to carry out careful structure-function relationship studies to understand key structural features that influence reactivity in natural binuclear zinc hydrolases, as it is the first report of a binuclear model system in a protein scaffold.

Designing heterotropically activated allosteric conformational switches using supercharging

Heterotropic allosteric activation of protein function, in which binding of one ligand thermodynamically activates the binding of another, different ligand or substrate, is a fundamental control mechanism in metabolism and as such has been a long-aspired capability in protein design. Here we show that greatly increasing the magnitude of a protein's net charge using surface supercharging transforms that protein into an allosteric ligand- and counterion-gated conformational molecular switch. To demonstrate this we first modified the designed helical bundle hemoprotein H4, creating a highly charged protein which both unfolds reversibly at low ionic strength and undergoes the ligand-induced folding transition commonly observed in signal transduction by intrinsically disordered proteins in biology. As a result of the high surface-charge density, ligand binding to this protein is allosterically activated up to 1,300-fold by low concentrations of divalent cations and the polyamine spermine. To extend this process further using a natural protein, we similarly modified Escherichia coli cytochrome b ₅₆₂ and the resulting protein behaves in a like manner. These simple model systems not only establish a set of general engineering principles which can be used to convert natural and designed soluble proteins into allosteric molecular switches useful in biodesign, sensing, and synthetic biology, the behavior we have demonstrated--functional activation of supercharged intrinsically disordered proteins by low concentrations of multivalent ions--may be a control mechanism utilized by Nature which has yet to be appreciated.

De novo protein design, a retrospective

Proteins are molecular machines whose function depends on their ability to achieve complex folds with precisely defined structural and dynamic properties. The rational design of proteins from first-principles, or de novo, was once considered to be impossible, but today proteins with a variety of folds and functions have been realized. We review the evolution of the field from its earliest days, placing particular emphasis on how this endeavor has illuminated our understanding of the principles underlying the folding and function of natural proteins, and is informing the design of macromolecules with unprecedented structures and properties. An initial set of milestones in de novo protein design focused on the construction of sequences that folded in water and membranes to adopt folded conformations. The first proteins were designed from first-principles using very simple physical models. As computers became more powerful, the use of the rotamer approximation allowed one to discover amino acid sequences that stabilize the desired fold. As the crystallographic database of protein structures expanded in subsequent years, it became possible to construct proteins by assembling short backbone fragments that frequently recur in Nature. The second set of milestones in de novo design involves the discovery of complex functions. Proteins have been designed to bind a variety of metals, porphyrins, and other cofactors. The design of proteins that catalyze hydrolysis and oxygen-dependent reactions has progressed significantly. However, de novo design of catalysts for energetically demanding reactions, or even proteins that bind with high affinity and specificity to highly functionalized complex polar molecules remains an importnant challenge that is now being achieved. Finally, the protein design contributed significantly to our understanding of membrane protein folding and transport of ions across membranes. The area of membrane protein design, or more generally of biomimetic polymers that function in mixed or non-aqueous environments, is now becoming increasingly possible.

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.

Amino acid modifications for conformationally constraining naturally occurring and engineered peptide backbones: Insights from the Protein Data Bank

Extensive efforts invested in understanding the rules of protein folding are now being applied, with good effect, in de novo design of proteins/peptides. For proteins containing standard α-amino acids alone, knowledge derived from experimentally determined three-dimensional (3D) structures of proteins and biologically active peptides are available from the Protein Data Bank (PDB), and the Cambridge Structural Database (CSD). These help predict and design protein structures, with reasonable confidence. However, our knowledge of 3D structures of biomolecules containing backbone modified amino acids is still evolving. A major challenge in de novo protein/peptide design concerns the engineering of conformationally constrained molecules with specific structural elements and chemical groups appropriately positioned for biological activity. This review explores four classes of amino acid modifications that constrain protein/peptide backbone structure. Systematic analysis of peptidic molecule structures (eg, bioactive peptides, inhibitors, antibiotics, and designed molecules), containing these backbone-modified amino acids, found in the PDB and CSD are discussed. The review aims to provide structure-function insights that will guide future design of proteins/peptides.

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.

Fast, cheap and out of control--Insights into thermodynamic and informatic constraints on natural protein sequences from de novo protein design

The accumulated results of thirty years of rational and computational de novo protein design have taught us important lessons about the stability, information content, and evolution of natural proteins. First, de novo protein design has complicated the assertion that biological function is equivalent to biological structure - demonstrating the capacity to abstract active sites from natural contexts and paste them into non-native topologies without loss of function. The structure-function relationship has thus been revealed to be either a generality or strictly true only in a local sense. Second, the simplification to "maquette" topologies carried out by rational protein design also has demonstrated that even sophisticated functions such as conformational switching, cooperative ligand binding, and light-activated electron transfer can be achieved with low-information design approaches. This is because for simple topologies the functional footprint in sequence space is enormous and easily exceeds the number of structures which could have possibly existed in the history of life on Earth. Finally, the pervasiveness of extraordinary stability in designed proteins challenges accepted models for the "marginal stability" of natural proteins, suggesting that there must be a selection pressure against highly stable proteins. This can be explained using recent theories which relate non-equilibrium thermodynamics and self-replication. This article is part of a Special Issue entitled Biodesign for Bioenergetics--The design and engineering of electronc transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Structural principles for computational and de novo design of 4Fe-4S metalloproteins

Iron-sulfur centers in metalloproteins can access multiple oxidation states over a broad range of potentials, allowing them to participate in a variety of electron transfer reactions and serving as catalysts for high-energy redox processes. The nitrogenase FeMoCO cluster converts di-nitrogen to ammonia in an eight-electron transfer step. The 2(Fe4S4) containing bacterial ferredoxin is an evolutionarily ancient metalloprotein fold and is thought to be a primordial progenitor of extant oxidoreductases. Controlling chemical transformations mediated by iron-sulfur centers such as nitrogen fixation, hydrogen production as well as electron transfer reactions involved in photosynthesis are of tremendous importance for sustainable chemistry and energy production initiatives. As such, there is significant interest in the design of iron-sulfur proteins as minimal models to gain fundamental understanding of complex natural systems and as lead-molecules for industrial and energy applications. Herein, we discuss salient structural characteristics of natural iron-sulfur proteins and how they guide principles for design. Model structures of past designs are analyzed in the context of these principles and potential directions for enhanced designs are presented, and new areas of iron-sulfur protein design are proposed. This article is part of a Special issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, protein networks, edited by Ronald L. Koder and J.L Ross Anderson.

Metal Stabilization of Collagen and de Novo Designed Mimetic Peptides

We explore the design of metal binding sites to modulate triple-helix stability of collagen and collagen-mimetic peptides. Globular proteins commonly utilize metals to connect tertiary structural elements that are well separated in sequence, constraining structure and enhancing stability. It is more challenging to engineer structural metals into fibrous protein scaffolds, which lack the extensive tertiary contacts seen in globular proteins. In the collagen triple helix, the structural adjacency of the carboxy-termini of the three chains makes this region an attractive target for introducing metal binding sites. We engineered His3 sites based on structural modeling constraints into a series of designed homotrimeric and heterotrimeric peptides, assessing the capacity of metal binding to improve stability and in the case of heterotrimers, affect specificity of assembly. Notable enhancements in stability for both homo- and heteromeric systems were observed upon addition of zinc(II) and several other metal ions only when all three histidine ligands were present. Metal binding affinities were consistent with the expected Irving-Williams series for imidazole. Unlike other metals tested, copper(II) also bound to peptides lacking histidine ligands. Acetylation of the peptide N-termini prevented copper binding, indicating proline backbone amide metal-coordination at this site. Copper similarly stabilized animal extracted Type I collagen in a metal-specific fashion, highlighting the potential importance of metal homeostasis within the extracellular matrix.

Crystal Structure of the Peroxo-diiron(III) Intermediate of Deoxyhypusine Hydroxylase, an Oxygenase Involved in Hypusination

Deoxyhypusine hydroxylase (DOHH) is a non-heme diiron enzyme involved in the posttranslational modification of a critical lysine residue of eukaryotic translation initiation factor 5A (eIF-5A) to yield the unusual amino acid residue hypusine. This modification is essential for the role of eIF-5A in translation and in nuclear export of a group of specific mRNAs. The diiron center of human DOHH (hDOHH) forms a peroxo-diiron(III) intermediate (hDOHHperoxo) when its reduced form reacts with O2. hDOHHperoxo has a lifetime exceeding that of the peroxo intermediates of other diiron enzymes by several orders of magnitude. Here we report the 1.7-Å crystal structures of hDOHHperoxo and a complex with glycerol. The structure of hDOHHperoxo reveals the presence of a μ-1,2-peroxo-diiron(III) species at the active site. Augmented by UV/Vis and Mössbauer spectroscopic studies, the crystal structures offer explanations for the extreme longevity of hDOHHperoxo and illustrate how the enzyme specifically recognizes its only substrate, deoxyhypusine-eIF-5A.

De novo design of a transmembrane Zn²⁺-transporting four-helix bundle

The design of functional membrane proteins from first principles represents a grand challenge in chemistry and structural biology. Here, we report the design of a membrane-spanning, four-helical bundle that transports first-row transition metal ions Zn(2+) and Co(2+), but not Ca(2+), across membranes. The conduction path was designed to contain two di-metal binding sites that bind with negative cooperativity. X-ray crystallography and solid-state and solution nuclear magnetic resonance indicate that the overall helical bundle is formed from two tightly interacting pairs of helices, which form individual domains that interact weakly along a more dynamic interface. Vesicle flux experiments show that as Zn(2+) ions diffuse down their concentration gradients, protons are antiported. These experiments illustrate the feasibility of designing membrane proteins with predefined structural and dynamic properties.

Cys-scanning disulfide crosslinking and bayesian modeling probe the transmembrane signaling mechanism of the histidine kinase, PhoQ

Bacteria transduce signals across the membrane using two-component systems (TCSs), consisting of a membrane-spanning sensor histidine kinase and a cytoplasmic response regulator. In gram-negative bacteria, the PhoPQ TCS senses cations and antimicrobial peptides, yet little is known about the structural changes involved in transmembrane signaling. We construct a model of PhoQ signal transduction using Bayesian inference, based on disulfide crosslinking data and homologous crystal structures. The data are incompatible with a single conformation but are instead consistent with two interconverting structures. These states differ in membrane depth of the periplasmic acidic patch and the reciprocal displacement of diagonal helices along the dimer interface. Studies of multiple histidine kinases suggest this repacking might be a common mode of signal transduction in sensor His-kinase receptors. Because a similar scissors model has been ruled out in CheA-linked chemoreceptors, the evidence suggests that sensor His-kinase and CheA-linked receptors possess different signaling mechanisms.

Computational design of metalloproteins

A number of design strategies exist for the development of novel metalloproteins. These strategies often exploit the inherent symmetry of metal coordination and local topology. Computational design of metal binding sites in flexible regions of proteins is challenging as the number of conformational degrees of freedom is significantly increased. Additionally, without pre-organization of the primary shell ligands by the protein fold, metal binding sites can rearrange according to the coordination constraints of the metal center. Examples of metal incorporation into existing folds, full fold design exploiting symmetry, and fold design in asymmetric scaffolds are presented.

Computational approaches for rational design of proteins with novel functionalities

Proteins are the most multifaceted macromolecules in living systems and have various important functions, including structural, catalytic, sensory, and regulatory functions. Rational design of enzymes is a great challenge to our understanding of protein structure and physical chemistry and has numerous potential applications. Protein design algorithms have been applied to design or engineer proteins that fold, fold faster, catalyze, catalyze faster, signal, and adopt preferred conformational states. The field of de novo protein design, although only a few decades old, is beginning to produce exciting results. Developments in this field are already having a significant impact on biotechnology and chemical biology. The application of powerful computational methods for functional protein designing has recently succeeded at engineering target activities. Here, we review recently reported de novo functional proteins that were developed using various protein design approaches, including rational design, computational optimization, and selection from combinatorial libraries, highlighting recent advances and successes.

Non-heme manganese catalase--the 'other' catalase

Non-heme manganese catalases are widely distributed over microbial life and represent an environmentally important alternative to heme-containing catalases in antioxidant defense. Manganese catalases contain a binuclear manganese complex as their catalytic active site rather than a heme, and cycle between Mn(2)(II,II) and Mn(2)(III,III) states during turnover. X-ray crystallography has revealed the key structural elements of the binuclear manganese active site complex that can serve as the starting point for computational studies on the protein. Four manganese catalase enzymes have been isolated and characterized, and the enzyme appears to have a broad phylogenetic distribution including both bacteria and archae. More than 100 manganese catalase genes have been annotated in genomic databases, although the assignment of many of these putative manganese catalases needs to be experimentally verified. Iron limitation, exposure to low levels of peroxide stress, thermostability and cyanide resistance may provide the biological and environmental context for the occurrence of manganese catalases.

Salt bridges: geometrically specific, designable interactions

Salt bridges occur frequently in proteins, providing conformational specificity and contributing to molecular recognition and catalysis. We present a comprehensive analysis of these interactions in protein structures by surveying a large database of protein structures. Salt bridges between Asp or Glu and His, Arg, or Lys display extremely well-defined geometric preferences. Several previously observed preferences are confirmed, and others that were previously unrecognized are discovered. Salt bridges are explored for their preferences for different separations in sequence and in space, geometric preferences within proteins and at protein-protein interfaces, co-operativity in networked salt bridges, inclusion within metal-binding sites, preference for acidic electrons, apparent conformational side chain entropy reduction on formation, and degree of burial. Salt bridges occur far more frequently between residues at close than distant sequence separations, but, at close distances, there remain strong preferences for salt bridges at specific separations. Specific types of complex salt bridges, involving three or more members, are also discovered. As we observe a strong relationship between the propensity to form a salt bridge and the placement of salt-bridging residues in protein sequences, we discuss the role that salt bridges might play in kinetically influencing protein folding and thermodynamically stabilizing the native conformation. We also develop a quantitative method to select appropriate crystal structure resolution and B-factor cutoffs. Detailed knowledge of these geometric and sequence dependences should aid de novo design and prediction algorithms.

Probing designability via a generalized model of helical bundle geometry

Because the space of folded protein structures is highly degenerate, with recurring secondary and tertiary motifs, methods for representing protein structure in terms of collective physically relevant coordinates are of great interest. By collapsing structural diversity to a handful of parameters, such methods can be used to delineate the space of designable structures (i.e., conformations that can be stabilized with a large number of sequences)-a crucial task for de novo protein design. We first demonstrate this on natural α-helical coiled coils using the Crick parameterization. We show that over 95% of known coiled-coil structures are within  1-Å C(α) root mean square deviation of a Crick-ideal backbone. Derived parameters show that natural geometric space of coiled coils is highly restricted and can be represented by "allowed" conformations amidst a potential continuum of conformers. Allowed structures have (1) restricted axial offsets between helices, which differ starkly between parallel and anti-parallel structures; (2) preferred superhelical radii, which depend linearly on the oligomerization state; (3) pronounced radius-dependent a- and d-position amino acid propensities; and (4) discrete angles of rotation of helices about their axes, which are surprisingly independent of oligomerization state or orientation. In all, we estimate the space of designable coiled-coil structures to be reduced at least 160-fold relative to the space of geometrically feasible structures. To extend the benefits of structural parameterization to other systems, we developed a general mathematical framework for parameterizing arbitrary helical structures, which reduces to the Crick parameterization as a special case. The method is successfully validated on a set of non-coiled-coil helical bundles, frequent in channels and transporter proteins, which show significant helix bending but not supercoiling. Programs for coiled-coil parameter fitting and structure generation are provided via a web interface at http://www.gevorggrigoryan.com/cccp/, and code for generalized helical parameterization is available upon request.

Designing artificial enzymes by intuition and computation

The rational design of artificial enzymes, either by applying physico-chemical intuition of protein structure and function or with the aid of computational methods, is a promising area of research with the potential to tremendously impact medicine, industrial chemistry and energy production. Designed proteins also provide a powerful platform for dissecting enzyme mechanisms of natural systems. Artificial enzymes have come a long way from simple α-helical peptide catalysts to proteins that facilitate multistep chemical reactions designed by state-of-the-art computational methods. Looking forward, we examine strategies employed by natural enzymes that could be used to improve the speed and selectivity of artificial catalysts.

Geometric constraints for porphyrin binding in helical protein binding sites

Helical bundles which bind heme and porphyrin cofactors have been popular targets for cofactor-containing de novo protein design. By analyzing a highly nonredundant subset of the protein databank we have determined a rotamer distribution for helical histidines bound to heme cofactors. Analysis of the entire nonredundant database for helical sequence preferences near the ligand histidine demonstrated little preference for amino acid side chain identity, size, or charge. Analysis of the database subdivided by ligand histidine rotamer, however, reveals strong preferences in each case, and computational modeling illuminates the structural basis for some of these findings. The majority of the rotamer distribution matches that predicted by molecular simulation of a single porphyrin-bound histidine residue placed in the center of an all-alanine helix, and the deviations explain two prominent features of natural heme protein binding sites: heme distortion in the case of the cytochromes C in the m166 histidine rotamer, and a highly prevalent glycine residue in the t73 histidine rotamer. These preferences permit derivation of helical consensus sequence templates which predict optimal side chain-cofactor packing interactions for each rotamer. These findings thus promise to guide future design endeavors not only in the creation of higher affinity heme and porphyrin binding sites, but also in the direction of bound cofactor geometry.

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.

Statistical theory for protein ensembles with designed energy landscapes

Combinatorial protein libraries provide a promising route to investigate the determinants and features of protein folding and to identify novel folding amino acid sequences. A library of sequences based on a pool of different monomer types are screened for folding molecules, consistent with a particular foldability criterion. The number of sequences grows exponentially with the length of the polymer, making both experimental and computational tabulations of sequences infeasible. Herein a statistical theory is extended to specify the properties of sequences having particular values of global energetic quantities that specify their energy landscape. The theory yields the site-specific monomer probabilities. A foldability criterion is derived that characterizes the properties of sequences by quantifying the energetic separation of the target state from low-energy states in the unfolded ensemble and the fluctuations of the energies in the unfolded state ensemble. For a simple lattice model of proteins, excellent agreement is observed between the theory and the results of exact enumeration. The theory may be used to provide a quantitative framework for the design and interpretation of combinatorial experiments.

Metalloprotein and metallo-DNA/RNAzyme design: current approaches, success measures, and future challenges

Specific metal-binding sites have been found in not only proteins but also DNA and RNA molecules. Together these metalloenzymes consist of a major portion of the enzyme family and can catalyze some of the most difficult biological reactions. Designing these metalloenzymes can be both challenging and rewarding because it can provide deeper insights into the structure and function of proteins and cheaper and more stable alternatives for biochemical and biotechnological applications. Toward this goal, both rational and combinatorial approaches have been used. The rational approach is good for designing metalloenzymes that are well characterized, such as heme proteins, while the combinatorial approach is better at designing those whose structures are poorly understood, such as metallo-DNA/RNAzymes. Among the rational approaches, de novo design is at its best when metal-binding sites reside in a scaffold whose structure has been designed de novo (e.g., alpha-helical bundles). Otherwise, design using native scaffolds can be equally effective, allowing more choices of scaffolds whose structural stability is often more resistant to multiple mutations. In addition, computational and empirical designs have both enjoyed successes. Because of the limitation in defining structural parameters for metal-binding sites, a computational approach is restricted to mostly metal-binding sites that are well defined, such as mono- or homonuclear centers. An empirical approach, even though it is less restrictive in the metal-binding sites to be designed, depends heavily on one's knowledge and choice of templates and targets. An emerging approach is a combination of both computational and empirical approaches. The success of these approaches can be measured not only by three-dimensional structural comparison between the designed and target enzymes but also by the total amount of insight obtained from the design process and studies of the designed enzymes. One of the biggest advantages of designed metalloenzymes is the potential of placing two different metal-binding sites in the same protein framework for comparison. A final measure of success is how one can utilize the insight gained from the intellectual exercise to design new metalloenzymes, including those with unprecedented structures and functions. Future challenges include designing more complex metalloenzymes such as heteronuclear metal centers with strong nanomolar or better affinities. A key to meeting this challenge is to focus on the design of not only primary but also secondary coordination spheres using a combination of improved computer programs, experimental design, and high-resolution crystallography.

Computer-based design of novel protein structures

Over the past 10 years there has been tremendous success in the area of computational protein design. Protein design software has been used to stabilize proteins, solubilize membrane proteins, design intermolecular interactions, and design new protein structures. A key motivation for these studies is that they test our understanding of protein energetics and structure. De novo design of novel structures is a particularly rigorous test because the protein backbone must be designed in addition to the amino acid side chains. A priori it is not guaranteed that the target backbone is even designable. To address this issue, researchers have developed a variety of methods for generating protein-like scaffolds and for optimizing the protein backbone in conjunction with the amino acid sequence. These protocols have been used to design proteins from scratch and to explore sequence space for naturally occurring protein folds.

Intelligent design: the de novo engineering of proteins with specified functions

One of the principal successes of de novo protein design has been the creation of small, robust protein-cofactor complexes which can serve as simplified models, or maquettes, of more complicated multicofactor protein complexes commonly found in nature. Different maquettes, generated by us and others, recreate a variety of aspects, or functional elements, recognized as parts of natural enzyme function. The current challenge is to both expand the palette of functional elements and combine and/or integrate them in recreating familiar enzyme activities or generating novel catalysis in the simplest protein scaffolds.

Artificial diiron proteins: from structure to function

De novo protein design provides an attractive approach for the construction of models to probe the features required for the function of complex metalloproteins. These minimal models contain the essential elements believed necessary for activity of the protein. In this article, we summarize the design, structure determination, and functional properties of a family of artificial diiron proteins.

Complementary nucleobase interaction enhances peptide-peptide recognition and self-replicating catalysis

The availability of the complementary interaction of nucleobases for influencing the formation of peptide architectures was explored. Nucleobases were incorporated as additional recognition elements in coiled-coil peptides by employing nucleobase amino acids (NBAs), which are artificial L-alpha-amino gamma-nucleobase-butyric acids. The effect of the base-pair interaction on intermolecular recognition between peptides was evaluated through a self-replication reaction. The self-replication reactions of the peptides with complementary base pairs such as thymine-adenine or guanine-cytosine at the g-g' heptad positions were accelerated in comparison with those of the peptides with mismatched base pairs or without nucleobases. However, thymine-adenine pairs at the e-e' positions did not enhance the self-replication. In the presence of a denaturant, the enhancement effects of complementary base pairs on the reaction disappeared. Thermal denaturation studies showed that the thymine-adenine pairs contributed to stabilization of the coiled-coil structure and that the pairs at the g-g' positions were more effective than those at the e-e' positions. The peptide-peptide interaction was reinforced by complementary nucleobase interactions appropriately arranged in the peptide structure; these led to acceleration of the self-replication reactions.

Designing proteins from the inside out

Globular proteins are characterized by the specific and tight packing of hydrophobic side-chains in the so-called "hydrophobic core." Formation of the core is key in folding, stabilization, and conformational specificity. The critical role of hydrophobic cores in maintaining the highly ordered structures present in natural proteins justifies the tremendous efforts devoted to their redesign. Both experimental and computational combinatorial-based approaches have been reported in the last years as powerful protein design tools. These manage to explore large regions of the sequence/conformational space, allowing the search for alternative protein core arrangements displaying native-like properties. The overall results obtained from core design projects have contributed significantly to our present knowledge of protein folding and function. In addition, core design has worked as a benchmark for the development of ambitious protein design projects that nowadays are allowing the de novo design of novel protein structures and functions.

Computational design and characterization of a monomeric helical dinuclear metalloprotein

The de novo design of di-iron proteins is an important step towards understanding the diversity of function among this complex family of metalloenzymes. Previous designs of due ferro (DF) proteins have resulted in tetrameric and dimeric four-helix bundles having crystallographically well-defined structures and active-site geometries. Here, the design and characterization of DFsc, a 114 residue monomeric four-helix bundle, is presented. The backbone was modeled using previous oligomeric structures and appropriate inter-helical turns. The identities of 26 residues were predetermined, including the primary and secondary ligands in the active site, residues involved in active site accessibility, and the gamma beta gamma beta turn between helices 2 and 3. The remaining 88 amino acid residues were determined using statistical computer aided design, which is based upon a recent statistical theory of protein sequences. Rather than sampling sequences, the theory directly provides the site-specific amino acid probabilities, which are then used to guide sequence design. The resulting sequence (DFsc) expresses well in Escherichia coli and is highly soluble. Sedimentation studies confirm that the protein is monomeric in solution. Circular dichroism spectra are consistent with the helical content of the target structure. The protein is structured in both the apo and the holo forms, with the metal-bound form exhibiting increased stability. DFsc stoichiometrically binds a variety of divalent metal ions, including Zn(II), Co(II), Fe(II), and Mn(II), with micromolar affinities. 15N HSQC NMR spectra of both the apo and Zn(II) proteins reveal excellent dispersion with evidence of a significant structural change upon metal binding. DFsc is then a realization of complete de novo design, where backbone structure, activity, and sequence are specified in the design process.

De novo designed cyclic-peptide heme complexes

The structural characterization of de novo designed metalloproteins together with determination of chemical reactivity can provide a detailed understanding of the relationship between protein structure and functional properties. Toward this goal, we have prepared a series of cyclic peptides that bind to water-soluble metalloporphyrins (FeIII and CoIII). Neutral and positively charged histidine-containing peptides bind with a high affinity, whereas anionic peptides bind only weakly to the negatively charged metalloporphyrin. Additionally, it was found that the peptide becomes helical only in the presence of the metalloporphyrin. CD experiments confirm that the metalloporphyrin binds specific cyclic peptides with high affinity and with isodichroic behavior. Thermal unfolding experiments show that the complex has "native-like" properties. Finally, NMR spectroscopy produced well dispersed spectra and experimental restraints that provide a high-resolution solution structure of the complexed peptide.

Evidence linking the Pseudomonas oleovorans alkane omega-hydroxylase, an integral membrane diiron enzyme, and the fatty acid desaturase family

Pseudomonas oleovorans alkane omega-hydroxylase (AlkB) is an integral membrane diiron enzyme that shares a requirement for iron and oxygen for activity in a manner similar to that of the non-heme integral membrane desaturases, epoxidases, acetylenases, conjugases, ketolases, decarbonylase and methyl oxidases. No overall sequence similarity is detected between AlkB and these desaturase-like enzymes by computer algorithms; however, they do contain a series of histidine residues in a similar relative positioning with respect to hydrophobic regions thought to be transmembrane domains. To test whether these conserved histidine residues are functionally equivalent to those of the desaturase-like enzymes we used scanning alanine mutagenesis to test if they are essential for activity of AlkB. These experiments show that alanine substitution of any of the eight conserved histidines results in complete inactivation, whereas replacement of three non-conserved histidines in close proximity to the conserved residues, results in only partial inactivation. These data provide the first experimental support for the hypotheses: (i) that the histidine motif in AlkB is equivalent to that in the desaturase-like enzymes and (ii) that the conserved histidine residues play a vital role such as coordinating the Fe ions comprising the diiron active site.

Preorganization of molecular binding sites in designed diiron proteins

De novo protein design provides an attractive approach to critically test the features that are required for metalloprotein structure and function. Previously we designed and crystallographically characterized an idealized dimeric model for the four-helix bundle class of diiron and dimanganese proteins [Dueferri 1 (DF1)]. Although the protein bound metal ions in the expected manner, access to its active site was blocked by large bulky hydrophobic residues. Subsequently, a substrate-access channel was introduced proximal to the metal-binding center, resulting in a protein with properties more closely resembling those of natural enzymes. Here we delineate the energetic and structural consequences associated with the introduction of these binding sites. To determine the extent to which the binding site was preorganized in the absence of metal ions, the apo structure of DF1 in solution was solved by NMR and compared with the crystal structure of the di-Zn(II) derivative. The overall fold of the apo protein was highly similar to that of the di-Zn(II) derivative, although there was a rotation of one of the helices. We also examined the thermodynamic consequences associated with building a small molecule-binding site within the protein. The protein exists in an equilibrium between folded dimers and unfolded monomers. DF1 is a highly stable protein (K(diss) = 0.001 fM), but the dissociation constant increases to 0.6 nM (deltadeltaG = 5.4 kcalmol monomer) as the active-site cavity is increased to accommodate small molecules.

Probing the configurational space of a metalloprotein core: an ab initio molecular dynamics study of Duo Ferro 1 binuclear Zn cofactor

We present three theoretical models of various degree of completeness to explore the chemical phase space available to the Glu4His2Zn2 cofactor found in the four-helix bundle of de novo designed metalloprotein Duo Ferro 1. We have found that the planewave DFT geometry optimization of 94-atom Model I, which contains both the protein scaffold constraints as well as the second shell hydrogen bonding network, reproduces the crystal structure bonding with remarkable accuracy (0.34 A). Surprisingly, the geometry optimization of 66-atom Model II (lacking the second shell hydrogen bonding) and 48-atom Model III (being also free of the protein scaffold constraints) still result in the fidelity with the crystallographic structure (RMSDs 0.29 and 0.34 A, respectively). To examine whether these structures are close to the global minimum as well as to investigate various conformational transitions to which the di-Zn cofactor may be susceptible to, we have carried out a 10 ps Car-Parrinello Molecular Dynamics (CPMD) simulation of Model III. We suggest that weak hydrogen bonds between imidazole hydrogens and carboxylate oxygens modulate the dynamical behavior of the system. One part of the molecule was found to be rigid due to the particular H(imidazole)-O(carboxylate) interaction restricting both the motion of the imidazole ring as well as the terminal carboxylate conformational mobility. The second half of the system was very flexible demonstrating a coupling of a transient formation of H(imidazole)-O(carboxylate) bonds with the spinning of the imidazole ring and syn-anti isomerization of the terminal carboxylate group. In addition, two low-energy snapshots from the 10 ps CPMD run were quenched, and their geometries were optimized, leading to two new isomers 48 kJ/mol lower in energy than the one associated with the crystal structure. We suggest that periodic quenching of the CPMD simulation snapshots of a minimalist model may be used as an efficient method to generate a large number of competitive local minima, which may be consequently pruned by imposing the protein scaffold constraints as well as further tuned by the second shell hydrogen bonding network.

Computational de novo design, and characterization of an A(2)B(2) diiron protein

Diiron proteins are found throughout nature and have a diverse range of functions; proteins in this class include methane monooxygenase, ribonucleotide reductase, Delta(9)-acyl carrier protein desaturase, rubrerythrin, hemerythrin, and the ferritins. Although each of these proteins has a very different overall fold, in every case the diiron active site is situated within a four-helix bundle. Additionally, nearly all of these proteins have a conserved Glu-Xxx-Xxx-His motif on two of the four helices with the Glu and His residues ligating the iron atoms. Intriguingly, subtle differences in the active site can result in a wide variety of functions. To probe the structural basis for this diversity, we designed an A(2)B(2) heterotetrameric four-helix bundle with an active site similar to those found in the naturally occurring diiron proteins. A novel computational approach was developed for the design, which considers the energy of not only the desired fold but also alternatively folded structures. Circular dichroism spectroscopy, analytical ultracentrifugation, and thermal unfolding studies indicate that the A and B peptides specifically associate to form an A(2)B(2) heterotetramer. Further, the protein binds Zn(II) and Co(II) in the expected manner and shows ferroxidase activity under single turnover conditions.

A hierarchic approach to the design of hexameric helical barrels

The design of large macromolecular assemblies is an endeavor with implications for protein engineering as well as nanotechnology. A hierarchic approach was used to design an antiparallel hexameric, tubular assembly of helices. In previous studies, a domain-swapped, dimeric three-helix bundle was designed from first principles. In the crystal lattice, three dimers associate around a 3-fold rotational axis to form a hexameric assembly. Although this hexameric assembly was not observed in solution, it was possible to stabilize its formation by changing three polar residues per monomer to hydrophobic (two Phe and one Trp) residues. Molecular models based on the crystallographic coordinates of DSD (PDB accession code 1G6U) show that these side-chains pack in the central cavity (the "supercore") of the hexameric bundle. Analytical ultracentrifugation, fluorescence spectroscopy, CD spectroscopy, and guanidine-HCl denaturation were used to determine the assembly of the hexamer. To probe the requirements for stabilizing the hexamer, we systematically varied the polarity and steric bulk of one of the Phe residues in the supercore of the hexamer. Depending on the nature of this side-chain, it is possible to modulate the stability of the hexamer in a predictable manner. This family of hexameric proteins may provide a useful framework for the construction of proteins that change their oligomeric states in response to binding of small molecules.

D(n)-symmetrical tertiary templates for the design of tubular proteins

Antiparallel helical bundles are found in a wide range of proteins. Often, four-helical bundles form tube-like structures, with binding sites for substrates or cofactors near their centers. For example, a transmembrane four-helical bundle in cytochrome bc(1) binds a pair of porphyrins in an elongated central cavity running down the center of the structure. Antiparallel helical barrels with larger diameters are found in the crystal structures of TolC and DSD, which form antiparallel 12-helical and six-helical bundles, respectively. The backbone geometries of the helical bundles of cytochrome bc(1), TolC, and DSD are well described using a simple D(n)-symmetrical model with only eight adjustable parameters. This parameterization provides an excellent starting point for construction of minimal models of these proteins as well as the de novo design of proteins with novel functions.

Crystal structure of manganese catalase from Lactobacillus plantarum

BACKGROUND

Catalases are important antioxidant metalloenzymes that catalyze disproportionation of hydrogen peroxide, forming dioxygen and water. Two families of catalases are known, one having a heme cofactor, and the other, a structurally distinct family containing nonheme manganese. We have solved the structure of the mesophilic manganese catalase from Lactobacillus plantarum and its azide-inhibited complex.

RESULTS

The crystal structure of the native enzyme has been solved at 1.8 A resolution by molecular replacement, and the azide complex of the native protein has been solved at 1.4 A resolution. The hexameric structure of the holoenzyme is stabilized by extensive intersubunit contacts, including a beta zipper and a structural calcium ion crosslinking neighboring subunits. Each subunit contains a dimanganese active site, accessed by a single substrate channel lined by charged residues. The manganese ions are linked by a mu1,3-bridging glutamate carboxylate and two mu-bridging solvent oxygens that electronically couple the metal centers. The active site region includes two residues (Arg147 and Glu178) that appear to be unique to the Lactobacillus plantarum catalase.

CONCLUSIONS

A comparison of L. plantarum and T. thermophilus catalase structures reveals the existence of two distinct structural classes, differing in monomer design and the organization of their active sites, within the manganese catalase family. These differences have important implications for catalysis and may reflect distinct biological functions for the two enzymes, with the L. plantarum enzyme serving as a catalase, while the T. thermophilus enzyme may function as a catalase/peroxidase.

De novo design of helical bundles as models for understanding protein folding and function

De novo protein design has proven to be a powerful tool for understanding protein folding, structure, and function. In this Account, we highlight aspects of our research on the design of dimeric, four-helix bundles. Dimeric, four-helix bundles are found throughout nature, and the history of their design in our laboratory illustrates our hierarchic approach to protein design. This approach has been successfully applied to create a completely native-like protein. Structural and mutational analysis allowed us to explore the determinants of native protein structure. These determinants were then applied to the design of a dinuclear metal-binding protein that can now serve as a model for this important class of proteins.

Miniaturized metalloproteins: application to iron-sulfur proteins

The miniaturization process applied to rubredoxins generated a class of peptide-based metalloprotein models, named METP (miniaturized electron transfer protein). The crystal structure of Desulfovibrio vulgaris rubredoxin was selected as a template for the construction of a tetrahedral (S(gamma)-Cys)(4) iron-binding site. Analysis of the structure showed that a sphere of 17 A in diameter, centered on the metal, circumscribes two unconnected approximately C(2) symmetry related beta-hairpins, each containing the -Cys-(Aaa)(2)-Cys- sequence. These observations provided a starting point for the design of an undecapeptide, which self assembles in the presence of tetrahedrally coordinating metal ions. The METP peptide was synthesized in good yield by standard methodologies. Successful assembly of the METP peptide with Co(II), Zn(II), Fe(II/III), in the expected 2:1 stoichiometry, was proven by UV-visible and circular dichroism spectroscopies. UV-visible analysis of the metal complexes indicated the four Cys ligands tetrahedrally arrange around the metal ion, as designed. Circular dichroism measurements on both the free and metal-bound forms revealed that the metal coordination drives the peptide chain to fold into a turned conformation. NMR characterization of the Zn(II)-METP complex fully supported the structure of the designed model. These results prove that METP reproduces the main features of rubredoxin.