Towards A Molecular Understanding of The Cannabinoid Related Orphan Receptor GPR18: A Focus on Its Constitutive Activity

New insights into sperm physiology regulation: Enlightenment from G-protein-coupled receptors

BACKGROUND

G-protein-coupled receptors are critical in many physiological and pathological processes in various organs. Serving as the control panel for sensing extracellular stimuli, G-protein-coupled receptors recognise various ligands, including light, temperature, odours, pheromones, hormones, neurotransmitters, chemokines, etc. Most recently, G-protein-coupled receptors residing in spermatozoa have been found to be indispensable for sperm function.

OBJECTIVE

Here, we have summarised cutting-edge findings on the functional mechanisms of G-protein-coupled receptors that are known to be associated with sperm functions and the activation of their downstream effectors, providing new insights into the roles of G-protein-coupled receptors in sperm physiology.

RESULTS

Emerging studies hint that alterations in G-protein-coupled receptors could affect sperm function, implicating their role in fertility, but solid evidence needs to be continuing excavated with various means. Several members of the G-protein-coupled receptor superfamily, including olfactory receptors, opsins, orphan G-protein-coupled receptors, CXC chemokine receptor 4, CC chemokine receptor 5 and CC chemokine receptor 6 as well as their downstream effector β-arrestins, etc., were suggested to be essential for sperm motility, capacitation, thermotaxis, chemotaxis, Ca2+ influx through CatSper channel and fertilisation capacity.

CONCLUSION

The present review provides a comprehensive overview of studies describing G-protein-coupled receptors and their potential action in sperm function. We also present a critical discussion of these issues, and a possible framework for future investigations on the diverse ligands, biological functions and cell signalling of G-protein-coupled receptors in spermatozoa. Here, the G-protein-coupled receptors and their related G proteins that specifically were identified in spermatozoa were summarised, and provided references valuable for further illumination, despite the evidence that is not overwhelming in most cases.

N-arachidonylglycine is a caloric state-dependent circulating metabolite which regulates human CD4+T cell responsiveness

Caloric deprivation interventions such as intermittent fasting and caloric restriction ameliorate metabolic and inflammatory disease. As a human model of caloric deprivation, a 24-h fast blunts innate and adaptive immune cell responsiveness relative to the refed state. Isolated serum at these time points confers these same immunomodulatory effects on transformed cell lines. To identify serum mediators orchestrating this, metabolomic and lipidomic analysis was performed on serum extracted after a 24-h fast and re-feeding. Bioinformatic integration with concurrent peripheral blood mononuclear cells RNA-seq analysis implicated key metabolite-sensing GPCRs in fasting-mediated immunomodulation. The putative GPR18 ligand N-arachidonylglycine (NAGly) was elevated during fasting and attenuated CD4⁺T cell responsiveness via GPR18 MTORC1 signaling. In parallel, NAGly reduced inflammatory Th1 and Th17 cytokines levels in CD4⁺T cells isolated from obese subjects, identifying a fasting-responsive metabolic intermediate that may contribute to the regulation of nutrient-level dependent inflammation associated with metabolic disease.

Structure Prediction, Evaluation, and Validation of GPR18 Lipid Receptor Using Free Programs

The GPR18 receptor, often referred to as the N-arachidonylglycine receptor, although assigned (along with GPR55 and GPR119) to the new class A GPCR subfamily-lipid receptors, officially still has the status of a class A GPCR orphan. While its signaling pathways and biological significance have not yet been fully elucidated, increasing evidence points to the therapeutic potential of GPR18 in relation to immune, neurodegenerative, and cancer processes to name a few. Therefore, it is necessary to understand the interactions of potential ligands with the receptor and the influence of particular structural elements on their activity. Thus, given the lack of an experimentally solved structure, the goal of the present study was to obtain a homology model of the GPR18 receptor in the inactive state, meeting all requirements in terms of protein structure quality and recognition of active ligands. To increase the reliability and precision of the predictions, different contemporary protein structure prediction methods and software were used and compared herein. To test the usability of the resulting models, we optimized and compared the selected structures followed by the assessment of the ability to recognize known, active ligands. The stability of the predicted poses was then evaluated by means of molecular dynamics simulations. On the other hand, most of the best-ranking contemporary CADD software/platforms for its full usability require rather expensive licenses. To overcome this down-to-earth obstacle, the overarching goal of these studies was to test whether it is possible to perform the thorough CADD experiments with high scientific confidence while using only license-free/academic software and online platforms. The obtained results indicate that a wide range of freely available software and/or academic licenses allow us to carry out meaningful molecular modelling/docking studies.

Targeting CB2 and TRPV1: Computational Approaches for the Identification of Dual Modulators

Both metabotropic (CBRs) and ionotropic cannabinoid receptors (ICRs) have implications in a range of neurological disorders. The metabotropic canonical CBRs CB1 and CB2 are highly implicated in these pathological events. However, selective targeting at CB2 versus CB1 offers optimized pharmacology due to the absence of psychoactive outcomes. The ICR transient receptor potential vanilloid type 1 (TRPV1) has also been reported to play a role in CNS disorders. Thus, activation of both targets, CB2 and TRPV1, offers a promising polypharmacological strategy for the treatment of neurological events including analgesia and neuroprotection. This brief research report aims to identify chemotypes with a potential dual CB2/TRPV1 profile. For this purpose, we have rationalized key structural features for activation and performed virtual screening at both targets using curated chemical libraries.

Therapeutic Exploitation of GPR18: Beyond the Cannabinoids?

GPR18 is a G-protein-coupled receptor that belongs to the orphan class A family. Even though it shares low sequence homology with the cannabinoid receptors CB₁R and CB₂R, a growing body of research suggests its relationship with the endocannabinoid system, not only because it is able to recognize cannabinoid ligands but also because of its expression and ability to heteromerize with CBRs. In this review, we aim to analyze the biological relevance, reported modulators, and structural features of GPR18. In order to guide future drug design in this field, highlights from molecular modeling of GPR18 will be provided.

Computational Investigations on the Binding Mode of Ligands for the Cannabinoid-Activated G Protein-Coupled Receptor GPR18

GPR18 is an orphan G protein-coupled receptor (GPCR) expressed in cells of the immune system. It is activated by the cannabinoid receptor (CB) agonist ∆⁹-tetrahydrocannabinol (THC). Several further lipids have been proposed to act as GPR18 agonists, but these results still require unambiguous confirmation. In the present study, we constructed a homology model of the human GPR18 based on an ensemble of three GPCR crystal structures to investigate the binding modes of the agonist THC and the recently reported antagonists which feature an imidazothiazinone core to which a (substituted) phenyl ring is connected via a lipophilic linker. Docking and molecular dynamics simulation studies were performed. As a result, a hydrophobic binding pocket is predicted to accommodate the imidazothiazinone core, while the terminal phenyl ring projects towards an aromatic pocket. Hydrophobic interaction of Cys251 with substituents on the phenyl ring could explain the high potency of the most potent derivatives. Molecular dynamics simulation studies suggest that the binding of imidazothiazinone antagonists stabilizes transmembrane regions TM1, TM6 and TM7 of the receptor through a salt bridge between Asp118 and Lys133. The agonist THC is presumed to bind differently to GPR18 than to the distantly related CB receptors. This study provides insights into the binding mode of GPR18 agonists and antagonists which will facilitate future drug design for this promising potential drug target.

Exploring immune-related genes with prognostic value in microenvironment of breast cancer from TCGA database

Breast cancer is one of the most common malignancies in women worldwide. Many studies have shown that tumor microenvironment cells, immune cells, and stromal cell infiltration have an important impact on prognosis, so it is important to identify biomarkers for achieving better treatment and prognosis.To better understand the relationship between immune and stromal cell-related genes and prognosis, we screened patients with breast cancer in The Cancer Genome Atlas (TCGA) database and divided them into high and low groups based on immune/stromal scores. We next identified differentially expressed immune-related genes that are significantly associated with the prognosis of patients with breast cancer for functional enrichment analysis and protein-protein interaction networks, respectively. Finally, we selected a separate breast cancer cohort in gene expression synthesis (GEO) for validation.Both immune scores and stromal scores are meaningful in the correlation of subtype classification. Disease-free survival of cases with the high score group of immune scores is statistically longer than the cases in the low score group. Differentially expressed immune-related genes extracted from the comparison can effectively evaluate the prognosis of patients with breast cancer and these genes are primarily involved in immune responses, extracellular matrix, and chemokine activity. At last, we obtained a series of verified tumor immune-related genes that predict the prognosis of patients with breast cancer.Combining the Estimation of Stromal and Immune Cells in Malignant Tumor Tissues using Expression database and the TCGA database to extract the list of tumor microenvironment related genes which may help to outline the prognosis of patients with breast cancer. Some previously overlooked genes have the potential to become additional biomarkers for breast cancer. Further research on these genes can reveal a new understanding of the potential relationship between tumor microenvironment and breast cancer prognosis.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.