Construction of a family of Cys2His2 zinc binding sites in the hydrophobic core of thioredoxin by structure-based design

Switchable Zinc(II)-Responsive Globular β-Sheet Peptide

The natural function of many proteins depends on their ability to switch their conformation driven by environmental changes. In this work, we present a small, monomeric β-sheet peptide that switches between a molten globule and a folded state through Zn(II) binding. The solvent-exposed hydrophobic core on the β-sheet surface was substituted by a His₃-site, whereas the internal hydrophobic core was left intact. Zn(II) is specifically recognized by the peptide relative to other divalent metal ions, binds in the lower micromolar range, and can be removed and re-added without denaturation of the peptide. In addition, the peptide is fully pH-switchable, has a pK_a of about 6, and survives several cycles of acidification and neutralization. In-depth structural characterization of the switch was achieved by concerted application of circular dichroism (CD) and multinuclear NMR spectroscopy. Thus, this study represents a viable approach toward a globular β-sheet Zn(II) mini-receptor prototype.

Designing hydrolytic zinc metalloenzymes

Zinc is an essential element required for the function of more than 300 enzymes spanning all classes. Despite years of dedicated study, questions regarding the connections between primary and secondary metal ligands and protein structure and function remain unanswered, despite numerous mechanistic, structural, biochemical, and synthetic model studies. Protein design is a powerful strategy for reproducing native metal sites that may be applied to answering some of these questions and subsequently generating novel zinc enzymes. From examination of the earliest design studies introducing simple Zn(II)-binding sites into de novo and natural protein scaffolds to current studies involving the preparation of efficient hydrolytic zinc sites, it is increasingly likely that protein design will achieve reaction rates previously thought possible only for native enzymes. This Current Topic will review the design and redesign of Zn(II)-binding sites in de novo-designed proteins and native protein scaffolds toward the preparation of catalytic hydrolytic sites. After discussing the preparation of Zn(II)-binding sites in various scaffolds, we will describe relevant examples for reengineering existing zinc sites to generate new or altered catalytic activities. Then, we will describe our work on the preparation of a de novo-designed hydrolytic zinc site in detail and present comparisons to related designed zinc sites. Collectively, these studies demonstrate the significant progress being made toward building zinc metalloenzymes from the bottom up.

Prediction of structures of zinc-binding proteins through explicit modeling of metal coordination geometry

Metal ions play an essential role in stabilizing protein structures and contributing to protein function. Ions such as zinc have well-defined coordination geometries, but it has not been easy to take advantage of this knowledge in protein structure prediction efforts. Here, we present a computational method to predict structures of zinc-binding proteins given knowledge of the positions of zinc-coordinating residues in the amino acid sequence. The method takes advantage of the "atom-tree" representation of molecular systems and modular architecture of the Rosetta3 software suite to incorporate explicit metal ion coordination geometry into previously developed de novo prediction and loop modeling protocols. Zinc cofactors are tethered to their interacting residues based on coordination geometries observed in natural zinc-binding proteins. The incorporation of explicit zinc atoms and their coordination geometry in both de novo structure prediction and loop modeling significantly improves sampling near the native conformation. The method can be readily extended to predict protein structures bound to other metal and/or small chemical cofactors with well-defined coordination or ligation geometry.

OptGraft: A computational procedure for transferring a binding site onto an existing protein scaffold

One of the many challenging tasks of protein design is the introduction of a completely new function into an existing protein scaffold. In this study, we introduce a new computational procedure OptGraft for placing a novel binding pocket onto a protein structure so as its geometry is minimally perturbed. This is accomplished by introducing a two-level procedure where we first identify where are the most appropriate locations to graft the new binding pocket into the protein fold by minimizing the departure from a set of geometric restraints using mixed-integer linear optimization. On identifying the suitable locations that can accommodate the new binding pocket, CHARMM energy calculations are employed to identify what mutations in the neighboring residues, if any, are needed to ensure that the minimum energy conformation of the binding pocket conserves the desired geometry. This computational framework is benchmarked against the results available in the literature for engineering a copper binding site into thioredoxin protein. Subsequently, OptGraft is used to guide the transfer of a calcium-binding pocket from thermitase protein (PDB: 1thm) into the first domain of CD2 protein (PDB:1hng). Experimental characterization of three de novo redesigned proteins with grafted calcium-binding centers demonstrated that they all exhibit high affinities for terbium (Kd) approximately 22, 38, and 55 microM) and can selectively bind calcium over magnesium.

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.

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.

Escherichia coli thioredoxin inhibition by cadmium: two mutually exclusive binding sites involving Cys32 and Asp26

Observations of thioredoxin inhibition by cadmium and of a positive role for thioredoxin in protection from Cd(2+) led us to investigate the thioredoxin-cadmium interaction properties. We used calorimetric and spectroscopic methods at different pH values to explore the relative contribution of putative binding residues (Cys32, Cys35, Trp28, Trp31 and Asp26) within or near the active site. At pH 8 or 7.5 two binding sites were identified by isothermal titration calorimetry with affinity constants of 10 x 10(6) m(-1) and 1 x 10(6) m(-1). For both sites, a proton was released upon Cd(2+) binding. One mole of Cd(2+) per mole of reduced thioredoxin was measured by mass spectrometry at these pH values, demonstrating that the two binding sites were partially occupied and mutually exclusive. Cd(2+) binding at either site totally inhibited the thiol-disulfide transferase activity of Trx. The absence of Cd(2+) interaction detected for oxidized or alkylated Trx and the inhibition of the enzymatic activity of thioredoxin by Cd(2+) supported the role of Cys32 at the first site. The fluorescence profile of Cd(2+)-bound thioredoxin differed, however, from that of oxidized thioredoxin, indicating that Cd(2+) was not coordinated with Cys32 and Cys35. From FTIR spectroscopy, we inferred that the second site might involve Asp26, a buried residue that deprotonates at a rather high and unusual pK(a) for a carboxylate (7.5/9.2). The pK(a) of the two residues Cys32 and Asp26 have been shown to be interdependent [Chivers, T. P. (1997) Biochemistry36, 14985-14991]. A mechanism is proposed in which Cd(2+) binding at the solvent-accessible thiolate group of Cys32 induces a decrease of the pK(a) of Asp26 and its deprotonation. Conversely, interaction between the carboxylate group of Asp26 and Cd(2+) at a second binding site induces Cys32 deprotonation and thioredoxin inhibition, so that Cd(2+) inhibits thioredoxin activity not only by binding at the Cys32 but also by interacting with Asp26.

Thioredoxin 2, an oxidative stress-induced protein, contains a high affinity zinc binding site

Two thioredoxins have been described in Escherichia coli, TrxA and Trx2. Both thioredoxins are capable of reducing disulfide bonds using a conserved pair of cysteine residues present in a WCGPC motif. A number of unique structural and regulatory features distinguish the Trx2 subfamily from the much larger TrxA family. The Trx2 subfamily has an additional N-terminal domain of +/- 30 residues, which contains two additional conserved CXXC motifs. Moreover, the gene coding for Trx2 is under control of the oxidative stress transcription factor OxyR in E. coli. This suggests that Trx2 may play a role in the cellular defense against oxidative stress. We show here that Trx2 contains zinc in a 1:1 stoichiometry, making it the first identified zinc-binding thioredoxin. The zinc atom is coordinated by the four cysteines of the two N-terminal CXXC motifs. The zinc center of Trx2 binds zinc with a very high affinity (Ka of >1018 m-1). We show that in vitro oxidation of the zinc binding cysteines by H2O2 releases the zinc and induces a conformational change. The zinc-free protein conserves its reductase activity. Altogether, our results suggest that the zinc center might play the role of a redox switch, changing a yet to be identified activity.

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.

NMR structure of the single QALGGH zinc finger domain from the Arabidopsis thaliana SUPERMAN protein

Zinc finger domains of the classical type represent the most abundant DNA binding domains in eukaryotic transcription factors. Plant proteins contain from one to four zinc finger domains, which are characterized by high conservation of the sequence QALGGH, shown to be critical for DNA-binding activity. The Arabidopsis thaliana SUPERMAN protein, which contains a single QALGGH zinc finger, is necessary for proper spatial development of reproductive floral tissues and has been shown to specifically bind to DNA. Here, we report the synthesis and UV and NMR spectroscopic structural characterization of a 37 amino acid SUPERMAN region complexed to a Zn(2+) ion (Zn-SUP37) and present the first high-resolution structure of a classical zinc finger domain from a plant protein. The NMR structure of the SUPERMAN zinc finger domain consists of a very well-defined betabetaalpha motif, typical of all other Cys(2)-His(2) zinc fingers structurally characterized. As a consequence, the highly conserved QALGGH sequence is located at the N terminus of the alpha helix. This region of the domain of animal zinc finger proteins consists of hypervariable residues that are responsible for recognizing the DNA bases. Therefore, we propose a peculiar DNA recognition code for the QALGGH zinc finger domain that includes all or some of the amino acid residues at positions -1, 2, and 3 (numbered relative to the N terminus of the helix) and possibly others at the C-terminal end of the recognition helix. This study further confirms that the zinc finger domain, though very simple, is an extremely versatile DNA binding motif.

Synthetic peptides mimicking the interleukin-6/gp130 interaction: a two-helix bundle system. Design and conformational studies

The objective of our study was to mimic in a typical reductionist approach the molecular interactions observed at the interface between the gp130 receptor and interleukin-6 during formation of their complex. A peptide system obtained by reproducing some of the interleukin-6/gp130 molecular interactions into a two-helix bundle structure was investigated. The solution conformational features of this system were determined by CD and NMR techniques. The CD titration experiments demonstrated that the interaction between the designed peptides is specific and based on a well-defined stoichiometry. The NMR data confirmed some of the structural features of the binding mechanism as predicted by the rational design and indicated that under our conditions the recognition specificity and affinity can be explained by the formation of a two-helix bundle. Thus, the data reported herein represent a promising indication on how to develop new peptides able to interfere with formation of the interleukin-6/gp130 complex.

Metal-dependent stabilization of an active HMG protein

Using a cloned single domain of the high mobility group protein 1 (HMGB1), we evaluated the effect of introducing metal binding site(s) on protein stability and function. An HMG domain is a conserved sequence of approximately 80 amino acids rich in basic, aromatic and proline residues that is active in binding DNA in a sequence- or structure-specific manner. The design strategy focuses on anchoring selected regions of the protein, specifically loops and turns in the molecule, using His-metal ligands. Changes in secondary structure, thermostability and DNA binding properties of a series of such mutants were evaluated. The two most stable mutant constructs contain three surface histidine replacements (two metal binding sites) in the regions encompassing both turns of the molecule. On ligation with the divalent nickel cation, the stability of these two triple histidine mutants (I38H/N51H/D55H and G39H/N51H/D55H) increases by 1.3 and 1.6 kcal/mol, respectively, relative to the wild-type protein, although the creation of binding sites per se destabilizes the protein. The DNA-binding properties of the modified proteins are not impaired by the introduction of the metal binding motifs. These results indicate that it is feasible to stabilize protein tertiary structure using appropriate placement of surface His-metal bonds without loss of function.

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 new zinc binding fold underlines the versatility of zinc binding modules in protein evolution

Many different zinc binding modules have been identified. Their abundance and variety suggests that the formation of zinc binding folds might be relatively common. We have determined the structure of CH1(1), a 27-residue peptide derived from the first cysteine/histidine-rich region (CH1) of CREB binding protein (CBP). This peptide forms a highly ordered zinc-dependent fold that is distinct from known folds. The structure differs from a subsequently determined structure of a larger region from the CH3 region of CBP, and the CH1(1) fold probably represents a nonphysiologically active form. Despite this, the fold is thermostable and tolerant to both multiple alanine mutations and changes in the zinc-ligand spacing. Our data support the idea that zinc binding domains may arise frequently. Additionally, such structures may prove useful as scaffolds for protein design, given their stability and robustness.

The metal site as a template for the metalloprotein structure formation

Achieving a thorough explanation of the behavior of metal sites in the formation of native metalloprotein structures is an exciting challenge in the biochemistry of metallobiomacromolecules. This study presents a personal insight into the subject. It is proposed that a metal center and its exogenous ligand compose a template. A template may impose a clear stereochemical preference on the loose peptide chains, and organize them into natural stereospecificity via the metal-ligand interaction, a long-range and strong interaction. Therefore, the stable peptide conformation induced by the template effect surrounding a template polyhedron could be called a template-mediated structural motif (TMSM).

Surface grafting onto template-assembled synthetic protein scaffolds in molecular recognition

Creating functional biological molecules de novo requires a detailed understanding of the intimate relationship between primary sequence, folding mechanism, and packing topology, and remains up to now a most challenging goal in protein design and mimicry. As a consequence, the use of well-defined robust macromolecules as scaffolds for the introduction of function by grafting surface residues has become a major objective in protein engineering and de novo design. In this article, the concept of scaffolds is demonstrated on some selected examples, illustrating that novel types of functional molecules can be generated. Reengineered proteins and, most notably, de novo designed peptide scaffolds exhibiting molecular function, are ideal tools for structure-function studies and as leads in drug design.

Conversion of a maltose receptor into a zinc biosensor by computational design

We have demonstrated that it is possible to radically change the specificity of maltose binding protein by converting it into a zinc sensor using a rational design approach. In this new molecular sensor, zinc binding is transduced into a readily detected fluorescence signal by use of an engineered conformational coupling mechanism linking ligand binding to reporter group response. An iterative progressive design strategy led to the construction of variants with increased zinc affinity by combining binding sites, optimizing the primary coordination sphere, and exploiting conformational equilibria. Intermediates in the design series show that the adaptive process involves both introduction and optimization of new functions and removal of adverse vestigial interactions. The latter demonstrates the importance of the rational design approach in uncovering cryptic phenomena in protein function, which cannot be revealed by the study of naturally evolved systems.

Nanohedra: using symmetry to design self assembling protein cages, layers, crystals, and filaments

A general strategy is described for designing proteins that self assemble into large symmetrical nanomaterials, including molecular cages, filaments, layers, and porous materials. In this strategy, one molecule of protein A, which naturally forms a self-assembling oligomer, A(n), is fused rigidly to one molecule of protein B, which forms another self-assembling oligomer, B(m). The result is a fusion protein, A-B, which self assembles with other identical copies of itself into a designed nanohedral particle or material, (A-B)(p). The strategy is demonstrated through the design, production, and characterization of two fusion proteins: a 49-kDa protein designed to assemble into a cage approximately 15 nm across, and a 44-kDa protein designed to assemble into long filaments approximately 4 nm wide. The strategy opens a way to create a wide variety of potentially useful protein-based materials, some of which share similar features with natural biological assemblies.

De novo design and structural characterization of proteins and metalloproteins

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

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

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

Rational design of nascent metalloenzymes

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

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

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

Hydrophobic core packing and protein design

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

The development of new biotechnologies using metalloprotein design

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