De novo heme proteins from designed combinatorial libraries.

We previously reported the design of a library of de novo amino acid sequences targeted to fold into four-helix bundles. The design of these sequences was based on a "binary code" strategy, in which the patterning of polar and nonpolar amino acids is specified explicitly, but the exact identities of the side chains is varied extensively (Kamtekar S, Schiffer JM, Xiong H, Babik JM, Hecht MH, 1993, Science 262:1680-1685). Because of this variability, the resulting collection of amino acid sequences may include de novo proteins capable of binding biologically important cofactors. To probe for such binding, the de novo sequences were screened for their ability to bind the heme cofactor. Among an initial collection of 30 binary code sequences, 15 are shown to bind heme and form bright red complexes. Characterization of several of these de novo heme proteins demonstrated that their absorption spectra and resonance Raman spectra resemble those of natural cytochromes. Because the design of these sequences is based on global features of polar/ nonpolar patterning, the finding that half of them bind heme highlights the power of the binary code strategy, and demonstrates that isolating de novo heme proteins does not require explicit design of the cofactor binding site. Because bound heme plays a key role in the functions of many natural proteins, these results suggest that binary code sequences may serve as initial prototypes for the development of large collections of functionally active de novo proteins.     

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.     

Design of proteins with hydrophobic and polar amino acids

A two amino acid (hydrophobic and polar) scheme is used to perform the design on target conformations corresponding to the native states of 20 single chain proteins. Strikingly, the percentage of successful identification of the nature of the residues benchmarked against naturally occurring proteins and their homologues is around 75%, independent of the complexity of the design procedure. Typically, the lowest success rate occurs for residues such as alanine that have a high secondary structure functionality. Using a simple lattice model, we argue that one possible shortcoming of the model studied may involve the coarse-graining of the 20 kinds of amino acids into just two effective types. Proteins 32:80–87, 1998. © 1998 Wiley-Liss, Inc.     

Cooperative thermal denaturation of proteins designed by binary patterning of polar and nonpolar amino acids

We previously reported a combinatorial strategy for designing alpha-helical proteins by assigning only the binary patterning of polar or nonpolar residues [Kamtekar, S., Schiffer, J. M., Xiong, H. Y., Babik, J. M., and Hecht, M. H. (1993) Science 262, 1680-1685]. Here we describe the finding that approximately half of the proteins in the original collection display some level of cooperativity in their thermal denaturation profiles. Many are monomeric in solution, demonstrating that the observed cooperativity is not merely a consequence of oligomerization. These findings demonstrate that although the combinatorial nature of the design strategy precludes explicit design of side-chain packing, binary patterning incorporates sufficient sequence information to generate de novo proteins with cooperatively folded structures. As binary partitioning of polar and nonpolar amino acids is an intrinsic part of the genetic code, these findings may bear on the early evolution of native proteins.     

Cold hardening induces transfer of fatty acids between polar and nonpolar lipid pools in the Arctic collembollan Megaphorura arctica

Cold hardiness in the Arctic Collembola Megaphorura arctica (Tullberg), formerly Onychiurus arcticus, has been the subject of extensive studies over the last decade. This species employs an unusual strategy known as cryoprotective dehydration to survive winter temperatures as low as ?25 °C. To expand knowledge of cryoprotective dehydration in M. arctica, the present study investigates how a reduction in ambient temperature affects the fatty acid composition of the total body lipid content along with polar (mainly membrane phospholipids) and nonpolar (mainly triacylglycerols) lipids. Most ectothermic animals compensate for changes in fluidity by regulating fatty acid composition, a process often described as homeoviscous adaptation. In M. arctica, changes in the fatty acid composition of total body lipid content during cold treatment are only moderate, with no clear pattern emerging. However, the levels of unsaturated fatty acids in the polar lipids increase with cold exposure, largely attributable to 16 : 1(n? 7), 18 : 1(n? 9), 18 : 3(n? 6) and 18 : 3(n? 3), whereas unsaturated fatty acid levels in the nonpolar lipids correspondingly decrease. These results suggest a reallocation of fatty acids between the two lipid pools as a response to a temperature reduction of 6 °C. Because of hypometabolism, a characteristic of cold adaptation, such a mechanism could be less energy demanding than de novo synthesis of fatty acids and may comprise part of an adaptive homeostatic response.     

Modern and prebiotic amino acids support distinct structural profiles in proteins

The earliest proteins had to rely on amino acids available on early Earth before the biosynthetic pathways for more complex amino acids evolved. In extant proteins, a significant fraction of the ‘late’ amino acids (such as Arg, Lys, His, Cys, Trp and Tyr) belong to essential catalytic and structure-stabilizing residues. How (or if) early proteins could sustain an early biosphere has been a major puzzle. Here, we analysed two combinatorial protein libraries representing proxies of the available sequence space at two different evolutionary stages. The first is composed of the entire alphabet of 20 amino acids while the second one consists of only 10 residues (ASDGLIPTEV) representing a consensus view of plausibly available amino acids through prebiotic chemistry. We show that compact conformations resistant to proteolysis are surprisingly similarly abundant in both libraries. In addition, the early alphabet proteins are inherently more soluble and refoldable, independent of the general Hsp70 chaperone activity. By contrast, chaperones significantly increase the otherwise poor solubility of the modern alphabet proteins suggesting their coevolution with the amino acid repertoire. Our work indicates that while both early and modern amino acids are predisposed to supporting protein structure, they do so with different biophysical properties and via different mechanisms.     

Stably folded de novo proteins from a designed combinatorial library

Binary patterning of polar and nonpolar amino acids has been used as the key design feature for constructing large combinatorial libraries of de novo proteins. Each position in a binary patterned sequence is designed explicitly to be either polar or nonpolar; however, the precise identities of these amino acids are varied extensively. The combinatorial underpinnings of the "binary code" strategy preclude explicit design of particular side chains at specified positions. Therefore, packing interactions cannot be specified a priori. To assess whether the binary code strategy can nonetheless produce well-folded de novo proteins, we constructed a second-generation library based upon a new structural scaffold designed to fold into 102-residue four-helix bundles. Characterization of five proteins chosen arbitrarily from this new library revealed that (1) all are alpha-helical and quite stable; (2) four of the five contain an abundance of tertiary interactions indicative of well-ordered structures; and (3) one protein forms a well-folded structure with native-like features. The proteins from this new 102-residue library are substantially more stable and dramatically more native-like than those from an earlier binary patterned library of 74-residue sequences. These findings demonstrate that chain length is a crucial determinant of structural order in libraries of de novo four-helix bundles. Moreover, these results show that the binary code strategy--if applied to an appropriately designed structural scaffold--can generate large collections of stably folded and/or native-like proteins.     

Comparison of the structures and stabilities of coiled‐coil proteins containing hexafluoroleucine and t‐butylalanine provides insight into the stabilizing effects of highly fluorinated amino acid side‐chains

Highly fluorinated analogs of hydrophobic amino acids are well known to increase the stability of proteins toward thermal unfolding and chemical denaturation, but there is very little data on the structural consequences of fluorination. We have determined the structures and folding energies of three variants of a de novo designed 4‐helix bundle protein whose hydrophobic cores contain either hexafluoroleucine (hFLeu) or t‐butylalanine (tBAla). Although the buried hydrophobic surface area is the same for all three proteins, the incorporation of tBAla causes a rearrangement of the core packing, resulting in the formation of a destabilizing hydrophobic cavity at the center of the protein. In contrast, incorporation of hFLeu, causes no changes in core packing with respect to the structure of the nonfluorinated parent protein which contains only leucine in the core. These results support the idea that fluorinated residues are especially effective at stabilizing proteins because they closely mimic the shape of the natural residues they replace while increasing buried hydrophobic surface area.     

Divergent evolution of a bifunctional de novo protein

Primordial proteins, the evolutionary ancestors of modern sequences, are presumed to have been minimally active and nonspecific. Following eons of selective pressure, these early progenitors evolved into highly active and specific proteins. While evolutionary trajectories from poorly active and multifunctional generalists toward highly active specialists likely occurred many times in evolutionary history, such pathways are difficult to reconstruct in natural systems, where primordial sequences are lost to time. To test the hypothesis that selection for enhanced activity leads to a loss of promiscuity, we evolved a de novo designed bifunctional protein. The parental protein, denoted Syn‐IF, was chosen from a library of binary patterned 4‐helix bundles. Syn‐IF was shown previously to rescue two different auxotrophic strains of E. coli: ΔilvA and Δfes. These two strains contain deletions for proteins with very different biochemical functions; IlvA is involved in isoleucine biosynthesis, while Fes is involved in iron assimilation. In two separate experiments, Syn‐IF, was evolved for faster rescue of either ΔilvA or Δfes. Following multiple rounds of mutagenesis, two new proteins were selected, each capable of rescuing the selected function significantly faster than the parental protein. In each case, the evolved protein also lost the ability to rescue the unselected function. In both evolutionary trajectories, the original bifunctional generalist was evolved into a monofunctional specialist with enhanced activity.     

Design, synthesis, and characterization of a novel hemoprotein

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

What limits the primary sequence space of natural proteins?

Abstract

Number of naturally occurring primary sequences of proteins is an infinitesimally small subset of the possible number of primary sequences that can be synthesized using 20 amino acids. Prevailing views ascribe this to slow and incremental mutational/selection evolutionary mechanisms. However, considering the large number of avenues available in form of diversity of emerging/evolving and/or disappearing living systems for exploring the primary sequence space over the evolutionary time scale of ～3.5 billion years, this remains a conjecture. Therefore, to investigate primary sequence space limitations, we carried out a systematic study for finding primary sequences absent in nature. We report the discovery of the smallest peptide sequence “Cysteine-Glutamine-Tryptophan-Tryptophan” that is not found in over half-a-million curated protein sequences in the Uniprot (Swiss-Prot) database. Additionally, we report a library of 83605 pentapeptides that are not found in any of the known protein sequences. Compositional analyses of these absent primary sequences yield a remarkably strong power relationship between the percentage occurrence of individual amino acids in all known protein sequences and their respective frequency of occurrence in the absent peptides, regardless of their specific position in the sequences. If random evolutionary mechanisms were responsible for limitations to the primary sequence space, then one would not expect any relationship between compositions of available and absent primary sequences. Thus, we conclusively show that stoichiometric constraints on amino acids limit the primary sequence space of proteins in nature. We discuss the possibly profound implications of our findings in both evolutionary and synthetic biology.

Communicated by Ramaswamy H. Sarma     

Cyclic oligomer design with de novo αβ‐proteins

We have previously shown that monomeric globular αβ‐proteins can be designed de novo with considerable control over topology, size, and shape. In this paper, we investigate the design of cyclic homo‐oligomers from these starting points. We experimented with both keeping the original monomer backbones fixed during the cyclic docking and design process, and allowing the backbone of the monomer to conform to that of adjacent subunits in the homo‐oligomer. The latter flexible backbone protocol generated designs with shape complementarity approaching that of native homo‐oligomers, but experimental characterization showed that the fixed backbone designs were more stable and less aggregation prone. Designed C2 oligomers with β‐strand backbone interactions were structurally confirmed through x‐ray crystallography and small‐angle X‐ray scattering (SAXS). In contrast, C3‐C5 designed homo‐oligomers with primarily nonpolar residues at interfaces all formed a range of oligomeric states. Taken together, our results suggest that for homo‐oligomers formed from globular building blocks, improved structural specificity will be better achieved using monomers with increased shape complementarity and with more polar interfaces.     

Structure and dynamics of de novo proteins from a designed superfamily of 4-helix bundles

Libraries of de novo proteins provide an opportunity to explore the structural and functional potential of biological molecules that have not been biased by billions of years of evolutionary selection. Given the enormity of sequence space, a rational approach to library design is likely to yield a higher fraction of folded and functional proteins than a stochastic sampling of random sequences. We previously investigated the potential of library design by binary patterning of hydrophobic and hydrophilic amino acids. The structure of the most stable protein from a binary patterned library of de novo 4-helix bundles was solved previously and shown to be consistent with the design. One structure, however, cannot fully assess the potential of the design strategy, nor can it account for differences in the stabilities of individual proteins. To more fully probe the quality of the library, we now report the NMR structure of a second protein, S-836. Protein S-836 proved to be a 4-helix bundle, consistent with design. The similarity between the two solved structures reinforces previous evidence that binary patterning can encode stable, 4-helix bundles. Despite their global similarities, the two proteins have cores that are packed at different degrees of tightness. The relationship between packing and dynamics was probed using the Modelfree approach, which showed that regions containing a high frequency of chemical exchange coincide with less well-packed side chains. These studies show (1) that binary patterning can drive folding into a particular topology without the explicit design of residue-by-residue packing, and (2) that within a superfamily of binary patterned proteins, the structures and dynamics of individual proteins are modulated by the identity and packing of residues in the hydrophobic core.     

Restricted sidechain plasticity in the structures of native proteins and complexes

Protein-design methodology can now generate models of protein structures and interfaces with computed energies in the range of those of naturally occurring structures. Comparison of the properties of native structures and complexes to isoenergetic design models can provide insight into the properties of the former that reflect selection pressure for factors beyond the energy of the native state. We report here that sidechains in native structures and interfaces are significantly more constrained than designed interfaces and structures with equal computed binding energy or stability, which may reflect selection against potentially deleterious non-native interactions.     

Step‐by‐step design of proteins for small molecule interaction: A review on recent milestones

Protein design is the field of synthetic biology that aims at developing de novo custom‐made proteins and peptides for specific applications. Despite exploring an ambitious goal, recent computational advances in both hardware and software technologies have paved the way to high‐throughput screening and detailed design of novel folds and improved functionalities. Modern advances in the field of protein design for small molecule targeting are described in this review, organized in a step‐by‐step fashion: from the conception of a new or upgraded active binding site, to scaffold design, sequence optimization, and experimental expression of the custom protein. In each step, contemporary examples are described, and state‐of‐the‐art software is briefly explored.     

Histidine placement in de novo-designed heme proteins.

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

Rational design of TNFα binding proteins based on the de novo designed protein DS119

De novo protein design offers templates for engineering tailor‐made protein functions and orthogonal protein interaction networks for synthetic biology research. Various computational methods have been developed to introduce functional sites in known protein structures. De novo designed protein scaffolds provide further opportunities for functional protein design. Here we demonstrate the rational design of novel tumor necrosis factor alpha (TNFα) binding proteins using a home‐made grafting program AutoMatch. We grafted three key residues from a virus 2L protein to a de novo designed small protein, DS119, with consideration of backbone flexibility. The designed proteins bind to TNFα with micromolar affinities. We further optimized the interface residues with RosettaDesign and significantly improved the binding capacity of one protein Tbab1‐4. These designed proteins inhibit the activity of TNFα in cellular luciferase assays. Our work illustrates the potential application of the de novo designed protein DS119 in protein engineering, biomedical research, and protein sequence‐structure‐function studies.     

Dissecting contact potentials for proteins: relative contributions of individual amino acids

Mimotopes mimic the three-dimensional topology of an antigen epitope, and are frequently recognized by antibodies with affinities comparable to those obtained for the original antibody-antigen interaction. Peptides and anti-idiotypic antibodies are two classes of protein mimotopes that mimic the topology (but not necessarily the sequence) of the parental antigen. In this study, we combine these two classes by selecting mimotopes based on single domain IgNAR antibodies, which display exceptionally long CDR3 loop regions (analogous to a constrained peptide library) presented in the context of an immunoglobulin framework with adjacent and supporting CDR1 loops. By screening an in vitro phage-display library of IgNAR variable domains (V(NAR)s) against the target antigen monoclonal antibody MAb5G8, we obtained four potential mimotopes. MAb5G8 targets a linear tripeptide epitope (AYP) in the flexible signal sequence of the Plasmodium falciparum Apical Membrane Antigen-1 (AMA1), and this or similar motifs were detected in the CDR loops of all four V(NAR)s. The V(NAR)s, 1-A-2, -7, -11, and -14, were demonstrated to bind specifically to this paratope by competition studies with an artificial peptide and all showed enhanced affinities (3-46 nM) compared to the parental antigen (175 nM). Crystallographic studies of recombinant proteins 1-A-7 and 1-A-11 showed that the SYP motifs on these V(NAR)s presented at the tip of the exposed CDR3 loops, ideally positioned within bulge-like structures to make contact with the MAb5G8 antibody. These loops, in particular in 1-A-11, were further stabilized by inter- and intra- loop disulphide bridges, hydrogen bonds, electrostatic interactions, and aromatic residue packing. We rationalize the higher affinity of the V(NAR)s compared to the parental antigen by suggesting that adjacent CDR1 and framework residues contribute to binding affinity, through interactions with other CDR regions on the antibody, though of course definitive support of this hypothesis will rely on co-crystallographic studies. Alternatively, the selection of mimotopes from a large (<4 x 10(8)) constrained library may have allowed selection of variants with even more favorable epitope topologies than present in the original antigenic structure, illustrating the power of in vivo selection of mimotopes from phage-displayed molecular libraries.     

Predicting the side‐chain dihedral angle distributions of nonpolar,aromatic, and polar amino acids using hard sphere models

The side‐chain dihedral angle distributions of all amino acids have been measured from myriad high‐resolution protein crystal structures. However, we do not yet know the dominant interactions that determine these distributions. Here, we explore to what extent the defining features of the side‐chain dihedral angle distributions of different amino acids can be captured by a simple physical model. We find that a hard‐sphere model for a dipeptide mimetic that includes only steric interactions plus stereochemical constraints is able to recapitulate the key features of the back‐bone dependent observed amino acid side‐chain dihedral angle distributions of Ser, Cys, Thr, Val, Ile, Leu, Phe, Tyr, and Trp. We find that for certain amino acids, performing the calculations with the amino acid of interest in the central position of a short α‐helical segment improves the match between the predicted and observed distributions. We also identify the atomic interactions that give rise to the differences between the predicted distributions for the hard‐sphere model of the dipeptide and that of the α‐helical segment. Finally, we point out a case where the hard‐sphere plus stereochemical constraint model is insufficient to recapitulate the observed side‐chain dihedral angle distribution, namely the distribution P(χ₃) for Met. Proteins 2014; 82:2574–2584. © 2014 Wiley Periodicals, Inc.