Abstract
Activation of κ opioid receptors (KORs) produces analgesia and aversion via distinct intracellular signaling pathways, but whether G protein-biased KOR agonists can be designed to have clinical utility will depend on a better understanding of the signaling mechanisms involved. We found that KOR activation produced conditioned place aversion and potentiated CPP for cocaine in male and female C57BL/6N mice. Consistent with this, males and females both showed arrestin-mediated increases in phospho-p38 MAPK following KOR activation. Unlike in males, however, KOR activation had inconsistent analgesic effects in females and KOR increased Gβγ-mediated ERK phosphorylation in males, but not females. KOR desensitization was not responsible for the lack of response in females because neither Grk3 nor Pdyn gene knock-out enhanced analgesia. Instead, responsiveness was estrous cycle dependent because KOR analgesia was evident during low estrogen phases of the cycle and in ovariectomized (OVX) females. Estradiol treatment of OVX females suppressed KOR-mediated analgesia, demonstrating that estradiol was sufficient to blunt Gβγ-mediated KOR signals. G protein-coupled receptor kinase 2 (GRK2) is known to regulate ERK activation, and we found that the inhibitory, phosphorylated form of GRK2 was significantly higher in intact females. GRK2/3 inhibition by CMPD101 increased KOR stimulation of phospho-ERK in females, decreased sex differences in KOR-mediated inhibition of dopamine release, and enhanced mu opioid receptor and KOR-mediated analgesia in females. In OVX females, estradiol increased the association between GRK2 and Gβγ. These studies suggest that estradiol, through increased phosphorylation of GRK2 and possible sequestration of Gβγ by GRK2, blunts G protein-mediated signals.
SIGNIFICANCE STATEMENT Chronic pain disorders are more prevalent in females than males, but opioid receptor agonists show inconsistent analgesic efficacy in females. κ opioid receptor (KOR) agonists have been tested in clinical trials for treating pain disorders based on their analgesic properties and low addictive potential. However, the molecular mechanisms underlying sex differences in KOR actions were previously unknown. Our studies identify an intracellular mechanism involving estradiol regulation of G protein-coupled receptor kinase 2 that is responsible for sexually dimorphic analgesic responses following opioid receptor activation. Understanding this mechanism will be critical for developing effective nonaddictive opioid analgesics for use in women and characterizing sexually dimorphic effects in other inhibitory G protein-coupled receptor signaling responses.
Introduction
The need for safer treatments for chronic pain is clear from the current crisis in opioid abuse and the growing number of overdose deaths (Volkow and Collins, 2017). G-biased κ opioid receptor (KOR) agonists are potentially safer analgesics because they lack the addictive properties of mu opioid receptor (MOR) agonists and lack the dysphoric effects of nonbiased KOR agonists (Pfeiffer et al., 1986; Shippenberg and Herz, 1986; Bruchas and Chavkin, 2010; Brust et al., 2016; Chavkin and Koob, 2016). The analgesic effects of KOR agonists occur via G protein-mediated signaling (Land et al., 2009), whereas the dysphoric properties of KOR agonists are dependent on G protein-coupled receptor kinase 3 (GRK3)/p38 MAPK and β-arrestin activation (Land et al., 2009; Ehrich et al., 2015). This dissociation between the analgesic and dysphoric properties of KOR actions has stimulated the development of G protein-biased KOR agonists (Brust et al., 2016) and some are in clinical trials for pain and itch disorders (Eisenach et al., 2003; Delvaux et al., 2004).
Chronic pain disorders are more prevalent in women (Berkley, 1997), but KOR agonists have low or inconsistent analgesic efficacy in female humans and rodents (Mogil et al., 2003; Stoffel et al., 2005; Craft, 2008). Female C57BL/6J and DBA/2J mice are significantly less sensitive to the analgesic effects of KOR agonists compared with males (Mogil et al., 2003). Dysphoric effects of prototypical KOR agonists (i.e., U50488) are present in male rodents (Chartoff and Mavrikaki, 2015), but females have dose- and species-dependent differences in reactivity to the anhedonic and aversive effects of KOR agonists. Female rats show decreased KOR agonist-mediated suppression of intracranial self-stimulation (Russell et al., 2014), whereas female California mice show conditioned place aversion to a lower dose of U50488 compared with males (Robles et al., 2014). These affective changes produced by KOR agonists are hypothesized to be estrous cycle independent (Russell et al., 2014), but KOR-mediated analgesia has been shown to be modulated by estradiol (Mogil et al., 2003). The present study was designed to identify the intracellular signaling mechanisms responsible for sexually dimorphic analgesic responses to opioid activation with the goal of providing insights to guide the development of better analgesics.
Here, we report that KOR activation produced conditioned place aversion and potentiated cocaine conditioned place preference (CPP) in both male and female C57BL/6N mice, but analgesic effects of KOR activation were attenuated in female mice. The analgesic effects of KOR activation were estradiol sensitive in female mice. We then demonstrated that activation of the estradiol-sensitive GRK2 in females increased sequestration of Gβγ and reduced effects of KOR activation on intracellular signaling, inhibition of dopamine release, and agonist-mediated analgesia. Together, these studies identify an estradiol-mediated intracellular signaling mechanism leading to sexually dimorphic responses that may generalize to other G protein-coupled receptors (GPCRs) acting through Gαi/o protein signaling mechanisms.
Materials and Methods
Subjects
Male and female C57BL/6N (B6) mice (n = 674) ranging from 8–16 weeks of age were used in these experiments. Breeding stocks were maintained in the University of Washington Vivarium and some additional WT animals were obtained from Charles River Laboratories. All experimental procedures were approved by the University of Washington Institutional Animal Use and Care Committee and were conducted in accordance with National Institutes of Health's “Principles of Laboratory Animal Care” (NIH Publication No. 86–23, revised 1985). All testing was during the light phase of the 12 h light/dark cycle. Prodynorphin (Sharifi et al., 2001) or GRK3 (Peppel et al., 1997) heterozygotes were bred to generate GRK3 (Grk3−/−) or PDYN (Pdyn−/−) knock-out mice and WT littermates. Mice were genotyped as described previously (McLaughlin et al., 2004; Terman et al., 2004). GRK3 or dynorphin knock-out mice were phenotypically indistinguishable from WT littermates in general locomotor behavior, weight, and lifespan. For estrous cycle determination, vaginal lavage was administered following behavioral testing and cells were placed on glass slides for cytology assessments (McLean et al., 2012).
Drugs
U50488 (2.5 or 10 mg/kg), norbinaltorphimine (10 mg/kg), morphine (3 mg/kg), cocaine (7.5 or 15 mg/kg), and 5′GNTI (0.03 mg/kg) were provided by the National Institute of Drug Abuse (NIDA) Drug Supply Program (Bethesda, MD) and were dissolved in saline for intraperitoneal administration in a volume of 10 ml/kg. U69593 was purchased from Santa Cruz Biotechnology. Nalfurafine (50 μg/kg; NIDA Drug Supply Program) was dissolved in saline and administered subcutaneously in a volume of 10 ml/kg. The GRK2/3 inhibitor CMPD101 (Tocris Bioscience) was dissolved in 10% ethanol/10% Cremaphor EL (Sigma-Aldrich)/80% saline (15 mg/kg) and administered intraperitoneally in a volume of 10 ml/kg. Estradiol (50 μg/kg; Cayman Chemical) was dissolved in 0.1% ethanol/0.1% Cremaphor EL/99% saline and administered intraperitoneally in a volume of 10 ml/kg.
Procedures
Conditioned place aversion/preference.
Mice were tested for aversion to U50488, a KOR selective agonist (Von Voigtlander and Lewis, 1982), in a two-chamber apparatus with distinct visual and tactile cues as described previously (Ehrich et al., 2015). All conditioning and testing sessions lasted 30 min and were recorded on video for analysis in Ethovision version 3.0 (Noldus). On day 1 (pretest), mice freely explored each side of the apparatus. Total time on each side was calculated and mice were then conditioned with U50488 (2.5 mg/kg) paired with their preferred side on subsequent days. On days 2 and 3 (conditioning), mice were confined to one side with saline treatment and, >4 h later, confined to the other side following U50488 administration. On day 4, mice were allowed to freely explore each side of the apparatus and time spent on the drug-paired floor during the test was measured. To test cocaine CPP, a three-chamber apparatus was used (McLaughlin et al., 2003). Mice were given a pretest on day 1 and then cocaine was paired with the less preferred side during conditioning on days 2 and 3. For repeated forced swim stress (rFSS), mice were given a 15 min swim stress (30°C) following the pretest session and four 6-min swims (separated by 6 min each) <10 min before the first cocaine conditioning session, as described previously (McLaughlin et al., 2003). During day 4 (posttest), mice were allowed to freely explore the apparatus. Preference score was determined by subtracting time on the drug-paired compartment during posttest from time on the drug-paired compartment during pretest (post − pre).
Warm water tail withdrawal test.
Mice were tested for latency (in seconds) to withdraw tail (flick) from 52.5°C water to assess anti-nociceptive effects of KOR activation (Bruchas et al., 2007a). Mice were tested for basal latency to flick and treated with U50488 (10 mg/kg), morphine (3 mg/kg), or saline and retested for latency to flick 30 min later. Doses for U50488 analgesia were selected based on Schattauer et al. (2017). In males, 1.0–7.5 mg/kg U50488 did not significantly change tail flick latency. Change in latency to flick from the pretest to the posttest was calculated. For experiments with the selective GRK2/3 inhibitor CMPD101 (Thal et al., 2011), mice were tested for basal latency, treated with CMPD101 (1.5 mg/kg), tested 30 min later, treated with U50488 (10 mg/kg) or morphine (3 mg/kg), and retested 30 min later. Pregnant mice were tested with U50488 before placement with a male and then retested with U50488 at 6 and 13 d following cohabitation. Investigators were blinded to treatment groups during testing.
Hot plate response test.
Mice were tested for latency to paw lick or jump following placement on the hotplate (57.5°C) immediately before treatment with saline or U50488 (10 mg/kg) and mice were retested 30 min later, as described previously (Terman et al., 2004). Changes in latency to respond were calculated as the difference between pretest and posttest.
Pruritis.
Antipruritic activity of nalfurafine was assayed as described previously (Schattauer et al., 2017). Mice were placed individually in observation boxes with a metal grid floor and were given 1 h to acclimate. Mice received either saline or nalfurafine (50 μg/kg, i.p.). Twenty minutes later, mice were injected subcutaneously on the midline of the back of the neck with saline or 5′-GNTI (0.03 mg/kg) and behavior was recorded on video for 15 min. Hind leg scratches directed to the back of the neck were counted. Data analyses were done by an investigator blinded to sex and drug treatment.
Ovariectomy.
Ovaries were removed from female mice (OVX) under isoflurane anesthesia (Smith et al., 2005) and complete removal was visually confirmed. Using mass spectrometry in serum from a subset of OVX mice, we confirmed that circulating estradiol was absent following ovariectomy. Mice were allowed 3 weeks for recovery from surgery before behavioral or biochemical experiments. Sham controls receiving surgery without ovary removal were not used in these studies because baseline tail flick scores in ovariectomized females were not different from WT females. For estradiol replacement experiments, mice received a priming injection of 50 μg/kg estradiol 3 d before testing and then a second injection of 50 μg/kg estradiol 60 min before behavioral testing or tissue collection (Green et al., 1970).
Isolation of proteins.
Mice were decapitated at the specified times following drug administration. Ventral striata, periaqueductal gray (midbrain), and spinal cord were dissected and homogenized with a Dounce homogenizer on ice in MAPK lysis buffer composed of 50 mm Tris-HCl, 300 mm NaCl, 1 mm EDTA, 1 mm Na3VO4, 1 mm NaF, 10% glycerol, protease inhibitor cocktail (Sigma-Aldrich), and phosphatase inhibitor cocktail (Millipore). For analysis of GRK2 protein levels, lysis buffer was supplemented with 1% Triton X-100. Samples were centrifuged 25 min at 30,000 × g at 4°C. Protein concentrations in supernatant were determined by bicinchoninic colorimetric assay (Thermo Fisher Scientific), mixed with 4× LDS buffer (10% lithium dodecyl sulfate, 25% glycerol, pH 8.4; Thermo Fisher Scientific) supplemented with 8% β-mercaptoethanol and stored at −20°C until analysis by western blot.
Western blot analysis.
Samples were loaded with 30 μg of protein (phospho-ERK1/2) or 40 μg of protein (phospho-p38 MAPK) onto 10% Bis-Tris precast gels (Thermo Fisher Scientific) and run at 100 V for 2 h. Blots were transferred to nitrocellulose (Thermo Fisher Scientific) for 1.5 h at 30 V. The nitrocellulose was blocked with 5% BSA-TBS 1 h at room temperature and stained overnight at 4°C for phospho-ERK1/2 (Cell Signaling Technology #cs9101, RRID:AB_331646, lot 28, 1:1000) and ERK2 (Santa Cruz Biotechnology #sc1647, RRID:AB_627547, 1:2000), phospho-p38 (Cell Signaling Technology #cs9211, RRID:AB_331641, lot 23, 1:1000) and p38α (Cell Signaling Technology #cs9228, RRID:AB_490886, lot 2, 1:1000), or phospho-GRK2(S670) (Sigma-Aldrich #MABS155, lot 2792488, 1:250) and GRK2 (Sigma-Aldrich #SAB4500592, lot 218278, 1:500) in 5% BSA-TBS. Blots were incubated in IRdye secondary (Li-Cor Biosciences #926-68070 lot C50721-05 and #926-32211 lot C506602-05, 1:10,000) in 1:1 Odyssey buffer (Li-Cor) and 5% milk-TBS 1 h at room temperature and then scanned on the Odyssey Infrared Imaging System (Li-Cor Biosciences). Band intensity was measured using the Odyssey software and expressed as phospho-immunoreactive protein intensity over the homologous protein band intensity. For total protein quantification, protein was normalized to Revert protein stain (Thermo Fisher Scientific). Data were normalized to percentage of control sample (saline, 100%). The Gβ antibody detected a single band at 37 kDa. ERK1/2 and phospho-ERK1/2 antibodies each detected two bands at 42 and 44 kDa. The antibody for p38α identified a single band at the predicted weight of ∼42 kDa. In addition to the expected band, the phospho-p38 antibody also detected nonspecific bands at ∼45 kDa; only the band overlapping with the p38α band was quantified (Bruchas et al., 2007b). The GRK2 antibody detected bands at the predicted weight of ∼80 kDa and at ∼55 kDa, whereas the phospho-GRK2 antibody detected multiple nonspecific bands, but a single band at ∼80 kDa; the ∼80 kDa band was identified as the specific band by GRK2 overexpression in HEK293 cells combined with phorbol-12-myristate-13-acetate treatments (Pitcher et al., 1999).
Coimmunoprecipitation.
Mice were decapitated at the specified times following drug administration. Ventral striata and spinal cord were dissected and homogenized with a Dounce homogenizer on ice in coimmunoprecipitation buffer (20 mm Tris, 135 mm NaCl, 0.5 mm EDTA, 10 mm sodium fluoride, 1 mm sodium orthovanadate, 1% Triton X-100, 10% glycerol, protease inhibitor cocktail (Sigma-Aldrich), and phosphatase inhibitor cocktail (Millipore), pH 7.5. Next, 700 μg of protein was incubated overnight at 4°C with 5 μl GRK2 antibody (Cell Signaling Technology #cs3982, RRID:AB_330680, lot 1) before incubation for 2 h with 25 μl washed magnetic protein A/G beads (Thermo Fisher Scientific) at 4°C. Samples were washed twice with buffer before elution at room temperature with 40 μl of 4× LDS. Samples were loaded onto 4–12% Bis-Tris precast gels (Thermo Fisher Scientific) and analyzed as described above, blotting for GRK2 (as above) and Gβ (Santa Cruz Biotechnology #sc-378, RRID:AB_2159943, lot C0515, 1:500) band intensity was normalized to GRK2 band intensity.
qRT-PCR analysis of KOR and GRK2 mRNA levels.
Mice were decapitated and ventral striata were dissected on ice. RNA was isolated using the Qiashredder and RNeasy Plus Mini Kits (Qiagen). RNA was amplified on a Stratagene Mx3000P (Agilent) using the Thermo Power SYBR Green RNA-to-CT 1-Step Kit (Thermo Fisher Scientific) according to the manufacturer's instructions, with the standard cycle settings: 50°C 30 min, 95°C 5 min, 40× (95°C 15 s, 62°C 1 min, followed by a dissociation curve ramping from 95°C to 55°C and back to 95°C. For KOR, the following primers were used: GTGGGCTTAGTGGGCAATTCT, AGTGGTAGTAACCAAAGCATCTG. For GRK2, the following primers were used: GCGCCAGCAAGAAGATCCT, GCAGAAGTCCCGGAAAAGCA. PPIA was selected as a reference gene after experiments demonstrated no difference in expression between male and female mice using the following primers: CACCGTGTTCTTCGACATCA and CAGTGCTCAGAGCTCGAAAGT. All samples were run in duplicate (20 ng of RNA) with a standard curve on each plate. Data were analyzed with the MxPro version 3.20 software (Agilent Technologies). Transcript levels were calculated from the standard curve, with KOR and GRK2 transcript levels then normalized to PPIA transcript levels.
Voltammetry.
Mice were decapitated and the head placed in preoxygenated, cold modified artificial CSF (aCSF) (248 mm sucrose was substituted for NaCl). The brain was removed and 250 μm coronal slices containing nucleus accumbens were prepared. To measure evoked dopamine release, nucleus accumbens slices were transferred to standard oxygenated aCSF and incubated for 1 h at 37°C before holding at room temperature before recordings. Slices were placed in a recording chamber and perfused with oxygenated aCSF (124 mm NaCl, 2.5 mm KCl, 1.25 mm NaH2PO4, 1.25 mm MgSO4, 2 mm CaCl2, 10 mm dextrose, 25 mm NaHCO3) at 31–33°C throughout recording. Carbon-fiber working electrodes were fabricated as described previously (Clark et al., 2010; Lemos et al., 2012) hand-cut to 150–200 μm past the capillary tip. A single electric pulse (120–200 μA, 2 ms per phase) was applied to a parallel bipolar stimulating electrode to evoke dopamine release and the potential at the working electrode was held at −0.4 V versus Ag/AgCl and ramped to +1.3 V and back to −0.4V at a rate of 10 Hz. Baseline dopamine release was measured every 2 min for 20 min before measuring release after stepwise increasing concentrations of U69593 (20 min per concentration). For CMPD101 experiments, aCSF with 30 μm CMPD101 was washed onto the slice for 10 min before washing on both 30 μm CMPD101 and increasing doses of U69593. Waveform generation, data acquisition, and analysis were performed using two PCI multifunction data acquisition cards and software written in LabVIEW version 7.1 (National Instruments) as described previously (Clark et al., 2010; Lemos et al., 2012). Responses in the last 8 min (4 recordings) were used for generating U69593 concentration curves in GraphPad Prism software and fitted using a three-parameter least-squares nonlinear regression with the top constrained to 100% and the Hill slope constrained to −1.
Experimental design and statistical analysis
Animal numbers were as follows: conditioned place aversion, n = 13–16 per group; CPP, n = 7–19 per group; tail flick, n = 3–6 per group for males and n = 7–19 per group for females; time course for tail flick, n = 6; hot plate, n = 6–9 per group; pruritis, n = 10–15 per group; dynorphin/GRK3 knock-out tail flick, n = 5–6 per group; OVX/estradiol replacement tail flick, n = 8–9 per group; striatal pERK, n = 8 per group; striatal p38, n = 8–11 per group; midbrain pERK, n = 4–5 per group; striatal GRK2 expression, n = 8 per group; phospho-GRK2, n = 11–13 per group; spinal cord pERK/p38, n = 7–9; GRK2/Gβ coimmunoprecipitation, n = 8–10 per group; CMPD101/phospho-ERK, n = 3–4 per group; dopamine voltammetry, n = 5–9 per group; CMPD101/U50488 tail flick, n = 13; CMPD101/Grk3−/−, n = 9–14 per group; and CMPD101/morphine tail flick, n = 8 per group. Data are expressed as mean ± SEM and were analyzed with Prism 7 software (GraphPad). For EC50, 95% confidence intervals (95% CIs) are reported. Group differences were determined using t tests, ANOVA, or repeated-measures ANOVA as described in the Results. Post hoc comparisons were analyzed with Dunnett's or Sidak's test. For all statistical tests, α was set to 0.05.
Results
KOR activation in females produces aversion but not analgesia
In male B6 mice, KOR activation stimulates Gαi/Gβγ signaling to produce analgesia and decrease itch, whereas arrestin-dependent signaling produces aversion and increases drug-seeking behaviors (Ehrich et al., 2015; Schattauer et al., 2017). We first tested for sex differences in arrestin-mediated behaviors following KOR activation. U50488, an unbiased KOR agonist (Von Voigtlander and Lewis, 1982; Schattauer et al., 2012), produced conditioned place aversion in male and female B6 mice (2.5 mg/kg; Fig. 1A). There was a main effect of drug (F(1,52) = 12.32, p = 0.009) and no main effect of sex or interaction between sex and drug. Planned comparisons between the saline- and U50488-treated mice showed a significant effect of U50488 in males (p = 0.027) and females (p = 0.02). As previously reported, females were more sensitive to the CPP effects of cocaine than males (Russo et al., 2003); significant cocaine CPP was evident in female mice at 7.5 mg/kg, whereas 15 mg/kg was required for males (Fig. 1B). This was confirmed with a one-sample t test showing that 7.5 mg/kg cocaine in unstressed males was not significantly different from a preference score of 0, whereas unstressed females (7.5 or 15 mg/kg; p < 0.0001) and unstressed males (15 mg/kg; p = 0.009) were significantly different from 0. rFSS has been shown to produce a KOR-dependent increase in cocaine-seeking behaviors (McLaughlin et al., 2003) and these effects are GRK3/arrestin dependent (Schindler et al., 2012). Stress potentiated cocaine CPP in both males and females (Fig. 1B). A 3-way ANOVA showed that there was no main effect of sex, but there was a main effect of cocaine dose (F(1,115) = 79.9, p < 0.001) and stress (F(1,115) = 7.14, p = 0.009). There was also a significant interaction between sex and dose (F(1,115) = 33.9, p < 0.001), with both males (p < 0.001) and females (p = 0.011) showing significant differences in preference scores between 15 mg/kg and 7.5 mg/kg of cocaine. These conditioned place aversion and stress-potentiated cocaine CPP studies indicated that arrestin-mediated behavioral effects of KOR activation were present in both males and females.
We then tested Gβγ-mediated behaviors and found that KOR activation produced significant antinociceptive responses in the tail flick assay in males, but not females. In the tail flick assay (Fig. 1C), there was a significant interaction between sex and drug treatment (F(1,43) = 4.71, p = 0.036). Post hoc analyses showed that males treated with KOR agonist had a significant increase in latency to flick compared with saline-treated males (p = 0.008) and were significantly different from females treated with KOR agonist (p = 0.049). Repeated testing of tail flick analgesia in females at 10 min intervals for 60 min following U50488 administration did not reveal any significant effect of drug treatment (Fig. 1D). Higher doses of U50488 were not tested in females because doses of up to 30 mg/kg do not produce significant analgesia in female C57BL/6 mice (Mogil et al., 2003). To interrogate alternative antinociceptive circuits, we next used the hot plate assay and confirmed that KOR-mediated analgesia was also absent in female mice in this assay of analgesia (Fig. 1E). There was a significant interaction between drug treatment and sex (F(1,26) = 12.2, p = 0.002) on hot plate response latency. Males treated with U50488 (10 mg/kg) had significantly different response latencies from saline-treated males (p = 0.0002) and females treated with U50488 (p = 0.001) in the hot plate assay.
To determine whether females showed antipruritic (anti-itch) effects of KOR activation, another hypothesized G protein-mediated behavior (Schattauer et al., 2017), we tested the G-biased KOR agonist nalfurafine (50 μg/kg) in males and females following application of 5′-GNTI, a KOR antagonist that induces compulsive scratching behaviors (Inan et al., 2011). We found that 5′-GNTI produced equivalent baseline levels of scratching behaviors in males and females and that significant antipruritic effects of nalfurafine were observed in both male and female mice (Fig. 1F). There was a significant interaction between sex and nalfurafine treatment (F(1,46) = 4.05, p = 0.05). Post hoc analyses showed that nalfurafine significantly decreased scratching in males (p < 0.0001) and females (p = 0.0005) compared with vehicle controls. These studies demonstrate that KOR produces arrestin-dependent behavioral effects in male and female B6 mice, but some G protein-mediated effects are differentially regulated in males and females.
Lack of antinociceptive response to KOR activation in females is estrogen sensitive and is not attributable to KOR desensitization
We tested whether KOR activation in females had greater intrinsic bias toward arrestin-mediated signaling by global deletion of the GRK3 (Bruchas et al., 2006). We hypothesized that deletion of arrestin-mediated KOR signals would uncover analgesic responses if basal receptor desensitization was elevated in females. However, GRK3 deletion (Grk3−/−) did not alter analgesic responses in female mice (Fig. 2A). We then tested whether altered basal dynorphin tone may contribute to the sexually dimorphic behavioral effects of KOR activation through receptor desensitization and found that global deletion of dynorphin (Pdyn−/−) did not produce KOR-mediated analgesic responses in female mice (Fig. 2A).
The sensitivity to KOR activation in the tail flick assay was estrous-cycle-dependent (Fig. 2B). Females were divided into estrus (low estrogen) or nonestrus groups and females in estrus showed significantly greater analgesia (t(23) = 2.14, p = 0.043) in response to KOR activation (Fig. 2B). To confirm the necessity of circulating gonadal hormones for KOR-mediated analgesia (Mogil et al., 2003), we assessed analgesia 3 weeks following OVX and found a significant effect of OVX on KOR-mediated analgesia in the tail flick assay (F(3,51) = 5.05, p = 0.004). Compared with intact females treated with KOR agonist, OVX females treated with KOR agonist showed a significant increase in tail flick analgesia (p = 0.009). Systemic administration of 17β-estradiol in OVX females impaired KOR-mediated analgesia compared with OVX females (p = 0.02) (Fig. 2C). There was no significant effect of OVX or 17β-estradiol on baseline latency to respond. These studies demonstrated that estrogen was necessary and sufficient to blunt KOR-mediated analgesia. To assess KOR-mediated analgesia during pregnancy, a period of altered hormone circulation (McCormack and Greenwald, 1974), we tested analgesia at multiple time points immediately before and following pair housing females with males. Thirteen days following pairing with a male, KOR-mediated analgesia was observed in the tail flick assay in female mice. A repeated-measures ANOVA confirmed this observation with a main effect of treatment (F(1.66,19.9) = 6.02, p = 0.01) and post hoc analyses showed that day 13 was significantly different (p = 0.024) from day 0 females (Fig. 2D).
KOR activation produces p38 MAPK phosphorylation in males and females but not ERK 1/2 phosphorylation in females
To determine whether differences in KOR signaling could be due to altered basal expression of KOR, we measured levels of KOR mRNA in males and females using qRT-PCR in ventral striatal extracts. We found no significant differences in male and female KOR mRNA in the striatum (Fig. 3A). Based on our behavioral data and Robles et al. (2014), we hypothesized that KOR activation should produce arrestin-mediated signals in males and females, but G protein-mediated signals should not be present in females. Western blot analysis demonstrated that there was a significant effect of KOR activation on phosphorylation of striatal p38 MAPK (F(1,34) = 12.35, p = 0.001), with no effect of sex and no interaction (Fig. 3B). In contrast, there was no main effect of KOR agonist on striatal ERK1/2 phosphorylation (Fig. 3C), but there was an interaction between sex and KOR agonist (F(3,37) = 4.98, p = 0.005). There was no significant difference in total ERK1/2 (males = 102 ± 9%; females = 98 ± 7% of average expression) or total p38 (males = 112 ± 17%; females = 88 ± 17% of average expression) protein levels. Post hoc analyses showed that ERK1/2 phosphorylation in males treated with U50488 were significantly different from saline-treated males (p = 0.007) and from females treated with KOR agonist (p = 0.0003), whereas U50488 treatment did not significantly increase striatal pERK in females (Fig. 3C). Representative Western blots for Figure 3, A–C, are shown in Figure 3, D and E. Representative blots for Figure 3, G and H, are shown in Figure 3F. Western blot analysis of spinal cord extracts 30 min following systemic U50488 (10 mg/kg) administration showed increased p38 phosphorylation in spinal cord in males and females (F(1,30) = 5.97, p = 0.021), with no significant effect of sex or interaction between sex and treatment (Fig. 3G). U50488 significantly increased pERK in males but not females at 30 min, with a significant interaction between sex and treatment (F(1,29) = 14.38, p = 0.019). Post hoc analysis showed that saline-treated males were significantly different from U50488-treated males (p = 0.038). There was no significant effect on pERK at 15 min in males or females, suggesting an altered time course for KOR signaling in the spinal cord compared with brain tissue. In midbrain samples targeted to include the periaqueductal gray, a region known to mediate opioid analgesia, p38 phosphorylation (n = 7–9) at 30 min was nonsignificantly increased in males (132 ± 22% of average expression) compared with females (104 ± 20% of average expression). Similar to p38 phosphorylation, pERK (n = 4–5) at 15 min was also nonsignificantly increased in males (144 ± 21% of average expression) compared with females (74 ± 20% of average expression) in the midbrain.
GRK2 sequesters Gβ/γ to decrease KOR-mediated ERK phosphorylation and dopamine release inhibition in females
G protein-coupled receptor kinase 2 (GRK2) can be modulated by estradiol via phosphorylation of the S670 GRK2 site (Dominguez et al., 2009) and GRK2 can sequester components of signaling pathways upstream from ERK1/2 phosphorylation, including Gβγ and MEK (Daaka et al., 1997; Pitcher et al., 1999; Raveh et al., 2010). GRK2 activity can blunt KOR-mediated effects in cardiac tissue (Chen et al., 2017), suggesting that there may be direct interactions between GRK2 and KOR in the brain. Using qRT-PCR in striatal extracts, we determined that total GRK2 levels were not significantly different between males and females (Fig. 4A). However, phosphorylation of the estrogen target S670 site of GRK2 was significantly greater in females (t(22) = 2.64, p = 0.034) compared with males (Fig. 4B) and there was an estradiol-sensitive increase in GRK2 association with Gβ in ovariectomized females (Fig. 4C). There was a significant effect of GRK2 coimmunoprecipitation compared with input (F(1,31) = 10.01, p = 0.004) and a significant interaction between sample group and estradiol treatment (F(1,31) = 10.74, p = 0.003). Post hoc analysis showed that estradiol treatment significantly increased (p = 0.005) coimmunoprecipitation between GRK2 and Gβγ in the immunoprecipitate, but not in the input sample. A similar effect was observed in spinal cord (n = 3–5), with an increase in GRK2 association with Gβγ in the immunoprecipitate following estradiol treatment (142 ± 19% of average expression) compared with the input sample treated with estradiol (104 ± 19% of average expression), suggesting that GRK2 modulates KOR signaling in the spinal cord as well.
Gβγ regulates ERK phosphorylation via PI3Kγ and Ras activity (Lopez-Ilasaca et al., 1997; Belcheva et al., 1998) and we hypothesized that inhibition of GRK2 would reveal, G protein-mediated effects of KOR activation in females. We first tested whether the selective GRK2/3 inhibitor CMPD101 (Thal et al., 2011; Lowe et al., 2015) (1.25 mg/kg; Tocris Bioscience) could unmask KOR-mediated ERK phosphorylation in females (Fig. 4D). There was a significant effect of sex (F(1,11) = 8.64, p = 0.014) and an interaction between sex and drug treatment (F(1,11) = 8.64, p = 0.0135). KOR activation in females treated 30 min after CMPD101 significantly increased ERK phosphorylation (p = 0.002) compared with females pretreated with saline. These results suggested that KOR/Gβγ-mediated intracellular signaling effects in the striatum were evident in females following GRK2/3 inhibition.
We next tested whether presynaptic inhibition of transmitter release by KOR activation (Chefer et al., 2005) was different in males and females using dopamine voltammetry in the striatum. Concentration–response curves were generated for the KOR agonist U69593 as described previously (Ehrich et al., 2015). Treatment with U69593 inhibited evoked dopamine release with an EC50 of 17 nm (95% CI: 11–27 nm) males and 63 nm (95% CI: 27–148 nm) in females. U69593 was significantly more potent in males than females across multiple doses (nonlinear regression, F(4,81) = 12.9, p < 0.0001) (Fig. 4E,F). There was a significant difference (F(2,19) = 5.199, p = 0.0158) in the logEC50 between groups and post hoc analysis showed that there was a nonsignificant trend (p = 0.087) toward a difference between males and females. Slices prepared from females incubated with 30 μm CMPD101 showed a significantly (p = 0.0013) lower U69593 logEC50 (EC50: 10 nm; 95% CI: 6–16 nm) compared with female slices in the absence of CMPD101. There was no effect of CMPD101 treatment alone on dopamine release. To further assess dose effects of CMPD101, we examined a U69593 dose close to the EC50 of all three groups at a 100 nm concentration (Fig. 4F) and found that there was a significant effect of group (F(2,19) = 9.69, p = 0.0013). Post hoc analysis confirmed the potency of U69593 was significantly lower in female brain slices than in either slices of male brain or slices of female brain in the presence of CPMD101 (p = 0.011). Together, these studies demonstrated that pharmacological inhibition of GRK2 in females could restore KOR-mediated ERK phosphorylation and dopamine release inhibition that was quantitatively similar to the effects produced in males.
Pharmacological inhibition of GRK2/3 enhances opioid-mediated analgesia in females
Our biochemical and voltammetric data showed that GRK2/3 inhibition with CMPD101 enables Gβγ-mediated KOR intracellular and neural circuit level signals in females, suggesting that KOR-mediated analgesia may also be regulated by estrogen-mediated phosphorylation of GRK2. We tested whether GRK2/3 blockade would unmask KOR analgesia in the warm water tail flick assay. CMPD101 pretreatment did not significantly alter baseline analgesic response; however, KOR activation in CMPD101 pretreated females produced a significant increase in tail flick latency (F(1,12.5) = 31.7, p < 0.0001; post hoc: p = 0.0003). To assess the specificity of CMPD101, which inhibits both GRK2 and GRK3, we determined the necessity of GRK3 for the effects of CMPD101 on KOR-mediated analgesia and there was a significant difference between groups (F(2,45) = 6.31, p = 0.004). U50488 in GRK3 knock-outs and littermate controls did not produce a significant analgesic response, but U50488 did increase tail withdrawal latency in female GRK3 knock-out mice pretreated with CMPD101. The response in the presence of CMPD101 was significantly different from WT controls (p = 0.013) and GRK3 knock-out mice (p = 0.007).
Based on our observations with KOR responses in females, we predicted that analgesic responses to MOR activation with morphine may also be attenuated by GRK2 activity in females. There was a significant interaction between drug treatment and sex (F(1,27) = 5.21, p = 0.031) and post hoc analyses showed that females pretreated with saline showed a significantly different response to morphine compared with females pretreated with CMPD101 (p = 0.022). We observed that blockade of GRK2 enhances KOR- and MOR-mediated analgesic responses in females. These studies demonstrate that estradiol stimulates GRK2 phosphorylation to increase Gβγ sequestration and block G protein-mediated biochemical, circuit, and behavioral effects of opioid receptor activation (for schematic, see Fig. 5D).
Discussion
These studies demonstrate that females show blunted intracellular, neural circuit, and behavioral responses to the G protein-mediated actions of KOR agonists. In particular, elevated estradiol was found to underlie the reduction in KOR function in females. Furthermore, we specifically identified estradiol upregulation of GRK2 activity to mediate the observed sex differences in opioid analgesic responses. GRK2-dependent blunting of the analgesic effects of MOR agonist-induced analgesia was also observed in females, suggesting that estradiol-stimulated GRK2 may regulate other GPCR systems.
Sex differences in KOR function
Whereas the G protein-mediated analgesic effects were subject to estradiol regulation, dysphoric effects of KOR activation were observed in both male and female mice. Although there are cycle-independent effects of KOR activation on aspects of reward-seeking behaviors in females (Russell et al., 2014), we did not observe significant sex effects on conditioned place aversion. These aversive effects are mediated by activation of the p38 MAPK (Bruchas et al., 2007b; Land et al., 2009; Bruchas et al., 2011; Ehrich et al., 2015, Schindler et al., 2012) and we found that p38 MAPK phosphorylation was increased in both B6 males and females following KOR activation. In contrast, G protein-mediated molecular or cellular effects, including ERK1/2 phosphorylation (Belcheva et al., 2005) and dopamine release inhibition (Ehrich et al., 2015), were blunted following KOR activation. The attenuation of G protein-mediated signals through estradiol/GRK2 regulation also decreased analgesic behavioral responses in B6 females following KOR activation. An alternative explanation for sex differences in KOR activation is that there are sex-dependent regional differences in kinase recruitment (Liu et al., 2018); however, this hypothesis cannot account for the CMPD101-sensitive differences in KOR-mediated ERK1/2 activation or inhibition of dopamine release.
Few selective KOR agonists have been tested for sex differences in analgesic efficacy in humans (Coffin et al., 1996, Delvaux et al., 1999), but nonselective MOR/KOR agonists have been reported to produce greater analgesia in women than in men (Gear et al., 1996). The discrepancy between rodent and human studies may be due to innate signaling differences in KOR function between rodents and humans (Schattauer et al., 2012; Broad et al., 2016) or to synergistic interactions between MOR/KOR. In female rats, increased MOR/KOR heterodimerization in the spinal cord during periods of high estrogen or progesterone is associated with increased analgesic responses to intrathecal morphine administration (Chakrabarti et al., 2010; Liu et al., 2011, 2017). In addition to species differences, there may be differences in spinal and supraspinal KOR responses (Lawson et al., 2010) that produce inconsistencies in KOR-mediated analgesia. Reviews of the literature describing analgesic effects of selective KOR agonists in rodents concluded that there were mixed effects of KOR agonists in females (Craft, 2003; Fillingim and Gear, 2004; Rasakham and Liu-Chen, 2011), but the mechanisms underlying this variability were previously unknown.
Estradiol blunts KOR-mediated analgesia
Mogil et al. (2003) showed that ovariectomy restored KOR-mediated analgesia and this effect was blocked by estradiol replacement. In agreement with Mogil et al. (2003), our results show that estradiol blunts KOR-mediated analgesia and suggests that reported inconsistencies in KOR-mediated analgesia in females may be due to differences in estrogen levels between study subjects. Progesterone may also alter KOR function in female rats because systemic progesterone administration in gonadectomized or progesterone receptor antagonism alters analgesic responses to U50488 (Stoffel et al., 2005; Liu et al., 2011). Changes in either or both of these hormones may underlie the analgesic effects we observed in pregnant female mice. The multiple sites of action underlying sex differences in KOR analgesia (Mogil et al., 2003; Chakrabarti et al., 2010) indicates that there are complex interactions among agonist dose, administration route, rodent species, and behavioral assays for pain that could be explored in future studies to determine optimal treatment strategies in humans.
Estradiol can produce complex effects on mood and behavior through estrogen receptor activation or interactions with GPCR signaling (Gillies and MacArthur, 2010; Martinez et al., 2016). Based on the observation that KOR activation did not produce analgesia or ERK phosphorylation, we hypothesized that an estradiol-regulated signaling system that controls G protein-biased signals was likely to be a mechanistic target for sexually dimorphic responses. GRK2 is phosphorylated following estradiol treatment, likely through estrogen receptor α (Dominguez et al., 2009), sequesters Gβγ subunits (Lodowski et al., 2003), and decreases ERK activation (Pitcher et al., 1999). GRK2 activity also desensitizes KOR activity in cardiac tissue (Chen et al., 2017), suggesting that estradiol modulation of GRK2 activity could underlie sex differences in KOR activity. Our results showed that, in females, GRK2 phosphorylation was increased and there was no change in GRK2 mRNA levels. We also observed increased association of GRK2 with Gβγ subunits following estradiol treatment in ovariectomized females, indicating that estradiol increased GRK2 sequestration of Gβγ to blunt G protein-mediated signals.
GRK2 regulates opioid analgesia in females
We found that pharmacological inhibition of GRK2 with CMPD101 attenuated sex differences following KOR activation. CMPD101 is a selective inhibitor for members of the GRK2 subfamily (GRK2 and GRK3), with little activity at GRK1 or GRK5 (Thal et al., 2011). To ensure that the effects that we observed were selective to GRK2, we tested the effect of CMPD101 in female Grk3−/− mice and found that the CMPD101-mediated increase in KOR analgesia was not GRK3 dependent. Global genetic deletion of GRK2 in mice causes death during embryonic development (Jaber et al., 1996), preventing an assessment of CMPD101 in GRK2 knock-out mice. In addition to sequestration of Gβγ, GRK2 can also interact with MEK to inhibit phosphorylation of ERK, suggesting an additional level of regulation over GPCR signaling that may contribute to our observed GRK2-mediated effects. GRK2 has also been shown to desensitize melanocortin-1 receptor (Sánchez-Más et al., 2005), which may contribute to sex-dependent analgesic effects of KOR agonists (Mogil et al., 2003). Together, these studies showed that GRK2 inhibition is sufficient to decrease KOR-mediated sex differences at a biochemical, neural circuit, and behavioral level in mice.
Whereas G protein-mediated analgesic effects were decreased, anti-pruritic effects of KOR activation remained evident in females. Although nalfurafine is a G-biased KOR compound (Schattauer et al., 2017), our results suggest that not all G protein-mediated effects may be subject to GRK2 regulation. For example, KOR activation can also disrupt cognition through G protein-mediated actions, but cognitive disruptions following KOR activation are not sexually dimorphic (Abraham et al., 2018). These findings suggest that there are some KOR/G protein-mediated behaviors that may not be regulated by GRK2 or may not require Gβγ activity. GRK2 interacts with Gβ1 and Gβ2, but not Gβ3 (Belcheva et al., 1998), indicating that there may be differential regulation of Gβ-mediated behaviors depending on the Gβ subunits or cell populations required for behavioral effects. Additionally, anti-pruritic and cognitive effects may be regulated by inhibitory interactions between Gαi and downstream effectors such as Ras, JNK, or protein kinase A (Edamatsu et al., 1998; Dessauer et al., 2002; Pierre et al., 2009). Munanairi et al. (2018) reported that KOR-mediated decreases in itch are sex independent and may not require G protein activation, instead using phospholipase C, supporting our observation of the anti-pruritic effects of KOR activation in males and females. The reported actions of KOR on diuresis (Leander, 1983; Craft et al., 2000) and locomotor activity (Robles et al., 2014) are also not sex dependent, suggesting that these behaviors are likely regulated by cellular and molecular mechanisms that are distinct from KOR-mediated analgesia. Additionally, KOR antagonists have also been reported to differ in efficacy between males and females (Laman-Maharg et al., 2018). To increase the efficacy and specificity of KOR agonists and antagonists, future studies could aim to dissociate the molecular signals or cellular populations producing analgesia from those producing diuretic, antipruritic, sedative, and cognitive effects of KOR activation.
Similar to the analgesic effects observed with KOR, blockade of GRK2 also promoted analgesic responses to morphine in females, a MOR-selective opiate (Goldstein and Naidu, 1989). GRK2 overexpression with MOR in cell culture increases activation of arrestin-dependent signals (Zhang et al., 1998), suggesting that GRK2 decreases MOR G protein-mediated signals, similar to GRK2-KOR interactions. Sex differences in morphine analgesia at low doses have been hypothesized to occur through regulation of G protein-coupled inwardly rectifying potassium channels (Mitrovic et al., 2003), which are also known to directly bind Gβγ for channel activation (Huang et al., 1995). Our studies demonstrated that sex differences in opioid receptor activity were the result of GRK2 activity and suggest that multiple inhibitory GPCRs could be subject to a similar sex-dependent regulation. Future studies could assess the role of GRK2 in decreasing Gαi protein-mediated behaviors in other GPCR systems.
Together, our studies demonstrate that KOR activation produces distinct analgesic effects in males and females in an estradiol- and GRK2-regulated manner. G protein-biased KOR analgesic compounds that are in clinical development may have sex-dependent effects and careful consideration of estradiol interactions with opioid receptor actions is required for the implementation of biased KOR agonists in clinical populations. Further understanding of the mechanisms underlying sexually dimorphic effects of opioid receptor activation could guide the development of highly efficacious nonaddictive opioid drugs for analgesia in males and females.
Footnotes
This work was supported by the National Institutes of Health (Grants P50 MH106428 and R01 DA030074 to C.C., Grant T32 DA07278 to A.D.A. and K.L.R., Grant T32 NS099578 to K.L.R., and Grant T32 GM007750 to J.H.C. and H.M.F.) and by the Brain & Behavior Research Foundation (NARSAD Young Investigator Award to B.B.L.). We thank Ashley A. Krull for assistance with ovariectomies.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Charles Chavkin, Department of Pharmacology, University of Washington, Box 357280, Seattle, WA 98195-7280. cchavkin{at}uw.edu