TCRmodel2: high-resolution modeling of T cell receptor recognition using deep learning.

The cellular immune system, which is a critical component of human immunity, uses T cell receptors (TCRs) to recognize antigenic proteins in the form of peptides presented by major histocompatibility complex (MHC) proteins. Accurate definition of the structural basis of TCRs and their engagement of peptide-MHCs can provide major insights into normal and aberrant immunity, and can help guide the design of vaccines and immunotherapeutics. Given the limited amount of experimentally determined TCR-peptide-MHC structures and the vast amount of TCRs within each individual as well as antigenic targets, accurate computational modeling approaches are needed. Here, we report a major update to our web server, TCRmodel, which was originally developed to model unbound TCRs from sequence, to now model TCR-peptide-MHC complexes from sequence, utilizing several adaptations of AlphaFold. This method, named TCRmodel2, allows users to submit sequences through an easy-to-use interface and shows similar or greater accuracy than AlphaFold and other methods to model TCR-peptide-MHC complexes based on benchmarking. It can generate models of complexes in 15 minutes, and output models are provided with confidence scores and an integrated molecular viewer. TCRmodel2 is available at https://tcrmodel.ibbr.umd.edu.

Rationally designed inhibitors of the Musashi protein-RNA interaction by hotspot mimicry.

RNA-binding proteins (RBPs) are key post-transcriptional regulators of gene expression, and thus underlie many important biological processes. Here, we developed a strategy that entails extracting a "hotspot pharmacophore" from the structure of a protein-RNA complex, to create a template for designing small-molecule inhibitors and for exploring the selectivity of the resulting inhibitors. We demonstrate this approach by designing inhibitors of Musashi proteins MSI1 and MSI2, key regulators of mRNA stability and translation that are upregulated in many cancers. We report this novel series of MSI1/MSI2 inhibitors is specific and active in biochemical, biophysical, and cellular assays. This study extends the paradigm of "hotspots" from protein-protein complexes to protein-RNA complexes, supports the "druggability" of RNA-binding protein surfaces, and represents one of the first rationally-designed inhibitors of non-enzymatic RNA-binding proteins. Owing to its simplicity and generality, we anticipate that this approach may also be used to develop inhibitors of many other RNA-binding proteins; we also consider the prospects of identifying potential off-target interactions by searching for other RBPs that recognize their cognate RNAs using similar interaction geometries. Beyond inhibitors, we also expect that compounds designed using this approach can serve as warheads for new PROTACs that selectively degrade RNA-binding proteins.

Rationally designed inhibitors of the Musashi protein-RNA interaction by hotspot mimicry.

RNA-binding proteins (RBPs) are key post-transcriptional regulators of gene expression, and thus underlie many important biological processes. Here, we developed a strategy that entails extracting a "hotspot pharmacophore" from the structure of a protein-RNA complex, to create a template for designing small-molecule inhibitors and for exploring the selectivity of the resulting inhibitors. We demonstrate this approach by designing inhibitors of Musashi proteins MSI1 and MSI2, key regulators of mRNA stability and translation that are upregulated in many cancers. We report this novel series of MSI1/MSI2 inhibitors is specific and active in biochemical, biophysical, and cellular assays. This study extends the paradigm of "hotspots" from protein-protein complexes to protein-RNA complexes, supports the "druggability" of RNA-binding protein surfaces, and represents one of the first rationally-designed inhibitors of non-enzymatic RNA-binding proteins. Owing to its simplicity and generality, we anticipate that this approach may also be used to develop inhibitors of many other RNA-binding proteins; we also consider the prospects of identifying potential off-target interactions by searching for other RBPs that recognize their cognate RNAs using similar interaction geometries. Beyond inhibitors, we also expect that compounds designed using this approach can serve as warheads for new PROTACs that selectively degrade RNA-binding proteins.

Structural assessment of HLA-A2-restricted SARS-CoV-2 spike epitopes recognized by public and private T-cell receptors.

T cells play a vital role in combatting SARS-CoV-2 and forming long-term memory responses. Whereas extensive structural information is available on neutralizing antibodies against SARS-CoV-2, such information on SARS-CoV-2-specific T-cell receptors (TCRs) bound to their peptide-MHC targets is lacking. Here we determine the structures of a public and a private TCR from COVID-19 convalescent patients in complex with HLA-A2 and two SARS-CoV-2 spike protein epitopes (YLQ and RLQ). The structures reveal the basis for selection of particular TRAV and TRBV germline genes by the public but not the private TCR, and for the ability of the TCRs to recognize natural variants of RLQ but not YLQ. Neither TCR recognizes homologous epitopes from human seasonal coronaviruses. By elucidating the mechanism for TCR recognition of an immunodominant yet variable epitope (YLQ) and a conserved but less commonly targeted epitope (RLQ), this study can inform prospective efforts to design vaccines to elicit pan-coronavirus immunity.

Identification and Validation of an Aspergillus nidulans Secondary Metabolite Derivative as an Inhibitor of the Musashi-RNA Interaction.

RNA-binding protein Musashi-1 (MSI1) is a key regulator of several stem cell populations. MSI1 is involved in tumor proliferation and maintenance, and it regulates target mRNAs at the translational level. The known mRNA targets of MSI1 include Numb, APC, and P21WAF-1, key regulators of Notch/Wnt signaling and cell cycle progression, respectively. In this study, we aim to identify small molecule inhibitors of MSI1-mRNA interactions, which could block the growth of cancer cells with high levels of MSI1. Using a fluorescence polarization (FP) assay, we screened small molecules from several chemical libraries for those that disrupt the binding of MSI1 to its consensus RNA. One cluster of hit compounds is the derivatives of secondary metabolites from Aspergillus nidulans. One of the top hits, Aza-9, from this cluster was further validated by surface plasmon resonance and nuclear magnetic resonance spectroscopy, which demonstrated that Aza-9 binds directly to MSI1, and the binding is at the RNA binding pocket. We also show that Aza-9 binds to Musashi-2 (MSI2) as well. To test whether Aza-9 has anti-cancer potential, we used liposomes to facilitate Aza-9 cellular uptake. Aza-9-liposome inhibits proliferation, induces apoptosis and autophagy, and down-regulates Notch and Wnt signaling in colon cancer cell lines. In conclusion, we identified a series of potential lead compounds for inhibiting MSI1/2 function, while establishing a framework for identifying small molecule inhibitors of RNA binding proteins using FP-based screening methodology.

Structural basis for oligoclonal T cell recognition of a shared p53 cancer neoantigen.

Adoptive cell therapy (ACT) with tumor-specific T cells can mediate cancer regression. The main target of tumor-specific T cells are neoantigens arising from mutations in self-proteins. Although the majority of cancer neoantigens are unique to each patient, and therefore not broadly useful for ACT, some are shared. We studied oligoclonal T-cell receptors (TCRs) that recognize a shared neoepitope arising from a driver mutation in the p53 oncogene (p53R175H) presented by HLA-A2. Here we report structures of wild-type and mutant p53-HLA-A2 ligands, as well as structures of three tumor-specific TCRs bound to p53R175H-HLA-A2. These structures reveal how a driver mutation in p53 rendered a self-peptide visible to T cells. The TCRs employ structurally distinct strategies that are highly focused on the mutation to discriminate between mutant and wild-type p53. The TCR-p53R175H-HLA-A2 complexes provide a framework for designing TCRs to improve potency for ACT without sacrificing specificity.

Targeting the interaction between RNA-binding protein HuR and FOXQ1 suppresses breast cancer invasion and metastasis.

Patients diagnosed with metastatic breast cancer have a dismal 5-year survival rate of only 24%. The RNA-binding protein Hu antigen R (HuR) is upregulated in breast cancer, and elevated cytoplasmic HuR correlates with high-grade tumors and poor clinical outcome of breast cancer. HuR promotes tumorigenesis by regulating numerous proto-oncogenes, growth factors, and cytokines that support major tumor hallmarks including invasion and metastasis. Here, we report a HuR inhibitor KH-3, which potently suppresses breast cancer cell growth and invasion. Furthermore, KH-3 inhibits breast cancer experimental lung metastasis, improves mouse survival, and reduces orthotopic tumor growth. Mechanistically, we identify FOXQ1 as a direct target of HuR. KH-3 disrupts HuR-FOXQ1 mRNA interaction, leading to inhibition of breast cancer invasion. Our study suggests that inhibiting HuR is a promising therapeutic strategy for lethal metastatic breast cancer.

Cellular origins and genetic landscape of cutaneous gamma delta T cell lymphomas.

Primary cutaneous γδ T cell lymphomas (PCGDTLs) represent a heterogeneous group of uncommon but aggressive cancers. Herein, we perform genome-wide DNA, RNA, and T cell receptor (TCR) sequencing on 29 cutaneous γδ lymphomas. We find that PCGDTLs are not uniformly derived from Vδ2 cells. Instead, the cell-of-origin depends on the tissue compartment from which the lymphomas are derived. Lymphomas arising from the outer layer of skin are derived from Vδ1 cells, the predominant γδ cell in the epidermis and dermis. In contrast, panniculitic lymphomas arise from Vδ2 cells, the predominant γδ T cell in the fat. We also show that TCR chain usage is non-random, suggesting common antigens for Vδ1 and Vδ2 lymphomas respectively. In addition, Vδ1 and Vδ2 PCGDTLs harbor similar genomic landscapes with potentially targetable oncogenic mutations in the JAK/STAT, MAPK, MYC, and chromatin modification pathways. Collectively, these findings suggest a paradigm for classifying, staging, and treating these diseases.

TCR3d: The T cell receptor structural repertoire database.

SUMMARY: T cell receptors (TCRs) are critical molecules of the adaptive immune system, capable of recognizing diverse antigens, including peptides, lipids and small molecules, and represent a rapidly growing class of therapeutics. Determining the structural and mechanistic basis of TCR targeting of antigens is a major challenge, as each individual has a vast and diverse repertoire of TCRs. Despite shared general recognition modes, diversity in TCR sequence and recognition represents a challenge to predictive modeling and computational techniques being developed to predict antigen specificity and mechanistic basis of TCR targeting. To this end, we have developed the TCR3d database, a resource containing all known TCR structures, with a particular focus on antigen recognition. TCR3d provides key information on antigen binding mode, interface features, loop sequences and germline gene usage. Users can interactively view TCR complex structures, search sequences of interest against known structures and sequences, and download curated datasets of structurally characterized TCR complexes. This database is updated on a weekly basis, and can serve the community as a centralized resource for those studying T cell receptors and their recognition. AVAILABILITY AND IMPLEMENTATION: The TCR3d database is available at https://tcr3d.ibbr.umd.edu/.

Natural product derivative Gossypolone inhibits Musashi family of RNA-binding proteins.

BACKGROUND: The Musashi (MSI) family of RNA-binding proteins is best known for the role in post-transcriptional regulation of target mRNAs. Elevated MSI1 levels in a variety of human cancer are associated with up-regulation of Notch/Wnt signaling. MSI1 binds to and negatively regulates translation of Numb and APC (adenomatous polyposis coli), negative regulators of Notch and Wnt signaling respectively. METHODS: Previously, we have shown that the natural product (-)-gossypol as the first known small molecule inhibitor of MSI1 that down-regulates Notch/Wnt signaling and inhibits tumor xenograft growth in vivo. Using a fluorescence polarization (FP) competition assay, we identified gossypolone (Gn) with a > 20-fold increase in Ki value compared to (-)-gossypol. We validated Gn binding to MSI1 using surface plasmon resonance, nuclear magnetic resonance, and cellular thermal shift assay, and tested the effects of Gn on colon cancer cells and colon cancer DLD-1 xenografts in nude mice. RESULTS: In colon cancer cells, Gn reduced Notch/Wnt signaling and induced apoptosis. Compared to (-)-gossypol, the same concentration of Gn is less active in all the cell assays tested. To increase Gn bioavailability, we used PEGylated liposomes in our in vivo studies. Gn-lip via tail vein injection inhibited the growth of human colon cancer DLD-1 xenografts in nude mice, as compared to the untreated control (P < 0.01, n = 10). CONCLUSION: Our data suggest that PEGylation improved the bioavailability of Gn as well as achieved tumor-targeted delivery and controlled release of Gn, which enhanced its overall biocompatibility and drug efficacy in vivo. This provides proof of concept for the development of Gn-lip as a molecular therapy for colon cancer with MSI1/MSI2 overexpression.

Peptide-MHC (pMHC) binding to a human antiviral T cell receptor induces long-range allosteric communication between pMHC- and CD3-binding sites.

T cells generate adaptive immune responses mediated by the T cell receptor (TCR)-CD3 complex comprising an αß TCR heterodimer noncovalently associated with three CD3 dimers. In early T cell activation, αß TCR engagement by peptide-major histocompatibility complex (pMHC) is first communicated to the CD3 signaling apparatus of the TCR-CD3 complex, but the underlying mechanism is incompletely understood. It is possible that pMHC binding induces allosteric changes in TCR conformation or dynamics that are then relayed to CD3. Here, we carried out NMR analysis and molecular dynamics (MD) simulations of both the α and ß chains of a human antiviral TCR (A6) that recognizes the Tax antigen from human T cell lymphotropic virus-1 bound to the MHC class I molecule HLA-A2. We observed pMHC-induced NMR signal perturbations in the TCR variable (V) domains that propagated to three distinct sites in the constant (C) domains: 1) the Cß FG loop projecting from the Vß/Cß interface; 2) a cluster of Cß residues near the Cß αA helix, a region involved in interactions with CD3; and 3) the Cα AB loop at the membrane-proximal base of the TCR. A biological role for each of these allosteric sites is supported by previous mutational and functional studies of TCR signaling. Moreover, the pattern of long-range, ligand-induced changes in TCR A6 revealed by NMR was broadly similar to that predicted by the MD simulations. We propose that the unique structure of the TCR ß chain enables allosteric communication between the TCR-binding sites for pMHC and CD3.

TCRmodel: high resolution modeling of T cell receptors from sequence.

T cell receptors (TCRs), along with antibodies, are responsible for specific antigen recognition in the adaptive immune response, and millions of unique TCRs are estimated to be present in each individual. Understanding the structural basis of TCR targeting has implications in vaccine design, autoimmunity, as well as T cell therapies for cancer. Given advances in deep sequencing leading to immune repertoire-level TCR sequence data, fast and accurate modeling methods are needed to elucidate shared and unique 3D structural features of these molecules which lead to their antigen targeting and cross-reactivity. We developed a new algorithm in the program Rosetta to model TCRs from sequence, and implemented this functionality in a web server, TCRmodel. This web server provides an easy to use interface, and models are generated quickly that users can investigate in the browser and download. Benchmarking of this method using a set of nonredundant recently released TCR crystal structures shows that models are accurate and compare favorably to models from another available modeling method. This server enables the community to obtain insights into TCRs of interest, and can be combined with methods to model and design TCR recognition of antigens. The TCRmodel server is available at: http://tcrmodel.ibbr.umd.edu/.

DARC: Mapping Surface Topography by Ray-Casting for Effective Virtual Screening at Protein Interaction Sites.

Protein-protein interactions represent an exciting and challenging target class for therapeutic intervention using small molecules. Protein interaction sites are often devoid of the deep surface pockets presented by "traditional" drug targets, and crystal structures reveal that inhibitors typically engage these sites using very shallow binding modes. As a consequence, modern virtual screening tools developed to identify inhibitors of traditional drug targets do not perform as well when they are instead deployed at protein interaction sites. To address the need for novel inhibitors of important protein interactions, here we introduce an alternate docking strategy specifically designed for this regime. Our method, termed DARC (Docking Approach using Ray-Casting), matches the topography of a surface pocket "observed" from within the protein to the topography "observed" when viewing a potential ligand from the same vantage point. We applied DARC to carry out a virtual screen against the protein interaction site of human antiapoptotic protein Mcl-1 and found that four of the top-scoring 21 compounds showed clear inhibition in a biochemical assay. The Ki values for these compounds ranged from 1.2 to 21 µM, and each had ligand efficiency comparable to promising small-molecule inhibitors of other protein-protein interactions. These hit compounds do not resemble the natural (protein) binding partner of Mcl-1, nor do they resemble any known inhibitors of Mcl-1. Our results thus demonstrate the utility of DARC for identifying novel inhibitors of protein-protein interactions.

Natural product (-)-gossypol inhibits colon cancer cell growth by targeting RNA-binding protein Musashi-1.

Musashi-1 (MSI1) is an RNA-binding protein that acts as a translation activator or repressor of target mRNAs. The best-characterized MSI1 target is Numb mRNA, whose encoded protein negatively regulates Notch signaling. Additional MSI1 targets include the mRNAs for the tumor suppressor protein APC that regulates Wnt signaling and the cyclin-dependent kinase inhibitor P21(WAF-1). We hypothesized that increased expression of NUMB, P21 and APC, through inhibition of MSI1 RNA-binding activity might be an effective way to simultaneously downregulate Wnt and Notch signaling, thus blocking the growth of a broad range of cancer cells. We used a fluorescence polarization assay to screen for small molecules that disrupt the binding of MSI1 to its consensus RNA binding site. One of the top hits was (-)-gossypol (Ki = 476 ± 273 nM), a natural product from cottonseed, known to have potent anti-tumor activity and which has recently completed Phase IIb clinical trials for prostate cancer. Surface plasmon resonance and nuclear magnetic resonance studies demonstrate a direct interaction of (-)-gossypol with the RNA binding pocket of MSI1. We further showed that (-)-gossypol reduces Notch/Wnt signaling in several colon cancer cell lines having high levels of MSI1, with reduced SURVIVIN expression and increased apoptosis/autophagy. Finally, we showed that orally administered (-)-gossypol inhibits colon cancer growth in a mouse xenograft model. Our study identifies (-)-gossypol as a potential small molecule inhibitor of MSI1-RNA interaction, and suggests that inhibition of MSI1's RNA binding activity may be an effective anti-cancer strategy.

Sequence and structure signatures of cancer mutation hotspots in protein kinases.

Protein kinases are the most common protein domains implicated in cancer, where somatically acquired mutations are known to be functionally linked to a variety of cancers. Resequencing studies of protein kinase coding regions have emphasized the importance of sequence and structure determinants of cancer-causing kinase mutations in understanding of the mutation-dependent activation process. We have developed an integrated bioinformatics resource, which consolidated and mapped all currently available information on genetic modifications in protein kinase genes with sequence, structure and functional data. The integration of diverse data types provided a convenient framework for kinome-wide study of sequence-based and structure-based signatures of cancer mutations. The database-driven analysis has revealed a differential enrichment of SNPs categories in functional regions of the kinase domain, demonstrating that a significant number of cancer mutations could fall at structurally equivalent positions (mutational hotspots) within the catalytic core. We have also found that structurally conserved mutational hotspots can be shared by multiple kinase genes and are often enriched by cancer driver mutations with high oncogenic activity. Structural modeling and energetic analysis of the mutational hotspots have suggested a common molecular mechanism of kinase activation by cancer mutations, and have allowed to reconcile the experimental data. According to a proposed mechanism, structural effect of kinase mutations with a high oncogenic potential may manifest in a significant destabilization of the autoinhibited kinase form, which is likely to drive tumorigenesis at some level. Structure-based functional annotation and prediction of cancer mutation effects in protein kinases can facilitate an understanding of the mutation-dependent activation process and inform experimental studies exploring molecular pathology of tumorigenesis.