Preclinical proof of principle for orally delivered Th17 antagonist miniproteins.

Interleukin (IL)-23 and IL-17 are well-validated therapeutic targets in autoinflammatory diseases. Antibodies targeting IL-23 and IL-17 have shown clinical efficacy but are limited by high costs, safety risks, lack of sustained efficacy, and poor patient convenience as they require parenteral administration. Here, we present designed miniproteins inhibiting IL-23R and IL-17 with antibody-like, low picomolar affinities at a fraction of the molecular size. The minibinders potently block cell signaling in vitro and are extremely stable, enabling oral administration and low-cost manufacturing. The orally administered IL-23R minibinder shows efficacy better than a clinical anti-IL-23 antibody in mouse colitis and has a favorable pharmacokinetics (PK) and biodistribution profile in rats. This work demonstrates that orally administered de novo-designed minibinders can reach a therapeutic target past the gut epithelial barrier. With high potency, gut stability, and straightforward manufacturability, de novo-designed minibinders are a promising modality for oral biologics.

Accurate de novo design of membrane-traversing macrocycles.

We use computational design coupled with experimental characterization to systematically investigate the design principles for macrocycle membrane permeability and oral bioavailability. We designed 184 6-12 residue macrocycles with a wide range of predicted structures containing noncanonical backbone modifications and experimentally determined structures of 35; 29 are very close to the computational models. With such control, we show that membrane permeability can be systematically achieved by ensuring all amide (NH) groups are engaged in internal hydrogen bonding interactions. 84 designs over the 6-12 residue size range cross membranes with an apparent permeability greater than 1 × 10-6 cm/s. Designs with exposed NH groups can be made membrane permeable through the design of an alternative isoenergetic fully hydrogen-bonded state favored in the lipid membrane. The ability to robustly design membrane-permeable and orally bioavailable peptides with high structural accuracy should contribute to the next generation of designed macrocycle therapeutics.

Blueprinting extendable nanomaterials with standardized protein blocks.

A wooden house frame consists of many different lumber pieces, but because of the regularity of these building blocks, the structure can be designed using straightforward geometrical principles. The design of multicomponent protein assemblies, in comparison, has been much more complex, largely owing to the irregular shapes of protein structures1. Here we describe extendable linear, curved and angled protein building blocks, as well as inter-block interactions, that conform to specified geometric standards; assemblies designed using these blocks inherit their extendability and regular interaction surfaces, enabling them to be expanded or contracted by varying the number of modules, and reinforced with secondary struts. Using X-ray crystallography and electron microscopy, we validate nanomaterial designs ranging from simple polygonal and circular oligomers that can be concentrically nested, up to large polyhedral nanocages and unbounded straight 'train track' assemblies with reconfigurable sizes and geometries that can be readily blueprinted. Because of the complexity of protein structures and sequence-structure relationships, it has not previously been possible to build up large protein assemblies by deliberate placement of protein backbones onto a blank three-dimensional canvas; the simplicity and geometric regularity of our design platform now enables construction of protein nanomaterials according to 'back of an envelope' architectural blueprints.

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.

Design of protein-binding proteins from the target structure alone.

The design of proteins that bind to a specific site on the surface of a target protein using no information other than the three-dimensional structure of the target remains a challenge1-5. Here we describe a general solution to this problem that starts with a broad exploration of the vast space of possible binding modes to a selected region of a protein surface, and then intensifies the search in the vicinity of the most promising binding modes. We demonstrate the broad applicability of this approach through the de novo design of binding proteins to 12 diverse protein targets with different shapes and surface properties. Biophysical characterization shows that the binders, which are all smaller than 65 amino acids, are hyperstable and, following experimental optimization, bind their targets with nanomolar to picomolar affinities. We succeeded in solving crystal structures of five of the binder-target complexes, and all five closely match the corresponding computational design models. Experimental data on nearly half a million computational designs and hundreds of thousands of point mutants provide detailed feedback on the strengths and limitations of the method and of our current understanding of protein-protein interactions, and should guide improvements of both. Our approach enables the targeted design of binders to sites of interest on a wide variety of proteins for therapeutic and diagnostic applications.

De novo protein design by deep network hallucination.

There has been considerable recent progress in protein structure prediction using deep neural networks to predict inter-residue distances from amino acid sequences1-3. Here we investigate whether the information captured by such networks is sufficiently rich to generate new folded proteins with sequences unrelated to those of the naturally occurring proteins used in training the models. We generate random amino acid sequences, and input them into the trRosetta structure prediction network to predict starting residue-residue distance maps, which, as expected, are quite featureless. We then carry out Monte Carlo sampling in amino acid sequence space, optimizing the contrast (Kullback-Leibler divergence) between the inter-residue distance distributions predicted by the network and background distributions averaged over all proteins. Optimization from different random starting points resulted in novel proteins spanning a wide range of sequences and predicted structures. We obtained synthetic genes encoding 129 of the network-'hallucinated' sequences, and expressed and purified the proteins in Escherichia coli; 27 of the proteins yielded monodisperse species with circular dichroism spectra consistent with the hallucinated structures. We determined the three-dimensional structures of three of the hallucinated proteins, two by X-ray crystallography and one by NMR, and these closely matched the hallucinated models. Thus, deep networks trained to predict native protein structures from their sequences can be inverted to design new proteins, and such networks and methods should contribute alongside traditional physics-based models to the de novo design of proteins with new functions.

Rapid and automated design of two-component protein nanomaterials using ProteinMPNN.

The design of protein-protein interfaces using physics-based design methods such as Rosetta requires substantial computational resources and manual refinement by expert structural biologists. Deep learning methods promise to simplify protein-protein interface design and enable its application to a wide variety of problems by researchers from various scientific disciplines. Here, we test the ability of a deep learning method for protein sequence design, ProteinMPNN, to design two-component tetrahedral protein nanomaterials and benchmark its performance against Rosetta. ProteinMPNN had a similar success rate to Rosetta, yielding 13 new experimentally confirmed assemblies, but required orders of magnitude less computation and no manual refinement. The interfaces designed by ProteinMPNN were substantially more polar than those designed by Rosetta, which facilitated in vitro assembly of the designed nanomaterials from independently purified components. Crystal structures of several of the assemblies confirmed the accuracy of the design method at high resolution. Our results showcase the potential of deep learning-based methods to unlock the widespread application of designed protein-protein interfaces and self-assembling protein nanomaterials in biotechnology.

Design of amyloidogenic peptide traps.

Segments of proteins with high ß-strand propensity can self-associate to form amyloid fibrils implicated in many diseases. We describe a general approach to bind such segments in ß-strand and ß-hairpin conformations using de novo designed scaffolds that contain deep peptide-binding clefts. The designs bind their cognate peptides in vitro with nanomolar affinities. The crystal structure of a designed protein-peptide complex is close to the design model, and NMR characterization reveals how the peptide-binding cleft is protected in the apo state. We use the approach to design binders to the amyloid-forming proteins transthyretin, tau, serum amyloid A1 and amyloid ß1-42 (Aß42). The Aß binders block the assembly of Aß fibrils as effectively as the most potent of the clinically tested antibodies to date and protect cells from toxic Aß42 species.

De novo design of proteins housing excitonically coupled chlorophyll special pairs.

Natural photosystems couple light harvesting to charge separation using a 'special pair' of chlorophyll molecules that accepts excitation energy from the antenna and initiates an electron-transfer cascade. To investigate the photophysics of special pairs independently of the complexities of native photosynthetic proteins, and as a first step toward creating synthetic photosystems for new energy conversion technologies, we designed C2-symmetric proteins that hold two chlorophyll molecules in closely juxtaposed arrangements. X-ray crystallography confirmed that one designed protein binds two chlorophylls in the same orientation as native special pairs, whereas a second designed protein positions them in a previously unseen geometry. Spectroscopy revealed that the chlorophylls are excitonically coupled, and fluorescence lifetime imaging demonstrated energy transfer. The cryo-electron microscopy structure of a designed 24-chlorophyll octahedral nanocage with a special pair on each edge closely matched the design model. The results suggest that the de novo design of artificial photosynthetic systems is within reach of current computational methods.

De novo design of small beta barrel proteins.

Small beta barrel proteins are attractive targets for computational design because of their considerable functional diversity despite their very small size (<70 amino acids). However, there are considerable challenges to designing such structures, and there has been little success thus far. Because of the small size, the hydrophobic core stabilizing the fold is necessarily very small, and the conformational strain of barrel closure can oppose folding; also intermolecular aggregation through free beta strand edges can compete with proper monomer folding. Here, we explore the de novo design of small beta barrel topologies using both Rosetta energy-based methods and deep learning approaches to design four small beta barrel folds: Src homology 3 (SH3) and oligonucleotide/oligosaccharide-binding (OB) topologies found in nature and five and six up-and-down-stranded barrels rarely if ever seen in nature. Both approaches yielded successful designs with high thermal stability and experimentally determined structures with less than 2.4 Å rmsd from the designed models. Using deep learning for backbone generation and Rosetta for sequence design yielded higher design success rates and increased structural diversity than Rosetta alone. The ability to design a large and structurally diverse set of small beta barrel proteins greatly increases the protein shape space available for designing binders to protein targets of interest.

De novo designed ice-binding proteins from twist-constrained helices.

Attaining molecular-level control over solidification processes is a crucial aspect of materials science. To control ice formation, organisms have evolved bewildering arrays of ice-binding proteins (IBPs), but these have poorly understood structure-activity relationships. We propose that reverse engineering using de novo computational protein design can shed light on structure-activity relationships of IBPs. We hypothesized that the model alpha-helical winter flounder antifreeze protein uses an unusual undertwisting of its alpha-helix to align its putative ice-binding threonine residues in exactly the same direction. We test this hypothesis by designing a series of straight three-helix bundles with an ice-binding helix projecting threonines and two supporting helices constraining the twist of the ice-binding helix. Our findings show that ice-recrystallization inhibition by the designed proteins increases with the degree of designed undertwisting, thus validating our hypothesis, and opening up avenues for the computational design of IBPs.

Improving the secretion of designed protein assemblies through negative design of cryptic transmembrane domains.

Computationally designed protein nanoparticles have recently emerged as a promising platform for the development of new vaccines and biologics. For many applications, secretion of designed nanoparticles from eukaryotic cells would be advantageous, but in practice, they often secrete poorly. Here we show that designed hydrophobic interfaces that drive nanoparticle assembly are often predicted to form cryptic transmembrane domains, suggesting that interaction with the membrane insertion machinery could limit efficient secretion. We develop a general computational protocol, the Degreaser, to design away cryptic transmembrane domains without sacrificing protein stability. The retroactive application of the Degreaser to previously designed nanoparticle components and nanoparticles considerably improves secretion, and modular integration of the Degreaser into design pipelines results in new nanoparticles that secrete as robustly as naturally occurring protein assemblies. Both the Degreaser protocol and the nanoparticles we describe may be broadly useful in biotechnological applications.

De novo design of protein homodimers containing tunable symmetric protein pockets.

Function follows form in biology, and the binding of small molecules requires proteins with pockets that match the shape of the ligand. For design of binding to symmetric ligands, protein homo-oligomers with matching symmetry are advantageous as each protein subunit can make identical interactions with the ligand. Here, we describe a general approach to designing hyperstable C2 symmetric proteins with pockets of diverse size and shape. We first designed repeat proteins that sample a continuum of curvatures but have low helical rise, then docked these into C2 symmetric homodimers to generate an extensive range of C2 symmetric cavities. We used this approach to design thousands of C2 symmetric homodimers, and characterized 101 of them experimentally. Of these, the geometry of 31 were confirmed by small angle X-ray scattering and 2 were shown by crystallographic analyses to be in close agreement with the computational design models. These scaffolds provide a rich set of starting points for binding a wide range of C2 symmetric compounds.

Improving Protein Expression, Stability, and Function with ProteinMPNN.

Natural proteins are highly optimized for function but are often difficult to produce at a scale suitable for biotechnological applications due to poor expression in heterologous systems, limited solubility, and sensitivity to temperature. Thus, a general method that improves the physical properties of native proteins while maintaining function could have wide utility for protein-based technologies. Here, we show that the deep neural network ProteinMPNN, together with evolutionary and structural information, provides a route to increasing protein expression, stability, and function. For both myoglobin and tobacco etch virus (TEV) protease, we generated designs with improved expression, elevated melting temperatures, and improved function. For TEV protease, we identified multiple designs with improved catalytic activity as compared to the parent sequence and previously reported TEV variants. Our approach should be broadly useful for improving the expression, stability, and function of biotechnologically important proteins.

Accurate computational design of three-dimensional protein crystals.

Protein crystallization plays a central role in structural biology. Despite this, the process of crystallization remains poorly understood and highly empirical, with crystal contacts, lattice packing arrangements and space group preferences being largely unpredictable. Programming protein crystallization through precisely engineered side-chain-side-chain interactions across protein-protein interfaces is an outstanding challenge. Here we develop a general computational approach for designing three-dimensional protein crystals with prespecified lattice architectures at atomic accuracy that hierarchically constrains the overall number of degrees of freedom of the system. We design three pairs of oligomers that can be individually purified, and upon mixing, spontaneously self-assemble into >100 µm three-dimensional crystals. The structures of these crystals are nearly identical to the computational design models, closely corresponding in both overall architecture and the specific protein-protein interactions. The dimensions of the crystal unit cell can be systematically redesigned while retaining the space group symmetry and overall architecture, and the crystals are extremely porous and highly stable. Our approach enables the computational design of protein crystals with high accuracy, and the designed protein crystals, which have both structural and assembly information encoded in their primary sequences, provide a powerful platform for biological materials engineering.

The nonreceptor tyrosine kinase SRMS inhibits autophagy and promotes tumor growth by phosphorylating the scaffolding protein FKBP51.

Nutrient-responsive protein kinases control the balance between anabolic growth and catabolic processes such as autophagy. Aberrant regulation of these kinases is a major cause of human disease. We report here that the vertebrate nonreceptor tyrosine kinase Src-related kinase lacking C-terminal regulatory tyrosine and N-terminal myristylation sites (SRMS) inhibits autophagy and promotes growth in a nutrient-responsive manner. Under nutrient-replete conditions, SRMS phosphorylates the PHLPP scaffold FK506-binding protein 51 (FKBP51), disrupts the FKBP51-PHLPP complex, and promotes FKBP51 degradation through the ubiquitin-proteasome pathway. This prevents PHLPP-mediated dephosphorylation of AKT, causing sustained AKT activation that promotes growth and inhibits autophagy. SRMS is amplified and overexpressed in human cancers where it drives unrestrained AKT signaling in a kinase-dependent manner. SRMS kinase inhibition activates autophagy, inhibits cancer growth, and can be accomplished using the FDA-approved tyrosine kinase inhibitor ibrutinib. This illuminates SRMS as a targetable vulnerability in human cancers and as a new target for pharmacological induction of autophagy in vertebrates.

Design of Diverse Asymmetric Pockets in De Novo Homo-oligomeric Proteins.

A challenge for design of protein-small-molecule recognition is that incorporation of cavities with size, shape, and composition suitable for specific recognition can considerably destabilize protein monomers. This challenge can be overcome through binding pockets formed at homo-oligomeric interfaces between folded monomers. Interfaces surrounding the central homo-oligomer symmetry axes necessarily have the same symmetry and so may not be well suited to binding asymmetric molecules. To enable general recognition of arbitrary asymmetric substrates and small molecules, we developed an approach to designing asymmetric interfaces at off-axis sites on homo-oligomers, analogous to those found in native homo-oligomeric proteins such as glutamine synthetase. We symmetrically dock curved helical repeat proteins such that they form pockets at the asymmetric interface of the oligomer with sizes ranging from several angstroms, appropriate for binding a single ion, to up to more than 20 Å across. Of the 133 proteins tested, 84 had soluble expression in E. coli, 47 had correct oligomeric states in solution, 35 had small-angle X-ray scattering (SAXS) data largely consistent with design models, and 8 had negative-stain electron microscopy (nsEM) 2D class averages showing the structures coming together as designed. Both an X-ray crystal structure and a cryogenic electron microscopy (cryoEM) structure are close to the computational design models. The nature of these proteins as homo-oligomers allows them to be readily built into higher-order structures such as nanocages, and the asymmetric pockets of these structures open rich possibilities for small-molecule binder design free from the constraints associated with monomer destabilization.

An enumerative algorithm for de novo design of proteins with diverse pocket structures.

To create new enzymes and biosensors from scratch, precise control over the structure of small-molecule binding sites is of paramount importance, but systematically designing arbitrary protein pocket shapes and sizes remains an outstanding challenge. Using the NTF2-like structural superfamily as a model system, we developed an enumerative algorithm for creating a virtually unlimited number of de novo proteins supporting diverse pocket structures. The enumerative algorithm was tested and refined through feedback from two rounds of large-scale experimental testing, involving in total the assembly of synthetic genes encoding 7,896 designs and assessment of their stability on yeast cell surface, detailed biophysical characterization of 64 designs, and crystal structures of 5 designs. The refined algorithm generates proteins that remain folded at high temperatures and exhibit more pocket diversity than naturally occurring NTF2-like proteins. We expect this approach to transform the design of small-molecule sensors and enzymes by enabling the creation of binding and active site geometries much more optimal for specific design challenges than is accessible by repurposing the limited number of naturally occurring NTF2-like proteins.

Structural basis of the atypical activation mechanism of KRASV14I.

RAS regulation and signaling are largely accomplished by direct protein-protein interactions, making RAS protein dynamics a critical determinant of RAS function. Here, we report a crystal structure of GDP-bound KRASV14I, a mutated KRAS variant associated with the developmental RASopathy disorder Noonan syndrome (NS), at 1.5-1.6 Å resolution. The structure is notable for revealing a marked extension of switch 1 away from the G-domain and nucleotide-binding site of the KRAS protein. We found that this extension is associated with a loss of the magnesium ion and a tilt in the position of the guanine base because of the additional carbon introduced by the isoleucine substitution. Hydrogen-deuterium exchange MS analysis confirmed that this conformation occurs in solution, but also disclosed a difference in kinetics when compared with KRASA146T, another RAS mutant that displays a nearly identical conformation in previously reported crystal structures. This conformational change contributed to a high rate of guanine nucleotide-exchange factor (GEF)-dependent and -independent nucleotide exchange and to an increase in affinity for SOS Ras/Rac GEF 1 (SOS1), which appears to be the major mode of activation for this RAS variant. These results highlight a mechanistic connection between KRASA146T and KRASV14I that may have implications for the regulation of these variants and for the development of therapeutic strategies to manage KRAS variant-associated disorders.

KRAS Switch Mutants D33E and A59G Crystallize in the State 1 Conformation.

KRAS switch loop movements play a crucial role in regulating RAS signaling, and alteration of these sensitive dynamics is a principal mechanism through which disease-associated RAS mutations lead to aberrant RAS activation. Prior studies suggest that despite a high degree of sequence similarity, the switches in KRAS are more dynamic than those in HRAS. We determined X-ray crystal structures of the rare tumorigenic KRAS mutants KRASD33E, in switch 1 (SW1), and KRASA59G, in switch 2 (SW2), bound to GDP and found these adopt nearly identical, open SW1 conformations as well as altered SW2 conformations. KRASA59G bound to a GTP analogue crystallizes in the same conformation. This open conformation is consistent with the inactive "state 1" previously observed for HRAS bound to GTP. For KRASA59G, switch rearrangements may be regulated by increased flexibility in the 57DXXGQ61 motif at codon 59. However, loss of interactions between side chains at codons 33 and 35 in the SW1 33DPT35 motif drives changes for KRASD33E. The 33DPT35 motif is conserved for multiple members of the RAS subfamily but is not found in RAB, RHO, ARF, or Gα families, suggesting that dynamics mediated by this motif may be important for determining the selectivity of RAS-effector interactions. Biochemically, the consequence of altered switch dynamics is the same, showing impaired interaction with the guanine exchange factor SOS and loss of GAP-dependent GTPase activity. However, interactions with the RBD of RAF are preserved. Overall, these observations add to a body of evidence suggesting that HRAS and KRAS show meaningful differences in functionality stemming from differential protein dynamics independent of the hypervariable region.