{'fraction_non_alphanumeric': 0.03670903757508032, 'fraction_numerical': 0.021343763095404387, 'mean_word_length': 5.208116545265349, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Differential Regulation of GPCRs—Are GRK Expression Levels the Key?

G protein-coupled receptors (GPCRs) comprise the largest family of transmembrane receptors and their signal transduction is tightly regulated by GPCR kinases (GRKs) and β-arrestins. In this review, we discuss novel aspects of the regulatory GRK/β-arrestin system. Therefore, we briefly revise the origin of the “barcode” hypothesis for GPCR/β-arrestin interactions, which states that β-arrestins recognize different receptor phosphorylation states to induce specific functions. We emphasize two important parameters which may influence resulting GPCR phosphorylation patterns: (A) direct GPCR–GRK interactions and (B) tissue-specific expression and availability of GRKs and β-arrestins. In most studies that focus on the molecular mechanisms of GPCR regulation, these expression profiles are underappreciated. Hence we analyzed expression data for GRKs and β-arrestins in 61 tissues annotated in the Human Protein Atlas. We present our analysis in the context of pathophysiological dysregulation of the GPCR/GRK/β-arrestin system. This tissue-specific point of view might be the key to unraveling the individual impact of different GRK isoforms on GPCR regulation.


INTRODUCTION
G protein-coupled receptors (GPCRs) constitute a family of over 800 membrane-localized receptors. They respond to a large variety of extracellular stimuli, among them, photons, odors, hormones, or neurotransmitters, to induce specific intracellular signaling (Marinissen and Gutkind, 2001). This is achieved by a vast diversity of ligand binding domains. Nevertheless, GPCRs share a seven-transmembrane architecture that undergoes large conformational changes during receptor activation in order to activate a common set of intracellular signaling proteins (Nygaard et al., 2013;Latorraca et al., 2017). Hence, G proteins, GPCR kinases (GRKs) and arrestins, as most prominent interaction partners of GPCRs, engage active receptors at their opened intracellular cavity in a similar fashion (Nygaard et al., 2013;Flock et al., 2017). This process usually involves the insertion of small loop structures or alpha-helical domains into the GPCR cavity. The similarities between the C-terminal alpha helix of G α subunits, the N-terminal domain of GRKs, and the finger loop region (FLR) of arrestins, which enable or enhance the interaction with active GPCRs, are highlighted in Figure 1A.
For the main signaling transducers, the trimeric G proteins, this interaction leads to a guanosine diphosphate (GDP)guanosine triphosphate (GTP) exchange followed by dissociation of G α and G βγ subunits (Oldham and Hamm, 2008;Flock et al., 2015). The now activated G protein subunits are able to individually regulate levels of second messengers [e.g., cyclic adenosine monophosphate (cAMP), calcium, or diacylglycerol (DAG)] to induce a cellular response. Subsequently, intracellular peptide stretches of active GPCRs are phosphorylated by GRKs. In turn, this accumulation of negative charges enables high affinity binding of arrestins (Gurevich and Gurevich, 2019), initiating the desensitization and internalization of receptors. As arrestins and G proteins utilize at least overlapping binding interfaces (DeWire et al., 2007), arrestin-bound receptors are canonically unable to further induce their primary signaling. Moreover, arrestins have been shown to serve as scaffolds for more than 100 intracellular proteins (Xiao et al., 2007;Crepieux et al., 2017), that enable the formation of specific effectorhubs, regulating intracellular trafficking and signaling of active GPCRs. In this review, we want to discuss the current state of research regarding the phosphorylation-dependent processes that underlie GPCR regulation. Moreover, we want to highlight the potential influence of tissue specific expression levels of GPCR-regulating genes on signaling outcomes.

ARRESTINS AND GRKs FACILITATE TARGETED DOWNSTREAM FUNCTIONS FOR HUNDREDS OF GPCRs
Human physiology features a sizeable amount of G α and G βγ subunits (Milligan and Kostenis, 2006). Thus, the diversity of primary GPCR signaling is adequately explained as different receptors preferably couple to specific combinations of G protein trimers (Inoue et al., 2019). However, the downregulation of most GPCRs is tightly controlled by only four ubiquitously expressed GRKs (GRK2, 3, 5, and 6) and two arrestin isoforms, namely β-arrestin1 and β-arrestin2. Still, the processes enabled by these proteins are highly diverse and seem specific for each GPCR. For some receptors, the interactions with GRKs and arrestins lead to desensitization and immediate recycling, redirecting the receptor back to the membrane after initial internalization (Claing et al., 2002). In contrast, certain GPCRs exhibit prolonged intracellular trafficking which localizes the receptors to specific intracellular compartments and may give rise to a second wave of endosomal signaling (Godbole et al., 2017).
GRKs have been shown to be allosterically activated via binding to active GPCRs (Palczewski et al., 1991;Chen et al., 1993;Huang and Tesmer, 2011). This binding mechanism has not been fully understood yet, but possibly features the insertion of a N-terminal α-helix into the cytoplasmic cavity of the GPCR. Although structural evidence is not necessarily conclusive (Cato et al., 2021), this mode of GRK-binding is highly attractive, as G proteins and arrestins probe for active GPCR conformations in a similar fashion (Figure 1A). In a cellular context, GRKbinding leads to the phosphorylation of active GPCRs at their intracellular sites. Notably, GRKs have also been shown to phosphorylate non-GPCR substrates (Palczewski et al., 1991;McCarthy and Akhtar, 2002), albeit with higher efficiency in the presence of active GPCRs. Thus, GRKs most likely also regulate other cellular processes in a phosphorylation-dependent manner, but in this review, we will predominantly discuss their impact on GPCR signaling. Non-visual GRKs are classified into two families Mushegian et al., 2012;Homan and Tesmer, 2014). GRK2 and GRK3 constitute the GRK2 family and are expressed in the cytosol. Subsequent to GPCR activation, GRK2 and 3 are recruited to the membrane, facilitated by GPCR complex formation and stabilizing interactions with G βγsubunits (Tesmer et al., 2005). In contrast, GRK4 family kinases, namely GRK4, 5, and 6, are generally membrane-associated. In this review, we will further focus on effects of GRK2, 3, 5, and 6. Some of these GRK isoforms have been shown to preferentially phosphorylate different residues at the intracellular side of GPCRs (Nobles et al., 2011), to induce receptor internalization and desensitization.
Upon binding to the active and phosphorylated GPCR, arrestins undergo conformational changes that involve the disruption of the polar core and three element interaction site, the two main auto-inhibitory intramolecular interactions. This renders the arrestin C-terminus and phosphate-sensing N-domain solvent-exposed, accompanied by a ∼18 • interdomain rotation. Especially, since arrestins have no enzymatic function, these conformational changes can be seen as hallmarks of arrestin activation. The release of the arrestin C-terminus is furthermore hypothesized to play a central role in the mediation of arrestindependent downstream functions. It harbors binding motifs for the adaptor protein 2 (AP2) complex and clathrin (Goodman et al., 1996;Krupnick et al., 1997;Laporte et al., 2000), in addition to a mitogen-activated protein kinase kinase (MEK) phosphorylation site (Cassier et al., 2017) that enables scaffolding of mitogen-activated protein kinases (MAPKs). Hence, arrestins are able to facilitate clathrin-dependent GPCR internalization and enhance G protein-induced MAPK signaling. In recent years arrestins have been shown to assume distinct conformational states, accommodating not only the active structure of a GPCR but also its specific intracellular phosphorylation. Depending on the overall geometry of the resulting GPCR-arrestin complex, a certain set of effector proteins may then be recruited to orchestrate specific functions.
Crystal structures (Shukla et al., 2013;Kang et al., 2015) and cryo-electron microscopy (Thomsen et al., 2016;Huang et al., 2020;Lee et al., 2020;Staus et al., 2020) studies have shown that these GPCR-arrestin complexes can occur in different configurations. Although they are most probably not mutually exclusive but rather present in a certain equilibrium, different GPCRs make use of distinct binding interfaces when coupling to arrestins. The two main interaction sites on the receptor are constituted by the opened intracellular cavity of the active GPCR and phosphorylated peptide stretches like the C-terminus or intracellular loop 3 (IL3). Arrestins bind to these phosphorylated regions via positive charges buried in their N-domain. Subsequently, the active GPCR cavity is engaged by the arrestin FLR, which is inserted into the receptor transmembrane helix bundle and might assume an alpha-helical structure to stabilize this interaction (Kang et al., 2015). GPCRarrestin complexes that make use of both binding interfaces were termed either "core, " "tight, " or "snuggly" and are usually characterized by high affinity binding and uncoupling of G proteins. Recently, GPCR-arrestin complexes were discovered that only rely on the interaction between the arrestin N-domain and the phosphorylated GPCR C-terminus (Thomsen et al., 2016;Nguyen et al., 2019). This complex configuration is independent of the FLR and does not utilize the transmembrane helix bundle binding interface, therefore still allowing further activation of G proteins ( Figure 1B). Moreover, arrestins that associate with GPCRs in this "hanging" configuration can still assume active conformations and have been shown to functionally increase receptor internalization (Kumari et al., 2017). Thus, GRK-mediated receptor phosphorylation is crucial for the formation of "core" and "hanging" GPCR-arrestin complexes. Phosphorylation is often also hypothesized to be the starting point of arrestin complex formation, however, the precise determination of succession of these binding events is still occluded, as arrestins also have an affinity for active, yet unphosphorylated GPCRs (Gurevich and Gurevich, 2006;Haider et al., 2019;Drube et al., 2021).
Differential spacing of negative charges at the receptor C-terminus has been shown to induce specific conformational changes in arrestins (Lee et al., 2016;Nuber et al., 2016;Mayer et al., 2019). Furthermore, these conformational states have been linked with distinct functional outcomes (Yang et al., 2015;Lee et al., 2016). As these findings suggest that every GPCRarrestin complex is formed in a specific configuration, this could explain how only two β-arrestin isoforms are able to mediate targeted processes for more than 800 different GPCRs. Based on this argumentation, the "barcode" hypothesis was put forward, stating that the arrestin N-domain is capable of recognizing a plethora of different GPCR phosphorylation states. Different phosphorylation patterns ("barcodes") would then only induce certain conformational changes that dictate arrestin functions for the interaction with a given GPCR (Figure 2A).

HOW ARRESTINS INTERPRET DIFFERENT PHOSPHORYLATION PATTERNS: THE "BARCODE" HYPOTHESIS
In its most straightforward interpretation, the "barcode" hypothesis states that arrestins react to different phosphorylation patterns via specific conformational changes in order to fulfill targeted functions (Figure 2A). This adequately explains how different GPCRs can experience divergent arrestinmediated regulation, and constitutes a solid foundation for the investigation of these phosphorylation-dependent processes. In line with this hypothesis, arrestins have been shown to undergo specific conformational changes for the coupling with different GPCRs (Lee et al., 2016;Nuber et al., 2016).
Multiple studies showed that different GRK isoforms preferentially phosphorylate specific sites of the same GPCR (Nobles et al., 2011;Doll et al., 2012;Miess et al., 2018). These findings expand the "barcode" hypothesis, as they suggest FIGURE 2 | Differential and GRK-specific GPCR phosphorylation induces distinct β-arrestin mediated functions. (A) The most straightforward way to interpret the "barcode" hypothesis, as different GPCRs feature different C-terminal phosphorylation patterns to induce distinct β-arrestin functions. (B) Individual GRK isoforms or families (GRK2/3 or GRK5/6) have been shown to preferentially phosphorylate specific sites at different GPCR C-termini (Nobles et al., 2011;Doll et al., 2012;Miess et al., 2018). Depending on the availability of kinases in a cellular system, the same GPCR could be phosphorylated by GRK2/3 or GRK5/6 only, to induce specific functions, or by all GRK isoforms to achieve the activation of all possible β-arrestin functions. (C) Certain GPCRs have been shown to be functionally phosphorylated by GRK2/3 only, or GRK2/3/5/6 (Drube et al., 2021). This might constitute another layer of coupling preference at the foundation of the "barcode" hypothesis. (D) Second messenger kinases, like PKC and PKA are activated by the primary G protein signaling and have been shown to phosphorylate GPCRs directly. Additionally, they are able to modulate the activity of certain GRK isoforms or families (Chuang et al., 1995;Winstel et al., 1996;Pronin and Benovic, 1997). that one receptor may feature different phosphorylation states depending on the cellular context and the availability of kinases. For example, the β2 adrenergic receptor (β2ADR) has been shown to be differentially phosphorylated by GRK2 or GRK6, resulting in kinase-specific C-terminal phosphorylation patterns (Nobles et al., 2011). From these results, a "site-specific barcode" hypothesis emerged, which suggests that GRK2/3 or GRK5/6 phosphorylate the receptor at different sites to induce divergent functions ( Figure 2B). Thus, depending on the available kinases, a GPCR could be phosphorylated at GRK2/3-or GRK5/6-specific sites only, or fully phosphorylated by all four GRK isoforms to induce all possible arrestin-mediated functions.
Indeed, there is evidence that supports this hypothesis as specific phosphorylation patterns have been linked with distinct conformational changes in β-arrestins and downstream functions (Yang et al., 2015). Interestingly, GRK2/3 phosphorylation was proposed to be the driver of receptor internalization, whereas GRK5/6-mediated GPCR phosphorylation was linked with increased ERK signaling Ren et al., 2005;Yang et al., 2015). In contrast to these reports, overlapping or even opposing effects for individual GRK isoforms were identified, depending on the used cellular system and the investigated receptor (Tran et al., 2004;Zhu et al., 2013).
The mentioned studies rely on siRNA/shRNA approaches or GRK inhibitors to investigate the impact of individual GRK isoforms on GPCR regulation. These methods bear the risk of co-analyzing a remaining expression of targeted GRK(s) in knockdown approaches, or potential off-target effects of pharmacological intervention. Furthermore, the impact of these methods depends on the initial endogenous GRK expression levels, which were not assessed in these studies. As an example, the knockdown or inhibition of GRK2 would have less pronounced effects in a cellular system that genuinely features a low expression of GRK2. Non-visual GRKs are usually thought of as ubiquitously expressed and their actual tissue distribution is underappreciated in most studies that focus on molecular mechanisms of GPCR regulation. Additionally, no clear consensus sequences have been identified for specific GRK isoforms, although efforts were made to fill this gap (Pinna and Ruzzene, 1996;Asai et al., 2014;Kang et al., 2020).
Recent studies which utilize the CRISPR/Cas9 technology to achieve a partial (Moller et al., 2020) or complete genetic ablation (Drube et al., 2021) of GRK2, 3, 5, and/or 6 suggest that GPCR-specific GRK-coupling preferences might determine which isoforms regulate a given receptor ( Figure 2C). Using β-arrestin recruitment as a read-out for GRK-mediated receptor regulation, two subsets of GPCRs have been identified (Drube et al., 2021): receptors that are functionally phosphorylated by GRK2, 3, 5, and 6 and those for which arrestin recruitment could only be mediated by GRK2 and 3. By analysis of the β2ADR, this study shows that even though GRK2 and GRK6 preferentially phosphorylate distinct C-terminal sites (Nobles et al., 2011), the individual overexpression of either kinase mediates β-arrestin recruitment to the same extent. These findings indicate that different GRK isoforms might be able to induce identical GPCR regulation on a molecular level, but specific contributions to these processes are ultimately defined by the relative tissue expression of GRK2, 3, 5, and 6.
GPCR phosphorylation patterns are also influenced by second messenger kinases like protein kinase A (PKA) or protein kinase C (PKC) (Figure 2D). Those kinases are activated via the primary G s or G q signaling pathways, respectively, and have been shown to phosphorylate GPCRs directly. Interestingly, PKC also phosphorylates GRKs and is able to modulate their activity (Chuang et al., 1995;Winstel et al., 1996;Pronin and Benovic, 1997). Thus, the resulting phosphorylation "barcode" of a GPCR might be changed by direct phosphorylation or via increasing or decreasing the activity of specific GRK isoforms, depending on the individual G protein-coupling preference. This cross-talk between GPCR regulating kinases is largely underappreciated in recent literature and needs more elaboration to complete our understanding of phosphorylation-dependent GPCR regulation.
Additionally, there are more unanswered biological questions at the foundation of the "barcode" hypothesis. Given that GRK isoforms preferably phosphorylate different sites, how is it that specific GRK consensus sequences are still elusive? Can a receptor molecule be phosphorylated by more than one GRK? If so, does the sequence in which a GPCR is phosphorylated by multiple GRK isoforms change the resulting phosphorylation pattern? These questions still need to be answered by future experiments in order to unravel the intricate details of GPCR regulation.

HOW THE "BARCODE" HYPOTHESIS CAN BE INTERPRETED STRUCTURALLY: INTRINSICALLY DISORDERED REGIONS
One possible extension to explain the "barcode" hypothesis structurally, which goes beyond pure electrostatic interactions of negatively charged phosphate groups on the receptor with basic amino acid side chains of arrestin, might be intrinsically disordered regions (IDRs) of the GPCR itself. IDRs are longer protein regions which do not show a persistent traditional secondary structure of an α-helix or β-sheet (van der Lee et al., 2014;Shammas et al., 2016). Such disordered regions are frequently found in proteins which are involved in signaling cascades (Wright and Dyson, 2015). Intriguingly, IDRs can form different secondary structures when interacting with specific binding partners. An impressive example is the protein p53 which was crystallized with 14 different binding partners and depending on the complex partner, the IDRs of p53 exhibited very different structures (Oldfield and Dunker, 2014). The analysis of GPCR sequences identified IDRs with >50 amino acids in three major receptor regions, namely the N-terminus, the third intracellular loop (IL3), and the receptor C-terminus (Jaakola et al., 2005;Venkatakrishnan et al., 2014). Not surprisingly, the two intracellular regions are well known to be involved in the signal transduction of GPCRs. Due to their flexibility, they frequently need to be truncated or substituted to increase receptor stability in structural biology approaches (Fonin et al., 2019). Furthermore, IDRs are frequently subject to posttranslational modifications which help to support structural stabilization of such regions (Venkatakrishnan et al., 2014). Most commonly IDRs are stabilized by phosphorylation, followed by less common ubiquitination (Bah and Forman-Kay, 2016). Both are well known post-translational modifications for GPCRs occurring within IL3 and the C-terminus of the receptor (Patwardhan et al., 2021). Since IDRs are characterized by a lack of persistent structure (Shammas et al., 2016), their folding state may greatly influence the kinetics of interactions with other partners. For example, increasing the proportion of IDRs with a structure that resembles the bound state might enhance the binding affinity for the partner protein (Shammas et al., 2016). This might be due to effects either on the binding on-rate (k on ) or off-rate (k off ) of the complex.
If we now carefully consider these possibilities, we can envision that a given GPCR interacts with a GRK and depending on their relative complex geometry, this event will add the first phosphate group to the receptor stretch which is closest to the active site of the GRK. Early experimental evidence for such a scenario was demonstrated for rhodopsin 30 years ago when even exogenous peptides in the vicinity of GRK1 were phosphorylated (Palczewski et al., 1991). This initial phosphorylation could have local structural consequences and allow or disallow certain residues of the receptor to be phosphorylated next. Depending on the GRK subtype, this can have different consequences for the phosphorylation pattern of a given GPCR. In the case that a GPCR is phosphorylated by more than one GRK, even the relative sequence of GRKs phosphorylating the receptor might have differential consequences. This relative order could be dominated by either different GRK expression levels or accessibility of the GPCR. Such a scenario could help to explain the apparent lack of consensus sequences for GRKs and account for altered GPCR signaling when certain GRKs are up-or downregulated under pathophysiological conditions. GRKs and β-arrestins are often stated to be ubiquitiously expressed (Nogues et al., 2018). However, a detailed comparative analysis of the tissue and cell-type specific expression pattern of β-arrestins or GRKs is currently not available (Nogues et al., 2018). Therefore, to understand the GRK/arrestin regulatory system in more detail, we analyzed the Human Protein Atlas (HPA) 1 (Thul et al., 2017) for reported expression levels of GRK2, 3, 5, 6, and the two β-arrestins. Furthermore, we included five human arrestin-domain-containing (ARRDC) proteins, also called α-arrestins, based on similarities in mechanistic substrate recognition (Aubry et al., 2009;Kang et al., 2014): Their yeast homolog proteins named ART (arrestin related trafficking adaptors as synonym for ARRDCs in yeast) are reported to use a basic patch in their arrestin domain to recognize the exposed acidic sorting motive of their substrate, for instance a nutrient transporter. To be recognized, the transporter must exist in a conformation that exposes the acidic sorting motif. This exposure occurs during the substrate transport process (active protein state) and is further assisted by phosphorylation (Kahlhofer et al., 2021). Interestingly, ARRDCs lack the auto-inhibitory polar core region seen in visual and β-arrestins and might therefore resemble more of an active arrestin state. Although little is currently known on the function of human ARRDCs, these proteins were reported to interact with GPCRs (Tian et al., 2016). In combination with the mechanistic similarities from their yeast homologs, this observation encouraged us to assemble the information on ARRDC expression besides the β-arrestins.

THE TISSUE PERSPECTIVE: ARE GRK EXPRESSION LEVELS THE KEY?
To evaluate the composition of GPCR-regulating systems for different tissues, we accessed the HPA and analyzed the relative tissue-specific expression levels for various GRK and arrestin isoforms. The HPA is a largescale project, aiming to elucidate human gene expression and localization in cells, tissues, and organs (Uhlen et al., 2015). Since its first publication in 2005, the website has been updated multiple times to include an 1 http://proteinatlas.org/humanproteome/cell increasing amount of data generated by different techniques and to combine information from various sources. To compare expression levels of the four ubiquitously expressed GRK isoforms, the two β-arrestins, and ARRDC1-5, we utilized the consensus transcriptomics data of the HPA, the Genotype-Tissue Expression (GTEx), and the Functional Annotation of Mammalian Genomes 5 (FANTOM5) project, made available on the HPA website 2 (HPA version 20.1, Ensembl version 92.38, last accessed March 10th, 2021). Although mRNA expression levels do not always equate to protein levels in the cells, we nevertheless assume that the mRNA levels somewhat reflect the resulting protein levels. Therefore, we used the available mRNA expression data for our analysis, as it is more detailed than the existing protein expression data. For each gene, the consensus normalized RNA expression (NX) value is calculated via normalization to the maximum expression value found in the three sources ( Table 1). By comparing the consensus NX values, different expression patterns within distinct tissues can be identified.
With this approach, we found tissues that predominantly express one GRK, with all other isoforms being comparatively lower expressed [e.g., GRK5 in heart muscle (23) or GRK2 in skin (46)]. The database also reveals tissues in which two GRK isoforms are comparably high expressed [e.g., GRK2 and 6 in bone marrow (7) or GRK2 and 5 in gallbladder (21)] or tissues with similar NX values for all GRKs [e.g., smooth muscle (48)]. Interestingly, some functional groups of tissues, categorized according to the HPA, share common expression patterns. For example, GRK2 is the predominant isoform expressed in all assessed tissues of the brain, whereas GRK3 is the most abundant isoform in adipose tissue. Bone marrow and lymphoid tissues feature high expression levels of GRK2 and GRK6. Furthermore, GRK6 is highly expressed in all assessed blood cells. Some of them express GRK6 and GRK2 at similar levels [dendritic cells (14) and total peripheral blood mononuclear cells (PBMC; 59)], or feature GRK2 as the second highest expressed isoform [e.g., granulocytes (22) and monocytes (31)]. In contrast, B-(5) and T-cells (52) show similar expression levels of GRK6 and GRK5.
These different GRK expression patterns occur alongside distinct expression levels of β-arrestins. Some tissues express the two β-arrestin isoforms at similar levels [e.g., colon (12) or lung (28)], while other tissues feature a predominant expression of one isoform [e.g., β-arrestin2 in bone marrow (7) or β-arrestin1 in pancreas (35)]. Considering the expression levels of ARRDC1-5 adds another layer of complexity to this system of GRKmediated GPCR regulation. To visualize the respective proteinspecific expression profiles for all listed tissues in the HPA, we prepared radar plots of β-arrestin1 and 2, GRK2, 3, 5, and 6 (Figure 3). Using a clustering heatmap [generated with R package pheatmap (Kolde, 2013. pheatmap: Pretty Heatmaps. R package version 1.0.12 3 )], we analyzed the relative expression of these genes, normalized to the respective maximal expression (Figure 4A). The clustering algorithm identified the highest degree of similarity for the relative expression profiles of GRK2, 6, and β-arrestin2, according to the Euclidean distance. Following  this analysis, we depicted the relative expression data for these three genes as an overlay radar chart ( Figure 4B). This overlay reveals stunningly similar tissue expression patterns for these three proteins. It is tempting to speculate that GRK2, 6, and β-arrestin2 constitute an intricate system, in which disbalance is unfavorable and might lead to dysfunctional GPCR regulation under pathological conditions.

PATHOPHYSIOLOGICAL EFFECTS OF DYSREGULATED GRK EXPRESSION CHANGES
As every cell of the human body expresses GPCRs, the regulated expression levels of GRKs and β-arrestins are crucial to maintain healthy cellular and organ functions. In the following section, we highlight selected examples where a dysregulation of this delicate regulatory system might contribute to the development or progress of different pathological conditions. The role of dysregulated GRK expression in the development of tumors was subject to extensive work and we refer to excellent reviews for further reading (Nogues et al., 2017;Nogues et al., 2018;Sun et al., 2018;Yu et al., 2018). All non-visual GRKs have been found to be dysregulated in at least one tumor model where they can act either as oncogenes or as tumor suppressors. As one example, the mean mRNA expression level of GRK5 for all analyzed tissues (Figure 3 and Table 1) is 9.7 NX. In brain tissues (indices 3, 6, 9, 10, 13, 24, 25, 30, 33, 39, 49, and 54 of Table 1), the expression levels of GRK5 range from 2.5 to 7.8 NX. Similarly low expression levels are also seen in prostate ( Table 1, index 40) with a relative expression of 7.5 NX. In these naturally low GRK5 expressing tissues, upregulation of GRK5 is increasing aggressiveness of glioma (Kaur et al., 2013) and is associated with increased proliferation of prostate cancer (Kim et al., 2012;Chakraborty et al., 2014;DeRita et al., 2017). In contrast, downregulation of GRK5 expression in colon ( Table 1, index 12), a tissue with a high expression of 18.8 NX, leads to promoted proliferation in colorectal cancer (Wu et al., 2011). GRKs cannot be generally classified as either tumor suppressors or promotors, since their influence on tumor progression is highly specific for individual cancer entities or tested cell lines . It is tempting to speculate that upregulation of GRK levels in tissues that naturally feature a low expression level of that given GRK, or downregulation of GRK levels in high expressing tissues might allow a prediction of the impact on cancer progression. This again strengthens the idea that the balance of different players in the given cellular context is the key for physiological regulation of cell growth.
GRKs are important regulators of cell migration, which is crucial for the formation of metastases. Hence, dysregulated GRK levels influence the migratory potential of cancer cells. Changes in GRK2 expression lead to different outcomes depending on the used stimuli and cell type and were extensively discussed elsewhere (Penela et al., 2014). Again, a general association of up-or downregulation of GRK2 with reduction or promotion of migration cannot be made. GRK3 regulates CXCR4-mediated migration and metastasis in breast cancer cell models (Billard et al., 2016). It was shown that shRNA mediated downregulation of GRK3 in breast cancer cell lines led to an increased migration toward CXCL12, whereas overexpression of GRK3 diminished the chemotaxis.
A study using GRK6 knockout (GRK6 −/− ) mice showed that the absence of GRK6 led to increased growth of subcutaneously injected Lewis lung cancer cells, and an increased formation of metastases formed by tail vain injected Lewis lung cancer cells (Raghuwanshi et al., 2013). In this model, CXCR2-mediated promotion of metastasis is regulated by GRK6, and the loss of this negative regulator promotes the malignant phenotype.
Besides the involvement of GRKs in cancer biology, the role of GRK2 in the cardiovascular system is also well studied Schumacher and Koch, 2017;Murga et al., 2019). The importance of GRK2 in the heart is highlighted by the fact, that homozygous GRK2 −/− mouse embryos exhibited a more than 70% decreased cardiac ejection fraction (Jaber et al., 1996), whereas heterozygous GRK2 +/− mice showed increased contractile function compared to wild type mice . This again indicates that the balanced expression is important for physiological function of GRKs and that any change in this delicate system often lead to unpredictable outcomes.
Besides adaptive dysregulation by pathophysiological conditions, gene mutations can also lead to altered expression levels. Mutations in GRK2 were detected in patients suffering FIGURE 3 | Visualization of the tissue-specific expression levels of GRK2, 3, 5, and 6, β-arrestin1 and -2. The consensus expression data of  from Jeune syndrome (Bosakova et al., 2020). In one patient a mutation was identified to cause a functional loss of GRK2. Interestingly, this did not lead to expected embryonic lethality as seen in mice (Jaber et al., 1996), as the patient was born alive, but passed away 5 days after birth. GRK2 was identified as an essential regulator of skeletogenesis (Bosakova et al., 2020). The patient had a very small chest and suffered from pulmonary insufficiency, but did not show gross abnormalities in the central nervous system. Functional analyses in the same study revealed an impairment of Hedgehog and canonical Wnt signaling leading to the observed phenotype.
All described examples so far pointed out the importance to maintain physiological GRK expression levels to prevent pathophysiological conditions. Although this is not the immediate focus of this review, the GPCR-GRK-β-arrestin system is also influenced by changes of β-arrestin expression levels. In 60% of patients suffering from Sezary Syndrome (a rare cutaneous T cell lymphoma), a mono-allelic loss of the β-arrestin2 gene was found (Cristofoletti et al., 2019). Cell culture experiments showed that downregulation of β-arrestin2 led to an impaired internalization of CXCR4 after CXCL12 stimulation, and it was hypothesized that this would lead to an increased migration toward high CXCL12 levels in skin. Another study found that, β-arrestin2 deficiency in dendritic cells promotes migration and cytokine production which contributes to autoimmune encephalomyelitis (Cai et al., 2019). The dysregulated expression of β-arrestin1 was found to be important in context of maternal-fetal tolerance in human pregnancies (Liu et al., 2021) where a strongly reduced mRNA expression of β-arrestin1 was found in villous samples of missed abortion.

CONCLUDING REMARKS
Taken together, there is surmounting evidence that the expression levels of GRKs, arrestins, and GPCRs play a crucial role in the development of pathological conditions. Literature suggests that the regulatory system of GPCRs is a common, yet finetuned machinery which is vital for the maintenance of healthy cellular functions. As different tissues express specific sets of GPCRs to properly react to extracellular stimuli, this regulatory system is adjusted via differential expression of GRKs and arrestins to service this exact set of GPCRs. Disturbance of this equilibrated regulation can then have differential consequences, especially considering that malignancies can also feature the overexpression or downregulation of GPCRs. This is highlighted by the seemingly unpredictable behavior of key players, as in cancer, they can act as both, tumor suppressors or oncogenes, depending on the pathological and cellular context. More work has to be done on mapping functional sets of GPCRs expressed by a given cell and understanding the individual impact of different GRK isoforms on their regulation. This tissue-specific point of view, in combination with further development and elaboration of the "barcode" hypothesis might be the key to unraveling the intricate details of GPCR regulation.

AUTHOR CONTRIBUTIONS
EM compiled and visualized the tissue expression data. RH and EM illustrated all figures. EM, RH, CH, and JD wrote the manuscript. All authors contributed to the article and approved the submitted version.