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.

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.

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.

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.

Cellular ERK phospho-form profiles with conserved preference for a switch-like pattern.

ERK is a member of the MAPK pathway with essential functions in cell proliferation, differentiation, and survival. Complete ERK activation by the kinase MEK requires dual phosphorylation at T and Y within the activation motif TEY. We show that exposure of primary mouse hepatocytes to hepatocyte growth factor (HGF) results in phosphorylation at the activation motif, but not of other residues nearby. To determine the relative abundances of unphosphorylated ERK and the three ERK phospho-forms pT, pY, and pTpY, we employed an extended one-source peptide/phosphopeptide standard method in combination with nanoUPLC-MS. This method enabled us to determine the abundances of phospho-forms with a relative variability of ≤5% (SD). We observed a switch-like preference of ERK phospho-form abundances toward the active, doubly phosphorylated and the inactive, unphosphorylated form. Interestingly, ERK phospho-form profiles were similar upon growth factor and cytokine stimulation. A screening of several murine and human cell systems revealed that the balance between TY- and pTpY-ERK is conserved while the abundances of pT- and pY-ERK are more variable within cell types. We show that the phospho-form profiles do not change by blocking MEK activity suggesting that cellular phosphatases determine the ERK phospho-form distribution. This study provides novel quantitative insights into multisite phosphorylation.

Enzyme-cleavable tandem peptides for quantitative studies in MS-based proteomics.

A novel type of peptide standard is introduced that consists of two peptides combined in one synthetic molecule and separated by a proteolytic cleavage site. Upon enzymatic digestion, the two peptides are released in a molar one-to-one ratio. This method enables the generation of exact equimolar mixtures of two peptides of any nature and origin, thereby providing a valuable tool for the investigation of fundamental phenomena in MS. The applicability of the method is exemplified by the analysis of the effect of peptide sequence variations on the relative ionization efficiency in ESI- and MALDI-MS.

Site-specific degree of phosphorylation in proteins measured by liquid chromatography-electrospray mass spectrometry.

This review focuses on quantitative protein phosphorylation analysis based on coverage of both the phosphorylated and nonphosphorylated forms. In this way, site-specific data on the degree of phosphorylation can be measured, generating the most detailed level of phosphorylation status analysis of proteins. To highlight the experimental challenges in this type of quantitative protein phosphorylation analysis, we discuss the typical workflows for mass spectrometry-based proteomics with a focus on the quantitative analysis of peptide/phosphopeptide ratios. We review workflows for measuring site-specific degrees of phosphorylation including the label-free approach, differential stable isotope labeling of analytes, and methods based on the addition of stable isotope labeled peptide/phosphopeptide pairs as internal standards. The discussion also includes the determination of phosphopeptide isoform abundance data for multiply phosphorylated motifs that contain information about the connectivity of phosphorylation events. The review closes with a prospective on the use of intact stable isotope labeled proteins as internal standards and a summarizing discussion of the typical accuracies of the individual methods.

Improving the precision of quantitative bottom-up proteomics based on stable isotope-labeled proteins.

Stable isotope dilution-based quantitative proteomics with intact labeled proteins as internal standards in combination with a bottom-up approach, i.e., with quantification on the peptide level, is an established method. To explore the technical precision of this approach, calmodulin-like protein 3 was prepared in non-labeled (light) and SILAC-type labeled (heavy) form by cell-free synthesis, mixed, digested with trypsin, and analyzed by UPLC-ESI-MS. In total, 16 light/heavy peptide pair ratios were determined. Pair-wise comparison of ratios of 12 peptides selected according to S/N ratios >50 revealed that the majority exhibited ratios, which were different at a high level of statistical significance (p < 0.001). HPLC-MALDI-MS ratio data confirmed this observation, thus excluding the ionization method as a source of the observed ratio differences. Variation of the digestion time from 0.25 to 4 h showed that the light/heavy ratios of most peptides decrease with time, indicating a kinetic isotope effect leading to preferred cleavage of light calmodulin-like protein 3. The subset of peptides with statistically identical ratios resulted in an average ratio with a RSD of 1.0 %. The light/heavy ratio calculated on the basis of these peptides probably provides the most accurate molar protein ratio.

Impaired neurotransmission in ether lipid-deficient nerve terminals.

Isolated defects of ether lipid (EL) biosynthesis in humans cause rhizomelic chondrodysplasia punctata type 2 and type 3, serious peroxisomal disorders. Using a previously described mouse model [Rodemer, C., Thai, T.P., Brugger, B., Kaercher, T., Werner, H., Nave, K.A., Wieland, F., Gorgas, K., and Just, W.W. (2003) Inactivation of ether lipid biosynthesis causes male infertility, defects in eye development and optic nerve hypoplasia in mice. Hum. Mol. Genet., 12, 1881-1895], we investigated the effect of EL deficiency in isolated murine nerve terminals (synaptosomes) on the pre-synaptic release of the neurotransmitters (NTs) glutamate and acetylcholine. Both Ca(2+)-dependent exocytosis and Ca(2+)-independent efflux of the transmitters were affected. EL-deficient synaptosomes respire at a reduced rate and exhibit a lowered adenosin-5'-triphosphate/adenosine diphosphate (ATP/ADP) ratio. Consequently, ATP-driven processes, such as synaptic vesicle cycling and maintenance of Na(+), K(+) and Ca(2+) homeostasis, might be disturbed. Analyzing reactive oxygen species in EL-deficient neural and non-neural tissues revealed that plasmalogens (PLs), the most abundant EL species in mammalian central nervous system, considerably contribute to the generation of the lipid peroxidation product malondialdehyde. Although EL-deficient tissue contains less lipid peroxidation products, fibroblasts lacking ELs are more susceptible to induced oxidative stress. In summary, these results suggest that due to the reduced energy state of EL-deficient tissue, the Ca(2+)-independent efflux of NTs increases while the Ca(2+)-dependent release declines. Furthermore, lack of PLs is mainly compensated for by an increase in the concentration of phosphatidylethanolamine and results in a significantly lowered level of lipid peroxidation products in the brain cortex and cerebellum.