Quantitative label-free phosphoproteomics strategy for multifaceted experimental designs

Advances in mass spectrometry-based phosphoproteomics for elucidating abscisic acid signaling and plant responses to abiotic stress

Abiotic stresses have significant impacts on crop yield and quality. Even though significant efforts during the past decade have been devoted to uncovering the core signaling pathways associated with the phytohormone abscisic acid (ABA) and abiotic stress in plants, abiotic stress signaling mechanisms in most crops remain largely unclear. The core components of the ABA signaling pathway, including early events in the osmotic stress-induced phosphorylation network, have recently been elucidated in Arabidopsis with the aid of phosphoproteomics technologies. We now know that SNF1-related kinases 2 (SnRK2s) are not only inhibited by the clade A type 2C protein phosphatases (PP2Cs) through dephosphorylation, but also phosphorylated and activated by upstream mitogen-activated protein kinase kinase kinases (MAP3Ks). Through describing the course of studies to elucidate abiotic stress and ABA signaling, we will discuss how we can take advantage of the latest innovations in mass-spectrometry-based phosphoproteomics and structural proteomics to boost our investigation of plant regulation and responses to ABA and abiotic stress.

Gel-Based and Gel-Free Phosphoproteomics to Measure and Characterize Mitochondrial Phosphoproteins

The mitochondrion is a key intracellular organelle regulating metabolic processes, oxidative stress, energy production, calcium homeostasis, and cell survival. Protein phosphorylation plays an important role in regulating mitochondrial functions and cellular signaling pathways. Dysregulation of protein phosphorylation status can cause protein malfunction and abnormal signal transduction, leading to organ dysfunction and disease. Investigating the mitochondrial phosphoproteins is therefore crucial to better understand the molecular and pathogenic mechanisms of many metabolic disorders. Conventional analyses of phosphoproteins, for instance, via western blotting, can be done only for proteins for which specific antibodies to their phosphorylated forms are available. Moreover, such an approach is not suitable for large-scale study of phosphoproteins. Currently, proteomics represents an important tool for large-scale analysis of proteins and their post-translational modifications, including phosphorylation. Here, we provide step-by-step protocols for the proteomics analysis of mitochondrial phosphoproteins (the phosphoproteome), using renal tubular cells as an example. These protocols include methods to effectively isolate mitochondria and to validate the efficacy of mitochondrial enrichment as well as its purity. We also provide detailed protocols for performing both gel-based and gel-free phosphoproteome analyses. The gel-based analysis involves two-dimensional gel electrophoresis and phosphoprotein-specific staining, followed by protein identification via mass spectrometry, whereas the gel-free approach is based on in-solution mass spectrometric identification of specific phosphorylation sites and residues. In all, these approaches allow large-scale analyses of mitochondrial phosphoproteins that can be applied to other cells and tissues of interest. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Mitochondrial isolation/purification from renal tubular cells Support Protocol: Validation of enrichment efficacy and purity of mitochondrial isolation Basic Protocol 2: Gel-based phosphoproteome analysis Basic Protocol 3: Gel-free phosphoproteome analysis.

Juvenile hormone membrane signaling phosphorylates USP and thus potentiates 20-hydroxyecdysone action in Drosophila

Juvenile hormone (JH) and 20-hydroxyecdysone (20E) coordinately regulate development and metamorphosis in insects. Two JH intracellular receptors, methoprene-tolerant (Met) and germ-cell expressed (Gce), have been identified in the fruit fly Drosophila melanogaster. To investigate JH membrane signaling pathway without the interference from JH intracellular signaling, we characterized phosphoproteome profiles of the Met gce double mutant in the absence or presence of JH in both chronic and acute phases. Functioning through a potential receptor tyrosine kinase and phospholipase C pathway, JH membrane signaling activated protein kinase C (PKC) which phosphorylated ultraspiracle (USP) at Ser35, the PKC phosphorylation site required for the maximal action of 20E through its nuclear receptor complex EcR-USP. The usp^S35A mutant, in which Ser was replaced with Ala at position 35 by genome editing, showed decreased expression of Halloween genes that are responsible for ecdysone biosynthesis and thus attenuated 20E signaling that delayed developmental timing. The usp^S35A mutant also showed lower Yorkie activity that reduced body size. Altogether, JH membrane signaling phosphorylates USP at Ser35 and thus potentiates 20E action that regulates the normal fly development. This study helps better understand the complex JH signaling network.

Phosphoproteomic analysis of lung tissue from patients with pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) is a rare disorder associated with high morbidity and mortality despite currently available treatments. We compared the phosphoproteome of lung tissue from subjects with idiopathic PAH (iPAH) obtained at the time of lung transplant with control lung tissue. The mass spectrometry-based analysis found 60,428 phosphopeptide features from which 6622 proteins were identified. Within the subset of identified proteins there were 1234 phosphopeptides with q < 0.05, many of which are involved in immune regulation, angiogenesis, and cell proliferation. Most notably there was a marked relative increase in phosphorylated (S378) IKZF3 (Aiolos), a zinc finger transcription factor that plays a key role in lymphocyte regulation. In vitro phosphorylation assays indicated that GSK3 alpha and/or GSK3 beta could phosphorylate IKZF3 at S378. Western blot analysis demonstrated increased pIKZF3 in iPAH lungs compared to controls. Immunohistochemistry demonstrated phosphorylated IKZF3 in lymphocytes surrounding severely hypertrophied pulmonary arterioles. In situ hybrization showed gene expression in lymphocyte aggregates in PAH samples. A BCL2 reporter assay showed that IKZF3 increased BCL2 promoter activity and demonstrated the potential role of phosphorylation of IKZF3 in the regulation of BCL mediated transcription. Kinase network analysis demonstrated potentially important regulatory roles of casein kinase 2, cyclin-dependent kinase 1 (CDK1), mitogen-associated protein kinases (MAPKs), and protein kinases (PRKs) in iPAH. Bioinformatic analysis demonstrated enrichment of RhoGTPase signaling and the potential importance of cGMP-dependent protein kinase 1 (PRKG). In conclusion, this unbiased phosphoproteomic analysis demonstrated several novel targets regulated by kinase networks in iPAH, and reinforced the potential role of immune regulation in the pathogenesis of iPAH. The identified up- and down-regulated phosphoproteins have potential to serve as biomarkers for PAH and to provide new insights for therapeutic strategies.

Quantitative Phosphoproteomic Comparison of Lens Proteins in Highly Myopic Cataract and Age-Related Cataract

Purpose

To investigate and compare the lens phosphoproteomes in patients with highly myopic cataract (HMC) or age-related cataract (ARC).

Methods

In this study, we undertook a comparative phosphoproteome analysis of the lenses from patients with HMC or ARC. Intact lenses from ARC and HMC patients were separated into the cortex and nucleus. After protein digestion, the phosphopeptides were quantitatively analyzed with TiO₂ enrichment and liquid chromatography-mass spectrometry. The potential functions of different phosphopeptides were assessed by Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis.

Results

In total, 522 phosphorylation sites in 164 phosphoproteins were identified. The number of phosphorylation sites was significantly higher in the cortex than in the nucleus, in both ARC and HMC lenses. The differentially phosphorylated peptides in the lens cortex and nucleus in HMC eyes were significantly involved in the glutathione metabolism pathway. The KEGG pathway enrichment analysis indicated that the differences in phosphosignaling mediators between the ARC and HMC lenses were associated with glycolysis and the level of phosphorylated phosphoglycerate kinase 1 was lower in HMC lenses than in ARC lenses.

Conclusions

We provide an overview of the differential phosphoproteomes of HMC and ARC lenses that can be used to clarify the molecular mechanisms underlying their different phenotypes.

Label-free quantitative proteomic analysis identifies CTNNB1 as a direct target of FOXP3 in gastric cancer cells

Forkhead box protein 3 (FOXP3) is expressed in numerous types of tumor cell and is associated with tumor progression and prognosis. A previous study reported that FOXP3 inhibited cellular proliferation and induced apoptosis of gastric cancer (GC) cells by activating the apoptosis signaling pathway. In the present study, label-free quantitative proteomic analysis and chromatin immunoprecipitation-polymerase chain reaction (ChIP-PCR) was performed to investigate the mechanism by which the anticancer role of FOXP3 was mediated and the proteins that with which it may interact. Label-free quantitative proteomic analysis was used to screen for proteins differentially expressed between FOXP3-overexpressing GC (AF) and vector (ANC) cells. Catenin β1 (CTNNB1) was one of the proteins that exhibited the greatest difference between AF and ANC among 3,313 proteins identified by liquid chromatography with tandem mass spectrometry analysis. The expression of CTNNB1 was evaluated by reverse transcription-quantitative PCR and western blotting. The association between FOXP3 and CTNNB1 was confirmed by ChIP-PCR in AGS cells. The changes in expression of epithelial-mesenchymal transition-associated proteins were analyzed by western blotting. The level of FOXP3 expression was positively associated with CTNNB1 and E-cadherin expression, but not with vimentin and N-cadherin expression. FOXP3 positively regulates CTNNB1 and binds to it directly. Along with the upregulation of glycogen synthase kinase 3β (GSK3β), which was also a protein whose expression was found to change significantly in proteomic analysis and has a key role in the Wnt pathway. This association is an attractive and novel hypothesis for the mechanism by which FOXP3 inhibits the invasion and metastasis of GC cells.

Acylation of Superoxide Dismutase 1 (SOD1) at K122 Governs SOD1-Mediated Inhibition of Mitochondrial Respiration

In this study, we employed proteomics to identify mechanisms of posttranslational regulation on cell survival signaling proteins. We focused on Cu-Zn superoxide dismutase (SOD1), which protects cells from oxidative stress. We found that acylation of K122 on SOD1, while not impacting SOD1 catalytic activity, suppressed the ability of SOD1 to inhibit mitochondrial metabolism at respiratory complex I. We found that deacylase depletion increased K122 acylation on SOD1, which blocked the suppression of respiration in a K122-dependent manner. In addition, we found that acyl-mimicking mutations at K122 decreased SOD1 accumulation in mitochondria, initially hinting that SOD1 may inhibit respiration directly within the intermembrane space (IMS). However, surprisingly, we found that forcing the K122 acyl mutants into the mitochondria with an IMS-targeting tag did not recover their ability to suppress respiration. Moreover, we found that suppressing or boosting respiration levels toggled SOD1 in or out of the mitochondria, respectively. These findings place SOD1-mediated inhibition of respiration upstream of its mitochondrial localization. Lastly, deletion-rescue experiments show that a respiration-defective mutant of SOD1 is also impaired in its ability to rescue cells from toxicity caused by SOD1 deletion. Together, these data suggest a previously unknown interplay between SOD1 acylation, metabolic regulation, and SOD1-mediated cell survival.

Analyses of PDE-regulated phosphoproteomes reveal unique and specific cAMP-signaling modules in T cells

Specific functions for different cyclic nucleotide phosphodiesterases (PDEs) have not yet been identified in most cell types. Conventional approaches to study PDE function typically rely on measurements of global cAMP, general increases in cAMP-dependent protein kinase (PKA), or the activity of exchange protein activated by cAMP (EPAC). Although newer approaches using subcellularly targeted FRET reporter sensors have helped define more compartmentalized regulation of cAMP, PKA, and EPAC, they have limited ability to link this regulation to downstream effector molecules and biological functions. To address this problem, we have begun to use an unbiased mass spectrometry-based approach coupled with treatment using PDE isozyme-selective inhibitors to characterize the phosphoproteomes of the functional pools of cAMP/PKA/EPAC that are regulated by specific cAMP-PDEs (the PDE-regulated phosphoproteomes). In Jurkat cells we find multiple, distinct PDE-regulated phosphoproteomes that can be defined by their responses to different PDE inhibitors. We also find that little phosphorylation occurs unless at least two different PDEs are concurrently inhibited in these cells. Moreover, bioinformatics analyses of these phosphoproteomes provide insight into the unique functional roles, mechanisms of action, and synergistic relationships among the different PDEs that coordinate cAMP-signaling cascades in these cells. The data strongly suggest that the phosphorylation of many different substrates contributes to cAMP-dependent regulation of these cells. The findings further suggest that the approach of using selective, inhibitor-dependent phosphoproteome analysis can provide a generalized methodology for understanding the roles of different PDEs in the regulation of cyclic nucleotide signaling.

An overview of experimental designs in HPLC method development and validation

Chemometric approaches have been increasingly viewed as precious complements to high performance liquid chromatographic practices, since a large number of variables can be simultaneously controlled to achieve the desired separations. Moreover, their applications may efficiently identify and optimize the significant factors to accomplish competent results through limited experimental trials. The present manuscript discusses usefulness of various chemometric approaches in high and ultra performance liquid chromatography for (i) methods development from dissolution studies and sample preparation to detection, considering the progressive substitution of traditional detectors with tandem mass spectrometry instruments and the importance of stability indicating assays (ii) method validation through screening and optimization designs. Choice of appropriate types of experimental designs so as to either screen the most influential factors or optimize the selected factors' combination and the mathematical models in chemometry have been briefly recalled and the advantages of chemometric approaches have been emphasized. The evolution of the design of experiments to the Quality by Design paradigm for method development has been reviewed and the Six Sigma practice as a quality indicator in chromatography has been explained. Chemometric applications and various strategies in chromatographic separations have been described.

The analysis of heterotaxy patients reveals new loss-of-function variants of GRK5

G protein-coupled receptor kinase 5 (GRK5) is a regulator of cardiac performance and a potential therapeutic target in heart failure in the adult. Additionally, we have previously classified GRK5 as a determinant of left-right asymmetry and proper heart development using zebrafish. We thus aimed to identify GRK5 variants of functional significance by analysing 187 individuals with laterality defects (heterotaxy) that were associated with a congenital heart defect (CHD). Using Sanger sequencing we identified two moderately frequent variants in GRK5 with minor allele frequencies <10%, and seven very rare polymorphisms with minor allele frequencies <1%, two of which are novel variants. Given their evolutionarily conserved position in zebrafish, in-depth functional characterisation of four variants (p.Q41L, p.G298S, p.R304C and p.T425M) was performed. We tested the effects of these variants on normal subcellular localisation and the ability to desensitise receptor signalling as well as their ability to correct the left-right asymmetry defect upon Grk5l knockdown in zebrafish. While p.Q41L, p.R304C and p.T425M responded normally in the first two aspects, neither p.Q41L nor p.R304C were capable of rescuing the lateralisation phenotype. The fourth variant, p.G298S was identified as a complete loss-of-function variant in all assays and provides insight into the functions of GRK5.

The G Protein-Coupled Estrogen Receptor Agonist G-1 Inhibits Nuclear Estrogen Receptor Activity and Stimulates Novel Phosphoproteomic Signatures

Estrogen exerts cellular effects through both nuclear (ESR1 and ESR2) and membrane-bound estrogen receptors (G-protein coupled estrogen receptor, GPER); however, it is unclear if they act independently or engage in crosstalk to influence hormonal responses. To investigate each receptor's role in proliferation, transcriptional activation, and protein phosphorylation in breast cancer cells (MCF-7), we employed selective agonists for ESR1 propyl-pyrazole-triol (PPT), ESR2 diarylpropionitrile (DPN), and GPER (G-1) and also determined the impact of xenoestrogens bisphenol-A (BPA) and genistein on these effects. As anticipated, 17β-estradiol (E2), PPT, DPN, BPA, and genistein each enhanced proliferation and activation of an ERE-driven reporter gene whereas G-1 had no significant impact. However, G-1 significantly reduced E2-, PPT-, DPN-, BPA-, and genistein-induced proliferation and ERE activation at doses greater than 500 nM indicating that G-1 mediated inhibition is not ESR isotype specific. As membrane receptors initiate cascades of phosphorylation events, we performed a global phosphoproteomic analysis on cells exposed to E2 or G-1 to identify potential targets of receptor crosstalk via downstream protein phosphorylation targets. Of the 211 phosphorylated proteins identified, 40 and 13 phosphoproteins were specifically modified by E2 and G-1, respectively. Subnetwork enrichment analysis revealed several processes related to cell cycle were specifically enriched by G-1 compared with E2. Further there existed a number of newly identified proteins that were specifically phosphorylated by G-1. These phosphorylation networks highlight specific proteins that may modulate the inhibitory effects of G-1 and suggest a novel role for interference with nuclear receptor activity driven by E2 and xenoestrogens.

S-Nitrosylation of Sarcomeric Proteins Depresses Myofilament Ca2+)Sensitivity in Intact Cardiomyocytes

AIMS

The heart responds to physiological and pathophysiological stress factors by increasing its production of nitric oxide (NO), which reacts with intracellular glutathione to form S-nitrosoglutathione (GSNO), a protein S-nitrosylating agent. Although S-nitrosylation protects some cardiac proteins against oxidative stress, direct effects on myofilament performance are unknown. We hypothesize that S-nitrosylation of sarcomeric proteins will modulate the performance of cardiac myofilaments.

RESULTS

Incubation of intact mouse cardiomyocytes with S-nitrosocysteine (CysNO, a cell-permeable low-molecular-weight nitrosothiol) significantly decreased myofilament Ca(2+) sensitivity. In demembranated (skinned) fibers, S-nitrosylation with 1 μM GSNO also decreased Ca(2+) sensitivity of contraction and 10 μM reduced maximal isometric force, while inhibition of relaxation and myofibrillar ATPase required higher concentrations (≥ 100 μM). Reducing S-nitrosylation with ascorbate partially reversed the effects on Ca(2+) sensitivity and ATPase activity. In live cardiomyocytes treated with CysNO, resin-assisted capture of S-nitrosylated protein thiols was combined with label-free liquid chromatography-tandem mass spectrometry to quantify S-nitrosylation and determine the susceptible cysteine sites on myosin, actin, myosin-binding protein C, troponin C and I, tropomyosin, and titin. The ability of sarcomere proteins to form S-NO from 10-500 μM CysNO in intact cardiomyocytes was further determined by immunoblot, with actin, myosin, myosin-binding protein C, and troponin C being the more susceptible sarcomeric proteins.

INNOVATION AND CONCLUSIONS

Thus, specific physiological effects are associated with S-nitrosylation of a limited number of cysteine residues in sarcomeric proteins, which also offer potential targets for interventions in pathophysiological situations.

Label-free quantitative phosphoproteomics with novel pairwise abundance normalization reveals synergistic RAS and CIP2A signaling

Hyperactivated RAS drives progression of many human malignancies. However, oncogenic activity of RAS is dependent on simultaneous inactivation of protein phosphatase 2A (PP2A) activity. Although PP2A is known to regulate some of the RAS effector pathways, it has not been systematically assessed how these proteins functionally interact. Here we have analyzed phosphoproteomes regulated by either RAS or PP2A, by phosphopeptide enrichment followed by mass-spectrometry-based label-free quantification. To allow data normalization in situations where depletion of RAS or PP2A inhibitor CIP2A causes a large uni-directional change in the phosphopeptide abundance, we developed a novel normalization strategy, named pairwise normalization. This normalization is based on adjusting phosphopeptide abundances measured before and after the enrichment. The superior performance of the pairwise normalization was verified by various independent methods. Additionally, we demonstrate how the selected normalization method influences the downstream analyses and interpretation of pathway activities. Consequently, bioinformatics analysis of RAS and CIP2A regulated phosphoproteomes revealed a significant overlap in their functional pathways. This is most likely biologically meaningful as we observed a synergistic survival effect between CIP2A and RAS expression as well as KRAS activating mutations in TCGA pan-cancer data set, and synergistic relationship between CIP2A and KRAS depletion in colony growth assays.

Quantifying reversible oxidation of protein thiols in photosynthetic organisms

Photosynthetic organisms use dynamic post-translational modifications to survive and adapt, which include reversible oxidative modifications of protein thiols that regulate protein structure, function, and activity. Efforts to quantify thiol modifications on a global scale have relied upon peptide derivatization, typically using isobaric tags such as TMT, ICAT, or iTRAQ that are more expensive, less accurate, and provide less proteome coverage than label-free approaches--suggesting the need for improved experimental designs for studies requiring maximal coverage and precision. Herein, we present the coverage and precision of resin-assisted thiol enrichment coupled to label-free quantitation for the characterization of reversible oxidative modifications on protein thiols. Using C. reinhardtii and Arabidopsis as model systems for algae and plants, we quantified 3662 and 1641 unique cysteinyl peptides, respectively, with median coefficient of variation (CV) of 13% and 16%. Further, our method is extendable for the detection of protein abundance changes and stoichiometries of cysteine oxidation. Finally, we demonstrate proof-of-principle for our method, and reveal that exogenous hydrogen peroxide treatment regulates the C. reinhardtii redox proteome by increasing or decreasing the level of oxidation of 501 or 67 peptides, respectively. As protein activity and function is controlled by oxidative modifications on protein thiols, resin-assisted thiol enrichment coupled to label-free quantitation can reveal how intracellular and environmental stimuli affect plant survival and fitness through oxidative stress.

Targeted quantitative phosphoproteomic analysis of erythrocyte membranes during blood bank storage

One of the hallmarks of blood bank stored red blood cells (RBCs) is the irreversible transition from a discoid to a spherocyte-like morphology with membrane perturbation and cytoskeleton disorders. Therefore, identification of the storage-associated modifications in the protein-protein interactions between the cytoskeleton and the lipid bilayer may contribute to enlighten the molecular mechanisms involved in the alterations of mechanical properties of stored RBCs. Here we report the results obtained analyzing RBCs after 0, 21 and 35 days of storage under standard blood banking conditions by label free mass spectrometry (MS)-based experiments. We could quantitatively measure changes in the phosphorylation level of crucial phosphopeptides belonging to β-spectrin, ankyrin-1, α-adducin, dematin, glycophorin A and glycophorin C proteins. Data have been validated by both western blotting and pseudo-Multiple Reaction Monitoring (MRM). Although each phosphopeptide showed a distinctive trend, a sharp increase in the phosphorylation level during the storage duration was observed. Phosphopeptide mapping and structural modeling analysis indicated that the phosphorylated residues localize in protein functional domains fundamental for the maintenance of membrane structural integrity. Along with previous morphological evidence acquired by electron microscopy, our results seem to indicate that 21-day storage may represent a key point for the molecular processes leading to the erythrocyte deformability reduction observed during blood storage. These findings could therefore be helpful in understanding and preventing the morphology-linked mechanisms responsible for the post-transfusion survival of preserved RBCs.

Phosphoproteomic profiling of human myocardial tissues distinguishes ischemic from non-ischemic end stage heart failure

The molecular differences between ischemic (IF) and non-ischemic (NIF) heart failure are poorly defined. A better understanding of the molecular differences between these two heart failure etiologies may lead to the development of more effective heart failure therapeutics. In this study extensive proteomic and phosphoproteomic profiles of myocardial tissue from patients diagnosed with IF or NIF were assembled and compared. Proteins extracted from left ventricular sections were proteolyzed and phosphopeptides were enriched using titanium dioxide resin. Gel- and label-free nanoscale capillary liquid chromatography coupled to high resolution accuracy mass tandem mass spectrometry allowed for the quantification of 4,436 peptides (corresponding to 450 proteins) and 823 phosphopeptides (corresponding to 400 proteins) from the unenriched and phospho-enriched fractions, respectively. Protein abundance did not distinguish NIF from IF. In contrast, 37 peptides (corresponding to 26 proteins) exhibited a ≥ 2-fold alteration in phosphorylation state (p<0.05) when comparing IF and NIF. The degree of protein phosphorylation at these 37 sites was specifically dependent upon the heart failure etiology examined. Proteins exhibiting phosphorylation alterations were grouped into functional categories: transcriptional activation/RNA processing; cytoskeleton structure/function; molecular chaperones; cell adhesion/signaling; apoptosis; and energetic/metabolism. Phosphoproteomic analysis demonstrated profound post-translational differences in proteins that are involved in multiple cellular processes between different heart failure phenotypes. Understanding the roles these phosphorylation alterations play in the development of NIF and IF has the potential to generate etiology-specific heart failure therapeutics, which could be more effective than current therapeutics in addressing the growing concern of heart failure.

Assessment of BAK1 activity in different plant receptor-like kinase complexes by quantitative profiling of phosphorylation patterns

Plant receptor-like kinases (RLKs) constitute a large family of receptors coordinating developmental programs with adaptation to environmental stresses including immune defenses. BRI1-ASSOCIATED KINASE 1 (BAK1), a member of the plant RLK family, forms receptor complexes with multiple RLK proteins including BRI1, FLS2, EFR and BIK1 to regulate responses to growth hormones or PAMPs. RLK activation and signal initiation involve protein complex formation and phosphorylation/dephosphorylation between BAK1 and its interacting partners. To gain new insight into how phosphorylation contributes to BAK1-mediated signaling specificity, we first mapped the phosphorylation patterns of BAK1 associated with different RLK partners (BRI1, FLS2, EFR and BIK1). Quantitative phospho-pattern profiling by label-free mass spectrometry revealed that differential phosphorylation patterns of RLK partners resulted from altered BAK1 phosphorylation status. More interestingly, the study of two BAK1 mutants (T450A and C408Y) both showing severe defect in immune defense yet normal growth phenotype suggested that varied phosphorylation patterns of RLK partners by BAK1 could be the molecular basis for selective regulation of multiple BAK1-dependent pathways. Taken together, this phospho-pattern profiling strategy allowed for explicit assessment of BAK1 kinase activity in different RLK complexes, which would facilitate elucidation of BAK1 diverse functions in plant development, defense, and adaptation.

BIOLOGICAL SIGNIFICANCE

BAK1 is a functionally important co-receptor known to interact with different receptor-like kinases (RLKs) to coordinate plant development and immune defenses. Our study first mapped the phosphorylation patterns of BAK1 associated with four RLK partners (BRI1, FLS2, EFR and BIK1), and further revealed that differential phosphorylation patterns of multiple RLK partners resulted from altered BAK1 phosphorylation status. More interestingly, the study of two BAK1 mutants suggested that varied phosphorylation patterns of RLK partners by BAK1 could be the basis for selective regulation of signaling pathways. Taken together, this phospho-pattern profiling strategy allowed for explicit assessment of BAK1 kinase activity in different RLK complexes, which would facilitate elucidation of BAK1 diverse functions in plant development, defense, and adaptation.

Single-step enrichment by Ti4+-IMAC and label-free quantitation enables in-depth monitoring of phosphorylation dynamics with high reproducibility and temporal resolution

Quantitative phosphoproteomics workflows traditionally involve additional sample labeling and fractionation steps for accurate and in-depth analysis. Here we report a high-throughput, straightforward, and comprehensive label-free phosphoproteomics approach using the highly selective, reproducible, and sensitive Ti(4+)-IMAC phosphopeptide enrichment method. We demonstrate the applicability of this approach by monitoring the phosphoproteome dynamics of Jurkat T cells stimulated by prostaglandin E2 (PGE2) over six different time points, measuring in total 108 snapshots of the phosphoproteome. In total, we quantitatively monitored 12,799 unique phosphosites over all time points with very high quantitative reproducibility (average r > 0.9 over 100 measurements and a median cv < 0.2). PGE2 is known to increase cellular cAMP levels, thereby activating PKA. The in-depth analysis revealed temporal regulation of a wide variety of phosphosites associated not only with PKA, but also with a variety of other classes of kinases. Following PGE2 stimulation, several pathways became only transiently activated, revealing that in-depth dynamic profiling requires techniques with high temporal resolution. Moreover, the large publicly available dataset provides a valuable resource for downstream PGE2 signaling dynamics in T cells, and cAMP-mediated signaling in particular. More generally, our method enables in-depth, quantitative, high-throughput phosphoproteome screening on any system, requiring very little sample, sample preparation, and analysis time.

Evaluation of quantitative performance of sequential immobilized metal affinity chromatographic enrichment for phosphopeptides

We evaluated a sequential elution protocol from immobilized metal affinity chromatography (SIMAC) employing gallium-based immobilized metal affinity chromatography (IMAC) in conjunction with titanium dioxide-based metal oxide affinity chromatography (MOAC). The quantitative performance of this SIMAC enrichment approach, assessed in terms of repeatability, dynamic range, and linearity, was evaluated using a mixture composed of tryptic peptides from caseins, bovine serum albumin, and phosphopeptide standards. Although our data demonstrate the overall consistent performance of the SIMAC approach under various loading conditions, the results also revealed that the method had limited repeatability and linearity for most phosphopeptides tested, and different phosphopeptides were found to have different linear ranges. These data suggest that, unless additional strategies are used, SIMAC should be regarded as a semiquantitative method when used in large-scale phosphoproteomics studies in complex backgrounds.

Automated, reproducible, titania-based phosphopeptide enrichment strategy for label-free quantitative phosphoproteomics

An automated phosphopeptide enrichment strategy is described using titanium dioxide (TiO2)-packed, fused silica capillaries for use with liquid chromatography (LC)-mass spectrometry (MS)/MS-based, label-free proteomics workflows. To correlate an optimum peptide:TiO2 loading ratio between different particle types, the ratio of phenyl phosphate-binding capacities was used. The optimum loading for the column was then verified through replicate enrichments of a range of quantities of digested rat brain tissue cell lysate. Fractions were taken during sample loading, multiple wash steps, and the elution steps and analyzed by LC-MS/MS to gauge the efficiency and reproducibility of the enrichment. Greater than 96% of the total phosphopeptides were detected in the elution fractions, indicating efficient trapping of the phosphopeptides on the first pass of enrichment. The quantitative reproducibility of the automated setup was also improved greatly with phosphopeptide intensities from replicate enrichments exhibiting a median coefficient of variation (CV) of 5.8%, and 80% of the identified phosphopeptides had CVs below 11.1%, while maintaining >85% specificity. By providing this high degree of analytical reproducibility, this method allows for label-free phosphoproteomics over large sample sets with complex experimental designs (multiple biological conditions, multiple biological replicates, multiple time-points, etc.), including large-scale clinical cohorts.

Longitudinal study of differential protein expression in an Alzheimer's mouse model lacking inducible nitric oxide synthase

Alzheimer's disease (AD) is a complex neurodegenerative process that involves altered brain immune, neuronal and metabolic functions. Understanding the underlying mechanisms has relied on mouse models that mimic components of AD pathology. We used gel-free, label-free LC-MS/MS to quantify protein and phosphopeptide levels in brains of APPSwDI/NOS2-/- (CVN-AD) mice. CVN-AD mice show a full spectrum of AD-like pathology, including amyloid deposition, hyperphosphorylated and aggregated tau, and neuronal loss that worsens with age. Tryptic digests, with or without phosphopeptide enrichment on an automated titanium dioxide LC system, were separated by online two-dimensional LC and analyzed on a Waters Synapt G2 HDMS, yielding relative expression for >950 proteins and >1100 phosphopeptides. Among differentially expressed proteins were known markers found in humans with AD, including GFAP and C1Q. Phosphorylation of connexin 43, not previously described in AD, was increased at 42 weeks, consistent with dysregulation of gap junctions and activation of astrocytes. Additional alterations in phosphoproteins suggests dysregulation of mitochondria, synaptic transmission, vesicle trafficking, and innate immune pathways. These data validate the CVN-AD mouse model of AD, identify novel disease and age-related changes in the brain during disease progression, and demonstrate the utility of integrating unbiased and phosphoproteomics for understanding disease processes in AD.

Identification and functional analysis of phosphorylation residues of the Arabidopsis BOTRYTIS-INDUCED KINASE1

Arabidopsis BOTRYTIS-INDUCED KINASE1 (BIK1) is a receptor-like cytoplasmic kinase acting early in multiple signaling pathways important for plant growth and innate immunity. It is known to form a signaling complex with a cell-surface receptor FLS2 and a co-receptor kinase BAK1 to transduce signals upon perception of pathogen-associated molecular patterns (PAMPs). Although site-specific phosphorylation is speculated to mediate the activation and function of BIK1, few studies have been devoted to complete profiling of BIK1 phosphorylation residues. Here, we identified nineteen in vitro autophosphorylation sites of BIK1 including three phosphotyrosine sites, thereby proving BIK1 is a dual-specificity kinase for the first time. The kinase activity of BIK1 substitution mutants were explicitly assessed using quantitative mass spectrometry (MS). Thr-237, Thr-242 and Tyr-250 were found to most significantly affect BIK1 activity in autophosphorylation and phosphorylation of BAK1 in vitro. A structural model of BIK1 was built to further illustrate the molecular functions of specific phosphorylation residues. We also mapped new sites of FLS2 phosphorylation by BIK1, which are different from those by BAK1. These in vitro results could provide new hypotheses for more in-depth in vivo studies leading to deeper understanding of how phosphorylation contributes to BIK1 activation and mediates downstream signaling specificity.

Proteomic analysis of human bronchoalveolar lavage fluid after subsgemental exposure

The analysis of airway fluid, as sampled by bronchoalveolar lavage (BAL), provides a minimally invasive route to interrogate lung biology in health and disease. Here, we used immunodepletion, coupled with gel- and label-free LC-MS/MS, for quantitation of the BAL fluid (BALF) proteome in samples recovered from human subjects following bronchoscopic instillation of saline, lipopolysaccharide (LPS) or house dust mite antigen into three distinct lung subsegments. Among more than 200 unique proteins quantified across nine samples, neutrophil granule-derived and acute phase proteins were most highly enriched in the LPS-exposed lobes. Of these, peptidoglycan response protein 1 was validated and confirmed as a novel marker of neutrophilic inflammation. Compared to a prior transcriptomic analysis of airway cells in this same cohort, the BALF proteome revealed a novel set of response factors. Independent of exposure, the enrichment of tracheal-expressed proteins in right lower lung lobes suggests a potential for constitutive intralobar variability in the BALF proteome; sampling of multiple lung subsegments also appears to aid in the identification of protein signatures that differentiate individuals at baseline. Collectively, this proof-of-concept study validates a robust workflow for BALF proteomics and demonstrates the complementary nature of proteomic and genomic techniques for investigating airway (patho)physiology.

Biochemical and biophysical changes underlie the mechanisms of basement membrane disruptions in a mouse model of dystroglycanopathy

Mutations in glycosyltransferases, such as protein O-mannose N-acetylglucosaminyltransferase 1 (POMGnT1), causes disruptions of basement membranes (BMs) that results in neuronal ectopias and muscular dystrophy. While the mutations diminish dystroglycan-mediated cell-ECM interactions, the cause and mechanism of BM disruptions remain unclear. In this study, we established an in vitro model to measure BM assembly on the surface of neural stem cells. Compared to control cells, the rate of BM assembly on POMGnT1 knockout neural stem cells was significantly reduced. Further, immunofluorescence staining and quantitative proteomic analysis of the inner limiting membrane (ILM), a BM of the retina, revealed that laminin-111 and nidogen-1 were reduced in POMGnT1 knockout mice. Finally, atomic force microscopy showed that the ILM from POMGnT1 knockout mice was thinner with an altered surface topography. The results combined demonstrate that reduced levels of key BM components cause physical changes that weaken the BM in POMGnT1 knockout mice. These changes are caused by a reduced rate of BM assembly during the developmental expansion of the neural tissue.

Circadian acetylome reveals regulation of mitochondrial metabolic pathways

The circadian clock is constituted by a complex molecular network that integrates a number of regulatory cues needed to maintain organismal homeostasis. To this effect, posttranslational modifications of clock proteins modulate circadian rhythms and are thought to convert physiological signals into changes in protein regulatory function. To explore reversible lysine acetylation that is dependent on the clock, we have characterized the circadian acetylome in WT and Clock-deficient (Clock(-/-)) mouse liver by quantitative mass spectrometry. Our analysis revealed that a number of mitochondrial proteins involved in metabolic pathways are heavily influenced by clock-driven acetylation. Pathways such as glycolysis/gluconeogenesis, citric acid cycle, amino acid metabolism, and fatty acid metabolism were found to be highly enriched hits. The significant number of metabolic pathways whose protein acetylation profile is altered in Clock(-/-) mice prompted us to link the acetylome to the circadian metabolome previously characterized in our laboratory. Changes in enzyme acetylation over the circadian cycle and the link to metabolite levels are discussed, revealing biological implications connecting the circadian clock to cellular metabolic state.

Proteomic analysis of ERK1/2-mediated human sickle red blood cell membrane protein phosphorylation

BACKGROUND

In sickle cell disease (SCD), the mitogen-activated protein kinase (MAPK) ERK1/2 is constitutively active and can be inducible by agonist-stimulation only in sickle but not in normal human red blood cells (RBCs). ERK1/2 is involved in activation of ICAM-4-mediated sickle RBC adhesion to the endothelium. However, other effects of the ERK1/2 activation in sickle RBCs leading to the complex SCD pathophysiology, such as alteration of RBC hemorheology are unknown.

RESULTS

To further characterize global ERK1/2-induced changes in membrane protein phosphorylation within human RBCs, a label-free quantitative phosphoproteomic analysis was applied to sickle and normal RBC membrane ghosts pre-treated with U0126, a specific inhibitor of MEK1/2, the upstream kinase of ERK1/2, in the presence or absence of recombinant active ERK2. Across eight unique treatment groups, 375 phosphopeptides from 155 phosphoproteins were quantified with an average technical coefficient of variation in peak intensity of 19.8%. Sickle RBC treatment with U0126 decreased thirty-six phosphopeptides from twenty-one phosphoproteins involved in regulation of not only RBC shape, flexibility, cell morphology maintenance and adhesion, but also glucose and glutamate transport, cAMP production, degradation of misfolded proteins and receptor ubiquitination. Glycophorin A was the most affected protein in sickle RBCs by this ERK1/2 pathway, which contained 12 unique phosphorylated peptides, suggesting that in addition to its effect on sickle RBC adhesion, increased glycophorin A phosphorylation via the ERK1/2 pathway may also affect glycophorin A interactions with band 3, which could result in decreases in both anion transport by band 3 and band 3 trafficking. The abundance of twelve of the thirty-six phosphopeptides were subsequently increased in normal RBCs co-incubated with recombinant ERK2 and therefore represent specific MEK1/2 phospho-inhibitory targets mediated via ERK2.

CONCLUSIONS

These findings expand upon the current model for the involvement of ERK1/2 signaling in RBCs. These findings also identify additional protein targets of this pathway other than the RBC adhesion molecule ICAM-4 and enhance the understanding of the mechanism of small molecule inhibitors of MEK/1/2/ERK1/2, which could be effective in ameliorating RBC hemorheology and adhesion, the hallmarks of SCD.

Identification of deacetylase substrates with the biotin switch approach

The identification of lysine-acetylated proteins and deacetylase substrates has primarily relied on protein immune-affinity techniques with antibodies that recognize acetylated lysine residues (Kac antibodies). While these antibody-based techniques are continuously improving, they can be limited by the narrow and many times unknown epitope specificity of Kac antibodies. An alternative approach is the biotin switch capture of deacetylated proteins. Similar in part to other biotin switch methodologies, this technique relies on the blocking of native lysine residues and removal of the modification of interest in vitro, after which the newly deacetylated proteins can be captured and identified by mass spectrometry (MS). In this chapter, we cover the essential steps of the procedure, highlight key points in the assay to reduce false positive protein identification, and discuss the quantitative MS methods useful for identifying the captured deacetylase substrates. We also discuss potential strategies and future improvements to overcome current limitations of the assay.

Solid-phase capture for the detection and relative quantification of S-nitrosoproteins by mass spectrometry

The proteomic analysis of S-nitrosylated protein (SNO-proteins) has long depended on the biotin switch technique (BST), which requires blocking of free thiols, ascorbate-based denitrosylation of SNO-Cys, biotinylation of nascent thiol and avidin-based affinity isolation. A more recent development is resin assisted-capture of SNO-proteins (SNO-RAC), which substitutes thiopropyl Sepharose (TPS) for biotin-avidin, thus reducing the number of steps required for enrichment of S-nitrosylated proteins. In addition, SNO-RAC facilitates on-resin proteolytic digestion following SNO-protein capture, greatly simplifying the purification of peptides containing sites of S-nitrosylation ("SNO-sites"). This resin-based approach has also now been applied to detection of alternative Cys-based modifications, including S-palmitoylation (Acyl-RAC) and S-oxidation (Ox-RAC). Here, we review the important steps to minimize false-positive identification of SNO-proteins, give detailed methods for processing of protein-bound TPS for mass spectrometry (MS) based analysis, and discuss the various quantitative MS methods that are compatible with SNO-RAC. We also discuss strategies to overcome the current limitations surrounding MS-based SNO-site localization in peptides containing more than one potential target Cys residue. This article therefore serves as a starting point and guide for the MS-focused exploration of SNO-proteomes by SNO-RAC.

Protein modification by deamidation indicates variations in joint extracellular matrix turnover

As extracellular proteins age, they undergo and accumulate nonenzymatic post-translational modifications that cannot be repaired. We hypothesized that these could be used to systemically monitor loss of extracellular matrix due to chronic arthritic diseases such as osteoarthritis (OA). To test this, we predicted sites of deamidation in cartilage oligomeric matrix protein (COMP) and confirmed, by mass spectroscopy, the presence of deamidated (Asp(64)) and native (Asn(64)) COMP epitopes (mean 0.95% deamidated COMP (D-COMP) relative to native COMP) in cartilage. An Asp(64), D-COMP-specific ELISA was developed using a newly created monoclonal antibody 6-1A12. In a joint replacement study, serum D-COMP (p = 0.017), but not total COMP (p = 0.5), declined significantly after replacement demonstrating a joint tissue source for D-COMP. In analyses of 450 participants from the Johnston County Osteoarthritis Project controlled for age, gender, and race, D-COMP was associated with radiographic hip (p < 0.0001) but not knee (p = 0.95) OA severity. In contrast, total COMP was associated with radiographic knee (p < 0.0001) but not hip (p = 0.47) OA severity. D-COMP was higher in soluble proteins extracted from hip cartilage proximal to OA lesions compared with remote from lesions (p = 0.007) or lesional and remote OA knee (p < 0.01) cartilage. Total COMP in cartilage did not vary by joint site or proximity to the lesion. This study demonstrates the presence of D-COMP in articular cartilage and the systemic circulation, and to our knowledge, it is the first biomarker to show specificity for a particular joint site. We believe that enrichment of deamidated epitope in hip OA cartilage indicates a lesser repair response of hip OA compared with knee OA cartilage.

Phosphoproteome exploration reveals a reformatting of cellular processes in response to low sterol biosynthetic capacity in Arabidopsis

Sterols are membrane-bound isoprenoid lipids that are required for cell viability and growth. In plants, it is generally assumed that 3-hydroxy-3-methylglutaryl-CoA-reductase (HMGR) is a key element of their biosynthesis, but the molecular regulation of that pathway is largely unknown. In an attempt to identify regulators of the biosynthetic flux from acyl-CoA toward phytosterols, we compared the membrane phosphoproteome of wild-type Arabidopsis thaliana and of a mutant being deficient in HMGR1. We performed a N-terminal labeling of microsomal peptides with a trimethoxyphenyl phosphonium (TMPP) derivative, followed by a quantitative assessment of phosphopeptides with a spectral counting method. TMPP derivatization of peptides resulted in an improved LC-MS/MS detection due to increased hydrophobicity in chromatography and ionization efficiency in electrospray. The phosphoproteome coverage was 40% higher with this methodology. We further found that 31 proteins were in a different phosphorylation state in the hmgr1-1 mutant as compared with the wild-type. One-third of these proteins were identified based on novel phosphopeptides. This approach revealed that phosphorylation changes in the Arabidopsis membrane proteome targets major cellular processes such as transports, calcium homeostasis, photomorphogenesis, and carbohydrate synthesis. A reformatting of these processes appears to be a response of a genetically reduced sterol biosynthesis.

Tunnel frit: a nonmetallic in-capillary frit for nanoflow ultra high-performance liquid chromatography-mass spectrometryapplications

In this study, an easy method to fabricate a durable in-capillary frit was developed for use in nanoflow liquid chromatography (nanoLC). A small orifice was tunneled into the sol-gel frit during the polymerization process resulting in the simple fabrication of a tunnel frit. A short packing tunnel frit column (2 cm, C(18) particles) was able to sustain over 10,000 psi continuous liquid flow for 10 days without observation of particle loss, and back pressure variation was less than 5%. The tunnel frit was successfully applied to the fabrication of nanoflow ultra high-performance liquid chromatography (nano-UHPLC) trap and analytical columns. In the analysis of tryptic peptides, the tunnel frit trap and analytical columns were demonstrated to have high separation efficiency and sensitivity. In analysis of phosphopeptides, the use of the nonmetallic tunnel frit column showed better sensitivity than the metallic frit column. This design can facilitate the preparation of nano-HPLC and nano-UHPLC columns and the packing material can easily be refilled when the column is severely contaminated or clogged.