Characteristics of a de novo designed protein

Evolution, folding, and design of TIM barrels and related proteins

Proteins are chief actors in life that perform a myriad of exquisite functions. This diversity has been enabled through the evolution and diversification of protein folds. Analysis of sequences and structures strongly suggest that numerous protein pieces have been reused as building blocks and propagated to many modern folds. This information can be traced to understand how the protein world has diversified. In this review, we discuss the latest advances in the analysis of protein evolutionary units, and we use as a model system one of the most abundant and versatile topologies, the TIM-barrel fold, to highlight the existing common principles that interconnect protein evolution, structure, folding, function, and design.

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 a four-fold symmetric TIM-barrel protein with atomic-level accuracy

Despite efforts for over 25 years, de novo protein design has not succeeded in achieving the TIM-barrel fold. Here we describe the computational design of four-fold symmetrical (β/α)8 barrels guided by geometrical and chemical principles. Experimental characterization of 33 designs revealed the importance of side chain-backbone hydrogen bonds for defining the strand register between repeat units. The X-ray crystal structure of a designed thermostable 184-residue protein is nearly identical to that of the designed TIM-barrel model. PSI-BLAST searches do not identify sequence similarities to known TIM-barrel proteins, and sensitive profile-profile searches indicate that the design sequence is distant from other naturally occurring TIM-barrel superfamilies, suggesting that Nature has sampled only a subset of the sequence space available to the TIM-barrel fold. The ability to design TIM barrels de novo opens new possibilities for custom-made enzymes.

A successful strategy for the recovering of active P21, an insoluble recombinant protein of Trypanosoma cruzi

Structural studies of proteins normally require large quantities of pure material that can only be obtained through heterologous expression systems and recombinant technique. In these procedures, large amounts of expressed protein are often found in the insoluble fraction, making protein purification from the soluble fraction inefficient, laborious, and costly. Usually, protein refolding is avoided due to a lack of experimental assays that can validate correct folding and that can compare the conformational population to that of the soluble fraction. Herein, we propose a validation method using simple and rapid 1D ¹H nuclear magnetic resonance (NMR) spectra that can efficiently compare protein samples, including individual information of the environment of each proton in the structure.

Octarellin VI: using rosetta to design a putative artificial (β/α)8 protein

The computational protein design protocol Rosetta has been applied successfully to a wide variety of protein engineering problems. Here the aim was to test its ability to design de novo a protein adopting the TIM-barrel fold, whose formation requires about twice as many residues as in the largest proteins successfully designed de novo to date. The designed protein, Octarellin VI, contains 216 residues. Its amino acid composition is similar to that of natural TIM-barrel proteins. When produced and purified, it showed a far-UV circular dichroism spectrum characteristic of folded proteins, with α-helical and β-sheet secondary structure. Its stable tertiary structure was confirmed by both tryptophan fluorescence and circular dichroism in the near UV. It proved heat stable up to 70°C. Dynamic light scattering experiments revealed a unique population of particles averaging 4 nm in diameter, in good agreement with our model. Although these data suggest the successful creation of an artificial α/β protein of more than 200 amino acids, Octarellin VI shows an apparent noncooperative chemical unfolding and low solubility.

Clusters of branched aliphatic side chains serve as cores of stability in the native state of the HisF TIM barrel protein

Imidazole-3-glycerol phosphate synthase is a heterodimeric allosteric enzyme that catalyzes consecutive reactions in imidazole biosynthesis through its HisF and HisH subunits. The unusually slow unfolding reaction of the isolated HisF TIM barrel domain from the thermophilic bacteria, Thermotoga maritima, enabled an NMR-based site-specific analysis of the main-chain hydrogen bonds that stabilize its native conformation. Very strong protection against exchange with solvent deuterium in the native state was found in a subset of buried positions in α-helices and pervasively in the underlying β-strands associated with a pair of large clusters of isoleucine, leucine and valine (ILV) side chains located in the α7(βα)8(βα)1-2 and α2(βα)3-6β7 segments of the (βα)8 barrel. The most densely packed region of the large cluster, α3(βα)4-6β7, correlates closely with the core of stability previously observed in computational, protein engineering and NMR dynamics studies, demonstrating a key role for this cluster in determining the thermodynamic and structural properties of the native state of HisF. When considered with the results of previous studies where ILV clusters were found to stabilize the hydrogen-bonded networks in folding intermediates for other TIM barrel proteins, it appears that clusters of branched aliphatic side chains can serve as cores of stability across the entire folding reaction coordinate of one of the most common motifs in biology.

Optimization of aromatic side chain size complementarity in the hydrophobic core of a designed coiled-coil

The coiled-coil structure plays an important roles, especially in protein assembly. Previously we constructed AAB-type heterotrimeric coiled-coils by manipulating the packing in the hydrophobic core using Trp and Ala residues, where one Trp and two Ala residues were placed in the hydrophobic core instead of three Ile residues. To optimize the packing complementarity in the hydrophobic core, we investigated the effects of introducing various aromatic amino acids on the formation of an AAB-type heterotrimeric coiled-coil, by circular dichroism, thermal stability, and nuclear magnetic resonance (NMR) studies. We found that the Phe residue was more suitable for heterotrimeric coiled-coil formation than the Trp residue, when combined with two Ala residues, whereas the Tyr and His residues did not induce the coiled-coil structure efficiently.

Two-metal ion, Ni(II) and Cu(II), binding alpha-helical coiled coil peptide

Metalloproteins are an attractive target for de novo design. Usually, natural proteins incorporate two or more (hetero- or homo-) metal ions into their frameworks to perform their functions, but the design of multiple metal-binding sites is usually difficult to achieve. Here, we undertook the de novo engineering of heterometal-binding sites, Ni(II) and Cu(II), into a designed coiled coil structure based on an isoleucine zipper (IZ) peptide. Previously, we described two peptides, IZ-3adH and IZ-3aH. The former has two His residues and forms a triple-stranded coiled coil after binding Ni(II), Zn(II), or Cu(II). The latter has one His residue, which allowed binding with Cu(II) and Zn(II), but not with Ni(II). On the basis of these properties, we newly designed IZ(5)-2a3adH as a heterometal-binding peptide. This peptide can bind Cu(II) and Ni(II) simultaneously in the hydrophobic core of the triple-stranded coiled coil. The first metal ion binding induced the folding of the peptide into the triple-stranded coiled coil, thereby promoting the second metal ion binding. This is the first example of a peptide that can bind two different metal ions. This construction should provide valuable insights for the de novo design of metalloproteins.

Folding of β/α-unit scrambled forms of S. cerevisiae triosephosphate isomerase: Evidence for autonomy of substructure formation and plasticity of hydrophobic and hydrogen bonding interactions in core of (β/α)8-barrel

The (beta/alpha)(8)-barrel domain consists of eight topologically equivalent supersecondary structural motifs known as beta/alpha-units. Each unit consists of a single beta-strand, an alpha-helix, and two loops. Evidence collected in recent years indicates that the (beta/alpha)(8)-barrel motif may not be a single, autonomously-folding domain, as was previously assumed. Segments of some (beta/alpha)(8)-barrels appear to fold autonomously. However, the extent to which this is true of various (beta/alpha)(8)-barrel domains remains to be explored. In this study, we have scrambled (reshuffled) the native order of beta/alpha-units (1-2-3-4-5-6-7-8) comprising the polypeptide chain of a model (beta/alpha)(8)-barrel from S. cerevisiae, triosephosphate isomerase (TIM). Total scrambling was effected in order to examine whether folding can still occur to yield beta/alpha-structures in spite of a global 'destruction' of native hydrophobic and hydrogen bonding interactions among beta/alpha-units, while still allowing the occurrence of native interactions within individual units. Our results demonstrate that scrambled full-barrel forms (2-4-6-8-1-3-5-7 and 1-3-5-7-2-4-6-8), as well as half-barrel (2-4-6-8) and quarter-barrel (1-3) forms of TIM fold into beta/alpha-structures that sustain tertiary and quaternary structural interactions. In particular, one variant (2-4-6-8-1-3-5-7) was found to fold and form a stable dimer with native-like structural content and other characteristics. Our results demonstrate that (beta/alpha)(8)-barrels can tolerate profound alterations of both strand-strand interactions responsible for the creation of the beta-barrel and the geometry of presentation of nonpolar sidechains into the hydrophobic core of the beta-barrel by individual beta-strands. These findings lend support to our recent proposal1 that a hierarchy of interactions probably regulates structure formation and stability in (beta/alpha)(8)-barrels, where folding proceeds successively through three stages: (i) the tentative formation of individual beta/alpha-units which associate through 'near-neighbor' diffusion-collision interactions into (ii) curved assemblies of multiple beta/alpha-units through sequence-independent hydrogen bonding of strands of neighboring units, leading finally to (iii) the association of curved (quarter/half-barrel) assemblies around a common hydrophobic core through packing interactions that remain plastic and amenable to change.

De novo backbone and sequence design of an idealized alpha/beta-barrel protein: evidence of stable tertiary structure

We have designed, synthesized, and characterized a 216 amino acid residue sequence encoding a putative idealized alpha/beta-barrel protein. The design was elaborated in two steps. First, the idealized backbone was defined with geometric parameters representing our target fold: a central eight parallel-stranded beta-sheet surrounded by eight parallel alpha-helices, connected together with short structural turns on both sides of the barrel. An automated sequence selection algorithm, based on the dead-end elimination theorem, was used to find the optimal amino acid sequence fitting the target structure. A synthetic gene coding for the designed sequence was constructed and the recombinant artificial protein was expressed in bacteria, purified and characterized. Far-UV CD spectra with prominent bands at 222nm and 208nm revealed the presence of alpha-helix secondary structures (50%) in fairly good agreement with the model. A pronounced absorption band in the near-UV CD region, arising from immobilized aromatic side-chains, showed that the artificial protein is folded in solution. Chemical unfolding monitored by tryptophan fluorescence revealed a conformational stability (DeltaG(H2O)) of 35kJ/mol. Thermal unfolding monitored by near-UV CD revealed a cooperative transition with an apparent T(m) of 65 degrees C. Moreover, the artificial protein did not exhibit any affinity for the hydrophobic fluorescent probe 1-anilinonaphthalene-8-sulfonic acid (ANS), providing additional evidence that the artificial barrel is not in the molten globule state, contrary to previously designed artificial alpha/beta-barrels. Finally, 1H NMR spectra of the folded and unfolded proteins provided evidence for specific interactions in the folded protein. Taken together, the results indicate that the de novo designed alpha/beta-barrel protein adopts a stable three-dimensional structure in solution. These encouraging results show that de novo design of an idealized protein structure of more than 200 amino acid residues is now possible, from construction of a particular backbone conformation to determination of an amino acid sequence with an automated sequence selection algorithm.

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.

Stability, catalytic versatility and evolution of the (beta alpha)(8)-barrel fold

The (beta alpha)(8)-barrel is a versatile single-domain protein fold that is adopted by a large number of enzymes. The (beta alpha)(8)-barrel fold has been used as a model to elucidate the structural basis of protein thermostability and in studies to interconvert catalytic activities or substrate specificities by rational design or directed evolution. Recently, the (beta alpha)(4)-half-barrel was identified as a possible structural subdomain.

Reverse engineering the (beta/alpha )8 barrel fold

The (beta/alpha)(8) barrel is the most commonly occurring fold among protein catalysts. To lay a groundwork for engineering novel barrel proteins, we investigated the amino acid sequence restrictions at 182 structural positions of the prototypical (beta/alpha)(8) barrel enzyme triosephosphate isomerase. Using combinatorial mutagenesis and functional selection, we find that turn sequences, alpha-helix capping and stop motifs, and residues that pack the interface between beta-strands and alpha-helices are highly mutable. Conversely, any mutation of residues in the central core of the beta-barrel, beta-strand stop motifs, and a single buried salt bridge between amino acids R189 and D227 substantially reduces catalytic activity. Four positions are effectively immutable: conservative single substitutions at these four positions prevent the mutant protein from complementing a triosephosphate isomerase knockout in Escherichia coli. At 142 of the 182 positions, mutation to at least one amino acid of a seven-letter amino acid alphabet produces a triosephosphate isomerase with wild-type activity. Consequently, it seems likely that (beta/alpha)(8) barrel structures can be encoded with a subset of the 20 amino acids. Such simplification would greatly decrease the computational burden of (beta/alpha)(8) barrel design.

Triose-phosphate isomerase (TIM) of the psychrophilic bacterium Vibrio marinus. Kinetic and structural properties

The purification and characterization of triose-phosphate isomerase from the psychrophilic bacterium Vibrio marinus (vTIM) is described. Crystal structures of the vTIM-sulfate complex and the vTIM-2-phosphoglycolate complex (at a 2.7-A resolution) are also presented. The optimal growth temperature of Vibrio marinus is 15 degrees C. Stability studies show that vTIM is an unstable protein with a half-life of only 10 min at 25 degrees C. The vTIM sequence is most closely related to the sequence of Escherichia coli TIM (eTIM) (66% identity), and several unique structural features described for eTIM are also seen in vTIM, but eTIM is considerably more stable. The Td values of vTIM and eTIM, determined by calorimetric studies, are 41 and 54 degrees C, respectively. Amino acid sequence comparison reveals that vTIM has an alanine in loop 8 (at position 238), whereas all other TIM sequences known to date have a serine. The vTIM mutant A238S was produced and characterized. Compared with wild type, the catalytic efficiency of the A238S mutant is somewhat reduced, and its stability is considerably increased.

New approaches in molecular structure prediction

In the past years, much effort has been put on the development of new methodologies and algorithms for the prediction of protein secondary and tertiary structures from (sequence) data; this is reviewed in detail. New approaches for these predictions such as neural network methods, genetic algorithms, machine learning, and graph theoretical methods are discussed. Secondary structure prediction algorithms were improved mostly by considering families of related proteins; however, for the reliable tertiary structure modeling of proteins, knowledge-based techniques are still preferred. Methods and examples with more or less successful results are described. Also, programs and parameterizations for energy minimisations, molecular dynamics, and electrostatic interactions have been improved, especially with respect to their former limits of applicability. Other topics discussed in this review include the use of traditional and on-line databases, the docking problem and surface properties of biomolecules, packing of protein cores, de novo design and protein engineering, prediction of membrane protein structures, the verification and reliability of model structures, and progress made with currently available software and computer hardware. In summary, the prediction of the structure, function, and other properties of a protein is still possible only within limits, but these limits continue to be moved.

Crystallization of a designed peptide from a molten globule ensemble

BACKGROUND

The design of amino acid sequences that adopt a desired three-dimensional fold has been of keen interest over the past decade. However, the design of proteins that adopt unique conformations is still a considerable problem. Until very recently, all of the designed proteins that have been extensively characterized possess the hallmarks of the molten globular state. Molten globular intermediates have been observed in both equilibrium and kinetic protein folding/stability studies, and understanding the forces that determine compact non-native states is critical for a comprehensive understanding of proteins. This paper describes the solution and early solid state characterization of peptides that form molten globular ensembles.

RESULTS & CONCLUSIONS

Crystals diffracting to 3.5 A resolution have been grown of a 16-residue peptide (alpha 1A) designed to form a tetramer of alpha-helices. In addition, a closely related peptide, alpha 1, has previously been shown to yield crystals that diffract to 1.2 A resolution. The solution properties of these two peptides were examined to determine whether their well defined crystalline conformations were retained in solution. On the basis of an examination of their NMR spectra, sedimentation equilibria, thermal unfolding, and ANS binding, it is concluded that the peptides form alpha-helical aggregates with properties similar to those of the molten globule state. Thus, for these peptides, the process of crystallization bears many similarities to models of protein folding. Upon dissolution, the peptides rapidly assume compact molten globular states similar to the molten globule like intermediates that are formed at short times after refolding is initiated. Following a rate-determining nucleation step, the peptides crystallize into a single or a small number of conformations in a process that mimics the formation of native structure in proteins.

Protein design: a hierarchic approach

The de novo design of peptides and proteins has recently emerged as an approach for investigating protein structure and function. Designed, helical peptides provide model systems for dissecting and quantifying the multiple interactions that stabilize secondary structure formation. De novo design is also useful for exploring the features that specify the stoichiometry and stability of alpha-helical coiled coils and for defining the requirements for folding into structures that resemble native, functional proteins. The design process often occurs in a series of discrete steps. Such steps reflect the hierarchy of forces required for stabilizing tertiary structures, beginning with hydrophobic forces and adding more specific interactions as required to achieve a unique, functional protein.

New strategies in protein design

Initially, it was hoped that very simple rules could be sued to design proteins that embody all the characteristics of natural proteins. Indeed, with single-domain proteins as targets, it has been possible to design proteins that adopt the desired global fold. Yet, designed proteins with well defined structures and properties that mimic those of natural proteins remain elusive. Recent efforts in protein design have been directed toward addressing the basis for non-native characteristics in most protein designs. Although it is clear that specific tertiary interactions between all residues in a protein contribute to the final folded state, much attention has been placed on optimizing the packing of side chains in the hydrophobic core, with substantial success.

A quantitative methodology for the de novo design of proteins

We have developed a general quantitative methodology for designing proteins de novo, which automatically produces sequences for any given plausible protein structure. The method incorporates statistical information, a theoretical description of protein structure, and motifs described in the literature. A model system embodying a portion of the quantitative methodology has been used to design many protein sequences for the phage 434 Cro and fibronectin type III domain folds, as well as several other structures. Residue sequences selected by this prototype share no significant identity with any natural protein. Nonetheless, 3-dimensional models of the designed sequences appear generally plausible. When examined using secondary structure prediction methods and profile analysis, the designed sequences generally score considerably better than the natural ones. The designed sequences are also in reasonable agreement with a sequence template. This quantitative methodology is likely to be capable of successfully designing new proteins and yielding fundamental insights about the determinants of protein structure.