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

Designing Artificial Metalloenzymes by Tuning of the Environment beyond the Primary Coordination Sphere

Metalloenzymes catalyze a variety of reactions using a limited number of natural amino acids and metallocofactors. Therefore, the environment beyond the primary coordination sphere must play an important role in both conferring and tuning their phenomenal catalytic properties, enabling active sites with otherwise similar primary coordination environments to perform a diverse array of biological functions. However, since the interactions beyond the primary coordination sphere are numerous and weak, it has been difficult to pinpoint structural features responsible for the tuning of activities of native enzymes. Designing artificial metalloenzymes (ArMs) offers an excellent basis to elucidate the roles of these interactions and to further develop practical biological catalysts. In this review, we highlight how the secondary coordination spheres of ArMs influence metal binding and catalysis, with particular focus on the use of native protein scaffolds as templates for the design of ArMs by either rational design aided by computational modeling, directed evolution, or a combination of both approaches. In describing successes in designing heme, nonheme Fe, and Cu metalloenzymes, heteronuclear metalloenzymes containing heme, and those ArMs containing other metal centers (including those with non-native metal ions and metallocofactors), we have summarized insights gained on how careful controls of the interactions in the secondary coordination sphere, including hydrophobic and hydrogen bonding interactions, allow the generation and tuning of these respective systems to approach, rival, and, in a few cases, exceed those of native enzymes. We have also provided an outlook on the remaining challenges in the field and future directions that will allow for a deeper understanding of the secondary coordination sphere a deeper understanding of the secondary coordintion sphere to be gained, and in turn to guide the design of a broader and more efficient variety of ArMs.

Rational De Novo Design of a Cu Metalloenzyme for Superoxide Dismutation

Superoxide dismutases (SODs) are highly efficient enzymes for superoxide dismutation and the first line of defense against oxidative stress. These metalloproteins contain a redox-active metal ion in their active site (Mn, Cu, Fe, Ni) with a tightly controlled reduction potential found in a close range around the optimal value of 0.36 V versus the normal hydrogen electrode (NHE). Rationally designed proteins with well-defined three-dimensional structures offer new opportunities for obtaining functional SOD mimics. Here, we explore four different copper-binding scaffolds: H3 (His3 ), H4 (His4 ), H2 DH (His3 Asp with two His and one Asp in the same plane) and H3 D (His3 Asp with three His in the same plane) by using the scaffold of the de novo protein GRα3 D. EPR and XAS analysis of the resulting copper complexes demonstrates that they are good CuII -bound structural mimics of Cu-only SODs. Furthermore, all the complexes exhibit SOD activity, though three orders of magnitude slower than the native enzyme, making them the first de novo copper SOD mimics.

Design of catalytic polypeptides and proteins

The first part of this review article lists examples of complete, empirical de novo design that made important contributions to the development of the field and initiated challenging projects. The second part of this article deals with computational design of novel enzymes in native protein scaffolds; active designs were refined through random and site-directed mutagenesis producing artificial enzymes with nearly native enzyme- like activities against a number of non-natural substrates. Combining aspects of de novo design and biological evolution of nature's enzymes has started and will accelerate the development of novel enzyme activities.

A triple-amplification SPR electrochemiluminescence assay for chloramphenicol based on polymer enzyme-linked nanotracers and exonuclease-assisted target recycling

The present study aimed to explore a novel triple-amplification electrochemiluminescence (ECL) assay for detecting of chloramphenicol (CAP). This strategy was based on single-stranded DNA-binding protein (SSB) and horseradish peroxidase (HRP) enzyme-linked polymer (EnVision reagent, EV) labeled on Au nanoparticles (EV-Au-SSB) as nanotracer and exonuclease-assisted target recycling. The composite probes were prepared via immunoreactions between the CdS nanocrystal (CdS NC)-functionalized partial complementary DNA and aptamer (CdSNCs/Apt-ssDNA1) as capture probes, and EV-Au-SSB as nanotracer. When the composite probe solution co-existed with CAP and Exo I, the aptamer on the capture probes preferentially combined with CAP, and then CAP-Apt and nanotracer complex were released into the solution. Subsequently, Exo I in the solution could further digest the CAP-Apt from the 3'-end of the aptamer and release CAP, which could participate in further reaction with the probes. It was worth mentioning that EV contained a large number of HRPs on its dendritic chain. In the EV-Au-SSB, Au could enhance ECL intensity of CdS NCs by surface plasmon resonance. What's more, HRPs on EV could catalyze the reaction of H2O2, which could obviously enhance ECL intensity of CdS NCs. This study demonstrated excellent performance of the triple-amplification ECL assay, which makes this aptasensor system suitable and promising for the practical application of CAP residues in fish samples. Moreover, the assay might provide a promising avenue to develop efficient aptasensors to determine small-molecule harmful substances in environmental monitoring and food safety.

Various strategies for obtaining oxidative artificial hemoproteins with a catalytic oxidative activity: from "Hemoabzymes" to "Hemozymes"?

The design of artificial hemoproteins that could lead to new biocatalysts for selective oxidation reactions using clean oxidants such as O 2 or H 2 O 2 under ecocompatible conditions constitutes a really promising challenge for a wide range of industrial applications. In vivo, such reactions are performed by heme-thiolate proteins, cytochromes P450, that catalyze the oxidation of drugs by dioxygen in the presence of electrons delivered from NADPH by cytochrome P450 reductase. Several strategies were used to design new artificial hemoproteins to mimic these enzymes, that associate synthetic metalloporphyrin derivatives to a protein that is supposed to induce a selectivity in the catalyzed reaction. A first generation of artificial hemoproteins or "hemoabzymes" was obtained by the non-covalent association of synthetic hemes such as N-methyl-mesoporphyrin IX, Fe(III) -α3β-tetra-o-carboxyphenylporphyrin or microperoxidase 8 with monoclonal antibodies raised against these cofactors. The obtained antibody-metalloporphyrin complexes displayed a peroxidase activity and some of them catalyzed the regio-selective nitration of phenols by H 2 O 2/ NO 2 and the stereo-selective oxidation of sulphides by H 2 O 2. A second generation of artificial hemoproteins or "hemozymes", was obtained by the non-covalent association of non-relevant proteins with metalloporphyrin derivatives. Several strategies were used, the most successful of which, named "host-guest" strategy involved the non-covalent incorporation of metalloporphyrin derivatives into easily affordable proteins. The artificial hemoproteins obtained were found to be able to perform efficiently the stereoselective oxidation of organic compounds such as sulphides and alkenes by H 2 O 2 and KHSO 5.

De novo design of stereochemically-bent sixteen-residue β-hairpin as a hydrolase mimic

Stepwise design of sixteen-residue β-hairpin as a hydrolase mimic involving fold design by stereochemical mutation followed by inverse-design of sequence.

From “hemoabzymes” to “hemozymes”: towards new biocatalysts for selective oxidations

Two generations of artificial hemoproteins have been obtained: “hemoabzymes”, by non-covalent association of synthetic hemes with monoclonal antibodies raised against these cofactors and “hemozymes”, by non-covalent association of non-relevant proteins with metalloporphyrin derivatives. A review of the different strategies employed as well as their structural and catalytic properties is presented here.

Redesigning the blue copper azurin into a redox-active mononuclear nonheme iron protein: preparation and study of Fe(II)-M121E azurin

Much progress has been made in designing heme and dinuclear nonheme iron enzymes. In contrast, engineering mononuclear nonheme iron enzymes is lagging, even though these enzymes belong to a large class that catalyzes quite diverse reactions. Herein we report spectroscopic and X-ray crystallographic studies of Fe(II)-M121E azurin (Az), by replacing the axial Met121 and Cu(II) in wild-type azurin (wtAz) with Glu and Fe(II), respectively. In contrast to the redox inactive Fe(II)-wtAz, the Fe(II)-M121EAz mutant can be readily oxidized by Na2IrCl6, and interestingly, the protein exhibits superoxide scavenging activity. Mössbauer and EPR spectroscopies, along with X-ray structural comparisons, revealed similarities and differences between Fe(II)-M121EAz, Fe(II)-wtAz, and superoxide reductase (SOR) and allowed design of the second generation mutant, Fe(II)-M121EM44KAz, that exhibits increased superoxide scavenging activity by 2 orders of magnitude. This finding demonstrates the importance of noncovalent secondary coordination sphere interactions in fine-tuning enzymatic activity.

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.

Identification of protein scaffolds for enzyme design using scaffold selection

The identification of suitable protein structures that can serve as scaffolds for the introduction of catalytic residues is crucial for the design of new enzymes. Here we describe how the automated and rapid scaffold search program ScaffoldSelection can be used to find the best starting points, namely protein structures that are most likely to tolerate the introduction and promote the proper formation of a specific catalytic motif.

Design of protein catalysts

Diverse engineering strategies have been developed to create enzymes with novel catalytic activities. Among these, computational approaches hold particular promise. Enzymes have been computationally designed to promote several nonbiological reactions, including a Diels-Alder cycloaddition, proton transfer, multistep retroaldol transformations, and metal-dependent hydrolysis of phosphotriesters. Although their efficiencies (kcat/KM = 0.1-100 M(-1) s(-1)) are typically low compared with those of the best natural enzymes (10(6)-10(8) M(-1) s(-1)), these catalysts are excellent starting points for laboratory evolution. This review surveys recent progress in combining computational and evolutionary approaches to enzyme design, together with insights into enzyme function gained from studies of the engineered catalysts.

Systematic optimization model and algorithm for binding sequence selection in computational enzyme design

A systematic optimization model for binding sequence selection in computational enzyme design was developed based on the transition state theory of enzyme catalysis and graph-theoretical modeling. The saddle point on the free energy surface of the reaction system was represented by catalytic geometrical constraints, and the binding energy between the active site and transition state was minimized to reduce the activation energy barrier. The resulting hyperscale combinatorial optimization problem was tackled using a novel heuristic global optimization algorithm, which was inspired and tested by the protein core sequence selection problem. The sequence recapitulation tests on native active sites for two enzyme catalyzed hydrolytic reactions were applied to evaluate the predictive power of the design methodology. The results of the calculation show that most of the native binding sites can be successfully identified if the catalytic geometrical constraints and the structural motifs of the substrate are taken into account. Reliably predicting active site sequences may have significant implications for the creation of novel enzymes that are capable of catalyzing targeted chemical reactions.

Energy functions in de novo protein design: current challenges and future prospects

In the past decade, a concerted effort to successfully capture specific tertiary packing interactions produced specific three-dimensional structures for many de novo designed proteins that are validated by nuclear magnetic resonance and/or X-ray crystallographic techniques. However, the success rate of computational design remains low. In this review, we provide an overview of experimentally validated, de novo designed proteins and compare four available programs, RosettaDesign, EGAD, Liang-Grishin, and RosettaDesign-SR, by assessing designed sequences computationally. Computational assessment includes the recovery of native sequences, the calculation of sizes of hydrophobic patches and total solvent-accessible surface area, and the prediction of structural properties such as intrinsic disorder, secondary structures, and three-dimensional structures. This computational assessment, together with a recent community-wide experiment in assessing scoring functions for interface design, suggests that the next-generation protein-design scoring function will come from the right balance of complementary interaction terms. Such balance may be found when more negative experimental data become available as part of a training set.

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.

A matching algorithm for catalytic residue site selection in computational enzyme design

A loop closure-based sequential algorithm, PRODA_MATCH, was developed to match catalytic residues onto a scaffold for enzyme design in silico. The computational complexity of this algorithm is polynomial with respect to the number of active sites, the number of catalytic residues, and the maximal iteration number of cyclic coordinate descent steps. This matching algorithm is independent of a rotamer library that enables the catalytic residue to take any required conformation during the reaction coordinate. The catalytic geometric parameters defined between functional groups of transition state (TS) and the catalytic residues are continuously optimized to identify the accurate position of the TS. Pseudo-spheres are introduced for surrounding residues, which make the algorithm take binding into account as early as during the matching process. Recapitulation of native catalytic residue sites was used as a benchmark to evaluate the novel algorithm. The calculation results for the test set show that the native catalytic residue sites were successfully identified and ranked within the top 10 designs for 7 of the 10 chemical reactions. This indicates that the matching algorithm has the potential to be used for designing industrial enzymes for desired reactions.

Improving computational protein design by using structure-derived sequence profile

Designing a protein sequence that will fold into a predefined structure is of both practical and fundamental interest. Many successful, computational designs in the last decade resulted from improved understanding of hydrophobic and polar interactions between side chains of amino acid residues in stabilizing protein tertiary structures. However, the coupling between main-chain backbone structure and local sequence has yet to be fully addressed. Here, we attempt to account for such coupling by using a sequence profile derived from the sequences of five residue fragments in a fragment library that are structurally matched to the five-residue segments contained in a target structure. We further introduced a term to reduce low complexity regions of designed sequences. These two terms together with optimized reference states for amino-acid residues were implemented in the RosettaDesign program. The new method, called RosettaDesign-SR, makes a 12% increase (from 34 to 46%) in fraction of proteins whose designed sequences are more than 35% identical to wild-type sequences. Meanwhile, it reduces 8% (from 22% to 14%) to the number of designed sequences that are not homologous to any known protein sequences according to psi-blast. More importantly, the sequences designed by RosettaDesign-SR have 2-3% more polar residues at the surface and core regions of proteins and these surface and core polar residues have about 4% higher sequence identity to wild-type sequences than by RosettaDesign. Thus, the proteins designed by RosettaDesign-SR should be less likely to aggregate and more likely to have unique structures due to more specific polar interactions.

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.

Motif-directed flexible backbone design of functional interactions

Computational protein design relies on a number of approximations to efficiently search the huge sequence space available to proteins. The fixed backbone and rotamer approximations in particular are important for formulating protein design as a discrete combinatorial optimization problem. However, the resulting coarse-grained sampling of possible side-chain terminal positions is problematic for the design of protein function, which depends on precise positioning of side-chain atoms. Although backbone flexibility can greatly increase the conformation freedom of side-chain functional groups, it is not obvious which backbone movements will generate the critical constellation of atoms responsible for protein function. Here, we report an automated method for identifying protein backbone movements that can give rise to any specified set of desired side-chain atomic placements and interactions, using protein-DNA interfaces as a model system. We use a library of previously observed protein-DNA interactions (motifs) and a rotamer-based description of side-chain conformation freedom to identify placements for the protein backbone that can give rise to a favorable side-chain interaction with DNA. We describe a tree-search algorithm for identifying those combinations of interactions from the library that can be realized with minimal perturbation of the protein backbone. We compare the efficiency of this method with the alternative approach of building and screening alternate backbone conformations.

Catalytic mechanism and performance of computationally designed enzymes for Kemp elimination

A series of enzymes for Kemp elimination of 5-nitrobenzisoxazole has been recently designed and tested. In conjunction with the design process, extensive computational analyses were carried out to evaluate the potential performance of four of the designs, as presented here. The enzyme-catalyzed reactions were modeled using mixed quantum and molecular mechanics (QM/MM) calculations in the context of Monte Carlo (MC) statistical mechanics simulations. Free-energy perturbation (FEP) calculations were used to characterize the free-energy surfaces for the catalyzed reactions as well as for reference processes in water. The simulations yielded detailed information about the catalytic mechanisms, activation barriers, and structural evolution of the active sites over the course of the reactions. The catalytic mechanism for the designed enzymes KE07, KE10(V131N), and KE15 was found to be concerted with proton transfer, generally more advanced in the transition state than breaking of the isoxazolyl N-O bond. On the basis of the free-energy results, all three enzymes were anticipated to be active. Ideas for further improvement of the enzyme designs also emerged. On the technical side, the synergy of parallel QM/MM and experimental efforts in the design of artificial enzymes is well illustrated.

Origin of the activity drop with the E50D variant of catalytic antibody 34E4 for Kemp elimination

In enzymes, multiple structural effects cooperatively lead to the high catalytic activity, while individually these effects can be small. The design of artificial enzymes requires the understanding and ability to manipulate such subtle effects. The 34E4 catalytic antibody, catalyzing Kemp elimination of 5-nitrobenzisoxazole, and its Glu50Asp (E50D) variant are the subject of the present investigation. This removal of only a methylene group yields an approximately 30-fold reduction in the rate for the catalyzed Kemp elimination. Here, the aim is to understand this difference in the catalytic performance. The mechanism of Kemp elimination catalyzed by 34E4 and the E50D mutant is elucidated using QM/MM Monte Carlo simulations and free energy perturbation theory. In both proteins, the reaction is shown to follow a single-step, concerted mechanism. In the mutant, the activation barrier rises by 2.4 kcal/mol, which corresponds to a 62-fold rate deceleration, which is in good agreement with the experimental data. The positions and functionality of the residues in the active site are monitored throughout the reaction. It is concluded that the looser contact with the base, shorter base-Asn58 contact, less favorable pi-stacking with Trp91 in the transition state of the reaction, and different solvation pattern all contribute to the reduction of the reaction rate in the E50D variant of 34E4.

TransCent: computational enzyme design by transferring active sites and considering constraints relevant for catalysis

Background

Computational enzyme design is far from being applicable for the general case. Due to computational complexity and limited knowledge of the structure-function interplay, heuristic methods have to be used.

Results

We have developed TransCent, a computational enzyme design method supporting the transfer of active sites from one enzyme to an alternative scaffold. In an optimization process, it balances requirements originating from four constraints. These are 1) protein stability, 2) ligand binding, 3) pK_a values of active site residues, and 4) structural features of the active site. Each constraint is handled by an individual software module. Modules processing the first three constraints are based on state-of-the-art concepts, i.e. RosettaDesign, DrugScore, and PROPKA. To account for the fourth constraint, knowledge-based potentials are utilized. The contribution of modules to the performance of TransCent was evaluated by means of a recapitulation test. The redesign of oxidoreductase cytochrome P450 was analyzed in detail. As a first application, we present and discuss models for the transfer of active sites in enzymes sharing the frequently encountered triosephosphate isomerase fold.

Conclusion

A recapitulation test on native enzymes showed that TransCent proposes active sites that resemble the native enzyme more than those generated by RosettaDesign alone. Additional tests demonstrated that each module contributes to the overall performance in a statistically significant manner.

Rational design of a novel calcium-binding site adjacent to the ligand-binding site on CD2 increases its CD48 affinity

Electrostatic interactions are important for molecular recognition processes including Ca2+-binding and cell adhesion. To understand these processes, we have successfully introduced a novel Ca2+-binding site into the non-Ca2+-dependent cell adhesion protein CD2 using our criteria that are specifically tailored to the structural and functional properties of the protein environment and charged adhesion surface. This designed site with ligand residues exclusively from the beta-sheets selectively binds to Ca2+ and Ln3+ over other mono- and divalent cations. While Ca2+ and Ln3+ binding specifically alters the local environment of the designed Ca2+-binding site, the designed protein undergoes a significantly smaller conformation change compared with those observed in naturally occurring Ca2+-binding sites that are composed of at least part of the flexible loop and helical regions. In addition, the CD2-CD48-binding affinity increased approximately threefold after protein engineering, suggesting that the cell adhesion of CD2 can be modulated by altering the local electrostatic environment. The study provides site-specific information for regulating cell adhesion within CD2 and gives insight into the structural factors required for Ca2+-modulated biological processes.

Minimalist active-site redesign: teaching old enzymes new tricks

Although nature evolves its catalysts over millions of years, enzyme engineers try to do it a bit faster. Enzyme active sites provide highly optimized microenvironments for the catalysis of biologically useful chemical transformations. Consequently, changes at these centers can have large effects on enzyme activity. The prediction and control of these effects provides a promising way to access new functions. The development of methods and strategies to explore the untapped catalytic potential of natural enzyme scaffolds has been pushed by the increasing demand for industrial biocatalysts. This Review describes the use of minimal modifications at enzyme active sites to expand their catalytic repertoires, including targeted mutagenesis and the addition of new reactive functionalities. Often, a novel activity can be obtained with only a single point mutation. The many successful examples of active-site engineering through minimal mutations give useful insights into enzyme evolution and open new avenues in biocatalyst research.

Novel enzymes through design and evolution

The generation of enzymes with new catalytic activities remains a major challenge. So far, several different strategies have been developed to tackle this problem, including site-directed mutagenesis, random mutagenesis (directed evolution), antibody catalysis, computational redesign, and de novo methods. Using these techniques, a broad array of novel enzymes has been created (aldolases, decarboxylases, dehydratases, isomerases, oxidases, reductases, and others), although their low efficiencies (10 to 100 M(-1) s(-l)) compared to those of the best natural enzymes (10(6) to 10(8) M(-1) s(-1)) remains a significant concern. Whereas rational design might be the most promising and versatile approach to generating new activities, directed evolution seems to be the best way to optimize the catalytic properties of novel enzymes. Indeed, impressive successes in enzyme engineering have resulted from a combination of rational and random design.

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.

Protein fabrication automation

Facile "writing" of DNA fragments that encode entire gene sequences potentially has widespread applications in biological analysis and engineering. Rapid writing of open reading frames (ORFs) for expressed proteins could transform protein engineering and production for protein design, synthetic biology, and structural analysis. Here we present a process, protein fabrication automation (PFA), which facilitates the rapid de novo construction of any desired ORF from oligonucleotides with low effort, high speed, and little human interaction. PFA comprises software for sequence design, data management, and the generation of instruction sets for liquid-handling robotics, a liquid-handling robot, a robust PCR scheme for gene assembly from synthetic oligonucleotides, and a genetic selection system to enrich correctly assembled full-length synthetic ORFs. The process is robust and scalable.

Local encoding of computationally designed enzyme activity

One aim of computational protein design is to introduce novel enzyme activity into proteins of known structure by predicting mutations that stabilize transition states. Previously, we showed that it is possible to introduce triose phosphate isomerase activity into the ribose-binding protein of Escherichia coli by constructing 17 mutations in the first two layers of residues that surround the wild-type ligand-binding site. Here, we report that these mutations can be "transplanted" into a homologous ribose-binding protein, isolated from the hyperthermophilic bacterium Thermoanaerobacter tengcongensis, with retention of catalytic activity, substrate affinity, and reaction pH dependence. The observed 10(5)-10(6)-fold rate enhancement corresponds to 70% of the maximally known transition-state binding energy. The wild-type sequences in these two homologues are almost perfectly conserved in the vicinity of their ribose-binding sites, but diverge significantly at increasing distance from these sites. The results demonstrate that the computationally designed mutations are sufficient to encode the observed enzyme activity, that all the observed activity is encoded locally within the layer of residues directly in contact with the substrate and that, in this case, at least 70% of transition state stabilization energy can be achieved using straightforward considerations of stereochemical complementarity between enzyme and reactants.

Of Folding and Function: Understanding Active-Site Context through Metalloenzyme Design

In the emerging field of biomolecular design, the introduction of metal-binding sites into loop or turn regions of known protein scaffolds has been utilized to create unique metalloprotein and metallopeptide systems for study. This Forum Article highlights examples of the modular-turn-substitution approach to design and the range of structural and mechanistic questions to which this tool can be applied. Examples from the authors' laboratory are given to show that lanthanide-binding metallopeptides, and now a full metallohomeodomain, can be generated by modular substitution of a Ca-binding EF-hand loop into the unrelated scaffold, the engrailed helix-turn-helix motif. We have previously shown that these peptides bind trivalent Ln(III) ions and promote DNA and phosphate hydrolysis, the targeted function. Here, a series of chimeric peptides are presented that differ only in the ninth loop position [given in parentheses; Peptides P3N (Asn), P3E (Glu), P3A (Ala), and P3W(D) (Asp]. This residue, a putative second-shell ligand stabilizing a coordinated water, was found to influence not only metal affinity but also peptide folding. The affinity for Tb(III) was determined by Trp-Tb fluorescence resonance energy transfer and followed the order Ka = P3W(D) > P3A approximately P3E > P3N. However, circular dichroism (CD) titrations with EuCl3 showed that only P3W(D) and P3N folded to any extent upon metal binding, indicating that the Asp/Asn side chains stabilize the central loop structure and thus propagate folding of the peripheral helices, whereas neither Ala nor Glu appears to be interacting with the metal to organize the loop. Finally, we investigated the longer range context of a given loop substitution by cloning and expressing a lanthanide-binding homeodomain (C2), whose loop insertion sequence is analogous to that of peptide P3W(D). We find by CD that apo-C2 has a significant helical structure (approximately 25% alphahelicity), which increases further upon the addition of Tb(III) (approximately 32% alpha helicity). The protein's Tb(III) affinity is similar to that of the chimeric peptides. However, unlike previously reported metallopeptides, we find that EuC2 does not appreciably promote phosphate or DNA cleavage, which suggests a difference in metal accessibility in the context of the full domain. We have demonstrated that substituting beta turns with metal-binding turns does not necessarily require homologous parental scaffolds or small flexible peptides but rather relies on the structural similarity of the motifs flanking the turn.

Solution structure of HndAc: a thioredoxin-like domain involved in the NADP-reducing hydrogenase complex

The NADP-reducing hydrogenase complex from Desulfovibrio fructosovorans is a heterotetramer encoded by the hndABCD operon. Sequence analysis indicates that the HndC subunit (52 kDa) corresponds to the NADP-reducing unit, and the HndD subunit (63.5 kDa) is homologous to Clostridium pasteurianum hydrogenase. The role of HndA and HndB subunits (18.8 kDa and 13.8 kDa, respectively) in the complex remains unknown. The HndA subunit belongs to the [2Fe-2S] ferredoxin family typified by C. pasteurianum ferredoxin. HndA is organized into two independent structural domains, and we report in the present work the NMR structure of its C-terminal domain, HndAc. HndAc has a thioredoxin-like fold consisting in four beta-strands and two relatively long helices. The [2Fe-2S] cluster is located near the surface of the protein and bound to four cysteine residues particularly well conserved in this class of proteins. Electron exchange between the HndD N-terminal [2Fe-2S] domain (HndDN) and HndAc has been previously evidenced, and in the present studies we have mapped the binding site of the HndDN domain on HndAc. A structural analysis of HndB indicates that it is a FeS subunit with 41% similarity with HndAc and it contains a possible thioredoxin-like fold. Our data let us propose that HndAc and HndB can form a heterodimeric intermediate in the electron transfer between the hydrogenase (HndD) active site and the NADP reduction site in HndC.

De novo proteins from designed combinatorial libraries

Combinatorial libraries of de novo amino acid sequences can provide a rich source of diversity for the discovery of novel proteins with interesting and important activities. Randomly generated sequences, however, rarely fold into well-ordered proteinlike structures. To enhance the quality of a library, features of rational design must be used to focus sequence diversity into those regions of sequence space that are most likely to yield folded structures. This review describes how focused libraries can be constructed by designing the binary pattern of polar and nonpolar amino acids to favor proteins that contain abundant secondary structure, while simultaneously burying hydrophobic side chains and exposing hydrophilic side chains to solvent. The "binary code" for protein design was used to construct several libraries of de novo proteins, including both alpha-helical and beta-sheet structures. The recently determined solution structure of a binary patterned four-helix bundle is well ordered, thereby demonstrating that sequences that have neither been selected by evolution (in vivo or in vitro) nor designed by computer can form nativelike proteins. Examples are presented demonstrating how binary patterned libraries have successfully produced well-ordered structures, cofactor binding, catalytic activity, self-assembled monolayers, amyloid-like nanofibrils, and protein-based biomaterials.

Computational design of a Zn2+ receptor that controls bacterial gene expression

The control of cellular physiology and gene expression in response to extracellular signals is a basic property of living systems. We have constructed a synthetic bacterial signal transduction pathway in which gene expression is controlled by extracellular Zn2+. In this system a computationally designed Zn2+-binding periplasmic receptor senses the extracellular solute and triggers a two-component signal transduction pathway via a chimeric transmembrane protein, resulting in transcriptional up-regulation of a beta-galactosidase reporter gene. The Zn2+-binding site in the designed receptor is based on a four-coordinate, tetrahedral primary coordination sphere consisting of histidines and glutamates. In addition, mutations were introduced in a secondary coordination sphere to satisfy the residual hydrogen-bonding potential of the histidines coordinated to the metal. The importance of the secondary shell interactions is demonstrated by their effect on metal affinity and selectivity, as well as protein stability. Three designed protein sequences, comprising two distinct metal-binding positions, were all shown to bind Zn2+ and to function in the cell-based assay, indicating the generality of the design methodology. These experiments demonstrate that biological systems can be manipulated with computationally designed proteins that have drastically altered ligand-binding specificities, thereby extending the repertoire of genetic control by extracellular signals.

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.