Quantitative profiling of posttranslational modifications of pathological tau via sarkosyl fractionation and mass spectrometry.

Tau protein aggregation is associated with posttranslational modifications (PTMs) in 75% of all dementia cases. The distribution of tau pathology and the presence of specific tau phosphorylation sites of interest are typically visualized and measured using antibodies. However, previous knowledge of the target epitopes is required. Additionally, antibodies can be used in a semi-quantitative manner but cannot be used to determine the absolute amount of tau or the extent of the modifications at specific sites or domains. Here we present a discovery assay that characterizes the global qualitative and quantitative tau modification landscape of a sample without a priori knowledge. Our workflow uses sarkosyl fractionation to extract the pathological tau species from sample-limited brain specimens, followed by mass spectrometry (MS) to characterize and quantify tau PTMs. The two-step MS-based proteomics approach includes an exploratory tau PTM analysis and a targeted full-length expressed stable isotope-labeled tau assay, which monitors specific unmodified tau peptides using a heavy isotope-labeled internal standard as a reference. This enables the absolute quantification of the respective tau peptides and the total tau amount in the sample, thus providing the modification extent of tau PTMs. This approach provides precise, comprehensive, qualitative and quantitative tau PTM profiling of the sample. It also enables the detailed molecular comparison of tau across multiple experiments, including a comparison between neurodegenerative diseases, stages of the disease, human patient heterogeneity and characterization of animal models. The approach is useful for studying the molecular features of pathological tau in neurodegeneration. The procedure requires 7-8 d and is suitable for users with expertise in targeted and untargeted MS-based protein analysis.

Integrative systems biology characterizes immune-mediated neurodevelopmental changes in murine Zika virus microcephaly.

Characterizing perturbation of molecular pathways in congenital Zika virus (ZIKV) infection is critical for improved therapeutic approaches. Leveraging integrative systems biology, proteomics, and RNA-seq, we analyzed embryonic brain tissues from an immunocompetent, wild-type congenital ZIKV infection mouse model. ZIKV induced a robust immune response accompanied by the downregulation of critical neurodevelopmental gene programs. We identified a negative correlation between ZIKV polyprotein abundance and host cell cycle-inducing proteins. We further captured the downregulation of genes/proteins, many of which are known to be causative for human microcephaly, including Eomesodermin/T-box Brain Protein 2 (EOMES/TBR2) and Neuronal Differentiation 2 (NEUROD2). Disturbances of distinct molecular pathways in neural progenitors and post-mitotic neurons may contribute to complex brain phenotype of congenital ZIKV infection. Overall, this report on protein- and transcript-level dynamics enhances understanding of the ZIKV immunopathological landscape through characterization of fetal immune response in the developing brain.

Urine Proteomics for Noninvasive Monitoring of Biomarkers in Bronchopulmonary Dysplasia.

INTRODUCTION: Current techniques to diagnose and/or monitor critically ill neonates with bronchopulmonary dysplasia (BPD) require invasive sampling of body fluids, which is suboptimal in these frail neonates. We tested our hypothesis that it is feasible to use noninvasively collected urine samples for proteomics from extremely low gestational age newborns (ELGANs) at risk for BPD to confirm previously identified proteins and biomarkers associated with BPD. METHODS: We developed a robust high-throughput urine proteomics methodology that requires only 50 µL of urine. We utilized the methodology with a proof-of-concept study validating proteins previously identified in invasively collected sample types such as blood and/or tracheal aspirates on urine collected within 72 h of birth from ELGANs (gestational age [26 ± 1.2] weeks) who were admitted to a single Neonatal Intensive Care Unit (NICU), half of whom eventually developed BPD (n = 21), while the other half served as controls (n = 21). RESULTS: Our high-throughput urine proteomics approach clearly identified several BPD-associated changes in the urine proteome recapitulating expected blood proteome changes, and several urinary proteins predicted BPD risk. Interestingly, 16 of the identified urinary proteins are known targets of drugs approved by the Food and Drug Administration. CONCLUSION: In addition to validating numerous proteins, previously found in invasively collected blood, tracheal aspirate, and bronchoalveolar lavage, that have been implicated in BPD pathophysiology, urine proteomics also suggested novel potential therapeutic targets. Ease of access to urine could allow for sequential proteomic evaluations for longitudinal monitoring of disease progression and impact of therapeutic intervention in future studies.

multiFLEX-LF: A Computational Approach to Quantify the Modification Stoichiometries in Label-Free Proteomics Data Sets.

In liquid-chromatography-tandem-mass-spectrometry-based proteomics, information about the presence and stoichiometry of protein modifications is not readily available. To overcome this problem, we developed multiFLEX-LF, a computational tool that builds upon FLEXIQuant, which detects modified peptide precursors and quantifies their modification extent by monitoring the differences between observed and expected intensities of the unmodified precursors. multiFLEX-LF relies on robust linear regression to calculate the modification extent of a given precursor relative to a within-study reference. multiFLEX-LF can analyze entire label-free discovery proteomics data sets in a precursor-centric manner without preselecting a protein of interest. To analyze modification dynamics and coregulated modifications, we hierarchically clustered the precursors of all proteins based on their computed relative modification scores. We applied multiFLEX-LF to a data-independent-acquisition-based data set acquired using the anaphase-promoting complex/cyclosome (APC/C) isolated at various time points during mitosis. The clustering of the precursors allows for identifying varying modification dynamics and ordering the modification events. Overall, multiFLEX-LF enables the fast identification of potentially differentially modified peptide precursors and the quantification of their differential modification extent in large data sets using a personal computer. Additionally, multiFLEX-LF can drive the large-scale investigation of the modification dynamics of peptide precursors in time-series and case-control studies. multiFLEX-LF is available at https://gitlab.com/SteenOmicsLab/multiflex-lf.

Phosphorylation-dependent control of Activity-regulated cytoskeleton-associated protein (Arc) protein by TNIK.

Activity-regulated cytoskeleton-associated protein (Arc) is an immediate early gene product that support neuroplastic changes important for cognitive function and memory formation. As a protein with homology to the retroviral Gag protein, a particular characteristic of Arc is its capacity to self-assemble into virus-like capsids that can package mRNAs and transfer those transcripts to other cells. Although a lot has been uncovered about the contributions of Arc to neuron biology and behavior, very little is known about how different functions of Arc are coordinately regulated both temporally and spatially in neurons. The answer to this question we hypothesized must involve the occurrence of different protein post-translational modifications acting to confer specificity. In this study, we used mass spectrometry and sequence prediction strategies to map novel Arc phosphorylation sites. Our approach led us to recognize serine 67 (S67) and threonine 278 (T278) as residues that can be modified by TNIK, which is a kinase abundantly expressed in neurons that shares many functional overlaps with Arc and has, along with its interacting proteins such as the NMDA receptor, and been implicated as a risk factor for psychiatric disorders. Furthermore, characterization of each residue using site-directed mutagenesis to create S67 and T278 mutant variants revealed that TNIK action at those amino acids can strongly influence Arc's subcellular distribution and self-assembly as capsids. Together, our findings reveal an unsuspected connection between Arc and TNIK. Better understanding of the interplay between these two proteins in neuronal cells could lead to new insights about apparition and progression of psychiatric disorders. Cover Image for this issue: https://doi.org/10.1111/jnc.15077.

FLEXIQuant-LF to quantify protein modification extent in label-free proteomics data.

Improvements in LC-MS/MS methods and technology have enabled the identification of thousands of modified peptides in a single experiment. However, protein regulation by post-translational modifications (PTMs) is not binary, making methods to quantify the modification extent crucial to understanding the role of PTMs. Here, we introduce FLEXIQuant-LF, a software tool for large-scale identification of differentially modified peptides and quantification of their modification extent without knowledge of the types of modifications involved. We developed FLEXIQuant-LF using label-free quantification of unmodified peptides and robust linear regression to quantify the modification extent of peptides. As proof of concept, we applied FLEXIQuant-LF to data-independent-acquisition (DIA) data of the anaphase promoting complex/cyclosome (APC/C) during mitosis. The unbiased FLEXIQuant-LF approach to assess the modification extent in quantitative proteomics data provides a better understanding of the function and regulation of PTMs. The software is available at https://github.com/SteenOmicsLab/FLEXIQuantLF.

Tau molecular diversity contributes to clinical heterogeneity in Alzheimer's disease.

Alzheimer's disease (AD) causes unrelenting, progressive cognitive impairments, but its course is heterogeneous, with a broad range of rates of cognitive decline1. The spread of tau aggregates (neurofibrillary tangles) across the cerebral cortex parallels symptom severity2,3. We hypothesized that the kinetics of tau spread may vary if the properties of the propagating tau proteins vary across individuals. We carried out biochemical, biophysical, MS and both cell- and animal-based-bioactivity assays to characterize tau in 32 patients with AD. We found striking patient-to-patient heterogeneity in the hyperphosphorylated species of soluble, oligomeric, seed-competent tau. Tau seeding activity correlates with the aggressiveness of the clinical disease, and some post-translational modification (PTM) sites appear to be associated with both enhanced seeding activity and worse clinical outcomes, whereas others are not. These data suggest that different individuals with 'typical' AD may have distinct biochemical features of tau. These data are consistent with the possibility that individuals with AD, much like people with cancer, may have multiple molecular drivers of an otherwise common phenotype, and emphasize the potential for personalized therapeutic approaches for slowing clinical progression of AD.

A Fast and Quantitative Method for Post-translational Modification and Variant Enabled Mapping of Peptides to Genomes.

Cross-talk between genes, transcripts, and proteins is the key to cellular responses; hence, analysis of molecular levels as distinct entities is slowly being extended to integrative studies to enhance the understanding of molecular dynamics within cells. Current tools for the visualization and integration of proteomics with other omics datasets are inadequate for large-scale studies. Furthermore, they only capture basic sequence identify, discarding post-translational modifications and quantitation. To address these issues, we developed PoGo to map peptides with associated post-translational modifications and quantification to reference genome annotation. In addition, the tool was developed to enable the mapping of peptides identified from customized sequence databases incorporating single amino acid variants. While PoGo is a command line tool, the graphical interface PoGoGUI enables non-bioinformatics researchers to easily map peptides to 25 species supported by Ensembl genome annotation. The generated output borrows file formats from the genomics field and, therefore, visualization is supported in most genome browsers. For large-scale studies, PoGo is supported by TrackHubGenerator to create web-accessible repositories of data mapped to genomes that also enable an easy sharing of proteogenomics data. With little effort, this tool can map millions of peptides to reference genomes within only a few minutes, outperforming other available sequence-identity based tools. This protocol demonstrates the best approaches for proteogenomics mapping through PoGo with publicly available datasets of quantitative and phosphoproteomics, as well as large-scale studies.

A Cost-Effective High-Throughput Plasma and Serum Proteomics Workflow Enables Mapping of the Molecular Impact of Total Pancreatectomy with Islet Autotransplantation.

Blood is an ideal body fluid for the discovery or monitoring of diagnostic and prognostic protein biomarkers. However, discovering robust biomarkers requires the analysis of large numbers of samples to appropriately represent interindividual variability. To address this analytical challenge, we established a high-throughput and cost-effective proteomics workflow for accurate and comprehensive proteomics at an analytical depth applicable for clinical studies. For validation, we processed 1 µL each from 62 plasma samples in 96-well plates and analyzed the product by quantitative data-independent acquisition liquid chromatography/mass spectrometry; the data were queried using feature quantification with Spectronaut. To show the applicability of our workflow to serum, we analyzed a unique set of samples from 48 chronic pancreatitis patients, pre and post total pancreatectomy with islet autotransplantation (TPIAT) surgery. We identified 16 serum proteins with statistically significant abundance alterations, which represent a molecular signature distinct from that of chronic pancreatitis. In summary, we established a cost-efficient high-throughput workflow for comprehensive proteomics using PVDF-membrane-based digestion that is robust, automatable, and applicable to small plasma and serum volumes, e.g., finger stick. Application of this plasma/serum proteomics workflow resulted in the first mapping of the molecular implications of TPIAT on the serum proteome.

Phosphorylation-Induced Motor Shedding Is Required at Mitosis for Proper Distribution and Passive Inheritance of Mitochondria.

While interphase mitochondria associate with microtubules, mitotic mitochondria dissociate from spindle microtubules and localize in the cell periphery. Here, we show that this redistribution is not mediated by mitochondrial active transport or tethering to the cytoskeleton. Instead, kinesin and dynein, which link mitochondria to microtubules, are shed from the mitochondrial surface. Shedding is driven by phosphorylation of mitochondrial and cytoplasmic targets by CDK1 and Aurora A. Forced recruitment of motor proteins to mitotic mitochondria to override this shedding prevents their proper symmetrical distribution and disrupts the balanced inheritance of mitochondria to daughter cells. Moreover, when mitochondria with bound dynein bind to the mitotic spindle, they arrest cell-cycle progression and produce binucleate cells. Thus, our results show that the regulated release of motor proteins from the mitochondrial surface is a critical mitotic event.

MStern Blotting-High Throughput Polyvinylidene Fluoride (PVDF) Membrane-Based Proteomic Sample Preparation for 96-Well Plates.

We describe a 96-well plate compatible membrane-based proteomic sample processing method, which enables the complete processing of 96 samples (or multiples thereof) within a single workday. This method uses a large-pore hydrophobic PVDF membrane that efficiently adsorbs proteins, resulting in fast liquid transfer through the membrane and significantly reduced sample processing times. Low liquid transfer speeds have prevented the useful 96-well plate implementation of FASP as a widely used membrane-based proteomic sample processing method. We validated our approach on whole-cell lysate and urine and cerebrospinal fluid as clinically relevant body fluids. Without compromising peptide and protein identification, our method uses a vacuum manifold and circumvents the need for digest desalting, making our processing method compatible with standard liquid handling robots. In summary, our new method maintains the strengths of FASP and simultaneously overcomes one of the major limitations of FASP without compromising protein identification and quantification.

Quantitative comparison of a human cancer cell surface proteome between interphase and mitosis.

The cell surface is the cellular compartment responsible for communication with the environment. The interior of mammalian cells undergoes dramatic reorganization when cells enter mitosis. These changes are triggered by activation of the CDK1 kinase and have been studied extensively. In contrast, very little is known of the cell surface changes during cell division. We undertook a quantitative proteomic comparison of cell surface-exposed proteins in human cancer cells that were tightly synchronized in mitosis or interphase. Six hundred and twenty-eight surface and surface-associated proteins in HeLa cells were identified; of these, 27 were significantly enriched at the cell surface in mitosis and 37 in interphase. Using imaging techniques, we confirmed the mitosis-selective cell surface localization of protocadherin PCDH7, a member of a family with anti-adhesive roles in embryos. We show that PCDH7 is required for development of full mitotic rounding pressure at the onset of mitosis. Our analysis provided basic information on how cell cycle progression affects the cell surface. It also provides potential pharmacodynamic biomarkers for anti-mitotic cancer chemotherapy.

Quantitative profiling of peptides from RNAs classified as noncoding.

Only a small fraction of the mammalian genome codes for messenger RNAs destined to be translated into proteins, and it is generally assumed that a large portion of transcribed sequences--including introns and several classes of noncoding RNAs (ncRNAs)--do not give rise to peptide products. A systematic examination of translation and physiological regulation of ncRNAs has not been conducted. Here we use computational methods to identify the products of non-canonical translation in mouse neurons by analysing unannotated transcripts in combination with proteomic data. This study supports the existence of non-canonical translation products from both intragenic and extragenic genomic regions, including peptides derived from antisense transcripts and introns. Moreover, the studied novel translation products exhibit temporal regulation similar to that of proteins known to be involved in neuronal activity processes. These observations highlight a potentially large and complex set of biologically regulated translational events from transcripts formerly thought to lack coding potential.

Basal activity of a PARP1-NuA4 complex varies dramatically across cancer cell lines.

Poly(ADP-ribose) polymerases (PARPs) catalyze poly(ADP-ribose) addition onto proteins, an important posttranslational modification involved in transcription, DNA damage repair, and stem cell identity. Previous studies established the activation of PARP1 in response to DNA damage, but little is known about PARP1 regulation outside of DNA repair. We developed an assay for measuring PARP activity in cell lysates and found that the basal activity of PARP1 was highly variable across breast cancer cell lines, independent of DNA damage. Sucrose gradient fractionation demonstrated that PARP1 existed in at least three biochemically distinct states in both high- and low-activity lines. A discovered complex containing the NuA4 chromatin-remodeling complex and PARP1 was responsible for high basal PARP1 activity, and NuA4 subunits were required for this activity. These findings present a pathway for PARP1 activation and a direct link between PARP1 and chromatin remodeling outside of the DNA damage response.

Gas2l3, a novel constriction site-associated protein whose regulation is mediated by the APC/C Cdh1 complex.

Growth arrest-specific 2-like protein 3 (Gas2l3) was recently identified as an Actin/Tubulin cross-linker protein that regulates cytokinesis. Using cell-free systems from both frog eggs and human cells, we show that the Gas2l3 protein is targeted for ubiquitin-mediated proteolysis by the APC/C(Cdh1) complex, but not by the APC/C(Cdc20) complex, and is phosphorylated by Cdk1 in mitosis. Moreover, late in cytokinesis, Gas2l3 is exclusively localized to the constriction sites, which are the narrowest parts of the intercellular bridge connecting the two daughter cells. Overexpression of Gas2l3 specifically interferes with cell abscission, which is the final stage of cell division, when the cutting of the intercellular bridge at the constriction sites occurs. We therefore suggest that Gas2l3 is part of the cellular mechanism that terminates cell division.

A non-parametric cutout index for robust evaluation of identified proteins.

This paper proposes a novel, automated method for evaluating sets of proteins identified using mass spectrometry. The remaining peptide-spectrum match score distributions of protein sets are compared to an empirical absent peptide-spectrum match score distribution, and a Bayesian non-parametric method reminiscent of the Dirichlet process is presented to accurately perform this comparison. Thus, for a given protein set, the process computes the likelihood that the proteins identified are correctly identified. First, the method is used to evaluate protein sets chosen using different protein-level false discovery rate (FDR) thresholds, assigning each protein set a likelihood. The protein set assigned the highest likelihood is used to choose a non-arbitrary protein-level FDR threshold. Because the method can be used to evaluate any protein identification strategy (and is not limited to mere comparisons of different FDR thresholds), we subsequently use the method to compare and evaluate multiple simple methods for merging peptide evidence over replicate experiments. The general statistical approach can be applied to other types of data (e.g. RNA sequencing) and generalizes to multivariate problems.

Plk1 negatively regulates PRC1 to prevent premature midzone formation before cytokinesis.

To achieve mitosis and cytokinesis, microtubules must assemble into distinct structures at different stages of cell division-mitotic spindles to segregate the chromosomes before anaphase and midzones to keep sister genomes apart and guide the cleavage furrow after anaphase. This temporal regulation is believed to involve Cdk1 kinase, which is inactivated in a switch-like way after anaphase. We found that inhibiting Plk1 caused premature assembly of midzones in cells still in metaphase, breaking the temporal regulation of microtubules. The antiparallel microtubule-bundling protein PRC1 plays a key role in organizing the midzone complex. We found that Plk1 negatively regulates PRC1 through phosphorylation of a single site, Thr-602, near the C-terminus of PRC1. We also found that microtubules stimulated Thr-602 phosphorylation by Plk1. This creates a potential negative feedback loop controlling PRC1 activity. It also made the extent of Thr-602 phosphorylation during mitotic arrest dependent on the mechanism of the arresting drug. Unexpectedly, we could not detect a preanaphase regulatory role for Cdk1 sites on PRC1. We suggest that PRC1 is regulated by Plk1, rather than Cdk1 as previously proposed, because its activity must be spatiotemporally regulated both preanaphase and postanaphase, and Cdk1 activity is too binary for this purpose.

Pleiotropic effects of miR-183~96~182 converge to regulate cell survival, proliferation and migration in medulloblastoma.

Medulloblastomas are the most common malignant brain tumors in children. Several large-scale genomic studies have detailed their heterogeneity, defining multiple subtypes with unique molecular profiles and clinical behavior. Increased expression of the miR-183~96~182 cluster of microRNAs has been noted in several subgroups, including the most clinically aggressive subgroup associated with genetic amplification of MYC. To understand the contribution of miR-183~96~182 to the pathogenesis of this aggressive subtype of medulloblastoma, we analyzed global gene expression and proteomic changes that occur upon modulation of miRNAs in this cluster individually and as a group in MYC-amplified medulloblastoma cells. Knockdown of the full miR-183~96~182 cluster results in enrichment of genes associated with apoptosis and dysregulation of the PI3K/AKT/mTOR signaling axis. Conversely, there is a relative enrichment of pathways associated with migration, metastasis and epithelial to mesenchymal transition, as well as pathways associated with dysfunction of DNA repair in cells with preserved miR-183 cluster expression. Immunocytochemistry and FACS analysis confirm induction of apoptosis upon knockdown of the miR-183 cluster. Importantly, cell-based migration and invasion assays verify the positive regulation of cell motility/migration by the miR-183 cluster, which is largely mediated by miR-182. We show that the effects on cell migration induced by the miR-183 cluster are coupled to the PI3K/AKT/mTOR pathway through differential regulation of AKT1 and AKT2 isoforms. Furthermore, we show that rapamycin inhibits cell motility/migration in medulloblastoma cells and phenocopies miR-183 cluster knockdown. Thus, the miR-183 cluster regulates multiple biological programs that converge to support the maintenance and metastatic potential of medulloblastoma.

Measuring phosphorylation-specific changes in response to kinase inhibitors in mammalian cells using quantitative proteomics.

Many cancers have been associated with the deregulation of kinases, and thus, kinases have become a prime target for the development of cancer treatments. This focus on kinases has resulted in the approval of several small-molecule kinase inhibitors for cancer treatments. Further, the use of these inhibitors as tools to study cancer has provided valuable information about biological mechanisms. However, to date, not much is known about the global effects of kinases on the proteome or phosphoproteome. In this protocol, we describe methodology to study the impact of kinase inhibitors on the proteome and phosphoproteome using mass spectrometry-based quantitative proteomics. More specifically, we focus on the effects of Aurora B kinase inhibitors on the proteome, cytoskeleton proteome, the phosphoproteome, and the cytoskeleton phosphoproteome during cell cycle. This methodology is easily extended to other biological studies whose aim is to study the global proteomic effects of a kinase inhibitor.

Interaction of survival of motor neuron (SMN) and HuD proteins with mRNA cpg15 rescues motor neuron axonal deficits.

Spinal muscular atrophy (SMA), caused by the deletion of the SMN1 gene, is the leading genetic cause of infant mortality. SMN protein is present at high levels in both axons and growth cones, and loss of its function disrupts axonal extension and pathfinding. SMN is known to associate with the RNA-binding protein hnRNP-R, and together they are responsible for the transport and/or local translation of ß-actin mRNA in the growth cones of motor neurons. However, the full complement of SMN-interacting proteins in neurons remains unknown. Here we used mass spectrometry to identify HuD as a novel neuronal SMN-interacting partner. HuD is a neuron-specific RNA-binding protein that interacts with mRNAs, including candidate plasticity-related gene 15 (cpg15). We show that SMN and HuD form a complex in spinal motor axons, and that both interact with cpg15 mRNA in neurons. CPG15 is highly expressed in the developing ventral spinal cord and can promote motor axon branching and neuromuscular synapse formation, suggesting a crucial role in the development of motor axons and neuromuscular junctions. Cpg15 mRNA previously has been shown to localize into axonal processes. Here we show that SMN deficiency reduces cpg15 mRNA levels in neurons, and, more importantly, cpg15 overexpression partially rescues the SMN-deficiency phenotype in zebrafish. Our results provide insight into the function of SMN protein in axons and also identify potential targets for the study of mechanisms that lead to the SMA pathology and related neuromuscular diseases.