We present an integrated approach that predicts and validates novel anti-cancer drug targets. We first built a classifier that integrates a variety of genomic and systematic datasets to prioritize drug targets specific for breast, pancreatic and ovarian cancer. We then devised strategies to inhibit these anti-cancer drug targets and selected a set of targets that are amenable to inhibition by small molecules, antibodies and synthetic peptides. We validated the predicted drug targets by showing strong anti-proliferative effects of both synthetic peptide and small molecule inhibitors against our predicted targets.

We report the design, construction, and validation of a highly diverse phage-displayed naïve ubiquitin variant (Ubv) library. We first conducted a mutation tolerance scan of 27 residues and confirmed that 24 of these could be substituted by chemically diverse amino acids without compromising the display of Ubvs on phage. Subsequently, we constructed a library containing 6.8×1010 unique members, in which these 24 positions were diversified with a degenerate codon that encodes for 6 aa that are prevalent in protein interaction sites. To ensure the optimal structural stability of the Ubvs, we constructed the library in a two-step process, whereby 12 positions were randomized first, and following the selection for displayed Ubvs, the resulting pool was further diversified at the other 12 positions. The resulting library was validated by conducting binding selections against a panel of 40 diverse protein antigens and was found to be as functional as a highly validated synthetic antibody library, yielding binders against 30 of the antigens. Detailed characterization of an Ubv that bound to the cell-surface receptor human epidermal growth factor receptor 3 revealed tight binding in the single-digit nanomolar range. Moreover, Ubvs that bound to two distinct sites on the intracellular adapter Grb2 could be combined to generate a potent inhibitor that functioned in cells. These results validate ubiquitin as a robust scaffold for the construction of naïve libraries that can be used to generate Ubvs that target signaling networks both outside and inside the cells.

The T cell antigen receptor (TCR)-peptide-major histocompatibility complex (MHC) interface is composed of conserved and diverse regions, yet the relative contribution of each in shaping recognition by T cells remains unclear. Here we isolated cross-reactive peptides with limited homology, which allowed us to compare the structural properties of nine peptides for a single TCR-MHC pair. The TCR's cross-reactivity was rooted in highly similar recognition of an apical 'hot-spot' position in the peptide with tolerance of sequence variation at ancillary positions. Furthermore, we found a striking structural convergence onto a germline-mediated interaction between the TCR CDR1α region and the MHC α2 helix in twelve TCR-peptide-MHC complexes. Our studies suggest that TCR-MHC germline-mediated constraints, together with a focus on a small peptide hot spot, might place limits on peptide antigen cross-reactivity.

Post-transcriptional regulation of mRNAs plays an essential role in the control of gene expression. mRNAs are regulated in ribonucleoprotein (RNP) complexes by RNA-binding proteins (RBPs) along with associated protein and noncoding RNA (ncRNA) cofactors. A global understanding of post-transcriptional control in any cell type requires identification of the components of all of its RNP complexes. We have previously shown that these complexes can be purified by immunoprecipitation using anti-RBP synthetic antibodies produced by phage display. To develop the large number of synthetic antibodies required for a global analysis of RNP complex composition, we have established a pipeline that combines (i) a computationally aided strategy for design of antigens located outside of annotated domains, (ii) high-throughput antigen expression and purification in Escherichia coli, and (iii) high-throughput antibody selection and screening. Using this pipeline, we have produced 279 antibodies against 61 different protein components of Drosophila melanogaster RNPs. Together with those produced in our low-throughput efforts, we have a panel of 311 antibodies for 67 RNP complex proteins. Tests of a subset of our antibodies demonstrated that 89% immunoprecipitate their endogenous target from embryo lysate. This panel of antibodies will serve as a resource for global studies of RNP complexes in Drosophila. Furthermore, our high-throughput pipeline permits efficient production of synthetic antibodies against any large set of proteins.

A synthetic phage-displayed antibody repertoire was constructed with equivalent chemical diversity in the third complementarity-determining regions of the heavy (CDR-H3) and light (CDR-L3) chains, which contrasts with natural antibodies in which CDR-H3 is much more diverse than CDR-L3 due to the genetic mechanisms that generate antibody encoding genes. Surprisingly, the synthetic repertoire yielded numerous functional antibodies that contained mutated CDR-L3 sequences but a fixed CDR-H3 sequence. Alanine-scanning analysis of antibodies that recognized 10 different antigens but contained a common CDR-H3 loop showed that, in most cases, the fixed CDR-H3 sequence was able to contribute favorably to antigen recognition, but in some cases, the loop was functionally inert. Structural analysis of one such antibody in complex with antigen showed that the inert CDR-H3 loop was nonetheless highly buried at the antibody-antigen interface. Taken together, these results show that CDR-H3 diversity is not necessarily required for the generation of antibodies that recognize diverse protein antigens with high affinity and specificity, and if given the chance, CDR-L3 readily assumes the dominant role for antigen recognition. These results contrast with the commonly accepted view of antigen recognition derived from the analysis of natural antibodies, in which CDR-H3 is presumed to be dominant and CDR-L3 is presumed to play an auxiliary role. Furthermore, the results show that natural antibody function is genetically constrained, and it should be possible to develop more functional synthetic antibody libraries by expanding the diversity of CDR-L3 beyond what is observed in nature.

Peptide recognition domains and transcription factors play crucial roles in cellular signaling. They bind linear stretches of amino acids or nucleotides, respectively, with high specificity. Experimental techniques that assess the binding specificity of these domains, such as microarrays or phage display, can retrieve thousands of distinct ligands, providing detailed insight into binding specificity. In particular, the advent of next-generation sequencing has recently increased the throughput of such methods by several orders of magnitude. These advances have helped reveal the presence of distinct binding specificity classes that co-exist within a set of ligands interacting with the same target. Here, we introduce a software system called MUSI that can rapidly analyze very large data sets of binding sequences to determine the relevant binding specificity patterns. Our pipeline provides two major advances. First, it can detect previously unrecognized multiple specificity patterns in any data set. Second, it offers integrated processing of very large data sets from next-generation sequencing machines. The results are visualized as multiple sequence logos describing the different binding preferences of the protein under investigation. We demonstrate the performance of MUSI by analyzing recent phage display data for human SH3 domains as well as microarray data for mouse transcription factors.

We applied a yeast-two-hybrid (Y2H) analysis to screen for ubiquitin variant (UbV) inhibitors of a human deubiquitinase (DUB), ubiquitin-specific protease 2 (USP2). The Y2H screen used USP2 as the bait and a prey library consisting of UbVs randomized at four specific positions, which were known to interact with USP2 from phage display analysis. The screen yielded numerous UbVs that bound to USP2 both as a Y2H interaction in vivo and as purified proteins in vitro. The Y2H-derived UbVs inhibited the catalytic activity of USP2 in vitro with nanomolar-range potencies, and they bound and inhibited USP2 in human cells. Mutational and structural analysis showed that potent and selective inhibition could be achieved by just two substitutions in a UbV, which exhibited improved hydrophobic and hydrophilic contacts compared to the wild-type ubiquitin interaction with USP2. Our results establish Y2H as an effective platform for the development of UbV inhibitors of DUBs in vivo, providing an alternative strategy for the analysis of DUBs that are recalcitrant to phage display and other in vitro methods.

Ubiquitin (Ub) is a small protein that post-translationally modifies a variety of substrates in eukaryotic cells to modulate substrate function. The ability of Ub to interact with numerous protein domains makes Ub an attractive scaffold for engineering ubiquitin variants (UbVs) with high target specificity. Previously, we identified a UbV that formed a non-covalent stable dimer via a β-strand exchange, and in the current work we identified and characterized the minimal substitutions in the primary sequence of Ub required to form a higher ordered complex. Using solution angle scattering and X-ray crystallography, we show that a single substitution of residue Gly10 to either Ala or Val is sufficient to convert Ub from a monomer to a dimer. We also investigate contributions to dimer formation by the residues in the surrounding sequence. These results can be used to develop next-generation phage-display libraries of UbVs to engineer new interfaces for protein recognition.

Beta-arrestins (βarrs) critically mediate desensitization, endocytosis and signalling of G protein-coupled receptors (GPCRs), and they scaffold a large number of interaction partners. However, allosteric modulation of their scaffolding abilities and direct targeting of their interaction interfaces to modulate GPCR functions selectively have not been fully explored yet. Here we identified a series of synthetic antibody fragments (Fabs) against different conformations of βarrs from phage display libraries. Several of these Fabs allosterically and selectively modulated the interaction of βarrs with clathrin and ERK MAP kinase. Interestingly, one of these Fabs selectively disrupted βarr-clathrin interaction, and when expressed as an intrabody, it robustly inhibited agonist-induced endocytosis of a broad set of GPCRs without affecting ERK MAP kinase activation. Our data therefore demonstrate the feasibility of selectively targeting βarr interactions using intrabodies and provide a novel framework for fine-tuning GPCR functions with potential therapeutic implications.

The recent Middle East respiratory syndrome coronavirus (MERS-CoV), Ebola and Zika virus outbreaks exemplify the continued threat of (re-)emerging viruses to human health, and our inability to rapidly develop effective therapeutic countermeasures. Many viruses, including MERS-CoV and the Crimean-Congo hemorrhagic fever virus (CCHFV) encode deubiquitinating (DUB) enzymes that are critical for viral replication and pathogenicity. They bind and remove ubiquitin (Ub) and interferon stimulated gene 15 (ISG15) from cellular proteins to suppress host antiviral innate immune responses. A variety of viral DUBs (vDUBs), including the MERS-CoV papain-like protease, are responsible for cleaving the viral replicase polyproteins during replication, and are thereby critical components of the viral replication cycle. Together, this makes vDUBs highly attractive antiviral drug targets. However, structural similarity between the catalytic cores of vDUBs and human DUBs complicates the development of selective small molecule vDUB inhibitors. We have thus developed an alternative strategy to target the vDUB activity through a rational protein design approach. Here, we report the use of phage-displayed ubiquitin variant (UbV) libraries to rapidly identify potent and highly selective protein-based inhibitors targeting the DUB domains of MERS-CoV and CCHFV. UbVs bound the vDUBs with high affinity and specificity to inhibit deubiquitination, deISGylation and in the case of MERS-CoV also viral replicative polyprotein processing. Co-crystallization studies further revealed critical molecular interactions between UbVs and MERS-CoV or CCHFV vDUBs, accounting for the observed binding specificity and high affinity. Finally, expression of UbVs during MERS-CoV infection reduced infectious progeny titers by more than four orders of magnitude, demonstrating the remarkable potency of UbVs as antiviral agents. Our results thereby establish a strategy to produce protein-based inhibitors that could protect against a diverse range of viruses by providing UbVs via mRNA or protein delivery technologies or through transgenic techniques.

Programmable nucleases, such as Cas9, are used for precise genome editing by homology-dependent repair (HDR). However, HDR efficiency is constrained by competition from other double-strand break (DSB) repair pathways, including non-homologous end-joining (NHEJ). We report the discovery of a genetically encoded inhibitor of 53BP1 that increases the efficiency of HDR-dependent genome editing in human and mouse cells. 53BP1 is a key regulator of DSB repair pathway choice in eukaryotic cells and functions to favor NHEJ over HDR by suppressing end resection, which is the rate-limiting step in the initiation of HDR. We screened an existing combinatorial library of engineered ubiquitin variants for inhibitors of 53BP1. Expression of one variant, named i53 (inhibitor of 53BP1), in human and mouse cells, blocked accumulation of 53BP1 at sites of DNA damage and improved gene targeting and chromosomal gene conversion with either double-stranded DNA or single-stranded oligonucleotide donors by up to 5.6-fold. Inhibition of 53BP1 is a robust method to increase efficiency of HDR-based precise genome editing.

RING and U-box E3 ubiquitin ligases regulate diverse eukaryotic processes and have been implicated in numerous diseases, but targeting these enzymes remains a major challenge. We report the development of three ubiquitin variants (UbVs), each binding selectively to the RING or U-box domain of a distinct E3 ligase: monomeric UBE4B, phosphorylated active CBL, or dimeric XIAP. Structural and biochemical analyses revealed that UbVs specifically inhibited the activity of UBE4B or phosphorylated CBL by blocking the E2∼Ub binding site. Surprisingly, the UbV selective for dimeric XIAP formed a dimer to stimulate E3 activity by stabilizing the closed E2∼Ub conformation. We further verified the inhibitory and stimulatory functions of UbVs in cells. Our work provides a general strategy to inhibit or activate RING/U-box E3 ligases and provides a resource for the research community to modulate these enzymes.

SH3 domains are protein modules that mediate protein-protein interactions in many eukaryotic signal transduction pathways. The majority of SH3 domains studied thus far act by binding to proline-rich sequences in partner proteins, but a growing number of studies have revealed alternative recognition mechanisms. We have comprehensively surveyed the specificity landscape of human SH3 domains in an unbiased manner using peptide-phage display and deep sequencing. Based on ∼70,000 unique binding peptides, we obtained 154 specificity profiles for 115 SH3 domains, which reveal that roughly half of the SH3 domains exhibit non-canonical specificities and collectively recognize a wide variety of peptide motifs, most of which were previously unknown. Crystal structures of SH3 domains with two distinct non-canonical specificities revealed novel peptide-binding modes through an extended surface outside of the canonical proline-binding site. Our results constitute a significant contribution toward a complete understanding of the mechanisms underlying SH3-mediated cellular responses.

Natural evolution encodes rich information about the structure and function of biomolecules in the genetic record. Previously, statistical analysis of co-variation patterns in natural protein families has enabled the accurate computation of 3D structures. Here, we explored generating similar information by experimental evolution, starting from a single gene and performing multiple cycles of in vitro mutagenesis and functional selection in Escherichia coli. We evolved two antibiotic resistance proteins, β-lactamase PSE1 and acetyltransferase AAC6, and obtained hundreds of thousands of diverse functional sequences. Using evolutionary coupling analysis, we inferred residue interaction constraints that were in agreement with contacts in known 3D structures, confirming genetic encoding of structural constraints in the selected sequences. Computational protein folding with interaction constraints then yielded 3D structures with the same fold as natural relatives. This work lays the foundation for a new experimental method (3Dseq) for protein structure determination, combining evolution experiments with inference of residue interactions from sequence information. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

Ebola (EBOV) and Marburg (MARV) are members of the Filoviridae family, which continue to emerge and cause sporadic outbreaks of hemorrhagic fever with high mortality rates. Filoviruses utilize their VP40 matrix protein to drive virion assembly and budding, in part, by recruitment of specific WW-domain-bearing host proteins via its conserved PPxY Late (L) domain motif. Here, we screened an array of 115 mammalian, bacterially expressed and purified WW-domains using a PPxY-containing peptide from MARV VP40 (mVP40) to identify novel host interactors. Using this unbiased approach, we identified Yes Associated Protein (YAP) and Transcriptional co-Activator with PDZ-binding motif (TAZ) as novel mVP40 PPxY interactors. YAP and TAZ function as downstream transcriptional effectors of the Hippo signaling pathway that regulates cell proliferation, migration and apoptosis. We demonstrate that ectopic expression of YAP or TAZ along with mVP40 leads to significant inhibition of budding of mVP40 VLPs in a WW-domain/PPxY dependent manner. Moreover, YAP colocalized with mVP40 in the cytoplasm, and inhibition of mVP40 VLP budding was more pronounced when YAP was localized predominantly in the cytoplasm rather than in the nucleus. A key regulator of YAP nuclear/cytoplasmic localization and function is angiomotin (Amot); a multi-PPxY containing protein that strongly interacts with YAP WW-domains. Interestingly, we found that expression of PPxY-containing Amot rescued mVP40 VLP egress from either YAP- or TAZ-mediated inhibition in a PPxY-dependent manner. Importantly, using a stable Amot-knockdown cell line, we found that expression of Amot was critical for efficient egress of mVP40 VLPs as well as egress and spread of authentic MARV in infected cell cultures. In sum, we identified novel negative (YAP/TAZ) and positive (Amot) regulators of MARV VP40-mediated egress, that likely function in part, via competition between host and viral PPxY motifs binding to modular host WW-domains. These findings not only impact our mechanistic understanding of virus budding and spread, but also may impact the development of new antiviral strategies.

Interleukin-4-induced-1 (IL4i1) is an amino acid oxidase secreted from immune cells. Recent observations have suggested that IL4i1 is pro-tumorigenic via unknown mechanisms. As IL4i1 has homologs in snake venoms (L-amino acid oxidases [LAAO]), we used comparative approaches to gain insight into the mechanistic basis of how conserved amino acid oxidases regulate cell fate and function. Using mammalian expressed recombinant proteins, we found that venom LAAO kills cells via hydrogen peroxide generation. By contrast, mammalian IL4i1 is non-cytotoxic and instead elicits a cell protective gene expression program inhibiting ferroptotic redox death by generating indole-3-pyruvate (I3P) from tryptophan. I3P suppresses ferroptosis by direct free radical scavenging and through the activation of an anti-oxidative gene expression program. Thus, the pro-tumor effects of IL4i1 are likely mediated by local anti-ferroptotic pathways via aromatic amino acid metabolism, arguing that an IL4i1 inhibitor may modulate tumor cell death pathways.

Synthetic antibody (Ab) technologies are efficient and cost-effective platforms for the generation of monoclonal Abs against human antigens. Yet, they typically depend on purified proteins, which exclude integral membrane proteins that require the lipid bilayers to support their native structure and function. Here, we present an Ab discovery strategy, termed CellectSeq, for targeting integral membrane proteins on native cells in complex environment. As proof of concept, we targeted three transmembrane proteins linked to cancer, tetraspanin CD151, carbonic anhydrase 9, and integrin-α11. First, we performed in situ cell-based selections to enrich phage-displayed synthetic Ab pools for antigen-specific binders. Then, we designed next-generation sequencing procedures to explore Ab diversities and abundances. Finally, we developed motif-based scoring and sequencing error-filtering algorithms for the comprehensive interrogation of next-generation sequencing pools to identify Abs with high diversities and specificities, even at extremely low abundances, which are very difficult to identify using manual sampling or sequence abundances.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in a global pandemic causing over 195 million infections and more than 4 million fatalities as of July 2021.To date, it has been demonstrated that a number of mutations in the spike glycoprotein (S protein) of SARS-CoV-2 variants of concern abrogate or reduce the neutralization potency of several therapeutic antibodies and vaccine-elicited antibodies. Therefore, the development of additional vaccine platforms with improved supply and logistic profile remains a pressing need. In this work, we have validated the applicability of a peptide-based strategy focused on a preventive as well as a therapeutic purpose. On the basis of the involvement of the dipeptidyl peptidase 4 (DPP4), in addition to the angiotensin converting enzyme 2 (ACE2) receptor in the mechanism of virus entry, we analyzed peptides bearing DPP4 sequences by protein-protein docking and assessed their ability to block pseudovirus infection in vitro. In parallel, we have selected and synthetized peptide sequences located within the highly conserved receptor-binding domain (RBD) of the S protein, and we found that RBD-based vaccines could better promote elicitation of high titers of neutralizing antibodies specific against the regions of interest, as confirmed by immunoinformatic methodologies and in vivo studies. These findings unveil a key antigenic site targeted by broadly neutralizing antibodies and pave the way to the design of pan-coronavirus vaccines.

Rift Valley fever virus (RVFV) is a zoonotic pathogen with pandemic potential. RVFV entry is mediated by the viral glycoprotein (Gn), but host entry factors remain poorly defined. Our genome-wide CRISPR screen identified low-density lipoprotein receptor-related protein 1 (mouse Lrp1/human LRP1), heat shock protein (Grp94), and receptor-associated protein (RAP) as critical host factors for RVFV infection. RVFV Gn directly binds to specific Lrp1 clusters and is glycosylation independent. Exogenous addition of murine RAP domain 3 (mRAPD3) and anti-Lrp1 antibodies neutralizes RVFV infection in taxonomically diverse cell lines. Mice treated with mRAPD3 and infected with pathogenic RVFV are protected from disease and death. A mutant mRAPD3 that binds Lrp1 weakly failed to protect from RVFV infection. Together, these data support Lrp1 as a host entry factor for RVFV infection and define a new target to limit RVFV infections.

Suppressor of cytokine signaling (SOCS)2 protein is a key negative regulator of the growth hormone (GH) and Janus kinase (JAK)-Signal Transducers and Activators of Transcription (STAT) signaling cascade. The central SOCS2-Src homology 2 (SH2) domain is characteristic of the SOCS family proteins and is an important module that facilitates recognition of targets bearing phosphorylated tyrosine (pTyr) residues. Here we identify an exosite on the SOCS2-SH2 domain which, when bound to a non-phosphorylated peptide (F3), enhances SH2 affinity for canonical phosphorylated ligands. Solution of the SOCS2/F3 crystal structure reveals F3 as an α-helix which binds on the opposite side of the SH2 domain to the phosphopeptide binding site. F3:exosite binding appears to stabilise the SOCS2-SH2 domain, resulting in slower dissociation of phosphorylated ligands and consequently, enhances binding affinity. This biophysical enhancement of SH2:pTyr binding affinity translates to increase SOCS2 inhibition of GH signaling.