Molecular recognition of morphine and fentanyl by the human µ-opioid receptor.

Morphine and fentanyl are among the most used opioid drugs that confer analgesia and unwanted side effects through both G protein and arrestin signaling pathways of µ-opioid receptor (µOR). Here, we report structures of the human µOR-G protein complexes bound to morphine and fentanyl, which uncover key differences in how they bind the receptor. We also report structures of µOR bound to TRV130, PZM21, and SR17018, which reveal preferential interactions of these agonists with TM3 side of the ligand-binding pocket rather than TM6/7 side. In contrast, morphine and fentanyl form dual interactions with both TM3 and TM6/7 regions. Mutations at the TM6/7 interface abolish arrestin recruitment of µOR promoted by morphine and fentanyl. Ligands designed to reduce TM6/7 interactions display preferential G protein signaling. Our results provide crucial insights into fentanyl recognition and signaling of µOR, which may facilitate rational design of next-generation analgesics.

Dimerization and antidepressant recognition at noradrenaline transporter.

The noradrenaline transporter has a pivotal role in regulating neurotransmitter balance and is crucial for normal physiology and neurobiology1. Dysfunction of noradrenaline transporter has been implicated in numerous neuropsychiatric diseases, including depression and attention deficit hyperactivity disorder2. Here we report cryo-electron microscopy structures of noradrenaline transporter in apo and substrate-bound forms, and as complexes with six antidepressants. The structures reveal a noradrenaline transporter dimer interface that is mediated predominantly by cholesterol and lipid molecules. The substrate noradrenaline binds deep in the central binding pocket, and its amine group interacts with a conserved aspartate residue. Our structures also provide insight into antidepressant recognition and monoamine transporter selectivity. Together, these findings advance our understanding of noradrenaline transporter regulation and inhibition, and provide templates for designing improved antidepressants to treat neuropsychiatric disorders.

GPCR activation and GRK2 assembly by a biased intracellular agonist.

Phosphorylation of G-protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) desensitizes G-protein signalling and promotes arrestin signalling, which is also modulated by biased ligands1-6. The molecular assembly of GRKs on GPCRs and the basis of GRK-mediated biased signalling remain largely unknown owing to the weak GPCR-GRK interactions. Here we report the complex structure of neurotensin receptor 1 (NTSR1) bound to GRK2, Gαq and the arrestin-biased ligand SBI-5537. The density map reveals the arrangement of the intact GRK2 with the receptor, with the N-terminal helix of GRK2 docking into the open cytoplasmic pocket formed by the outward movement of the receptor transmembrane helix 6, analogous to the binding of the G protein to the receptor. SBI-553 binds at the interface between GRK2 and NTSR1 to enhance GRK2 binding. The binding mode of SBI-553 is compatible with arrestin binding but clashes with the binding of Gαq protein, thus providing a mechanism for its arrestin-biased signalling capability. In sum, our structure provides a rational model for understanding the details of GPCR-GRK interactions and GRK2-mediated biased signalling.

Molecular Basis for Hormone Recognition and Activation of Corticotropin-Releasing Factor Receptors.

Corticotropin-releasing factor (CRF) and the three related peptides urocortins 1-3 (UCN1-UCN3) are endocrine hormones that control the stress responses by activating CRF1R and CRF2R, two members of class B G-protein-coupled receptors (GPCRs). Here, we present two cryoelectron microscopy (cryo-EM) structures of UCN1-bound CRF1R and CRF2R with the stimulatory G protein. In both structures, UCN1 adopts a single straight helix with its N terminus dipped into the receptor transmembrane bundle. Although the peptide-binding residues in CRF1R and CRF2R are different from other members of class B GPCRs, the residues involved in receptor activation and G protein coupling are conserved. In addition, both structures reveal bound cholesterol molecules to the receptor transmembrane helices. Our structures define the basis of ligand-binding specificity in the CRF receptor-hormone system, establish a common mechanism of class B GPCR activation and G protein coupling, and provide a paradigm for studying membrane protein-lipid interactions for class B GPCRs.

Structures of Gi-bound metabotropic glutamate receptors mGlu2 and mGlu4.

The metabotropic glutamate receptors (mGlus) have key roles in modulating cell excitability and synaptic transmission in response to glutamate (the main excitatory neurotransmitter in the central nervous system)1. It has previously been suggested that only one receptor subunit within an mGlu homodimer is responsible for coupling to G protein during receptor activation2. However, the molecular mechanism that underlies the asymmetric signalling of mGlus remains unknown. Here we report two cryo-electron microscopy structures of human mGlu2 and mGlu4 bound to heterotrimeric Gi protein. The structures reveal a G-protein-binding site formed by three intracellular loops and helices III and IV that is distinct from the corresponding binding site in all of the other G-protein-coupled receptor (GPCR) structures. Furthermore, we observed an asymmetric dimer interface of the transmembrane domain of the receptor in the two mGlu-Gi structures. We confirmed that the asymmetric dimerization is crucial for receptor activation, which was supported by functional data; this dimerization may provide a molecular basis for the asymmetric signal transduction of mGlus. These findings offer insights into receptor signalling of class C GPCRs.

Structure of GPR101-Gs enables identification of ligands with rejuvenating potential.

GPR101 is an orphan G protein-coupled receptor actively participating in energy homeostasis. Here we report the cryo-electron microscopy structure of GPR101 constitutively coupled to Gs heterotrimer, which reveals unique features of GPR101, including the interaction of extracellular loop 2 within the 7TM bundle, a hydrophobic chain packing-mediated activation mechanism and the structural basis of disease-related mutants. Importantly, a side pocket is identified in GPR101 that facilitates in silico screening to identify four small-molecule agonists, including AA-14. The structure of AA-14-GPR101-Gs provides direct evidence of the AA-14 binding at the side pocket. Functionally, AA-14 partially restores the functions of GH/IGF-1 axis and exhibits several rejuvenating effects in wild-type mice, which are abrogated in Gpr101-deficient mice. In summary, we provide a structural basis for the constitutive activity of GPR101. The structure-facilitated identification of GPR101 agonists and functional analysis suggest that targeting this orphan receptor has rejuvenating potential.

Activation of the GLP-1 receptor by a non-peptidic agonist.

Class B G-protein-coupled receptors are major targets for the treatment of chronic diseases, including diabetes and obesity1. Structures of active receptors reveal peptide agonists engage deep within the receptor core, leading to an outward movement of extracellular loop 3 and the tops of transmembrane helices 6 and 7, an inward movement of transmembrane helix 1, reorganization of extracellular loop 2 and outward movement of the intracellular side of transmembrane helix 6, resulting in G-protein interaction and activation2-6. Here we solved the structure of a non-peptide agonist, TT-OAD2, bound to the glucagon-like peptide-1 (GLP-1) receptor. Our structure identified an unpredicted non-peptide agonist-binding pocket in which reorganization of extracellular loop 3 and transmembrane helices 6 and 7 manifests independently of direct ligand interaction within the deep transmembrane domain pocket. TT-OAD2 exhibits biased agonism, and kinetics of G-protein activation and signalling that are distinct from peptide agonists. Within the structure, TT-OAD2 protrudes beyond the receptor core to interact with the lipid or detergent, providing an explanation for the distinct activation kinetics that may contribute to the clinical efficacy of this compound series. This work alters our understanding of the events that drive the activation of class B receptors.

Structural analysis of the dual agonism at GLP-1R and GCGR.

Glucagon-like peptide-1 receptor (GLP-1R) and glucagon receptor (GCGR), two members of class B1 G protein-coupled receptors, play important roles in glucose homeostasis and energy metabolism. They share a high degree of sequence homology but have different functionalities. Unimolecular dual agonists of both receptors developed recently displayed better clinical efficacies than that of monotherapy. To study the underlying molecular mechanisms, we determined high-resolution cryo-electron microscopy structures of GLP-1R or GCGR in complex with heterotrimeric Gs protein and three GLP-1R/GCGR dual agonists including peptide 15, MEDI0382 (cotadutide) and SAR425899 with variable activating profiles at GLP-1R versus GCGR. Compared with related structures reported previously and supported by our published pharmacological data, key residues responsible for ligand recognition and dual agonism were identified. Analyses of peptide conformational features revealed a difference in side chain orientations within the first three residues, indicating that distinct engagements in the deep binding pocket are required to achieve receptor selectivity. The middle region recognizes extracellular loop 1 (ECL1), ECL2, and the top of transmembrane helix 1 (TM1) resulting in specific conformational changes of both ligand and receptor, especially the dual agonists reshaped ECL1 conformation of GLP-1R relative to that of GCGR, suggesting an important role of ECL1 interaction in executing dual agonism. Structural investigation of lipid modification showed a better interaction between lipid moiety of MEDI0382 and TM1-TM2 cleft, in line with its increased potency at GCGR than SAR425899. Together, the results provide insightful information for the design and development of improved therapeutics targeting these two receptors simultaneously.

Structural insights into the interaction of three Y-shaped ligands with PI3Kα.

Class IA phosphoinositide 3-kinase alpha (PI3Kα) is an important drug target because it is one of the most frequently mutated proteins in human cancers. However, small molecule inhibitors currently on the market or under development have safety concerns due to a lack of selectivity. Therefore, other chemical scaffolds or unique mechanisms of catalytic kinase inhibition are needed. Here, we report the cryo-electron microscopy structures of wild-type PI3Kα, the dimer of p110α and p85α, in complex with three Y-shaped ligands [cpd16 (compound 16), cpd17 (compound 17), and cpd18 (compound 18)] of different affinities and no inhibitory effect on the kinase activity. Unlike ATP-competitive inhibitors, cpd17 adopts a Y-shaped conformation with one arm inserted into a binding pocket formed by R770 and W780 and the other arm lodged in the ATP-binding pocket at an angle that is different from that of the ATP phosphate tail. Such a special interaction induces a conformation of PI3Kα resembling that of the unliganded protein. These observations were confirmed with two isomers (cpd16 and cpd18). Further analysis of these Y-shaped ligands revealed the structural basis of differential binding affinities caused by stereo- or regiochemical modifications. Our results may offer a different direction toward the design of therapeutic agents against PI3Kα.

Molecular basis of signal transduction mediated by the human GIPR splice variants.

Glucose-dependent insulinotropic polypeptide receptor (GIPR) is a potential drug target for metabolic disorders. It works with glucagon-like peptide-1 receptor and glucagon receptor in humans to maintain glucose homeostasis. Unlike the other two receptors, GIPR has at least 13 reported splice variants (SVs), more than half of which have sequence variations at either C or N terminus. To explore their roles in endogenous peptide-mediated GIPR signaling, we determined the cryoelectron microscopy (cryo-EM) structures of the two N terminus-altered SVs (referred as GIPR-202 and GIPR-209 in the Ensembl database, SV1 and SV2 here, respectively) and investigated the outcome of coexpressing each of them in question with GIPR in HEK293T cells with respect to ligand binding, receptor expression, cAMP (adenosine 3,5-cyclic monophosphate) accumulation, ß-arrestin recruitment, and cell surface localization. It was found that while both N terminus-altered SVs of GIPR neither bound to the hormone nor elicited signal transduction per se, they suppressed ligand binding and cAMP accumulation of GIPR. Meanwhile, SV1 reduced GIPR-mediated ß-arrestin 2 responses. The cryo-EM structures of SV1 and SV2 showed that they reorganized the extracellular halves of transmembrane helices 1, 6, and 7 and extracellular loops 2 and 3 to adopt a ligand-binding pocket-occupied conformation, thereby losing binding ability to the peptide. The results suggest a form of signal bias that is constitutive and ligand-independent, thus expanding our knowledge of biased signaling beyond pharmacological manipulation (i.e., ligand specific) as well as constitutive and ligand-independent (e.g., SV1 of the growth hormone-releasing hormone receptor).

Single hormone or synthetic agonist induces Gs/Gi coupling selectivity of EP receptors via distinct binding modes and propagating paths.

To accomplish concerted physiological reactions, nature has diversified functions of a single hormone at at least two primary levels: 1) Different receptors recognize the same hormone, and 2) different cellular effectors couple to the same hormone-receptor pair [R.P. Xiao, Sci STKE 2001, re15 (2001); L. Hein, J. D. Altman, B.K. Kobilka, Nature 402, 181-184 (1999); Y. Daaka, L. M. Luttrell, R. J. Lefkowitz, Nature 390, 88-91 (1997)]. Not only these questions lie in the heart of hormone actions and receptor signaling but also dissecting mechanisms underlying these questions could offer therapeutic routes for refractory diseases, such as kidney injury (KI) or X-linked nephrogenic diabetes insipidus (NDI). Here, we identified that Gs-biased signaling, but not Gi activation downstream of EP4, showed beneficial effects for both KI and NDI treatments. Notably, by solving Cryo-electron microscope (cryo-EM) structures of EP3-Gi, EP4-Gs, and EP4-Gi in complex with endogenous prostaglandin E2 (PGE2)or two synthetic agonists and comparing with PGE2-EP2-Gs structures, we found that unique primary sequences of prostaglandin E2 receptor (EP) receptors and distinct conformational states of the EP4 ligand pocket govern the Gs/Gi transducer coupling selectivity through different structural propagation paths, especially via TM6 and TM7, to generate selective cytoplasmic structural features. In particular, the orientation of the PGE2 ω-chain and two distinct pockets encompassing agonist L902688 of EP4 were differentiated by their Gs/Gi coupling ability. Further, we identified common and distinct features of cytoplasmic side of EP receptors for Gs/Gi coupling and provide a structural basis for selective and biased agonist design of EP4 with therapeutic potential.

Cryo-EM structures of cancer-specific helical and kinase domain mutations of PI3Kα.

Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that perform multiple and important cellular functions. The protein investigated here belongs to class IA of the PI3Ks; it is a dimer consisting of a catalytic subunit, p110α, and a regulatory subunit, p85α, and is referred to as PI3Kα. The catalytic subunit p110α is frequently mutated in cancer. The mutations induce a gain of function and constitute a driving force in cancer development. About 80% of these mutations lead to single-amino-acid substitutions in one of three sites of p110α: two in the helical domain of the protein (E542K and E545K) and one at the C-terminus of the kinase domain (H1047R). Here, we report the cryo-electron microscopy structures of these mutants in complex with the p110α-specific inhibitor BYL-719. The H1047R mutant rotates its sidechain to a new position and weakens the kα11 activation loop interaction, thereby reducing the inhibitory effect of p85α on p110α. E542K and E545K completely abolish the tight interaction between the helical domain of p110α and the N-terminal SH2 domain of p85α and lead to the disruption of all p85α binding and a dramatic increase in flexibility of the adaptor-binding domain (ABD) in p110α. Yet, the dimerization of PI3Kα is preserved through the ABD-p85α interaction. The local and global structural features induced by these mutations provide molecular insights into the activation of PI3Kα, deepen our understanding of the oncogenic mechanism of this important signaling molecule, and may facilitate the development of mutant-specific inhibitors.

Structural basis of peptidomimetic agonism revealed by small- molecule GLP-1R agonists Boc5 and WB4-24.

Glucagon-like peptide-1 receptor (GLP-1R) agonists are effective in treating type 2 diabetes and obesity with proven cardiovascular benefits. However, most of these agonists are peptides and require subcutaneous injection except for orally available semaglutide. Boc5 was identified as the first orthosteric nonpeptidic agonist of GLP-1R that mimics a broad spectrum of bioactivities of GLP-1 in vitro and in vivo. Here, we report the cryoelectron microscopy structures of Boc5 and its analog WB4-24 in complex with the human GLP-1R and Gs protein. Bound to the extracellular domain, extracellular loop 2, and transmembrane (TM) helices 1, 2, 3, and 7, one arm of both compounds was inserted deeply into the bottom of the orthosteric binding pocket that is usually accessible by peptidic agonists, thereby partially overlapping with the residues A8 to D15 in GLP-1. The other three arms, meanwhile, extended to the TM1-TM7, TM1-TM2, and TM2-TM3 clefts, showing an interaction feature substantially similar to the previously known small-molecule agonist LY3502970. Such a unique binding mode creates a distinct conformation that confers both peptidomimetic agonism and biased signaling induced by nonpeptidic modulators at GLP-1R. Further, the conformational difference between Boc5 and WB4-24, two closed related compounds, provides a structural framework for fine-tuning of pharmacological efficacy in the development of future small-molecule therapeutics targeting GLP-1R.

Nanobodies and chemical cross-links advance the structural and functional analysis of PI3Kα.

Nanobodies and chemical cross-linking were used to gain information on the identity and positions of flexible domains of PI3Kα. The application of chemical cross-linking mass spectrometry (CXMS) facilitated the identification of the p85 domains BH, cSH2, and SH3 as well as their docking positions on the PI3Kα catalytic core. Binding of individual nanobodies to PI3Kα induced activation or inhibition of enzyme activity and caused conformational changes that could be correlated with enzyme function. Binding of nanobody Nb3-126 to the BH domain of p85α substantially improved resolution for parts of the PI3Kα complex, and binding of nanobody Nb3-159 induced a conformation of PI3Kα that is distinct from known PI3Kα structures. The analysis of CXMS data also provided mechanistic insights into the molecular underpinning of the flexibility of PI3Kα.

Abnormal global alternative RNA splicing in COVID-19 patients.

Viral infections can alter host transcriptomes by manipulating host splicing machinery. Despite intensive transcriptomic studies on SARS-CoV-2, a systematic analysis of alternative splicing (AS) in severe COVID-19 patients remains largely elusive. Here we integrated proteomic and transcriptomic sequencing data to study AS changes in COVID-19 patients. We discovered that RNA splicing is among the major down-regulated proteomic signatures in COVID-19 patients. The transcriptome analysis showed that SARS-CoV-2 infection induces widespread dysregulation of transcript usage and expression, affecting blood coagulation, neutrophil activation, and cytokine production. Notably, CD74 and LRRFIP1 had increased skipping of an exon in COVID-19 patients that disrupts a functional domain, which correlated with reduced antiviral immunity. Furthermore, the dysregulation of transcripts was strongly correlated with clinical severity of COVID-19, and splice-variants may contribute to unexpected therapeutic activity. In summary, our data highlight that a better understanding of the AS landscape may aid in COVID-19 diagnosis and therapy.

Highly Luminescent Chiral Double π-Helical Nanoribbons.

Unveiling the mechanism behind chirality propagation and dissymmetry amplification at the molecular level is of significance for the development of chiral systems with comprehensively outstanding chiroptical performances. Herein, we have presented a straightforward Cu-mediated Ullmann homocoupling approach to synthesize perylene diimide-entwined double π-helical nanoribbons encompassing dimer, trimer, and tetramer while producing homochiral or heterochiral linking of chiral centers. A significant dissymmetry amplification was achieved, with absorption dissymmetry factors (|gabs|) increasing from 0.009 to 0.017 and further to 0.019, and luminescence dissymmetry factors (|glum|) rising from 0.007 to 0.013 and eventually to 0.015 for homochiral double π-helical oligomers. The disparity of magnetic transition dipole moment (m) densities in homochiral and heterochiral tetramers by time-dependent density functional theory calculations confirmed that homochiral oligomerization can maximize the total m, which is favorable for achieving ever-increasing g factors. Notably, these double π-helices exhibited exceptional photoluminescence quantum yields (ΦPL) ranging from 83 to 95%. The circularly polarized luminescence brightness (BCPL) eventually reached a remarkable 575 M-1 cm-1 for the homochiral tetramer, which is among the highest values reported for chiral small molecules. This kind of linearly extended double π-helices offers a platform for a comprehensive understanding of the mechanism behind chirality propagation and dissymmetry amplification.

Phase-plate cryo-EM structure of a biased agonist-bound human GLP-1 receptor-Gs complex.

The class B glucagon-like peptide-1 (GLP-1) G protein-coupled receptor is a major target for the treatment of type 2 diabetes and obesity. Endogenous and mimetic GLP-1 peptides exhibit biased agonism-a difference in functional selectivity-that may provide improved therapeutic outcomes. Here we describe the structure of the human GLP-1 receptor in complex with the G protein-biased peptide exendin-P5 and a Gαs heterotrimer, determined at a global resolution of 3.3 Å. At the extracellular surface, the organization of extracellular loop 3 and proximal transmembrane segments differs between our exendin-P5-bound structure and previous GLP-1-bound GLP-1 receptor structure. At the intracellular face, there was a six-degree difference in the angle of the Gαs-α5 helix engagement between structures, which was propagated across the G protein heterotrimer. In addition, the structures differed in the rate and extent of conformational reorganization of the Gαs protein. Our structure provides insights into the molecular basis of biased agonism.

Structure of the glucagon receptor in complex with a glucagon analogue.

Class B G-protein-coupled receptors (GPCRs), which consist of an extracellular domain (ECD) and a transmembrane domain (TMD), respond to secretin peptides to play a key part in hormonal homeostasis, and are important therapeutic targets for a variety of diseases. Previous work has suggested that peptide ligands bind to class B GPCRs according to a two-domain binding model, in which the C-terminal region of the peptide targets the ECD and the N-terminal region of the peptide binds to the TMD binding pocket. Recently, three structures of class B GPCRs in complex with peptide ligands have been solved. These structures provide essential insights into peptide ligand recognition by class B GPCRs. However, owing to resolution limitations, the specific molecular interactions for peptide binding to class B GPCRs remain ambiguous. Moreover, these previously solved structures have different ECD conformations relative to the TMD, which introduces questions regarding inter-domain conformational flexibility and the changes required for receptor activation. Here we report the 3.0 Å-resolution crystal structure of the full-length human glucagon receptor (GCGR) in complex with a glucagon analogue and partial agonist, NNC1702. This structure provides molecular details of the interactions between GCGR and the peptide ligand. It reveals a marked change in the relative orientation between the ECD and TMD of GCGR compared to the previously solved structure of the inactive GCGR-NNC0640-mAb1 complex. Notably, the stalk region and the first extracellular loop undergo major conformational changes in secondary structure during peptide binding, forming key interactions with the peptide. We further propose a dual-binding-site trigger model for GCGR activation-which requires conformational changes of the stalk, first extracellular loop and TMD-that extends our understanding of the previously established two-domain peptide-binding model of class B GPCRs.

CCDC92 promotes podocyte injury by regulating PA28α/ABCA1/cholesterol efflux axis in type 2 diabetic mice.

Podocyte lipotoxicity mediated by impaired cellular cholesterol efflux plays a crucial role in the development of diabetic kidney disease (DKD), and the identification of potential therapeutic targets that regulate podocyte cholesterol homeostasis has clinical significance. Coiled-coil domain containing 92 (CCDC92) is a novel molecule related to metabolic disorders and insulin resistance. However, whether the expression level of CCDC92 is changed in kidney parenchymal cells and the role of CCDC92 in podocytes remain unclear. In this study, we found that Ccdc92 was significantly induced in glomeruli from type 2 diabetic mice, especially in podocytes. Importantly, upregulation of Ccdc92 in glomeruli was positively correlated with an increased urine albumin-to-creatinine ratio (UACR) and podocyte loss. Functionally, podocyte-specific deletion of Ccdc92 attenuated proteinuria, glomerular expansion and podocyte injury in mice with DKD. We further demonstrated that Ccdc92 contributed to lipid accumulation by inhibiting cholesterol efflux, finally promoting podocyte injury. Mechanistically, Ccdc92 promoted the degradation of ABCA1 by regulating PA28α-mediated proteasome activity and then reduced cholesterol efflux. Thus, our studies indicate that Ccdc92 contributes to podocyte injury by regulating the PA28α/ABCA1/cholesterol efflux axis in DKD.

Activation of PTH1R alleviates epididymitis and orchitis through Gq and ß-arrestin-1 pathways.

Inflammation in the epididymis and testis contributes significantly to male infertility. Alternative therapeutic avenues treating epididymitis and orchitis are expected since current therapies using antibiotics have limitations associated to side effects and are commonly ineffective for inflammation due to nonbacterial causes. Here, we demonstrated that type 1 parathyroid hormone receptor (PTH1R) and its endogenous agonists, parathyroid hormone (PTH) and PTH-related protein (PTHrP), were mainly expressed in the Leydig cells of testis as well as epididymal epithelial cells. Screening the secretin family G protein-coupled receptor identified that PTH1R in the epididymis and testis was down-regulated in mumps virus (MuV)- or lipopolysaccharide (LPS)-induced inflammation. Remarkably, activation of PTH1R by abaloparatide (ABL), a Food and Drug Administration-approved treatment for postmenopausal osteoporosis, alleviated MuV- or LPS-induced inflammatory responses in both testis and epididymis and significantly improved sperm functions in both mouse model and human samples. The anti-inflammatory effects of ABL were shown to be regulated mainly through the Gq and ß-arrestin-1 pathway downstream of PTH1R as supported by the application of ABL in Gnaq± and Arrb1-/- mouse models. Taken together, our results identified an important immunoregulatory role for PTH1R signaling in the epididymis and testis. Targeting to PTH1R might have a therapeutic effect for the treatment of epididymitis and orchitis or other inflammatory disease in the male reproductive system.