Quantitative Phosphoproteomics Dissection of Seven-transmembrane Receptor Signaling Using Full and Biased Agonists*

Seven-transmembrane receptors (7TMRs) signal through the well described heterotrimeric G proteins but can also activate G protein-independent signaling pathways of which the impact and complexity are less understood. The angiotensin II type 1 receptor (AT1R) is a prototypical 7TMR and an important drug target in cardiovascular diseases. “Biased agonists” with intrinsic “functional selectivity” that simultaneously blocks Gαq protein activity and activates G protein-independent pathways of the AT1R confer important perspectives in treatment of cardiovascular diseases. In this study, we performed a global quantitative phosphoproteomics analysis of the AT1R signaling network. We analyzed ligand-stimulated SILAC (stable isotope labeling by amino acids in cell culture) cells by high resolution (LTQ-Orbitrap) MS and compared the phosphoproteomes of the AT1R agonist angiotensin II and the biased agonist [Sar1,Ile4,Ile8]angiotensin II (SII angiotensin II), which only activates the Gαq protein-independent signaling. We quantified more than 10,000 phosphorylation sites of which 1183 were regulated by angiotensin II or its analogue SII angiotensin II. 36% of the AT1R-regulated phosphorylations were regulated by SII angiotensin II. Analysis of phosphorylation site patterns showed a striking distinction between protein kinases activated by Gαq protein-dependent and -independent mechanisms, and we now place protein kinase D as a key protein involved in both Gαq-dependent and -independent AT1R signaling. This study provides substantial novel insight into angiotensin II signal transduction and is the first study dissecting the differences between a full agonist and a biased agonist from a 7TMR on a systems-wide scale. Importantly, it reveals a previously unappreciated diversity and quantity of Gαq protein-independent signaling and uncovers novel signaling pathways. We foresee that the amount and diversity of G protein-independent signaling may be more pronounced than previously recognized for other 7TMRs as well. Quantitative mass spectrometry is a promising tool for evaluation of the signaling properties of biased agonists to other receptors in the future.

Seven-transmembrane receptors (7TMRs) signal through the well described heterotrimeric G proteins but can also activate G protein-independent signaling pathways of which the impact and complexity are less understood. The angiotensin II type 1 receptor (AT 1 R) is a prototypical 7TMR and an important drug target in cardiovascular diseases. "Biased agonists" with intrinsic "functional selectivity" that simultaneously blocks G␣ q protein activity and activates G protein-independent pathways of the AT 1 R confer important perspectives in treatment of cardiovascular diseases. In this study, we performed a global quantitative phosphoproteomics analysis of the AT 1 R signaling network. We analyzed ligand-stimulated SILAC (stable isotope labeling by amino acids in cell culture) cells by high resolution (LTQ-Orbitrap) MS and compared the phosphoproteomes of the AT 1 R agonist angiotensin II and the biased agonist [Sar 1 ,Ile 4 ,Ile 8 ]angiotensin II (SII angiotensin II), which only activates the G␣ q protein-independent signaling. We quantified more than 10,000 phosphorylation sites of which 1183 were regulated by angiotensin II or its analogue SII angiotensin II. 36% of the AT 1 R-regulated phosphorylations were regulated by SII angiotensin II. Analysis of phosphorylation site patterns showed a striking distinction between protein kinases activated by G␣ q protein-dependent and -independent mechanisms, and we now place protein kinase D as a key protein involved in both G␣ q -dependent and -independent AT 1 R signaling. This study provides substantial novel insight into angiotensin II signal transduction and is the first study dissecting the differences between a full agonist and a biased agonist from a 7TMR on a systems-wide scale. Importantly, it reveals a previously unappreciated diversity and quantity of G␣ q protein-independent signaling and uncovers novel signaling pathways. We foresee that the amount and diversity of G protein-independent signaling may be more pronounced than previously recognized for other 7TMRs as well. Quantitative mass spectrometry is a promising tool for evaluation of the signaling properties of biased agonists to other receptors in the future.

Molecular & Cellular Proteomics 9:1540 -1553, 2010.
Seven-transmembrane receptors (7TMRs) 1 constitute the largest family of plasma membrane receptors. These receptors regulate a variety of biological processes and are currently the most prominent therapeutic targets, highlighted by an abundance of clinically available drugs (1)(2)(3)(4)(5). Their cellular responses were until recently thought to depend primarily on G protein activation and generation of second messengers such as inositol 1,4,5-trisphosphate, diacylglycerol, and cAMP, reflected by the term G protein-coupled receptors. However, the 7TMR structure also confers activation of G protein-independent signaling initiated by direct interaction with ␤-arrestin or tyrosine kinases (1-3, 6, 7). The discovery of alternative signaling pathways to the G protein-initiated pathways has also led to the discovery of regulatory peptides that can have agonistic effects on one pathway while simultaneously antagonizing another. Such pharmacological separation of signaling pathways has an enormous clinical potential, which is described by pharmacological concepts such as "functional selectivity" or "biased agonism" (3,8,9). A number of antagonists have only been evaluated by G protein-dependent readouts and may therefore turn out to be biased agonists. An example of such biased agonism is the retrospective discovery that the widely used inverse ␤ 2 -adrenergic receptor agonist carvedilol, which has been superior to other ␤-blockers in the prevention of heart failure, functions as a biased agonist. The reason for the increased survival of patients treated with carvedilol is unclear, but this drug acts as a biased agonist by activating the Erk1/2 mitogen-activated protein kinases while inhibiting G␣ s -mediated effects (10). This could be the explanation behind its clinical success because the G protein-independent signaling could confer the effects that improve survival (10). Similar mechanisms are likely to underlie the molecular mechanism behind other successful 7TMR drugs. Comparable observations that G protein-independent signaling may provide beneficial effects have been made in model systems for angiotensin II (Ang II) signaling via the Ang II type 1 receptor (AT 1 R). The AT 1 R is commonly used as a model system to study G␣ q -dependent versus -independent signaling because of the availability of excellent and well defined biological tools for this receptor.
The AT 1 R serves as a key regulator of cardiovascular physiology, maintaining salt and fluid homeostasis and blood pressure (11,12). The AT 1 R is also involved in a number of medical conditions including endothelial dysfunction and atheroma, cardiac hypertrophy and failure, atrial fibrillation, nephropathy, insulin resistance, and cancer (13) and is therefore a prominent drug target in cardiovascular diseases, reflected by the common use of AT 1 R blockers such as losartan and angiotensin-converting enzyme inhibitors, for example ramipril. AT 1 R signaling is particularly interesting with respect to functional selectivity because selective activation of G␣ q proteinindependent signaling has proven relevant in primary cells, in isolated organs, and in vivo (3,14,15). AT 1 R signaling pathways encompass a complex network of signaling molecules including heterotrimeric G proteins, protein kinases, and scaffold proteins (11,12). The biased agonist [Sar 1 ,Ile 4 ,Ile 8 ]Ang II (SII Ang II) is a well established AT 1 R ligand that provides the opportunity to study G␣ q -independent effects (16). SII Ang II selectively inhibits the G␣ q protein signaling branch while still activating the remaining network that is activated from the AT 1 R such as tyrosine kinase signaling (6, 7), ␤-arrestin signaling (16), and the possibly minor contribution from other trimeric G proteins such as G 12/13 , G i , and G o (17)(18)(19) that collectively is called G␣ q -independent signaling (Fig. 1A). Activation of Erk1/2 has been thoroughly investigated with respect to G␣ q protein-dependent versus ␤-arrestin-dependent signaling (20 -22). These studies have established that Erk1/2 is activated by both pathways and furthermore that there is a striking pathway-dependent difference in the function and localization of activated Erk1/2 (20). G␣ q -dependent activation of Erk1/2 results in nuclear trans-location and gene transcription, whereas G␣ q -independent activation leads to a tight association with ␤-arrestin, which results in cytosolic sequestration of Erk1/2 (20 -22). This difference in Erk1/2 localization has been related to several important phenotypic differences. The ␤-arrestin part of G␣ qindependent signaling has been linked to cellular survival and renewal as well as regulation of cell migration (23,24). Conversely, the G␣ q protein-mediated component of AT 1 R signaling has been associated with cellular death and fibrosis leading to cardiac hypertrophy and progression to heart failure (3). These divergent roles of the two signaling pathways indicate that full therapeutic blockade of the AT 1 R may have drawbacks in some situations due to blocking desirable physiological functions. Biased agonism provides a novel pharmacological potential as the use of biased agonists may allow selective blocking of undesired effects while maintaining desired effects (2,3).
Global and site-specific analysis of in vivo phosphorylation sites by quantitative mass spectrometry has emerged as the method of choice to investigate cell signaling pathways in an unbiased fashion (25). In this study, we characterized the AT 1 R signaling network by high resolution quantitative phosphoproteomics and thereby delineated the differences between the full agonist Ang II and the biased agonist SII Ang II. To characterize the AT 1 R signaling network, we used stable isotope labeling by amino acids in cell culture (SILAC) in combination with phosphopeptide enrichment and high performance mass spectrometry (26). This approach has become increasingly powerful in the analysis of signaling pathways because it allows for a high throughput and quantitative mapping of changes in cellular protein phosphorylation on a systems-wide scale. The majority of the identified proteins in our study have not previously been associated with AT 1 R signaling, and we found G␣ q protein-independent signaling to play an unexpected large role in AT 1 R signaling.

PKD Kinase Assays
AT 1 R-HEK cells were stimulated as described above and lysed on ice for 15min in lysis buffer (0.5% Triton X-100, 150 mM NaCl, 50 mM Tris-HCl, pH 7.4, 1 mM Na 3 VO 4 , 5 mM EDTA, 50 mM NaF, 10 nM calyculin A, 10 M leupeptin, 5 M pepstatin, 1 g/ml aprotinin). Lysates were cleared for 15 min at 13,000 rpm. Immunoprecipitation was performed by addition of 1.5 g of PKD antibody (Santa Cruz Biotechnology, Inc.) for 1 h followed by 1-h addition of Protein G-agarose beads. The beads were washed three times with lysis buffer and twice in kinase buffer (45 mM Tris-HCl, pH 7.4, 15 mM MgCl 2 , 1.5 mM dithiothreitol) and resuspended in 20 l of kinase buffer. The kinase reaction was initiated by addition of 10 l of reaction buffer (450 M ATP, 0.25 g/l syntide 2 (PLARTLSVA-GLPGKK), 0.2 Ci of [␥-32 P]ATP), and phosphorylation was allowed under vigorous shaking for 10 min at 30°C. 20 l of the kinase reaction were spotted onto phosphocellulose paper (Whatman p81), and the paper was washed with 150 mM orthophosphoric acid. Phosphorylation of syntide 2 was determined by autoradiography on a STORM TM phosphorimaging system using ImageQuant TM software (GE Healthcare).

In-solution Digestion
The acetone-precipitated lysate was resolubilized in an 8 M urea buffer (6 M urea ϩ 2 M thiourea in 10 mM HEPES, pH 8). The soluble proteins were reduced for 30 min at room temperature with 1 mM DTT and alkylated for 15 min by 5.5 mM iodoacetamide. Endoproteinase Lys-C (Wako) was added at 1:100 (w/w), and the lysates were digested for 4 h at room temperature. The resulting peptide mixtures were diluted 4-fold with deionized water to achieve a final urea concentration below 2 M. Trypsin (modified sequencing grade, Promega) was added at 1:100 (w/w), and the sample was digested overnight. Trypsin and Lys-C activity was quenched by acidification of the reaction mixtures with a 10% TFA solution to pH ϳ2. The peptide mixture was desalted and concentrated on a C 18 SepPak cartridge (Waters) and eluted with 2 ϫ 1 ml of 60% ACN in 0.3% TFA.
Proteins from the insoluble chromatin pellet were extracted by repeated incubations in a mixture of 8 M urea and Benzonase (Merck). The extracted proteins were treated, digested, and purified by the same approach as the acetone precipitate above.

Strong Cation Exchange Chromatography
Peptide fractionation by strong cation exchange (SCX) chromatography was performed to reduce the complexity of the peptide mixture.
The in-solution digests were loaded onto an SCX column (1-ml Resource S, GE Healthcare) using an Ä kta FPLC system (GE Healthcare). Peptides were separated by a linear gradient for 30 min from 100% buffer A (5 mM KH 2 PO 4 , pH 2.7, 30% ACN) to 30% buffer B (350 mM KCl, 5 mM KH 2 PO 4 , pH 2.7, 30% ACN) with a flow rate of 1 ml/min. Fractions were collected every minute using a fraction collector. In total, 16 SCX fractions including the column flow-through were collected for further analysis.

Phosphopeptide Enrichment
Phosphopeptides were enriched using Titansphere chromatography as described (26). Briefly, phosphopeptides were enriched on titanium dioxide beads (10-m Titansphere, GL Sciences) precoated with 2,5-dihydroxybenzoic acid. 5 g of coated-beads were added to each SCX fraction and incubated with rotation for 15 min. The beads were washed once with 200 l of 30% ACN in 3% TFA and once with 200 l of 60% ACN in 0.3% TFA and transferred in 200 l of 75% acetonitrile in 0.3% TFA on top of a reversed phase C 8 StageTip. The bound phosphopeptides were eluted directly into a 96-well plate by 2 volumes of 20 l of 25% ACN in 15% NH 4 OH, water (pH Ͼ 11). The eluates were immediately dried in a SpeedVac at 60°C.

LC-MS/MS
The dried phosphopeptide mixtures were acidified with 5% acetonitrile in 0.3% TFA to an end volume of 8 l and analyzed by on-line nanoflow LC-MS/MS as described previously (26) with a few modifications. Briefly, all nano-LC-MS/MS experiments were performed on an EASY-nLC TM system (Proxeon Biosystems, Odense, Denmark) connected to the LTQ-Orbitrap XL (Thermo Electron, Bremen, Germany) through a nanoelectrospray ion source.
The tryptic peptide and phosphopeptide mixtures were separated in a 15-cm analytical column (75-m inner diameter) in housepacked with 3-m C 18 beads (Reprosil-AQ Pur, Dr. Maisch) with a 90-min gradient from 5 to 30% acetonitrile in 0.5% acetic acid. The effluent from the HPLC was directly electrosprayed into the mass spectrometer.
The MS instrument was operated in data-dependent mode to automatically switch between full-scan MS and MS/MS acquisition. Survey full-scan MS spectra (from m/z 300 to 2000) were acquired in the orbitrap with resolution R ϭ 60,000 at m/z 400 (after accumulation to a "target value" of 1e6 in the linear ion trap). The 10 most intense peptide ions with charge states Ն2 were sequentially isolated to a target value of 5e3 and fragmented in the linear ion trap by multistage activation (MSA; or pseudo-MS3) (26,29) All fragment ion spectra were recorded with the LTQ detectors. For all measurements with the orbitrap detector, a lock mass ion from ambient air (m/z 445.120024) was used as an internal calibrant as described (30). Typical mass spectrometric conditions were as follows: spray voltage, 2.4 kV; no sheath and auxiliary gas flow; heated capillary temperature, 150°C; normalized collision energy, 35% for MSA in the LTQ with wideband activation enabled. The ion selection threshold was 500 counts for MS2, and the maximum allowed ion time was 500 ms for full scans and 150 ms for MS/MS. An activation q of 0.25 and activation time of 30 ms were used. For MSA, neutral losses of m/z 32.66, 48.99, and 97.97 were activated.

Raw MS Data Analysis
Peptide Identification by Mascot and MaxQuant-Raw Orbitrap full-scan MS and ion trap MSA spectra were processed by MaxQuant as described (31). In brief, all identified SILAC triplets were quantified, accurate precursor masses were determined using the entire LC elution profiles, and MS/MS spectra were filtered to contain at most six peaks per 100-Da interval and merged into peak list files (*.msm).
Peptides and proteins were identified by Mascot (Matrix Science, London, UK) via automated database matching of all tandem mass spectra against an in-house curated concatenated target/decoy database, a forward and reversed version of the human International Protein Index sequence database (version 3.37, 138,632 forward and reversed protein sequences from the European Bioinformatics Institute) supplemented with common contaminants such as human keratins, bovine serum proteins, and porcine trypsin. Tandem mass spectra were initially matched with a mass tolerance of 7 ppm on precursor masses and 0.5 Da for fragment ions and strict trypsin specificity, allowing for up to three missed tryptic cleavage sites. N-Pyroglutamine (Gln ϩ17.026549 Da), oxidized methionine (ϩ15.994915 Da), acetylation of protein (N terminus ϩ42.010565 Da), and phosphorylation of serine, threonine, and tyrosine (Ser/Thr/Tyr ϩ79.966331 Da) were searched as variable modifications.
Peptide Filtering, Quantitation, and Phosphosite Localization by MaxQuant-The resulting Mascot result files (*.dat) were loaded into the MaxQuant software suite version 1.0.12.35 for further processing. By automatically filtering on peptide length, mass error precision estimates, and Mascot score of all forward and reversed peptide identifications, we fixed the estimated false discovery rate (FDR) of all peptide and protein identifications at 1% using a Bayesian approach as detailed previously (32). The FDR is a framework to address the multiple hypothesis testing problem by restricting the number of false positives to a fixed percentage as opposed to e.g. a fixed Mascot score threshold. The database setup and FDR estimation are equivalent to those in Ref. 33; additional variable modifications in the search algorithm were used to increase the search space. Ratio determination by MaxQuant was done by robust linear regression of the two-dimensional centroid values of the intensity/time elution profiles for all isotope distributions that were present in the different SILAC states to determine the overall ratio difference. To pinpoint the actual phosphorylated amino acid residue(s) within all identified phosphopeptide sequences in an unbiased manner, MaxQuant calculated the localization probabilities of all putative serine, threonine, and tyrosine phosphorylation sites using the post-translation modification score algorithm as described (26). For final ratio calculation, MaxQuant calculates the median of all peptide ratios associated with a specific phosphorylation site including variants such as oxidation of methionine and missed cleavages. As a last filter, only sites with a localization probability Ͼ75% and a posttranslation modification delta score of at least 5 were included in the analysis. The extensively annotated spectra showing multiple water and phosphopeptide losses (34) explain more than 70% of the fragment ions (supplemental Fig. X). All spectra can be found in supplemental Data IX.
The determined regulation of a specific phosphorylation site identified in this study is most likely not related to the Mascot score as ratio information is not included in the identification step. To verify this, we plotted the Mascot score distributions and localization probability distributions from the phosphopeptides we report as regulated and compared them with those not reported as regulated (supplemental Fig. XI). As illustrated, the distributions overlapped almost completely, indicating that there indeed is no relation between identification and showing a significant ratio.
Cutoffs, Phosphorylation Site Annotation, and Analysis-The two biological replicates and the two time points were considered as four independent experiments. Phosphopeptides were required to have at least two independent ratios above the cutoff to be considered reproducibly regulated. This method ensures that peptides identified across all four conditions are also more accurately measured in the individual experiments because of a higher peptide count. The cutoff was set at 3 times the median standard deviation of the log-transformed ratios of the phosphorylated peptides. The ratios of the nonphosphorylated peptides were used to estimate the proportion of false positives as shown in Fig. 2B. The final cutoff was set at a ratio of 1.5 observed in at least two experiments, which provides an estimated 0.3% false positive rate.
Regulated phosphosites were annotated with UniProt accession numbers and compared with UniProt, PHOSIDA, and Phospho.ELM databases of known phosphosites (supplemental Tables I and II). We evaluated whether a phosphorylated protein was previously observed to be regulated by Ang II by searching for reports in PubMed (supplemental Tables I and II). Functional grouping of regulated proteins was done manually. All regulated sites are shown in supplemental Table IV, and all the detected phosphosites are shown in supplemental Table V. All identified proteins are listed in supplemental Table VII, and all modified peptides are listed in supplemental Table VIII.
Phosphosite Sequence Motif Logo Plots-To identify enriched sequence motifs in the phosphorylation site data set, we made use of an in-house algorithm that iteratively tests for position-specific overrepresentation of any amino acid groups in a prealigned list of the regulated sequences compared with the average occurrence of the amino acid group in a list of sequences surrounding non-regulated phosphosites. In each round of iteration, the most significant amino acid group is excluded in a position-specific fashion from both lists (35). Statistically significant over-representation is calculated using the R (a programming language and software environment for statistical computing and graphics) implementation of Fisher's exact test. Grouping of amino acids was done on the basis of related chemical properties (acidic, basic, aromatic, aliphatic, hydrophilic, amide, polar, and cyclic). The Perl package of WebLogo was used internally in the algorithm to visualize the enriched sequence motifs as logo plots. For the kinase motif analysis, we used a sequence window of Ϯ6 amino acids surrounding the central phosphorylated serine, threonine, or tyrosine residue. All up-regulated phosphorylation site sequences were compared with non-regulated phosphorylation site sequences from the results from either Ang II or SII Ang II, respectively. We considered a motif significant if it fulfilled our conservative cutoff of p Ͻ 0.001 on the Bonferroni correction-adjusted p values.

Dissecting AT 1 R Signaling by Mass Spectrometry-To
identify novel proteins in the AT 1 R signaling network and determine the weight and nature of G␣ q -dependent and -independent signaling, we used SILAC combined with phosphopeptide enrichment and high resolution mass spectrometry. This approach simultaneously enables high throughput and quantitative analysis of changes in multiple protein phosphorylation sites (26). We compared AT 1 R-HEK cells incubated with and without Ang II and SII Ang II. Ang II stimulation fully activates the AT 1 R, whereas SII Ang II stimulation only activates the ␤-arrestin-dependent and other G␣ q protein-independent signaling events (16,20). As can be seen from Fig. 1B, SII Ang II was unable to induce inositol accumulation (a readout for G␣ q activity) but could still activate Erk1/2, which has previously been shown to be ␤-arrestin-dependent (20) (Fig. 1B).
Two pools of triple SILAC-labeled cells were stimulated with vehicle, Ang II, and SII Ang II for 3 or 15 min, and these experiments were repeated as biological, independent repli-cates. Lysates from vehicle-, Ang II-, and SII Ang II-treated cells were mixed, the proteins were enzymatically digested with trypsin and separated into 16 fractions, and phosphopeptides were enriched on TiO 2 columns (26) (Fig. 1B). Each fraction was independently analyzed (the last four fractions were combined) by on-line reversed phase nanoflow LC-MS on a hybrid linear ion trap orbitrap mass spectrometer (LTQ-Orbitrap), and phosphopeptides were fragmented by multistage activation ("pseudo-MS3") (26). All SILAC-encoded peptides eluted as characteristic triplets, and their intensities, reflecting the relative phosphorylation state of the peptide, were measured by high resolution orbitrap mass spectrometry (Fig. 1B). To validate the reliability of the method, the relative phosphorylation intensities from mass spectrometry analysis were compared with standard Western blotting analysis for selected proteins. As exemplified with Erk1/2, the changes of the peak intensities were comparable with what was observed using phosphosite-specific antibodies (Fig. 1C).
AT 1 Receptor Signaling Affects Multiple Cellular Processes: 36% Is Regulated by G␣ q -independent Signaling-With this approach, we identified 10,965 unique phosphosites with more than 99% certainty of correct peptide identification and 75% certainty of phosphosite localization. These peptides were found on 2661 different proteins, and the distribution of Tyr(P), Thr(P), and Ser(P) was 1.2, 12.0, and 86.8%, respectively, which is comparable with previous observations (26) (Fig. 2, A and C, and supplemental Table III). The distribution of all up-regulated Tyr(P), Thr(P), and Ser(P) sites was 0.6, 11.1, and 88.3%, respectively (Fig. 2C and supplemental Tables I and II), Ang II SII Ang II  1. Experimental setup. A, illustration showing the principle signaling pathways from the AT 1 R. The AT 1 R associates with heterotrimeric G proteins but can also signal through G protein-independent mechanisms, which include direct binding to tyrosine kinases or signaling complexes scaffolded by ␤-arrestin (␤-ARR). Regulation of Erk1/2 is well described in this system and characterized by differential localization upon G protein-dependent and -independent activation. B, SII Ang II mediates G␣ q -independent Erk1/2 activation in AT1R-HEK cells. Upper panel, dose-response curves for agonistinduced inositol phosphate (IP) accumulations are shown as a mean of three independent experiments performed in triplicate. Data are depicted in cpm/well (ϮS.E.). The concentration of SII Ang II was 187 times the concentration of Ang II due to equally lower binding affinity (16). Lower panel, Western blot analysis of Ang II-and SII Ang IIinduced Erk1/2 activation as a function of time. Data are presented as the mean and error bars reflect Ϯ S.E. of four experiments. C, the workflow of the quantitative phosphoproteomics approach by SILAC. Three populations of AT 1 R-HEK cells were labeled with normal and stable isotope-substituted arginine and lysine, creating proteins and peptides distinguishable by mass. The medium and heavy labeled cell pools were stimulated with 100 nM Ang II or 18.7 M SII Ang II, respectively (for 3 and 15 min), the cells were lysed, and lysates were mixed 1:1:1. Proteins were digested by trypsin and endoproteinase Lys-C and fractionated by SCX, and phosphopeptides were enriched on TiO 2 columns. The resulting high resolution orbitrap full-scan MS spectra reveal SILAC triplets of each peptide, allowing quantification of hormone-induced changes in the phosphorylation level between the three groups. The 3-and 15-min samples were treated as individual experiments. A phospho-Erk (pERK) spectrum is shown to illustrate the high resolution MS spectrum. For comparison, a Western blot of phospho-Erk1/2 on the same lysate is illustrated in the top right corner. Ctrl, control; (ph), phosphorylation. and a similar distribution was seen for the Ang II-and SII Ang II-regulated peptides. The number of phosphorylation sites identified is relatively high compared with other proteomics studies, and therefore the analysis uncovered many novel phosphorylation sites (36 -38).
In the following bioinformatics study of our data set, we only analyzed the phosphorylation sites that were up-regulated as our experience is that down-regulation is less reproducible and can be due to both activation of phosphatases and general degradation of protein, which can change the absolute protein concentration much faster than synthesis. The cutoff was set at 3 times the median standard deviation of the log-transformed ratios of the phosphorylated peptides, equal to 1.5-fold. The ratios of all non-phosphorylated peptides detected in the study, which were expected to be unaffected by the short term treatment (ratio of 1:1:1), were used to estimate the proportion of false positives (0.3%) in our regulated set of phosphopeptides as shown on Fig. 2B.
1183 phosphosites were found to be regulated; 36% of these or 427 phosphosites were regulated by SII Ang II. This group represents phosphosites that are affected by G␣ qindependent signaling. A smaller group of 58 sites was found to be regulated by SII Ang II but not Ang II; many of these sites are probably also regulated by Ang II but did not reach the threshold.
In general, we found SII Ang II to be a weaker inducer of common phosphorylation sites than Ang II alone. This is not surprising as we and others have previously shown that SII Ang II is only a partial agonist on ␤-arrestin-dependent signaling as reflected by a weaker recruitment of ␤-arrestin (8,39). In this study, SII Ang II was also a weaker inducer than Ang II of almost all phosphorylations (generally 20 -50% of the Ang II response). Despite this, we used a 1.5-fold ratio increase as the cutoff for both Ang II and SII Ang II. However, the impact of AT 1 Rmediated G␣ q -independent signaling may be underestimated because phosphosites that are only regulated 1.5-3-fold with Ang II might also be regulated by SII Ang II but fall under the cutoff. Using these parameters, we found 1183 up-regulated phosphorylation sites on 527 phosphoproteins. For some phosphosites, SII Ang II stimulation yielded 50% of Ang II stimulation and on some sites only 20%. This difference could indicate that some of the Ang II-regulated phosphosites are affected by both G␣ q -dependent and -independent signaling, whereas others are strictly G␣ q -independent. An example of this is the well described Erk1/2 kinases discussed above in which activation loop phosphorylation sites were regulated about 20- (2nd highest ratio )  2. Distribution of phosphopeptides. A, we identified 10,965 unique phosphosites. 1183 of these peptides were reproducibly regulated 1.5-fold by Ang II or SII Ang II after 3-and/or 15-min stimulation. 36% of the AT 1 R-regulated phosphorylation sites were regulated by SII Ang II. B, the cutoff was determined as three standard deviations from the mean, equal to 1.5-fold. Comparing the distribution of non-phosphorylated peptides with that of phosphorylated peptides allowed us to determine a false positive identification rate of 0.3%. C, the distribution of all identified serine, threonine, and tyrosine phosphorylation sites compared with the distribution of regulated sites (REG). D, the 1183 phosphorylation sites were manually grouped in functional groups according to the main known cellular function of the protein on which they reside. The outer circle indicates the fraction of phosphorylation sites regulated by Ang II only (gray) and by SII Ang II (black). fold with Ang II but only 4-fold with SII Ang II. Erk1/2 are in fact regulated by both G␣ q -dependent and -independent pathways. Thus, SII Ang II only affects the G␣ q -independent pathway and moreover is a partial agonist on this pathway, and thus SII Ang II only yielded 20% of the regulation observed with Ang II. For other kinases such as CAMK2D, SII Ang II up-regulated phosphorylation sites to around 50% of Ang II induction, indicating that G␣ q -independent signaling is dominant for regulation of this protein.
The present study adds substantially to uncovering the Ang II signaling network as more than 80% of these phosphoproteins have not previously been linked to Ang II signaling (supplemental Tables I and II). Moreover, the analysis uncovered many novel phosphorylation sites as 32% of the regulated phosphorylation sites are not reported in public databases of known phosphosites (UniProt, PHOSIDA, and Phospho.ELM) (supplemental Tables I and II). We identified 3866 phosphorylation sites that are not previously reported in UniProt. Detailed information on all identified phosphosites is provided in supplemental Table III. As depicted in Fig. 2D, the phosphoproteins regulated by Ang II affect a variety of different cellular functions. We manually assigned the phosphoproteins into functional groups based on their associated gene ontology terms, and the inner circle of Fig. 2D shows the overall distribution of phosphosites, whereas the outer ring of the circle illustrates the fraction of phosphosites regulated by Ang II only (gray) and phosphosites regulated by SII Ang II (black). These data demonstrate a previously unappreciated versatility of G␣ q -independent signaling and participation in Ang II-regulated cellular processes. According to this observation, the description "G protein-coupled receptors" can be misleading, whereas "7TM receptors" may be more accurate.
G␣ q -dependent and -independent Signaling Activates Different Kinases-The protein kinases are the backbone in cellular signaling pathways by transmitting signals from the plasma membrane to intracellular proteins through site-specific phosphorylation cascades. We identified 100 reproducibly up-regulated phosphorylation sites positioned on 42 kinases (Table I and see Fig. 5). Of these kinases, 23 have not previously been described in Ang II signaling. SII Ang II stimulation increased phosphorylation levels of 19 kinases. Although including mostly well known Ang II-regulated kinases, 14 of these kinases were not known to be regulated independently of G␣ q protein. To evaluate the reproducibility of phosphorylation levels of phosphosites found in both biological replicates, we performed a Pearson's analysis of correlation. Regulation of phosphorylation sites on kinases was highly reproducible, showing a correlation of r 2 ϭ 0.90 between the two biological replicates performed (Fig. 3A).
Comparing linear sequence motifs surrounding the regulated phosphorylation sites revealed differences in protein kinases activated by Ang II only (G protein-dependent) and Ang II plus SII Ang II (G protein-independent). Both Ang II signaling through G␣ q -dependent and -independent pathways preferentially stimulated phosphorylation of a basophilic AGC/CAM kinase motif ((R/K)XXp(S/T) where p represents phosphorylation) (Fig. 3B, upper panel). Indeed, among the 19 kinases phosphorylated by SII Ang II, 10 of them belong to the AGC/CAM kinase family, and seven of the 23 Ang II kinases regulated by Ang II only are AGC/CAM kinases. Interestingly, the phosphorylation of a proline-directed Erk-like phosphorylation motif (PXp(S/T)P) was only significant among the peptides that were regulated by Ang II in a G␣ q -dependent manner (Fig. 3B, middle panel). This remarkable observation of divergence could reflect that Erk1/2 translocates very fast to the nucleus upon G␣ q -dependent activation to reach its nuclear substrates, whereas it is tied up in the cytosol by ␤-arrestin and activates AGC-and CAM-like kinases after G␣ qindependent activation. As we discuss in detail below, the PKD consensus motif ((L/V/I/M)XRXXp(S/T)) was by far the most predominant motif within the group of SII Ang II-regulated peptides and was also over-represented among the Ang II-regulated peptides, indicating substantial PKD activity. The AGC/CAM kinases CAMK2␦ and PKD have previously been found in association with ␤-arrestin 2, supporting their role in G␣ q -independent Ang II signaling (40).
To validate that the Mascot Score for peptide identification did not influence the data interpretation, we investigated whether the motif analysis was related to certainty of peptide identification by applying a filter on the sites identified based on a Mascot score cutoff of 48.9, equal to the upper 25% quantile and reapplying the motif analysis. The most significant feature was still the AGC/CAM motif for both sites regulated uniquely by Ang II (33% or 70 of 212 sites) and for sites regulated by both SII and Ang II (42% or 42 of 100 sites). This compares with a non-regulated background where this motif was found in only 14.8% or 331 of 2229 sites. The subset of these sites containing the PKD substrate motif was still significant for both uniquely Ang II-regulated sites (11% or 24 of 212 sites) and sites regulated by both SII and Ang II (14% or 14 of 100 sites) compared with a background of non-regulated sites (2.7% or 60 of 2229 sites). However, the Erk motif was no longer significant; although by closer analysis, the same trend was seen (8.9% of sites uniquely regulated by Ang II compared with a background of 4.8%). This analysis convincingly demonstrates that the Mascot score of a given peptide does not influence the likelihood of it being determined as regulated. The results also show that the motif analysis is limited in sensitivity given fewer sites but a differing Mascot score does not affect the motif extracted.
Transactivation of Membrane Receptors by Ang II-Ang II is known to influence activation and signal transduction of several other plasma membrane receptors, for example the EGF receptor (41,42). Our data show reproducible transphosphorylation of 10 membrane receptors, six of which depend on G␣ q protein activation, whereas four were phosphorylated independently of the G␣ q signaling (Table II). Ang II stimulation induced phosphorylation of both the insulin and IGF2 receptors. This is interesting, because several pieces of clinical and biological evidence indicate cross-talk between AT 1 R and insulin receptor signaling (43,44), but to our knowledge a direct link at the receptor level has not been reported. Furthermore, phosphorylation independent of G␣ q protein activation was also seen on three previously undiscovered FGF receptor phosphorylation sites. Ang II also resulted in phosphorylation of other 7TMRs including the ␤ 2 -adrenergic receptor, which is an AT 1 R interaction partner on a functional level (45).
Protein Kinase D Is a Key Protein Kinase in G␣ q Proteinindependent Ang II Signaling-The group of kinases that were phosphorylated after SII Ang II stimulation includes all three isoforms of PKD (Table I). Furthermore, we found 52 SII Ang II-up-regulated phosphorylation sites, which all conform to the PKD consensus phosphorylation motif (V/L/I)XRXX(S/T) Phosphoproteomics of AT 1 R Signaling (46), pointing to a substantial activation of PKD and phosphorylation of its substrates (Fig. 3B). The 52 putative PKD phosphorylation sites and protein substrates are depicted in supplemental Table VI. Because of the potentially large impact of PKD activation, we further explored how PKDs are regulated by G␣ q -independent signaling. First we validated whether the observed phosphorylation of PKD isoforms leads to their activation by performing Western blot analysis on lysates from Ang II-and SII Ang II-treated cells. We used a phosphospecific antibody against the autophosphorylation site, serine 910, which is known to correlate with the kinase activity of PKD1. As depicted in Fig. 4A, both Ang II and SII Ang II stimulation induced a robust phosphorylation of serine 910. Furthermore, the autophosphorylation of PKD1 and the activity of PKD1 measured in a PKD1 kinase activity assay (Fig. 4C) showed that both Ang II and SII Ang II rapidly and robustly activated the PKD1. PKD activation has been described to be mediated by novel PKC family members (47,48). PKCs are central Ang II-regulated kinases. It was previously shown that Ang II activates PKD through PKC␦ (47), but until now, this activation has been considered G␣ q -dependent. The SII Ang II-mediated PKD activation is also dependent on PKC␦ as it is was sensitive to the PKC inhibitor Gö 6983 and the PKC␦-specific inhibitor rottlerin (Fig. 4B). PKC␦ was activated by both Ang II and SII Ang II, which were evaluated by immunoblotting against phospho-Thr 507 (data not shown), which is the kinase domain activation site. In contrast to SII Ang II activation of PKD, activation by Ang II appears to be dependent on other PKC isotypes as it was more effectively blocked with the broader PKC inhibitor Gö 6983 than with the PKC␦-specific rottlerin (Fig. 4B).
A well known pathway leading to PKC␦ activation is activation of Rho and ROCK (49,50). Considering the large A, the regulation level of phosphorylation on protein kinases was very reproducible between the two biological replicates. A -fold change plot of the Ang II versus control ratio from the first biological replicate against the second reveals an impressive reproducibility resulting in a correlation of r 2 ϭ 0.9. B, G␣ q -dependently and -independently activated kinases phosphorylate distinct motifs. Top panel, among both the Ang II-and the SII Ang II-regulated peptides more than a third of the phosphorylation sites observed to be regulated were on the AGC/CAM kinase motif RXXp(S/T). Middle panel, the Erk1/2 kinase motif PXp(S/T)P is specific for the peptides that are only regulated by Ang II; i.e. of the G protein-dependently phosphorylated peptides, 83 of the peptides in the Ang II only group were phosphorylated on an Erk1/2 motif. Bottom panel, a further analysis of the SII Ang II-regulated peptides matching the AGC/CAM kinase motif showed that the dominant phosphorylation site within this group matched the PKD consensus motif ((V/L/I)XRXX(S/T)). * indicates statistically significant over-representation compared with non-regulated phosphopeptides. Statistically significant over-representation was calculated using the R implementation of Fisher's exact test. Ctrl, control; ANG, Ang II SII, SII Ang. amount of phosphorylated GEFs and GTPase-activating proteins in our data, we believe that the observed PKC␦ and PKD activation could be a result of ROCK activity. The mechanism of small GTPase activation by the AT 1 R has been poorly elucidated. The anchoring protein AKAP13, which is also highly phosphorylated by SII Ang II, has been shown to function as a GEF for Rho and simultaneously to interact with both PKC␦ and PKD, perhaps as a scaffold to bring the two kinases in close proximity (51). To test whether the observed PKD activation was ROCK-dependent, we measured PKD activation by the ability of PKD to phosphorylate syntide 2 in the presence of the ROCKspecific inhibitor Y27632. As depicted in Fig. 4C, the ROCK inhibitor inhibited 94% of the SII Ang II-mediated PKD activation and only 45% of the Ang II-mediated PKD activation. These data suggest a model where Rho/ROCK/PKC␦/PKD is activated independently of G␣ q activation (Fig. 4D).

DISCUSSION
In this study, we identified more than 500 phosphoproteins regulated by Ang II of which the majority has not previously been implicated in Ang II signaling (supplemental Tables I and II and Fig. 2). Importantly, our data reveal a previously unexpected amount and versatility of G␣ q -independent signaling because SII Ang II regulated 36% of the AT 1 R-regulated phosphosites and affected the same cellular functions as the G␣ q -dependent part of the signaling (Fig. 2). Furthermore, G␣ q -dependent and -independent AT 1 R signaling pathways are very distinct as reflected by activation of different kinds of kinases. This was observed as both phosphorylation of the actual kinases and "fingerprints" that activated kinases leave on other phosphorylated substrate proteins due to their peptide motif specificity (Table I and Fig. 3B).
Quite unexpectedly, we observed that Ang II signaling almost exclusively occurs through serine (88.3%) and threonine (11.1%) phosphorylation ( Fig. 2B and supplemental Tables I and II). We observed very little tyrosine phosphorylation (0.6%) and only identified regulation of five tyrosine kinases (EPHA2, FYN, insulin receptor, RET, and FGFR3) (Table I). However, there is no doubt that other tyrosine kinases such as the Jak and Src kinases as well as EGF receptor tyrosine kinase are important mediators of AT 1 R signaling in certain cellular conditions (6,11,13,(52)(53)(54). In a similar large scale phosphoproteomics analysis of the EGF receptor, we found that 6% of the regulated phosphosites were Tyr(P) (26). This provides assurance that our method is able to detect tyrosine phosphorylation sites, and therefore our data indicate that EGF receptor transactivation and tyrosine kinase signaling may play a minor role in AT 1 R signaling in HEK293 cells that might to some extent be due to low expression of EGF receptors.
This unbiased mass spectrometric analysis uncovered unexpected regulation of specific phosphorylation patterns and proteins in the two AT 1 R signaling pathways. We showed that PKD was activated by SII Ang II through PKC␦-and ROCKdependent signal transduction, whereas Ang II activation of PKD via G␣ q protein may also depend on other PKCs (Fig.  4D). These data challenge the understanding of PKCs as strictly G␣ q -dependent regulated kinases in the Ang II system and place PKD as a very central hub in G␣ q -independent mediated signal transduction from the AT 1 R. PKD has not attracted much attention in cardiovascular signaling compared with its relatives in the PKC family. However, emerging data point to a central role for PKDs in AT 1 R-regulated signal transduction, resulting in the induction of cardiac hypertrophy as well as in protection against apoptosis, migration, remodeling, and proliferation (55,56). The regulation of these very different cellular outcomes is presumably related to distinct cellular localization of the PKD isoforms and organization into different signalosomes (56). Interestingly, this study showed that PKD can also be regulated by G␣ q -independent mechanisms via the Rho/ROCK pathway. A substantial fraction of the Ang II-regulated proteins displayed here is involved in small GTPase activation, plasma membrane remodeling, and cytoskeletal rearrangements. In general, we found many GEFs and GTPase-activating proteins in both groups, indicating a substantial activity of small GTPase signaling. We believe that a focus on small GTPases may reveal a novel branch of the G protein-independent signaling network.
In recent years, research on Ang II has been focused on the cellular phenotypes resulting from G protein-dependent versus ␤-arrestin-dependent signal transduction. Cellular outcomes such as gene expression changes, proliferation, migration, apoptosis, and hypertrophy have been linked to Ang II signaling (22,24,40,57). A study of Ang II-mediated transcriptional changes concluded that transcriptional regulation is dependent on G␣ q activation and inhibited by ␤-arrestins (58). In support of this observation, we found that Ang II up-regulated phosphorylation of several transcription factors, for example c-JUN, HOXA3, and EP400 in a strictly G␣ q -dependent manner (supplemental Table I). On the other hand, we also measured an extensive SII Ang II-stimulated phosphorylation site regulation of nuclear proteins, strongly indicating that the G␣ q -independent signal is somehow propagated to the nucleus (supplemental Table II). As the ␤-arrestindependent signal to the nucleus does presumably not induce transcriptional activation, it could instead mediate transcriptional repression. Indeed, we observed SII Ang II-regulated phosphorylation of transcriptional regulators such as LARP, MIER, and TRIM33 that are known to mediate transcriptional repression. SII Ang II also regulated the programmed cell death protein 4 (PDCD4) at several sites including Ser 67 . It has been shown that upon Ser 67 phosphorylation PDCD4 translocates to the nucleus and represses c-JUN-and ␤-cateninmediated transcription (59). PDCD4 has also been proposed to form a complex with ␤-arrestins 1 and 2 (40). The lack of transcriptional activity upon ␤-arrestin signaling may also be due to nuclear exclusion of proteins central for progression of the transcriptional process. One example is CAMK2␦, which is involved in both growth responses and transcriptional regulation. CAMK2␦ is phosphorylated on multiple sites by SII Ang II stimulation in a cluster (Ser 330 -Thr 337 ) adjacent to the nuclear localization signal, possibly inhibiting nuclear translocation of CAMK2␦ (60).
SII Ang II stimulates protein synthesis in both HEK293 and vascular smooth muscle cells by activation of the translation initiator EIF4A (61). In agreement with this finding, we observed phosphorylation on EIF4B and EIF5B on multiple sites by both Ang II and SII Ang II. We also observed Ang IIstimulated phosphorylation on 12 proteins involved in the ubiquitination-mediated protein degradation process (supplemental Table I). Only three of these proteins (calpastatin and the ubiquitin ligases TRIM33 and PRAJA1) were affected by SII Ang II stimulation (supplemental Table II).   The statistical difference between Ang II and SII Ang II with and without inhibitors was tested with a one-tailed paired Student's t test. *, p Ͻ 0.05. Data are presented as the mean and error bars reflect Ϯ S.E. of three experiments. C, kinase assay showing phosphorylation of the PKD substrate peptide syntide 2 in the presence and absence of 30 M ROCK inhibitor Y27632. AT 1 R-HEK cells were starved overnight and pretreated with Y27632 for 30 min before stimulation with 10 nM Ang II and 1.87 M SII Ang II for 3 min. Lysates were immunoprecipitated with anti-PKD antibody, and syntide 2 peptide was subjected to phosphorylation for 10 min in the presence of [ 32 P]ATP. Syntide 2 was spotted to Whatmann paper and washed thoroughly before autoradiography. The statistical difference between Ang II and SII Ang II with and without inhibitors was tested with a one-tailed paired Student's t test. *, p Ͻ 0.05. Data are presented as the mean and error bars reflect Ϯ S.E. of three experiments. D, schematic indicating the G protein-independent pathway to PKD activation going through Rho/ ROCK/PKC␦. ANG, Ang II; PMA, phorbol 12-myristate 13-acetate; max, maximum.
In the Ang II signaling system, G␣ q -dependent mechanisms have been reported to increase apoptosis, whereas ␤-arrestin-mediated signaling may protect against apoptosis (57). We found ligand-regulated phosphorylation of both pro-and antiapoptotic proteins by Ang II (e.g. proapoptotic protein reticulon 4 (NOGO) and the antiapoptotic BCL-2 family member BIS) as well as SII Ang II (e.g. the well known antiapoptotic kinases RSK2 and p70S6K and the antiapoptotic protein PEA-15 as well as the proapoptotic DIDO-1).
We found a variety of proteins involved in cell cycle regulation to be involved in both G␣ q -dependent and -independent signaling. DNA polymerase 1 (POLA1) was phosphorylated in a strictly G␣ q -dependent manner on Ser 192 , a previous unidentified phosphorylation site. The nuclear interaction partner of anaplastic lymphoma kinase (NIPA) was phosphorylated by SII Ang II on a cluster of 7 residues including Ser 354 and Ser 359 , which most likely results in degradation of the nuclear interaction partner of anaplastic lymphoma kinase and accumulation of cyclin B in late G 2 , which in theory leads to cell cycle progression (27). However, SII Ang II did not result in proliferation in the AT 1 R-HEK cells (data not shown) as was observed in cardiac myocytes and NIH3T3 cells (22,39).
Upon AT 1 R activation, the formation of clathrin-coated pits and receptor endocytosis are initiated, which require a considerable rearrangement of the cytoskeleton. Here, we observed that dynamin, filamin A, and vimentin, which are cen-tral filamentous proteins involved in endocytosis, F-actin polymerization, and migration, were regulated strictly by Ang II. Contradictory, cellular migration toward an attracting chemical factor (chemotaxis) has been directly related to ␤-arrestin 2-dependent AT 1 R signaling (24). Numerous proteins involved in the endocytosis and cytoskeletal reorganization were regulated by SII Ang II including GIT1, REPS1, and PACS1. CONCLUSIONS This is the first study dissecting the differences between a full agonist and a biased agonist of a 7TMR by mass spectrometry. Using a high resolution and quantitative mass spectrometric approach, we identified more than 500 phosphoproteins regulated by Ang II of which the majority has not previously been implicated in Ang II signaling. Our data reveal a previously unexpected amount and versatility of G␣ q -independent signaling from this 7TMR because SII Ang II affected 36% of the identified AT 1 R signaling network and influenced all Ang II-regulated cellular events. In addition, we also established that the two signaling pathways are in fact very distinct in their nature, which is reflected by activation of different protein kinases. The shotgun mass spectrometric approach did not uncover all proteins in the network but provided multiple pieces of the puzzle. For many proteins such as PKD, this unbiased analysis uncovered unexpected regulation patterns The illustration depicts the novel AT 1 R-regulated kinases identified in this study. Note that the kinases illustrated as G protein-independent may be activated by both pathways, whereas the kinases regulated strictly by Ang II are considered as G protein-dependent. ␤-ARR, ␤-arrestin. and exposed proteins not previously considered to be a part of the Ang II signaling network (Fig. 5). We believe that this mapping of the Ang II signalosome will serve as a basis for future studies on the regulation of Ang II signaling molecules as well as provide a broader and more in-depth overview of the effect on cellular outcomes mediated by different AT 1 R signaling pathways, which could be important for the pharmacological use of biased agonists in cardiovascular diseases. Looking at AT 1 R signaling as a model system for 7TMR signaling in general, the amount and diversity of G␣ q -independent signaling may be larger than previously acknowledged in other systems as well. Moreover, our study holds promise for quantitative mass spectrometry as a universal tool in the evaluation and analysis of biased agonists to other receptors in the future. The raw data associated with this manuscript may be downloaded from ProteomeCommons.org Tranche, https://proteomecommons. org/tranche/, using the following hash codes: all spectra: AkG/ I6DDb4N6PFouQySWFB7f0vNPdfzVN0sQ/wDFPxdBDYschQpϩ 0fvbHAd6C9MWlHaϩe1hOm9ZK3X1itkM7m/B1MywAAAAAAABVtwϭϭ; pass phrase, ANG2gitte; and annotated spectra: 8dTCR-89xDMxbL7j4fbhHDxs/6463EzqsuTfURzr/NBCG3b/3d0Sb8UWDFm-7khvibmcuhEV1sZYrufaPR7uxmDwY68mYAAAAAAAACzwϭϭ; pass phrase, ANG2gitte.
□ S This article contains the following supplemental data: Table I,   phosphopeptides regulated by Ang II alone; Table II, phosphopeptides regulated by ANG II plus SII Ang II; Table III, phosphopeptides  regulated by SII Ang II alone; Table IV, all regulated phosphopeptides;  Table V, all peptides; Table VI, protein kinase D motif peptides; Table  VII, all proteins; Table VIII,  ʈ Both authors contributed equally to this work.