From synthetic coiled coils to functional proteins: automated design of a receptor for the calmodulin-binding domain of calcineurin

De novo design of high-affinity binders of bioactive helical peptides

Many peptide hormones form an α-helix on binding their receptors1-4, and sensitive methods for their detection could contribute to better clinical management of disease5. De novo protein design can now generate binders with high affinity and specificity to structured proteins6,7. However, the design of interactions between proteins and short peptides with helical propensity is an unmet challenge. Here we describe parametric generation and deep learning-based methods for designing proteins to address this challenge. We show that by extending RFdiffusion8 to enable binder design to flexible targets, and to refining input structure models by successive noising and denoising (partial diffusion), picomolar-affinity binders can be generated to helical peptide targets by either refining designs generated with other methods, or completely de novo starting from random noise distributions without any subsequent experimental optimization. The RFdiffusion designs enable the enrichment and subsequent detection of parathyroid hormone and glucagon by mass spectrometry, and the construction of bioluminescence-based protein biosensors. The ability to design binders to conformationally variable targets, and to optimize by partial diffusion both natural and designed proteins, should be broadly useful.

Chirality manifestation in elastic coupling between the layers of double-walled carbon nanotubes

The search for new relatively easy physicochemical methods for the structural identification of carbon nanotubes represents a key challenge. Here, analyzing the experimental data on double-walled carbon nanotubes (DWCNTs) obtained by us and taken from the literature, we have expressed the magnitude of elastic coupling between two tubular walls forming a DWCNT as a simple function dependent not only on DWCNT diameters but also on the difference between the chirality angles of the constituent nanotubes. To get this quite unexpected result, which allows us to relate more precisely the structural parameters of a DWCNT with frequencies of its radial breathing-like modes (RBLM), we have developed a new model for the RBLM dynamics that takes into account a possible deposition of water molecules from ambient air onto the DWCNT surface. The model constructed allows us to predict the higher prevalence of DWCNTs comprising two walls with identical handedness. The application of the results obtained for the identification of DWCNTs is also considered.

De Novo Design of Four-Helix Bundle Metalloproteins: One Scaffold, Diverse Reactivities

De novo protein design represents an attractive approach for testing and extending our understanding of metalloprotein structure and function. Here, we describe our work on the design of DF (Due Ferri or two-iron in Italian), a minimalist model for the active sites of much larger and more complex natural diiron and dimanganese proteins. In nature, diiron and dimanganese proteins protypically bind their ions in 4-Glu, 2-His environments, and they catalyze diverse reactions, ranging from hydrolysis, to O₂-dependent chemistry, to decarbonylation of aldehydes. In the design of DF, the position of each atom-including the backbone, the first-shell ligands, the second-shell hydrogen-bonded groups, and the well-packed hydrophobic core-was bespoke using precise mathematical equations and chemical principles. The first member of the DF family was designed to be of minimal size and complexity and yet to display the quintessential elements required for binding the dimetal cofactor. After thoroughly characterizing its structural, dynamic, spectroscopic, and functional properties, we added additional complexity in a rational stepwise manner to achieve increasingly sophisticated catalytic functions, ultimately demonstrating substrate-gated four-electron reduction of O₂ to water. We also briefly describe the extension of these studies to the design of proteins that bind nonbiological metal cofactors (a synthetic porphyrin and a tetranuclear cluster), and a Zn²⁺/proton antiporting membrane protein. Together these studies demonstrate a successful and generally applicable strategy for de novo metalloprotein design, which might indeed mimic the process by which primordial metalloproteins evolved. We began the design process with a highly symmetrical backbone and binding site, by using point-group symmetry to assemble the secondary structures that position the amino acid side chains required for binding. The resulting models provided a rough starting point and initial parameters for the subsequent precise design of the final protein using modern methods of computational protein design. Unless the desired site is itself symmetrical, this process requires reduction of the symmetry or lifting it altogether. Nevertheless, the initial symmetrical structure can be helpful to restrain the search space during assembly of the backbone. Finally, the methods described here should be generally applicable to the design of highly stable and robust catalysts and sensors. There is considerable potential in combining the efficiency and knowledge base associated with homogeneous metal catalysis with the programmability, biocompatibility, and versatility of proteins. While the work reported here focuses on testing and learning the principles of natural metalloproteins by designing and studying proteins one at a time, there is also considerable potential for using designed proteins that incorporate both biological and nonbiological metal ion cofactors for the evolution of novel catalysts.

Comparing side chain packing in soluble proteins, protein-protein interfaces, and transmembrane proteins

We compare side chain prediction and packing of core and non-core regions of soluble proteins, protein-protein interfaces, and transmembrane proteins. We first identified or created comparable databases of high-resolution crystal structures of these 3 protein classes. We show that the solvent-inaccessible cores of the 3 classes of proteins are equally densely packed. As a result, the side chains of core residues at protein-protein interfaces and in the membrane-exposed regions of transmembrane proteins can be predicted by the hard-sphere plus stereochemical constraint model with the same high prediction accuracies (>90%) as core residues in soluble proteins. We also find that for all 3 classes of proteins, as one moves away from the solvent-inaccessible core, the packing fraction decreases as the solvent accessibility increases. However, the side chain predictability remains high (80% within 30°) up to a relative solvent accessibility, rSASA≲0.3, for all 3 protein classes. Our results show that ≈40% of the interface regions in protein complexes are "core", that is, densely packed with side chain conformations that can be accurately predicted using the hard-sphere model. We propose packing fraction as a metric that can be used to distinguish real protein-protein interactions from designed, non-binding, decoys. Our results also show that cores of membrane proteins are the same as cores of soluble proteins. Thus, the computational methods we are developing for the analysis of the effect of hydrophobic core mutations in soluble proteins will be equally applicable to analyses of mutations in membrane proteins.

Empirical and computational design of iron-sulfur cluster proteins

Here, we compare two approaches of protein design. A computational approach was used in the design of the coiled-coil iron-sulfur protein, CCIS, as a four helix bundle binding an iron-sulfur cluster within its hydrophobic core. An empirical approach was used for designing the redox-chain maquette, RCM as a four-helix bundle assembling iron-sulfur clusters within loops and one heme in the middle of its hydrophobic core. We demonstrate that both ways of design yielded the desired proteins in terms of secondary structure and cofactors assembly. Both approaches, however, still have much to improve in predicting conformational changes in the presence of bound cofactors, controlling oligomerization tendency and stabilizing the bound iron-sulfur clusters in the reduced state. Lessons from both ways of design and future directions of development are discussed. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Computational design of virus-like protein assemblies on carbon nanotube surfaces

There is a general need for the engineering of protein-like molecules that organize into geometrically specific superstructures on molecular surfaces, directing further functionalization to create richly textured, multilayered assemblies. Here we describe a computational approach whereby the surface properties and symmetry of a targeted surface define the sequence and superstructure of surface-organizing peptides. Computational design proceeds in a series of steps that encode both surface recognition and favorable intersubunit packing interactions. This procedure is exemplified in the design of peptides that assemble into a tubular structure surrounding single-walled carbon nanotubes (SWNTs). The geometrically defined, virus-like coating created by these peptides converts the smooth surfaces of SWNTs into highly textured assemblies with long-scale order, capable of directing the assembly of gold nanoparticles into helical arrays along the SWNT axis.

Specificity in computational protein design

A long-standing goal of computational protein design is to create proteins similar to those found in Nature. One motivation is to harness the exquisite functional capabilities of proteins for our own purposes. The extent of similarity between designed and natural proteins also reports on how faithfully our models represent the selective pressures that determine protein sequences. As the field of protein design shifts emphasis from reproducing native-like protein structure to function, it has become important that these models treat the notion of specificity in molecular interactions. Although specificity may, in some cases, be achieved by optimization of a desired protein in isolation, methods have been developed to address directly the desire for proteins that exhibit specific functions and interactions.

A rational route to probing membrane proteins

Computationally designed peptides against transmembrane helices point to new approaches to studying membrane protein function.

A recent report describes the design of short peptides that bind specifically to transmembrane regions of integrins, providing an exciting tool for probing the biology of membrane proteins.

Computational design of peptides that target transmembrane helices

A variety of methods exist for the design or selection of antibodies and other proteins that recognize the water-soluble regions of proteins; however, companion methods for targeting transmembrane (TM) regions are not available. Here, we describe a method for the computational design of peptides that target TM helices in a sequence-specific manner. To illustrate the method, peptides were designed that specifically recognize the TM helices of two closely related integrins (alphaIIbbeta3 and alphavbeta3) in micelles, bacterial membranes, and mammalian cells. These data show that sequence-specific recognition of helices in TM proteins can be achieved through optimization of the geometric complementarity of the target-host complex.

A novel fluorescent probe for protein binding and folding studies:p-cyano-phenylalanine

Recently, it is has been shown that the C=N stretching vibration of a non-natural amino acid, p-cyano-phenylalanine (PheCN), could be used as an infrared reporter of local environment. Here, we further showed that the fluorescence emission of PheCN is also sensitive to solvent and, therefore, could be used as a novel optical probe for protein binding and folding studies. Moreover, we found that the fluorescence quantum yield of PheCN is nearly five times larger than that of phenylalanine and, more importantly, can be selectively excited even when other aromatic amino acids are present, thus making it a more versatile fluorophore. To test the feasibility of using PheCN as a practical fluorescent probe, we studied the binding of calmodulin (CaM) to a peptide derived from the CaM-binding domain of skeletal muscle myosin light chain kinase (MLCK). The peptide (MLCK3CN) contains a single PheCN residue and has been shown to bind to CaM with high affinity. As expected, addition of CaM into a MLCK3CN solution resulted in quenching of the PheCN fluorescence. A series of stochiometric titrations further allowed us to determine the binding affinity (Kd) of this peptide to CaM. Taken together, these results indicated that the PheCN fluorescence is sensitive to environment and could be applicable to a wide variety of biological problems.

Design of a heterospecific, tetrameric, 21-residue miniprotein with mixed alpha/beta structure

The study of short, autonomously folding peptides, or "miniproteins," is important for advancing our understanding of protein stability and folding specificity. Although many examples of synthetic alpha-helical structures are known, relatively few mixed alpha/beta structures have been successfully designed. Only one mixed-secondary structure oligomer, an alpha/beta homotetramer, has been reported thus far. In this report, we use structural analysis and computational design to convert this homotetramer into the smallest known alpha/beta-heterotetramer. Computational screening of many possible sequence/structure combinations led efficiently to the design of short, 21-residue peptides that fold cooperatively and autonomously into a specific complex in solution. A 1.95 A crystal structure reveals how steric complementarity and charge patterning encode heterospecificity. The first- and second-generation heterotetrameric miniproteins described here will be useful as simple models for the analysis of protein-protein interaction specificity and as structural platforms for the further elaboration of folding and function.

Dual surface selection methodology for the identification of thrombin binding epitopes from hotspot biased phage-display libraries

Protein libraries biased towards amino-acid residues found at so-called 'hotspots' were incorporated into the beta-sheet region of the thermostable variant (HTB1) of the B1 domain of the immunoglobulin (IgG) binding protein G and expressed as gene 3 fusions on M13 bacteriophage. The HTB1 library (2.2 x 10(9)) variants with a minimal 12 amino acid basis set were selected for binding IgG, to ensure structural conservation, and subsequently to thrombin to evolve a thrombin-binding function. We believe that this dual surface selection strategy will have great utility in evolving new bi-functional proteins without compromising structure. Furthermore the discrete beta-sheet epitopes identified by our methodology will lend itself to small-molecule mimicry of beta-sheets.

Computational design of protein–protein interactions

Computational protein design strategies have been developed to reengineer protein-protein interfaces in an automated, generalizable fashion. In the past two years, these methods have been successfully applied to generate chimeric proteins and protein pairs with specificities different from naturally occurring protein-protein interactions. Although there are shortcomings in current approaches, both in the way conformational space is sampled and in the energy functions used to evaluate designed conformations, the successes suggest we are now entering an era in which computational methods can be used to modulate, reengineer and design protein-protein interaction networks in living cells.

Metal-Chelating Amino Acids As Building Blocks For Synthetic Receptors Sensing Metal Ions And Histidine-Tagged Proteins

Protein structure and function rely on a still not fully understood interplay of energetic and entropic constraints defined by the permutation of the twenty genetically encoded amino acids. Many attempts have been undertaken to design peptide-peptide interaction pairs and synthetic receptors de novo by using this limited number of building blocks. We describe a rational approach to creating a building block based on a tailored metal-chelating amino acid. Nepsilon,Nepsilon-bis(carboxymethyl)-L-lysine can be flexibly introduced into peptides by 9-fluorenylmethoxycarbonyl solid-phase chemistry. The corresponding metal-chelating peptides act as metal sensors and synthetic receptors for histidine-tagged proteins. These biochemical tweezers will open new ways to control protein-protein interactions, to design peptide-based interaction pairs, or to generate switchable protein function.

Protein grafting of an HIV-1-inhibiting epitope

Protein grafting, the transfer of a binding epitope of one ligand onto the surface of another protein, is a potentially powerful technique for presenting peptides in preformed and active three-dimensional conformations. Its utility, however, has been limited by low biological activity of the designed ligands and low tolerance of the protein scaffolds to surface substitutions. Here, we graft the complete binding epitope (19 nonconsecutive amino acids with a solvent-accessible surface area of >2,000 A2) of an HIV-1 C-peptide, which is derived from the C-terminal region of HIV-1 gp41 and potently inhibits HIV-1 entry into cells, onto the surface of a GCN4 leucine zipper. The designed peptide, named C34coil, displays a potent antiviral activity approaching that of the native ligand. Moreover, whereas the linear C-peptide is unstructured and sensitive to degradation by proteases, C34coil is well structured, conformationally stable, and exhibits increased resistance to proteolytic degradation compared with the linear peptide. In addition to being a structured antiviral inhibitor, C34coil may also serve as the basis for the development of an alternative class of immunogens. This study demonstrates that "one-shot" protein grafting, without subsequent rounds of optimization, can be used to create ligands with structural conformations and improved biomedical properties.

Preorganization of molecular binding sites in designed diiron proteins

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

Modulating calmodulin binding specificity through computational protein design

We report the computational redesign of the protein-binding interface of calmodulin (CaM), a small, ubiquitous Ca(2+)-binding protein that is known to bind to and regulate a variety of functionally and structurally diverse proteins. The CaM binding interface was optimized to improve binding specificity towards one of its natural targets, smooth muscle myosin light chain kinase (smMLCK). The optimization was performed using optimization of rotamers by iterative techniques (ORBIT), a protein design program that utilizes a physically based force-field and the Dead-End Elimination theorem to compute sequences that are optimal for a given protein scaffold. Starting from the structure of the CaM-smMLCK complex, the program considered 10(22) amino acid residue sequences to obtain the lowest-energy CaM sequence. The resulting eightfold mutant, CaM_8, was constructed and tested for binding to a set of seven CaM target peptides. CaM_8 displayed high binding affinity to the smMLCK peptide (1.3nM), similar to that of the wild-type protein (1.8nM). The affinity of CaM_8 to six other target peptides was reduced, as intended, by 1.5-fold to 86-fold. Hence, CaM_8 exhibited increased binding specificity, preferring the smMLCK peptide to the other targets. Studies of this type may increase our understanding of the origins of binding specificity in protein-ligand complexes and may provide valuable information that can be used in the design of novel protein receptors and/or ligands.

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

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

A hierarchic approach to the design of hexameric helical barrels

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

A designed protein with packing between left-handed and right-handed helices

A common motif in protein structures is the assembly of alpha-helices. Natural alpha-helical assemblies, such as helical bundles and coiled coils, consist of multiple right-handed alpha-helices. Here we design a protein complex containing both left-handed and right-handed helices, with peptides of D- and L-amino acids, respectively. The two peptides, D-Acid and L-Base, feature hydrophobic heptad repeats and are designed to pack against each other in a "knobs-into-holes" manner. In solution, the peptides form a stable, helical heterotetramer with tight packing in the most solvent-protected core. This motif may be useful for designing protease-resistant, helical D-peptide ligands against biological protein targets.

Two structures of cyclophilin 40: folding and fidelity in the TPR domains

BACKGROUND

The "large immunophilin" family consists of domains of cyclophilin or FK506 binding protein linked to a tetratricopeptide (TPR) domain. They are intimately associated with steroid receptor complexes and bind to the C-terminal domain of Hsp90 via the TPR domain. The competitive binding of specific large immunophilins and other TPR-Hsp90 proteins provides a regulatory mechanism for Hsp90 chaperone activity.

RESULTS

We have solved the X-ray structures of monoclinic and tetragonal forms of Cyp40. In the monoclinic form, the TPR domain consists of seven helices of variable length incorporating three TPR motifs, which provide a convincing binding surface for the Hsp90 C-terminal MEEVD sequence. The C-terminal residues of Cyp40 protrude out beyond the body of the TPR domain to form a charged helix-the putative calmodulin binding site. However, in the tetragonal form, two of the TPR helices have straightened out to form one extended helix, providing a dramatically different conformation of the molecule.

CONCLUSIONS

The X-ray structures are consistent with the role of Cyclophilin 40 as a multifunctional signaling protein involved in a variety of protein-protein interactions. The intermolecular helix-helix interactions in the tetragonal form mimic the intramolecular interactions found in the fully folded monoclinic form. These conserved intra- and intermolecular TPR-TPR interactions are illustrative of a high-fidelity recognition mechanism. The two structures also open up the possibility that partially folded forms of TPR may be important in domain swapping and protein recognition.

Review: protein design--where we were, where we are, where we're going

Protein design has become a powerful approach for understanding the relationship between amino acid sequence and 3-dimensional structure. In the past 5 years, there have been many breakthroughs in the development of computational methods that allow the selection of novel sequences given the structure of a protein backbone. Successful design of protein scaffolds has now paved the way for new endeavors to design function. The ability to design sequences compatible with a fold may also be useful in structural and functional genomics by expanding the range of proteins used for fold recognition and for the identification of functionally important domains from multiple sequence alignments.

Design of three-dimensional domain-swapped dimers and fibrous oligomers

Three-dimensional (3D) domain-swapped proteins are intermolecularly folded analogs of monomeric proteins; both are stabilized by the identical interactions, but the individual domains interact intramolecularly in monomeric proteins, whereas they form intermolecular interactions in 3D domain-swapped structures. The structures and conditions of formation of several domain-swapped dimers and trimers are known, but the formation of higher order 3D domain-swapped oligomers has been less thoroughly studied. Here we contrast the structural consequences of domain swapping from two designed three-helix bundles: one with an up-down-up topology, and the other with an up-down-down topology. The up-down-up topology gives rise to a domain-swapped dimer whose structure has been determined to 1.5 A resolution by x-ray crystallography. In contrast, the domain-swapped protein with an up-down-down topology forms fibrils as shown by electron microscopy and dynamic light scattering. This demonstrates that design principles can predict the oligomeric state of 3D domain-swapped molecules, which should aid in the design of domain-swapped proteins and biomaterials.

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.

Retrostructural analysis of metalloproteins: application to the design of a minimal model for diiron proteins

De novo protein design provides an attractive approach for the construction of models to probe the features required for function of complex metalloproteins. The metal-binding sites of many metalloproteins lie between multiple elements of secondary structure, inviting a retrostructural approach to constructing minimal models of their active sites. The backbone geometries comprising the metal-binding sites of zinc fingers, diiron proteins, and rubredoxins may be described to within approximately 1 A rms deviation by using a simple geometric model with only six adjustable parameters. These geometric models provide excellent starting points for the design of metalloproteins, as illustrated in the construction of Due Ferro 1 (DF1), a minimal model for the Glu-Xxx-Xxx-His class of dinuclear metalloproteins. This protein was synthesized and structurally characterized as the di-Zn(II) complex by x-ray crystallography, by using data that extend to 2.5 A. This four-helix bundle protein is comprised of two noncovalently associated helix-loop-helix motifs. The dinuclear center is formed by two bridging Glu and two chelating Glu side chains, as well as two monodentate His ligands. The primary ligands are mostly buried in the protein interior, and their geometries are stabilized by a network of hydrogen bonds to second-shell ligands. In particular, a Tyr residue forms a hydrogen bond to a chelating Glu ligand, similar to a motif found in the diiron-containing R2 subunit of Escherichia coli ribonucleotide reductase and the ferritins. DF1 also binds cobalt and iron ions and should provide an attractive model for a variety of diiron proteins that use oxygen for processes including iron storage, radical formation, and hydrocarbon oxidation.