Elicitation of Potent Neutralizing Antibody Responses by Designed Protein Nanoparticle Vaccines for SARS-CoV-2

A safe, effective, and scalable vaccine is needed to halt the ongoing SARS-CoV-2 pandemic. We describe the structure-based design of self-assembling protein nanoparticle immunogens that elicit potent and protective antibody responses against SARS-CoV-2 in mice. The nanoparticle vaccines display 60 SARS-CoV-2 spike receptor-binding domains (RBDs) in a highly immunogenic array and induce neutralizing antibody titers 10-fold higher than the prefusion-stabilized spike despite a 5-fold lower dose. Antibodies elicited by the RBD nanoparticles target multiple distinct epitopes, suggesting they may not be easily susceptible to escape mutations, and exhibit a lower binding:neutralizing ratio than convalescent human sera, which may minimize the risk of vaccine-associated enhanced respiratory disease. The high yield and stability of the assembled nanoparticles suggest that manufacture of the nanoparticle vaccines will be highly scalable. These results highlight the utility of robust antigen display platforms and have launched cGMP manufacturing efforts to advance the SARS-CoV-2-RBD nanoparticle vaccine into the clinic.

Induction of Potent Neutralizing Antibody Responses by a Designed Protein Nanoparticle Vaccine for Respiratory Syncytial Virus

Respiratory syncytial virus (RSV) is a worldwide public health concern for which no vaccine is available. Elucidation of the prefusion structure of the RSV F glycoprotein and its identification as the main target of neutralizing antibodies have provided new opportunities for development of an effective vaccine. Here, we describe the structure-based design of a self-assembling protein nanoparticle presenting a prefusion-stabilized variant of the F glycoprotein trimer (DS-Cav1) in a repetitive array on the nanoparticle exterior. The two-component nature of the nanoparticle scaffold enabled the production of highly ordered, monodisperse immunogens that display DS-Cav1 at controllable density. In mice and nonhuman primates, the full-valency nanoparticle immunogen displaying 20 DS-Cav1 trimers induced neutralizing antibody responses ∼10-fold higher than trimeric DS-Cav1. These results motivate continued development of this promising nanoparticle RSV vaccine candidate and establish computationally designed two-component nanoparticles as a robust and customizable platform for structure-based vaccine design.

Computational protein design enables a novel one-carbon assimilation pathway

We describe a computationally designed enzyme, formolase (FLS), which catalyzes the carboligation of three one-carbon formaldehyde molecules into one three-carbon dihydroxyacetone molecule. The existence of FLS enables the design of a new carbon fixation pathway, the formolase pathway, consisting of a small number of thermodynamically favorable chemical transformations that convert formate into a three-carbon sugar in central metabolism. The formolase pathway is predicted to use carbon more efficiently and with less backward flux than any naturally occurring one-carbon assimilation pathway. When supplemented with enzymes carrying out the other steps in the pathway, FLS converts formate into dihydroxyacetone phosphate and other central metabolites in vitro. These results demonstrate how modern protein engineering and design tools can facilitate the construction of a completely new biosynthetic pathway.

Principles of Protein Stability and Their Application in Computational Design

Proteins are increasingly used in basic and applied biomedical research. Many proteins, however, are only marginally stable and can be expressed in limited amounts, thus hampering research and applications. Research has revealed the thermodynamic, cellular, and evolutionary principles and mechanisms that underlie marginal stability. With this growing understanding, computational stability design methods have advanced over the past two decades starting from methods that selectively addressed only some aspects of marginal stability. Current methods are more general and, by combining phylogenetic analysis with atomistic design, have shown drastic improvements in solubility, thermal stability, and aggregation resistance while maintaining the protein's primary molecular activity. Stability design is opening the way to rational engineering of improved enzymes, therapeutics, and vaccines and to the application of protein design methodology to large proteins and molecular activities that have proven challenging in the past.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Recent advances in de novo protein design: Principles, methods, and applications

The computational de novo protein design is increasingly applied to address a number of key challenges in biomedicine and biological engineering. Successes in expanding applications are driven by advances in design principles and methods over several decades. Here, we review recent innovations in major aspects of the de novo protein design and include how these advances were informed by principles of protein architecture and interactions derived from the wealth of structures in the Protein Data Bank. We describe developments in de novo generation of designable backbone structures, optimization of sequences, design scoring functions, and the design of the function. The advances not only highlight design goals reachable now but also point to the challenges and opportunities for the future of the field.

Alteration of enzyme specificity by computational loop remodeling and design

Altering the specificity of an enzyme requires precise positioning of side-chain functional groups that interact with the modified groups of the new substrate. This requires not only sequence changes that introduce the new functional groups but also sequence changes that remodel the structure of the protein backbone so that the functional groups are properly positioned. We describe a computational design method for introducing specific enzyme-substrate interactions by directed remodeling of loops near the active site. Benchmark tests on 8 native protein-ligand complexes show that the method can recover native loop lengths and, often, native loop conformations. We then use the method to redesign a critical loop in human guanine deaminase such that a key side-chain interaction is made with the substrate ammelide. The redesigned enzyme is 100-fold more active on ammelide and 2.5e4-fold less active on guanine than wild-type enzyme: The net change in specificity is 2.5e6-fold. The structure of the designed protein was confirmed by X-ray crystallographic analysis: The remodeled loop adopts a conformation that is within 1-A Calpha RMSD of the computational model.

Perturbing the energy landscape for improved packing during computational protein design

The FastDesign protocol in the molecular modeling program Rosetta iterates between sequence optimization and structure refinement to stabilize de novo designed protein structures and complexes. FastDesign has been used previously to design novel protein folds and assemblies with important applications in research and medicine. To promote sampling of alternative conformations and sequences, FastDesign includes stages where the energy landscape is smoothened by reducing repulsive forces. Here, we discover that this process disfavors larger amino acids in the protein core because the protein compresses in the early stages of refinement. By testing alternative ramping strategies for the repulsive weight, we arrive at a scheme that produces lower energy designs with more native-like sequence composition in the protein core. We further validate the protocol by designing and experimentally characterizing over 4000 proteins and show that the new protocol produces higher stability proteins.

Computational design of affinity and specificity at protein-protein interfaces

The computer-based design of protein-protein interactions is a rigorous test of our understanding of molecular recognition and an attractive approach for creating novel tools for cell and molecular research. Considerable attention has been placed on redesigning the affinity and specificity of naturally occurring interactions. Several studies have shown that reducing the desolvation costs for binding while preserving shape complimentarity and hydrogen bonding is an effective strategy for improving binding affinities. In favorable cases specificity has been designed by focusing only on interactions with the target protein, while in cases with closely related off-target proteins it has been necessary to explicitly disfavor unwanted binding partners. The rational design of protein-protein interactions from scratch is still an unsolved problem, but recent developments in flexible backbone design and energy functions hold promise for the future.

Computational design of a self-assembling symmetrical β-propeller protein

The modular structure of many protein families, such as β-propeller proteins, strongly implies that duplication played an important role in their evolution, leading to highly symmetrical intermediate forms. Previous attempts to create perfectly symmetrical propeller proteins have failed, however. We have therefore developed a new and rapid computational approach to design such proteins. As a test case, we have created a sixfold symmetrical β-propeller protein and experimentally validated the structure using X-ray crystallography. Each blade consists of 42 residues. Proteins carrying 2-10 identical blades were also expressed and purified. Two or three tandem blades assemble to recreate the highly stable sixfold symmetrical architecture, consistent with the duplication and fusion theory. The other proteins produce different monodisperse complexes, up to 42 blades (180 kDa) in size, which self-assemble according to simple symmetry rules. Our procedure is suitable for creating nano-building blocks from different protein templates of desired symmetry.

High-resolution design of a protein loop

Despite having irregular structure, protein loops often adopt specific conformations that are critical to protein function. Most studies in de novo protein design have focused on creating proteins with regular elements of secondary structure connected by very short loops or turns. To design longer protein loops that adopt specific conformations, we have developed a protocol within the Rosetta molecular modeling program that iterates between optimizing the sequence and conformation of a loop in search of low-energy sequence-structure pairs. We have tested the procedure by designing 10-residue loops for the connection between the second and third strand in the beta-sandwich protein tenascin. Three low-energy designs from 7,200 flexible backbone trajectories were selected for experimental characterization. All three designs, called LoopA, LoopB, and LoopC, adopt stable folded structures. High-resolution crystal structures of LoopA and LoopB have been solved. LoopB adopts a structure very similar to the design model (0.46 A rmsd), and all but one of the side chains are modeled in the correct rotamers. LoopA crystallized at low pH in a structure that differs dramatically from our design model. It forms a strand-swapped dimer mediated by hydrogen bonds to protonated glutamic acids. Gel filtration indicates that the protein is not a dimer at neutral pH. These results suggest that the high-resolution design of protein loops is possible; however, they also highlight how small changes in protein energetics can dramatically perturb the low free energy structure of a protein.

Combining computational and experimental screening for rapid optimization of protein properties

We present a combined computational and experimental method for the rapid optimization of proteins. Using beta-lactamase as a test case, we redesigned the active site region using our Protein Design Automation technology as a computational screen to search the entire sequence space. By eliminating sequences incompatible with the protein fold, Protein Design Automation rapidly reduced the number of sequences to a size amenable to experimental screening, resulting in a library of approximately equal 200,000 mutants. These were then constructed and experimentally screened to select for variants with improved resistance to the antibiotic cefotaxime. In a single round, we obtained variants exhibiting a 1,280-fold increase in resistance. To our knowledge, all of the mutations were novel, i.e., they have not been identified as beneficial by random mutagenesis or DNA shuffling or seen in any of the naturally occurring TEM beta-lactamases, the most prevalent type of Gram-negative beta-lactamases. This combined approach allows for the rapid improvement of any property that can be screened experimentally and provides a powerful broadly applicable tool for protein engineering.

Structure-based protocol for identifying mutations that enhance protein-protein binding affinities

The ability to manipulate protein binding affinities is important for the development of proteins as biosensors, industrial reagents, and therapeutics. We have developed a structure-based method to rationally predict single mutations at protein-protein interfaces that enhance binding affinities. The protocol is based on the premise that increasing buried hydrophobic surface area and/or reducing buried hydrophilic surface area will generally lead to enhanced affinity if large steric clashes are not introduced and buried polar groups are not left without a hydrogen bond partner. The procedure selects affinity enhancing point mutations at the protein-protein interface using three criteria: (1) the mutation must be from a polar amino acid to a non-polar amino acid or from a non-polar amino acid to a larger non-polar amino acid, (2) the free energy of binding as calculated with the Rosetta protein modeling program should be more favorable than the free energy of binding calculated for the wild-type complex and (3) the mutation should not be predicted to significantly destabilize the monomers. The performance of the computational protocol was experimentally tested on two separate protein complexes; Galpha(i1) from the heterotrimeric G-protein system bound to the RGS14 GoLoco motif, and the E2, UbcH7, bound to the E3, E6AP from the ubiquitin pathway. Twelve single-site mutations that were predicted to be stabilizing were synthesized and characterized in the laboratory. Nine of the 12 mutations successfully increased binding affinity with five of these increasing binding by over 1.0 kcal/mol. To further assess our approach we searched the literature for point mutations that pass our criteria and have experimentally determined binding affinities. Of the eight mutations identified, five were accurately predicted to increase binding affinity, further validating the method as a useful tool to increase protein-protein binding affinities.

Discovery of Novel Gain-of-Function Mutations Guided by Structure-Based Deep Learning

Despite the promise of deep learning accelerated protein engineering, examples of such improved proteins are scarce. Here we report that a 3D convolutional neural network trained to associate amino acids with neighboring chemical microenvironments can guide identification of novel gain-of-function mutations that are not predicted by energetics-based approaches. Amalgamation of these mutations improved protein function in vivo across three diverse proteins by at least 5-fold. Furthermore, this model provides a means to interrogate the chemical space within protein microenvironments and identify specific chemical interactions that contribute to the gain-of-function phenotypes resulting from individual mutations.

Generation of longer emission wavelength red fluorescent proteins using computationally designed libraries

The longer emission wavelengths of red fluorescent proteins (RFPs) make them attractive for whole-animal imaging because cells are more transparent to red light. Although several useful RFPs have been developed using directed evolution, the quest for further red-shifted and improved RFPs continues. Herein, we report a structure-based rational design approach to red-shift the fluorescence emission of RFPs. We applied a combined computational and experimental approach that uses computational protein design as an in silico prescreen to generate focused combinatorial libraries of mCherry mutants. The computational procedure helped us identify residues that could fulfill interactions hypothesized to cause red-shifts without destabilizing the protein fold. These interactions include stabilization of the excited state through H-bonding to the acylimine oxygen atom, destabilization of the ground state by hydrophobic packing around the charged phenolate, and stabilization of the excited state by a π-stacking interaction. Our methodology allowed us to identify three mCherry mutants (mRojoA, mRojoB, and mRouge) that display emission wavelengths > 630 nm, representing red-shifts of 20-26 nm. Moreover, our approach required the experimental screening of a total of ∼5,000 clones, a number several orders of magnitude smaller than those previously used to achieve comparable red-shifts. Additionally, crystal structures of mRojoA and mRouge allowed us to verify fulfillment of the interactions hypothesized to cause red-shifts, supporting their contribution to the observed red-shifts.

AbDesign: An algorithm for combinatorial backbone design guided by natural conformations and sequences

Computational design of protein function has made substantial progress, generating new enzymes, binders, inhibitors, and nanomaterials not previously seen in nature. However, the ability to design new protein backbones for function--essential to exert control over all polypeptide degrees of freedom--remains a critical challenge. Most previous attempts to design new backbones computed the mainchain from scratch. Here, instead, we describe a combinatorial backbone and sequence optimization algorithm called AbDesign, which leverages the large number of sequences and experimentally determined molecular structures of antibodies to construct new antibody models, dock them against target surfaces and optimize their sequence and backbone conformation for high stability and binding affinity. We used the algorithm to produce antibody designs that target the same molecular surfaces as nine natural, high-affinity antibodies; in five cases interface sequence identity is above 30%, and in four of those the backbone conformation at the core of the antibody binding surface is within 1 Å root-mean square deviation from the natural antibodies. Designs recapitulate polar interaction networks observed in natural complexes, and amino acid sidechain rigidity at the designed binding surface, which is likely important for affinity and specificity, is high compared to previous design studies. In designed anti-lysozyme antibodies, complementarity-determining regions (CDRs) at the periphery of the interface, such as L1 and H2, show greater backbone conformation diversity than the CDRs at the core of the interface, and increase the binding surface area compared to the natural antibody, potentially enhancing affinity and specificity.

Engineering Transcriptional Regulator Effector Specificity Using Computational Design and In Vitro Rapid Prototyping: Developing a Vanillin Sensor

The pursuit of circuits and metabolic pathways of increasing complexity and robustness in synthetic biology will require engineering new regulatory tools. Feedback control based on relevant molecules, including toxic intermediates and environmental signals, would enable genetic circuits to react appropriately to changing conditions. In this work, variants of qacR, a tetR family repressor, were generated by computational protein design and screened in a cell-free transcription-translation (TX-TL) system for responsiveness to a new targeted effector. The modified repressors target vanillin, a growth-inhibiting small molecule found in lignocellulosic hydrolysates and other industrial processes. Promising candidates from the in vitro screen were further characterized in vitro and in vivo in a gene circuit. The screen yielded two qacR mutants that respond to vanillin both in vitro and in vivo. While the mutants exhibit some toxicity to cells, presumably due to off-target effects, they are prime starting points for directed evolution toward vanillin sensors with the specifications required for use in a dynamic control loop. We believe this process, a combination of the generation of variants coupled with in vitro screening, can serve as a framework for designing new sensors for other target compounds.

Computational design of a PAK1 binding protein

We describe a computational protocol, called DDMI, for redesigning scaffold proteins to bind to a specified region on a target protein. The DDMI protocol is implemented within the Rosetta molecular modeling program and uses rigid-body docking, sequence design, and gradient-based minimization of backbone and side-chain torsion angles to design low-energy interfaces between the scaffold and target protein. Iterative rounds of sequence design and conformational optimization were needed to produce models that have calculated binding energies that are similar to binding energies calculated for native complexes. We also show that additional conformation sampling with molecular dynamics can be iterated with sequence design to further lower the computed energy of the designed complexes. To experimentally test the DDMI protocol, we redesigned the human hyperplastic discs protein to bind to the kinase domain of p21-activated kinase 1 (PAK1). Six designs were experimentally characterized. Two of the designs aggregated and were not characterized further. Of the remaining four designs, three bound to the PAK1 with affinities tighter than 350 muM. The tightest binding design, named Spider Roll, bound with an affinity of 100 muM. NMR-based structure prediction of Spider Roll based on backbone and (13)C(beta) chemical shifts using the program CS-ROSETTA indicated that the architecture of human hyperplastic discs protein is preserved. Mutagenesis studies confirmed that Spider Roll binds the target patch on PAK1. Additionally, Spider Roll binds to full-length PAK1 in its activated state but does not bind PAK1 when it forms an auto-inhibited conformation that blocks the Spider Roll target site. Subsequent NMR characterization of the binding of Spider Roll to PAK1 revealed a comparably small binding 'on-rate' constant (<<10(5) M(-1) s(-1)). The ability to rationally design the site of novel protein-protein interactions is an important step towards creating new proteins that are useful as therapeutics or molecular probes.

Computationally designed high specificity inhibitors delineate the roles of BCL2 family proteins in cancer

Many cancers overexpress one or more of the six human pro-survival BCL2 family proteins to evade apoptosis. To determine which BCL2 protein or proteins block apoptosis in different cancers, we computationally designed three-helix bundle protein inhibitors specific for each BCL2 pro-survival protein. Following in vitro optimization, each inhibitor binds its target with high picomolar to low nanomolar affinity and at least 300-fold specificity. Expression of the designed inhibitors in human cancer cell lines revealed unique dependencies on BCL2 proteins for survival which could not be inferred from other BCL2 profiling methods. Our results show that designed inhibitors can be generated for each member of a closely-knit protein family to probe the importance of specific protein-protein interactions in complex biological processes.

Development of High Affinity and High Specificity Inhibitors of Matrix Metalloproteinase 14 through Computational Design and Directed Evolution

Degradation of the extracellular matrices in the human body is controlled by matrix metalloproteinases (MMPs), a family of more than 20 homologous enzymes. Imbalance in MMP activity can result in many diseases, such as arthritis, cardiovascular diseases, neurological disorders, fibrosis, and cancers. Thus, MMPs present attractive targets for drug design and have been a focus for inhibitor design for as long as 3 decades. Yet, to date, all MMP inhibitors have failed in clinical trials because of their broad activity against numerous MMP family members and the serious side effects of the proposed treatment. In this study, we integrated a computational method and a yeast surface display technique to obtain highly specific inhibitors of MMP-14 by modifying the natural non-specific broad MMP inhibitor protein N-TIMP2 to interact optimally with MMP-14. We identified an N-TIMP2 mutant, with five mutations in its interface, that has an MMP-14 inhibition constant (K_i ) of 0.9 pm, the strongest MMP-14 inhibitor reported so far. Compared with wild-type N-TIMP2, this variant displays ∼900-fold improved affinity toward MMP-14 and up to 16,000-fold greater specificity toward MMP-14 relative to other MMPs. In an in vitro and cell-based model of MMP-dependent breast cancer cellular invasiveness, this N-TIMP2 mutant acted as a functional inhibitor. Thus, our study demonstrates the enormous potential of a combined computational/directed evolution approach to protein engineering. Furthermore, it offers fundamental clues into the molecular basis of MMP regulation by N-TIMP2 and identifies a promising MMP-14 inhibitor as a starting point for the development of protein-based anticancer therapeutics.

Rapid search for tertiary fragments reveals protein sequence-structure relationships

Finding backbone substructures from the Protein Data Bank that match an arbitrary query structural motif, composed of multiple disjoint segments, is a problem of growing relevance in structure prediction and protein design. Although numerous protein structure search approaches have been proposed, methods that address this specific task without additional restrictions and on practical time scales are generally lacking. Here, we propose a solution, dubbed MASTER, that is both rapid, enabling searches over the Protein Data Bank in a matter of seconds, and provably correct, finding all matches below a user-specified root-mean-square deviation cutoff. We show that despite the potentially exponential time complexity of the problem, running times in practice are modest even for queries with many segments. The ability to explore naturally plausible structural and sequence variations around a given motif has the potential to synthesize its design principles in an automated manner; so we go on to illustrate the utility of MASTER to protein structural biology. We demonstrate its capacity to rapidly establish structure-sequence relationships, uncover the native designability landscapes of tertiary structural motifs, identify structural signatures of binding, and automatically rewire protein topologies. Given the broad utility of protein tertiary fragment searches, we hope that providing MASTER in an open-source format will enable novel advances in understanding, predicting, and designing protein structure.

Antibody humanization by structure-based computational protein design

Antibodies derived from non-human sources must be modified for therapeutic use so as to mitigate undesirable immune responses. While complementarity-determining region (CDR) grafting-based humanization techniques have been successfully applied in many cases, it remains challenging to maintain the desired stability and antigen binding affinity upon grafting. We developed an alternative humanization approach called CoDAH ("Computationally-Driven Antibody Humanization") in which computational protein design methods directly select sets of amino acids to incorporate from human germline sequences to increase humanness while maintaining structural stability. Retrospective studies show that CoDAH is able to identify variants deemed beneficial according to both humanness and structural stability criteria, even for targets lacking crystal structures. Prospective application to TZ47, a murine anti-human B7H6 antibody, demonstrates the approach. Four diverse humanized variants were designed, and all possible unique VH/VL combinations were produced as full-length IgG1 antibodies. Soluble and cell surface expressed antigen binding assays showed that 75% (6 of 8) of the computationally designed VH/VL variants were successfully expressed and competed with the murine TZ47 for binding to B7H6 antigen. Furthermore, 4 of the 6 bound with an estimated KD within an order of magnitude of the original TZ47 antibody. In contrast, a traditional CDR-grafted variant could not be expressed. These results suggest that the computational protein design approach described here can be used to efficiently generate functional humanized antibodies and provide humanized templates for further affinity maturation.

Increasing sequence diversity with flexible backbone protein design: the complete redesign of a protein hydrophobic core

Protein design tests our understanding of protein stability and structure. Successful design methods should allow the exploration of sequence space not found in nature. However, when redesigning naturally occurring protein structures, most fixed backbone design algorithms return amino acid sequences that share strong sequence identity with wild-type sequences, especially in the protein core. This behavior places a restriction on functional space that can be explored and is not consistent with observations from nature, where sequences of low identity have similar structures. Here, we allow backbone flexibility during design to mutate every position in the core (38 residues) of a four-helix bundle protein. Only small perturbations to the backbone, 1-2 Å, were needed to entirely mutate the core. The redesigned protein, DRNN, is exceptionally stable (melting point >140°C). An NMR and X-ray crystal structure show that the side chains and backbone were accurately modeled (all-atom RMSD = 1.3 Å).

Shaping up the protein folding funnel by local interaction: lesson from a structure prediction study

Predicting protein tertiary structure by folding-like simulations is one of the most stringent tests of how much we understand the principle of protein folding. Currently, the most successful method for folding-based structure prediction is the fragment assembly (FA) method. Here, we address why the FA method is so successful and its lesson for the folding problem. To do so, using the FA method, we designed a structure prediction test of "chimera proteins." In the chimera proteins, local structural preference is specific to the target sequences, whereas nonlocal interactions are only sequence-independent compaction forces. We find that these chimera proteins can find the native folds of the intact sequences with high probability indicating dominant roles of the local interactions. We further explore roles of local structural preference by exact calculation of the HP lattice model of proteins. From these results, we suggest principles of protein folding: For small proteins, compact structures that are fully compatible with local structural preference are few, one of which is the native fold. These local biases shape up the funnel-like energy landscape.

Computational design of a leucine-rich repeat protein with a predefined geometry

Structure-based protein design offers a possibility of optimizing the overall shape of engineered binding scaffolds to match their targets better. We developed a computational approach for the structure-based design of repeat proteins that allows for adjustment of geometrical features like length, curvature, and helical twist. By combining sequence optimization of existing repeats and de novo design of capping structures, we designed leucine-rich repeats (LRRs) from the ribonuclease inhibitor (RI) family that assemble into structures with a predefined geometry. The repeat proteins were built from self-compatible LRRs that are designed to interact to form highly curved and planar assemblies. We validated the geometrical design approach by engineering a ring structure constructed from 10 self-compatible repeats. Protein design can also be used to increase our structural understanding of repeat proteins. We use our design constructs to demonstrate that buried Cys play a central role for stability and folding cooperativity in RI-type LRR proteins. The computational procedure presented here may be used to develop repeat proteins with various geometrical shapes for applications where greater control of the interface geometry is desired.