Elaboration of the Homer1 Recognition Landscape Reveals Incomplete Divergence of Paralogous EVH1 Domains.

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of Short Linear Motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. In doing so, we expanded current understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

Marginal specificity in protein interactions constrains evolution of a paralogous family.

The evolution of novel functions in biology relies heavily on gene duplication and divergence, creating large paralogous protein families. Selective pressure to avoid detrimental cross-talk often results in paralogs that exhibit exquisite specificity for their interaction partners. But how robust or sensitive is this specificity to mutation? Here, using deep mutational scanning, we demonstrate that a paralogous family of bacterial signaling proteins exhibits marginal specificity, such that many individual substitutions give rise to substantial cross-talk between normally insulated pathways. Our results indicate that sequence space is locally crowded despite overall sparseness, and we provide evidence that this crowding has constrained the evolution of bacterial signaling proteins. These findings underscore how evolution selects for "good enough" rather than optimized phenotypes, leading to restrictions on the subsequent evolution of paralogs.

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.

Neural network-derived Potts models for structure-based protein design using backbone atomic coordinates and tertiary motifs.

Designing novel proteins to perform desired functions, such as binding or catalysis, is a major goal in synthetic biology. A variety of computational approaches can aid in this task. An energy-based framework rooted in the sequence-structure statistics of tertiary motifs (TERMs) can be used for sequence design on predefined backbones. Neural network models that use backbone coordinate-derived features provide another way to design new proteins. In this work, we combine the two methods to make neural structure-based models more suitable for protein design. Specifically, we supplement backbone-coordinate features with TERM-derived data, as inputs, and we generate energy functions as outputs. We present two architectures that generate Potts models over the sequence space: TERMinator, which uses both TERM-based and coordinate-based information, and COORDinator, which uses only coordinate-based information. Using these two models, we demonstrate that TERMs can be utilized to improve native sequence recovery performance of neural models. Furthermore, we demonstrate that sequences designed by TERMinator are predicted to fold to their target structures by AlphaFold. Finally, we show that both TERMinator and COORDinator learn notions of energetics, and these methods can be fine-tuned on experimental data to improve predictions. Our results suggest that using TERM-based and coordinate-based features together may be beneficial for protein design and that structure-based neural models that produce Potts energy tables have utility for flexible applications in protein science.

High-throughput computational discovery of inhibitory protein fragments with AlphaFold.

Peptides can bind to specific sites on larger proteins and thereby function as inhibitors and regulatory elements. Peptide fragments of larger proteins are particularly attractive for achieving these functions due to their inherent potential to form native-like binding interactions. Recently developed experimental approaches allow for high-throughput measurement of protein fragment inhibitory activity in living cells. However, it has thus far not been possible to predict de novo which of the many possible protein fragments bind their protein targets, let alone act as inhibitors. We have developed a computational method, FragFold, that employs AlphaFold to predict protein fragment binding to full-length protein targets in a high-throughput manner. Applying FragFold to thousands of fragments tiling across diverse proteins revealed peaks of predicted binding along each protein sequence. These predictions were compared with experimentally measured peaks of inhibitory activity in E. coli. We establish that our approach is a sensitive predictor of protein fragment function: Evaluating inhibitory fragments derived from known protein-protein interaction interfaces, we found 87% were predicted by FragFold to bind in a native-like mode. Across full protein sequences, 68% of FragFold-predicted binding peaks match experimentally measured inhibitory peaks. This is true even when the underlying inhibitory mechanism is unclear from existing structural data, and we find FragFold is able to predict novel binding modes for inhibitory fragments of unknown structure, explaining previous genetic and biochemical data for these fragments. The success rate of FragFold demonstrates that this computational approach should be broadly applicable for discovering inhibitory protein fragments across proteomes.

Molecular determinants of TRAF6 binding specificity suggest that native interaction partners are not optimized for affinity.

TRAF6 is an adaptor protein involved in signaling pathways that are essential for development and the immune system. It participates in many protein-protein interactions, some of which are mediated by the C-terminal MATH domain, which binds to short peptide segments containing the motif PxExx[FYWHDE], where x is any amino acid. Blocking MATH domain interactions is associated with favorable effects in various disease models. To better define TRAF6 MATH domain binding preferences, we screened a combinatorial library using bacterial cell-surface peptide display. We identified 236 of the best TRAF6-interacting peptides and a set of 1,200 peptides that match the sequence PxE but do not bind TRAF6 MATH. The peptides that were most enriched in the screen bound TRAF6 tighter than previously measured native peptides. To better understand the structural basis for TRAF6 interaction preferences, we built all-atom structural models of the MATH domain in complex with high-affinity binders and nonbinders identified in the screen. We identified favorable interactions for motif features in binders as well as negative design elements distributed across the motif that can disfavor or preclude binding. Searching the human proteome revealed that the most biologically relevant TRAF6 motif matches occupy a different sequence space from the best hits discovered in combinatorial library screening, suggesting that native interactions are not optimized for affinity. Our experimentally determined binding preferences and structural models support the design of peptide-based interaction inhibitors with higher affinities than endogenous TRAF6 ligands.

Tertiary motifs as building blocks for the design of protein-binding peptides.

Despite advances in protein engineering, the de novo design of small proteins or peptides that bind to a desired target remains a difficult task. Most computational methods search for binder structures in a library of candidate scaffolds, which can lead to designs with poor target complementarity and low success rates. Instead of choosing from pre-defined scaffolds, we propose that custom peptide structures can be constructed to complement a target surface. Our method mines tertiary motifs (TERMs) from known structures to identify surface-complementing fragments or "seeds." We combine seeds that satisfy geometric overlap criteria to generate peptide backbones and score the backbones to identify the most likely binding structures. We found that TERM-based seeds can describe known binding structures with high resolution: the vast majority of peptide binders from 486 peptide-protein complexes can be covered by seeds generated from single-chain structures. Furthermore, we demonstrate that known peptide structures can be reconstructed with high accuracy from peptide-covering seeds. As a proof of concept, we used our method to design 100 peptide binders of TRAF6, seven of which were predicted by Rosetta to form higher-quality interfaces than a native binder. The designed peptides interact with distinct sites on TRAF6, including the native peptide-binding site. These results demonstrate that known peptide-binding structures can be constructed from TERMs in single-chain structures and suggest that TERM information can be applied to efficiently design novel target-complementing binders.

Native proline-rich motifs exploit sequence context to target actin-remodeling Ena/VASP protein ENAH.

The human proteome is replete with short linear motifs (SLiMs) of four to six residues that are critical for protein-protein interactions, yet the importance of the sequence surrounding such motifs is underexplored. We devised a proteomic screen to examine the influence of SLiM sequence context on protein-protein interactions. Focusing on the EVH1 domain of human ENAH, an actin regulator that is highly expressed in invasive cancers, we screened 36-residue proteome-derived peptides and discovered new interaction partners of ENAH and diverse mechanisms by which context influences binding. A pocket on the ENAH EVH1 domain that has diverged from other Ena/VASP paralogs recognizes extended SLiMs and favors motif-flanking proline residues. Many high-affinity ENAH binders that contain two proline-rich SLiMs use a noncanonical site on the EVH1 domain for binding and display a thermodynamic signature consistent with the two-motif chain engaging a single domain. We also found that photoreceptor cilium actin regulator (PCARE) uses an extended 23-residue region to obtain a higher affinity than any known ENAH EVH1-binding motif. Our screen provides a way to uncover the effects of proteomic context on motif-mediated binding, revealing diverse mechanisms of control over EVH1 interactions and establishing that SLiMs can't be fully understood outside of their native context.

A distributed residue network permits conformational binding specificity in a conserved family of actin remodelers.

Metazoan proteomes contain many paralogous proteins that have evolved distinct functions. The Ena/VASP family of actin regulators consists of three members that share an EVH1 interaction domain with a 100 % conserved binding site. A proteome-wide screen revealed photoreceptor cilium actin regulator (PCARE) as a high-affinity ligand for ENAH EVH1. Here, we report the surprising observation that PCARE is ~100-fold specific for ENAH over paralogs VASP and EVL and can selectively bind ENAH and inhibit ENAH-dependent adhesion in cells. Specificity arises from a mechanism whereby PCARE stabilizes a conformation of the ENAH EVH1 domain that is inaccessible to family members VASP and EVL. Structure-based modeling rapidly identified seven residues distributed throughout EVL that are sufficient to differentiate binding by ENAH vs. EVL. By exploiting the ENAH-specific conformation, we rationally designed the tightest and most selective ENAH binder to date. Our work uncovers a conformational mechanism of interaction specificity that distinguishes highly similar paralogs and establishes tools for dissecting specific Ena/VASP functions in processes including cancer cell invasion.

Data-driven computational protein design.

Computational protein design can generate proteins not found in nature that adopt desired structures and perform novel functions. Although proteins could, in theory, be designed with ab initio methods, practical success has come from using large amounts of data that describe the sequences, structures, and functions of existing proteins and their variants. We present recent creative uses of multiple-sequence alignments, protein structures, and high-throughput functional assays in computational protein design. Approaches range from enhancing structure-based design with experimental data to building regression models to training deep neural nets that generate novel sequences. Looking ahead, deep learning will be increasingly important for maximizing the value of data for protein design.

Spatial Multiplexing of Fluorescent Reporters for Imaging Signaling Network Dynamics.

In order to analyze how a signal transduction network converts cellular inputs into cellular outputs, ideally one would measure the dynamics of many signals within the network simultaneously. We found that, by fusing a fluorescent reporter to a pair of self-assembling peptides, it could be stably clustered within cells at random points, distant enough to be resolved by a microscope but close enough to spatially sample the relevant biology. Because such clusters, which we call signaling reporter islands (SiRIs), can be modularly designed, they permit a set of fluorescent reporters to be efficiently adapted for simultaneous measurement of multiple nodes of a signal transduction network within single cells. We created SiRIs for indicators of second messengers and kinases and used them, in hippocampal neurons in culture and intact brain slices, to discover relationships between the speed of calcium signaling, and the amplitude of PKA signaling, upon receiving a cAMP-driving stimulus.

Precision Calcium Imaging of Dense Neural Populations via a Cell-Body-Targeted Calcium Indicator.

Methods for one-photon fluorescent imaging of calcium dynamics can capture the activity of hundreds of neurons across large fields of view at a low equipment complexity and cost. In contrast to two-photon methods, however, one-photon methods suffer from higher levels of crosstalk from neuropil, resulting in a decreased signal-to-noise ratio and artifactual correlations of neural activity. We address this problem by engineering cell-body-targeted variants of the fluorescent calcium indicators GCaMP6f and GCaMP7f. We screened fusions of GCaMP to natural, as well as artificial, peptides and identified fusions that localized GCaMP to within 50 µm of the cell body of neurons in mice and larval zebrafish. One-photon imaging of soma-targeted GCaMP in dense neural circuits reported fewer artifactual spikes from neuropil, an increased signal-to-noise ratio, and decreased artifactual correlation across neurons. Thus, soma-targeting of fluorescent calcium indicators facilitates usage of simple, powerful, one-photon methods for imaging neural calcium dynamics.

Tertiary Structural Motif Sequence Statistics Enable Facile Prediction and Design of Peptides that Bind Anti-apoptotic Bfl-1 and Mcl-1.

Understanding the relationship between protein sequence and structure well enough to design new proteins with desired functions is a longstanding goal in protein science. Here, we show that recurring tertiary structural motifs (TERMs) in the PDB provide rich information for protein-peptide interaction prediction and design. TERM statistics can be used to predict peptide binding energies for Bcl-2 family proteins as accurately as widely used structure-based tools. Furthermore, design using TERM energies (dTERMen) rapidly and reliably generates high-affinity peptide binders of anti-apoptotic proteins Bfl-1 and Mcl-1 with just 15%-38% sequence identity to any known native Bcl-2 family protein ligand. High-resolution structures of four designed peptides bound to their targets provide opportunities to analyze the strengths and limitations of the computational design method. Our results support dTERMen as a powerful approach that can complement existing tools for protein engineering.

Peptide design by optimization on a data-parameterized protein interaction landscape.

Many applications in protein engineering require optimizing multiple protein properties simultaneously, such as binding one target but not others or binding a target while maintaining stability. Such multistate design problems require navigating a high-dimensional space to find proteins with desired characteristics. A model that relates protein sequence to functional attributes can guide design to solutions that would be hard to discover via screening. In this work, we measured thousands of protein-peptide binding affinities with the high-throughput interaction assay amped SORTCERY and used the data to parameterize a model of the alpha-helical peptide-binding landscape for three members of the Bcl-2 family of proteins: Bcl-xL, Mcl-1, and Bfl-1. We applied optimization protocols to explore extremes in this landscape to discover peptides with desired interaction profiles. Computational design generated 36 peptides, all of which bound with high affinity and specificity to just one of Bcl-xL, Mcl-1, or Bfl-1, as intended. We designed additional peptides that bound selectively to two out of three of these proteins. The designed peptides were dissimilar to known Bcl-2-binding peptides, and high-resolution crystal structures confirmed that they engaged their targets as expected. Excellent results on this challenging problem demonstrate the power of a landscape modeling approach, and the designed peptides have potential uses as diagnostic tools or cancer therapeutics.

Iterative optimization yields Mcl-1-targeting stapled peptides with selective cytotoxicity to Mcl-1-dependent cancer cells.

Bcl-2 family proteins regulate apoptosis, and aberrant interactions of overexpressed antiapoptotic family members such as Mcl-1 promote cell transformation, cancer survival, and resistance to chemotherapy. Discovering potent and selective Mcl-1 inhibitors that can relieve apoptotic blockades is thus a high priority for cancer research. An attractive strategy for disabling Mcl-1 involves using designer peptides to competitively engage its binding groove, mimicking the structural mechanism of action of native sensitizer BH3-only proteins. We transformed Mcl-1-binding peptides into α-helical, cell-penetrating constructs that are selectively cytotoxic to Mcl-1-dependent cancer cells. Critical to the design of effective inhibitors was our introduction of an all-hydrocarbon cross-link or "staple" that stabilizes α-helical structure, increases target binding affinity, and independently confers binding specificity for Mcl-1 over related Bcl-2 family paralogs. Two crystal structures of complexes at 1.4 Å and 1.9 Å resolution demonstrate how the hydrophobic staple induces an unanticipated structural rearrangement in Mcl-1 upon binding. Systematic sampling of staple location and iterative optimization of peptide sequence in accordance with established design principles provided peptides that target intracellular Mcl-1. This work provides proof of concept for the development of potent, selective, and cell-permeable stapled peptides for therapeutic targeting of Mcl-1 in cancer, applying a design and validation workflow applicable to a host of challenging biomedical targets.

PixelDB: Protein-peptide complexes annotated with structural conservation of the peptide binding mode.

PixelDB, the Peptide Exosite Location Database, compiles 1966 non-redundant, high-resolution structures of protein-peptide complexes filtered to minimize the impact of crystal packing on peptide conformation. The database is organized to facilitate study of structurally conserved versus non-conserved elements of protein-peptide engagement. PixelDB clusters complexes based on the structural similarity of the peptide-binding protein, and by comparing complexes within a cluster highlights examples of domains that engage peptides using more than one binding mode. PixelDB also identifies conserved peptide core structural motifs characteristic of each binding mode. Peptide regions that flank core motifs often make non-structurally conserved interactions with the protein surface in regions we call exosites. Many examples establish that exosite contacts can be important for enhancing protein binding and interaction specificity. PixelDB provides a resource for computational and structural biologists to study, model, and predict core-motif and exosite-contacting peptide interactions. PixelDB is available to the community without restriction in a convenient flat-file format with accompanying visualization tools.

Modular assembly of a protein nanotriangle using orthogonally interacting coiled coils.

Synthetic protein assemblies that adopt programmed shapes would support many applications in nanotechnology. We used a rational design approach that exploits the modularity of orthogonally interacting coiled coils to create a self-assembled protein nanotriangle. Coiled coils have frequently been used to construct nanoassemblies and materials, but rarely with successful prior specification of the resulting structure. We designed a heterotrimer from three pairs of heterodimeric coiled coils that mediate specific interactions while avoiding undesired crosstalk. Non-associating pairs of coiled-coil units were strategically fused to generate three chains that were predicted to preferentially form the heterotrimer, and a rational annealing process led to the desired oligomer. Extensive biophysical characterization and modeling support the formation of a molecular triangle, which is a shape distinct from naturally occurring supramolecular nanostructures. Our approach can be extended to design more complex nanostructures using additional coiled-coil modules, other protein parts, or templated surfaces.

Epistatic mutations in PUMA BH3 drive an alternate binding mode to potently and selectively inhibit anti-apoptotic Bfl-1.

Overexpression of anti-apoptotic Bcl-2 family proteins contributes to cancer progression and confers resistance to chemotherapy. Small molecules that target Bcl-2 are used in the clinic to treat leukemia, but tight and selective inhibitors are not available for Bcl-2 paralog Bfl-1. Guided by computational analysis, we designed variants of the native BH3 motif PUMA that are > 150-fold selective for Bfl-1 binding. The designed peptides potently trigger disruption of the mitochondrial outer membrane in cells dependent on Bfl-1, but not in cells dependent on other anti-apoptotic homologs. High-resolution crystal structures show that designed peptide FS2 binds Bfl-1 in a shifted geometry, relative to PUMA and other binding partners, due to a set of epistatic mutations. FS2 modified with an electrophile reacts with a cysteine near the peptide-binding groove to augment specificity. Designed Bfl-1 binders provide reagents for cellular profiling and leads for developing enhanced and cell-permeable peptide or small-molecule inhibitors.

Combinatorial bZIP dimers display complex DNA-binding specificity landscapes.

How transcription factor dimerization impacts DNA-binding specificity is poorly understood. Guided by protein dimerization properties, we examined DNA binding specificities of 270 human bZIP pairs. DNA interactomes of 80 heterodimers and 22 homodimers revealed that 72% of heterodimer motifs correspond to conjoined half-sites preferred by partnering monomers. Remarkably, the remaining motifs are composed of variably-spaced half-sites (12%) or 'emergent' sites (16%) that cannot be readily inferred from half-site preferences of partnering monomers. These binding sites were biochemically validated by EMSA-FRET analysis and validated in vivo by ChIP-seq data from human cell lines. Focusing on ATF3, we observed distinct cognate site preferences conferred by different bZIP partners, and demonstrated that genome-wide binding of ATF3 is best explained by considering many dimers in which it participates. Importantly, our compendium of bZIP-DNA interactomes predicted bZIP binding to 156 disease associated SNPs, of which only 20 were previously annotated with known bZIP motifs.