Early Immune Changes Support Signet Ring Cell Dormancy in CDH1-Driven Hereditary Diffuse Gastric Carcinogenesis

Stage IA gastric adenocarcinoma, characterized by foci of intramucosal signet ring cells (SRC), is found in nearly all asymptomatic patients with germline pathogenic CDH1 variants and hereditary diffuse gastric cancer syndrome (HDGC). The molecular steps involved in initiating malignant transformation and promoting SRC dormancy in HDGC are unknown. Here, whole-exome bulk RNA sequencing (RNA-seq) of SRCs and adjacent non-SRC epithelium (NEP) was performed on laser-capture microdissected (LCM) regions of interest found in risk-reducing total gastrectomy specimens from patients with HDGC (Clinicaltrials.gov ID: NCT03030404). In total, 20 patients (6 male, 14 female) with confirmed HDGC were identified. Analysis of differentially expressed genes (DEG) demonstrated upregulation of certain individual EMT and proliferation genes. However, no oncogenic pathways were found to be upregulated in SRCs. Rather, SRC regions had significant enrichment in pathways involved in T-cell signaling. CIBERSORTx predicted significant increases in the presence of regulatory T cells (Treg) specific to SRC regions. IHC confirmed an increase in FOXP3+ cells in SRC foci, as well as elevations in CD4+ T cells and HLA-DR staining. In summary, the tumor immune microenvironment is microscopically inseparable from stage IA gastric SRCs using a granular isolation technique. An elevation in CD4+ T cells within SRC regions correlates with clinically observed SRC dormancy, while Treg upregulation represents a potential immune escape mechanism.

IMPLICATIONS

Characterization of the tumor-immune microenvironment in HDGC underscores the potential for the immune system to shape the transcriptional profile of the earliest tumors, which suggests immune-directed therapy as a potential cancer interception strategy in diffuse-type gastric cancer.

The Fanconi anemia core complex promotes CtIP-dependent end-resection to drive homologous recombination at DNA double-strand breaks

Homologous Recombination (HR) is a high-fidelity repair mechanism of DNA Double-Strand Breaks (DSBs), which are induced by irradiation, genotoxic chemicals or physiological DNA damaging processes. DSBs are also generated as intermediates during the repair of interstrand crosslinks (ICLs). In this context, the Fanconi anemia (FA) core complex, which is effectively recruited to ICLs, promotes HR-mediated DSB-repair. However, whether the FA core complex also promotes HR at ICL-independent DSBs remains controversial. Here, we identified the FA core complex members FANCL and Ube2T as HR-promoting factors in a CRISPR/Cas9-based screen with cells carrying the DSB-repair reporter DSB-Spectrum. Using isogenic cell-line models, we validated the HR-function of FANCL and Ube2T, and demonstrated a similar function for their ubiquitination-substrate FANCD2. We further show that FANCL and Ube2T are directly recruited to DSBs and are required for the accumulation of FANCD2 at these break sites. Mechanistically, we demonstrate that FANCL ubiquitin ligase activity is required for the accumulation of the nuclease CtIP at DSBs, and consequently for optimal end-resection and Rad51 loading. CtIP overexpression rescues HR in FANCL-deficient cells, validating that FANCL primarily regulates HR by promoting CtIP recruitment. Together, these data demonstrate that the FA core complex and FANCD2 have a dual genome maintenance function by promoting repair of DSBs as well as the repair of ICLs.

Biphasic JNK-Erk signaling separates the induction and maintenance of cell senescence after DNA damage induced by topoisomerase II inhibition

Genotoxic stress in mammalian cells, including those caused by anti-cancer chemotherapy, can induce temporary cell-cycle arrest, DNA damage-induced senescence (DDIS), or apoptotic cell death. Despite obvious clinical importance, it is unclear how the signals emerging from DNA damage are integrated together with other cellular signaling pathways monitoring the cell's environment and/or internal state to control different cell fates. Using single-cell-based signaling measurements combined with tensor partial least square regression (t-PLSR)/principal component analysis (PCA) analysis, we show that JNK and Erk MAPK signaling regulates the initiation of cell senescence through the transcription factor AP-1 at early times after doxorubicin-induced DNA damage and the senescence-associated secretory phenotype (SASP) at late times after damage. These results identify temporally distinct roles for signaling pathways beyond the classic DNA damage response (DDR) that control the cell senescence decision and modulate the tumor microenvironment and reveal fundamental similarities between signaling pathways responsible for oncogene-induced senescence (OIS) and senescence caused by topoisomerase II inhibition. A record of this paper's transparent peer review process is included in the supplemental information.

Neutrophil heterogeneity and emergence of a distinct population of CD11b/CD18-activated low-density neutrophils after trauma

INTRODUCTION

Multiple large clinical trauma trials have documented an increased susceptibility to infection after injury. Although neutrophils (polymorphonuclear leukocytes [PMNs]) were historically considered a homogeneous cell type, we hypothesized that injury could alter neutrophil heterogeneity and predispose to dysfunction. To explore whether trauma modifies PMN heterogeneity, we performed an observational mass-spectrometry-based cytometry study on total leukocytes and low-density PMNs found in the peripheral blood mononuclear cell fraction of leukocytes from healthy controls and trauma patients.

METHODS

A total of 74 samples from 12 trauma patients, each sampled at 1 or more time points, and matched controls were fractionated and profiled by mass-spectrometry-based cytometry using a panel of 44 distinct markers. After deconvolution and conservative gating on neutrophils, data were analyzed using Seurat, followed by clustering of principal components.

RESULTS

Eleven distinct neutrophil populations were resolved in control and trauma neutrophils based on differential protein surface marker expression. Trauma markedly altered the basal heterogeneity of neutrophil subgroups seen in the control samples, with loss of a dominant population of resting neutrophils marked by high expression of C3AR and low levels of CD63, CD64, and CD177 (cluster 1), and expansion of two alternative neutrophil populations, one of which is marked by high expression of CD177 with suppression of CD10, CD16, C3AR, CD63, and CD64 (cluster 6). Remarkably, following trauma, a substantially larger percentage of neutrophils sediment in the monocyte fraction. These low-density neutrophils bear markers of functional exhaustion and form a unique trauma-induced population (cluster 9) with markedly upregulated expression of active surface adhesion molecules (activated CD11b/CD18), with suppression of nearly all other surface markers, including receptors for formyl peptides, leukotrienes, chemokines, and complement.

CONCLUSION

Circulating neutrophils demonstrate considerable evidence of functional heterogeneity that is markedly altered by trauma. Trauma induces evolution of a novel, exhausted, low-density neutrophil population with immunosuppressive features.

Screening for CD19-specific chimaeric antigen receptors with enhanced signalling via a barcoded library of intracellular domains

The immunostimulatory intracellular domains (ICDs) of chimaeric antigen receptors (CARs) are essential for converting antigen recognition into antitumoural function. Although there are many possible combinations of ICDs, almost all current CARs rely on combinations of CD3𝛇, CD28 and 4-1BB. Here we show that a barcoded library of 700,000 unique CD19-specific CARs with diverse ICDs cloned into lentiviral vectors and transduced into Jurkat T cells can be screened at high throughput via cell sorting and next-generation sequencing to optimize CAR signalling for antitumoural functions. By using this screening approach, we identified CARs with new ICD combinations that, compared with clinically available CARs, endowed human primary T cells with comparable tumour control in mice and with improved proliferation, persistence, exhaustion and cytotoxicity after tumour rechallenge in vitro. The screening strategy can be adapted to other disease models, cell types and selection conditions, and could be used to improve adoptive cell therapies and to expand their utility to new disease indications.

Proteomics of Coagulopathy Following Injury Reveals Limitations of Using Laboratory Assessment to Define Trauma-Induced Coagulopathy to Predict Massive Transfusion

Objective

Trauma-induced coagulopathy (TIC) is provoked by multiple mechanisms and is perceived to be one driver of massive transfusions (MT). Single laboratory values using prothrombin time (INR) or thrombelastography (TEG) are used to clinically define this complex process. We used a proteomics approach to test whether current definitions of TIC (INR, TEG, or clinical judgement) are sufficient to capture the majority of protein changes associated with MT.

Methods

Eight level-I trauma centers contributed blood samples from patients available early after injury. TIC was defined as INR >1.5 (INR-TIC), TEG maximum amplitude <50mm (TEG-TIC), or clinical judgement (Clin-TIC) by the trauma surgeon. MT was defined as > 10 units of red blood cells in 24 hours or > 4 units RBC/hour during the first 4 hr. SomaLogic proteomic analysis of 1,305 proteins was performed. Pathways associated with proteins dysregulated in patients with each TIC definition and MT were identified.

Results

Patients (n=211) had a mean injury severity score of 24, with a MT and mortality rate of 22% and 12%, respectively. We identified 578 SOMAscan analytes dysregulated among MT patients, of which INR-TIC, TEG-TIC, and Clin-TIC patients showed dysregulation only in 25%, 3%, and 4% of these, respectively. TIC definitions jointly failed to show changes in 73% of the protein levels associated with MT, and failed to identify 26% of patients that received a massive transfusion. INR-TIC and TEG-TIC patients showed dysregulation of proteins significantly associated with complement activity. Proteins dysregulated in Clin-TIC or massive transfusion patients were not significantly associated with any pathway.

Conclusion

These data indicate there are unexplored opportunities to identify patients at risk for massive bleeding. Only a small subset of proteins that are dysregulated in patients receiving MT are statistically significantly dysregulated among patients whose TIC is defined based solely on laboratory measurements or clinical assessment.

Cancer-Associated Fibroblasts and Squamous Epithelial Cells Constitute a Unique Microenvironment in a Mouse Model of Inflammation-Induced Colon Cancer

The tumor microenvironment plays a key role in the pathogenesis of colorectal tumors and contains various cell types including epithelial, immune, and mesenchymal cells. Characterization of the interactions between these cell types is necessary for revealing the complex nature of tumors. In this study, we used single-cell RNA-seq (scRNA-seq) to compare the tumor microenvironments between a mouse model of sporadic colorectal adenoma (Lrig1^CreERT2/+;Apc^2lox14/+) and a mouse model of inflammation-driven colorectal cancer induced by azoxymethane and dextran sodium sulfate (AOM/DSS). While both models develop tumors in the distal colon, we found that the two tumor types have distinct microenvironments. AOM/DSS tumors have an increased abundance of two populations of cancer-associated fibroblasts (CAFs) compared with APC tumors, and we revealed their divergent spatial association with tumor cells using multiplex immunofluorescence (MxIF) imaging. We also identified a unique squamous cell population in AOM/DSS tumors, whose origins were distinct from anal squamous epithelial cells. These cells were in higher proportions upon administration of a chemotherapy regimen of 5-Fluorouracil/Irinotecan. We used computational inference algorithms to predict cell-cell communication mediated by ligand-receptor interactions and downstream pathway activation, and identified potential mechanistic connections between CAFs and tumor cells, as well as CAFs and squamous epithelial cells. This study provides important preclinical insight into the microenvironment of two distinct models of colorectal tumors and reveals unique roles for CAFs and squamous epithelial cells in the AOM/DSS model of inflammation-driven cancer.

An in vitro intestinal model captures immunomodulatory properties of the microbiota in inflammation

Considerable effort has been put forth to understand mechanisms by which the microbiota modulates and responds to inflammation. Here, we explored whether oxidation metabolites produced by the host during inflammation, sodium nitrate and trimethylamine oxide, impact the composition of a human stool bacterial population in a gut simulator. We then assessed whether an immune-competent in vitro intestinal model responded differently to spent medium from bacteria exposed to these cues compared to spent medium from a control bacterial population. The host-derived oxidation products were found to decrease levels of Bacteroidaceae and overall microbiota metabolic potential, while increasing levels of proinflammatory Enterobacteriaceae and lipopolysaccharide in bacterial cultures, reflecting shifts that occur in vivo in inflammation. Spent microbiota media induced elevated intracellular mucin levels and reduced intestinal monolayer integrity as reflected in transepithelial electrical resistance relative to fresh medium controls. However, multiplexed cytokine analysis revealed markedly different cytokine signatures from intestinal cultures exposed to spent medium with added oxidation products relative to spent control medium, while cytokine signatures of cultures exposed to fresh media were similar regardless of addition of host-derived cues. Further, the presence of immune cells in the intestinal model was required for this differentiation of cytokine signatures. This study indicates that simple in vitro immune-competent intestinal models can capture bacterial-mammalian cross-talk in response to host-derived oxidation products and supports utility of these systems for mechanistic studies of interactions between the gut microbiome and host in inflammation.

The injury response to DNA damage in live tumor cells promotes antitumor immunity

Although immune checkpoint blockade (ICB) has strong clinical benefit for treating some tumor types, it fails in others, indicating a need for additional modalities to enhance the ICB effect. Here, we identified one such modality by using DNA damage to create a live, injured tumor cell adjuvant. Using an optimized ex vivo coculture system, we found that treating tumor cells with specific concentrations of etoposide, mitoxantrone, or doxorubicin markedly enhanced dendritic cell–mediated T cell activation. These immune-enhancing effects of DNA damage did not correlate with immunogenic cell death markers or with the extent of apoptosis or necroptosis; instead, these effects were mediated by live injured cells with activation of the DNA-PK, ATR, NF-κB, p38 MAPK, and RIPK1 signaling pathways. In mice, intratumoral injection of ex vivo etoposide–treated tumor cells in combination with systemic ICB (by anti-PD-1 and anti-CTLA4 antibodies) increased the number of intratumoral CD103+ dendritic cells and circulating tumor-antigen–specific CD8+ T cells, decreased tumor growth, and improved survival. These effects were absent in Batf3−/− mice and in mice in which the DNA-damaging drug was injected directly into the tumor, due to DNA damage in the immune cells. The combination treatment induced complete tumor regression in a subset of mice that were then able to reject tumor rechallenge, indicating that the injured cell adjuvant treatment induced durable antitumor immunological memory. These results provide a strategy for enhancing the efficacy of immune checkpoint inhibition in tumor types that do not respond to this treatment modality by itself.

Transite: A Computational Motif-Based Analysis Platform That Identifies RNA-Binding Proteins Modulating Changes in Gene Expression

RNA-binding proteins (RBPs) play critical roles in regulating gene expression by modulating splicing, RNA stability, and protein translation. Stimulus-induced alterations in RBP function contribute to global changes in gene expression, but identifying which RBPs are responsible for the observed changes remains an unmet need. Here, we present Transite, a computational approach that systematically infers RBPs influencing gene expression through changes in RNA stability and degradation. As a proof of principle, we apply Transite to RNA expression data from human patients with non-small-cell lung cancer whose tumors were sampled at diagnosis or after recurrence following treatment with platinum-based chemotherapy. Transite implicates known RBP regulators of the DNA damage response and identifies hnRNPC as a new modulator of chemotherapeutic resistance, which we subsequently validated experimentally. Transite serves as a framework for the identification of RBPs that drive cell-state transitions and adds additional value to the vast collection of publicly available gene expression datasets.

Krismer et al. present a computational approach to identify RNA-binding proteins (RBPs) that modulate post-transcriptional control of gene expression using RNA expression data as inputs. By applying this approach to publicly available patient datasets, they identify and experimentally confirm that the RBP hnRNPC contributes to chemotherapy resistance in lung cancer.

Phosphoproteomics identifies microglial Siglec-F inflammatory response during neurodegeneration

Alzheimer's disease (AD) is characterized by the appearance of amyloid-β plaques, neurofibrillary tangles, and inflammation in brain regions involved in memory. Using mass spectrometry, we have quantified the phosphoproteome of the CK-p25, 5XFAD, and Tau P301S mouse models of neurodegeneration. We identified a shared response involving Siglec-F which was upregulated on a subset of reactive microglia. The human paralog Siglec-8 was also upregulated on microglia in AD. Siglec-F and Siglec-8 were upregulated following microglial activation with interferon gamma (IFNγ) in BV-2 cell line and human stem cell-derived microglia models. Siglec-F overexpression activates an endocytic and pyroptotic inflammatory response in BV-2 cells, dependent on its sialic acid substrates and immunoreceptor tyrosine-based inhibition motif (ITIM) phosphorylation sites. Related human Siglecs induced a similar response in BV-2 cells. Collectively, our results point to an important role for mouse Siglec-F and human Siglec-8 in regulating microglial activation during neurodegeneration.

Multiplexed Plasma Immune Mediator Signatures Can Differentiate Sepsis From NonInfective SIRS: American Surgical Association 2020 Annual Meeting Paper

OBJECTIVES

Sepsis and sterile both release "danger signals' that induce the systemic inflammatory response syndrome (SIRS). So differentiating infection from SIRS can be challenging. Precision diagnostic assays could limit unnecessary antibiotic use, improving outcomes.

METHODS

After surveying human leukocyte cytokine production responses to sterile damage-associated molecular patterns (DAMPs), bacterial pathogen-associated molecular patterns, and bacteria we created a multiplex assay for 31 cytokines. We then studied plasma from patients with bacteremia, septic shock, "severe sepsis," or trauma (ISS ≥15 with circulating DAMPs) as well as controls. Infections were adjudicated based on post-hospitalization review. Plasma was studied in infection and injury using univariate and multivariate means to determine how such multiplex assays could best distinguish infective from noninfective SIRS.

RESULTS

Infected patients had high plasma interleukin (IL)-6, IL-1α, and triggering receptor expressed on myeloid cells-1 (TREM-1) compared to controls [false discovery rates (FDR) <0.01, <0.01, <0.0001]. Conversely, injury suppressed many mediators including MDC (FDR <0.0001), TREM-1 (FDR <0.001), IP-10 (FDR <0.01), MCP-3 (FDR <0.05), FLT3L (FDR <0.05), Tweak, (FDR <0.05), GRO-α (FDR <0.05), and ENA-78 (FDR <0.05). In univariate studies, analyte overlap between clinical groups prevented clinical relevance. Multivariate models discriminated injury and infection much better, with the 2-group random-forest model classifying 11/11 injury and 28/29 infection patients correctly in out-of-bag validation.

CONCLUSIONS

Circulating cytokines in traumatic SIRS differ markedly from those in health or sepsis. Variability limits the accuracy of single-mediator assays but machine learning based on multiplexed plasma assays revealed distinct patterns in sepsis- and injury-related SIRS. Defining biomarker release patterns that distinguish specific SIRS populations might allow decreased antibiotic use in those clinical situations. Large prospective studies are needed to validate and operationalize this approach.

Shank3 mutation in a mouse model of autism leads to changes in the S-nitroso-proteome and affects key proteins involved in vesicle release and synaptic function

Mutation in the SHANK3 human gene leads to different neuropsychiatric diseases including Autism Spectrum Disorder (ASD), intellectual disabilities and Phelan-McDermid syndrome. Shank3 disruption in mice leads to dysfunction of synaptic transmission, behavior, and development. Protein S-nitrosylation, the nitric oxide (NO^•)-mediated posttranslational modification (PTM) of cysteine thiols (SNO), modulates the activity of proteins that regulate key signaling pathways. We tested the hypothesis that Shank3 mutation would generate downstream effects on PTM of critical proteins that lead to modification of synaptic functions. SNO-proteins in two ASD-related brain regions, cortex and striatum of young and adult InsG3680(+/+) mice (a human mutation-based Shank3 mouse model), were identified by an innovative mass spectrometric method, SNOTRAP. We found changes of the SNO-proteome in the mutant compared to WT in both ages. Pathway analysis showed enrichment of processes affected in ASD. SNO-Calcineurin in mutant led to a significant increase of phosphorylated Synapsin1 and CREB, which affect synaptic vesicle mobilization and gene transcription, respectively. A significant increase of 3-nitrotyrosine was found in the cortical regions of the adult mutant, signaling both oxidative and nitrosative stress. Neuronal NO^• Synthase (nNOS) was examined for levels and localization in neurons and no significant difference was found in WT vs. mutant. S-nitrosoglutathione concentrations were higher in mutant mice compared to WT. This is the first study on NO^•-related molecular changes and SNO-signaling in the brain of an ASD mouse model that allows the characterization and identification of key proteins, cellular pathways, and neurobiological mechanisms that might be affected in ASD.

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

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.

Substrate-based kinase activity inference identifies MK2 as driver of colitis

Inflammatory bowel disease (IBD) is a chronic and debilitating disorder that has few treatment options due to a lack of comprehensive understanding of its molecular pathogenesis. We used multiplexed mass spectrometry to collect high-content information on protein phosphorylation in two different mouse models of IBD. Because the biological function of the vast majority of phosphorylation sites remains unknown, we developed Substrate-based Kinase Activity Inference (SKAI), a methodology to infer kinase activity from phosphoproteomic data. This approach draws upon prior knowledge of kinase-substrate interactions to construct custom lists of kinases and their respective substrate sites, termed kinase-substrate sets that employ prior knowledge across organisms. This expansion as much as triples the amount of prior knowledge available. We then used these sets within the Gene Set Enrichment Analysis framework to infer kinase activity based on increased or decreased phosphorylation of its substrates in a dataset. When applied to the phosphoproteomic datasets from the two mouse models, SKAI predicted largely non-overlapping kinase activation profiles. These results suggest that chronic inflammation may arise through activation of largely divergent signaling networks. However, the one kinase inferred to be activated in both mouse models was mitogen-activated protein kinase-activated protein kinase 2 (MAPKAPK2 or MK2), a serine/threonine kinase that functions downstream of p38 stress-activated mitogen-activated protein kinase. Treatment of mice with active colitis with ATI450, an orally bioavailable small molecule inhibitor of the MK2 pathway, reduced inflammatory signaling in the colon and alleviated the clinical and histological features of inflammation. These studies establish MK2 as a therapeutic target in IBD and identify ATI450 as a potential therapy for the disease.

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

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Proteogenomic Network Analysis of Context-Specific KRAS Signaling in Mouse-to-Human Cross-Species Translation

The highest frequencies of KRAS mutations occur in colorectal carcinoma (CRC) and pancreatic ductal adenocarcinoma (PDAC). The ability to target downstream pathways mediating KRAS oncogenicity is limited by an incomplete understanding of the contextual cues modulating the signaling output of activated K-RAS. We performed mass spectrometry on mouse tissues expressing wild-type or mutant Kras to determine how tissue context and genetic background modulate oncogenic signaling. Mutant Kras dramatically altered the proteomes and phosphoproteomes of preneoplastic and neoplastic colons and pancreases in a context-specific manner. We developed an approach to statistically humanize the mouse networks with data from human cancer and identified genes within the humanized CRC and PDAC networks synthetically lethal with mutant KRAS. Our studies demonstrate the context-dependent plasticity of oncogenic signaling, identify non-canonical mediators of KRAS oncogenicity within the KRAS-regulated signaling network, and demonstrate how statistical integration of mouse and human datasets can reveal cross-species therapeutic insights.

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.

Tissue-Specific Oncogenic Activity of KRASA146T

KRAS is the most frequently mutated oncogene. The incidence of specific KRAS alleles varies between cancers from different sites, but it is unclear whether allelic selection results from biological selection for specific mutant KRAS proteins. We used a cross-disciplinary approach to compare KRASG12D, a common mutant form, and KRASA146T, a mutant that occurs only in selected cancers. Biochemical and structural studies demonstrated that KRASA146T exhibits a marked extension of switch 1 away from the protein body and nucleotide binding site, which activates KRAS by promoting a high rate of intrinsic and guanine nucleotide exchange factor-induced nucleotide exchange. Using mice genetically engineered to express either allele, we found that KRASG12D and KRASA146T exhibit distinct tissue-specific effects on homeostasis that mirror mutational frequencies in human cancers. These tissue-specific phenotypes result from allele-specific signaling properties, demonstrating that context-dependent variations in signaling downstream of different KRAS mutants drive the KRAS mutational pattern seen in cancer. SIGNIFICANCE: Although epidemiologic and clinical studies have suggested allele-specific behaviors for KRAS, experimental evidence for allele-specific biological properties is limited. We combined structural biology, mass spectrometry, and mouse modeling to demonstrate that the selection for specific KRAS mutants in human cancers from different tissues is due to their distinct signaling properties.See related commentary by Hobbs and Der, p. 696.This article is highlighted in the In This Issue feature, p. 681.

Comprehensive substrate specificity profiling of the human Nek kinome reveals unexpected signaling outputs

Human NimA-related kinases (Neks) have multiple mitotic and non-mitotic functions, but few substrates are known. We systematically determined the phosphorylation-site motifs for the entire Nek kinase family, except for Nek11. While all Nek kinases strongly select for hydrophobic residues in the −3 position, the family separates into four distinct groups based on specificity for a serine versus threonine phospho-acceptor, and preference for basic or acidic residues in other positions. Unlike Nek1-Nek9, Nek10 is a dual-specificity kinase that efficiently phosphorylates itself and peptide substrates on serine and tyrosine, and its activity is enhanced by tyrosine auto-phosphorylation. Nek10 dual-specificity depends on residues in the HRD+2 and APE-4 positions that are uncommon in either serine/threonine or tyrosine kinases. Finally, we show that the phosphorylation-site motifs for the mitotic kinases Nek6, Nek7 and Nek9 are essentially identical to that of their upstream activator Plk1, suggesting that Nek6/7/9 function as phospho-motif amplifiers of Plk1 signaling.

Acidification of Tumor at Stromal Boundaries Drives Transcriptome Alterations Associated with Aggressive Phenotypes

Acidosis is a fundamental feature of the tumor microenvironment, which directly regulates tumor cell invasion by affecting immune cell function, clonal cell evolution, and drug resistance. Despite the important association of tumor microenvironment acidosis with tumor cell invasion, relatively little is known regarding which areas within a tumor are acidic and how acidosis influences gene expression to promote invasion. Here, we injected a labeled pH-responsive peptide to mark acidic regions within tumors. Surprisingly, acidic regions were not restricted to hypoxic areas and overlapped with highly proliferative, invasive regions at the tumor-stroma interface, which were marked by increased expression of matrix metalloproteinases and degradation of the basement membrane. RNA-seq analysis of cells exposed to low pH conditions revealed a general rewiring of the transcriptome that involved RNA splicing and enriched for targets of RNA binding proteins with specificity for AU-rich motifs. Alternative splicing of Mena and CD44, which play important isoform-specific roles in metastasis and drug resistance, respectively, was sensitive to histone acetylation status. Strikingly, this program of alternative splicing was reversed in vitro and in vivo through neutralization experiments that mitigated acidic conditions. These findings highlight a previously underappreciated role for localized acidification of tumor microenvironment in the expression of an alternative splicing-dependent tumor invasion program. SIGNIFICANCE: This study expands our understanding of acidosis within the tumor microenvironment and indicates that acidosis induces potentially therapeutically actionable changes to alternative splicing.

ROS and Oxidative Stress Are Elevated in Mitosis during Asynchronous Cell Cycle Progression and Are Exacerbated by Mitotic Arrest

Although elevated levels of reactive oxygen species (ROS) have been observed in cancer cells and cancer cells aberrantly proliferate, it is not known whether the level of reactive oxygen species and the accumulation of oxidative damage to macromolecules vary across the cell cycle. Here, we measure the prevalence of reactive oxygen species and of biomolecule oxidation across the cell cycle in freely cycling cancer cells. We report that reactive oxygen species vary during the cell cycle and peak in mitosis, resulting in mitotic accumulation of oxidized protein cysteine residues. Prolonged mitotic arrest further increased the levels of ROS and the abundance of oxidatively damaged biomolecules, including cysteine-sulfenic-acid-containing proteins and 8-oxoguanine. These finding suggest that mitotic arrest agents may enhance the effects of ROS-dependent anticancer therapies.

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.

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.

Systems Approaches to Cancer Biology

Cancer systems biology aims to understand cancer as an integrated system of genes, proteins, networks, and interactions rather than an entity of isolated molecular and cellular components. The inaugural Systems Approaches to Cancer Biology Conference, cosponsored by the Association of Early Career Cancer Systems Biologists and the National Cancer Institute of the NIH, focused on the interdisciplinary field of cancer systems biology and the challenging cancer questions that are best addressed through the combination of experimental and computational analyses. Attendees found that elucidating the many molecular features of cancer inevitably reveals new forms of complexity and concluded that ensuring the reproducibility and impact of cancer systems biology studies will require widespread method and data sharing and, ultimately, the translation of important findings to the clinic. Cancer Res; 76(23); 6774-7. ©2016 AACR.

DNA Repair Capacity in Multiple Pathways Predicts Chemoresistance in Glioblastoma Multiforme

Cancer cells can resist the effects of DNA-damaging therapeutic agents via utilization of DNA repair pathways, suggesting that DNA repair capacity (DRC) measurements in cancer cells could be used to identify patients most likely to respond to treatment. However, the limitations of available technologies have so far precluded adoption of this approach in the clinic. We recently developed fluorescence-based multiplexed host cell reactivation (FM-HCR) assays to measure DRC in multiple pathways. Here we apply a mathematical model that uses DRC in multiple pathways to predict cellular resistance to killing by DNA-damaging agents. This model, developed using FM-HCR and drug sensitivity measurements in 24 human lymphoblastoid cell lines, was applied to a panel of 12 patient-derived xenograft (PDX) models of glioblastoma to predict glioblastoma response to treatment with the chemotherapeutic DNA-damaging agent temozolomide. This work showed that, in addition to changes in O⁶-methylguanine DNA methyltransferase (MGMT) activity, small changes in mismatch repair (MMR), nucleotide excision repair (NER), and homologous recombination (HR) capacity contributed to acquired temozolomide resistance in PDX models and led to reduced relative survival prolongation following temozolomide treatment of orthotopic mouse models in vivo Our data indicate that measuring the combined status of MMR, HR, NER, and MGMT provided a more robust prediction of temozolomide resistance than assessments of MGMT activity alone. Cancer Res; 77(1); 198-206. ©2016 AACR.

PTP1B-dependent regulation of receptor tyrosine kinase signaling by the actin-binding protein Mena

The actin-binding protein Mena regulates RTK signaling after growth factor stimulation in tumor cells by a novel mechanism. The alternatively spliced MenaINV isoform disrupts this attenuation to drive sensitivity to growth factors, resistance to targeted inhibitors, and ultimately tumor invasion and metastasis.

During breast cancer progression, alternative mRNA splicing produces functionally distinct isoforms of Mena, an actin regulator with roles in cell migration and metastasis. Aggressive tumor cell subpopulations express Mena^INV, which promotes tumor cell invasion by potentiating EGF responses. However, the mechanism by which this occurs is unknown. Here we report that Mena associates constitutively with the tyrosine phosphatase PTP1B and mediates a novel negative feedback mechanism that attenuates receptor tyrosine kinase signaling. On EGF stimulation, complexes containing Mena and PTP1B are recruited to the EGFR, causing receptor dephosphorylation and leading to decreased motility responses. Mena also interacts with the 5′ inositol phosphatase SHIP2, which is important for the recruitment of the Mena-PTP1B complex to the EGFR. When Mena^INV is expressed, PTP1B recruitment to the EGFR is impaired, providing a mechanism for growth factor sensitization to EGF, as well as HGF and IGF, and increased resistance to EGFR and Met inhibitors in signaling and motility assays. In sum, we demonstrate that Mena plays an important role in regulating growth factor–induced signaling. Disruption of this attenuation by Mena^INV sensitizes tumor cells to low–growth factor concentrations, thereby increasing the migration and invasion responses that contribute to aggressive, malignant cell phenotypes.

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 α-actinin-4 (ACTN4), Tyr(4) and Tyr(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(4) site. We found that introducing a phosphomimetic mutation at only Tyr(31) was sufficient to inhibit the interaction with actin in vitro. However, molecular dynamics simulations predicted that Tyr(31) is mostly buried and that phosphorylation of Tyr(4) would increase the solvent exposure and thus kinase accessibility of Tyr(31). In fibroblast cells, EGF stimulation increased tyrosine phosphorylation of a mutant form of ACTN4 with a phosphorylation-mimicking residue at Tyr(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(4) and Tyr(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(31) can be dictated by the kinase affinity for Tyr(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.

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

Meacham et al. performed a genome-scale shRNA screen for modulators of B-cell leukemia progression in vivo and revealed dramatic distinctions between the relative effects of shRNAs on the growth of tumor cells in culture versus in their native microenvironment. They identified “context-specific” regulators of leukemia development, including the gene encoding the zinc finger protein Phf6. While inactivating mutations in Phf6 are commonly observed in human myeloid and T-cell malignancies, they found that Phf6 suppression in B-cell malignancies impairs tumor progression.

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.

Protein kinases display minimal interpositional dependence on substrate sequence: potential implications for the evolution of signalling networks

Characterization of in vitro substrates of protein kinases by peptide library screening provides a wealth of information on the substrate specificity of kinases for amino acids at particular positions relative to the site of phosphorylation, but provides no information concerning interdependence among positions. High-throughput techniques have recently made it feasible to identify large numbers of in vivo kinase substrates. We used data from experiments on the kinases ATM/ATR and CDK1, and curated CK2 substrates to evaluate the prevalence of interactions between substrate positions within a motif and the utility of these interactions in predicting kinase substrates. Among these data, evidence of interpositional sequence dependencies is strikingly rare, and what dependency exists does little to aid in the prediction of novel kinase substrates. Significant increases in the ability of models to predict kinase-substrate specificity beyond position-independent models must come largely from inclusion of elements of biological and cellular context, rather than further analysis of substrate sequences alone. Our results suggest that, evolutionarily, kinase substrate fitness exists in a smooth energetic landscape. Taken with results from others indicating that phosphopeptide-binding domains do exhibit interpositional dependence, our data suggest that incorporation of new substrate molecules into phospho-signalling networks may be rate-limited by the evolution of suitability for binding by phosphopeptide-binding domains.

Quantitative phospho-proteomics to investigate the polo-like kinase 1-dependent phospho-proteome

Polo-like kinase 1 (PLK1) is a key regulator of mitotic progression and cell division, and small molecule inhibitors of PLK1 are undergoing clinical trials to evaluate their utility in cancer therapy. Despite this importance, current knowledge about the identity of PLK1 substrates is limited. Here we present the results of a proteome-wide analysis of PLK1-regulated phosphorylation sites in mitotic human cells. We compared phosphorylation sites in HeLa cells that were or were not treated with the PLK1-inhibitor BI 4834, by labeling peptides via methyl esterification, fractionation of peptides by strong cation exchange chromatography, and phosphopeptide enrichment via immobilized metal affinity chromatography. Analysis by quantitative mass spectrometry identified 4070 unique mitotic phosphorylation sites on 2069 proteins. Of these, 401 proteins contained one or multiple phosphorylation sites whose abundance was decreased by PLK1 inhibition. These include proteins implicated in PLK1-regulated processes such as DNA damage, mitotic spindle formation, spindle assembly checkpoint signaling, and chromosome segregation, but also numerous proteins that were not suspected to be regulated by PLK1. Analysis of amino acid sequence motifs among phosphorylation sites down-regulated under PLK1 inhibition in this data set identified two potential novel variants of the PLK1 consensus motif.

Spatial exclusivity combined with positive and negative selection of phosphorylation motifs is the basis for context-dependent mitotic signaling

The timing and localization of events during mitosis are controlled by the regulated phosphorylation of proteins by the mitotic kinases, which include Aurora A, Aurora B, Nek2 (never in mitosis kinase 2), Plk1 (Polo-like kinase 1), and the cyclin-dependent kinase complex Cdk1/cyclin B. Although mitotic kinases can have overlapping subcellular localizations, each kinase appears to phosphorylate its substrates on distinct sites. To gain insight into the relative importance of local sequence context in kinase selectivity, identify previously unknown substrates of these five mitotic kinases, and explore potential mechanisms for substrate discrimination, we determined the optimal substrate motifs of these major mitotic kinases by positional scanning oriented peptide library screening (PS-OPLS). We verified individual motifs with in vitro peptide kinetic studies and used structural modeling to rationalize the kinase-specific selection of key motif-determining residues at the molecular level. Cross comparisons among the phosphorylation site selectivity motifs of these kinases revealed an evolutionarily conserved mutual exclusion mechanism in which the positively and negatively selected portions of the phosphorylation motifs of mitotic kinases, together with their subcellular localizations, result in proper substrate targeting in a coordinated manner during mitosis.

PTMScout, a Web resource for analysis of high throughput post-translational proteomics studies

The rate of discovery of post-translational modification (PTM) sites is increasing rapidly and is significantly outpacing our biological understanding of the function and regulation of those modifications. To help meet this challenge, we have created PTMScout, a web-based interface for viewing, manipulating, and analyzing high throughput experimental measurements of PTMs in an effort to facilitate biological understanding of protein modifications in signaling networks. PTMScout is constructed around a custom database of PTM experiments and contains information from external protein and post-translational resources, including gene ontology annotations, Pfam domains, and Scansite predictions of kinase and phosphopeptide binding domain interactions. PTMScout functionality comprises data set comparison tools, data set summary views, and tools for protein assignments of peptides identified by mass spectrometry. Analysis tools in PTMScout focus on informed subset selection via common criteria and on automated hypothesis generation through subset labeling derived from identification of statistically significant enrichment of other annotations in the experiment. Subset selection can be applied through the PTMScout flexible query interface available for quantitative data measurements and data annotations as well as an interface for importing data set groupings by external means, such as unsupervised learning. We exemplify the various functions of PTMScout in application to data sets that contain relative quantitative measurements as well as data sets lacking quantitative measurements, producing a set of interesting biological hypotheses. PTMScout is designed to be a widely accessible tool, enabling generation of multiple types of biological hypotheses from high throughput PTM experiments and advancing functional assignment of novel PTM sites. PTMScout is available at http://ptmscout.mit.edu.

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.

Action-at-a-distance interactions enhance protein binding affinity

The identification of protein mutations that enhance binding affinity may be achieved by computational or experimental means, or by a combination of the two. Sources of affinity enhancement may include improvements to the net balance of binding interactions of residues forming intermolecular contacts at the binding interface, such as packing and hydrogen-bonding interactions. Here we identify noncontacting residues that make substantial contributions to binding affinity and that also provide opportunities for mutations that increase binding affinity of the TEM1 beta-lactamase (TEM1) to the beta-lactamase inhibitor protein (BLIP). A region of BLIP not on the direct TEM1-binding surface was identified for which changes in net charge result in particularly large increases in computed binding affinity. Some mutations to the region have previously been characterized, and our results are in good correspondence with this results of that study. In addition, we propose novel mutations to BLIP that were computed to improve binding significantly without contacting TEM1 directly. This class of noncontacting electrostatic interactions could have general utility in the design and tuning of binding interactions.

A computational method for the analysis and prediction of protein:phosphopeptide-binding sites

Phosphopeptide-binding domains, including the FHA, SH2, WW, WD40, MH2, and Polo-box domains, as well as the 14-3-3 proteins, exert control functions in important processes such as cell growth, division, differentiation, and apoptosis. Structures and mechanisms of phosphopeptide binding are generally diverse, revealing few general principles. A computational method for analysis of phosphopeptide-binding domains was therefore developed to elucidate the physical and chemical nature of phosphopeptide binding, given this lack of structural similarity. The surfaces of nine phosphopeptide-binding proteins, representing seven distinct classes of phosphopeptide-binding modules, were discretized, and encoded with information about amino acid identity, surface curvature, and electrostatic potential at every point on the surface in order to identify local surface properties enriched in phosphoresidue contact sites. Cross-validation indicated that propensities corresponding to this enrichment calculated from a subset of the training data could be used to predict the phosphoresidue contact site on proteins not used in training with no false negative results, and with few unconfirmed positive predictions. The locations of phosphoresidue contact sites were then predicted on the surfaces of the checkpoint kinase Chk1 and the BRCA1 BRCT repeat domain, and these predictions are consistent with recent experimental evidence.

Proteomic identification of 14-3-3zeta as a mitogen-activated protein kinase-activated protein kinase 2 substrate: role in dimer formation and ligand binding

Mitogen-activated protein kinase (MAPK)-activated protein kinase 2 (MAPKAPK2) mediates multiple p38 MAPK-dependent inflammatory responses. To define the signal transduction pathways activated by MAPKAPK2, we identified potential MAPKAPK2 substrates by using a functional proteomic approach consisting of in vitro phosphorylation of neutrophil lysate by active recombinant MAPKAPK2, protein separation by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and phosphoprotein identification by peptide mass fingerprinting with matrix-assisted laser desorption ionization mass spectrometry (MALDI-MS) and protein database analysis. One of the eight candidate MAPKAPK2 substrates identified was the adaptor protein, 14-3-3zeta. We confirmed that MAPKAPK2 interacted with and phosphorylated 14-3-3zeta in vitro and in HEK293 cells. The chemoattractant formyl-methionyl-leucyl-phenylalanine (fMLP) stimulated p38-MAPK-dependent phosphorylation of 14-3-3 proteins in human neutrophils. Mutation analysis showed that MAPKAPK2 phosphorylated 14-3-3zeta at Ser-58. Computational modeling and calculation of theoretical binding energies predicted that both phosphorylation at Ser-58 and mutation of Ser-58 to Asp (S58D) compromised the ability of 14-3-3zeta to dimerize. Experimentally, S58D mutation significantly impaired both 14-3-3zeta dimerization and binding to Raf-1. These data suggest that MAPKAPK2-mediated phosphorylation regulates 14-3-3zeta functions, and this MAPKAPK2 activity may represent a novel pathway mediating p38 MAPK-dependent inflammation.