Comprehensive Analysis of in Vivo Phosphoproteome of Mouse Liver

Oct 20, 2015 - In the KEGG pathway analysis, phosphorylated microsomal proteins were highly enriched in ribosome protein processing in ER and ribosome...
0 downloads 0 Views 1MB Size
Subscriber access provided by CMU Libraries - http://library.cmich.edu

Article

Comprehensive analysis of in vivo phosphoproteome of mouse liver microsomes Oh Kwang Kwon, Juhee Sim, Sun Ju Kim, Eunji Sung, Jin Young Kim, Tae Cheon Jeong, and Sangkyu Lee J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/acs.jproteome.5b00812 • Publication Date (Web): 20 Oct 2015 Downloaded from http://pubs.acs.org on October 29, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Proteome Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Comprehensive analysis of in vivo phosphoproteome of mouse liver microsomes

Oh Kwang Kwon †,¶, JuHee Sim †,¶, Sun Ju Kim †, Eunji Sung †, Jin Young Kim ‡, Tae Cheon Jeong §, Sangkyu Lee *,†



College of Pharmacy, Research Institute of Pharmaceutical Sciences, Kyungpook National

University, Daegu 41566, Republic of Korea ‡

Mass Spectrometry Research Center, Korea Basic Science Institute, Ochang, Chungbuk

28115, Republic of Korea §

College of Pharmacy, Yeungnam University, Gyeongsan 38541, Republic of Korea



Those authors are equally contributed.

*Correspondence: Sangkyu Lee, Ph.D., Professor, College of Pharmacy, Research Institute of Pharmaceutical Sciences, Kyungpook National University, 80 Daehak-ro, Buk-gu 702-701, Republic

of

Korea.

Tel.:

+82-53-950

8571;

Fax:

[email protected].

1

ACS Paragon Plus Environment

+82-53-950-8557;

E-mail:

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Protein phosphorylation at serine, threonine, and tyrosine residues are some of the most widespread reversible post-translational modifications. Microsomes are vesicle-like bodies, not ordinarily present within living cells, which form from pieces of the endoplasmic reticulum (ER), plasma membrane, mitochondria, or Golgi apparatus of broken eukaryotic cells. Here, we investigated the total phosphoproteome of mouse liver microsomes (MLM) using TiO2 enrichment of phosphopeptides coupled with in-line two-dimensional (2D)-LC-MS/MS. In total, 699 phosphorylation sites in 527 proteins were identified in MLMs. When compared with the current phosphoSitePlus database, 155 novel phosphoproteins were identified in MLM. The distributions of phosphosites were 89.4%, 8.0%, and 2.6% for phosphoserine, phosphotheronine, and phosphotyrosine, respectively. By Motif-X analysis, 8 Ser motifs and 1 Thr motif were found, and 5 acidic, 2 basophilic-, and 2 proline-directed motifs were assigned. The potential functions of phosphoproteins in MLM were assigned by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. In GO annotation, phosphorylated microsomal proteins were involved in mRNA processing, mRNA metabolic processes and RNA splicing. In the KEGG pathway analysis, phosphorylated microsomal proteins were highly enriched in ribosome protein processing in ER, ribosomes, and in RNA transport. We determined that 52 and 23 phosphoproteins were potential substrates of cAMP-dependent protein kinase A and casein kinase II, respectively, many of which are 40S/60S ribosomal proteins. Overall, our results provide an overview of features of protein phosphorylation in MLMs that should be a valuable resource for the future understanding of protein synthesis or translation involving phosphorylation.

Key words Microsomes; phosphorylation; ribosomal protein; cytochrome P450; protein processing 2

ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Introduction Phosphorylation of proteins at serine, threonine, and tyrosine residues is one of the most widespread reversible post-translational modifications (PTM)

1

. Generally, reversible

phosphorylation, a key mechanism of regulating protein activities, participates in diverse cellular processes including metabolism, cell communication, cell growth, and development 2. The balance between reversible phosphorylation and dephosphorylation is coordinated by the dynamic interplay between protein kinases, phosphatases, and phosphorus-binding proteins. Although more than 30% of human proteins are known to be phosphorylated, phosphorylation of regulatory proteins is a transient modification that shows low stoichiometry

3, 4

. Phosphorylation occurs only on particular proteins or at particular times,

resulting in a low steady state abundance of phosphorylated proteins. Therefore, specific and efficient phosphopeptide-enrichment methods are needed to capture and analyze the large scale of the phosphoproteome in complex samples. The enzymes found in liver microsomes normally contribute to various biological reactions, including metabolism of chemicals, the various functions of membrane bound enzymes, and lipid-protein interactions

5, 6

. Microsomes are spontaneously formed vesicle-

like bodies, not ordinarily present in living cells, derived from pieces of the endoplasmic reticulum (ER), plasma membrane, mitochondria, or Golgi apparatus of broken eukaryotic cells 7. Typically, microsomes are separated by differential centrifugation of homogenized cells or tissues. Using methods based on enrichment of ER fractions, prepared microsomes usually show a high level of metabolic activity due to large amounts of metabolic enzymes, such as cytochrome P450 (CYP) 8. CYPs are important components of xenobioticmetabolizing

monooxygenase,

which

accounts

biotransformation of clinically relevant drugs

for

approximately

9

75%

of

the

. Hepatic CYPs are ER-anchored

3

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 34

hemoproteins involved in the metabolism of numerous endo- and xenobiotics, including drugs, steroids, and carcinogens 10. During the past decades, global phosphoproteomic analyses have been conducted in H. sapiens and M. musculus using MS-based proteomic and phosphopeptide enrichment methods 1, 11, 12. Nevertheless, only proteomic analyses have been applied to microsomes, and the extent of phosphorylation of microsomal proteins remains undefined

6, 13, 14

. A better

understanding of protein phosphorylation in microsomes is needed to improve our knowledge of important microsomal processes. Here we investigated the global microsomal phosphoproteome using TiO2 enrichment of phosphopeptides coupled with in-line twodimensional (2D)-LC-MS/MS.

4

ACS Paragon Plus Environment

Page 5 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Experimental procedures Preparation of mouse liver microsomes Specific pathogen-free male ICR mice (28 to 30 g) were obtained from the Orient Co. (Seoul, Korea). The animals purchased at 5 weeks of age were acclimated for at least 2 weeks. They were randomly housed four per cage and strictly maintained at 23 ± 3 °C and 40–60% relative humidity. A 12 h light/12 h dark cycle was used with an intensity of 150–300 lx. After animals were subjected to necropsy, the livers were removed and homogenized with four volumes of ice-cold 0.1 M potassium phosphate buffer, pH 7.4. The liver homogenates were centrifuged at 9,000 ×g for 20 min at 4 °C in order to obtain an S9 fraction. The S9 fraction was centrifuged at 105,000 ×g for 1 h at 4 °C using an ultracentrifuge, and the supernatants were discarded. The liver microsomal pellets were washed by pipetting using ice-cold 0.1 M potassium phosphate buffer, pH 7.4, and were centrifuged once again under the same conditions. The microsomal pellets were ground on ice in 250 mM sucrose in 0.1 M potassium phosphate buffer by using a Dounce tissue grinder. The prepared micosome were confirmed by Western blotting (Supplemental Figure S1). The liver microsomal pellets isolated by differential centrifugation were stored at -80 °C until use

15, 16

. The concentration

of microsomal protein was determined using Bradford assay (Bio-Rad Laboratories, Inc., Hercules, CA). All animal procedures complied with the guidelines issued by the Society of Toxicology (USA; 1989). This study was approved by the Institutional Review Board of the Kyungpook National University (IRB-KNU 2014-0082-1).

Tryptic digestion of microsomal proteins To analyze the proteins, the microsome fractions were lysed using lysis buffer at 4 °C for 5

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

30 min and centrifuged at 4 °C for 10 min. Lysis buffer contained RIPA buffer (25 mM TrisHCl, pH 7.6, 150 mM NaCl, 1% NP-40, 1% sodium deoxycholate, 0.1% SDS, Thermo Scientific, Rockford, IL) with protease inhibitor cocktail (1 mM AEBSF, 0.8 µM aprotinin, 50 µM bestatin, 15 µM E-64, 20 µM leupeptin, and 10 µM pepstatin A; Calbiochem, Merck Millipore, Billerica, MA) and phosphatase inhibitor (10 mM sodium fluoride, 1 mM sodium orthovanadate, 2 mM β-glycerophosphate disodium salt hydrate, 2 mM sodium pyrophosphate decahydrate, Sigma Aldrich, St. Louis, MO). The supernatants were transferred to new tubes and precipitated with 10% trichloroacetic acid (TCA) at 4 °C for 4 h. The precipitated protein was centrifuged at 13,000 ×g, 4 °C for 10 min, and washed twice with cold acetone. Two milligrams of protein were reduced using 15 mM dithiothreitol in 25 mM ammonium bicarbonate at 56 °C for 30 min, and then alkylated using 60 mM iodoacetamide in 25 mM ammonium bicarbonate for 30 min at RT in the dark. Protein samples were then added to 75 mM cysteine in 25 mM ammonium bicarbonate at 25 °C for 30 min in the dark. After adding trypsin to the protein (at a ratio of 1:50 w/w), peptide samples were digested at 37 °C overnight. Digested peptides were dried using a centrifugal vacuum concentrator, and then additional trypsin was added at 1:50 (w/w). Finally, the mixed peptides were dried using a centrifugal vacuum concentrator.

Enrichment of phosphopeptides The phosphopeptides were captured using a phosphopeptide enrichment titanium dioxide (TiO2) kit (TitansphereTM, Phos-TiO kit, GL Sciences Inc., Japan). The TiO2 tips were activated with 200 µl of 100% ACN and equilibrated with 200 µl of solution B (1 M glycolic acid in 80% ACN with 2% TFA) using centrifugation at 3,000 ×g for 2 min. For adsorption 6

ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

of peptides, peptides were mixed with 1000 µl of solution B at RT, rinsed using solution B and then with solution A (80% ACN with 2% TFA), each at 1,000 ×g for 5 min. Phosphopeptides were eluted with 200 µl of 0.5% ammonium hydroxide solution, and with a 200 µl volume of 30% ACN for an additional elution. Enriched phosphopeptides were dried using a centrifugal vacuum concentrator. The samples were desalted using C18 Zip-tip (Thermo Pierce, Rockford, IL) according to manufacturer’s instructions and dried using a centrifugal vacuum concentrator. The sample was analyzed by 2D-LC-MS/MS performed on a nano ACQUITY Ultra Performance LC System (Waters; Milford, MA) and an LTQOrbitrap velos pro (Thermo Scientific; Rockford, IL) fitted with a nano-electrospray source. The details of the LC-MS/MS method are described in Supporting Information.

Western blotting The population of phosphoproteins in microsomes was screened by Western blotting. The microsomes (30 µg) were separated by 12% SDS-PAGE and transferred to a PVDF membrane (Roche, Germany) by wet transfer at 300 mA for 1 h on ice. After membranes were blocked with TBST (20 mM Tris, 137 mM sodium chloride, 0.1% Tween-20) containing 5% bovine serum albumin (BSA) for 3 h at RT, membranes were incubated with each primary

antibody;

anti-phosphoserine

(Abcam,

UK),

anti-phosphothreonine,

anti-

phosphotyrosine (Cell Signaling, Danvers, MA) by rotation at 4 °C overnight. Membranes were washed with TBST, exposed to horseradish peroxidase (HRP)-linked secondary antibody (Cell Signaling, Danvers, MA) with rotation for 2 h at RT, and washed with TBST. The membranes were developed using ECL reagent (GE Healthcare, UK) and detected using Image Quant LAS 4000 mini (GE Healthcare, UK).

7

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 34

PKA and CKII enzyme reactions For reaction with cAMP-dependent protein kinase A (PKA) using a PKA assay kit (Promega, Madison, WI), 500 µg of pooled CD-1 male mouse liver microsome (MLM, XenoTech, LLC., Lenexa, USA) were diluted in 5× PKA reaction buffer (100 mM Tris-HCl, pH 7.4, 50 mM MgCl2, 5 mM ATP), plus PKA activator 5× solution (5 µM cAMP in water), and preincubated at 30 °C for 1 min. To the MLM reaction samples were added PKA solution (2 µg/ml PKA in 350 mM K3PO4, 0.1 mM DTT) or buffer only without PKA solution (negative control, only added 350 mM K3PO4, 0.1 mM DTT), and incubation was continued at 30 °C for 2.5 h. For reaction with casein kinase II (CKII) using recombinant CKII (New England BioLabs Inc., Ipswich, MA), 500 µg of pooled MLM were diluted in CKII reaction buffer (20 mM Tris-HCl, pH 7.5, 50 mM KCl, 10 mM MgCl2) with 0.2 mM ATP solution in water and preincubated at 30 °C for 1 min. CKII (2500 U/5 µL) was added and incubation continued at 30 °C for 2.5 h. Both kinase reactions were incubated in the presence of phosphatase inhibitor cocktail (10 mM sodium fluoride, 1 mM sodium orthovanadate, 2 mM β-glycerophosphate disodium salt hydrate, 2 mM sodium pyrophosphate decahydrate). For immunoblotting, approximately 10 µg of the reaction samples were transferred to new tubes and 5× SDS sample buffer was added. Reaction samples (400 µg) were treated with lysis buffer to obtain trypsin-digested peptides.

Identification and validation of phosphorylation Phosphopeptides and proteins were identified using MaxQuant (version 1.4.1.2)

17

. The

resulting MS/MS spectral data were used to search a database containing 45,269 protein sequences

of

Mus

musculus

downloaded

from

UniProtKB

(http://www.uniprot.org/proteomes/UP000000589). Search parameters were as follows: 8

ACS Paragon Plus Environment

Page 9 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

trypsin digestion with up to 2 missed cleavages; fixed modification of carbamidomethylation (Cys); variable modification of oxidation (Met) and phosphorylation (Ser/Thr/Tyr). Maximum mass tolerance allowance was 5 ppm (monoisotopic) for MS/MS tolerance and 0.02 Da (monoisotopic) for fragment ions. MaxQuant search results were filtered with a MaxQuant score ≥40 and occupation of phosphorylation site probability ≥ 0.5, with reverse and contamination peptides removed from the results. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD003043 18.

Bioinformatic analysis To obtain information about the identified phosphoproteins, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and the annotated protein KEGG database were searched. For GO annotation, phosphoproteins were mapped to GO ID. If the identified phosphorylation substrates were not annotated in the UniProt or GOA databases, InterProScan software was used to annotate a protein’s GO function based on a protein sequence alignment method 19. The annotated phosphoproteins were further classified by GO annotation based on three categories: biological process, cellular component, and molecular function. The pathways associated with phosphoproteins in microsomes were annotated with KEGG online service tool KAAS. The annotation results were mapped onto the KEGG pathway database using the KEGG online service tool KEGG mapper

20

. WoLF

PSORT was used for subcellular localization prediction 21. For functional enrichment analysis, Fisher’s exact test was used to test for enrichment or depletion (right-tailed test) of specific annotation terms among members of the resulting 9

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

protein clusters. Derived p-values were further adjusted to address multiple hypothesis testing by the method proposed by Benjamini and Hochberg

22

. Terms having adjusted p-values

below 0.05 in any of the clusters were treated as significant. For motif analysis, Software motif-x was used to analyze the models of sequences composed of amino acids in specific positions within phospho-13-mers (6 amino acids upstream and downstream of the phosphorylation site) in all protein sequences 23. Further, all the database protein sequences were used as the background database parameter, with all other parameters set at default.

10

ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Results Global phosphoproteome in mouse liver microsomes The liver microsomal lysates were prepared from ICR male mice, with 3 individual livers pooled prior to sample digestion. The separated microsomes had been confirmed by Western blotting using CYP2B and cytochrome c oxidase subunit IV (COX IV) indicating specific microsomal and mitochondrial proteins, respectively (Supplemental Figure S1). All samples had strong blots for CYP2B, on the other hand COX IV was not showed in prepared samples. To confirm the distribution of microsomal phosphorylation, microsomal proteins were separated by 1D SDS-PAGE (Figure 1). The Coomassie blue-stained gel showed a closely migrating set of proteins in the range of 43–55 kDa, typical of the migration of CYP proteins 24

(Figure 1A). The distribution of phosphorylations in prepared microsomes was detected by

Western blot using anti-phospho-Ser, phospho-Thr, and phospho-Tyr antibodies (Figure 1B). Analogous patterns of immunoblots with well-characterized patterns were obtained from three individual microsome preparations. To identify the global phosphoproteome in mouse liver microsomes, microsomal lysates were prepared as solution tryptic digests (Figure 2A). After trypsin digestion, the peptide mixture was enriched for phosphopeptides using TiO2 chromatography. Enriched phosphopeptides were analyzed by nano-LC-MS/MS on Orbitrap velos pro. Phosphoproteins were identified by MaxQuant with the use of the current UniProt database (containing 45,269 protein sequences). To ensure the quality of the analysis, annotated spectra of all detected phosphopeptides were examined manually, and the mass tolerance for precursor ions was less than 5 ppm (Supplemental Figure S2A). Using orthogonal filtering criteria, we filtered the search results sequentially and established a phosphopeptide data set with a MaxQuant score ≥ 40, modified sites probability ≥ 0.5, and removal of reverse and contaminating peptides (Supplemental 11

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 34

Figure S2B). In total, 699 phosphorylation sites in 527 phosphorylated proteins were identified in mouse liver microsomal fractions (Supplemental Table S1). 372 phosphoproteins were identified in the phosphoSitePlus database

25

, and 155 (29%) novel phosphoproteins

were identified for the first time in mouse liver microsomes (Figure 2B). Detailed information of the novel phosphorylated peptides and proteins are summarized in Supplemental Table S2.

Characterization of phosphorylation in mouse liver microsomes The detected distributions of phosphosites were 89.4%, 8.0%, and 2.6% for phosphoserine, phosphothreonine, and phosphotyrosine, respectively (Figure 2C). The high ratio of phosphoserine was close to that observed in a large-scale phosphoproteomics study of a mouse liver cell line 12, 26.

In this study, the Motif-X algorithm was used to extract specific

microsomal phosphorylation motifs from the data (Figure 3). The kinase classification for these motifs was confirmed by reference to a previously reported motif database

27

. As a

result, 8 Ser motifs and 1 Thr motif were found, and 5 acidic (n = 178), 2 basophilic (n = 108) and 2 proline-directed motifs (n = 190) were assigned. The acidic motifs containing E or D at positions +1, +2, or +3 were dominant, consistent with reactions by CKII 4. We identified 2 basophilic motifs (RxxS), which are specific target sites of PKA or protein kinase

28

, and 2

proline-directed motifs (SP and TP) found as substrates of mitogen-activated protein kinase (MAPK) 4, 29. To understand potential functional implications of phosphorylation in the microsomal fractions, we analyzed enrichments of GO annotation and KEGG pathway analysis (Figure 4 and Supplemental Figure S3). GO analysis showed that microsomal phosphorylations were enriched for proteins that are involved in mRNA processing, mRNA metabolic processes, and 12

ACS Paragon Plus Environment

Page 13 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

RNA splicing in biological processes (Figure 4A). With respect to cellular compartment, the phosphorylations were mainly detected in membrane proteins associated with the ER. Phosphorylation

of

microsomal

proteins

is

therefore

functionally

important

for

oxidoreductase activity and RNA binding. In the KEGG pathway analysis, the most highly enriched categories for microsomal phosphorylation were ribosome protein processing in the ER, ribosomes, and RNA transport (Figure 4B). The greatest enrichment of microsomal phosphorylation was in proteins involved in processing in the ER and translation of mRNA (Supplemental Figures S4 and S5). A total of 27 phosphosites in 18 phosphoproteins related to protein processing was identified, and three phosphorylation sites in ribosome-binding protein 1 were revealed as novel phosphorylations in microsomal proteins (Table 1). Phosphorylations of CYP isoforms were previously identified in humans 24. Specifically, the phosphorylation at S478 of hCYP3A4 is well-known to enhance the ubiquitinationmediated degradation of CYP

30

. Previous results from studies that focused on the

identification of phosphorylation of human and murine CYP were poorly analyzed. In this study, we identified a total of 10 phosphorylation sites in 8 murine CYP proteins (Table 2). Four phosphorylation sites, S128, T275, S134, and S138 in CYP 2B10, 2D10, 3A11, and 3A13, respectively, had not been reported previously, and the S129 phosphorylation of CYP2E1 had been previously identified as a phosphorylation site in humans 24.

Identification of potential PKA and CKII kinase substrates by motif analysis Based on identified motifs, we hypothesized that PKA and CKII might be the main kinases regulating the microsomal phosphorylation. To investigate the substrates of PKA and CKII, mouse microsomal proteins were incubated with or without PKA or CKII kinase. We confirmed an increased phosphorylation of microsomal proteins using immunoblotting, 13

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

except Tyr-phosphorylation (Figure 5). The 52 and 23 phosphoproteins as potential substrates of PKA and CKII, respectively, were identified and are listed in Supplemental Tables S3 and S4. The 17 and 4 phosphopeptides in 40S/60S ribosomal proteins that were identified as substrates of PKA and CKII kinase in mouse liver microsomes are listed in Table 3. When compared with the PhosphoSite database, 8 phosphorylation sites were revealed as novel phosphorylation sites at positions S272, T18, T60, S9, S82, S31, S53, and S32 in 60S ribosomal proteins L5, S16, S20, L23, S25, S29, L35, and L36a, respectively.

14

ACS Paragon Plus Environment

Page 14 of 34

Page 15 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Discussion From a global phosphoproteome analysis of mouse liver microsomes, an enrichment of phosphorylation in ribosomal proteins and proteins involved in processing was revealed. Specifically, we identified phosphorylations at 27 sites in 18 proteins involved in protein processing in microsomal proteins. Among these, calnexin is a molecular chaperone of ER with folding-promoting functions and an abundant integral membrane phosphoprotein of the ER of eukaryotic cells

31, 32

. We verified the 3 phosphorylation sites on calnexin that were

previously identified as phosphorylated sites at Ser553, 563, and 582 33. Several kinases have been shown to have roles in the phosphorylation of calnexin. The phosphorylation at the Ser563 in calnexin by CK2 increases its association with ribosomes

34

. ERK1, MAPK was

also shown to phosphorylate Ser563 both in vitro and in vivo 35. In this study, the previously known 3 phosphorylation sites in calnexin were confirmed, indicating the reliability of the characterized phosphorylation sites in microsomal proteins we have identified. Here we also identified 3 phosphorylation sites at Ser154, 155, and 167 in ribosomebinding protein 1 (RRBP1). RRBP1 is a membrane-associated and ribosome-binding protein located on the rough ER that has essential roles in the transportation of intracellular proteins in mammalian cells

36

. The abnormal expression of RRBP1 is related to tumorigenesis and

progression of colorectal and lung cancers 37, 38. In breast cancer, RRBP1 is upregulated at the mRNA level and overexpressed in 84% (177/219) of the breast carcinoma cases tested 39. In addition, RRBP1 overexpression highly affects overall survival in patients with early-stage (I and II) tumors, indicating that RRBP1 is a potentially important marker for disease prognosis 40

. In general, the phosphorylated forms of proteins are typically associated with biologically

activated forms, as compared with non-phosphorylated proteins 2. Therefore, our findings suggest that the phosphorylation of RRBP1 is an interesting molecule that may be involved in 15

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 34

progression of cancer, and the potential role of phosphorylation in RRBP1 can be identified in further studies to serve as a cancer biomarker or prognostic factor. CYPs are categorized into 1A, 2A, 2B, 2C, 2D, 2E, and 3A subfamilies, which have important roles in drug metabolism

41

. There are 57 human and 102 murine putatively

functional genes 42. In 2003, PKA was identified as a major contributor of phosphorylation of CYP

43

. The phosphorylation of CYP leads to a marked decrease in the monooxygenase

activity associated with UBC7/gp78–mediated CYP3A4 and CYP2E1 ubiquitination

30, 44, 45

.

Its cell surface receptor, gp78/AMFR, is also localized to the endoplasmic reticulum, where it functions as an E3 ubiquitin ligase

46

. The turnover of CYP3A4 involves endoplasmic

reticulum-associated degradation via the ubiquitin (Ub)-dependent 26S proteasomal system that relies on two highly complementary E2 Ub-conjugating-E3 Ub-ligase complexes, as well as PKA and PKC

47

. It is therefore important to identify the phosphorylation sites of CYPs

since the phosphorylation of CYP is involved in the regulation of multiple metabolic activities. Here, we identified the 10 phosphorylation sites of murine cytochrome P450s using a global phosphoproteome approach in mouse liver microsomes. As shown in Table 2, the phosphorylation of CYP was highly similar in humans and mice. We identified a consensus phosphorylated sequence [RFS(ph)XXXX] in humans and mice for CYP2A12, 2B10, 2D26, and 2E1 24. The phosphorylation at S134 on CYP3A11 also showed the same result as seen in hCYP3A4 30. The phosphorylation in CYPs is significantly involved in regulation of catalytic activities, and the variation of CYP activities may seriously affect results of pharmacokinetic studies during drug development processes. In this manuscript, a high similarity of phosphorylation patterns between human and mouse CYP isoforms was revealed, suggesting that it will be helpful to search for correlations between results in human and mouse 16

ACS Paragon Plus Environment

Page 17 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

pharmacokinetic studies. In mouse liver microsomes, PKA mediated a total of 66 phosphorylations, and 17 phosphorylation sites found in 40S/60S ribosomal proteins that comprise the ribosome complex were mediated by PKA’s catalytic activity (Table 3). PKA is the archetypical phosphokinase, sharing a catalytic core with the entire protein kinase superfamily

48

.

Moreover, PKA contributes to peptide-chain elongation in mammalian cells by phosphorylating the eukaryotic translation elongation factor 2 (eEF2) kinase

49

. The

phosphorylation of eEF2 kinase causes phosphorylation and inactivation of eEF2, resulting in a decrease in the elongation of peptide chains

50

. Although PKA contributes to diverse

cellular biological functions, such as cellular signaling, cell cycle, proliferation, and metabolism

51, 52

, its interaction with ribosomal proteins was not well-known. Moreover, the

function of phosphorylation of ribosomal proteins in the ribosome complex is also not clear. Here we showed that the subunit proteins of ribosomes were widely phosphorylated by PKA, and the modified sites in phosphoproteins were detected in mouse liver microsomes. The modulation of PKA activity might be an essential regulating factor of proteins in ribosomes or ER membrane. The biological functions of the phosphorylations of ribosome subunits should be the subjects of further study. In present study, we identified global levels of phosphorylation in mouse liver microsomes that were focused on ribosomal proteins and CYPs. 155 phosphoproteins were newly identified, and 10 phosphorylation sites in CYPs and a large population of phosphorylation sites in ribosomal proteins were also investigated. Use of the mouse as a model mammal, and as a surrogate for human biology, assumes reasonable similarity between the two. It is therefore of interest to compare the similarities and differences of protein phosphorylations, and to clarify the limits of extrapolation from mouse to human. Therefore, 17

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 34

it will be very important to expand on in vivo experiments that generate information about phosphorylation in murine models, as was presented here.

Acknowledgements This study was supported by the Korean Health Technology R&D Project, Ministry of Health & Welfare, Republic of Korea (Grant No: A112026) and by the National Research Foundation

(NRF)

grant

funded

by

the

Korean

government

2012R1A4A1028835).

18

ACS Paragon Plus Environment

(MSIP)

(NRF-

Page 19 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

References 1.

Olsen, J. V.; Blagoev, B.; Gnad, F.; Macek, B.; Kumar, C.; Mortensen, P.; Mann, M.,

Global, in vivo, and site-specific phosphorylation dynamics in signaling networks. Cell 2006, 127, (3), 635-48. 2.

Kweon, H. K.; Andrews, P. C., Quantitative analysis of global phosphorylation

changes with high-resolution tandem mass spectrometry and stable isotopic labeling. Methods 2013, 61, (3), 251-9. 3.

Cohen, P., The origins of protein phosphorylation. Nat Cell Biol 2002, 4, (5), E127-

30. 4.

Olsen, J. V.; Vermeulen, M.; Santamaria, A.; Kumar, C.; Miller, M. L.; Jensen, L. J.;

Gnad, F.; Cox, J.; Jensen, T. S.; Nigg, E. A.; Brunak, S.; Mann, M., Quantitative phosphoproteomics reveals widespread full phosphorylation site occupancy during mitosis. Sci Signal 2010, 3, (104), ra3. 5.

Heinemann, F. S.; Ozols, J., Isolation and structural analysis of microsomal

membrane proteins. Front Biosci 1998, 3, d483-93. 6.

Peng, F.; Zhan, X.; Li, M. Y.; Fang, F.; Li, G.; Li, C.; Zhang, P. F.; Chen, Z.,

Proteomic and bioinformatics analyses of mouse liver microsomes. Int J Proteomics 2012, 2012, 832569. 7.

Mathias, R. A.; Chen, Y. S.; Kapp, E. A.; Greening, D. W.; Mathivanan, S.; Simpson,

R. J., Triton X-114 phase separation in the isolation and purification of mouse liver microsomal membrane proteins. Methods 2011, 54, (4), 396-406. 8.

Nelson, A. C.; Huang, W.; Moody, D. E., Variables in human liver microsome

preparation: impact on the kinetics of l-alpha-acetylmethadol (LAAM) n-demethylation and dextromethorphan O-demethylation. Drug Metab Dispos 2001, 29, (3), 319-25. 9.

Guengerich, F. P., Cytochrome p450 and chemical toxicology. Chem Res Toxicol

2008, 21, (1), 70-83. 10.

Hersman, E. M.; Bumpus, N. N., A targeted proteomics approach for profiling

murine cytochrome P450 expression. J Pharmacol Exp Ther 2014, 349, (2), 221-8. 11.

Wang, X.; Bian, Y.; Cheng, K.; Gu, L. F.; Ye, M.; Zou, H.; Sun, S. S.; He, J. X., A

large-scale protein phosphorylation analysis reveals novel phosphorylation motifs and phosphoregulatory networks in Arabidopsis. J Proteomics 2013, 78, 486-98. 12.

Pan, C.; Gnad, F.; Olsen, J. V.; Mann, M., Quantitative phosphoproteome analysis of 19

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

a mouse liver cell line reveals specificity of phosphatase inhibitors. Proteomics 2008, 8, (21), 4534-46. 13.

Kanaeva, I. P.; Petushkova, N. A.; Lisitsa, A. V.; Lokhov, P. G.; Zgoda, V. G.;

Karuzina, II; Archakov, A. I., Proteomic and biochemical analysis of the mouse liver microsomes. Toxicol In Vitro 2005, 19, (6), 805-12. 14.

Golizeh, M.; Sleno, L., Optimized proteomic analysis of rat liver microsomes using

dual enzyme digestion with 2D-LC-MS/MS. J Proteomics 2013, 82, 166-78. 15.

Clement, B.; Demesmaeker, M.; Linne, S., Microsomal Catalyzed N-Hydroxylation

of Guanabenz and Reduction of the N-Hydroxylated Metabolite: Characterization of the Two Reactions and Genotoxic Potential of Guanoxabenz1. Chemical research in toxicology 1996, 9, (4), 682-688. 16.

Lee, S. K.; Kim, J. H.; Seo, Y. M.; Kim, H. C. H.; Kang, M. J.; Jeong, H. G.; Lee, E.

S.; Jeong, T. C., In vitro characterization of the enzymes involved in the metabolism of 1furan-2-yl-3-pyridin-2-yl-propenone, an anti-inflammatory propenone compound. Archives of pharmacal research 2008, 31, (6), 764-770. 17.

Cox, J.; Mann, M., MaxQuant enables high peptide identification rates,

individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol 2008, 26, (12), 1367-72. 18.

Vizcaino, J. A.; Deutsch, E. W.; Wang, R.; Csordas, A.; Reisinger, F.; Rios, D.;

Dianes, J. A.; Sun, Z.; Farrah, T.; Bandeira, N.; Binz, P. A.; Xenarios, I.; Eisenacher, M.; Mayer, G.; Gatto, L.; Campos, A.; Chalkley, R. J.; Kraus, H. J.; Albar, J. P.; MartinezBartolome, S.; Apweiler, R.; Omenn, G. S.; Martens, L.; Jones, A. R.; Hermjakob, H., ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat Biotechnol 2014, 32, (3), 223-6. 19.

Jones, P.; Binns, D.; Chang, H. Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.;

Maslen, J.; Mitchell, A.; Nuka, G.; Pesseat, S.; Quinn, A. F.; Sangrador-Vegas, A.; Scheremetjew, M.; Yong, S. Y.; Lopez, R.; Hunter, S., InterProScan 5: genome-scale protein function classification. Bioinformatics 2014, 30, (9), 1236-40. 20.

Kanehisa, M.; Goto, S., KEGG: kyoto encyclopedia of genes and genomes. Nucleic

Acids Res 2000, 28, (1), 27-30. 21.

Horton, P.; Park, K. J.; Obayashi, T.; Fujita, N.; Harada, H.; Adams-Collier, C. J.;

Nakai, K., WoLF PSORT: protein localization predictor. Nucleic Acids Res 2007, 35, (Web Server issue), W585-7. 20

ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

22.

Klipper-Aurbach, Y.; Wasserman, M.; Braunspiegel-Weintrob, N.; Borstein, D.;

Peleg, S.; Assa, S.; Karp, M.; Benjamini, Y.; Hochberg, Y.; Laron, Z., Mathematical formulae for the prediction of the residual beta cell function during the first two years of disease in children and adolescents with insulin-dependent diabetes mellitus. Med Hypotheses 1995, 45, (5), 486-90. 23.

Chou, M. F.; Schwartz, D., Biological sequence motif discovery using motif-x. Curr

Protoc Bioinformatics 2011, Chapter 13, Unit 13 15-24. 24.

Redlich, G.; Zanger, U. M.; Riedmaier, S.; Bache, N.; Giessing, A. B.; Eisenacher,

M.; Stephan, C.; Meyer, H. E.; Jensen, O. N.; Marcus, K., Distinction between human cytochrome P450 (CYP) isoforms and identification of new phosphorylation sites by mass spectrometry. J Proteome Res 2008, 7, (11), 4678-88. 25.

Hornbeck, P. V.; Zhang, B.; Murray, B.; Kornhauser, J. M.; Latham, V.; Skrzypek, E.,

PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res 2015, 43, (Database issue), D512-20. 26.

Ding, S. J.; Qian, W. J.; Smith, R. D., Quantitative proteomic approaches for

studying phosphotyrosine signaling. Expert Rev Proteomics 2007, 4, (1), 13-23. 27.

Amanchy, R.; Periaswamy, B.; Mathivanan, S.; Reddy, R.; Tattikota, S. G.; Pandey,

A., A curated compendium of phosphorylation motifs. Nature biotechnology 2007, 25, (3), 285-286. 28.

Pearson, R. B.; Kemp, B. E., Protein kinase phosphorylation site sequences and

consensus specificity motifs: tabulations. Methods in enzymology 1991, 200, 62-81. 29.

Li, M.; Stefansson, B.; Wang, W.; Schaefer, E. M.; Brautigan, D. L., Phosphorylation

of the Pro-X-Thr-Pro site in phosphatase inhibitor-2 by cyclin-dependent protein kinase during M-phase of the cell cycle. Cellular signalling 2006, 18, (8), 1318-1326. 30.

Wang, Y.; Guan, S.; Acharya, P.; Liu, Y.; Thirumaran, R. K.; Brandman, R.; Schuetz,

E. G.; Burlingame, A. L.; Correia, M. A., Multisite phosphorylation of human liver cytochrome P450 3A4 enhances Its gp78- and CHIP-mediated ubiquitination: a pivotal role of its Ser-478 residue in the gp78-catalyzed reaction. Mol Cell Proteomics 2012, 11, (2), M111 010132. 31.

Ireland, B. S.; Niggemann, M.; Williams, D. B., In vitro assays of the functions of

calnexin and calreticulin, lectin chaperones of the endoplasmic reticulum. Methods Mol Biol 2006, 347, 331-42. 32.

Chevet, E.; Smirle, J.; Cameron, P. H.; Thomas, D. Y.; Bergeron, J. J., Calnexin 21

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

phosphorylation: linking cytoplasmic signalling to endoplasmic reticulum lumenal functions. Semin Cell Dev Biol 2010, 21, (5), 486-90. 33.

Trost, M.; English, L.; Lemieux, S.; Courcelles, M.; Desjardins, M.; Thibault, P., The

phagosomal proteome in interferon-gamma-activated macrophages. Immunity 2009, 30, (1), 143-54. 34.

Wong, H. N.; Ward, M. A.; Bell, A. W.; Chevet, E.; Bains, S.; Blackstock, W. P.;

Solari, R.; Thomas, D. Y.; Bergeron, J. J., Conserved in vivo phosphorylation of calnexin at casein kinase II sites as well as a protein kinase C/proline-directed kinase site. J Biol Chem 1998, 273, (27), 17227-35. 35.

Chevet, E.; Wong, H. N.; Gerber, D.; Cochet, C.; Fazel, A.; Cameron, P. H.; Gushue,

J. N.; Thomas, D. Y.; Bergeron, J. J., Phosphorylation by CK2 and MAPK enhances calnexin association with ribosomes. EMBO J 1999, 18, (13), 3655-66. 36.

Benyamini, P.; Webster, P.; Meyer, D. I., Knockdown of p180 eliminates the terminal

differentiation of a secretory cell line. Mol Biol Cell 2009, 20, (2), 732-44. 37.

Tsai, H. Y.; Yang, Y. F.; Wu, A. T.; Yang, C. J.; Liu, Y. P.; Jan, Y. H.; Lee, C. H.;

Hsiao, Y. W.; Yeh, C. T.; Shen, C. N.; Lu, P. J.; Huang, M. S.; Hsiao, M., Endoplasmic reticulum ribosome-binding protein 1 (RRBP1) overexpression is frequently found in lung cancer patients and alleviates intracellular stress-induced apoptosis through the enhancement of GRP78. Oncogene 2013, 32, (41), 4921-31. 38.

Krasnov, G. S.; Oparina, N.; Khankin, S. L.; Mashkova, T. D.; Ershov, A. N.;

Zatsepina, O. G.; Karpov, V. L.; Beresten, S. F., [Colorectal cancer 2D-proteomics: identification of altered protein expression]. Mol Biol (Mosk) 2009, 43, (2), 348-56. 39.

Telikicherla, D.; Marimuthu, A.; Kashyap, M. K.; Ramachandra, Y. L.; Mohan, S.;

Roa, J. C.; Maharudraiah, J.; Pandey, A., Overexpression of ribosome binding protein 1 (RRBP1) in breast cancer. Clin Proteomics 2012, 9, (1), 7. 40.

Liang, X.; Sun, S.; Zhang, X.; Wu, H.; Tao, W.; Liu, T.; Wei, W.; Geng, J.; Pang, D.,

Expression of ribosome-binding protein 1 correlates with shorter survival in Her-2 positive breast cancer. Cancer Sci 2015. 41.

Guengerich, F. P.; Cheng, Q., Orphans in the human cytochrome P450 superfamily:

approaches to discovering functions and relevance in pharmacology. Pharmacol Rev 2011, 63, (3), 684-99. 42.

Nelson, D. R.; Zeldin, D. C.; Hoffman, S. M.; Maltais, L. J.; Wain, H. M.; Nebert, D.

W., Comparison of cytochrome P450 (CYP) genes from the mouse and human genomes, 22

ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

including nomenclature recommendations for genes, pseudogenes and alternative-splice variants. Pharmacogenetics 2004, 14, (1), 1-18. 43.

Oesch-Bartlomowicz, B.; Oesch, F., Cytochrome-P450 phosphorylation as a

functional switch. Arch Biochem Biophys 2003, 409, (1), 228-34. 44.

Wang, Y.; Liao, M.; Hoe, N.; Acharya, P.; Deng, C.; Krutchinsky, A. N.; Correia, M.

A., A role for protein phosphorylation in cytochrome P450 3A4 ubiquitin-dependent proteasomal degradation. J Biol Chem 2009, 284, (9), 5671-84. 45.

Wang, Y.; Guan, S.; Acharya, P.; Koop, D. R.; Liu, Y.; Liao, M.; Burlingame, A. L.;

Correia, M. A., Ubiquitin-dependent proteasomal degradation of human liver cytochrome P450 2E1: identification of sites targeted for phosphorylation and ubiquitination. J Biol Chem 2011, 286, (11), 9443-56. 46.

Fairbank, M.; St-Pierre, P.; Nabi, I. R., The complex biology of autocrine motility

factor/phosphoglucose isomerase (AMF/PGI) and its receptor, the gp78/AMFR E3 ubiquitin ligase. Mol Biosyst 2009, 5, (8), 793-801. 47.

Wang, Y.; Kim, S. M.; Trnka, M. J.; Liu, Y.; Burlingame, A. L.; Correia, M. A.,

Human liver cytochrome P450 3A4 ubiquitination: molecular recognition by UBC7-gp78 autocrine motility factor receptor and UbcH5a-CHIP-Hsc70-Hsp40 E2-E3 ubiquitin ligase complexes. J Biol Chem 2015, 290, (6), 3308-32. 48.

Veglia, G.; Cembran, A., Role of conformational entropy in the activity and

regulation of the catalytic subunit of protein kinase A. FEBS J 2013, 280, (22), 5608-15. 49.

Redpath, N. T.; Proud, C. G., Cyclic AMP-dependent protein kinase phosphorylates

rabbit reticulocyte elongation factor-2 kinase and induces calcium-independent activity. Biochem J 1993, 293 ( Pt 1), 31-4. 50.

Browne, G. J.; Proud, C. G., Regulation of peptide-chain elongation in mammalian

cells. Eur J Biochem 2002, 269, (22), 5360-8. 51.

Kotani, S.; Tanaka, H.; Yasuda, H.; Todokoro, K., Regulation of APC activity by

phosphorylation and regulatory factors. J Cell Biol 1999, 146, (4), 791-800. 52.

Brushia, R. J.; Walsh, D. A., Phosphorylase kinase: the complexity of its regulation is

reflected in the complexity of its structure. Front Biosci 1999, 4, D618-41.

23

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

1

Page 24 of 34

Table 1. Selected list of phosphopeptides from proteins in processing pathways in endoplasmic reticulum from mouse liver microsomes Protein accession A2AVJ7

Protein names

Ribosome-binding protein 1

Gene names RRBP1

KEGG name P180

Position

Modified sequence

S167

_S(ph)AILEATPK_

S154

_VEPAVS(ph)SIVNSIQVLASK_

S155

_VEPAVSS(ph)IVNSIQVLASK_

O08795

Glucosidase 2 subunit beta

GLU2B

GIcII

S168

_S(ph)LEDQVETLR_

P07901

Heat shock protein (HSP) 90-alpha

HSP90A

Hap90

S263

_ESDDKPEIEDVGS(ph)DEEEEEKK_

P08113

Endoplasmin (Fragment)

F7C312

GRP94

S447

_GVVDSDDLPLNVS(ph)R_

P11499

Heat shock protein (HSP) 90-beta (Fragment)

E9Q3D6

Hap90

S255

_IEDVGS(ph)DEEDDSGKDKK_

P20029

78 kDa glucose-regulated protein

GRP78

BiP

T649

_LYGSGGPPPTGEEDT(ph)SEKDEL_

T474

_VYEGERPLT(ph)K_

S582

_AEEDEILNRS(ph)PR_

S553

_QKS(ph)DAEEDGVTGSQDEEDSKPK_

S563

_QKSDAEEDGVTGS(ph)QDEEDSKPK_

P35564 Q01853 Q3TDQ1

Calnexin Transitional endoplasmic reticulum ATPase Dolichyl-diphosphooligosaccharide--protein glycosyltransferase subunit STT3B

CALX

CNX

TERA

P97

S197

_EDEEES(ph)LNEVGYDDIGGCR_

STT3B

OSTs

S495

_ENPPVEDS(ph)SDEDDKRNPGNLYDK_

T375

_EELEQQT(ph)DGDCDEEDDDKDGEVPK_

S356

_SQHS(ph)SGNGNDFEMITK_ _GTENGVNGTVTSNGADS(ph)PR_

Q8BU14

Translocation protein SEC62

SEC62

Sce62/63

Q91V04

Translocating chain-associated membrane protein 1

TRAM1

TRAM

S365

Q922R8

Protein disulfide-isomerase A6

Q3TML0

PDIs

S428

_EPWDGKDGELPVEDDIDLS(ph)DVELDDLEKDEL_

S13

_PGPTPSGTNVGS(ph)SGR_

S17

_PGPTPSGTNVGSSGRS(ph)PSK_

S49

_TTS(ph)AGTGGMWR_

S268

_VEM(ox)GTSSQNDVDMSWIPQETLNQINKAS(ph)PR_

Q9CQS8 Q9CY50

Protein transport protein Sec61 subunit beta Translocon-associated protein subunit alpha

SC61B SSRA

Sec61 TRAP 24

ACS Paragon Plus Environment

Page 25 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Journal of Proteome Research

Q9DCF9

Translocon-associated protein subunit gamma

Q9JKR6

Hypoxia upregulated protein 1

Q9R049

E3 ubiquitin-protein ligase AMFR

Q9Z2B5

Eukaryotic translation initiation factor 2-alpha kinase 3

SSRG

TRAP

HYOU1

NEF

AMFR

gp78

E9QQ30

PERK

S11

_QQS(ph)EEDLLLQDFSR_

S567

_VESVFETLVEDS(ph)PEEESTLTK_

S552

_LEETLDFS(ph)EVELEPIEVEDFEAR_

S542

_VPLDLS(ph)PR_

S711

_EQIEVIAPS(ph)PER_

2

25

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

3

Page 26 of 34

Table 2. Selected list of phosphopeptides of murine cytochrome P450s in mouse liver microsomes

Human#

Mouse

4

Protein accession

Protein names

Position

Modified sequence

Score

Sequence

Protein

P56593

Cytochrome P450 2A12

S130

_RFS(ph)IATLR_

136.38

_RFS(ph)SIATLR_

CYP2A6

Q9WUD0

Cytochrome P450 2B10

S128*

_RFS(ph)LATMR_

53.567

_RFS(ph)VTTMR_

CYP2B6

P24456

Cytochrome P450 2D10

T275*

_NLT(ph)DAFLAEIEK_

211.26

-

-

Q8CIM7

Cytochrome P450 2D26

S138

_FS(ph)VSTLR_

152.23

-RFS(ph)VSTLR_

CYP2D6

Q05421

Cytochrome P450 2E1

S129

_RFS(ph)LSILR_

142.97

_RFS(ph)LTTLR_

CYP2E1

P33267

Cytochrome P450 2F2

S128

_RFS(ph)VQILR_

136.53

-

-

Q64464

Cytochrome P450 3A13

S134*

_ALLS(ph)PTFTSGR_

164.68

_SLLS(ph)PTFTSGK_

CYP3A4

Q64464

Cytochrome P450 3A13

S501

_VVSRDETVS(ph)DE_

79.089

-

-

Q64459

Cytochrome P450 3A11

T138*

_ALLSPTFT(ph)SGK_

162.48

-

-

Q64459

Cytochrome P450 3A11

S139

_ALLSPTFTS(ph)GK_

169.58

_RLRSLLSPTFTS(ph)GK_

CYP3A4

#

*, novel phosphorylation sites in mouse liver microsomes. , Described in literature as phosphorylation

5

26

ACS Paragon Plus Environment

24, 30

.

Page 27 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

6

Journal of Proteome Research

Table 3. The phosphorylation of ribosomal proteins by PKA and CK2 in mouse liver microsomes Kinase

PKA

CK2

7

Protein accession

Protein names

Score

Modified sequence

P27659

60S ribosomal protein L3

155.3

_HGS(ph)LGFLPR_

S13

Q9D8E6

60S ribosomal protein L4

55.8

_RNT(ph)ILRQAR_

T339

P47962

60S ribosomal protein L5

54.7

_MS(ph)LAQKKDR_

S272*

P47911

60S ribosomal protein L6

109.2

_AGS(ph)DAAASRPR_

S21

P14131

40S ribosomal protein S16

67.1

_T(ph)ATAVAHCK_

T18*

P62717

60S ribosomal protein L18a

86.9

_AHS(ph)IQIM(ox)K_

S123

P60867

40S ribosomal protein S20

62.5

_MPTKT(ph)LR_

T60*

P62830

60S ribosomal protein L23

88.1

_GGS(ph)SGAKFR_

S9*

Q8BP67

60S ribosomal protein L24

100.5

_AITGAS(ph)LADIMAK_

S86

P62852

40S ribosomal protein S25

117.5

_IRGS(ph)LAR_

S82*

P47915

60S ribosomal protein L29

57.4

_SQRYES(ph)LK_

S31*

50.2

_AVS(ph)ARAEAIK_

S66

P41105

60S ribosomal protein L28

122.96

_RAS(ph)AILR_

S115

Q9D1R9

60S ribosomal protein L34

99.5

_RLSYNT(ph)ASNK_

T15

62.861

_RLS(ph)YNTASNKTR_

S12

Q6ZWV7

60S ribosomal protein L35

72.1

_VVRKS(ph)IAR_

S53*

P83882

60S ribosomal protein L36a

69.5

_GKDS(ph)LYAQGK_

S32*

Q9D8E6

60S ribosomal protein L4

184.2

_ILKS(ph)PEIQR_

S295

P47963

60S ribosomal protein L13

155.6

_KGDS(ph)S(ph)AEELK_

S139

155.6

_KGDS(ph)S(ph)AEELK_

S140

P62270

40S ribosomal protein S18

173

_RAGELT(ph)EDEVER_

T60

*, novel phosphorylation site in ribosomal proteins. 27

ACS Paragon Plus Environment

Position

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

8

Figure legends

9

Figure 1. Detection of phosphoproteome in mouse liver microsomes. Electrophoretic

10

separation of microsomal proteins, as revealed by Coomassie staining (A) and

11

screening of phosphoproteome by Western blot (B). The individually prepared three

12

microsomal fractions (lanes 1–3, 30 µg of protein per lane) were separated on a 12%

13

Tris-glycine gel, followed by immunoblotting using anti-phospho-Ser, phospho-Thr,

14

and phospho-Tyr antibodies, respectively.

15

Figure 2. Scheme of the analytical workflow to identify phosphorylation sites in mouse

16

liver microsomes. Experimental workflow (A), area-proportional Venn diagram for

17

overlapping identified unique phosphoproteins (B), and frequency of phosphoserine,

18

phosphothreonine, and phosphotyrosine (C). The microsomes had been prepared from

19

three independent liver microsome preparations.

20

Figure 3. Motif analysis of Ser and Thr phosphorylation sites in mouse liver microsomes.

21

Acidic (A), basophilic (B) and proline-directed (D) motifs. Sites that were within six

22

residues of protein termini were not used in the motif logo. Ser and Thr

23

phosphorylation sites were identified in mouse liver microsomes by Motif X using the

24

13-mer sequences.

25

Figure 4. Characterization of mouse liver microsomal phosphoproteome. Phosphorylated

26

protein functional enrichment by biological processes, cellular component, and

27

molecular function (A) and phosphorylated protein functional enrichment by KEGG

28

pathway (B).

29

Figure 5. Immunoblots to detect phosphorylation in mouse liver microsomes after

30

reaction with PKA and CKII. Lane 1; reaction without PKA, lane 2: reaction with 28

ACS Paragon Plus Environment

Page 28 of 34

Page 29 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

31

Journal of Proteome Research

PKA, lane 3: reaction without CKII; lane 4: reaction with CKII.

29

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

32

Figure 1.

33

30

ACS Paragon Plus Environment

Page 30 of 34

Page 31 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

34

Journal of Proteome Research

Figure 2

35 36

31

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

37

Figure 3

38 39

32

ACS Paragon Plus Environment

Page 32 of 34

Page 33 of 34

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

40

Journal of Proteome Research

Figure 4

(A)

mRNA processing

4.84

Organic substance transport

3.91

Translation

Biological Process

3.74

Establishment of localization in cell

3.58

mRNA metabolic process

3.31

RNA processing

2.17

Cellular localization

2.12

RNA splicing

2.08

Ribosome

Cellular Component

5.98

Spliceosomal complex

4.32

Organelle membrane

4.18

Ribonucleoprotein complex

4.16

Endoplasmic reticulum part

3.85

Endoplasmic reticulum membrane

3.69 2.98

Endoplasmic reticulum Nuclear outer membrane-endoplasmic reticulum membrane network Endomembrane system

Molecular Function

2.77 2.07

Oxidoreductase activity

10.8

RNA binding

3.75

0

2

4

6

8

10

Fold enrichment

(B) (mmu03010) Ribosome

4,40

(mmu04141) Protein processing in endoplasmic reticulum

4.29

(mmu00140) Steroid hormone biosynthesis

3.87

(mmu03013) RNA transport

3.79

(mmu03040) Spliceosome

3.45

(mmu04910) Insulin signaling pathway

2.27

0

1

2

3

4

Fold enrichment

41 42

33

ACS Paragon Plus Environment

5

12

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

43

Figure 5

44

34

ACS Paragon Plus Environment

Page 34 of 34