Patch-DCA: improved protein interface prediction by utilizing structural information and clustering DCA scores.

MOTIVATION: Over the past decade, there have been impressive advances in determining the 3D structures of protein complexes. However, there are still many complexes with unknown structures, even when the structures of the individual proteins are known. The advent of protein sequence information provides an opportunity to leverage evolutionary information to enhance the accuracy of protein-protein interface prediction. To this end, several statistical and machine learning methods have been proposed. In particular, direct coupling analysis has recently emerged as a promising approach for identification of protein contact maps from sequential information. However, the ability of these methods to detect protein-protein inter-residue contacts remains relatively limited. RESULTS: In this work, we propose a method to integrate sequential and co-evolution information with structural and functional information to increase the performance of protein-protein interface prediction. Further, we present a post-processing clustering method that improves the average relative F1 score by 70% and 24% and the average relative precision by 80% and 36% in comparison with two state-of-the-art methods, PSICOV and GREMLIN. AVAILABILITY AND IMPLEMENTATION: https://github.com/BioMLBoston/PatchDCA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The germline factor DDX4 contributes to the chemoresistance of small cell lung cancer cells.

Human cancers often re-express germline factors, yet their mechanistic role in oncogenesis and cancer progression remains unknown. Here we demonstrate that DEAD-box helicase 4 (DDX4), a germline factor and RNA helicase conserved in all multicellular organisms, contributes to increased cell motility and cisplatin-mediated drug resistance in small cell lung cancer (SCLC) cells. Proteomic analysis suggests that DDX4 expression upregulates proteins related to DNA repair and immune/inflammatory response. Consistent with these trends in cell lines, DDX4 depletion compromised in vivo tumor development while its overexpression enhanced tumor growth even after cisplatin treatment in nude mice. Further, the relatively higher DDX4 expression in SCLC patients correlates with decreased survival and shows increased expression of immune/inflammatory response markers. Taken together, we propose that DDX4 increases SCLC cell survival, by increasing the DNA damage and immune response pathways, especially under challenging conditions such as cisplatin treatment.

Cancer-associated fibroblasts are the main contributors to epithelial-to-mesenchymal signatures in the tumor microenvironment.

Epithelial-to-mesenchymal transition (EMT) is associated with tumor initiation, metastasis, and drug resistance. However, the mechanisms underlying these associations are largely unknown. We studied several tumor types to identify the source of EMT gene expression signals and a potential mechanism of resistance to immuno-oncology treatment. Across tumor types, EMT-related gene expression was strongly associated with expression of stroma-related genes. Based on RNA sequencing of multiple patient-derived xenograft models, EMT-related gene expression was enriched in the stroma versus parenchyma. EMT-related markers were predominantly expressed by cancer-associated fibroblasts (CAFs), cells of mesenchymal origin which produce a variety of matrix proteins and growth factors. Scores derived from a 3-gene CAF transcriptional signature (COL1A1, COL1A2, COL3A1) were sufficient to reproduce association between EMT-related markers and disease prognosis. Our results suggest that CAFs are the primary source of EMT signaling and have potential roles as biomarkers and targets for immuno-oncology therapies.

Regulation of neuroendocrine plasticity by the RNA-binding protein ZFP36L1.

Some small cell lung cancers (SCLCs) are highly sensitive to inhibitors of the histone demethylase LSD1. LSD1 inhibitors are thought to induce their anti-proliferative effects by blocking neuroendocrine differentiation, but the mechanisms by which LSD1 controls the SCLC neuroendocrine phenotype are not well understood. To identify genes required for LSD1 inhibitor sensitivity in SCLC, we performed a positive selection genome-wide CRISPR/Cas9 loss of function screen and found that ZFP36L1, an mRNA-binding protein that destabilizes mRNAs, is required for LSD1 inhibitor sensitivity. LSD1 binds and represses ZFP36L1 and upon LSD1 inhibition, ZFP36L1 expression is restored, which is sufficient to block the SCLC neuroendocrine differentiation phenotype and induce a non-neuroendocrine "inflammatory" phenotype. Mechanistically, ZFP36L1 binds and destabilizes SOX2 and INSM1 mRNAs, two transcription factors that are required for SCLC neuroendocrine differentiation. This work identifies ZFP36L1 as an LSD1 target gene that controls the SCLC neuroendocrine phenotype and demonstrates that modulating mRNA stability of lineage transcription factors controls neuroendocrine to non-neuroendocrine plasticity.

Plasticity in the Absence of NOTCH Uncovers a RUNX2-Dependent Pathway in Small Cell Lung Cancer.

Neuroendocrine to nonneuroendocrine plasticity supports small cell lung cancer (SCLC) tumorigenesis and promotes immunogenicity. Approximately 20% to 25% of SCLCs harbor loss-of-function (LOF) NOTCH mutations. Previous studies demonstrated that NOTCH functions as a SCLC tumor suppressor, but can also drive nonneuroendocrine plasticity to support SCLC growth. Given the dual functionality of NOTCH, it is not understood why SCLCs select for LOF NOTCH mutations and how these mutations affect SCLC tumorigenesis. In a CRISPR-based genetically engineered mouse model of SCLC, genetic loss of Notch1 or Notch2 modestly accelerated SCLC tumorigenesis. Interestingly, Notch-mutant SCLCs still formed nonneuroendocrine subpopulations, and these Notch-independent, nonneuroendocrine subpopulations were driven by Runx2-mediated regulation of Rest. Notch2-mutant nonneuroendocrine cells highly express innate immune signaling genes including stimulator of interferon genes (STING) and were sensitive to STING agonists. This work identifies a Notch-independent mechanism to promote nonneuroendocrine plasticity and suggests that therapeutic approaches to activate STING could be selectively beneficial for SCLCs with NOTCH2 mutations. SIGNIFICANCE: A genetically engineered mouse model of NOTCH-mutant SCLC reveals that nonneuroendocrine plasticity persists in the absence of NOTCH, driven by a RUNX2-REST-dependent pathway and innate immune signaling.

MET-Induced CD73 Restrains STING-Mediated Immunogenicity of EGFR-Mutant Lung Cancer.

Immunotherapy has shown limited efficacy in patients with EGFR-mutated lung cancer. Efforts to enhance the immunogenicity of EGFR-mutated lung cancer have been unsuccessful to date. Here, we discover that MET amplification, the most common mechanism of resistance to third-generation EGFR tyrosine kinase inhibitors (TKI), activates tumor cell STING, an emerging determinant of cancer immunogenicity (1). However, STING activation was restrained by ectonucleosidase CD73, which is induced in MET-amplified, EGFR-TKI-resistant cells. Systematic genomic analyses and cell line studies confirmed upregulation of CD73 in MET-amplified and MET-activated lung cancer contexts, which depends on coinduction of FOSL1. Pemetrexed (PEM), which is commonly used following EGFR-TKI treatment failure, was identified as an effective potentiator of STING-dependent TBK1-IRF3-STAT1 signaling in MET-amplified, EGFR-TKI-resistant cells. However, PEM treatment also induced adenosine production, which inhibited T-cell responsiveness. In an allogenic humanized mouse model, CD73 deletion enhanced immunogenicity of MET-amplified, EGFR-TKI-resistant cells, and PEM treatment promoted robust responses regardless of CD73 status. Using a physiologic antigen recognition model, inactivation of CD73 significantly increased antigen-specific CD8+ T-cell immunogenicity following PEM treatment. These data reveal that combined PEM and CD73 inhibition can co-opt tumor cell STING induction in TKI-resistant EGFR-mutated lung cancers and promote immunogenicity. SIGNIFICANCE: MET amplification upregulates CD73 to suppress tumor cell STING induction and T-cell responsiveness in TKI-resistant, EGFR-mutated lung cancer, identifying a strategy to enhance immunogenicity and improve treatment.

Diminished Efficacy of Programmed Death-(Ligand)1 Inhibition in STK11- and KEAP1-Mutant Lung Adenocarcinoma Is Affected by KRAS Mutation Status.

INTRODUCTION: STK11 and KEAP1 mutations (STK11 mutant [STK11MUT] and KEAP1MUT) are among the most often mutated genes in lung adenocarcinoma (LUAD). Although STK11MUT has been associated with resistance to programmed death-(ligand)1 (PD-[L]1) inhibition in KRASMUT LUAD, its impact on immunotherapy efficacy in KRAS wild-type (KRASWT) LUAD is currently unknown. Whether KEAP1MUT differentially affects outcomes to PD-(L)1 inhibition in KRASMUT and KRASWT LUAD is also unknown. METHODS: Clinicopathologic and genomic data were collected from September 2013 to September 2020 from patients with advanced LUAD at the Dana-Farber Cancer Institute/Massachusetts General Hospital cohort and the Memorial Sloan Kettering Cancer Center/MD Anderson Cancer Center cohort. Clinical outcomes to PD-(L)1 inhibition were analyzed according to KRAS, STK11, and KEAP1 mutation status in two independent cohorts. The Cancer Genome Atlas transcriptomic data were interrogated to identify differences in tumor gene expression and tumor immune cell subsets, respectively, according to KRAS/STK11 and KRAS/KEAP1 comutation status. RESULTS: In the combined cohort (Dana-Farber Cancer Institute/Massachusetts General Hospital + Memorial Sloan Kettering Cancer Center/MD Anderson Cancer Center) of 1261 patients (median age = 61 y [range: 22-92], 708 women [56.1%], 1065 smokers [84.4%]), KRAS mutations were detected in 536 cases (42.5%), and deleterious STK11 and KEAP1 mutations were found in 20.6% (260 of 1261) and 19.2% (231 of 1202) of assessable cases, respectively. In each independent cohort and in the combined cohort, STK11 and KEAP1 mutations were associated with significantly worse progression-free (STK11 hazard ratio [HR] = 2.04, p < 0.0001; KEAP1 HR = 2.05, p < 0.0001) and overall (STK11 HR = 2.09, p < 0.0001; KEAP1 HR = 2.24, p < 0.0001) survival to immunotherapy uniquely among KRASMUT but not KRASWT LUADs. Gene expression ontology and immune cell enrichment analyses revealed that the presence of STK11 or KEAP1 mutations results in distinct immunophenotypes in KRASMUT, but not in KRASWT, lung cancers. CONCLUSIONS: STK11 and KEAP1 mutations confer worse outcomes to immunotherapy among patients with KRASMUT but not among KRASWT LUAD. Tumors harboring concurrent KRAS/STK11 and KRAS/KEAP1 mutations display distinct immune profiles in terms of gene expression and immune cell infiltration.

The Extracellular Milieu of Toxoplasma's Lytic Cycle Drives Lab Adaptation, Primarily by Transcriptional Reprogramming.

Evolve and resequencing (E&R) was applied to lab adaptation of Toxoplasma gondii for over 1,500 generations with the goal of mapping host-independent in vitro virulence traits. Phenotypic assessments of steps across the lytic cycle revealed that only traits needed in the extracellular milieu evolved. Nonsynonymous single-nucleotide polymorphisms (SNPs) in only one gene, a P4 flippase, fixated across two different evolving populations, whereas dramatic changes in the transcriptional signature of extracellular parasites were identified. Newly developed computational tools correlated phenotypes evolving at different rates with specific transcriptomic changes. A set of 300 phenotype-associated genes was mapped, of which nearly 50% is annotated as hypothetical. Validation of a select number of genes by knockouts confirmed their role in lab adaptation and highlights novel mechanisms underlying in vitro virulence traits. Further analyses of differentially expressed genes revealed the development of a "pro-tachyzoite" profile as well as the upregulation of the fatty acid biosynthesis (FASII) pathway. The latter aligned with the P4 flippase SNP and aligned with a low abundance of medium-chain fatty acids at low passage, indicating this is a limiting factor in extracellular parasites. In addition, partial overlap with the bradyzoite differentiation transcriptome in extracellular parasites indicated that stress pathways are involved in both situations. This was reflected in the partial overlap between the assembled ApiAP2 and Myb transcription factor network underlying the adapting extracellular state with the bradyzoite differentiation program. Overall, E&R is a new genomic tool successfully applied to map the development of polygenic traits underlying in vitro virulence of T. gondii. IMPORTANCE It has been well established that prolonged in vitro cultivation of Toxoplasma gondii augments progression of the lytic cycle. This lab adaptation results in increased capacities to divide, migrate, and survive outside a host cell, all of which are considered host-independent virulence factors. However, the mechanistic basis underlying these enhanced virulence features is unknown. Here, E&R was utilized to empirically characterize the phenotypic, genomic, and transcriptomic changes in the non-lab-adapted strain, GT1, during 2.5 years of lab adaptation. This identified the shutdown of stage differentiation and upregulation of lipid biosynthetic pathways as the key processes being modulated. Furthermore, lab adaptation was primarily driven by transcriptional reprogramming, which rejected the starting hypothesis that genetic mutations would drive lab adaptation. Overall, the work empirically shows that lab adaptation augments T. gondii's in vitro virulence by transcriptional reprogramming and that E&R is a powerful new tool to map multigenic traits.

Dynamic single-cell RNA sequencing identifies immunotherapy persister cells following PD-1 blockade.

Resistance to oncogene-targeted therapies involves discrete drug-tolerant persister cells, originally discovered through in vitro assays. Whether a similar phenomenon limits efficacy of programmed cell death 1 (PD-1) blockade is poorly understood. Here, we performed dynamic single-cell RNA-Seq of murine organotypic tumor spheroids undergoing PD-1 blockade, identifying a discrete subpopulation of immunotherapy persister cells (IPCs) that resisted CD8+ T cell-mediated killing. These cells expressed Snai1 and stem cell antigen 1 (Sca-1) and exhibited hybrid epithelial-mesenchymal features characteristic of a stem cell-like state. IPCs were expanded by IL-6 but were vulnerable to TNF-α-induced cytotoxicity, relying on baculoviral IAP repeat-containing protein 2 (Birc2) and Birc3 as survival factors. Combining PD-1 blockade with Birc2/3 antagonism in mice reduced IPCs and enhanced tumor cell killing in vivo, resulting in durable responsiveness that matched TNF cytotoxicity thresholds in vitro. Together, these data demonstrate the power of high-resolution functional ex vivo profiling to uncover fundamental mechanisms of immune escape from durable anti-PD-1 responses, while identifying IPCs as a cancer cell subpopulation targetable by specific therapeutic combinations.

Intrinsic Immunogenicity of Small Cell Lung Carcinoma Revealed by Its Cellular Plasticity.

Small cell lung carcinoma (SCLC) is highly mutated, yet durable response to immune checkpoint blockade (ICB) is rare. SCLC also exhibits cellular plasticity, which could influence its immunobiology. Here we discover that a distinct subset of SCLC uniquely upregulates MHC I, enriching for durable ICB benefit. In vitro modeling confirms epigenetic recovery of MHC I in SCLC following loss of neuroendocrine differentiation, which tracks with derepression of STING. Transient EZH2 inhibition expands these nonneuroendocrine cells, which display intrinsic innate immune signaling and basally restored antigen presentation. Consistent with these findings, murine nonneuroendocrine SCLC tumors are rejected in a syngeneic model, with clonal expansion of immunodominant effector CD8 T cells. Therapeutically, EZH2 inhibition followed by STING agonism enhances T-cell recognition and rejection of SCLC in mice. Together, these data identify MHC I as a novel biomarker of SCLC immune responsiveness and suggest novel immunotherapeutic approaches to co-opt SCLC's intrinsic immunogenicity. SIGNIFICANCE: SCLC is poorly immunogenic, displaying modest ICB responsiveness with rare durable activity. In profiling its plasticity, we uncover intrinsically immunogenic MHC Ihi subpopulations of nonneuroendocrine SCLC associated with durable ICB benefit. We also find that combined EZH2 inhibition and STING agonism uncovers this cell state, priming cells for immune rejection.This article is highlighted in the In This Issue feature, p. 1861.

Treatment-Induced Tumor Dormancy through YAP-Mediated Transcriptional Reprogramming of the Apoptotic Pathway.

Eradicating tumor dormancy that develops following epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor (TKI) treatment of EGFR-mutant non-small cell lung cancer, is an attractive therapeutic strategy but the mechanisms governing this process are poorly understood. Blockade of ERK1/2 reactivation following EGFR TKI treatment by combined EGFR/MEK inhibition uncovers cells that survive by entering a senescence-like dormant state characterized by high YAP/TEAD activity. YAP/TEAD engage the epithelial-to-mesenchymal transition transcription factor SLUG to directly repress pro-apoptotic BMF, limiting drug-induced apoptosis. Pharmacological co-inhibition of YAP and TEAD, or genetic deletion of YAP1, all deplete dormant cells by enhancing EGFR/MEK inhibition-induced apoptosis. Enhancing the initial efficacy of targeted therapies could ultimately lead to prolonged treatment responses in cancer patients.