Erythropoietin-driven dynamic proteome adaptations during erythropoiesis prevent iron overload in the developing embryo.

Erythropoietin (Epo) ensures survival and proliferation of colony-forming unit erythroid (CFU-E) progenitor cells and their differentiation to hemoglobin-containing mature erythrocytes. A lack of Epo-induced responses causes embryonic lethality, but mechanisms regulating the dynamic communication of cellular alterations to the organismal level remain unresolved. By time-resolved transcriptomics and proteomics, we show that Epo induces in CFU-E cells a gradual transition from proliferation signature proteins to proteins indicative for differentiation, including heme-synthesis enzymes. In the absence of the Epo receptor (EpoR) in embryos, we observe a lack of hemoglobin in CFU-E cells and massive iron overload of the fetal liver pointing to a miscommunication between liver and placenta. A reduction of iron-sulfur cluster-containing proteins involved in oxidative phosphorylation in these embryos leads to a metabolic shift toward glycolysis. This link connecting erythropoiesis with the regulation of iron homeostasis and metabolic reprogramming suggests that balancing these interactions is crucial for protection from iron intoxication and for survival.

Light contamination in stable isotope-labelled internal peptide standards is frequent and a potential source of false discovery and quantitation error in proteomics.

In mass spectrometry-based proteomics, heavy internal standards are used to validate target peptide detections and to calibrate peptide quantitation. Here, we report light contamination present in heavy labelled synthetic peptides of high isotopic enrichment. Application of such peptides as assay-internal standards potentially compromises the detection and quantitation especially of low abundant cellular peptides. Therefore, it is important to adopt guidelines to prevent false-positive identifications of endogenous light peptides as well as errors in their quantitation from biological samples.

MSPypeline: a python package for streamlined data analysis of mass spectrometry-based proteomics.

Summary: Mass spectrometry-based proteomics is increasingly employed in biology and medicine. To generate reliable information from large datasets and ensure comparability of results, it is crucial to implement and standardize the quality control of the raw data, the data processing steps and the statistical analyses. MSPypeline provides a platform for importing MaxQuant output tables, generating quality control reports, data preprocessing including normalization and performing exploratory analyses by statistical inference plots. These standardized steps assess data quality, provide customizable figures and enable the identification of differentially expressed proteins to reach biologically relevant conclusions. Availability and implementation: The source code is available under the MIT license at https://github.com/siheming/mspypeline with documentation at https://mspypeline.readthedocs.io. Benchmark mass spectrometry data are available on ProteomeXchange (PXD025792). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

A method for the quantitative determination of glycerophospholipid regioisomers by UPLC-ESI-MS/MS.

Diacyl glycerophospholipids (GPs) belong to the most abundant lipid species in living organisms and consist of a glycerol backbone with fatty acyl groups in sn-1 and sn-2 and a polar head group in the sn-3 position. Regioisomeric mixed diacyl GPs have the same fatty acyl composition but differ in their allocation to sn-1 or sn-2 of the glycerol unit. In-depth analysis of regioisomeric mixed diacyl GP species composed of fatty acyl moieties that are similar in length and degree of saturation typically requires either chemical derivatization or sophisticated analytical instrumentation, since these types of regioisomers are not well resolved under standard ultra-performance liquid chromatography (UPLC) conditions. Here, we introduce a simple and fast method for diacyl GP regioisomer analysis employing UPLC tandem mass spectrometry (MS/MS). This GP regioisomer analysis is based both on minor chromatographic retention time shifts and on major differences in relative abundances of the two fatty acyl anion fragments observed in MS/MS. To monitor these differences with optimal precision, MS/MS spectra are recorded continuously over the UPLC elution profile of the lipid species of interest. Quantification of relative abundances of the regioisomers was performed by algorithms that we have developed for this purpose. The method was applied to commercially available mixed diacyl GP standards and to total lipid extracts of Escherichia coli (E. coli) and bovine liver. To validate our results, we determined regioisomeric ratios of phosphatidylcholine (PC) standards using phospholipase A2-specific release of fatty acids from the sn-2 position of the glycerol backbone. Our results show that most analyzed mixed diacyl GPs of biological origin exhibit significantly higher regioisomeric purity than synthetic lipid standards. In summary, this method can be implemented in routine LC-MS/MS-based lipidomics workflows without the necessity for additional chemical additives, derivatizations, or instrumentation.

Resolving the Combinatorial Complexity of Smad Protein Complex Formation and Its Link to Gene Expression.

Upon stimulation of cells with transforming growth factor ß (TGF-ß), Smad proteins form trimeric complexes and activate a broad spectrum of target genes. It remains unresolved which of the possible Smad complexes are formed in cellular contexts and how these contribute to gene expression. By combining quantitative mass spectrometry with a computational selection strategy, we predict and provide experimental evidence for the three most relevant Smad complexes in the mouse hepatoma cell line Hepa1-6. Utilizing dynamic pathway modeling, we specify the contribution of each Smad complex to the expression of representative Smad target genes, and show that these contributions are conserved in human hepatoma cell lines and primary hepatocytes. We predict, based on gene expression data of patient samples, increased amounts of Smad2/3/4 proteins and Smad2 phosphorylation as hallmarks of hepatocellular carcinoma and experimentally verify this prediction. Our findings demonstrate that modeling approaches can disentangle the complexity of transcription factor complex formation and its impact on gene expression.

Model Based Targeting of IL-6-Induced Inflammatory Responses in Cultured Primary Hepatocytes to Improve Application of the JAK Inhibitor Ruxolitinib.

IL-6 is a central mediator of the immediate induction of hepatic acute phase proteins (APP) in the liver during infection and after injury, but increased IL-6 activity has been associated with multiple pathological conditions. In hepatocytes, IL-6 activates JAK1-STAT3 signaling that induces the negative feedback regulator SOCS3 and expression of APPs. While different inhibitors of IL-6-induced JAK1-STAT3-signaling have been developed, understanding their precise impact on signaling dynamics requires a systems biology approach. Here we present a mathematical model of IL-6-induced JAK1-STAT3 signaling that quantitatively links physiological IL-6 concentrations to the dynamics of IL-6-induced signal transduction and expression of target genes in hepatocytes. The mathematical model consists of coupled ordinary differential equations (ODE) and the model parameters were estimated by a maximum likelihood approach, whereas identifiability of the dynamic model parameters was ensured by the Profile Likelihood. Using model simulations coupled with experimental validation we could optimize the long-term impact of the JAK-inhibitor Ruxolitinib, a therapeutic compound that is quickly metabolized. Model-predicted doses and timing of treatments helps to improve the reduction of inflammatory APP gene expression in primary mouse hepatocytes close to levels observed during regenerative conditions. The concept of improved efficacy of the inhibitor through multiple treatments at optimized time intervals was confirmed in primary human hepatocytes. Thus, combining quantitative data generation with mathematical modeling suggests that repetitive treatment with Ruxolitinib is required to effectively target excessive inflammatory responses without exceeding doses recommended by the clinical guidelines.

IL-1ß-induced and p38MAPK-dependent activation of the mitogen-activated protein kinase-activated protein kinase 2 (MK2) in hepatocytes: Signal transduction with robust and concentration-independent signal amplification.

The IL-1ß induced activation of the p38MAPK/MAPK-activated protein kinase 2 (MK2) pathway in hepatocytes is important for control of the acute phase response and regulation of liver regeneration. Many aspects of the regulatory relevance of this pathway have been investigated in immune cells in the context of inflammation. However, very little is known about concentration-dependent activation kinetics and signal propagation in hepatocytes and the role of MK2. We established a mathematical model for IL-1ß-induced activation of the p38MAPK/MK2 pathway in hepatocytes that was calibrated to quantitative data on time- and IL-1ß concentration-dependent phosphorylation of p38MAPK and MK2 in primary mouse hepatocytes. This analysis showed that, in hepatocytes, signal transduction from IL-1ß via p38MAPK to MK2 is characterized by strong signal amplification. Quantification of p38MAPK and MK2 revealed that, in hepatocytes, at maximum, 11.3% of p38MAPK molecules and 36.5% of MK2 molecules are activated in response to IL-1ß. The mathematical model was experimentally validated by employing phosphatase inhibitors and the p38MAPK inhibitor SB203580. Model simulations predicted an IC50 of 1-1.2 µm for SB203580 in hepatocytes. In silico analyses and experimental validation demonstrated that the kinase activity of p38MAPK determines signal amplitude, whereas phosphatase activity affects both signal amplitude and duration. p38MAPK and MK2 concentrations and responsiveness toward IL-1ß were quantitatively compared between hepatocytes and macrophages. In macrophages, the absolute p38MAPK and MK2 concentration was significantly higher. Finally, in line with experimental observations, the mathematical model predicted a significantly higher half-maximal effective concentration for IL-1ß-induced pathway activation in macrophages compared with hepatocytes, underscoring the importance of cell type-specific differences in pathway regulation.

A timeline of stable isotopes and mass spectrometry in the life sciences.

This review retraces the role of stable isotopes and mass spectrometry in the life sciences. The timeline is divided into four segments covering the years 1920-1950, 1950-1980, 1980-2000, and 2000 until today. For each period methodic progress and typical applications are discussed. Application of stable isotopes is driven by improvements of mass spectrometry, chromatography, and related fields in sensitivity, mass accuracy, structural specificity, complex sample handling ability, data output, and data evaluation. We currently experience the vision of omics-type analyses, that is, the comprehensive identification and quantification of a complete compound class within one or a few analytical runs. This development is driven by stable isotopes without competition by radioisotopes. In metabolic studies as classic field of isotopic tracer experiments, stable isotopes and radioisotopes were competing solutions, with stable isotopes as the long-term junior partner. Since the 1990s the number of metabolic studies with radioisotopes decreases, whereas stable isotope studies retain their slow but stable upward tendency. Unique fields of stable isotopes are metabolic tests in newborns, metabolic experiments in healthy controls, newborn screening for inborn errors, quantification of drugs and drug metabolites in doping control, natural isotope fractionation in geology, ecology, food authentication, or doping control, and more recently the field of quantitative omics-type analyses. There, cells or whole organisms are systematically labeled with stable isotopes to study proteomic differences or specific responses to stimuli or genetic manipulation. The duo of stable isotopes and mass spectrometry will probably continue to grow in the life sciences, since it delivers reference-quality quantitative data with molecular specificity, often combined with informative isotope effects. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:58-85, 2017.

Identification of Cell Type-Specific Differences in Erythropoietin Receptor Signaling in Primary Erythroid and Lung Cancer Cells.

Lung cancer, with its most prevalent form non-small-cell lung carcinoma (NSCLC), is one of the leading causes of cancer-related deaths worldwide, and is commonly treated with chemotherapeutic drugs such as cisplatin. Lung cancer patients frequently suffer from chemotherapy-induced anemia, which can be treated with erythropoietin (EPO). However, studies have indicated that EPO not only promotes erythropoiesis in hematopoietic cells, but may also enhance survival of NSCLC cells. Here, we verified that the NSCLC cell line H838 expresses functional erythropoietin receptors (EPOR) and that treatment with EPO reduces cisplatin-induced apoptosis. To pinpoint differences in EPO-induced survival signaling in erythroid progenitor cells (CFU-E, colony forming unit-erythroid) and H838 cells, we combined mathematical modeling with a method for feature selection, the L1 regularization. Utilizing an example model and simulated data, we demonstrated that this approach enables the accurate identification and quantification of cell type-specific parameters. We applied our strategy to quantitative time-resolved data of EPO-induced JAK/STAT signaling generated by quantitative immunoblotting, mass spectrometry and quantitative real-time PCR (qRT-PCR) in CFU-E and H838 cells as well as H838 cells overexpressing human EPOR (H838-HA-hEPOR). The established parsimonious mathematical model was able to simultaneously describe the data sets of CFU-E, H838 and H838-HA-hEPOR cells. Seven cell type-specific parameters were identified that included for example parameters for nuclear translocation of STAT5 and target gene induction. Cell type-specific differences in target gene induction were experimentally validated by qRT-PCR experiments. The systematic identification of pathway differences and sensitivities of EPOR signaling in CFU-E and H838 cells revealed potential targets for intervention to selectively inhibit EPO-induced signaling in the tumor cells but leave the responses in erythroid progenitor cells unaffected. Thus, the proposed modeling strategy can be employed as a general procedure to identify cell type-specific parameters and to recommend treatment strategies for the selective targeting of specific cell types.

Microquantification of phospholipid classes by stable isotope dilution and nanoESI mass spectrometry.

A new method for microquantification of phospholipid classes by nanoelectrospray mass spectrometry and stable isotope dilution is presented. The method covers the sum of phosphatidylcholine and sphingomyelin and in addition selectively quantifies phosphatidylethanolamine, phosphatidylserine, and phosphatidylinositol. A phospholipid class is quantified together with its corresponding lyso-species due to the presence of a common head group. Phospholipids are extracted from tissue lysates, hydrolysed by hydrofluoric acid, and the liberated polar head groups choline, ethanolamine, serine, and inositol are quantified by nanoelectrospray mass spectrometry using deuterium-labeled analogs of the head groups as internal standards. The method is applied to tissue samples of a gastrointestinal tumor and of corresponding non-affected control tissue. In the tumor sample, the abovementioned phospholipids were found at roughly threefold elevated concentrations with a virtually unaltered relative abundance profile.

Context-specific flow through the MEK/ERK module produces cell- and ligand-specific patterns of ERK single and double phosphorylation.

The same pathway, such as the mitogen-activated protein kinase (MAPK) pathway, can produce different cellular responses, depending on stimulus or cell type. We examined the phosphorylation dynamics of the MAPK kinase MEK and its targets extracellular signal-regulated kinase 1 and 2 (ERK1/2) in primary hepatocytes and the transformed keratinocyte cell line HaCaT A5 exposed to either hepatocyte growth factor or interleukin-6. By combining quantitative mass spectrometry with dynamic modeling, we elucidated network structures for the reversible threonine and tyrosine phosphorylation of ERK in both cell types. In addition to differences in the phosphorylation and dephosphorylation reactions, the HaCaT network model required two feedback mechanisms, which, as the experimental data suggested, involved the induction of the dual-specificity phosphatase DUSP6 and the scaffold paxillin. We assayed and modeled the accumulation of the double-phosphorylated and active form of ERK1/2, as well as the dynamics of the changes in the monophosphorylated forms of ERK1/2. Modeling the differences in the dynamics of the changes in the distributions of the phosphorylated forms of ERK1/2 suggested that different amounts of MEK activity triggered context-specific responses, with primary hepatocytes favoring the formation of double-phosphorylated ERK1/2 and HaCaT A5 cells that produce both the threonine-phosphorylated and the double-phosphorylated form. These differences in phosphorylation distributions explained the threshold, sensitivity, and saturation of the ERK response. We extended the findings of differential ERK phosphorylation profiles to five additional cultured cell systems and matched liver tumor and normal tissue, which revealed context-specific patterns of the various forms of phosphorylated ERK.

Ca(2+) -complex stability of GAPAGPLIVPY peptide in gas and aqueous phase, investigated by affinity capillary electrophoresis and molecular dynamics simulations and compared to mass spectrometric results.

Strong, sequence-specific gas-phase bindings between proline-rich peptides and alkaline earth metal ions in nanoESI-MS experiments were reported by Lehmann et al. (Rapid Commun. Mass Spectrom. 2006, 20, 2404-2410), however its relevance for physiological-like aqueous phase is uncertain. Therefore, the complexes should also be studied in aqueous solution and the relevance of the MS method for binding studies be evaluated. A mobility shift ACE method was used for determining the binding between the small peptide GAPAGPLIVPY and various metal ions in aqueous solution. The findings were compared to the MS results and further explained using computational methods. While the MS data showed a strong alkaline earth ion binding, the ACE results showed nonsignificant binding. The proposed vacuum state complex also decomposed during a molecular dynamic simulation in aqueous solution. This study shows that the formed stable peptide-metal ion adducts in the gas phase by ESI-MS does not imply the existence of analogous adducts in the aqueous phase. Comparing peptide-metal ion interaction under the gaseous MS and aqueous ACE conditions showed huge difference in binding behavior.

Immunosuppression and Aberrant T Cell Development in the Absence of N-Myristoylation.

N-myristoylation refers to the attachment of myristic acid to the N-terminal glycine of proteins and substantially affects their intracellular targeting and functions. The thymus represents an organ with a prominent N-myristoylation activity. To elucidate the role of protein N-myristoylation for thymocyte development, we generated mice with a T cell lineage-specific deficiency in N-myristoyl transferase (Nmt)1 and 2. Depletion of Nmt activity in T cells led to a defective transmission of TCR signals, a developmental blockage of thymocytes at the transition from double-negative 3 to 4 stages, and a reduction of all the following stages. We could demonstrate that Lck and myristoylated alanine-rich C kinase substrate, two main myristoylated kinases in T cells, were mislocalized in the absence of Nmt activity. N-myristoylation was also indispensable for early and distal TCR signaling events such as CD3ζ, Zap70, and Erk activation and for release of cytokines such as IFN-γ and IL-2. As a consequence, the initiation and propagation of the TCR signaling cascade was severely impaired. Furthermore, we showed that the absence of myristoylation had an immunosuppressive effect on T cells in vivo after treatment with CpG and stimulation of the TCR with the staphylococcal enterotoxin B superantigen. Therefore, protein myristoylation is indispensable in T cell development and activation and its inhibition might offer a novel strategy to achieve immunosuppression.

Microquantification of inorganic and organic phosphate by negative ion electrospray tandem mass spectrometry.

A new method for inorganic phosphate microquantification is introduced based on negative ion electrospray tandem mass spectrometry and stable isotope dilution by (18)O4-labeled phosphate. Quantification is performed using the non-labeled and (18)O3-labeled [P(18)O3](-) fragment ions at m/z 79 and 85, respectively, formed by dissociation of the [H2PO4](-) ion at m/z 97 and 105, respectively, visible in negative ion electrospray ionization mass spectrometry (ESI-MS) spectra. Tandem mass spectrometry was selected to remove an overlap with the isobaric [HSO4](-) ion at m/z 97 of sulfate and to establish an optimal sensitivity of the quantification assay. It is demonstrated that the assay can also measure the sum of inorganic and phosphoryl phosphate by prior enzymatic hydrolysis of phosphoryl phosphate. The assay works with phosphate concentrations in the micromolar range and, in combination with nano-ESI, is capable to quantitate absolute amounts of phosphate in the low nanogram range from complex samples.

Identification of isoform-specific dynamics in phosphorylation-dependent STAT5 dimerization by quantitative mass spectrometry and mathematical modeling.

STAT5A and STAT5B are important transcription factors that dimerize and transduce activation signals of cytokine receptors directly to the nucleus. A typical cytokine that mediates STAT5 activation is erythropoietin (Epo). Differential functions of STAT5A and STAT5B have been reported. However, the extent to which phosphorylated STAT5A and STAT5B (pSTAT5A, pSTAT5B) form homo- or heterodimers is not understood, nor is how this might influence the signal transmission to the nucleus. To study this, we designed a concept to investigate the isoform-specific dimerization behavior of pSTAT5A and pSTAT5B that comprises isoform-specific immunoprecipitation (IP), measurement of the degree of phosphorylation, and isoform ratio determination between STAT5A and STAT5B. For the main analytical method, we employed quantitative label-free and -based mass spectrometry. For the cellular model system, we used Epo receptor (EpoR)-expressing BaF3 cells (BaF3-EpoR) stimulated with Epo. Three hypotheses of dimer formation between pSTAT5A and pSTAT5B were used to explain the analytical results by a static mathematical model: formation of (i) homodimers only, (ii) heterodimers only, and (iii) random formation of homo- and heterodimers. The best agreement between experimental data and model simulations was found for the last case. Dynamics of cytoplasmic STAT5 dimerization could be explained by distinct nuclear import rates and individual nuclear retention for homo- and heterodimers of phosphorylated STAT5.

Preparation of heteroelement-incorporated and stable isotope-labeled protein standards for quantitative proteomics.

A major obstacle for further development of quantitative proteomics is the lack of accurately quantified protein standards. The following protocol describes innovative methods for the production of stable isotope-labeled protein standards. Their production is achieved by cell-free protein synthesis, which enables simultaneous incorporation of selenomethionine and stable isotope-labeled amino acids. The selenium tag allows sensitive and accurate quantification by inductively coupled plasma mass spectrometry (ICP-MS). The stable isotope label allows internal standardization in mass spectrometry-based proteomics by electrospray ionization-tandem mass spectrometry (ESI-MS/MS). Both label types can be placed within a single protein RISQ standard (recombinant isotope-labeled and selenium quantified) or can be distributed over two types of related RSQ and RIQ standards for the same target protein (recombinant selenium quantified and recombinant isotope-labeled and quantified). The combination of cell-free synthesis as production method with ICP-MS and ESI-MS/MS as detection methods results in protein standards, which are quantified at an outstanding level of accuracy.

One-source peptide/phosphopeptide ratio standards for accurate and site-specific determination of the degree of phosphorylation.

Reversible protein phosphorylation is a key mediator for intracellular signal transduction. Here we describe an innovative method for the production of pairs of peptide standards designed for quantitative mass spectrometry. These standard pairs can be used for site-specific analysis of the degree of phosphorylation of proteins in a bottom-up approach. The method starts from an isotopically labeled phosphopeptide analogue of the analyte phosphopeptide and ends up with a labeled peptide/phosphopeptide ratio standard in which the molar ratio between the phosphorylated and the unphosphorylated form is exactly defined. The signals of the ratio standard are used to standardize the corresponding analyte signals. This compensates for differences in LC recovery or ionization efficiency between the phosphorylated and unphosphorylated forms. The method can also be extended to quantitative analysis of multisite phosphorylation in a single peptide, which is exemplified for the presence of two phosphorylation sites. Peptide/phosphopeptide ratio standards exhibit high ratio accuracy, since ratio adjustment is performed by volumetric operations only.

Phosphorylation of multifunctional galectins by protein kinases CK1, CK2, and PKA.

Phosphorylation is known to have a strong impact on protein functions. We analyzed members of the lectin family of multifunctional galectins as targets of the protein kinases CK1, CK2, and PKA. Galectins are potent growth regulators able to bind both glycan and peptide motifs at intra- and extracellular sites. Performing in vitro kinase assays, galectin phosphorylation was detected by phosphoprotein staining and autoradiography. The insertion of phosphoryl groups varied to a large extent depending on the type of kinase applied and the respective galectin substrate. Sites of phosphorylation observed in the recombinant galectins were determined by a strategic combination of phosphopeptide enrichment and nano-ultra-performance liquid chromatography tandem mass spectrometry (nanoUPLC-MS/MS). By in silico modeling, phosphorylation sites were visualized three-dimensionally. Our results reveal galectin-type-specific Ser-/Thr-dependent phosphorylation beyond the known example of galectin-3. These data are the basis for functional studies and also illustrate the analytical sensitivity of the applied methods for further work on human lectins.

Biliary phosphatidylcholine and lysophosphatidylcholine profiles in sclerosing cholangitis.

AIM: To analyze phospholipid profiles in intrahepatic bile from patients with primary sclerosing cholangitis (PSC) and secondary sclerosing cholangitis (SSC). METHODS: Intrahepatic bile specimens collected via endoscopic retrograde cholangiography from 41 patients were analyzed. Fourteen of these patients were diagnosed with PSC, 10 with SSC, 11 with choledocholithiasis or no identifiable biliary disease, and 6 with cholangiocellular carcinoma (CCC). Bile acid, cholesterol, protein, and bilirubin contents as well as pancreas lipase activity in bile were determined by biochemical methods. Phosphatidylcholine (PC) and lysophosphatidylcholine (LPC) species were quantified using nano-electrospray ionization tandem mass spectrometry. RESULTS: Bile from all the examined patient groups showed a remarkably similar PC and LPC species composition, with only minor statistical differences. Total biliary PC concentrations were highest in controls (8030 ± 1843 µmol/L) and lowest in patients with CCC (1969 ± 981 µmol/L) (P = 0.005, controls vs SSC and CCC, respectively, P < 0.05). LPC contents in bile were overall low (4.2% ± 1.8%). Biliary LPC/PC ratios and ratios of biliary PC to bilirubin, PC to cholesterol, PC to protein, and PC to bile acids showed no intergroup differences. CONCLUSION: PC and LPC profiles being similar in patients with or without sclerosing cholangitis, these phospholipids are likely not of major pathogenetic importance in this disease group.