Abstract
Water excretion by the kidney is regulated by the neurohypophyseal peptide hormone vasopressin through actions in renal collecting duct cells to regulate the water channel protein aquaporin-2. Vasopressin signaling is initiated by binding to a G-protein–coupled receptor called V2R, which signals through heterotrimeric G-protein subunit Gsα, adenylyl cyclase 6, and activation of the cAMP-regulated protein kinase (PKA). Signaling events coupling PKA activation and aquaporin-2 regulation were largely unknown until the advent of modern protein mass spectrometry techniques that allow proteome-wide quantification of protein phosphorylation changes (phosphoproteomics). This short review documents phosphoproteomic findings in collecting duct cells describing the response to V2R-selective vasopressin agonists and antagonists, the response to CRISPR-mediated deletion of PKA, results from in vitro phosphorylation studies using recombinant PKA, the response to the broad-spectrum kinase inhibitor H89 (N-[2-p-bromocinnamylamino-ethyl]-5-isoquinolinesulphonamide), and the responses underlying lithium-induced nephrogenic diabetes insipidus. These phosphoproteomic data sets have been made available online for modeling vasopressin signaling and signaling downstream from other G-protein-coupled receptors.
SIGNIFICANCE STATEMENT New developments in protein mass spectrometry are facilitating progress in identification of signaling networks. Using mass spectrometry, it is now possible to identify and quantify thousands of phosphorylation sites in a given cell type (phosphoproteomics). The authors describe the use of phosphoproteomics technology to identify signaling mechanisms downstream from a G-protein-coupled receptor, the vasopressin V2 subtype receptor, and its role of the regulation and dysregulation of water excretion in the kidney. Data from multiple phosphoproteomic data sets are provided as web-based resources.
Introduction
Renal water excretion is controlled by the peptide hormone vasopressin largely through regulation of the water channel aquaporin-2 (AQP2) in collecting duct principal cells (Fig. 1). The response is mediated by a G-protein–coupled receptor, the vasopressin V2 receptor (V2R, gene symbol: Avpr2). The V2R signals via the heterotrimeric G-protein stimulatory α subunit, Gsα (Erlenbach et al., 2001), which links to activation of adenylyl cyclase 6 in collecting duct cells (Roos et al., 2012), resulting in increased cAMP production. The V2R can also signal through binding of β arrestins, which act as scaffolds for elements of the ERK mitogen-activated protein kinase cascade, thereby increasing ERK activation in response to ligand binding (Tohgo et al., 2003) (See Effect of Vasopressin on the Phosphoproteome of the Renal Collecting Duct and V2 Receptor Antagonists and Hyponatremic Disorders: How Vaptans Affect the Collecting Duct Phosphoproteome.).
Renal collecting ducts can be microdissected from rodent or rabbit kidneys and perfused in vitro to study the water permeability response to vasopressin (Burg et al., 1966). Using this technique, it has been observed that addition of vasopressin results in a rapid increase in transepithelial water permeability (Grantham and Burg, 1966; Wall et al., 1992). The increase begins after an approximately 40-second delay and requires approximately 20 minutes for completion (Wall et al., 1992). If isolated perfused collecting ducts are exposed to cAMP analogs, the transepithelial water permeability increases with a time course similar to that seen with vasopressin (Grantham and Burg, 1966; Star et al., 1988; Wall et al., 1992), implying that the signaling processes that are essential for the water permeability response to vasopressin are cAMP-dependent. The two known effectors of cAMP in collecting duct cells are protein kinase A (PKA) and exchange protein directly activated by cAMP (EPAC) 1 and EPAC2 (coded by Rapgef3 and Rapgef4). Deletion of PKA in collecting duct cells abolishes nearly all of the protein phosphorylation changes associated with vasopressin stimulation, in addition to eliminating Aqp2 gene transcription (Isobe et al., 2017; Datta et al., 2020). On the other hand, EPAC gene deletions did not have a pronounced effect on AQP2 regulation when knocked out in mice (Cherezova et al., 2019), although EPAC has been implicated in vasopressin-induced calcium mobilization (Yip, 2006). Thus, in collecting duct cells, the major aspects of vasopressin signaling are mediated by PKA.
Knowledge of the upstream aspects of vasopressin signaling in the collecting duct through activation of PKA is well established (Fig. 1). Signaling downstream from PKA leading to regulation of AQP2 activity in collecting duct cells is less well understood. Discovery of the relevant signaling pathways in collecting duct cells and how they are changed in disorders of water balance is of interest because of the high prevalence of such disorders. A relatively new tool that is shedding light on vasopressin signaling in health and disease is protein mass spectrometry, especially the use of protein mass spectrometry to identify and quantify phosphorylation events in the cell (phosphoproteomics). The goal of this review is to summarize progress from phosphoproteomic analysis of AQP2-expressing collecting duct cells to identify signaling changes seen in response to vasopressin, signaling changes resulting from deletion of PKA in collecting duct cells, identification of PKA substrates, signaling changes seen in polyuric disorders, and signaling changes associated with water-retaining (hyponatremic) disorders. We also summarize essential features of phosphoproteomic methodology, including bioinformatic analysis and data integration.
Cellular Physiology of AQP2-Expressing Renal Collecting Duct Cells
Detailed studies over several decades have identified a number of cellular-level processes that are regulated in response to vasopressin in collecting duct cells (Table 1) (Ganote et al., 1968; Kirk et al., 1984; Star et al., 1988; Mishler et al., 1990; Champigneulle et al., 1993; Nielsen and Knepper, 1993; Nielsen et al., 1993, 1995; Simon et al., 1993; Naruse et al., 1995; Sabolić et al., 1995; Chou et al., 2000, 2004, 2008; Tamma et al., 2001; Brown, 2003; Yamaguchi et al., 2003; Nunes et al., 2008; Hasler et al., 2009; Nedvetsky et al., 2010; Khositseth et al., 2011; Miller et al., 2013; Sandoval et al., 2013, 2016; Loo et al., 2013). The most thoroughly studied of these processes are the actions of vasopressin to increase Aqp2 gene transcription (Hasler et al., 2002; Sandoval et al., 2016) and the action of vasopressin to regulate membrane trafficking of the AQP2 protein to increase its abundance in the plasma membrane of collecting duct cells (Nielsen et al., 1995). Detailed discussion of the cellular level responses in Table 1 is beyond the scope of this treatise but can be found in prior review articles (Knepper, 1997; Knepper and Inoue, 1997; Brown et al., 1998; Sasaki et al., 1998; Nielsen et al., 1999, 2002; Verkman, 1999; Klussmann and Rosenthal, 2001; Brown, 2003; Noda and Sasaki, 2005; Valenti et al., 2005; Bichet, 2006; Boone and Deen, 2008; Moeller and Fenton, 2012; Fenton et al., 2013; Jung and Kwon, 2016, 2019). One goal of this review is to map recently obtained phosphoproteomic data from both native and cultured collecting duct cells to the processes summarized in Table 1 to identify the signaling events responsible for the physiologic responses. The phosphoproteomic data sets discussed in this review are listed (with their hyperlink addresses) in Table 2. The data sets can also be interrogated at the Kidney Systems Biology Project website (https://hpcwebapps.cit.nih.gov/ESBL/Database/), which allows users to browse and search a large number of data sets acquired from proteomics and next-generation sequencing studies.
Phosphoproteomic Methodology and Bioinformatics
Full details of the methodology used for phosphoproteomic analysis are beyond the scope of this article. The underlying principles of the method are summarized in Fig. 2. The approach is a so-called bottom-up methodology in which proteins are digested with proteases having specific cleavage sites (most commonly trypsin, which cleaves at the carboxyl-terminal side of the amino acids lysine and arginine), and the resulting peptides are sequenced by tandem mass spectrometry coupled to a chromatographic method that stratifies peptide delivery to the mass spectrometer, so-called LC-MS/MS. For phosphoproteomics, additional chromatographic steps are included for phosphopeptide enrichment based on the negative charges of the added phosphate groups (Fig. 2). The mass spectrometer readout is a list of mass-to-charge ratio peaks and their intensities as a function of time. The amino acid sequence for a given peptide is typically determined by fragmentation (e.g., via high-energy collisions of the peptide ions with an inert gas), which produces a series of mass-to-charge ratio peaks that can be used to determine the sequence and phosphorylation site(s) by pattern matching. Many proprietary and open-source programs are available for this task.
Quantification is usually achieved using labeling techniques. For cells in culture, the most widely used method is SILAC (Ong et al., 2002). With this method, cells are grown in media containing amino acids in which different stable isotopes are incorporated, e.g., C12 and N14 in control cells and C13 and N15 in vasopressin-treated cells. This allows the same peptide to be seen by the mass spectrometer twice (light peptides for control cells, heavy peptides for vasopressin-treated cells in the above example). The relative abundance of a particular peptide in the two samples can be quantified by comparing the peak intensities.
A different quantification strategy, based on isobaric tagging, is usually used in phosphoproteomic analysis of native tissues (Cheng et al., 2016). An example of isobaric tagging methodology is TM (ThermoFisher Inc.), which allows simultaneous quantification of up to 11 multiplexed samples in a single LC-MS/MS run. With this method, chemical tags of equal molecular mass are covalently attached to the processed peptides, a different one for each sample, allowing the samples to be combined. The tagged peptides are fragmented in the mass spectrometer to generate a series of sample-specific reporter ion peaks whose heights give measures of the phosphopeptide abundance in each of the multiplexed input samples. This method is cost-effective because of multiplexing. It allows a high degree of precision and high throughput, allowing tens of thousands of phosphopeptides to be quantified in each LC-MS/MS run.
The tangible product from a quantitative phosphoproteomic analysis is a list of phosphopeptides with the following attributes: 1) parent protein indicated by official gene symbol, 2) the identified phosphorylation site(s), and 3) a quantitation value typically indicated as a ratio between experimental and control values. Added to this are statistical parameters that help to rank the responses with regard to the likelihood that a given response is not a false positive. Typically, three or more replicates are required for reliable statistical analysis.
The bioinformatic challenge is to map the phosphopeptide list to underlying biologic processes to generate hypotheses about molecular mechanisms. For this review, the relevant biologic processes are the ones listed in Table 1. One such process is vasopressin-mediated depolymerization of filamentous actin. The task, then, is to identify the phosphoproteins in a given data set that are involved in regulation of actin polymerization. This can be done with Gene Ontology (GO) biologic function terms. Biologic function terms for a set of regulated phosphoproteins can be identified by inputting the corresponding official gene symbols into computer programs such as the Automated Bioinformatics Extractor (https://helixweb.nih.gov/ESBL/ABE/). This program will provide all GO terms for each regulated phosphoprotein. After transferring this information to a spreadsheet, it is relatively easy to identify which phosphoproteins have GO terms involving actin polymerization. These proteins and their regulated phosphorylation sites can then be studied further, e.g., through use of genome editing techniques (CRISPR-Cas9) (Isobe et al., 2020) to block gene expression or to mutate the phosphorylated amino acid. Usually, however, prior to such reductionist approaches, it is wise to use Bayesian analysis (Bradford et al., 2014; Xue et al., 2017) to assess the strength of prior data pointing to a particular phosphoprotein as a candidate for a mechanistic role in a targeted biologic process.
Another important bioinformatic task is the identification of the protein kinases responsible for observed changes in phosphorylation. It is well known that many protein kinases have sequence preferences. For example, PKA prefers serines or threonines in a motif (R/K)-(R/K)-X-p(S/T)-, where (R/K) indicates either arginine or lysine, X indicates any amino acid, and p(S/T) indicates phosphorylated serine or threonine. Thus, the sequence surrounding a phosphorylation site can be used to narrow down the possible kinases that could phosphorylate that site. An equally important factor that determines kinase/target interactions is colocalization in the cell. Proteomic analysis of subcellular fractions (Bolger et al., 2012; Schenk et al., 2012; Yang et al., 2015; Pickering et al., 2016) can be used to find candidate protein kinases that could target particular phosphoproteins based on colocalization.
Often, regulated phosphoproteins are themselves protein kinases, and these kinases become candidates for mediation of regulated phosphorylation. In a minority of cases, the effect of a given phosphorylation event on kinase activity is known from prior research; these sites can often be identified using data available on two commercial websites, namely PhosphositePlus (https://www.phosphosite.org/) and PhosphoNet (http://www.phosphonet.ca/). When such prior data are available, the change of kinase activity in a particular experiment can be inferred.
To provide a tool for cross-comparing the different data sets described in this paper (Table 2), we used a database program, MySQL, and provided a web-based interrogation interface (https://big.nhlbi.nih.gov/). Users can enter official gene symbols in the submission box, and the output will be all of the quantified sites for that protein in any of the phosphoproteomic studies reported in this paper. Other types of -omics data are included for reference. For example, entering "Ctnnb1" (β catenin) revealed that Ser552 phosphorylation is increased by vasopressin in both native inner medullary collecting duct cells and cultured mouse mpkCCD cells, whereas the V2R antagonist decreases phosphorylation, and lithium treatment had no effect.
Effect of Vasopressin on the Phosphoproteome of the Renal Collecting Duct
The simplest experiment to identify vasopressin signaling pathways in the renal collecting duct is to expose collecting duct cells to vasopressin or its vehicle and carry out mass spectrometry–based quantitative phosphoproteomics. Experiments following this strategy have been done in both cultured collecting duct cells (mouse mpkCCD) (Rinschen et al., 2010; Datta et al., 2020) and suspensions of native inner medullary collecting duct cells from rats (Hoffert et al., 2006, 2012; Deshpande et al., 2019). The studies have used a V2R–selective vasopressin analog desmopressin, D-amino D-arginine vasopressin (dDAVP), which is used clinically for treatment of central diabetes insipidus (Christensen and Rittig, 2006; Oiso et al., 2013). Over the years, protein mass spectrometry has become more and more sensitive, resulting in the ability to quantify more and more phosphorylation sites (Table 3). Despite a relative lack of sensitivity in the earliest phosphoproteomic studies of the vasopressin response in collecting duct cells, they provided critical information that spurred progress, e.g., the identification of a cluster of four vasopressin-regulated phosphorylation sites in the COOH-terminal tail of AQP2 that are critical to regulation of AQP2 trafficking (Hoffert et al., 2006, 2008, 2012; Bansal et al., 2010; Rinschen et al., 2010). However, only in the last 2 or 3 years have we begun to approach comprehensive phosphorylation site quantification in phosphoproteomic analysis.
A general finding in all phosphoproteomic studies of the vasopressin response in collecting duct cells is that only a small fraction (typically around 2%) of phosphorylation sites among those quantified exhibit changes in response to vasopressin. This finding is compatible with the idea that control of physiologic processes via GPCR-dependent signaling is highly selective with regard to phosphorylation targets. Hence, comprehensive analysis of phosphorylation changes has the potential to uncover specific cellular processes that are selectively regulated. Motif analysis in these studies demonstrates that the phosphorylation sites increased in abundance by vasopressin are dominated by those that fit the motif –(R/K)-(R/K)-X-p(S/T), where (R/K) means either arginine or lysine, X means any amino acid, and p(S/T) means phosphorylation of either serine or threonine. This motif is similar to that associated with the action of PKA but is also consistent with a number of protein kinases, including protein kinase G (Prkg1 or Prkg2) (Miller et al., 2008), myotonic dystrophy protein kinase (Dmpk or related kinases) (Miller et al., 2008), and death-associated protein kinase (Dapk1, Dapk2, or Dapk3) (Douglass et al., 2012). In contrast, the phosphorylation sites decreased in abundance in response to vasopressin are dominated by those that fit the motif -X-p(S/T)-P, where P is proline, here in position +1 relative to the phosphorylated residue. Such phosphorylation sites are typically targets of members of two families of protein kinases: the cyclin-dependent kinases (CDKs) and the mitogen-activated kinases (MAPKs). This finding suggests that vasopressin signaling results in inactivation of one or more of these “proline-directed” kinases. Indeed, vasopressin has been found to reduce activation of both ERK1 and ERK2 (Mapk3 and Mapk1, respectively) in collecting duct cells (Pisitkun et al., 2008; Rinschen et al., 2010). The fall in ERK activity in response to vasopressin contrasts with β arrestin–mediated GPCR-dependent activation of the MAPK signaling seen with some GPCRs (Luttrell et al., 2001). Therefore, these results suggest that, with regard to vasopressin signaling in the collecting duct, cAMP-mediated signaling dominates over β arrestin signaling. There is also evidence that vasopressin increases the activity of a proline-directed protein phosphatase called PP2A, through PKA-mediated phosphorylation of a phosphatase inhibitor, cAMP-regulated phosphoprotein 19 (Arpp19) (Deshpande et al., 2019), which could explain some of the downregulated -X-p(S/T)-P sites seen in response to vasopressin.
Vasopressin V2 receptors are also present in two renal cell types that do not express AQP2, namely distal convoluted tubule cells and thick ascending limb cells (Lee et al., 2015), in which vasopressin regulates Na-Cl reabsorption from the tubule lumen to the blood. Large-scale phosphoproteomic studies in cultured distal convoluted tubule cells (Cheng et al., 2015) and suspensions of thick ascending limb cells (Gunaratne et al., 2010) showed that vasopressin elicits a response that is very similar to that seen in collecting duct cells, with increases at basic sites with the motif –(R/K)-(R/K)-X-p(S/T) and decreases in sites with the consensus proline-dependent motif -X-p(S/T)-P. Thus, the general characteristics of the response to vasopressin acting through the V2 receptor are independent of target cell type. Furthermore, parathyroid hormone binding to the parathyroid hormone receptor, a Gsα-coupled receptor, in MC3T3-E1 preosteoblast cells produced the same general pattern of phosphorylation changes (Williams et al., 2016), demonstrating that this response pattern is not unique to the V2 receptor or to the kidney.
Figure 3 shows some of the results from the most comprehensive vasopressin-response phosphoproteomic study in native collecting duct cells (Deshpande et al., 2019). Several of the phosphoprotein targets of vasopressin signaling mapped to cellular physiologic functions corresponding to those in Table 1. Many of those that had basophilic phosphorylation sites that increased with vasopressin are known targets of PKA (underlined, Fig. 3), cementing the conclusion that vasopressin signaling is largely mediated by PKA activation. Nevertheless, some elements of vasopressin signaling appear to be PKA-independent (Datta et al., 2020). This includes PKA-independent, vasopressin-induced phosphorylation changes associated with activation of one or more kinases of the AMPK-related kinase family (also known as SNF1-family kinases). It also includes increased phosphorylation of the protein kinase coded by Stk39 (commonly called "SPAK") at a site associated with its activation and decreased phosphorylation and decreased activation of atypical protein kinase C iota.
However, in the study of Deshpande et al. (2019), phosphorylation changes were seen in several other protein kinases at sites that mediate changes in catalytic activity. Specifically, protein kinase D1, PCTAIRE domain kinase 1 (Cdk16), PCTAIRE domain kinase 3 (Cdk18), and Src undergo phosphorylation changes consistent with increases in activity in response to vasopressin, whereas Araf and Raf1 undergo apparent decreases in activity in response to vasopressin. Protein kinase D1, Cdk16, and Cdk18 are involved in regulation of membrane trafficking, whereas Src, Araf, and Raf1 regulate MAP kinase activities (Deshpande et al., 2019). The full data set describing the phosphoproteomic response to vasopressin in native rat collecting ducts is provided on a publicly accessible website (https://hpcwebapps.cit.nih.gov/ESBL/Database/IMCD-Phos/).
Because protein kinases must come into physical contact with their targets to phosphorylate them, phosphoproteomic analysis can be construed as a form of proximity analysis, potentially identifying subcompartments in the cell that contain the regulated forms of the relevant protein kinases. For example, in the Deshpande et al. (2019) analysis, it was concluded that PKA phosphorylates its targets predominantly at membrane surfaces, particularly targeting plasma membrane proteins.
Effect of PKA Deletion on the Phosphoproteome of the Renal Collecting Duct
It is widely accepted that the downstream effects of activation of the V2 vasopressin receptor are mediated chiefly by PKA (Fig. 1). To identify direct and indirect phosphorylation targets of PKA in collecting duct cells, Isobe et al. (2017) used a combination of genome editing and phosphoproteomics. CRISPR-Cas9 was used to eliminate PKA expression (deletion of both PKA catalytic α and PKA catalytic β) in mouse mpkCCD cells, a cultured cell model of the kidney collecting duct principal cell. The wild-type mpkCCD cells have been previously shown to display both long-term (transcriptional) and short-term (membrane trafficking) regulation of aquaporin-2 in response to vasopressin (Yu et al., 2009). The CRISPR/phosphoproteomics experiments identified 229 direct PKA phosphorylation sites, many of which were novel (Isobe et al., 2017). These phosphoproteins are mapped to specific cellular structures or functions in Fig. 4. Note that the indicated protein groups correspond closely to the known cellular effects of vasopressin in collecting duct cells (Table 1).
The group 1 PKA targets (“Transcriptional Regulation,” Fig. 4) include two transcription factors, designated by official gene symbols Nfatc2 and Nfkb1. Members of the NFATc (Li et al., 2007) and NFκB (Hasler et al., 2008) transcription factor families have already been shown to regulate Aqp2 gene transcription. NFATc transcription factors are also regulated by PKA in another way (https://hpcwebapps.cit.nih.gov/ESBL/PKANetwork/Transcription.html), namely by increasing intracellular calcium through phosphorylation of inositol 1,4,5-triphosphate receptor 3 (Gene symbol: Itpr3), a calcium release channel. The vasopressin-mediated calcium mobilization in collecting duct cells is seen as calcium spikes that increase in frequency with addition of vasopressin (Yip, 2002; Pisitkun et al., 2008). The increased calcium activates the phosphatase calcineurin, which dephosphorylates NFATc, causing its translocation into the nucleus. Another group 1 PKA target, β catenin (Ctnnb1), functions as a transcriptional coregulator in the Wnt signaling pathway, which is important to collecting duct maturation and differentiation (Schmidt-Ott and Barasch, 2008).
Identification of Potential PKA Targets by In Vitro Phosphorylation and Mass Spectrometry
Another approach to the identification of potential PKA targets is to carry out in vitro incubation of dephosphorylated proteins with the recombinant active kinase, followed by mass spectrometry–based phosphoproteomics. This has been done in two studies aimed at identifying phosphorylation pathways involved in the regulation of aquaporin-2. One study used a mixture of dephosphorylated proteins from rat whole kidney, liver, brain, and small intestine (Douglass et al., 2012), whereas the other used dephosphorylated proteins isolated from rat inner medullary collecting duct (Bradford et al., 2014), incubating each with recombinant PKA catalytic α subunit. The data are available at https://hpcwebapps.cit.nih.gov/ESBL/Database/PKA-in-vitro/.
In general, these studies confirmed an R-R-X-p(S/T) substrate preference for PKA but also identified many PKA phosphorylation target sites that deviated from this motif, albeit with basic amino acids in at least one upstream position within four amino acids from the phosphorylated serine or threonine. It is recognized that, in the functioning cell, whether a kinase targets a particular substrate is dependent on two factors: kinase specificity and colocalization of the kinase with the substrate (Linding et al., 2008). Thus, in vitro phosphorylation of a particular phosphorylation site by PKA does not necessarily mean that that site will be phosphorylated by PKA in the intact cell.
Effect of Kinase Inhibitor H89 on Protein Phosphorylation in PKA-Null Collecting Duct Cells
N-[2-p-bromocinnamylamino-ethyl]-5-isoquinolinesulphonamide (H89) is often used as a PKA specific inhibitor to study the involvement of PKA in signaling pathways. However, evidence from cell-free experiments has suggested that H89 can also inhibit other protein kinases (Davies et al., 2000). Limbutara et al. (2019) used PKA-null and PKA-intact mouse cell lines derived from mpkCCD cells (see Effect of PKA Deletion on the Phosphoproteome of the Renal Collecting Duct) and quantitative phosphoproteomics to investigate the specificity of H89 over the range of concentrations commonly used in the literature. From a total of 14,507 phosphorylation sites quantified, the authors found that 402 phosphorylation sites were significantly changed in abundance in PKA-intact cells, and 217 sites were significantly changed in PKA-null cells. Concentration-response data are available at https://esbl.nhlbi.nih.gov/H89/. Analyses of sequence logos generated from significantly decreased phosphorylation sites in PKA-intact and PKA-null cells both revealed a preference for basic amino acids at position −3 and −2. Thus, H89 appears to inhibit basophilic kinases other than PKA in intact cells. Likely H89 targets include basophilic protein kinases such as protein kinase B, ribosomal S6 inase, AMP-activated protein kinase, and Rho kinase. Thus, the finding that a particular process is altered by H89 cannot be considered sufficient to conclude that PKA is involved in that process.
Polyuric Disorders: Phosphoproteomics of Lithium-Induced Nephrogenic Diabetes Insipidus
Several disorders of water balance are attributable to loss of vasopressin-mediated regulation of AQP2 in collecting duct cells. Phosphoproteomics can be used in animal models of these disorders to ascertain which signaling pathways are perturbed. In diabetes insipidus, patients excrete excessive amounts of water (polyuria) in an uncontrolled fashion. An example is lithium-induced nephrogenic diabetes insipidus (NDI). Lithium salts provide an effective treatment in many patients with psychiatric manic-depressive disorder (bipolar disorder), but the overall efficacy of lithium salts is limited by concomitant NDI (Knepper et al., 2015). Animal studies have revealed that the polyuria in lithium-induced NDI is due to a nearly complete loss of AQP2 protein from collecting duct cells (Marples et al., 1995) due to loss of Aqp2 gene transcription (Kortenoeven et al., 2012) associated with entry of the cells into the cell cycle with G2 arrest (de Groot et al., 2014). Because the protein kinase glycogen synthase kinase 3-beta (Gsk3β) is directly inhibited by lithium ions (Klein and Melton, 1996), prior studies have focused on the role of Gsk3β in the regulation of Aqp2 gene expression, showing that Gsk3β inhibition or deletion negatively regulates adenylyl cyclase activity (Rao et al., 2010). To carry out an unbiased assessment of signaling changes in collecting duct cells in response to short-term lithium administration in rats, we carried out quantitative phosphoproteomics in native inner medullary collecting ducts (Trepiccione et al., 2014). The study quantified 1093 unique phosphopeptides, including 152 that were increased and 56 that were decreased in response to lithium. The full data set is available online at https://helixweb.nih.gov/ESBL/Database/iPALT/. The phosphoprotein groups with phosphorylation sites that were altered by lithium treatment are shown in Fig. 5. A major finding was that active site phosphorylation in both ERK1 (MAPK3) and ERK2 (MAPK1) was increased in response to lithium, consistent with an increase in ERK activity, i.e., opposite to that seen in response to vasopressin. ERK has many targets, but an important set of targets are transcription factors in the ETS family and AP-1 family whose phosphorylation triggers an immediate early response (Frost et al., 1994; Gille et al., 1995; Kelly and Siebenlist, 1995; Murphy and Blenis, 2006). ERK-mediated phosphorylation of ETS and AP-1 transcription factors can lead to a series of transcriptional changes associated with entry into the cell cycle and dedifferentiation, potentially explaining the loss of AQP2 expression. Subsequent studies using RNA-Seq to identify mRNA changes in microdissected collecting ducts from lithium-treated rats were consistent with this model (Sung et al., 2019). One element of the response was a noncanonical activation of NFκB with selective upregulation of Nfkb2 and Relb transcripts. NFκB is a transcription factor ordinarily associated with inflammation, and its activation had previously been shown to repress Aqp2 gene expression (Hasler et al., 2008).
A long-standing hypothesis about the mechanism of lithium inhibition of water transport in the collecting duct is that it inhibits cAMP production (Anger et al., 1990). A detailed transcriptomic analysis using RNA-Seq in collecting ducts microdissected from lithium-treated rats indeed showed a decrease in mRNA coding for the V2 receptor after prolonged exposure, consistent with a role of decreased cAMP production late in the response to lithium (Sung et al., 2019). Thus, based on the combination of phosphoproteomic studies and RNA-Seq studies, the loss of expression of the Aqp2 gene in lithium-induced nephrogenic diabetes insipidus appears to be due to several signaling changes in collecting duct cells. In particular, there is ERK activation triggering an immediate early gene expression pattern followed by an NFκB-dependent inflammatory-like response and, ultimately, decreased V2 receptor expression, all contributing to a decrease in Aqp2 gene transcription.
V2 Receptor Antagonists and Hyponatremic Disorders: How Vaptans Affect the Collecting Duct Phosphoproteome
Water balance abnormalities also include hyponatremic disorders in which patients reabsorb inappropriately large amounts of water from their collecting ducts, resulting in dilution of body fluids and hyponatremia (Knepper et al., 2015). Such disorders are quite common in hospitalized patients, with a prevalence of hyponatremia sometimes greater than 30% in internal medicine services in tertiary care hospitals (Upadhyay et al., 2006). Hyponatremia is sometimes associated with congestive heart failure and hepatic cirrhosis but is most frequently seen with the syndrome of inappropriate antidiuretic hormone hypersecretion. In these syndromes, hyponatremia is the result of nonosmotic release of vasopressin or other antidiuretic substances into the circulation (Knepper et al., 2015). One treatment option is nonpeptide antagonists of vasopressin receptors, the so-called vaptans. Vaptans that antagonize V2 receptors in the kidney have been used clinically to treat chronic hyponatremia (Schrier et al., 2006). In addition, one of the vaptans, namely tolvaptan, has been shown to slow the progression of autosomal dominant polycystic kidney disease (Torres et al., 2012) and has recently been approved by the US Food and Drug Administration for this purpose (Chebib et al., 2018).
Satavaptan was first described as a highly potent and selective V2R antagonist that inhibited vasopressin-stimulated adenylyl cyclase activity and produced a marked aquaretic response in rats (Serradeil-Le Gal et al., 1996). More recently, satavaptan has been demonstrated to act as an inverse agonist of Gsα-mediated vasopressin signaling (i.e., it can decrease adenylyl cyclase activity even in the absence of ligand stimulation) (Azzi et al., 2003) and can act as a partial agonist in the β arrestin pathway (i.e., it can recruit β arrestin and stimulate the MAP kinase pathway) (Azzi et al., 2003).
Hoffert et al. (2014) investigated the effects of satavaptan on vasopressin signaling using quantitative phosphoproteomics in native rat inner medullary collecting ducts. In general, for phosphorylation sites that are altered by vasopressin signaling, satavaptan had opposite effects, consistent with a predominance of antagonism of Gsα-mediated vasopressin signaling (Fig. 6). The authors reported a limited number of phosphorylation sites in which the response to vasopressin was in the same direction as the response to satavaptan and proposed that these sites could be downstream from β arrestin. β Arrestin signaling is generally associated with an activation of the MAP kinase pathway (Daaka et al., 1998). However, active site phosphorylation of ERK1 and ERK2 showed no change in response to satavaptan.
Integration of Phosphoproteomic Results
An important goal in systems biology is to piece together information from various -omics studies to identify pathways and mechanisms involved in biologic responses. Ultimately, the objective is to derive causal models (directed graphs and associated differential equations) that can be used to organize mechanistic knowledge and make predictions. In our studies investigating vasopressin’s actions in the kidney, we want to create causal models that explain each of the vasopressin responses listed in Table 1. This endeavor can be viewed as consisting of two steps: 1) identifying the essential model elements (proteins) that will be “nodes” in the causal models and 2) identifying how these proteins interact (“edges” in the model). The phosphoproteomic studies described in this review focus mainly on the first step, identification of the elements of the model. It should be emphasized that protein phosphorylation and dephosphorylation are only a part of cell signaling, and other types of -omic data can fill in the information gap, including standard proteomics to detect proteins that change in abundance in response to vasopressin (Khositseth et al., 2011), proteomics of subcellular fractions that can identify translocation of proteins within cells in response to vasopressin (Bolger et al., 2012; Schenk et al., 2012; Yang et al., 2015; Pickering et al., 2016), quantification of other post-translational modifications such as acetylation (Hyndman et al., 2018), and the combination of chromatin precipitation coupled to mass spectrometry (Hwang et al., 2017) and chromatin precipitation coupled to DNA sequencing (Jung et al., 2018) that together can identify proteins that bind to chromatin and their sites of binding. The second step, identification of functional interactions between proteins, is more challenging. Often, functional interactions between two proteins are already documented in the literature and can be identified through fundamental tools like PubMed or more sophisticated tools such as STRING (https://string-db.org/) and Ingenuity Pathway Analysis (http://pages.ingenuity.com/rs/ingenuity/images/IPA_data_sheet.pdf), although the supporting data for these interactions are generally in different cell types than the cell of interest and need experimental confirmation. Ultimately, candidate signaling pathways and networks need to be tested experimentally, e.g., through protein overexpression or siRNA knockdowns/CRISPR deletions, etc., as well as time course studies that can test requirements for causality. Formal computational models can be used to predict the responses within complex networks. Ultimately, more formal tests of candidate models can be achieved by analyzing ensembles of phosphoproteomic data sets via structural equation modeling or path analysis (Stein et al., 2017).
Although the emphasis in this paper is on data sets describing responses to vasopressin in collecting duct cells and how vasopressin signaling may be perturbed in water balance disorders, we propose that the data will be useful for modeling PKA-dependent signaling pathways involving other receptors in other cell types. To facilitate the use of the data, all data sets can be browsed and/or downloaded at the Kidney Systems Biology Project website (https://hpcwebapps.cit.nih.gov/ESBL/Database/).
Authorship Contributions
Wrote or contributed to the writing of the manuscript: Salhadar, Matthews, Raghuram, Limbutara, Yang, Datta, Chou, Knepper.
Footnotes
- Received January 23, 2020.
- Accepted March 24, 2020.
↵1 Current affiliation: Yenepoya Research Center, Yenepoya (Deemed to be University), Mangalore, India.
Based on a talk given by M.A.K. preented as part of an ASPET symposium “Phosphoproteomic Analysis of G-Protein–Coupled Pathways” in Orlando Florida on April 7, 2019.
This research was supported by the Intramural Research Program of the National Heart, Lung, and Blood Institute [Projects ZIA-HL001285 and ZIA-HL006129].
Authors declare that there is no conflict of interest.
Abbreviations
- AMPK
- AMP-activated protein kinase
- AQP2
- aquaporin-2
- CDK
- cyclin-dependent kinase
- Dapk
- death-associated protein kinase
- dDAVP
- D-amino D-arginine vasopressin
- EPAC
- exchange protein directly activated by cAMP
- ERK
- extracellular signal–regulated kinase
- GO
- Gene Ontology
- GPCR
- G-protein coupled receptor
- Gsα
- heterotrimeric G-protein stimulatory alpha subunit
- H89
- N-[2-p-bromocinnamylamino-ethyl]-5-isoquinolinesulphonamide
- LC-MS/MS
- liquid chromatography coupled to tandem mass spectrometry
- MAPK
- mitogen-activated protein kinase
- mpkCCD
- a mouse epithelial cell line with characteristics of cortical collecting duct pricipal cells
- NDI
- nephrogenic diabetes insipidus
- PKA
- protein kinase A
- RNA-seq
- RNA-sequencing
- SILAC
- stable isotope labeling of amino acids in cultured cells
- V2R
- type 2 vasopressin receptor
- U.S. Government work not protected by U.S. copyright