Tandem phosphorylation within an intrinsically disordered region regulates ACTN4 function

Phosphorylated residues occur preferentially in the intrinsically disordered regions of eukaryotic proteins. In the disordered amino-terminal region of human a-actinin-4 (ACTN4), Tyr[superscript 4] and Tyr[superscript 31] are phosphorylated in cells stimulated with epidermal growth factor (EGF), and a mutant with phosphorylation-mimicking mutations of both tyrosines exhibits reduced interaction with actin in vitro. Cleavage of ACTN4 by m-calpain, a protease that in motile cells is predominantly activated at the rear, removes the Tyr[superscript 4] site. We found that introducing a phosphomimetic mutation at only Tyr[superscript 31] was sufficient to inhibit the interaction with actin in vitro. However, molecular dynamics simulations predicted that Tyr[superscript 31] is mostly buried and that phosphorylation of Tyr[superscript 4] would increase the solvent exposure and thus kinase accessibility of Tyr[superscript 31]. In fibroblast cells, EGF stimulation increased tyrosine phosphorylation of a mutant form of ACTN4 with a phosphorylation-mimicking residue at Tyr[superscript 4], whereas a truncated mutant representing the product of m-calpain cleavage exhibited EGF-stimulated tyrosine phosphorylation at a background amount similar to that observed for a double phosphomimetic mutant of Tyr[superscript 4] and Tyr[superscript 31]. We also found that inhibition of the receptor tyrosine kinases of the TAM family, such as AXL, blocked EGF-stimulated tyrosine phosphorylation of ACTN4. Mathematical modeling predicted that the kinetics of phosphorylation at Tyr[superscript 31] can be dictated by the kinase affinity for Tyr[superscript 4]. This study suggests that tandem-site phosphorylation within intrinsically disordered regions provides a mechanism for a site to function as a switch to reveal a nearby function-regulating site.National Institutes of Health (U.S.) (Grant R01 GM69668

An integrated comparative phosphoproteomic and bioinformatic approach reveals a novel class of MPM-2 motifs upregulated in EGFRvIII-expressing Glioblastoma Cells

Glioblastoma (GBM, WHO grade IV) is an aggressively proliferative and invasive brain tumor that carries a poor clinical prognosis with a median survival of 9 to 12 months. In a prior phosphoproteomic study performed in the U87MG glioblastoma cell line, we identified tyrosine phosphorylation events that are regulated as a result of titrating EGFRvIII, a constitutively active mutant of the epidermal growth factor receptor (EGFR) associated with poor prognosis in GBM patients. In the present study, we have used the phosphoserine/phosphothreonine-specific antibody MPM-2 (mitotic protein monoclonal #2) to quantify serine/threonine phosphorylation events in the same cell lines. By employing a bioinformatic tool to identify amino acid sequence motifs regulated in response to increasing oncogene levels, a set of previously undescribed MPM-2 epitope sequence motifs orthogonal to the canonical “pS/pT-P” motif was identified. These motifs contain acidic amino acids in combinations of the −5, −2, +1, +3, and +5 positions relative to the phosphorylated amino acid. Phosphopeptides containing these motifs are upregulated in cells expressing EGFRvIII, raising the possibility of a general role for a previously unrecognized acidophilic kinase (e.g. casein kinase II (CK2)) in cell proliferation downstream of EGFR signaling.National Cancer Institute (U.S.). Integrative Cancer Biology Program (grant U54-CA112967)National Cancer Institute (U.S.). Bioengineering Research Partnership (grant R01-CA96504)National Institutes of Health (U.S.) (grant R01-GM60594

Robustness and applicability of transcription factor and pathway analysis tools on single-cell RNA-seq data

Many functional analysis tools have been developed to extract functional and mechanistic insight from bulk transcriptome data. With the advent of single-cell RNA sequencing (scRNA-seq), it is in principle possible to do such an analysis for single cells. However, scRNA-seq data has characteristics such as drop-out events and low library sizes. It is thus not clear if functional TF and pathway analysis tools established for bulk sequencing can be applied to scRNA-seq in a meaningful way.To address this question, we perform benchmark studies on simulated and real scRNA-seq data. We include the bulk-RNA tools PROGENy, GO enrichment, and DoRothEA that estimate pathway and transcription factor (TF) activities, respectively, and compare them against the tools SCENIC/AUCell and metaVIPER, designed for scRNA-seq. For the in silico study, we simulate single cells from TF/pathway perturbation bulk RNA-seq experiments. We complement the simulated data with real scRNA-seq data upon CRISPR-mediated knock-out. Our benchmarks on simulated and real data reveal comparable performance to the original bulk data. Additionally, we show that the TF and pathway activities preserve cell type-specific variability by analyzing a mixture sample sequenced with 13 scRNA-seq protocols. We also provide the benchmark data for further use by the community.Our analyses suggest that bulk-based functional analysis tools that use manually curated footprint gene sets can be applied to scRNA-seq data, partially outperforming dedicated single-cell tools. Furthermore, we find that the performance of functional analysis tools is more sensitive to the gene sets than to the statistic used

A genome-scale in vivo loss-of-function screen identifies Phf6 as a lineage-specific regulator of leukemia cell growth

We performed a genome-scale shRNA screen for modulators of B-cell leukemia progression in vivo. Results from this work revealed dramatic distinctions between the relative effects of shRNAs on the growth of tumor cells in culture versus in their native microenvironment. Specifically, we identified many “context-specific” regulators of leukemia development. These included the gene encoding the zinc finger protein Phf6. While inactivating mutations in PHF6 are commonly observed in human myeloid and T-cell malignancies, we found that Phf6 suppression in B-cell malignancies impairs tumor progression. Thus, Phf6 is a “lineage-specific” cancer gene that plays opposing roles in developmentally distinct hematopoietic malignancies.Massachusetts Institute of Technology. Department of Biology (Training Grant)National Cancer Institute (U.S.). Integrative Cancer Biology Program (U54-CA112967-06)National Institutes of Health (U.S.) (RO1-CA128803-05

Integrated mapping of pharmacokinetics and pharmacodynamics in a patient-derived xenograft model of glioblastoma

Therapeutic options for the treatment of glioblastoma remain inadequate despite concerted research efforts in drug development. Therapeutic failure can result from poor permeability of the blood-brain barrier, heterogeneous drug distribution, and development of resistance. Elucidation of relationships among such parameters could enable the development of predictive models of drug response in patients and inform drug development. Complementary analyses were applied to a glioblastoma patient-derived xenograft model in order to quantitatively map distribution and resulting cellular response to the EGFR inhibitor erlotinib. Mass spectrometry images of erlotinib were registered to histology and magnetic resonance images in order to correlate drug distribution with tumor characteristics. Phosphoproteomics and immunohistochemistry were used to assess protein signaling in response to drug, and integrated with transcriptional response using mRNA sequencing. This comprehensive dataset provides simultaneous insight into pharmacokinetics and pharmacodynamics and indicates that erlotinib delivery to intracranial tumors is insufficient to inhibit EGFR tyrosine kinase signaling.National Institutes of Health (U.S.) (U54 CA210180)MIT/Mayo Physical Sciences Center for Drug Distribution and Drug Efficacy in Brain TumorsDana-Farber Cancer Institute (PLGA Fund)Lundbeck FoundationNovo Nordisk Foundatio

Mass balance of the Greenland Ice Sheet from 1992 to 2018

In recent decades, the Greenland Ice Sheet has been a major contributor to global sea-level rise1,2, and it is expected to be so in the future3. Although increases in glacier flow4–6 and surface melting7–9 have been driven by oceanic10–12 and atmospheric13,14 warming, the degree and trajectory of today’s imbalance remain uncertain. Here we compare and combine 26 individual satellite measurements of changes in the ice sheet’s volume, flow and gravitational potential to produce a reconciled estimate of its mass balance. Although the ice sheet was close to a state of balance in the 1990s, annual losses have risen since then, peaking at 335 ± 62 billion tonnes per year in 2011. In all, Greenland lost 3,800 ± 339 billion tonnes of ice between 1992 and 2018, causing the mean sea level to rise by 10.6 ± 0.9 millimetres. Using three regional climate models, we show that reduced surface mass balance has driven 1,971 ± 555 billion tonnes (52%) of the ice loss owing to increased meltwater runoff. The remaining 1,827 ± 538 billion tonnes (48%) of ice loss was due to increased glacier discharge, which rose from 41 ± 37 billion tonnes per year in the 1990s to 87 ± 25 billion tonnes per year since then. Between 2013 and 2017, the total rate of ice loss slowed to 217 ± 32 billion tonnes per year, on average, as atmospheric circulation favoured cooler conditions15 and as ocean temperatures fell at the terminus of Jakobshavn Isbræ16. Cumulative ice losses from Greenland as a whole have been close to the IPCC’s predicted rates for their high-end climate warming scenario17, which forecast an additional 50 to 120 millimetres of global sea-level rise by 2100 when compared to their central estimate

Mass balance of the Greenland and Antarctic ice sheets from 1992 to 2020

Ice losses from the Greenland and Antarctic ice sheets have accelerated since the 1990s, accounting for a significant increase in the global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume, and in Earth's gravity field. Between 1992 and 2020, the ice sheets contributed 21.0±1.9g€¯mm to global mean sea level, with the rate of mass loss rising from 105g€¯Gtg€¯yr-1 between 1992 and 1996 to 372g€¯Gtg€¯yr-1 between 2016 and 2020. In Greenland, the rate of mass loss is 169±9g€¯Gtg€¯yr-1 between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86g€¯Gtg€¯yr-1 in 2017 to 444g€¯Gtg€¯yr-1 in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (82±9g€¯Gtg€¯yr-1) and, to a lesser extent, from the Antarctic Peninsula (13±5g€¯Gtg€¯yr-1). East Antarctica remains close to a state of balance, with a small gain of 3±15g€¯Gtg€¯yr-1, but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at 10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021)

Multiplexed relative and absolute quantitative immunopeptidomics reveals MHC I repertoire alterations induced by CDK4/6 inhibition

© 2020, The Author(s). Peptides bound to class I major histocompatibility complexes (MHC) play a critical role in immune cell recognition and can trigger an antitumor immune response in cancer. Surface MHC levels can be modulated by anticancer agents, altering immunity. However, understanding the peptide repertoire’s response to treatment remains challenging and is limited by quantitative mass spectrometry-based strategies lacking normalization controls. We describe an experimental platform that leverages recombinant heavy isotope-coded peptide MHCs (hipMHCs) and multiplex isotope tagging to quantify peptide repertoire alterations using low sample input. HipMHCs improve quantitative accuracy of peptide repertoire changes by normalizing for variation across analyses and enable absolute quantification using internal calibrants to determine copies per cell of MHC antigens, which can inform immunotherapy design. Applying this platform in melanoma cell lines to profile the immunopeptidome response to CDK4/6 inhibition and interferon-γ — known modulators of antigen presentation — uncovers treatment-specific alterations, connecting the intracellular response to extracellular immune presentation

VISAGE Reveals a Targetable Mitotic Spindle Vulnerability in Cancer Cells

There is an unmet need for new antimitotic drug combinations that target cancer-specific vulnerabilities. Based on our finding of elevated biomolecule oxidation in mitotically arrested cancer cells, we combined Plk1 inhibitors with TH588, an MTH1 inhibitor that prevents detoxification of oxidized nucleotide triphosphates. This combination showed robust synergistic killing of cancer, but not normal, cells that, surprisingly, was MTH1-independent. To dissect the underlying synergistic mechanism, we developed VISAGE, a strategy integrating experimental synergy quantification with computational-pathway-based gene expression analysis. VISAGE predicted, and we experimentally confirmed, that this synergistic combination treatment targeted the mitotic spindle. Specifically, TH588 binding to β-tubulin impaired microtubule assembly, which when combined with Plk1 blockade, synergistically disrupted mitotic chromosome positioning to the spindle midzone. These findings identify a cancer-specific mitotic vulnerability that is targetable using Plk1 inhibitors with microtubule-destabilizing agents and highlight the general utility of the VISAGE approach to elucidate molecular mechanisms of drug synergy.National Institutes of Health (U.S.) (Grant R35-ES028374)National Institutes of Health (U.S.) (Grant R01-GM104047)National Institutes of Health (U.S.) (Grant R01-ES015339)National Institutes of Health (U.S.) (Grant U54-CA112967)National Institutes of Health (U.S.) (Grant U54-CA217377)American Cancer Society. Post-Doctoral Fellowship (PF-13-355-01-TBE)National Cancer Institute (U.S.) (Support Grant P30-ES002109

Inflammatory but not mitogenic contexts prime synovial fibroblasts for compensatory signaling responses to p38 inhibition

Rheumatoid arthritis (RA) is a chronic inflammatory disorder that causes joint pain, swelling, and loss of function. Development of effective new drugs has proven challenging in part because of the complexities and interconnected nature of intracellular signaling networks that complicate the effects of pharmacological interventions. We characterized the kinase signaling pathways that are activated in RA and evaluated the multivariate effects of targeted inhibitors. Synovial fluids from RA patients activated the kinase signaling pathways JAK, JNK, p38, and MEK in synovial fibroblasts (SFs), a stromal cell type that promotes RA progression. Kinase inhibitors enhanced signaling of “off-target” pathways in a manner dependent on stimulatory context. Inhibitors of p38, which have been widely explored in clinical trials for RA, resulted in undesirable increases in nuclear factor B (NF-B), JNK, and MEK signaling in SFs in inflammatory, but not mitogenic, contexts. This was mediated by the transcription factor CREB, which functions in part within a negative feedback loop in MAPK signaling. CREB activation was induced predominately by p38 in response to inflammatory stimuli, but by MEK in response to mitogenic stimuli; hence, the effects of drugs targeting p38 or MEK were markedly different in SFs cultured under mitogenic or inflammatory conditions. Together, these findings illustrate how stimulatory context can alter dominance in pathway cross-talk even for a fixed network topology, thereby providing a rationale for why p38 inhibitors deliver limited benefits in RA and demonstrating the need for careful consideration of p38-targeted drugs in inflammation-related disorders.National Institutes of Health (U.S.) (Grant U54HL127365)National Institutes of Health (U.S.) (Grant P50GM107618)United States. Army Research Office (Grant W911NF-09-0001)National Institutes of Health (U.S.) (Grant 5F32AR062931