All Publications


  • Leveraging nonstructural data to predict structures and affinities of protein-ligand complexes. Proceedings of the National Academy of Sciences of the United States of America Paggi, J. M., Belk, J. A., Hollingsworth, S. A., Villanueva, N., Powers, A. S., Clark, M. J., Chemparathy, A. G., Tynan, J. E., Lau, T. K., Sunahara, R. K., Dror, R. O. 1800; 118 (51)

    Abstract

    Over the past five decades, tremendous effort has been devoted to computational methods for predicting properties of ligands-i.e., molecules that bind macromolecular targets. Such methods, which are critical to rational drug design, fall into two categories: physics-based methods, which directly model ligand interactions with the target given the target's three-dimensional (3D) structure, and ligand-based methods, which predict ligand properties given experimental measurements for similar ligands. Here, we present a rigorous statistical framework to combine these two sources of information. We develop a method to predict a ligand's pose-the 3D structure of the ligand bound to its target-that leverages a widely available source of information: a list of other ligands that are known to bind the same target but for which no 3D structure is available. This combination of physics-based and ligand-based modeling improves pose prediction accuracy across all major families of drug targets. Using the same framework, we develop a method for virtual screening of drug candidates, which outperforms standard physics-based and ligand-based virtual screening methods. Our results suggest broad opportunities to improve prediction of various ligand properties by combining diverse sources of information through customized machine-learning approaches.

    View details for DOI 10.1073/pnas.2112621118

    View details for PubMedID 34921117

  • ecDNA hubs drive cooperative intermolecular oncogene expression. Nature Hung, K. L., Yost, K. E., Xie, L., Shi, Q., Helmsauer, K., Luebeck, J., Schopflin, R., Lange, J. T., Chamorro Gonzalez, R., Weiser, N. E., Chen, C., Valieva, M. E., Wong, I. T., Wu, S., Dehkordi, S. R., Duffy, C. V., Kraft, K., Tang, J., Belk, J. A., Rose, J. C., Corces, M. R., Granja, J. M., Li, R., Rajkumar, U., Friedlein, J., Bagchi, A., Satpathy, A. T., Tjian, R., Mundlos, S., Bafna, V., Henssen, A. G., Mischel, P. S., Liu, Z., Chang, H. Y. 2021

    Abstract

    Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs-clusters of around 10-100 ecDNAs within the nucleus-enable intermolecular enhancer-gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer-gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy.

    View details for DOI 10.1038/s41586-021-04116-8

    View details for PubMedID 34819668

  • A human mutation in STAT3 promotes type 1 diabetes through a defect in CD8+ T cell tolerance. The Journal of experimental medicine Warshauer, J. T., Belk, J. A., Chan, A. Y., Wang, J., Gupta, A. R., Shi, Q., Skartsis, N., Peng, Y., Phipps, J. D., Acenas, D., Smith, J. A., Tamaki, S. J., Tang, Q., Gardner, J. M., Satpathy, A. T., Anderson, M. S. 2021; 218 (8)

    Abstract

    Naturally occurring cases of monogenic type 1 diabetes (T1D) help establish direct mechanisms driving this complex autoimmune disease. A recently identified de novo germline gain-of-function (GOF) mutation in the transcriptional regulator STAT3 was found to cause neonatal T1D. We engineered a novel knock-in mouse incorporating this highly diabetogenic human STAT3 mutation (K392R) and found that these mice recapitulated the human autoimmune diabetes phenotype. Paired single-cell TCR and RNA sequencing revealed that STAT3-GOF drives proliferation and clonal expansion of effector CD8+ cells that resist terminal exhaustion. Single-cell ATAC-seq showed that these effector T cells are epigenetically distinct and have differential chromatin architecture induced by STAT3-GOF. Analysis of islet TCR clonotypes revealed a CD8+ cell reacting against known antigen IGRP, and STAT3-GOF in an IGRP-reactive TCR transgenic model demonstrated that STAT3-GOF intrinsic to CD8+ cells is sufficient to accelerate diabetes onset. Altogether, these findings reveal a diabetogenic CD8+ T cell response that is restrained in the presence of normal STAT3 activity and drives diabetes pathogenesis.

    View details for DOI 10.1084/jem.20210759

    View details for PubMedID 34115115

  • Differential usage of transcriptional repressor Zeb2 enhancers distinguishes adult and embryonic hematopoiesis. Immunity Huang, X., Ferris, S. T., Kim, S., Choudhary, M. N., Belk, J. A., Fan, C., Qi, Y., Sudan, R., Xia, Y., Desai, P., Chen, J., Ly, N., Shi, Q., Bagadia, P., Liu, T., Guilliams, M., Egawa, T., Colonna, M., Diamond, M. S., Murphy, T. L., Satpathy, A. T., Wang, T., Murphy, K. M. 2021

    Abstract

    The transcriptional repressor ZEB2 regulates development of many cell fates among somatic, neural, and hematopoietic lineages, but the basis for its requirement in these diverse lineages is unclear. Here, we identified a 400-basepair (bp) region located 165 kilobases (kb) upstream of the Zeb2 transcriptional start site (TSS) that binds the E proteins at several E-box motifs and was active in hematopoietic lineages. Germline deletion of this 400-bp region (Zeb2Delta-165mice) specifically prevented Zeb2 expression in hematopoietic stem cell (HSC)-derived lineages. Zeb2Delta-165 mice lacked development of plasmacytoid dendritic cells (pDCs), monocytes, and B cells. All macrophages in Zeb2Delta-165 mice were exclusively of embryonic origin. Using single-cell chromatin profiling, we identified a second Zeb2 enhancer located at+164-kb that was selectively active in embryonically derived lineages, but not HSC-derived ones. Thus, Zeb2 expression in adult, but not embryonic, hematopoiesis is selectively controlled by the -165-kb Zeb2 enhancer.

    View details for DOI 10.1016/j.immuni.2021.04.015

    View details for PubMedID 34004142

  • The Role of Aire in the selection of Regulatory T cells in Type 1 Diabetes Bridge, J., Balolong, J., Belk, J., Dong, S., Cramer, N., Mowery, C., Crawford, F., Ye, J., Satpathy, A., Kappler, J. W., Bluestone, J. A., Anderson, M. S. AMER ASSOC IMMUNOLOGISTS. 2021
  • Transient rest restores functionality in exhausted CAR-T cells through epigenetic remodeling. Science (New York, N.Y.) Weber, E. W., Parker, K. R., Sotillo, E., Lynn, R. C., Anbunathan, H., Lattin, J., Good, Z., Belk, J. A., Daniel, B., Klysz, D., Malipatlolla, M., Xu, P., Bashti, M., Heitzeneder, S., Labanieh, L., Vandris, P., Majzner, R. G., Qi, Y., Sandor, K., Chen, L., Prabhu, S., Gentles, A. J., Wandless, T. J., Satpathy, A. T., Chang, H. Y., Mackall, C. L. 2021; 372 (6537)

    Abstract

    T cell exhaustion limits immune responses against cancer and is a major cause of resistance to chimeric antigen receptor (CAR)-T cell therapeutics. Using murine xenograft models and an in vitro model wherein tonic CAR signaling induces hallmark features of exhaustion, we tested the effect of transient cessation of receptor signaling, or rest, on the development and maintenance of exhaustion. Induction of rest through enforced down-regulation of the CAR protein using a drug-regulatable system or treatment with the multikinase inhibitor dasatinib resulted in the acquisition of a memory-like phenotype, global transcriptional and epigenetic reprogramming, and restored antitumor functionality in exhausted CAR-T cells. This work demonstrates that rest can enhance CAR-T cell efficacy by preventing or reversing exhaustion, and it challenges the notion that exhaustion is an epigenetically fixed state.

    View details for DOI 10.1126/science.aba1786

    View details for PubMedID 33795428

  • Discovery and functional interrogation of SARS-CoV-2 RNA-host protein interactions. Cell Flynn, R. A., Belk, J. A., Qi, Y., Yasumoto, Y., Wei, J., Alfajaro, M. M., Shi, Q., Mumbach, M. R., Limaye, A., DeWeirdt, P. C., Schmitz, C. O., Parker, K. R., Woo, E., Chang, H. Y., Horvath, T. L., Carette, J. E., Bertozzi, C. R., Wilen, C. B., Satpathy, A. T. 2021

    Abstract

    SARS-CoV-2 is the cause of a pandemic with growing global mortality. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with ChIRP-MS data from three other RNA viruses defined viral specificity of RNA-host protein interactions. Targeted CRISPR screens revealed that the majority of functional RNA-binding proteins protect the host from virus-induced cell death, and comparative CRISPR screens across seven RNA viruses revealed shared and SARS-specific antiviral factors. Finally, by combining the RNA-centric approach and functional CRISPR screens, we demonstrated a physical and functional connection between SARS-CoV-2 and mitochondria, highlighting this organelle as a general platform for antiviral activity. Altogether, these data provide a comprehensive catalog of functional SARS-CoV-2 RNA-host protein interactions, which may inform studies to understand the host-virus interface and nominate host pathways that could be targeted for therapeutic benefit.

    View details for DOI 10.1016/j.cell.2021.03.012

    View details for PubMedID 33743211

  • Affinity-Restricted Memory B Cells Dominate Recall Responses to Heterologous Flaviviruses. Immunity Wong, R., Belk, J. A., Govero, J., Uhrlaub, J. L., Reinartz, D., Zhao, H., Errico, J. M., D'Souza, L., Ripperger, T. J., Nikolich-Zugich, J., Shlomchik, M. J., Satpathy, A. T., Fremont, D. H., Diamond, M. S., Bhattacharya, D. 2020

    Abstract

    Memory B cells (MBCs) can respond to heterologous antigens either by molding new specificities through secondary germinal centers (GCs) or by selecting preexisting clones without further affinity maturation. To distinguish these mechanisms in flavivirus infections and immunizations, we studied recall responses to envelope protein domain III (DIII). Conditional deletion of activation-induced cytidine deaminase (AID) between heterologous challenges of West Nile, Japanese encephalitis, Zika, and dengue viruses did not affect recall responses. DIII-specific MBCs were contained mostly within the plasma-cell-biased CD80+ subset, and few GCs arose following heterologous boosters, demonstrating that recall responses are confined by preexisting clonal diversity. Measurement of monoclonal antibody (mAb) binding affinity to DIII proteins, timed AID deletion, single-cell RNA sequencing, and lineage tracing experiments point to selection of relatively low-affinity MBCs as a mechanism to promote diversity. Engineering immunogens to avoid this MBC diversity may facilitate flavivirus-type-specific vaccines with minimized potential for infection enhancement.

    View details for DOI 10.1016/j.immuni.2020.09.001

    View details for PubMedID 33010224

  • Systematic discovery and functional interrogation of SARS-CoV-2 viral RNA-host protein interactions during infection. bioRxiv : the preprint server for biology Flynn, R. A., Belk, J. A., Qi, Y. n., Yasumoto, Y. n., Schmitz, C. O., Mumbach, M. R., Limaye, A. n., Wei, J. n., Alfajaro, M. M., Parker, K. R., Chang, H. Y., Horvath, T. L., Carette, J. E., Bertozzi, C. n., Wilen, C. B., Satpathy, A. T. 2020

    Abstract

    Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of a pandemic with growing global mortality. There is an urgent need to understand the molecular pathways required for host infection and anti-viral immunity. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with viral ChIRP-MS data from three other positive-sense RNA viruses defined pan-viral and SARS-CoV-2-specific host interactions. Functional interrogation of these factors with a genome-wide CRISPR screen revealed that the vast majority of viral RNA-binding proteins protect the host from virus-induced cell death, and we identified known and novel anti-viral proteins that regulate SARS-CoV-2 pathogenicity. Finally, our RNA-centric approach demonstrated a physical connection between SARS-CoV-2 RNA and host mitochondria, which we validated with functional and electron microscopy data, providing new insights into a more general virus-specific protein logic for mitochondrial interactions. Altogether, these data provide a comprehensive catalogue of SARS-CoV-2 RNA-host protein interactions, which may inform future studies to understand the mechanisms of viral pathogenesis, as well as nominate host pathways that could be targeted for therapeutic benefit.· ChIRP-MS of SARS-CoV-2 RNA identifies a comprehensive viral RNA-host protein interaction network during infection across two species· Comparison to RNA-protein interaction networks with Zika virus, dengue virus, and rhinovirus identify SARS-CoV-2-specific and pan-viral RNA protein complexes and highlights distinct intracellular trafficking pathways· Intersection of ChIRP-MS and genome-wide CRISPR screens identify novel SARS-CoV-2-binding proteins with pro- and anti-viral function· Viral RNA-RNA and RNA-protein interactions reveal specific SARS-CoV-2-mediated mitochondrial dysfunction during infection.

    View details for DOI 10.1101/2020.10.06.327445

    View details for PubMedID 33052334

    View details for PubMedCentralID PMC7553159