Rapid search for tertiary fragments reveals protein sequence-structure relationships

A Designed Zn2+ Sensor Domain Transmits Binding Information to Transmembrane Histidine Kinases

Generating stimulus-responsive allosteric signaling de novo is a significant challenge in protein design. In natural systems like bacterial histidine kinases (HKs), signal transduction occurs when ligand binding initiates a signal that is amplified across biological membranes over long distances to induce large-scale rearrangements and phosphorylation relays. Here, we ask whether our understanding of protein design and multidomain, intramolecular signaling has progressed sufficiently to enable engineering of a HK with tunable de novo components. We generated de novo metal-binding sensor domains and substituted them for the native sensor domain of a transmembrane HK, affording chimeras that transduce signals initiated from a de novo sensor. Signaling depended on the designed sensor's stability and the interdomain linker's phase and length. These results show the usefulness of de novo design to elucidate the biochemical mechanisms and principles of transmembrane signaling.

Surfactant-like peptide gels are based on cross-β amyloid fibrils

Surfactant-like peptides, in which hydrophilic and hydrophobic residues are encoded within different domains in the peptide sequence, undergo facile self-assembly in aqueous solution to form supramolecular hydrogels. These peptides have been explored extensively as substrates for the creation of functional materials since a wide variety of amphipathic sequences can be prepared from commonly available amino acid precursors. The self-assembly behavior of surfactant-like peptides has been compared to that observed for small molecule amphiphiles in which nanoscale phase separation of the hydrophobic domains drives the self-assembly of supramolecular structures. Here, we investigate the relationship between sequence and supramolecular structure for a pair of bola-amphiphilic peptides, Ac-KLIIIK-NH2 (L2) and Ac-KIIILK-NH2 (L5). Despite similar length, composition, and polar sequence pattern, L2 and L5 form morphologically distinct assemblies, nanosheets and nanotubes, respectively. Cryo-EM helical reconstruction was employed to determine the structure of the L5 nanotube at near-atomic resolution. Rather than displaying self-assembly behavior analogous to conventional amphiphiles, the packing arrangement of peptides in the L5 nanotube displayed steric zipper interfaces that resembled those observed in the structures of β-amyloid fibrils. Like amyloids, the supramolecular structures of the L2 and L5 assemblies were sensitive to conservative amino acid substitutions within an otherwise identical amphipathic sequence pattern. This study highlights the need to better understand the relationship between sequence and supramolecular structure to facilitate the development of functional peptide-based materials for biomaterials applications.

Designed peptide binders and nanobodies as PROTAC starting points for targeted degradation of PCNA and BCL6

The efficient degradation of pathogenic proteins, particularly proliferating cell nuclear antigen (PCNA) and B-cell lymphoma 6 protein (BCL6), is crucial for treating various diseases related to cancer. As key biological macromolecules, PCNA plays a critical role in DNA replication and repair, while BCL6 acts as a transcriptional repressor involved in B-cell lymphoma. To enhance the efficiency and specificity of protein degradation, we developed a RS80E-based bioPROTACs system that consists of truncated variants of Ring-B-boxed coiled-coil (RBCC) domains (RS80E) with improved degradation efficiency fused to an AI-driven binder/nanobody targeting specific antigens. Combining state-of-the-art methodologies such as ProteinMPNN, RFdiffusion, AlphaFold3, AlphaFold2, and HADDOCK, we identified binders for PCNA and predicted spatial interrelationships. Employing fragment-based and alanine scanning methods, we designed nanobodies targeting PCNA and BCL6 by combinatorially designing CDR3 and grafting them onto nanobody scaffolds. Significantly, our results demonstrate the utility of bioPROTACs in degrading PCNA and BCL6, thereby activating p53 and promoting apoptosis. This highlights the therapeutic potential of targeting PCNA and BCL6 degradation and lays the groundwork for developing PCNA and BCL6-degrading therapeutics. In summary, our system offers a modular and rapid pathway for exploration other intractable therapeutic targets, and emphasizes the importance of interdisciplinary methods in advancing therapeutic interventions.

Rational Design Principles for De Novo α-Helical Peptide Barrels with Dynamic Conductive Channels

Despite advances in peptide and protein design, the rational design of membrane-spanning peptides that form conducting channels remains challenging due to our imperfect understanding of the sequence-to-structure relationships that drive membrane insertion, assembly, and conductance. Here, we describe the design and computational and experimental characterization of a series of coiled coil-based peptides that form transmembrane α-helical barrels with conductive channels. Through a combination of rational and computational design, we obtain barrels with 5 to 7 helices, as characterized in detergent micelles. In lipid bilayers, these peptide assemblies exhibit two conductance states with relative populations dependent on the applied potential: (i) low-conductance states that correlate with variations in the designed amino-acid sequences and modeled coiled-coil barrel geometries, indicating stable transmembrane α-helical barrels; and (ii) high-conductance states in which single channels change size in discrete steps. Notably, the high-conductance states are similar for all peptides in contrast to the low-conductance states. This indicates the formation of large, dynamic channels, as observed in natural barrel-stave peptide channels. These findings establish rational routes to design and tune functional membrane-spanning peptide channels with specific conductance and geometry.

De Novo Design of Proteins That Bind Naphthalenediimides, Powerful Photooxidants with Tunable Photophysical Properties

De novo protein design provides a framework to test our understanding of protein function and build proteins with cofactors and functions not found in nature. Here, we report the design of proteins designed to bind powerful photooxidants and the evaluation of the use of these proteins to generate diffusible small-molecule reactive species. Because excited-state dynamics are influenced by the dynamics and hydration of a photooxidant's environment, it was important to not only design a binding site but also to evaluate its dynamic properties. Thus, we used computational design in conjunction with molecular dynamics (MD) simulations to design a protein, designated NBP (NDI Binding Protein), that held a naphthalenediimide (NDI), a powerful photooxidant, in a programmable molecular environment. Solution NMR confirmed the structure of the complex. We evaluated two NDI cofactors in this de novo protein using ultrafast pump-probe spectroscopy to evaluate light-triggered intra- and intermolecular electron transfer function. Moreover, we demonstrated the utility of this platform to activate multiple molecular probes for protein labeling.

Assessing Generative Model Coverage of Protein Structures with SHAPES

Recent advances in generative modeling enable efficient sampling of protein structures, but their tendency to optimize for designability imposes a bias toward idealized structures at the expense of loops and other complex structural motifs critical for function. We introduce SHAPES (Structural and Hierarchical Assessment of Proteins with Embedding Similarity) to evaluate five state-of-the-art generative models of protein structures. Using structural embeddings across multiple structural hierarchies, ranging from local geometries to global protein architectures, we reveal substantial undersampling of the observed protein structure space by these models. We use Fréchet Protein Distance (FPD) to quantify distributional coverage. Different models are distinct in their coverage behavior across different sampling noise scales and temperatures; the frequency of TERtiary Motifs (TERMs) further supports the observations. More robust sequence design and structure prediction methods are likely crucial in guiding the development of models with improved coverage of the designable protein space.

A designed Zn 2+ sensor domain transmits binding information to transmembrane histidine kinases

Generating stimulus-responsive, allosteric signaling de novo is a significant challenge in protein design. In natural systems like bacterial histidine kinases (HKs), signal transduction occurs when ligand binding initiates a signal that is amplified across biological membranes over long distances to induce large-scale rearrangements and phosphorylation relays. Here, we ask whether our understanding of protein design and multi-domain, intramolecular signaling has progressed sufficiently to enable engineering of a HK with tunable de novo components. We generated de novo metal-binding sensor domains and substituted them for the native sensor domain of a transmembrane HK, affording chimeras that transduce signals initiated from a de novo sensor. Signaling depended on the designed sensor's stability and the interdomain linker's phase and length. These results show the usefulness of de novo design to elucidate biochemical mechanisms and principles for design of new signaling systems.

Exploring "dark-matter" protein folds using deep learning

De novo protein design explores uncharted sequence and structure space to generate novel proteins not sampled by evolution. A main challenge in de novo design involves crafting "designable" structural templates to guide the sequence searches toward adopting target structures. We present a convolutional variational autoencoder that learns patterns of protein structure, dubbed Genesis. We coupled Genesis with trRosetta to design sequences for a set of protein folds and found that Genesis is capable of reconstructing native-like distance and angle distributions for five native folds and three novel, the so-called "dark-matter" folds as a demonstration of generalizability. We used a high-throughput assay to characterize the stability of the designs through protease resistance, obtaining encouraging success rates for folded proteins. Genesis enables exploration of the protein fold space within minutes, unrestricted by protein topologies. Our approach addresses the backbone designability problem, showing that small neural networks can efficiently learn structural patterns in proteins. A record of this paper's transparent peer review process is included in the supplemental information.

Unveiling the hidden interactome of CRBN molecular glues with chemoproteomics

Targeted protein degradation and induced proximity refer to strategies that leverage the recruitment of proteins to facilitate their modification, regulation or degradation. As prospective design of glues remains challenging, unbiased discovery methods are needed to unveil hidden chemical targets. Here we establish a high throughput affinity purification mass spectrometry workflow in cell lysates for the unbiased identification of molecular glue targets. By mapping the targets of 20 CRBN-binding molecular glues, we identify 298 protein targets and demonstrate the utility of enrichment methods for identifying novel targets overlooked using established methods. We use a computational workflow to estimate target confidence and perform a biochemical screen to identify a lead compound for the new non-ZF target PPIL4. Our study provides a comprehensive inventory of targets chemically recruited to CRBN and delivers a robust and scalable workflow for identifying new drug-induced protein interactions in cell lysates.

PB-GPT: An innovative GPT-based model for protein backbone generation

With advanced computational methods, it is now feasible to modify or design proteins for specific functions, a process with significant implications for disease treatment and other medical applications. Protein structures and functions are intrinsically linked to their backbones, making the design of these backbones a pivotal aspect of protein engineering. In this study, we focus on the task of unconditionally generating protein backbones. By means of codebook quantization and compression dictionaries, we convert protein backbone structures into a distinctive coded language and propose a GPT-based protein backbone generation model, PB-GPT. To validate the generalization performance of the model, we trained and evaluated the model on both public datasets and small protein datasets. The results demonstrate that our model has the capability to unconditionally generate elaborate, highly realistic protein backbones with structural patterns resembling those of natural proteins, thus showcasing the significant potential of large language models in protein structure design.

Protocol for Designing De Novo Noncanonical Peptide Binders in OSPREY

D-peptides, the mirror image of canonical L-peptides, offer numerous biological advantages that make them effective therapeutics. This article details how to use DexDesign, the newest OSPREY-based algorithm, for designing these D-peptides de novo. OSPREY physics-based models precisely mimic energy-equivariant reflection operations, enabling the generation of D-peptide scaffolds from L-peptide templates. Due to the scarcity of D-peptide:L-protein structural data, DexDesign calls a geometric hashing algorithm, Method of Accelerated Search for Tertiary Ensemble Representatives, as a subroutine to produce a synthetic structural dataset. DexDesign enables mixed-chirality designs with a new user interface and also reduces the conformation and sequence search space using three new design techniques: Minimum Flexible Set, Inverse Alanine Scanning, and K*-based Mutational Scanning.

Collagen-targeted protein nanomicelles for the imaging of non-alcoholic steatohepatitis

In vivo molecular imaging tools hold immense potential to drive transformative breakthroughs by enabling researchers to visualize cellular and molecular interactions in real-time and/or at high resolution. These advancements will facilitate a deeper understanding of fundamental biological processes and their dysregulation in disease states. Here, we develop and characterize a self-assembling protein nanomicelle called collagen type I binding - thermoresponsive assembled protein (Col1-TRAP) that binds tightly to type I collagen in vitro with nanomolar affinity. For ex vivo visualization, Col1-TRAP is labeled with a near-infrared fluorescent dye (NIR-Col1-TRAP). Both Col1-TRAP and NIR-Col1-TRAP display approximately a 3.8-fold greater binding to type I collagen compared to TRAP when measured by surface plasmon resonance (SPR). We present a proof-of-concept study using NIR-Col1-TRAP to detect fibrotic type I collagen deposition ex vivo in the livers of mice with non-alcoholic steatohepatitis (NASH). We show that NIR-Col1-TRAP demonstrates significantly decreased plasma recirculation time as well as increased liver accumulation in the NASH mice compared to mice without disease over 4 hours. As a result, NIR-Col1-TRAP shows potential as an imaging probe for NASH with in vivo targeting performance after injection in mice. STATEMENT OF SIGNIFICANCE: Direct molecular imaging of fibrosis in NASH patients enables the diagnosis and monitoring of disease progression with greater specificity and resolution than do elastography-based methods or blood tests. In addition, protein-based imaging probes are more advantageous than alternatives due to their biodegradability and scalable biosynthesis. With the aid of computational modeling, we have designed a self-assembled protein micelle that binds to fibrillar and monomeric collagen in vitro. After the protein was labeled with near-infrared fluorescent dye, we injected the compound into mice fed on a NASH diet. NIR-Col1-TRAP clears from the serum faster in these mice compared to control mice, and accumulates significantly more in fibrotic livers.This work advances the development of targeted protein probes for in vivo fibrosis imaging.

Rationally seeded computational protein design of ɑ-helical barrels

Computational protein design is advancing rapidly. Here we describe efficient routes starting from validated parallel and antiparallel peptide assemblies to design two families of α-helical barrel proteins with central channels that bind small molecules. Computational designs are seeded by the sequences and structures of defined de novo oligomeric barrel-forming peptides, and adjacent helices are connected by loop building. For targets with antiparallel helices, short loops are sufficient. However, targets with parallel helices require longer connectors; namely, an outer layer of helix-turn-helix-turn-helix motifs that are packed onto the barrels. Throughout these computational pipelines, residues that define open states of the barrels are maintained. This minimizes sequence sampling, accelerating the design process. For each of six targets, just two to six synthetic genes are made for expression in Escherichia coli. On average, 70% of these genes express to give soluble monomeric proteins that are fully characterized, including high-resolution structures for most targets that match the design models with high accuracy.

De novo design of peptides that bind specific conformers of α-synuclein

Insoluble amyloids rich in cross-β fibrils are observed in a number of neurodegenerative diseases. Depending on the clinicopathology, the amyloids can adopt distinct supramolecular assemblies, termed conformational strains. However, rapid methods to study amyloids in a conformationally specific manner are lacking. We introduce a novel computational method for de novo design of peptides that tile the surface of α-synuclein fibrils in a conformationally specific manner. Our method begins by identifying surfaces that are unique to the conformational strain of interest, which becomes a "target backbone" for the design of a peptide binder. Next, we interrogate structures in the PDB with high geometric complementarity to the target. Then, we identify secondary structural motifs that interact with this target backbone in a favorable, highly occurring geometry. This method produces monomeric helical motifs with a favorable geometry for interaction with the strands of the underlying amyloid. Each motif is then symmetrically replicated to form a monolayer that tiles the amyloid surface. Finally, amino acid sequences of the peptide binders are computed to provide a sequence with high geometric and physicochemical complementarity to the target amyloid. This method was applied to a conformational strain of α-synuclein fibrils, resulting in a peptide with high specificity for the target relative to other amyloids formed by α-synuclein, tau, or Aβ40. This designed peptide also markedly slowed the formation of α-synuclein amyloids. Overall, this method offers a new tool for examining conformational strains of amyloid proteins.

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.

DexDesign: A new OSPREY-based algorithm for designing de novo D-peptide inhibitors

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan, which enable exponential reductions in the size of the peptide sequence search space. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a new framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead therapeutic candidates. We provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

DexDesign: an OSPREY-based algorithm for designing de novo D-peptide inhibitors

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead inhibitors. We also provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small-molecule-binding proteins has seen exciting recent progress; however, the ability to achieve exquisite affinity for binding small molecules while tuning specificity has not yet been demonstrated directly from computation. Here, we develop a computational procedure that results in the highest affinity binders to date with predetermined relative affinities, targeting a series of PARP1 inhibitors. Two of four designed proteins bound with affinities ranging from < 5 nM to low μM, in a predictable manner. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free-energy calculations performed directly on the designed models are in excellent agreement with the experimentally measured affinities, suggesting that the de novo design of small-molecule-binding proteins with tuned interaction energies is now feasible entirely from computation. We expect these methods to open many opportunities in biomedicine, including rapid sensor development, antidote design, and drug delivery vehicles.

De novo Design of Peptides that Bind Specific Conformers of α-Synuclein

Insoluble amyloids rich in cross-β fibrils are observed in a number of neurodegenerative diseases. Depending on the clinicopathology, the amyloids can adopt distinct supramolecular assemblies, termed conformational strains. However, rapid methods to study amyloid in a conformationally specific manner are lacking. We introduce a novel computational method for de novo design of peptides that tile the surface of α-synuclein fibrils in a conformationally specific manner. Our method begins by identifying surfaces that are unique to the conformational strain of interest, which becomes a "target backbone" for the design of a peptide binder. Next, we interrogate structures in the PDB database with high geometric complementarity to the target. Then, we identify secondary structural motifs that interact with this target backbone in a favorable, highly occurring geometry. This method produces monomeric helical motifs with a favorable geometry for interaction with the strands of the underlying amyloid. Each motif is then symmetrically replicated to form a monolayer that tiles the amyloid surface. Finally, amino acid sequences of the peptide binders are computed to provide a sequence with high geometric and physicochemical complementarity to the target amyloid. This method was applied to a conformational strain of α-synuclein fibrils, resulting in a peptide with high specificity for the target relative to other amyloids formed by α-synuclein, tau, or Aβ40. This designed peptide also markedly slowed the formation of α-synuclein amyloids. Overall, this method offers a new tool for examining conformational strains of amyloid proteins.

A Suite of Designed Protein Cages Using Machine Learning Algorithms and Protein Fragment-Based Protocols

Designed protein cages and related materials provide unique opportunities for applications in biotechnology and medicine, while methods for their creation remain challenging and unpredictable. In the present study, we apply new computational approaches to design a suite of new tetrahedrally symmetric, self-assembling protein cages. For the generation of docked poses, we emphasize a protein fragment-based approach, while for de novo interface design, a comparison of computational protocols highlights the power and increased experimental success achieved using the machine learning program ProteinMPNN. In relating information from docking and design, we observe that agreement between fragment-based sequence preferences and ProteinMPNN sequence inference correlates with experimental success. Additional insights for designing polar interactions are highlighted by experimentally testing larger and more polar interfaces. In all, using X-ray crystallography and cryo-EM, we report five structures for seven protein cages, with atomic resolution in the best case reaching 2.0 Å. We also report structures of two incompletely assembled protein cages, providing unique insights into one type of assembly failure. The new set of designed cages and their structures add substantially to the body of available protein nanoparticles, and to methodologies for their creation.

De novo design of protein interactions with learned surface fingerprints

Physical interactions between proteins are essential for most biological processes governing life1. However, the molecular determinants of such interactions have been challenging to understand, even as genomic, proteomic and structural data increase. This knowledge gap has been a major obstacle for the comprehensive understanding of cellular protein-protein interaction networks and for the de novo design of protein binders that are crucial for synthetic biology and translational applications2-9. Here we use a geometric deep-learning framework operating on protein surfaces that generates fingerprints to describe geometric and chemical features that are critical to drive protein-protein interactions10. We hypothesized that these fingerprints capture the key aspects of molecular recognition that represent a new paradigm in the computational design of novel protein interactions. As a proof of principle, we computationally designed several de novo protein binders to engage four protein targets: SARS-CoV-2 spike, PD-1, PD-L1 and CTLA-4. Several designs were experimentally optimized, whereas others were generated purely in silico, reaching nanomolar affinity with structural and mutational characterization showing highly accurate predictions. Overall, our surface-centric approach captures the physical and chemical determinants of molecular recognition, enabling an approach for the de novo design of protein interactions and, more broadly, of artificial proteins with function.

Peptides from human BNIP5 and PXT1 and non-native binders of pro-apoptotic BAK can directly activate or inhibit BAK-mediated membrane permeabilization

Apoptosis is important for development and tissue homeostasis, and its dysregulation can lead to diseases, including cancer. As an apoptotic effector, BAK undergoes conformational changes that promote mitochondrial outer membrane disruption, leading to cell death. This is termed "activation" and can be induced by peptides from the human proteins BID, BIM, and PUMA. To identify additional peptides that can regulate BAK, we used computational protein design, yeast surface display screening, and structure-based energy scoring to identify 10 diverse new binders. We discovered peptides from the human proteins BNIP5 and PXT1 and three non-native peptides that activate BAK in liposome assays and induce cytochrome c release from mitochondria. Crystal structures and binding studies reveal a high degree of similarity among peptide activators and inhibitors, ruling out a simple function-determining property. Our results shed light on the vast peptide sequence space that can regulate BAK function and will guide the design of BAK-modulating tools and therapeutics.

PPI-Miner: A Structure and Sequence Motif Co-Driven Protein-Protein Interaction Mining and Modeling Computational Method

Protein-protein interactions (PPIs) play important roles in biological processes of life, and predicting PPIs becomes a critical scientific issue of concern. Most PPIs occur through small domains or motifs (fragments), which are challenging and laborious to map by standard biochemical approaches because they generally require the cloning of several truncation mutants. Here, we present a computational method, named as PPI-Miner, to fish potential protein interacting partners utilizing protein motifs as queries. In brief, this work first developed a motif-matching algorithm designed to identify the proteins that contain sequential or structural similar motifs with the given query motif. Being aligned to the query motif, the binding mode of the discovered motif and its receptor protein will be initially determined to be used to build PPI complexes accordingly. Eventually, a PPI complex structure could be built and optimized with a designed automatic protocol. Besides discovering PPIs, PPI-Miner can also be applied to other areas, i.e., the rational design of molecular glues and protein vaccines. In this work, PPI-Miner was employed to mine the potential cereblon (CRBN) substrates from human proteome. As a result, 1,739 candidates were predicted, and 16 of them have been experimentally validated in previous studies. The source code of PPI-Miner can be obtained from the GitHub repository (https://github.com/Wang-Lin-boop/PPI-Miner), the webserver is freely available for users (https://bailab.siais.shanghaitech.edu.cn/services/ppi-miner), and the database of predicted CRBN substrates is accessible at https://bailab.siais.shanghaitech.edu.cn/services/crbn-subslib.

Fragment-based computational design of antibodies targeting structured epitopes

De novo design methods hold the promise of reducing the time and cost of antibody discovery while enabling the facile and precise targeting of predetermined epitopes. Here, we describe a fragment-based method for the combinatorial design of antibody binding loops and their grafting onto antibody scaffolds. We designed and tested six single-domain antibodies targeting different epitopes on three antigens, including the receptor-binding domain of the SARS-CoV-2 spike protein. Biophysical characterization showed that all designs are stable and bind their intended targets with affinities in the nanomolar range without in vitro affinity maturation. We further discuss how a high-resolution input antigen structure is not required, as similar predictions are obtained when the input is a crystal structure or a computer-generated model. This computational procedure, which readily runs on a laptop, provides a starting point for the rapid generation of lead antibodies binding to preselected epitopes.

A generic framework for hierarchical de novo protein design

De novo protein design enables the exploration of novel sequences and structures absent from the natural protein universe. De novo design also stands as a stringent test for our understanding of the underlying physical principles of protein folding and may lead to the development of proteins with unmatched functional characteristics. The first fundamental challenge of de novo design is to devise "designable" structural templates leading to sequences that will adopt the predicted fold. Here, we built on the TopoBuilder (TB) de novo design method, to automatically assemble structural templates with native-like features starting from string descriptors that capture the overall topology of proteins. Our framework eliminates the dependency of hand-crafted and fold-specific rules through an iterative, data-driven approach that extracts geometrical parameters from structural tertiary motifs. We evaluated the TopoBuilder framework by designing sequences for a set of five protein folds and experimental characterization revealed that several sequences were folded and stable in solution. The TopoBuilder de novo design framework will be broadly useful to guide the generation of artificial proteins with customized geometries, enabling the exploration of the protein universe.

Catalysis and Electron Transfer in De Novo Designed Metalloproteins

One of the hallmark advances in our understanding of metalloprotein function is showcased in our ability to design new, non-native, catalytically active protein scaffolds. This review highlights progress and milestone achievements in the field of de novo metalloprotein design focused on reports from the past decade with special emphasis on de novo designs couched within common subfields of bioinorganic study: heme binding proteins, monometal- and dimetal-containing catalytic sites, and metal-containing electron transfer sites. Within each subfield, we highlight several of what we have identified as significant and important contributions to either our understanding of that subfield or de novo metalloprotein design as a discipline. These reports are placed in context both historically and scientifically. General suggestions for future directions that we feel will be important to advance our understanding or accelerate discovery are discussed.

Matching protein surface structural patches for high-resolution blind peptide docking

Modeling interactions between short peptides and their receptors is a challenging docking problem due to the peptide flexibility, resulting in a formidable sampling problem of peptide conformation in addition to its orientation. Alternatively, the peptide can be viewed as a piece that complements the receptor monomer structure. Here, we show that the peptide conformation can be determined based on the receptor backbone only and sampled using local structural motifs found in solved protein monomers and interfaces, independent of sequence similarity. This approach outperforms current peptide docking protocols and promotes new directions for peptide interface design.

Peptide docking can be perceived as a subproblem of protein–protein docking. However, due to the short length and flexible nature of peptides, many do not adopt one defined conformation prior to binding. Therefore, to tackle a peptide docking problem, not only the relative orientation, but also the bound conformation of the peptide needs to be modeled. Traditional peptide-centered approaches use information about peptide sequences to generate representative conformer ensembles, which can then be rigid-body docked to the receptor. Alternatively, one may look at this problem from the viewpoint of the receptor, namely, that the protein surface defines the peptide-bound conformation. Here, we present PatchMAN (Patch-Motif AligNments), a global peptide-docking approach that uses structural motifs to map the receptor surface with backbone scaffolds extracted from protein structures. On a nonredundant set of protein–peptide complexes, starting from free receptor structures, PatchMAN successfully models and identifies near-native peptide–protein complexes in 58%/84% within 2.5 Å/5 Å interface backbone RMSD, with corresponding sampling in 81%/100% of the cases, outperforming other approaches. PatchMAN leverages the observation that structural units of peptides with their binding pocket can be found not only within interfaces, but also within monomers. We show that the bound peptide conformation is sampled based on the structural context of the receptor only, without taking into account any sequence information. Beyond peptide docking, this approach opens exciting new avenues to study principles of peptide–protein association, and to the design of new peptide binders. PatchMAN is available as a server at https://furmanlab.cs.huji.ac.il/patchman/.

Structure‐conditioned amino‐acid couplings: how contact geometry affects pairwise sequence preferences

Relating a protein's sequence to its conformation is a central challenge for both structure prediction and sequence design. Statistical contact potentials, as well as their more descriptive versions that account for side‐chain orientation and other geometric descriptors, have served as simplistic but useful means of representing second‐order contributions in sequence–structure relationships. Here we ask what happens when a pairwise potential is conditioned on the fully defined geometry of interacting backbones fragments. We show that the resulting structure‐conditioned coupling energies more accurately reflect pair preferences as a function of structural contexts. These structure‐conditioned energies more reliably encode native sequence information and more highly correlate with experimentally determined coupling energies. Clustering a database of interaction motifs by structure results in ensembles of similar energies and clustering them by energy results in ensembles of similar structures. By comparing many pairs of interaction motifs and showing that structural similarity and energetic similarity go hand‐in‐hand, we provide a tangible link between modular sequence and structure elements. This link is applicable to structural modeling, and we show that scoring CASP models with structured‐conditioned energies results in substantially higher correlation with structural quality than scoring the same models with a contact potential. We conclude that structure‐conditioned coupling energies are a good way to model the impact of interaction geometry on second‐order sequence preferences.

De novo metalloprotein design

Natural metalloproteins perform many functions - ranging from sensing to electron transfer and catalysis - in which the position and property of each ligand and metal, is dictated by protein structure. De novo protein design aims to define an amino acid sequence that encodes a specific structure and function, providing a critical test of the hypothetical inner workings of (metallo)proteins. To date, de novo metalloproteins have used simple, symmetric tertiary structures - uncomplicated by the large size and evolutionary marks of natural proteins - to interrogate structure-function hypotheses. In this Review, we discuss de novo design applications, such as proteins that induce complex, increasingly asymmetric ligand geometries to achieve function, as well as the use of more canonical ligand geometries to achieve stability. De novo design has been used to explore how proteins fine-tune redox potentials and catalyse both oxidative and hydrolytic reactions. With an increased understanding of structure-function relationships, functional proteins including O2-dependent oxidases, fast hydrolases, and multi-proton/multi-electron reductases, have been created. In addition, proteins can now be designed using xeno-biological metals or cofactors and principles from inorganic chemistry to derive new-to-nature functions. These results and the advances in computational protein design suggest a bright future for the de novo design of diverse, functional metalloproteins.

Neuroserpin and transthyretin are extracellular chaperones that preferentially inhibit amyloid formation

Neuroserpin is a secreted protease inhibitor known to inhibit amyloid formation by the Alzheimer’s beta peptide (Aβ). To test whether this effect was constrained to Aβ, we used a range of in vitro assays to demonstrate that neuroserpin inhibits amyloid formation by several different proteins and protects against the associated cytotoxicity but, unlike other known chaperones, has a poor ability to inhibit amorphous protein aggregation. Collectively, these results suggest that neuroserpin has an unusual chaperone selectivity for intermediates on the amyloid-forming pathway. Bioinformatics analyses identified a highly conserved 14-residue region containing an α helix shared between neuroserpin and the thyroxine-transport protein transthyretin, and we subsequently demonstrated that transthyretin also preferentially inhibits amyloid formation. Last, we used rationally designed neuroserpin mutants to demonstrate a direct involvement of the conserved 14-mer region in its chaperone activity. Identification of this conserved region may prove useful in the future design of anti-amyloid reagents.

Deterministic chaos in the self-assembly of β sheet nanotubes from an amphipathic oligopeptide

The self-assembly of designed peptides into filaments and other higher-order structures has been the focus of intense interest because of the potential for creating new biomaterials and biomedical devices. These peptide assemblies have also been used as models for understanding biological processes, such as the pathological formation of amyloid. We investigate the assembly of an octapeptide sequence, Ac-FKFEFKFE-NH₂, motivated by prior studies that demonstrated that this amphipathic β strand peptide self-assembled into fibrils and biocompatible hydrogels. Using high-resolution cryoelectron microscopy (cryo-EM), we are able to determine the atomic structure for two different coexisting forms of the fibrils, containing four and five β sandwich protofilaments, respectively. Surprisingly, the inner walls in both forms are parallel β sheets, while the outer walls are antiparallel β sheets. Our results demonstrate the chaotic nature of peptide self-assembly and illustrate the importance of cryo-EM structural analysis to understand the complex phase behavior of these materials at near-atomic resolution.

Rossmann-toolbox: a deep learning-based protocol for the prediction and design of cofactor specificity in Rossmann fold proteins

The Rossmann fold enzymes are involved in essential biochemical pathways such as nucleotide and amino acid metabolism. Their functioning relies on interaction with cofactors, small nucleoside-based compounds specifically recognized by a conserved βαβ motif shared by all Rossmann fold proteins. While Rossmann methyltransferases recognize only a single cofactor type, the S-adenosylmethionine, the oxidoreductases, depending on the family, bind nicotinamide (nicotinamide adenine dinucleotide, nicotinamide adenine dinucleotide phosphate) or flavin-based (flavin adenine dinucleotide) cofactors. In this study, we showed that despite its short length, the βαβ motif unambiguously defines the specificity towards the cofactor. Following this observation, we trained two complementary deep learning models for the prediction of the cofactor specificity based on the sequence and structural features of the βαβ motif. A benchmark on two independent test sets, one containing βαβ motifs bearing no resemblance to those of the training set, and the other comprising 38 experimentally confirmed cases of rational design of the cofactor specificity, revealed the nearly perfect performance of the two methods. The Rossmann-toolbox protocols can be accessed via the webserver at https://lbs.cent.uw.edu.pl/rossmann-toolbox and are available as a Python package at https://github.com/labstructbioinf/rossmann-toolbox.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

A fragment-based protein interface design algorithm for symmetric assemblies

Theoretical and experimental advances in protein engineering have led to the creation of precisely defined, novel protein assemblies of great size and complexity, with diverse applications. One powerful approach involves designing a new attachment or binding interface between two simpler symmetric oligomeric protein components. The required methods of design, which present both similarities and key differences compared to problems in protein docking, remain challenging and are not yet routine. With the aim of more fully enabling this emerging area of protein material engineering, we developed a computer program, nanohedra, to introduce two key advances. First, we encoded in the program the construction rules (i.e. the search space parameters) that underlie all possible symmetric material constructions. Second, we developed algorithms for rapidly identifying favorable docking/interface arrangements based on tabulations of empirical patterns of known protein fragment-pair associations. As a result, the candidate poses that nanohedra generates for subsequent amino acid interface design appear highly native-like (at the protein backbone level), while simultaneously conforming to the exacting requirements for symmetry-based assembly. A retrospective computational analysis of successful vs failed experimental studies supports the expectation that this should improve the success rate for this challenging area of protein engineering.

Structural analysis of cross α-helical nanotubes provides insight into the designability of filamentous peptide nanomaterials

The exquisite structure-function correlations observed in filamentous protein assemblies provide a paradigm for the design of synthetic peptide-based nanomaterials. However, the plasticity of quaternary structure in sequence-space and the lability of helical symmetry present significant challenges to the de novo design and structural analysis of such filaments. Here, we describe a rational approach to design self-assembling peptide nanotubes based on controlling lateral interactions between protofilaments having an unusual cross-α supramolecular architecture. Near-atomic resolution cryo-EM structural analysis of seven designed nanotubes provides insight into the designability of interfaces within these synthetic peptide assemblies and identifies a non-native structural interaction based on a pair of arginine residues. This arginine clasp motif can robustly mediate cohesive interactions between protofilaments within the cross-α nanotubes. The structure of the resultant assemblies can be controlled through the sequence and length of the peptide subunits, which generates synthetic peptide filaments of similar dimensions to flagella and pili.

De Novo Design, Solution Characterization, and Crystallographic Structure of an Abiological Mn-Porphyrin-Binding Protein Capable of Stabilizing a Mn(V) Species

De novo protein design offers the opportunity to test our understanding of how metalloproteins perform difficult transformations. Attaining high-resolution structural information is critical to understanding how such designs function. There have been many successes in the design of porphyrin-binding proteins; however, crystallographic characterization has been elusive, limiting what can be learned from such studies as well as the extension to new functions. Moreover, formation of highly oxidizing high-valent intermediates poses design challenges that have not been previously implemented: (1) purposeful design of substrate/oxidant access to the binding site and (2) limiting deleterious oxidation of the protein scaffold. Here we report the first crystallographically characterized porphyrin-binding protein that was programmed to not only bind a synthetic Mn-porphyrin but also maintain binding site access to form high-valent oxidation states. We explicitly designed a binding site with accessibility to dioxygen units in the open coordination site of the Mn center. In solution, the protein is capable of accessing a high-valent Mn(V)-oxo species which can transfer an O atom to a thioether substrate. The crystallographic structure is within 0.6 Å of the design and indeed contained an aquo ligand with a second water molecule stabilized by hydrogen bonding to a Gln side chain in the active site, offering a structural explanation for the observed reactivity.

Allosteric cooperation in a de novo-designed two-domain protein

We describe the de novo design of an allosterically regulated protein, which comprises two tightly coupled domains. One domain is based on the DF (Due Ferri in Italian or two-iron in English) family of de novo proteins, which have a diiron cofactor that catalyzes a phenol oxidase reaction, while the second domain is based on PS1 (Porphyrin-binding Sequence), which binds a synthetic Zn-porphyrin (ZnP). The binding of ZnP to the original PS1 protein induces changes in structure and dynamics, which we expected to influence the catalytic rate of a fused DF domain when appropriately coupled. Both DF and PS1 are four-helix bundles, but they have distinct bundle architectures. To achieve tight coupling between the domains, they were connected by four helical linkers using a computational method to discover the most designable connections capable of spanning the two architectures. The resulting protein, DFP1 (Due Ferri Porphyrin), bound the two cofactors in the expected manner. The crystal structure of fully reconstituted DFP1 was also in excellent agreement with the design, and it showed the ZnP cofactor bound over 12 Å from the dimetal center. Next, a substrate-binding cleft leading to the diiron center was introduced into DFP1. The resulting protein acts as an allosterically modulated phenol oxidase. Its Michaelis-Menten parameters were strongly affected by the binding of ZnP, resulting in a fourfold tighter Km and a 7-fold decrease in kcat These studies establish the feasibility of designing allosterically regulated catalytic proteins, entirely from scratch.

Expanding the space of protein geometries by computational design of de novo fold families

Naturally occurring proteins vary the precise geometries of structural elements to create distinct shapes optimal for function. We present a computational design method, loop-helix-loop unit combinatorial sampling (LUCS), that mimics nature's ability to create families of proteins with the same overall fold but precisely tunable geometries. Through near-exhaustive sampling of loop-helix-loop elements, LUCS generates highly diverse geometries encompassing those found in nature but also surpassing known structure space. Biophysical characterization showed that 17 (38%) of 45 tested LUCS designs encompassing two different structural topologies were well folded, including 16 with designed non-native geometries. Four experimentally solved structures closely matched the designs. LUCS greatly expands the designable structure space and offers a new paradigm for designing proteins with tunable geometries that may be customizable for novel functions.

A defined structural unit enables de novo design of small-molecule-binding proteins

The de novo design of proteins that bind highly functionalized small molecules represents a great challenge. To enable computational design of binders, we developed a unit of protein structure-a van der Mer (vdM)-that maps the backbone of each amino acid to statistically preferred positions of interacting chemical groups. Using vdMs, we designed six de novo proteins to bind the drug apixaban; two bound with low and submicromolar affinity. X-ray crystallography and mutagenesis confirmed a structure with a precisely designed cavity that forms favorable interactions in the drug-protein complex. vdMs may enable design of functional proteins for applications in sensing, medicine, and catalysis.

Restriction of S-adenosylmethionine conformational freedom by knotted protein binding sites

S-adenosylmethionine (SAM) is one of the most important enzyme substrates. It is vital for the function of various proteins, including large group of methyltransferases (MTs). Intriguingly, some bacterial and eukaryotic MTs, while catalysing the same reaction, possess significantly different topologies, with the former being a knotted one. Here, we conducted a comprehensive analysis of SAM conformational space and factors that affect its vastness. We investigated SAM in two forms: free in water (via NMR studies and explicit solvent simulations) and bound to proteins (based on all data available in the PDB and on all-atom molecular dynamics simulations in water). We identified structural descriptors—angles which show the major differences in SAM conformation between unknotted and knotted methyltransferases. Moreover, we report that this is caused mainly by a characteristic for knotted MTs compact binding site formed by the knot and the presence of adenine-binding loop. Additionally, we elucidate conformational restrictions imposed on SAM molecules by other protein groups in comparison to conformational space in water.

The topology of a folded polypeptide chain has great impact on the resulting protein function and its interaction with ligands. Interestingly, topological constraints appear to affect binding of one of the most ubiquitous substrates in the cell, S-adenosylmethionine (SAM), to its target proteins. Here, we demonstrate how binding sites of specific proteins restrict SAM conformational freedom in comparison to its unbound state, with a special interest in proteins with non-trivial topology, including an exciting group of knotted methyltransferases. Using a vast array of computational methods combined with NMR experiments, we identify key structural features of knotted methyltransferases that impose unorthodox SAM conformations. We compare them with the characteristics of standard, unknotted SAM binding proteins. These results are significant for understanding differences between analogous, yet topologically different enzymes, as well as for future rational drug design.

De novo protein design enables the precise induction of RSV-neutralizing antibodies

De novo protein design has been successful in expanding the natural protein repertoire. However, most de novo proteins lack biological function, presenting a major methodological challenge. In vaccinology, the induction of precise antibody responses remains a cornerstone for next-generation vaccines. Here, we present a protein design algorithm called TopoBuilder, with which we engineered epitope-focused immunogens displaying complex structural motifs. In both mice and nonhuman primates, cocktails of three de novo-designed immunogens induced robust neutralizing responses against the respiratory syncytial virus. Furthermore, the immunogens refocused preexisting antibody responses toward defined neutralization epitopes. Overall, our design approach opens the possibility of targeting specific epitopes for the development of vaccines and therapeutic antibodies and, more generally, will be applicable to the design of de novo proteins displaying complex functional motifs.

ALEPH: a network-oriented approach for the generation of fragment-based libraries and for structure interpretation

The analysis of large structural databases reveals general features and relationships among proteins, providing useful insight. A different approach is required to characterize ubiquitous secondary-structure elements, where flexibility is essential in order to capture small local differences. The ALEPH software is optimized for the analysis and the extraction of small protein folds by relying on their geometry rather than on their sequence. The annotation of the structural variability of a given fold provides valuable information for fragment-based molecular-replacement methods, in which testing alternative model hypotheses can succeed in solving difficult structures when no homology models are available or are successful. ARCIMBOLDO_BORGES combines the use of composite secondary-structure elements as a search model with density modification and tracing to reveal the rest of the structure when both steps are successful. This phasing method relies on general fold libraries describing variations around a given pattern of β-sheets and helices extracted using ALEPH. The program introduces characteristic vectors defined from the main-chain atoms as a way to describe the geometrical properties of the structure. ALEPH encodes structural properties in a graph network, the exploration of which allows secondary-structure annotation, decomposition of a structure into small compact folds, generation of libraries of models representing a variation of a given fold and finally superposition of these folds onto a target structure. These functions are available through a graphical interface designed to interactively show the results of structure manipulation, annotation, fold decomposition, clustering and library generation. ALEPH can produce pictures of the graphs, structures and folds for publication purposes.

A general-purpose protein design framework based on mining sequence-structure relationships in known protein structures

Current state-of-the-art approaches to computational protein design (CPD) aim to capture the determinants of structure from physical principles. While this has led to many successful designs, it does have strong limitations associated with inaccuracies in physical modeling, such that a reliable general solution to CPD has yet to be found. Here, we propose a design framework-one based on identifying and applying patterns of sequence-structure compatibility found in known proteins, rather than approximating them from models of interatomic interactions. We carry out extensive computational analyses and an experimental validation for our method. Our results strongly argue that the Protein Data Bank is now sufficiently large to enable proteins to be designed by using only examples of structural motifs from unrelated proteins. Because our method is likely to have orthogonal strengths relative to existing techniques, it could represent an important step toward removing remaining barriers to robust CPD.

Dissociating STAT4 and STAT5 Signaling Inhibitory Functions of SOCS3: Effects on CD8 T Cell Responses

Cytokines are critical for guiding the differentiation of T lymphocytes to perform specialized tasks in the immune response. Developing strategies to manipulate cytokine-signaling pathways holds promise to program T cell differentiation toward the most therapeutically useful direction. Suppressor of cytokine signaling (SOCS) proteins are attractive targets, as they effectively inhibit undesirable cytokine signaling. However, these proteins target multiple signaling pathways, some of which we may need to remain uninhibited. SOCS3 inhibits IL-12 signaling but also inhibits the IL-2–signaling pathway. In this study, we use computational protein design based on SOCS3 and JAK crystal structures to engineer a mutant SOCS3 with altered specificity. We generated a mutant SOCS3 designed to ablate interactions with JAK1 but maintain interactions with JAK2. We show that this mutant does indeed ablate JAK1 inhibition, although, unexpectedly, it still coimmunoprecipitates with JAK1 and does so to a greater extent than with JAK2. When expressed in CD8 T cells, mutant SOCS3 preserved inhibition of JAK2-dependent STAT4 phosphorylation following IL-12 treatment. However, inhibition of STAT phosphorylation was ablated following stimulation with JAK1-dependent cytokines IL-2, IFN-α, and IL-21. Wild-type SOCS3 inhibited CD8 T cell expansion in vivo and induced a memory precursor phenotype. In vivo T cell expansion was restored by expression of the mutant SOCS3, and this also reverted the phenotype toward effector T cell differentiation. These data show that SOCS proteins can be engineered to fine-tune their specificity, and this can exert important changes to T cell biology.

Protein Design by Provable Algorithms

Protein design algorithms can leverage provable guarantees of accuracy to provide new insights and unique optimized molecules.

Packing of apolar side chains enables accurate design of highly stable membrane proteins

The features that stabilize the structures of membrane proteins remain poorly understood. Polar interactions contribute modestly, and the hydrophobic effect contributes little to the energetics of apolar side-chain packing in membranes. Disruption of steric packing can destabilize the native folds of membrane proteins, but is packing alone sufficient to drive folding in lipids? If so, then membrane proteins stabilized by this feature should be readily designed and structurally characterized-yet this has not been achieved. Through simulation of the natural protein phospholamban and redesign of variants, we define a steric packing code underlying its assembly. Synthetic membrane proteins designed using this code and stabilized entirely by apolar side chains conform to the intended fold. Although highly stable, the steric complementarity required for their folding is surprisingly stringent. Structural informatics shows that the designed packing motif recurs across the proteome, emphasizing a prominent role for precise apolar packing in membrane protein folding, stabilization, and evolution.

X-ray Crystal Structure of the Influenza A M2 Proton Channel S31N Mutant in Two Conformational States: An Open and Shut Case

The amantadine-resistant S31N mutant of the influenza A M2 proton channel has become prevalent in currently circulating viruses. Here, we have solved an X-ray crystal structure of M2(22-46) S31N that contains two distinct conformational states within its asymmetric unit. This structure reveals the mechanism of adamantane resistance in both conformational states of the M2 channel. In the Inward_open conformation, the mutant Asn31 side chain faces the channel pore and sterically blocks the adamantane binding site. In the Inward_closed conformation, Asn31 forms hydrogen bonds with carbonyls at the monomer-monomer interface, which twists the monomer helices and constricts the channel pore at the drug binding site. We also examine M2(19-49) WT and S31N using solution NMR spectroscopy and show that distribution of the two conformational states is dependent on both detergent choice and experimental pH.

Rosetta FunFolDes - A general framework for the computational design of functional proteins

The robust computational design of functional proteins has the potential to deeply impact translational research and broaden our understanding of the determinants of protein function and stability. The low success rates of computational design protocols and the extensive in vitro optimization often required, highlight the challenge of designing proteins that perform essential biochemical functions, such as binding or catalysis. One of the most simplistic approaches for the design of function is to adopt functional motifs in naturally occurring proteins and transplant them to computationally designed proteins. The structural complexity of the functional motif largely determines how readily one can find host protein structures that are "designable", meaning that are likely to present the functional motif in the desired conformation. One promising route to enhance the "designability" of protein structures is to allow backbone flexibility. Here, we present a computational approach that couples conformational folding with sequence design to embed functional motifs into heterologous proteins-Rosetta Functional Folding and Design (FunFolDes). We performed extensive computational benchmarks, where we observed that the enforcement of functional requirements resulted in designs distant from the global energetic minimum of the protein. An observation consistent with several experimental studies that have revealed function-stability tradeoffs. To test the design capabilities of FunFolDes we transplanted two viral epitopes into distant structural templates including one de novo "functionless" fold, which represent two typical challenges where the designability problem arises. The designed proteins were experimentally characterized showing high binding affinities to monoclonal antibodies, making them valuable candidates for vaccine design endeavors. Overall, we present an accessible strategy to repurpose old protein folds for new functions. This may lead to important improvements on the computational design of proteins, with structurally complex functional sites, that can perform elaborate biochemical functions related to binding and catalysis.

CATS (Coordinates of Atoms by Taylor Series): protein design with backbone flexibility in all locally feasible directions

Motivation

When proteins mutate or bind to ligands, their backbones often move significantly, especially in loop regions. Computational protein design algorithms must model these motions in order to accurately optimize protein stability and binding affinity. However, methods for backbone conformational search in design have been much more limited than for sidechain conformational search. This is especially true for combinatorial protein design algorithms, which aim to search a large sequence space efficiently and thus cannot rely on temporal simulation of each candidate sequence.

Results

We alleviate this difficulty with a new parameterization of backbone conformational space, which represents all degrees of freedom of a specified segment of protein chain that maintain valid bonding geometry (by maintaining the original bond lengths and angles and ω dihedrals). In order to search this space, we present an efficient algorithm, CATS, for computing atomic coordinates as a function of our new continuous backbone internal coordinates. CATS generalizes the iMinDEE and EPIC protein design algorithms, which model continuous flexibility in sidechain dihedrals, to model continuous, appropriately localized flexibility in the backbone dihedrals ϕ and ψ as well. We show using 81 test cases based on 29 different protein structures that CATS finds sequences and conformations that are significantly lower in energy than methods with less or no backbone flexibility do. In particular, we show that CATS can model the viability of an antibody mutation known experimentally to increase affinity, but that appears sterically infeasible when modeled with less or no backbone flexibility.

Availability and implementation

Our code is available as free software at https://github.com/donaldlab/OSPREY_refactor.

Supplementary information

Supplementary data are available at Bioinformatics online.