ORIGINAL RESEARCH article

Front. Neuroinform., 07 February 2024

Volume 18 - 2024 | https://doi.org/10.3389/fninf.2024.1294939

In silico analyses of the involvement of GPR55, CB1R and TRPV1: response to THC, contribution to temporal lobe epilepsy, structural modeling and updated evolution

  • 1. Worcester Biomedical Research Group, School of Science and the Environment, University of Worcester, Worcester, United Kingdom

  • 2. Sustainable Environments Research Group, School of Science and the Environment University of Worcester, Worcester, United Kingdom

  • 3. School of Science and the Environment University of Worcester, Worcester, United Kingdom

Article metrics

View details

6

Citations

7,9k

Views

2,3k

Downloads

Abstract

Introduction:

The endocannabinoid (eCB) system is named after the discovery that endogenous cannabinoids bind to the same receptors as the phytochemical compounds found in Cannabis. While endogenous cannabinoids include anandamide (AEA) and 2-arachidonoylglycerol (2-AG), exogenous phytocannabinoids include Δ-9 tetrahydrocannabinol (THC) and cannabidiol (CBD). These compounds finely tune neurotransmission following synapse activation, via retrograde signaling that activates cannabinoid receptor 1 (CB1R) and/or transient receptor potential cation channel subfamily V member 1 (TRPV1). Recently, the eCB system has been linked to several neurological diseases, such as neuro-ocular abnormalities, pain insensitivity, migraine, epilepsy, addiction and neurodevelopmental disorders. In the current study, we aim to: (i) highlight a potential link between the eCB system and neurological disorders, (ii) assess if THC exposure alters the expression of eCB-related genes, and (iii) identify evolutionary-conserved residues in CB1R or TRPV1 in light of their function.

Methods:

To address this, we used several bioinformatic approaches, such as transcriptomic (Gene Expression Omnibus), protein–protein (STRING), phylogenic (BLASTP, MEGA) and structural (Phyre2, AutoDock, Vina, PyMol) analyzes.

Results:

Using RNA sequencing datasets, we did not observe any dysregulation of eCB-related transcripts in major depressive disorders, bipolar disorder or schizophrenia in the anterior cingulate cortex, nucleus accumbens or dorsolateral striatum. Following in vivo THC exposure in adolescent mice, GPR55 was significantly upregulated in neurons from the ventral tegmental area, while other transcripts involved in the eCB system were not affected by THC exposure. Our results also suggest that THC likely induces neuroinflammation following in vitro application on mice microglia. Significant downregulation of TPRV1 occurred in the hippocampi of mice in which a model of temporal lobe epilepsy was induced, confirming previous observations. In addition, several transcriptomic dysregulations were observed in neurons of both epileptic mice and humans, which included transcripts involved in neuronal death. When scanning known interactions for transcripts involved in the eCB system (n = 12), we observed branching between the eCB system and neurophysiology, including proteins involved in the dopaminergic system. Our protein phylogenic analyzes revealed that CB1R forms a clade with CB2R, which is distinct from related paralogues such as sphingosine-1-phosphate, receptors, lysophosphatidic acid receptors and melanocortin receptors. As expected, several conserved residues were identified, which are crucial for CB1R receptor function. The anandamide-binding pocket seems to have appeared later in evolution. Similar results were observed for TRPV1, with conserved residues involved in receptor activation.

Conclusion:

The current study found that GPR55 is upregulated in neurons following THC exposure, while TRPV1 is downregulated in temporal lobe epilepsy. Caution is advised when interpreting the present results, as we have employed secondary analyzes. Common ancestors for CB1R and TRPV1 diverged from jawless vertebrates during the late Ordovician, 450 million years ago. Conserved residues are identified, which mediate crucial receptor functions.

1 Introduction

The endocannabinoid (eCB) system acts as a regulator during neurotransmission. Indeed, following synaptic activation, eCBs will be released in the post-synapse, which, in turn, will decrease neuronal activation (Lu and Mackie, 2016). Historically, anandamide (AEA) and 2-arachidonoylglycerol (2-AG) were discovered first and named after the receptors they involve, mimicking the effects observed after “cannabis” exposure (Lu and Mackie, 2016). Expression of the cannabinoid receptor type 1 (CB1R) is mostly confined to the central nervous system (Matsuda et al., 1990), while CB2R expression is considered to be mainly located at the periphery (Munro et al., 1993), although potential clinical applications of targeting CB2R in neurological disorders have been highlighted recently (Kibret et al., 2022), likely due to microglial expression of CB2R (Tanaka et al., 2020). Other known eCB receptors include transient receptor potential vanilloid 1 (TRPV1) and G protein-coupled receptor 55 (GPR55) (Kano et al., 2009). Furthermore, GPR18 (McHugh et al., 2010), GPR119 (Overton et al., 2006) and some PPARs [peroxisome proliferator-activated receptors, reviewed in 2016 (O’Sullivan, 2016)] have also been identified as potential eCB receptors.

Derived from the omega-6 arachidonic acid (C20:4 n-6), both AEA and 2-AG are endogenous ligands of CB1R and induce decreased neuronal excitation (Ameri et al., 1999; Wilson and Nicoll, 2001; Straiker et al., 2023). Both molecules seem to involve a retrograde mechanism, whereby synthesis of AEA or 2-AG is located in the post-synapse, followed by diffusion in the synaptic cleft to then act on pre-synaptic vesicular fusion (Wilson and Nicoll, 2001). However, non-retrograde signaling has also been observed, as well as neuron-astrocyte communication (Castillo et al., 2012). The biosynthesis of AEA involves activation of N-acyltransferase (NAT) and N-acylphosphatidylethanolamine-hydrolyzing phospholipase D (NAPE-PLD), while 2-AG requires both phospholipase C (PLC) and diacylglycerol lipase (DAGL) (Kano et al., 2009). Degradation of AEA is performed by activation of fatty acid amide hydrolase (FAAH) and 2-AG is broken down by monoacyl-glycerol lipase (MAGL) (Kano et al., 2009). Interestingly, other degradation pathways can be triggered through two other enzymes: cyclooxygenase and lipoxygenase (Kano et al., 2009). The biological activity and regional distributions of these enzymes, at both the organ- and cell-specific levels, have been summarized previously (Basavarajappa, 2007).

Recent evidence has pinpointed that dysfunctions of the eCB system may cause neuro-ocular abnormalities (Bainbridge et al., 2022), pain insensitivity (Habib et al., 2019), migraine (Juhasz et al., 2009), epilepsy (Ludányi et al., 2008), addiction (Herman et al., 2006; Agrawal et al., 2009; Proudnikov et al., 2010), metabolic disease (de Miguel-Yanes et al., 2011), obesity (Russo et al., 2007) and neurodevelopmental disorders (Lu et al., 2008; Jung et al., 2012; Siniscalco et al., 2013; Gomis-González et al., 2016; Servadio et al., 2016; Wei et al., 2016; Barchel et al., 2018; Karhson et al., 2018; Melancia et al., 2018; Aran et al., 2019; Zamberletti et al., 2019; Wu et al., 2020). Many clinical observations have been paralleled and explained by animal studies, often linking AEA and neurodevelopmental disorders. In a recent meta-analysis, medical cannabis, largely used in the form of either tetrahydrocannabinol (THC) or cannabidiol (CBD), was found to be efficient for the treatment of several neurological diseases, such as chronic pain, spasticity, substance use disorder, epilepsy, Tourette’s syndrome and Parkinson’s Disease (Bilbao and Spanagel, 2022). Thus, activation of the eCB by natural or synthetic agonists provides efficient clinical benefits, which could be explained by the negative feedback role of the eCB system in neurotransmission.

The most probable evolutionary origin of the eCB system was explained in 2006 (McPartland et al., 2006). In fact, different evolutionary origins have been discovered, depending on which eCB is scrutinized. TRPV1 and GPR55 arose rather recently, likely in mammals (McPartland et al., 2006). The oldest member of the eCB family, FAAH, was present in early eukaryotic cells, while NAPE-PLD arose in opisthokonts, which comprise both animals and fungi. Recent evidence has corroborated previous observations, as two FAAH isoforms have been studied in the legume Medicago truncatula (Arias-Gaguancela et al., 2023), in addition to the known expression in the angiosperm Arabidopsis thaliana (Haq and Kilaru, 2020). Furthermore, a duplication event likely induced the appearance of CB1R and CB2R in vertebrates from an invertebrate deuterostome ancestor (Elphick et al., 2003), in which the pre-CBR protein was responsible for axonal regulation (Egertová and Elphick, 2007).

In this article, we aim to:

  • (i) Assess the potential transcriptomic alterations of 13 specific eCB-related transcripts in epilepsy, major depressive disorder, bipolar disorder or schizophrenia.

  • (ii) Assess the impact of THC on these 13 specific eCB-related transcripts.

  • (iii) Identify paralogues of CB1R or TRPV1 to infer evolutionary ancestry.

  • (iv) Identify conserved residues in CB1R or TRPV1 that mediate receptor function.

2 Methods

2.1 In silico analysis of publicly available RNA sequencing datasets

Previously-published RNA sequencing datasets were retrieved from the Gene Expression Omnibus (GEO) repository (Edgar et al., 2002; Clough and Barrett, 2016). Briefly, GEO datasets were scanned for adequate output using specific search terms. Table 1 describes the datasets used in the current study, with corresponding accession links given in Supplementary Table S1. We focused on 13 specific genes from the eCB system (Supplementary Table S2). Supplementary material M1 summarizes all DESeq2 results of the present study.

Table 1

DatasetSpeciePublicationData typeNorm needed?SettingTissueRegionsIntervention groupControl groupBehavioral testing performed prior to RNA-seq
GSE80655HsRamaker et al. (2017)Raw aligned countsYesPost mortemBrainAntC, NAC, DLSMDD, BD or SCZHealthyn/a
GSE189821MmZuo et al. (2022)Raw aligned countsYesIn vivoBrainPFC, NAC, AntC, DLS, VTAIntraperitoneal THC,
10 mg/kg/day for 21 days, P30-44 start, washout 14 days
Vehicle ethanol:tween:saline,
1:1:18
Yes
GSE116813MmJouroukhin et al. (2019)Raw aligned countsyesIn vivobrainHCSubcutaneous THC,
8 mg/kg/day for 21 days, P30 start, washout 21 days
Vehicle ethanol:cremaphor:saline,
1:1:18
Yes
GSE70689MmJuknat et al. (2013)Raw aligned countsNoIn vitroMicroglia-THC, 10 μM, 2 h.None, as final
[ethanol] 0.1%
n/a
GSE77578MmSrivastava et al. (2018)Raw aligned countsnoIn vivoBrainHCEP induced by single intraperitoneal pilocarpine injectionNon-epileptic (pilocarpine-naïve)No
GSE190451HsChen et al. (2023)Raw aligned countsYesPost mortemBrainNeoCEP (temporal lobe)Non-epileptic patients (undergoing benign brain tumor surgery, healthy tissue collected)n/a
GSE74150RnDamasceno et al. (2020)Raw aligned countsYesIn vivoBrainInferior colliculusEP (audiogenic)Non-epilepticNo

Detailed description of all RNA sequencing datasets retrieved for the present study.

GSE, GEO (gene expression omnibus) series experient; Hs, Homo sapiens; Mm, Mus musculus; Rn, Rattus norvegicus; NAC, nucleus accumbens, AntC, anterior cingulate; DLS, dorsolateral striatum; HC, hippocampus; NeoC, neocortex; PFC, prefrontal cortex; VTA, ventral tegmental area; MDD, major depressive disorder; BD, bipolar disorder; SCZ, schizophrenia; EP, epilepsy; THC, tetrahydrocannabinol; Norm, normalization; n/a, non-applicable. Please see Supplementary Table S1 for accession links.

Differential gene expression analysis was performed in the freeware R Studio version 4.0.3 (R Core Team, 2013; RStudio Team, 2020). Following the workflow described by Sanchis et al. (2021), we first retrieved count data, using the GEOquery package. If necessary, count data were then normalized using the DESeq2 package. Differential expression analyzes were performed using the DESeq2 package. Gene annotations were performed with either the org.Hs.eg.bd (Homo sapiens), org.Mm.eg.db (Mus musculus) or AnnotationDbi packages, as appropriate. Volcano plots were drawn with the EnhancedVolcano package. Gene ontology (over-representation analysis, ORA) was performed with the GOplot package using Biological Processes. Gene set enrichment analysis (GSEA) was performed with clusterProfiler and enrichplot packages. Scatterplot graphs were drawn using the ggplot2 package. All of the above-mentioned packages are available under open access at the Comprehensive R Archive Network (CRAN) repository.1

2.2 Protein phylogenetics

Phylogenetic analyzes were performed in Mega X (Kumar et al., 2018) following alignment of FASTA sequences. Evolutional analysis was performed using the maximum likelihood method. Initial trees were obtained by applying Neighbor-Join (NJ) and BioNJ algorithms (Gascuel, 1997) to a matrix of pairwise distances estimated using the Jones-Taylor-Thornton (JTT)-based model (Jones et al., 1992). The final drawn trees are those with the highest log likelihoods. One hundred bootstraps (r = 100) were used to infer confidence in the output trees. BLASTP was used to scan different species with Homo sapiens as bait (reference sequence), available at the National Center for Biotechnology Information (NCBI).2

To search for conserved domains/residues, protein sequence (Homo sapiens) alignments was performed in Clustal Omega, available at the European Molecular Biology Laboratory – European Bioinformatics Institute (EMBL-EBI).3

2.3 Protein structural analyzes

Structural prediction of transmembrane protein topology (TMHMM, Transmembrane prediction using Hidden Markov Models) was performed from protein FASTA sequences using the platform provided by the Department of Health Technology at the Technical University of Denmark,4 which was made publicly available in 2001 (Krogh et al., 2001). Homology models were created using Phyre2 (Kelley and Sternberg, 2009; Kelley et al., 2015). For docking of anandamide into TRPV1, the Cryo-EM structure of human TRPV1 in complex with the analgesic drug SB-366791 (PDB: 8GFA) was used. The A chain of the multimer was used with ligand and phospholipids removed from the .pdb file. The anandamide ligand (CHEMBL15848) was retrieved from the ChEMBL webserver as a .csv file and converted to a .pdb file using the Online SMILES Translator and Structure File Generator and then to a .pdbqt file with merged non-polar hydrogens with AutoDockTools version 1.5.7. The ligand was docked into the receptor using AutoDock Vina 1.1.2 (Trott and Olson, 2010) with a docking grid 20 Å in length in the x, y, and z directions centered on the “tunnel” anandamide-binding pocket previously identified (Morales et al., 2022), flexible N438, Y555, Y554, Y487 and D708 side-chains and an exhaustiveness of 20. Structural analysis was performed by superposition in PyMol (Schrödinger and DeLano, 2020). A professional academic license was acquired for PyMol.

2.4 Protein–protein interaction networks

Protein–protein interaction networks (Szklarczyk et al., 2019) were predicted using the STRING (Search Tool for the Retrieval of Interacting Genes) database server,5 which contains both known and predicted protein interactions. Queries were pasted with either a single protein name or multiple names (as appropriate). When relevant, outer networks were searched for by expanding the nodes on display, using the built-in database tool, to reflect the immediate, intermediate and extended protein interactomes. For protein clustering within the networks, k mean clustering (based on centroids) was applied (Dalmaijer et al., 2022) directly in the server.

3 Results

3.1 Neurodevelopmental disorders and the eCB system: no change for 13 selected transcripts

We first examined dataset GSE80655 (Ramaker et al., 2017), which contains RNA sequencing of 281 postmortem samples extracted from patients with either major depressive disorder (MDD, n = 69), bipolar disorder (BP, n = 71), schizophrenia (SCZ, n = 71) or controls (n = 70), in an attempt to decipher the potential role of the eCB system in these pathologies (Table 1). Out of 13 transcripts of interest (CNR1, CNR2, GPR55, GPR110, TRPV1, MGLL, DAGLA, PLCB1, NAPEPLD, FAAH, ABHD6, PLA2G4A and ALOX8, Supplementary Table S2), none were considered significantly up- or down-regulated in MDD (Figure 1A), respectively defined as log2(fold change) >1 or < −1 and -log10(p) ≥ 1.3 (p = 0.05). Similar results were observed with BP (Figure 1B) and SCZ (Figure 1C). These results indicate that our 13 transcripts of interest are not dysregulated under these conditions.

Figure 1

3.2 Chronic THC induces significant up-regulation of GPR55 in the ventral tegmental area of mice

Next, we investigated whether we could detect transcriptional changes in our 13 transcripts of interest following acute or chronic THC exposure. To that end, we used GSE189821, GSE116813 and GSE70689 (Table 1). In GSE189821, 119 samples were analyzed from adolescent mice chronically exposed to THC (3 weeks at 10 mg/kg/day) or not (vehicle) and following a 14-day washout period, mimicking human behavior. In accordance with previous observations (Zuo et al., 2022), we found only minor transcriptomic dysregulations of our 13 eCB-related transcripts of interest following THC exposure in these mice (males and females pooled). This might be due to the long washout period that was applied. However, we noted 31 differentially expressed genes, with 13 transcripts down-regulated and 18 transcripts up-regulated, outside of our list of interest, when pooling all different brain regions together (Figure 2A). In the ventral tegmental area, GPR55 was found to be significantly upregulated after THC exposure, the only one from our list of interest and out of all brain regions examined. In GSE116813, 12 hippocampal samples were sequenced from mice exposed to chronic THC (21 days at 8 mg/kg/day) or not (Jouroukhin et al., 2019) and after 21 days of washout. After subsetting the dataset to only include wild-type mice, we observed 134 differentially expressed transcripts, 38 up-regulated and 96 down-regulated (Figure 2B). As observed with the previous dataset, none of these differentially-expressed transcripts were linked to our list of interest. The next dataset, GSE70689, was used to observe acute transcriptional changes in mice microglia exposed to in vitro THC (10 μM for 2 h) (Juknat et al., 2013). Here, we performed subsetting of microglia exposed to either THC or vehicle, excluding samples exposed to cannabidiol. Using such a method, 60 transcripts were found to be significantly differentially expressed, 2 up-regulated and 58 downregulated (Figure 2C). None of our transcripts of interest were found to be significantly dysregulated in such a dataset. A further analysis revealed that none of the previously differentially expressed genes (outside our 13 transcripts of interest) overlapped in the three datasets (Venn diagram in Figure 2D). In addition to these observations, we also ran ORA Gene Ontology (GO) analyzes for each dataset. GSE189821 and GSE116813 returned no GO enrichment, while GSE70689 revealed many GO enrichments. These included transcripts involved in immune response, such as T cell activation, lymphocyte migration, cytokine-mediated signaling and ERK (extracellular signal-regulated kinases) cascades (Figure 2D). These surprising results indicate that our 13 transcripts of interest are only minimally involved in the response to in vivo THC exposure, at least in mice. However, in vitro exposure to THC likely activates microglia, in an attempt to prevent THC-induced inflammation. Representations of the three datasets (GSE189821, GSE116813 and GSE70689) and our 13 transcripts of interest can be found on Figure 2E. Note that GPR55 is the only transcript of interest that was significantly upregulated after chronic in vivo THC exposure in mice, and only in the VTA (Figure 2E). This suggests that our 13 transcripts of interest were not affected by THC exposure, at least under these specific conditions, except for GPR55. GSEA analyzes did not reveal any biological processes of interest (not shown).

Figure 2

3.3 TPRV1 is significantly downregulated in the epileptic tissue of patients with temporal lobe epilepsy

Our next step was to investigate the eCB system with regard to temporal lobe epilepsy, in light of extensive evidence of the links between the two (Alger, 2004; Kwan Cheung et al., 2019; Bilbao and Spanagel, 2022; Martínez-Aguirre et al., 2022; Sugaya and Kano, 2022). Three datasets were thus investigated: GSE77578 (Srivastava et al., 2018), GSE190451 (Chen et al., 2023) and GSE74150 (Damasceno et al., 2020). In GSE77578, 56 hippocampal samples were compared from epileptic or control mice. Here, we only compared control mice to epileptic mice (temporal lobe epilepsy), excluding mice treated with an antiepileptic drug. We found 129 differentially expressed transcripts in such a model of temporal lobe epilepsy, with 43 transcripts up-regulated and 86 down-regulated (Figure 3A). In GSE190451, 6 neocortical samples were compared between patients with temporal lobe epilepsy (n = 3) and aged-matched non-epileptic controls (n = 3), an analysis that differs slightly from the published article, in which the authors analyzed a total of 8 samples (Chen et al., 2023). We found a total of 1873 differentially expressed transcripts in temporal lobe epilepsy, 833 up-regulated and 1,040 down-regulated (Figure 3B). In GSE74150, an animal model of genetic epilepsy was used (Wistar audiogenic rat). Here, we found 67 differentially expressed transcripts, 53 down-regulated and 14 up-regulated (Figure 3C). Up- and down-regulated transcripts in the three datasets presented some overlap (Figure 3D, Supplementary Table S3), showing that transcripts involved in, or responding to, epilepsy seem to be found across species and models. Several gene ontology enrichments (ORA) were detected in the first two datasets, ranging from regulation of neurophysiology to neuronal death. Interestingly, 37 biological process terms were found to be common across the two first datasets (Supplementary Table S4), such as synapse organization, neuron death, myelination and plasma membrane organization (Figure 3C). GSEA analyzes did not reveal any biological processes of interest (not shown). Out of our 13 transcripts of interest, TRPV1 was found to be significantly down-regulated in the neocortex of patients with temporal lobe epilepsy (Figure 3D). These results confirm, once again, that TRPV1 plays an important role in the pathophysiology of epilepsy, as observed previously (Jia et al., 2015; Saffarzadeh et al., 2015; Kong et al., 2019; Wang et al., 2019; Lazarini-Lopes et al., 2022).

Figure 3

3.4 The eCB system interactome

Using the STRING database (Szklarczyk et al., 2019), containing both known and predicted protein interactions, we first fed into the database single protein names (from our list of 13 transcripts/proteins). Analysis of the immediate interactomes revealed great disparities in terms of protein overlap in pair-wise comparisons (Figure 4A). Interestingly, great overlap was observed with CB1R, CB2R and ADGRF1 (GPR110). We also noted great protein overlap in the interactomes of DAGLA, NAPE-PLD and ABHD6. As we expanded the interactome searches, similar conclusions were reached (Figure 4A). Surprisingly, TRPV1 and ALOX15B presented the lowest overlap in all levels, whether immediate, intermediated or extended interactomes (Figure 4A). We next scanned the database using our 13 proteins of interest by bulk feeding all 13 proteins into the search equation. When examining the immediate interactome of these 13 proteins, only the same 13 proteins were returned (Figure 4B), providing no additional information of interest. However, when the interactome search was expanded to an intermediate level, new proteins appeared, including GNAI1, GNAI2, GNAQ, and TRMP8 (Figure 4B). When the search was expanded to another level (extended interactome), new protein interactions arose, such as CALM3, CALML3, CALML5, CALML6 and TRPA1 (Figure 4B). Our final step was to apply k means clustering using a fixed number of clusters (k = 3). Very surprisingly, 4 proteins clustered outside of the main cluster containing the proteins of interest: TRPV1, PLC1B, ALOX15B and PLA2G4A (Figure 4C). We also noted that TRPM8 and GATB clustered with the majority of our proteins of interest, via FAAH, acting as a link within the extended interactome (Figure 4C).

Figure 4

3.5 Paralogues of human CBRs

Our next step was to try to determine how CB1R (the protein produced from the CNR1 gene) evolved within the human proteome. Using GeneCards (Stelzer et al., 2016) with Homo sapiens CB1R as input, 18 paralogues were returned. As expected, the closest paralogue was CNR2 (Figure 5A). CNR1 and CNR2 formed a distinct clade, with the next most closely related G-coupled protein receptors (GPCRs) being lysophosphatidic acid receptors (LPARs) and sphingosine-1-phosphate receptors (S1PRs). The furthest clade included melanocortin receptors 1–4 (MCRs) and GPR119. Inspection of the sequences revealed three conserved regions (Figure 5B), which could be divided into residues conserved within the clade containing CNR, LPAR and S1PR families (Figure 5B, blue) and residues only conserved within the CNR clade (Figure 5B, orange). Three conserved motifs (CMs) were identified: CM1, CM2 and CM3 are located within transmembrane (TM) 2, 3 and 7, respectively (Figure 5C). We then examined the location and function of these residues in the Homo sapiens CB1R structure bound to an anandamide analog (Krishna Kumar et al., 2023) and in the lysophosphatidic acid receptor structure bound to its native lysophosphatidic acid ligand (PDB: 7TD1) (Liu et al., 2022). The residues conserved within all three CNR, LPAR and S1PR families are largely found in the lower part of the receptor and include five of the six residues identified as being important in receptor activation (Figure 5D), as observed before (Hua et al., 2020). The ligand binding pocket in both receptors is formed by a large pocket of hydrophobic residues. Most of the residues conserved only within the CNR clade are located in this area and include some hydrophobic residues and H178 which forms a hydrogen bond with the terminal hydroxyl of the ligand (Figures 5D,E).

Figure 5

Out of 13 deuterostome species spanning from mammals to echinoderms (plus insects as protostomes), CB1R orthologues were found in placental mammals, marsupials, monotremes, birds, reptiles, amphibians, teleosts, elasmobranchs and Agnatha (Supplementary Table S5), but absent in cephalochordates, tunicates, echinoderms and insects. Furthermore, CB1R orthologues were also not found in plants, bacteria, Cnidaria and Porifera, except in Klebsiella (gram-negative bacteria) which presented 100% homology with the bird CB1R sequence (not shown), thus likely accounting for a database error in GenBank. When the same experiment was conducted with CB2R, we observed the presence of orthologues in placental mammals, marsupials, monotremes, birds, reptiles, amphibians and teleosts, the latter presenting two orthologues (Supplementary Table S5). Within the BLASTP database, unknown CBRs were found in cephalochordates and tunicates. These results suggest that the classical endocannabinoid receptors (CB1R and CB2R) are confined to a subgroup of deuterostomes, containing urochordates, cephalochordates and chordates, but absent from other deuterostomes, other animals, plants and prokaryotes, thus partly corroborating previous observations (McPartland et al., 2006; Elphick, 2012). Furthermore, it is suggested that a duplication event has occurred in teleosts since they diverged from other chordates. A further phylogenetic analysis revealed a clear division between CB1R and CB2R into two distinct clades (Figure 6A). We noted that most taxa presented copies of both proteins. There are two explanations for the lack of CNR2 paralogues in agnathans and elasmobranchs: (i) a duplication may have occurred since the divergence of other chordates and elasmobranchs or (ii) CNR2 paralogues were lost by both agnathans and elasmobranchs, the first scenario being by far the more likely explanation. Finally, we noted that the CBRs of Cephalochordates and Tunicates are very different, making it difficult to distinguish these two proteins from other related GPCRs. To address this, we performed multiple alignments with CB1R and CB2R protein sequences from the above-mentioned species. Our first alignment included cephalochordates and tunicates and led to great disruptions of the CMs (Figure 6B). The second alignment, which was performed without cephalochordates and tunicates, produced more defined, but yet not optimal, CMs (not shown). Our final alignment was performed without cephalochordates, tunicates and Danio rerio. This final approach produced robust identification of the three CMs (Figure 6C), thus corroborating the above results. Subsequently, we created homology models of the agnathan CB1R and cephalochordate unclassified CBR. Both the agnathan (Petromyzon marinus: ModelArchive ma-yvfgy) and cephalochordate (Branchiostoma floridae: ModelArchive ma-mi5fb) receptors contain the six residues important for receptor activation. However, the cephalochordate receptor lacks the majority of the residues involved in the AEA-binding pocket, including a histidine equivalent (H178) needed for hydrogen-bond formation with the terminal hydroxyl (Figures 6D,E).

Figure 6

3.6 Paralogues of human TRPVs

Our final analyzes focused on TRPV1, as this is the second eCB receptor in the central nervous system (together with CB1R). Using similar approaches, we first determined the phylogeny of the six Homo sapiens paralogues, using TRPV1 as the reference sequence. We observed that TRPV1 and TRPV2 diverged recently, which arose from a duplication event (Figure 7A). Likewise, divergence of TPRV3, TRPV5 and TRPV6 arose earlier, with TPRV5 and TRPV6 being the most recent split (Figure 7A). Two small CMs were noted throughout those 6 different proteins, including Y375-G376-P377, L382-Y383-D384-L385, L547-G548-W549 and R557-G558 (Figure 7B). Using Homo sapiens TRPV1 sequence as bait, we found that TPRV5 or TRPV6 were present across the most diverse range of taxa (Supplementary Table S6). These receptors are present in Bilatera but absent from Porifera, the latter presenting weak protein homology. These observations suggest that the TRPV family is specific to the Bilatera. Furthermore, TPRV1 and TRPV2 resulted from later duplication events in the common ancestor of jawed vertebrates and amniotes, respectively (Figure 7C). In addition, our phylogenetic analysis highlighted that eCB-binding receptors (TRPV1 and TRPV2) form one clade, which includes taxa from elasmobranchs onwards, with a more recent duplication forming TRPV2, found in amniotes only (Figure 7C), as suggested above (Supplementary Table S6). Alignments of TPRV1 and TRPV2 from these species correctly identified conserved residues involved in channel gating and selectivity, as well as ligand binding [anandamide, vanilloid, double-knot toxin (DkTx), Figure 7D]. Conserved residues were often found throughout species, indicative of the importance of these residues for TRPVs function. Finally, out of these 28 functional residues in TRPV1, only 8 were conserved in TRPV5, whether fully conserved or present with conservative substitutions (Figure 7E). Taken together, these results pinpoint a divergence of eCB receptors between jawed and jawless vertebrates, around 450 million years ago during the late Ordovician. These observations are consistent with previous studies (Elphick et al., 2003; McPartland, 2004; McPartland et al., 2006).

Figure 7

We examined the structures of TRPV1 and TRPV4 to identify residues important for ligand binding. TRPV1 can be activated by endogenous ligands such as AEA and exogenous ligands such as capsaicin, the component of chili peppers which creates a sensation of burning (Lam et al., 2005). Structural studies have identified the capsaicin binding site, known as the vanilloid-binding pocket (Nadezhdin et al., 2021; Kwon et al., 2022). There are no experimentally determined structures of the TRPV1 receptor bound to AEA. However, using a model of human TRPV1 based on the experimental rat structure, molecular dynamics and docking studies have identified two possible binding sites for anandamide (Muller et al., 2020; Morales et al., 2022), represented in Figure 8A. One site was the vanilloid-binding pocket (VBP) but the preferred site was a tunnel (AEA tunnel) involving interactions with Y554, Y555, Y487, D708, and N438. We docked AEA into the AEA tunnel using the recently experimentally determined structure of human TRPV1 (PDB: 8GFA) (Neuberger et al., 2023). One of nine docking solutions had a similar orientation to that shown in Morales et alia (Morales et al., 2022), with Y555, N438 and Y554 forming hydrogen bonds with the amide group hydrogen and the terminal hydroxyl group oxygen (Figure 8B, ModelArchive ma-dol4y). The five interacting residues are well conserved across species in TRPV1 and TRPV2. AEA itself does not directly bind to TPRV4, however this receptor is activated by an AEA metabolite 5′,6′-epoxyeicosatrienoic acid (5′,6′-EET) (Watanabe et al., 2003). Homo sapiens TRPV4 (PDB: 8TD1) (Nadezhdin et al., 2021) is missing equivalents to TRPV1 Y847 and Y555 which stabilize AEA with hydrogen bonds in the AEA tunnel. Simulations have suggested that residues K535, F549, Q550, Y591 and R594 of TRPV4 form the 5′,6′-EET binding pocket (Berna-Erro et al., 2017), as displayed in Figure 8C. Of these, Y591 is equivalent to Y554 in TRPV1 which forms a hydrogen bond with the terminal hydroxyl of AEA, but the remaining pocket is very different. We also examined differences in the vanilloid-binding pocket by superimposing capsaicin from the rat (Rattus norvegicus) TRPV2 complex structure (PDB: 8SLY) (Gochman et al., 2023) onto the Homo sapiens TRPV1 receptor (PDB: 8GFA). Whilst some important hydrophobic residues are conserved in the equivalent pocket in TRPV4, there are no equivalent residues for the hydrogen-bonding Y511 and T550 which were previously observed to be involved in propagating ligand-activated conformational changes (Figure 8D) (Kwon et al., 2022).

Figure 8

4 Discussion

The location of CB1Rs is restricted to pre-synaptic terminals of both excitatory and inhibitory neurons (Mátyás et al., 2008). The physiological roles of CB1R are to mediate mood (Sbarski and Akirav, 2018), memory (Borgan et al., 2019a), stress (Arévalo et al., 2001; Barna et al., 2004) and locomotion (Han et al., 2023). These functions were reviewed in 2020 (Cristino et al., 2020). Activation of CB1Rs induces synaptic plasticity via long-term depression (Wu et al., 2015). Indeed, following synaptic activation, endogenous endocannabinoids (AEA and 2-AG) will be released from the post-synaptic element to then diffuse to the pre-synaptic element, where both molecules will act as agonists at CB1 receptors. Such activation of CB1 receptors will reduce neurotransmitter release, ending the cycle of synapse activation. This mechanism is often referred to as endocannabinoid-mediated negative feedback. Interestingly, CB1Rs were found to be expressed in the nucleus accumbens of vervet monkeys (Kucera et al., 2018), thus highlighting the role of these receptors in the reward circuitry (Méndez-Díaz et al., 2019; Braunscheidel et al., 2022; Ceccarini et al., 2022; Han et al., 2023). Dysfunction of CB1Rs has been linked to a wide array of neurological disorders such as Parkinson’s Disease, Alzheimer’s Disease, Huntington Disease and multiple sclerosis [reviewed in Cristino et al. (2020)].

In contrast to CB1R, TRPV1 is expressed in both presynaptic (Tominaga et al., 1998; Medvedeva et al., 2008) and postsynaptic (Bae et al., 2004; Grueter et al., 2010) elements. The physiological roles of TRPV1 intertwined those of CB1R, such as memory and locomotion, but also other functions, such as food intake (Cristino et al., 2020). Extensive evidence suggests that TRPV1 can be linked to epilepsy (Jia et al., 2015; Saffarzadeh et al., 2015; Kong et al., 2019; Wang et al., 2019; Lazarini-Lopes et al., 2022), although causal determination remains unclear. Indeed, in naïve rodents, chronic exposure to the TRPV1 agonist capsaicin induces spontaneous seizures (Jia et al., 2015). On the other hand, pre-treatment with the TRPV1 antagonist capsazepine decreased seizure severity and mortality in pentylenetetrazol-induced epilepsy (GABAA antagonist) (Jia et al., 2015). Another study has observed altered electrophysiological properties in microglial cells following capsaicin, which was hypothesized as driving hyperthermia-induced seizures (Kong et al., 2019), although direct experimental confirmation has not been performed. Increased expression of TRPV1 was also observed after the successful induction of temporal lope epilepsy through kindling (audiogenic epilepsy) in rats (Lazarini-Lopes et al., 2022). Similar results were observed in another study using a model based on pilocarpine exposure (Saffarzadeh et al., 2015). Another elegant confirmation of the link between epilepsy and TRPV1 was observed in animals lacking TRPV1 (knock-out). Indeed, compared to wild-type mice, mice lacking TRPV1 presented decreased seizure latency, duration and severity (Wang et al., 2019). In agreement with these observations, our current transcriptomic results report a downregulation of TPRV1 in temporal lobe epilepsy. Indeed, we found a significant downregulation of TRPV1 mRNA in the neocortex of patients with temporal lobe epilepsy compared to control tissues. As mentioned above, preclinical studies have shown that agonism of TRPV1, using capsaicin, can induce tonic–clonic seizures, while TPRV1 antagonism reduced mortality in animals in which seizures were induced (Jia et al., 2015). These results were also observed in knock-out animals for TRPV1 (Jia et al., 2015). However, contrasting results were observed in a 2020 study, in which capsaicin inhibited epileptiform activities in the prefrontal cortex (Pasierski and Szulczyk, 2020). The same authors also witnessed decreased neuronal excitation, thus explaining the above-mentioned effects (Pasierski and Szulczyk, 2020), which are mediated by decreased sodium currents. These were not isolated results, as other studies have also concluded that capsaicin possesses anti-convulsive properties in rats (Abdel-Salam et al., 2020) and mice (Lee et al., 2011; Pezzoli et al., 2014). Interestingly, hippocampal expression of CB1Rs was also found to be significantly decreased in patients suffering from epilepsy (Ludányi et al., 2008). Thus, in light of these previous results and those from the present study, TRPV1 could be of interest for developing new therapeutics for patients with refractory epilepsy.

The role of the eCB system in neurodevelopmental disorders has been described previously. Indeed, many studies observed a direct relationship between altered CB1R availability in schizophrenia (Dean et al., 2001; Zavitsanou et al., 2004; Newell et al., 2006; Eggan et al., 2008; Urigüen et al., 2009; Wong et al., 2010; Dalton et al., 2011; Jenko et al., 2012; Ceccarini et al., 2013; Ranganathan et al., 2016; Muguruza et al., 2019; Bloch Priel et al., 2023) or during psychosis (Borgan et al., 2019b, 2021), which is also corroborated by animal models of schizophrenia (Szűcs et al., 2016; Stark et al., 2022). Furthermore, a link between THC exposure and schizophrenia has been evidenced in humans (Hjorthøj et al., 2021) and animal models (Rodríguez et al., 2017). Compared to control patients, levels of AEA are increased in the blood (De Marchi et al., 2003; Ibarra-Lecue et al., 2022) and cerebrospinal fluid (Giuffrida et al., 2004; Leweke et al., 2007; Minichino et al., 2019) in patients with schizophrenia. In addition, exposure to THC alters AEA levels (Leweke et al., 2007; Seillier et al., 2020; Ibarra-Lecue et al., 2022). Furthermore, lower brain FAAH availability was observed in adolescent chronic cannabis users (Jacobson et al., 2021). Targeting FAAH has proven an efficient treatment to decrease cannabis use and withdrawal (D’Souza et al., 2019). During psychosis, increased levels of DAGL and NAPE are reported in mononuclear cells (peripheral blood) compared to control patients (Bioque et al., 2013), two enzymes controlling endocannabinoid production. In addition, such a study also reported increased FAAH levels during psychosis (Bioque et al., 2013). Altogether, these studies provide ample evidence for a link between the eCB system and neurodevelopmental disorders. Recent studies reviewed the relationship of the eCB in several neurological pathologies such as depression (Gallego-Landin et al., 2021), schizophrenia (Garani et al., 2021), autism (Su et al., 2021) and substance use disorder (Navarro et al., 2022). In the current study, our transcriptomic analysis found that GPR55 is upregulated in the ventral tegmental area of adolescent mice following THC exposure, even after a long washout period (2 weeks). These results can be explained by the fact that THC possesses psychoactive effects (Ameri et al., 1999) and that the ventral tegmental area is known to respond to THC (French et al., 1997; Lupica et al., 2004; Ostlund et al., 2022). However, what remains surprising is the fact that GPR55, and not CB1R, is upregulated after THC exposure, in light of two studies demonstrating loss of CB1R-dependent plasticity after chronic THC exposure (Friend et al., 2017; Ostlund et al., 2022). The answer appears to lie within the capacity of GPR55 to respond to THC (Ryberg et al., 2007; Lauckner et al., 2008), which has been the reason why GPR55 is classed as an endocannabinoid receptor. In fact, THC, AEA, methanandamide (stable metabolic analog of AEA) and JWH015 (cannabinoid agonist) all seem to elicit intracellular responses in GPR55-transfected HEK cells (Lauckner et al., 2008), thus proving a link between GPR55 and THC. Furthermore, CBD is also known to abolish GPR55-mediated synaptic events (Sylantyev et al., 2013), highlighting again a strong interaction between GPR55 and the eCB system.

The two phytocannabinoids THC and CBD are known to bind to CB1R (Jakowiecki and Filipek, 2016; Hua et al., 2017; Chung et al., 2019), since the only chemical difference between these two compounds is a bond disconnection in CBD (Shao et al., 2016). Functionally, CBD has a low affinity for the CB1R (McPartland et al., 2007) and is described as an antagonist at CB1R (Thomas et al., 2007). However, evidence suggests that CBD rather behaves as a negative allosteric modulator of THC (Laprairie et al., 2015; Chung et al., 2019). Interestingly, CBD also binds to TRPV1-4 channels, as observed recently (Costa et al., 2004; Qin et al., 2008; Iannotti et al., 2014; Muller et al., 2019; Starkus et al., 2019; Galaj et al., 2020; Etemad et al., 2022; Landucci et al., 2022). These observations were also made outside of the central nervous system (De Petrocellis et al., 2012; Nabissi et al., 2013; de la Harpe et al., 2022).

We also observed that the amino acid sequence of TRPV1 in Danio rerio is missing Y554 and Y555, which have been shown to allow AEA binding (Muller et al., 2020). Instead, TRPV1 in zebrafish possesses A554 and A555. The switch from polar amino acids (Tyr-Tyr) to non-polar amino acids (Ala-Ala) in the TRPV1 orthologue of zebrafish does not seem to impact heat detection, which is the main function of TRPV1 (Zheng and Wen, 2019). Despite these amino acids substitutions, in vitro and in vivo observations confirmed that TRPV1 is triggered above 25°C in zebrafish (Gau et al., 2013). In contrast, heat activation of TRPV1 in mammals occurs above 42°C (Caterina et al., 1997), which could be explained by slightly different structural conformations following amino acid substitutions. Interestingly, TRPV3 seems to possess opposite heat-detection functions between humans and reptiles (Xenopus tropicalis), as observed before (Saito et al., 2011), with the receptor activated by cold temperatures ( 16°C). Behaviorally, TRPA1 seems to be the receptor mediating avoidance of dangerous temperatures in frogs (Saito et al., 2022). Similar to our findings, Saito et alia also reported that, in vertebrates, TRPV channels form very distinct clades (Saito et al., 2011). Furthermore, TRPV1 and TRPA1 seem to have been co-expressed in an early vertebrate common ancestor (Saito et al., 2012).

The link between TRPV1 and heat can be experienced when eating hot chili peppers, which dose-dependently trigger receptor activation (Mohapatra and Nau, 2003; Voets et al., 2004; Neelands et al., 2005; Vetter et al., 2006; Vasquez et al., 2020; Matsuyama et al., 2021). THC can bind TRPV2-4, TRPA1 and TRPM8, but not TRPV1 (Muller et al., 2019). This seems very different from CBD, which has been shown to bind to TRPV1 (see above). What remains fascinating is the fact that TRPV channels are involved in sensory perception (chemical nociception, heat detection) in the peripheral nervous system, whilst TRPV1 in the central nervous system is involved in synaptic plasticity, neurite growth and signaling (Moran et al., 2004).

In terms of evolution, our results suggest that the eCB system appeared during the late Ordovician. McPartland and colleagues studied the evolution of all components of the eCB system in a study published in 2006 (McPartland et al., 2006). In their article, authors evidenced that TRPV1 and GPR55 are limited to mammals, whilst CB2R and DAGL were found amongst vertebrates (McPartland et al., 2006). Other components, such as MAGL, are found in chordates, while all eukaryotes express FAAH (McPartland et al., 2006). Thus, it can be safely hypothesized that the origins of the eCB appeared in a common eukaryotic ancestor. Such a study also suggested that ligand-metabolizing enzymes of the eCB system evolved prior to functional cannabinoid receptors (McPartland et al., 2006). In the absence of these receptors, such as CB1R, CB2R and TRPV1, endocannabinoids, acting as signaling molecules, such as AEA and 2-AG, might possess other functional properties (McPartland et al., 2006). A likely explanation is modulating physiological excitability in the nervous system (McPartland, 2004). 2-AG is expressed in insects and, apart from signaling properties, 2-AG might serve as a deterrent to predators, such as vertebrates, who do express functional receptors to 2-AG, such as TRPV1. This mechanism was likely driven by convergent evolution. Interestingly, this feeding-deterrent mechanism was also suggested for plants (McPartland et al., 2006), expressing endocannabinoid-like compounds but without eCB receptors. Intertwined evolution of receptors, enzymes and ligands of the eCB system was likely driven by convergence, divergence and parallel evolution (McPartland and Guy, 2004). To the best of our knowledge, pinpointing an exact timeframe for such mechanisms is currently impossible. However, appearance of the eCB in the late Ordovician seems to be the accepted consensus, and our current study further corroborates these previous suggestions.

Indeed, recent studies have observed brain-like structures in fossils from the Cambrian-Ordovician eras (Dong et al., 2022; Ortega-Hernández et al., 2022; Pates et al., 2022), which seem to support the idea of an earlier origin of the nervous system, likely occurring before the appearance of annelids (Parry and Caron, 2019). This pinpoints to the nervous system originating just before 500 million years ago, with the eCB system likely appearing 50 million years after the former. Further evidence has been gathered in invertebrates (annelids and insects), in which observations of eCB-like receptors were discovered (Stefano et al., 1997; Elphick, 1998, 2012; McPartland et al., 2006). Some components of the eCB system have been found in species much older than annelids. Indeed, FAAH, which is responsible for metabolism of AEA, has been identified in the moss Physcomitrella patens, with a high catabolic activity for AEA (Haq and Kilaru, 2020). Furthermore, 9 paralogues were identified in this species, with some closely related to rat FAAH, while others closely related to plant FAAH (Haq and Kilaru, 2020). It is now accepted that the nervous system appeared in bilaterians ancestors. Recovered fossils suggest that the nervous system was already present 540 million years ago (McPartland and Guy, 2004; Budd, 2013; Heger et al., 2020). Furthermore, several studies confirmed the presence of a nervous system in bilaterians (Adey, 1951; Holland et al., 2013; Yang et al., 2013, 2016; Martín-Durán et al., 2018), sponges (Leys, 2015; Musser et al., 2021; Yañez-Guerra et al., 2022) or ctenophores (Jékely et al., 2015; Moroz, 2015; Burkhardt et al., 2023). Interestingly, the oldest bilaterians fossils recovered so far range between 540 and 585 million years (Martin et al., 2000; Pecoits et al., 2012; Chen et al., 2019; Evans et al., 2020). Altogether, it can be safely assumed that the nervous system appeared before such a timeframe.

To conclude, we provide here new evidence that GPR55 is upregulated in the ventral tegmental area of mice following THC exposure and that TPRV1 is downregulated in the epileptogenic focus of patients with temporal lobe epilepsy. These results further support previous observations linking the eCB system to epilepsy, cannabis and neurodevelopmental disorders. We also report here an updated analysis of the evolution of CB1R and TRPV1, two major actors in the eCB system.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

AC: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing. MW: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing. KM: Formal analysis, Investigation, Writing – review & editing. MDM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Software, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the University of Worcester, grants SSE2206 (to MDM) and ROFF22234MDM (to MDM).

Acknowledgments

The authors would like to thank the University of Worcester for support with funding. In addition, the authors would like to thank Pr. Eleanor Bradley, Dr. Derek McDougall, and Dr. Steven Coles.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fninf.2024.1294939/full#supplementary-material

References

  • 1

    Abdel-SalamO. M. E.SleemA. A.SayedM. A. E. B. M.YounessE. R.ShaffieN. (2020). Capsaicin exerts anti-convulsant and neuroprotective effects in Pentylenetetrazole-induced seizures. Neurochem. Res.45, 10451061. doi: 10.1007/s11064-020-02979-3

  • 2

    AdeyW. R. (1951). The nervous system of the earthworm Megascolex. J. Comp. Neurol.94, 57103. doi: 10.1002/cne.900940104

  • 3

    AgrawalA.WetherillL.DickD. M.XueiX.HinrichsA.HesselbrockV.et al. (2009). Evidence for association between polymorphisms in the cannabinoid receptor 1 (CNR1) gene and cannabis dependence. Am. J. Med. Genet. B Neuropsychiatr. Genet.150B, 736740. doi: 10.1002/ajmg.b.30881

  • 4

    AlgerB. E. (2004). Endocannabinoids and their implications for epilepsy. Epilepsy Curr.4, 169173. doi: 10.1111/j.1535-7597.2004.04501.x

  • 5

    AmeriA.WilhelmA.SimmetT. (1999). Effects of the endogeneous cannabinoid, anandamide, on neuronal activity in rat hippocampal slices. Br. J. Pharmacol.126, 18311839. doi: 10.1038/sj.bjp.0702478

  • 6

    AranA.EylonM.HarelM.PolianskiL.NemirovskiA.TepperS.et al. (2019). Lower circulating endocannabinoid levels in children with autism spectrum disorder. Mol. Autism.10:2. doi: 10.1186/s13229-019-0256-6

  • 7

    ArévaloC.de MiguelR.Hernández-TristánR. (2001). Cannabinoid effects on anxiety-related behaviours and hypothalamic neurotransmitters. Pharmacol. Biochem. Behav.70, 123131. doi: 10.1016/s0091-3057(01)00578-0

  • 8

    Arias-GaguancelaO.HerrellE.AzizM.ChapmanK. D. (2023). Two legume fatty acid amide hydrolase isoforms with distinct preferences for microbial- and plant-derived acylamides. Sci. Rep.13:7486. doi: 10.1038/s41598-023-34754-z

  • 9

    BaeY. C.OhJ. M.HwangS. J.ShigenagaY.ValtschanoffJ. G. (2004). Expression of vanilloid receptor TRPV1 in the rat trigeminal sensory nuclei. J. Comp. Neurol.478, 6271. doi: 10.1002/cne.20272

  • 10

    BainbridgeM. N.MazumderA.OgasawaraD.Abou JamraR.BernardG.BertiniE.et al. (2022). Endocannabinoid dysfunction in neurological disease: neuro-ocular DAGLA-related syndrome. Brain145, 33833390. doi: 10.1093/brain/awac223

  • 11

    BarchelD.StolarO.De-HaanT.Ziv-BaranT.SabanN.FuchsD. O.et al. (2018). Oral Cannabidiol use in children with autism Spectrum disorder to treat related symptoms and co-morbidities. Front. Pharmacol.9:1521. doi: 10.3389/fphar.2018.01521

  • 12

    BarnaI.ZelenaD.ArszovszkiA. C.LedentC. (2004). The role of endogenous cannabinoids in the hypothalamo-pituitary-adrenal axis regulation: in vivo and in vitro studies in CB1 receptor knockout mice. Life Sci.75, 29592970. doi: 10.1016/j.lfs.2004.06.006

  • 13

    BasavarajappaB. S. (2007). Critical enzymes involved in endocannabinoid metabolism. Protein Pept. Lett.14, 237246. doi: 10.2174/092986607780090829

  • 14

    Berna-ErroA.Izquierdo-SerraM.SepúlvedaR. V.Rubio-MoscardoF.Doñate-MaciánP.SerraS. A.et al. (2017). Structural determinants of 5′,6′-epoxyeicosatrienoic acid binding to and activation of TRPV4 channel. Sci. Rep.7:10522. doi: 10.1038/s41598-017-11274-1

  • 15

    BilbaoA.SpanagelR. (2022). Medical cannabinoids: a pharmacology-based systematic review and meta-analysis for all relevant medical indications. BMC Med.20:259. doi: 10.1186/s12916-022-02459-1

  • 16

    BioqueM.García-BuenoB.MacDowellK. S.MeseguerA.SaizP. A.ParelladaM.et al. (2013). Peripheral endocannabinoid system dysregulation in first-episode psychosis. Neuropsychopharmacology38, 25682577. doi: 10.1038/npp.2013.165

  • 17

    Bloch PrielS.YitzhakyA.GurwitzD.HertzbergL. (2023). Cannabinoid receptor gene CNR1 is downregulated in subcortical brain samples and upregulated in blood samples of individuals with schizophrenia: a participant data systematic meta-analysis. Eur. J. Neurosci.58, 35403554. doi: 10.1111/ejn.16122

  • 18

    BorganF.BeckK.ButlerE.McCutcheonR.VeroneseM.VernonA.et al. (2019a). The effects of cannabinoid 1 receptor compounds on memory: a meta-analysis and systematic review across species. Psychopharmacology236, 32573270. doi: 10.1007/s00213-019-05283-3

  • 19

    BorganF.LaurikainenH.VeroneseM.MarquesT. R.Haaparanta-SolinM.SolinO.et al. (2019b). In vivo availability of cannabinoid 1 receptor levels in patients with first-episode psychosis. JAMA Psychiatry76, 10741084. doi: 10.1001/jamapsychiatry.2019.1427

  • 20

    BorganF.VeroneseM.Reis MarquesT.LythgoeD. J.HowesO. (2021). Association between cannabinoid 1 receptor availability and glutamate levels in healthy controls and drug-free patients with first episode psychosis: a multi-modal PET and 1H-MRS study. Eur. Arch. Psychiatry Clin. Neurosci.271, 677687. doi: 10.1007/s00406-020-01191-2

  • 21

    BraunscheidelK. M.OkasM. P.FlorescoS. B.WoodwardJ. J. (2022). Cannabinoid receptor type 1 antagonists alter aspects of risk/reward decision making independent of toluene-mediated effects. Psychopharmacology239, 13371347. doi: 10.1007/s00213-021-05914-8

  • 22

    BuddG. E. (2013). The Cambrian explosion: the reconstruction of animal biodiversity.— by Douglas H. Erwin and James W. Valentine. Syst. Biol.62, 915917. doi: 10.1093/sysbio/syt043

  • 23

    BurkhardtP.ColgrenJ.MedhusA.DigelL.NaumannB.Soto-AngelJ. J.et al. (2023). Syncytial nerve net in a ctenophore adds insights on the evolution of nervous systems. Science380, 293297. doi: 10.1126/science.ade5645

  • 24

    CastilloP. E.YountsT. J.ChávezA. E.HashimotodaniY. (2012). Endocannabinoid signaling and synaptic function. Neuron76, 7081. doi: 10.1016/j.neuron.2012.09.020

  • 25

    CaterinaM. J.SchumacherM. A.TominagaM.RosenT. A.LevineJ. D.JuliusD. (1997). The capsaicin receptor: a heat-activated ion channel in the pain pathway. Nature389, 816824. doi: 10.1038/39807

  • 26

    CeccariniJ.De HertM.Van WinkelR.PeuskensJ.BormansG.KranasterL.et al. (2013). Increased ventral striatal CB1 receptor binding is related to negative symptoms in drug-free patients with schizophrenia. NeuroImage79, 304312. doi: 10.1016/j.neuroimage.2013.04.052

  • 27

    CeccariniJ.KooleM.Van LaereK. (2022). Cannabinoid receptor availability modulates the magnitude of dopamine release in vivo in the human reward system: a preliminary multitracer positron emission tomography study. Addict. Biol.27:e13167. doi: 10.1111/adb.13167

  • 28

    ChenZ.-P.WangS.ZhaoX.FangW.WangZ.YeH.et al. (2023). Lipid-accumulated reactive astrocytes promote disease progression in epilepsy. Nat. Neurosci.26, 542554. doi: 10.1038/s41593-023-01288-6

  • 29

    ChenZ.ZhouC.YuanX.XiaoS. (2019). Death march of a segmented and trilobate bilaterian elucidates early animal evolution. Nature573, 412415. doi: 10.1038/s41586-019-1522-7

  • 30

    ChungH.FierroA.Pessoa-MahanaC. D. (2019). Cannabidiol binding and negative allosteric modulation at the cannabinoid type 1 receptor in the presence of delta-9-tetrahydrocannabinol: an in silico study. PLoS One14:e0220025. doi: 10.1371/journal.pone.0220025

  • 31

    CloughE.BarrettT. (2016). The gene expression omnibus database. Methods Mol. Biol.1418, 93110. doi: 10.1007/978-1-4939-3578-9_5

  • 32

    CostaB.GiagnoniG.FrankeC.TrovatoA. E.ColleoniM. (2004). Vanilloid TRPV1 receptor mediates the antihyperalgesic effect of the nonpsychoactive cannabinoid, cannabidiol, in a rat model of acute inflammation. Br. J. Pharmacol.143, 247250. doi: 10.1038/sj.bjp.0705920

  • 33

    CristinoL.BisognoT.Di MarzoV. (2020). Cannabinoids and the expanded endocannabinoid system in neurological disorders. Nat. Rev. Neurol.16, 929. doi: 10.1038/s41582-019-0284-z

  • 34

    D’SouzaD. C.Cortes-BrionesJ.CreaturaG.BluezG.ThurnauerH.DeasoE.et al. (2019). Efficacy and safety of a fatty acid amide hydrolase inhibitor (PF-04457845) in the treatment of cannabis withdrawal and dependence in men: a double-blind, placebo-controlled, parallel group, phase 2a single-site randomised controlled trial. Lancet Psychiatry6, 3545. doi: 10.1016/S2215-0366(18)30427-9

  • 35

    DalmaijerE. S.NordC. L.AstleD. E. (2022). Statistical power for cluster analysis. BMC Bioinformatics23:205. doi: 10.1186/s12859-022-04675-1

  • 36

    DaltonV. S.LongL. E.WeickertC. S.ZavitsanouK. (2011). Paranoid schizophrenia is characterized by increased CB1 receptor binding in the dorsolateral prefrontal cortex. Neuropsychopharmacology36, 16201630. doi: 10.1038/npp.2011.43

  • 37

    DamascenoS.Gómez-NietoR.Garcia-CairascoN.Herrero-TurriónM. J.MarínF.LopézD. E. (2020). Top common differentially expressed genes in the epileptogenic nucleus of two strains of rodents susceptible to Audiogenic seizures: WAR and GASH/Sal. Front. Neurol.11:33. doi: 10.3389/fneur.2020.00033

  • 38

    de la HarpeA.BeukesN.FrostC. L. (2022). CBD activation of TRPV1 induces oxidative signaling and subsequent ER stress in breast cancer cell lines. Biotechnol. Appl. Biochem.69, 420430. doi: 10.1002/bab.2119

  • 39

    De MarchiN.De PetrocellisL.OrlandoP.DanieleF.FezzaF.Di MarzoV. (2003). Endocannabinoid signalling in the blood of patients with schizophrenia. Lipids Health Dis.2:5. doi: 10.1186/1476-511X-2-5

  • 40

    de Miguel-YanesJ. M.ManningA. K.ShraderP.McAteerJ. B.GoelA.HamstenA.et al. (2011). Variants at the endocannabinoid receptor CB1 gene (CNR1) and insulin sensitivity, type 2 diabetes, and coronary heart disease. Obesity (Silver Spring)19, 20312037. doi: 10.1038/oby.2011.135

  • 41

    De PetrocellisL.OrlandoP.MorielloA. S.AvielloG.StottC.IzzoA. A.et al. (2012). Cannabinoid actions at TRPV channels: effects on TRPV3 and TRPV4 and their potential relevance to gastrointestinal inflammation. Acta Physiol (Oxf.)204, 255266. doi: 10.1111/j.1748-1716.2011.02338.x

  • 42

    DeanB.SundramS.BradburyR.ScarrE.CopolovD. (2001). Studies on [3H]CP-55940 binding in the human central nervous system: regional specific changes in density of cannabinoid-1 receptors associated with schizophrenia and cannabis use. Neuroscience103, 915. doi: 10.1016/S0306-4522(00)00552-2

  • 43

    DongX.DuanB.LiuJ.DonoghueP. C. J. (2022). Internal anatomy of a fossilized embryonic stage of the Cambrian-Ordovician scalidophoran Markuelia. R. Soc. Open Sci.9:220115. doi: 10.1098/rsos.220115

  • 44

    EdgarR.DomrachevM.LashA. E. (2002). Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res.30, 207210. doi: 10.1093/nar/30.1.207

  • 45

    EgertováM.ElphickM. R. (2007). Localization of CiCBR in the invertebrate chordate Ciona intestinalis: evidence of an ancient role for cannabinoid receptors as axonal regulators of neuronal signalling. J. Comp. Neurol.502, 660672. doi: 10.1002/cne.21331

  • 46

    EgganS. M.HashimotoT.LewisD. A. (2008). Reduced cortical cannabinoid 1 receptor messenger RNA and protein expression in schizophrenia. Arch. Gen. Psychiatry65, 772784. doi: 10.1001/archpsyc.65.7.772

  • 47

    ElphickM. R. (1998). An invertebrate G-protein coupled receptor is a chimeric cannabinoid/melanocortin receptor. Brain Res.780, 168171. doi: 10.1016/S0006-8993(97)01297-3

  • 48

    ElphickM. R. (2012). The evolution and comparative neurobiology of endocannabinoid signalling. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci.367, 32013215. doi: 10.1098/rstb.2011.0394

  • 49

    ElphickM. R.SatouY.SatohN. (2003). The invertebrate ancestry of endocannabinoid signalling: an orthologue of vertebrate cannabinoid receptors in the urochordate Ciona intestinalis. Gene302, 95101. doi: 10.1016/s0378-1119(02)01094-6

  • 50

    EtemadL.KarimiG.AlaviM. S.RoohbakhshA. (2022). Pharmacological effects of cannabidiol by transient receptor potential channels. Life Sci.300:120582. doi: 10.1016/j.lfs.2022.120582

  • 51

    EvansS. D.HughesI. V.GehlingJ. G.DroserM. L. (2020). Discovery of the oldest bilaterian from the Ediacaran of South Australia. Proc. Natl. Acad. Sci.117, 78457850. doi: 10.1073/pnas.2001045117

  • 52

    FrenchE. D.DillonK.WuX. (1997). Cannabinoids excite dopamine neurons in the ventral tegmentum and substantia nigra. Neuroreport8, 649652. doi: 10.1097/00001756-199702100-00014

  • 53

    FriendL.WeedJ.SandovalP.NuferT.OstlundI.EdwardsJ. G. (2017). CB1-dependent Long-term depression in ventral tegmental area GABA neurons: a novel target for marijuana. J. Neurosci.37, 1094310954. doi: 10.1523/JNEUROSCI.0190-17.2017

  • 54

    GalajE.BiG.-H.YangH.-J.XiZ.-X. (2020). Cannabidiol attenuates the rewarding effects of cocaine in rats by CB2, 5-HT1A and TRPV1 receptor mechanisms. Neuropharmacology167:107740. doi: 10.1016/j.neuropharm.2019.107740

  • 55

    Gallego-LandinI.García-BaosA.Castro-ZavalaA.ValverdeO. (2021). Reviewing the role of the endocannabinoid system in the pathophysiology of depression. Front. Pharmacol.12:762738. doi: 10.3389/fphar.2021.762738

  • 56

    GaraniR.WattsJ. J.MizrahiR. (2021). Endocannabinoid system in psychotic and mood disorders, a review of human studies. Prog. Neuro-Psychopharmacol. Biol. Psychiatry106:110096. doi: 10.1016/j.pnpbp.2020.110096

  • 57

    GascuelO. (1997). BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol. Biol. Evol.14, 685695. doi: 10.1093/oxfordjournals.molbev.a025808

  • 58

    GauP.PoonJ.Ufret-VincentyC.SnelsonC. D.GordonS. E.RaibleD. W.et al. (2013). The zebrafish Ortholog of TRPV1 is required for heat-induced locomotion. J. Neurosci.33, 52495260. doi: 10.1523/JNEUROSCI.5403-12.2013

  • 59

    GiuffridaA.LewekeF. M.GerthC. W.SchreiberD.KoetheD.FaulhaberJ.et al. (2004). Cerebrospinal anandamide levels are elevated in acute schizophrenia and are inversely correlated with psychotic symptoms. Neuropsychopharmacology29, 21082114. doi: 10.1038/sj.npp.1300558

  • 60

    GochmanA.TanX.-F.BaeC.ChenH.SwartzK. J.Jara-OsegueraA. (2023). Cannabidiol sensitizes TRPV2 channels to activation by 2-APB. elife12:e86166. doi: 10.7554/eLife.86166

  • 61

    Gomis-GonzálezM.Busquets-GarciaA.MatuteC.MaldonadoR.MatoS.OzaitaA. (2016). Possible therapeutic doses of cannabinoid type 1 receptor antagonist reverses key alterations in fragile X syndrome mouse model. Genes (Basel)7:56. doi: 10.3390/genes7090056

  • 62

    GrueterB. A.BrasnjoG.MalenkaR. C. (2010). Postsynaptic TRPV1 triggers cell type-specific long-term depression in the nucleus accumbens. Nat. Neurosci.13, 15191525. doi: 10.1038/nn.2685

  • 63

    HabibA. M.OkorokovA. L.HillM. N.BrasJ. T.LeeM.-C.LiS.et al. (2019). Microdeletion in a FAAH pseudogene identified in a patient with high anandamide concentrations and pain insensitivity. Br. J. Anaesth.123, e249e253. doi: 10.1016/j.bja.2019.02.019

  • 64

    HanX.LiangY.HempelB.JordanC. J.ShenH.BiG.-H.et al. (2023). Cannabinoid CB1 receptors are expressed in a subset of dopamine neurons and underlie cannabinoid-induced aversion, Hypoactivity, and anxiolytic effects in mice. J. Neurosci.43, 373385. doi: 10.1523/JNEUROSCI.1493-22.2022

  • 65

    HaqI.KilaruA. (2020). An endocannabinoid catabolic enzyme FAAH and its paralogs in an early land plant reveal evolutionary and functional relationship with eukaryotic orthologs. Sci. Rep.10:3115. doi: 10.1038/s41598-020-59948-7

  • 66

    HegerP.ZhengW.RottmannA.PanfilioK. A.WieheT. (2020). The genetic factors of bilaterian evolution. elife9:e45530. doi: 10.7554/eLife.45530

  • 67

    HermanA. I.KranzlerH. R.CubellsJ. F.GelernterJ.CovaultJ. (2006). Association study of the CNR1 gene exon 3 alternative promoter region polymorphisms and substance dependence. Am. J. Med. Genet. B Neuropsychiatr. Genet.141B, 499503. doi: 10.1002/ajmg.b.30325

  • 68

    HjorthøjC.PosseltC. M.NordentoftM. (2021). Development over time of the population-attributable risk fraction for Cannabis use disorder in schizophrenia in Denmark. JAMA Psychiatry78, 10131019. doi: 10.1001/jamapsychiatry.2021.1471

  • 69

    HollandL. Z.CarvalhoJ. E.EscrivaH.LaudetV.SchubertM.ShimeldS. M.et al. (2013). Evolution of bilaterian central nervous systems: a single origin?EvoDevo4:27. doi: 10.1186/2041-9139-4-27

  • 70

    HuaT.LiX.WuL.Iliopoulos-TsoutsouvasC.WangY.WuM.et al. (2020). Activation and signaling mechanism revealed by cannabinoid receptor-Gi complex structures. Cell180, 655665.e18. doi: 10.1016/j.cell.2020.01.008

  • 71

    HuaT.VemuriK.NikasS. P.LaprairieR. B.WuY.QuL.et al. (2017). Crystal structures of agonist-bound human cannabinoid receptor CB1. Nature547, 468471. doi: 10.1038/nature23272

  • 72

    IannottiF. A.HillC. L.LeoA.AlhusainiA.SoubraneC.MazzarellaE.et al. (2014). Nonpsychotropic plant cannabinoids, Cannabidivarin (CBDV) and Cannabidiol (CBD), activate and desensitize transient receptor potential Vanilloid 1 (TRPV1) channels in vitro: potential for the treatment of neuronal Hyperexcitability. ACS Chem. Neurosci.5, 11311141. doi: 10.1021/cn5000524

  • 73

    Ibarra-LecueI.Unzueta-LarrinagaP.Barrena-BarbadilloR.VillateA.HorrilloI.MendivilB.et al. (2022). Cannabis use selectively modulates circulating biomarkers in the blood of schizophrenia patients. Addict. Biol.27:e13233. doi: 10.1111/adb.13233

  • 74

    JacobsonM. R.WattsJ. J.Da SilvaT.TyndaleR. F.RusjanP. M.HouleS.et al. (2021). Fatty acid amide hydrolase is lower in young cannabis users. Addict. Biol.26:e12872. doi: 10.1111/adb.12872

  • 75

    JakowieckiJ.FilipekS. (2016). Hydrophobic ligand entry and exit pathways of the CB1 cannabinoid receptor. J. Chem. Inf. Model56, 24572466. doi: 10.1021/acs.jcim.6b00499

  • 76

    JékelyG.PapsJ.NielsenC. (2015). The phylogenetic position of ctenophores and the origin(s) of nervous systems. EvoDevo6:1. doi: 10.1186/2041-9139-6-1

  • 77

    JenkoK. J.HirvonenJ.HenterI. D.AndersonK. B.ZoghbiS. S.HydeT. M.et al. (2012). Binding of a tritiated inverse agonist to cannabinoid CB1 receptors is increased in patients with schizophrenia. Schizophr. Res.141, 185188. doi: 10.1016/j.schres.2012.07.021

  • 78

    JiaY.-F.LiY.-C.TangY.-P.CaoJ.WangL.-P.YangY.-X.et al. (2015). Interference of TRPV1 function altered the susceptibility of PTZ-induced seizures. Front. Cell Neurosci.9:20. doi: 10.3389/fncel.2015.00020

  • 79

    JonesD. T.TaylorW. R.ThorntonJ. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci.8, 275282. doi: 10.1093/bioinformatics/8.3.275

  • 80

    JouroukhinY.ZhuX.ShevelkinA. V.HasegawaY.AbazyanB.SaitoA.et al. (2019). Adolescent Δ9-tetrahydrocannabinol exposure and astrocyte-specific genetic vulnerability converge on nuclear factor-κB-Cyclooxygenase-2 signaling to impair memory in adulthood. Biol. Psychiatry85, 891903. doi: 10.1016/j.biopsych.2018.07.024

  • 81

    JuhaszG.LazaryJ.ChaseD.PeggE.DowneyD.TothZ. G.et al. (2009). Variations in the cannabinoid receptor 1 gene predispose to migraine. Neurosci. Lett.461, 116120. doi: 10.1016/j.neulet.2009.06.021

  • 82

    JuknatA.PietrM.KozelaE.RimmermanN.LevyR.GaoF.et al. (2013). Microarray and pathway analysis reveal distinct mechanisms underlying cannabinoid-mediated modulation of LPS-induced activation of BV-2 microglial cells. PLoS One8:e61462. doi: 10.1371/journal.pone.0061462

  • 83

    JungK.-M.SepersM.HenstridgeC. M.LassalleO.NeuhoferD.MartinH.et al. (2012). Uncoupling of the endocannabinoid signalling complex in a mouse model of fragile X syndrome. Nat. Commun.3:1080. doi: 10.1038/ncomms2045

  • 84

    KanoM.Ohno-ShosakuT.HashimotodaniY.UchigashimaM.WatanabeM. (2009). Endocannabinoid-mediated control of synaptic transmission. Physiol. Rev.89, 309380. doi: 10.1152/physrev.00019.2008

  • 85

    KarhsonD. S.KrasinskaK. M.DallaireJ. A.LiboveR. A.PhillipsJ. M.ChienA. S.et al. (2018). Plasma anandamide concentrations are lower in children with autism spectrum disorder. Mol. Autism.9:18. doi: 10.1186/s13229-018-0203-y

  • 86

    KelleyL. A.MezulisS.YatesC. M.WassM. N.SternbergM. J. E. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc.10, 845858. doi: 10.1038/nprot.2015.053

  • 87

    KelleyL. A.SternbergM. J. E. (2009). Protein structure prediction on the web: a case study using the Phyre server. Nat. Protoc.4, 363371. doi: 10.1038/nprot.2009.2

  • 88

    KibretB. G.IshiguroH.HoriuchiY.OnaiviE. S. (2022). New insights and potential therapeutic targeting of CB2 cannabinoid receptors in CNS disorders. Int. J. Mol. Sci.23:975. doi: 10.3390/ijms23020975

  • 89

    KongW.WangX.YangX.HuangW.HanS.YinJ.et al. (2019). Activation of TRPV1 contributes to recurrent febrile seizures via inhibiting the microglial M2 phenotype in the immature brain. Front. Cell. Neurosci.13:442. doi: 10.3389/fncel.2019.00442

  • 90

    Krishna KumarK.RobertsonM. J.ThadhaniE.WangH.SuomivuoriC.-M.PowersA. S.et al. (2023). Structural basis for activation of CB1 by an endocannabinoid analog. Nat. Commun.14:2672. doi: 10.1038/s41467-023-37864-4

  • 91

    KroghA.LarssonB.von HeijneG.SonnhammerE. L. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J. Mol. Biol.305, 567580. doi: 10.1006/jmbi.2000.4315

  • 92

    KuceraR.BouskilaJ.ElkriefL.Fink-JensenA.PalmourR.BouchardJ.-F.et al. (2018). Expression and localization of CB1R, NAPE-PLD, and FAAH in the vervet monkey nucleus accumbens. Sci. Rep.8:8689. doi: 10.1038/s41598-018-26826-2

  • 93

    KumarS.StecherG.LiM.KnyazC.TamuraK. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol.35, 15471549. doi: 10.1093/molbev/msy096

  • 94

    Kwan CheungK. A.PeirisH.WallaceG.HollandO. J.MitchellM. D. (2019). The interplay between the endocannabinoid system, epilepsy and cannabinoids. Int. J. Mol. Sci.20:6079. doi: 10.3390/ijms20236079

  • 95

    KwonD. H.ZhangF.FedorJ. G.SuoY.LeeS.-Y. (2022). Vanilloid-dependent TRPV1 opening trajectory from cryoEM ensemble analysis. Nat. Commun.13:2874. doi: 10.1038/s41467-022-30602-2

  • 96

    LamP. M. W.McDonaldJ.LambertD. G. (2005). Characterization and comparison of recombinant human and rat TRPV1 receptors: effects of exo- and endocannabinoids†. Br. J. Anaesth.94, 649656. doi: 10.1093/bja/aei098

  • 97

    LanducciE.Pellegrini-GiampietroD. E.GianoncelliA.RibaudoG. (2022). Cannabidiol preferentially binds TRPV2: a novel mechanism of action. Neural Regen. Res.17, 26932694. doi: 10.4103/1673-5374.335821

  • 98

    LaprairieR. B.BagherA. M.KellyM. E. M.Denovan-WrightE. M. (2015). Cannabidiol is a negative allosteric modulator of the cannabinoid CB1 receptor. Br. J. Pharmacol.172, 47904805. doi: 10.1111/bph.13250

  • 99

    LaucknerJ. E.JensenJ. B.ChenH.-Y.LuH.-C.HilleB.MackieK. (2008). GPR55 is a cannabinoid receptor that increases intracellular calcium and inhibits M current. Proc. Natl. Acad. Sci. USA105, 26992704. doi: 10.1073/pnas.0711278105

  • 100

    Lazarini-LopesW.Silva-CardosoG. K.Leite-PanissiC. R. A.Garcia-CairascoN. (2022). Increased TRPV1 channels and FosB protein expression are associated with chronic epileptic seizures and Anxiogenic-like behaviors in a preclinical model of temporal lobe epilepsy. Biomedicines10:416. doi: 10.3390/biomedicines10020416

  • 101

    LeeT.-H.LeeJ.-G.YonJ.-M.OhK.-W.BaekI.-J.NahmS.-S.et al. (2011). Capsaicin prevents kainic acid-induced epileptogenesis in mice. Neurochem. Int.58, 634640. doi: 10.1016/j.neuint.2011.01.027

  • 102

    LewekeF. M.GiuffridaA.KoetheD.SchreiberD.NoldenB. M.KranasterL.et al. (2007). Anandamide levels in cerebrospinal fluid of first-episode schizophrenic patients: impact of cannabis use. Schizophr. Res.94, 2936. doi: 10.1016/j.schres.2007.04.025

  • 103

    LeysS. P. (2015). Elements of a “nervous system” in sponges. J. Exp. Biol.218, 581591. doi: 10.1242/jeb.110817

  • 104

    LiuS.PaknejadN.ZhuL.KiharaY.RayM.ChunJ.et al. (2022). Differential activation mechanisms of lipid GPCRs by lysophosphatidic acid and sphingosine 1-phosphate. Nat. Commun.13:731. doi: 10.1038/s41467-022-28417-2

  • 105

    LuH.-C.MackieK. (2016). An introduction to the endogenous cannabinoid system. Biol. Psychiatry79, 516525. doi: 10.1016/j.biopsych.2015.07.028

  • 106

    LuA. T.OgdieM. N.JärvelinM.-R.MoilanenI. K.LooS. K.McCrackenJ. T.et al. (2008). Association of the cannabinoid receptor gene (CNR1) with ADHD and post-traumatic stress disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet.147B, 14881494. doi: 10.1002/ajmg.b.30693

  • 107

    LudányiA.ErossL.CzirjákS.VajdaJ.HalászP.WatanabeM.et al. (2008). Downregulation of the CB1 cannabinoid receptor and related molecular elements of the endocannabinoid system in epileptic human hippocampus. J. Neurosci.28, 29762990. doi: 10.1523/JNEUROSCI.4465-07.2008

  • 108

    LupicaC. R.RiegelA. C.HoffmanA. F. (2004). Marijuana and cannabinoid regulation of brain reward circuits. Br. J. Pharmacol.143, 227234. doi: 10.1038/sj.bjp.0705931

  • 109

    MartinM. W.GrazhdankinD. V.BowringS. A.EvansD. A.FedonkinM. A.KirschvinkJ. L. (2000). Age of Neoproterozoic bilatarian body and trace fossils, White Sea, Russia: implications for metazoan evolution. Science288, 841845. doi: 10.1126/science.288.5467.841

  • 110

    Martín-DuránJ. M.PangK.BørveA.H. S.FuruA.CannonJ. T.et al. (2018). Convergent evolution of bilaterian nerve cords. Nature553, 4550. doi: 10.1038/nature25030

  • 111

    Martínez-AguirreC.CinarR.RochaL. (2022). Targeting endocannabinoid system in epilepsy: for good or for bad. Neuroscience482, 172185. doi: 10.1016/j.neuroscience.2021.12.013

  • 112

    MatsudaL. A.LolaitS. J.BrownsteinM. J.YoungA. C.BonnerT. I. (1990). Structure of a cannabinoid receptor and functional expression of the cloned cDNA. Nature346, 561564. doi: 10.1038/346561a0

  • 113

    MatsuyamaM.TeradaY.Yamazaki-ItoT.ItoK. (2021). A luminescence-based human TRPV1 assay system for quantifying pungency in spicy foods. Food Secur.10:151. doi: 10.3390/foods10010151

  • 114

    MátyásF.UrbánG. M.WatanabeM.MackieK.ZimmerA.FreundT. F.et al. (2008). Identification of the sites of 2-arachidonoylglycerol synthesis and action imply retrograde endocannabinoid signaling at both GABAergic and glutamatergic synapses in the ventral tegmental area. Neuropharmacology54, 95107. doi: 10.1016/j.neuropharm.2007.05.028

  • 115

    McHughD.HuS. S. J.RimmermanN.JuknatA.VogelZ.WalkerJ. M.et al. (2010). N-arachidonoyl glycine, an abundant endogenous lipid, potently drives directed cellular migration through GPR18, the putative abnormal cannabidiol receptor. BMC Neurosci.11:44. doi: 10.1186/1471-2202-11-44

  • 116

    McPartlandJ. M. (2004). Phylogenomic and chemotaxonomic analysis of the endocannabinoid system. Brain Res. Brain Res. Rev.45, 1829. doi: 10.1016/j.brainresrev.2003.11.005

  • 117

    McPartlandJ. M.GlassM.PertweeR. G. (2007). Meta-analysis of cannabinoid ligand binding affinity and receptor distribution: interspecies differences. Br. J. Pharmacol.152, 583593. doi: 10.1038/sj.bjp.0707399

  • 118

    McPartlandJ. M.GuyG. (2004). “The evolution of Cannabis and coevolution with the cannabinoid receptor - a hypothesis” in The medicinal use of Cannabis. eds. GuyG. W.RobsonR.StrongK.WhittleB. (London: Royal Society of Pharmacists), 71102.

  • 119

    McPartlandJ. M.MatiasI.Di MarzoV.GlassM. (2006). Evolutionary origins of the endocannabinoid system. Gene370, 6474. doi: 10.1016/j.gene.2005.11.004

  • 120

    MedvedevaY. V.KimM.-S.UsachevY. M. (2008). Mechanisms of prolonged presynaptic Ca2+ signaling and glutamate release induced by TRPV1 activation in rat sensory neurons. J. Neurosci.28, 52955311. doi: 10.1523/JNEUROSCI.4810-07.2008

  • 121

    MelanciaF.SchiaviS.ServadioM.CartocciV.CampolongoP.PalmeryM.et al. (2018). Sex-specific autistic endophenotypes induced by prenatal exposure to valproic acid involve anandamide signalling. Br. J. Pharmacol.175, 36993712. doi: 10.1111/bph.14435

  • 122

    Méndez-DíazM.Amancio-BelmontO.Estrada-GonzálezV.Ruiz-ContrerasA. E.Prospéro-GarcíaO. (2019). CB1R mediates oleamide’s reward while 5HT2cR mediates aversion in the nucleus accumbens shell of rats. Neurosci. Lett.706, 189193. doi: 10.1016/j.neulet.2019.05.033

  • 123

    MinichinoA.SeniorM.BrondinoN.ZhangS. H.GodlewskaB. R.BurnetP. W. J.et al. (2019). Measuring disturbance of the endocannabinoid system in psychosis: a systematic review and Meta-analysis. JAMA Psychiatry76, 914923. doi: 10.1001/jamapsychiatry.2019.0970

  • 124

    MohapatraD. P.NauC. (2003). Desensitization of capsaicin-activated currents in the Vanilloid receptor TRPV1 is decreased by the cyclic AMP-dependent protein kinase pathway*. J. Biol. Chem.278, 5008050090. doi: 10.1074/jbc.M306619200

  • 125

    MoralesP.MullerC.JagerovicN.ReggioP. H. (2022). Targeting CB2 and TRPV1: computational approaches for the identification of dual modulators. Front. Mol. Biosci.9:841190. doi: 10.3389/fmolb.2022.841190

  • 126

    MoranM. M.XuH.ClaphamD. E. (2004). TRP ion channels in the nervous system. Curr. Opin. Neurobiol.14, 362369. doi: 10.1016/j.conb.2004.05.003

  • 127

    MorozL. L. (2015). Convergent evolution of neural systems in ctenophores. J. Exp. Biol.218, 598611. doi: 10.1242/jeb.110692

  • 128

    MuguruzaC.MorentinB.MeanaJ. J.AlexanderS. P.CalladoL. F. (2019). Endocannabinoid system imbalance in the postmortem prefrontal cortex of subjects with schizophrenia. J. Psychopharmacol.33, 11321140. doi: 10.1177/0269881119857205

  • 129

    MullerC.LynchD. L.HurstD. P.ReggioP. H. (2020). A closer look at anandamide interaction with TRPV1. Front. Mol. Biosci.7:144. doi: 10.3389/fmolb.2020.00144

  • 130

    MullerC.MoralesP.ReggioP. H. (2019). Cannabinoid ligands targeting TRP channels. Front. Mol. Neurosci.11:487. doi: 10.3389/fnmol.2018.00487

  • 131

    MunroS.ThomasK. L.Abu-ShaarM. (1993). Molecular characterization of a peripheral receptor for cannabinoids. Nature365, 6165. doi: 10.1038/365061a0

  • 132

    MusserJ. M.SchippersK. J.NickelM.MizzonG.KohnA. B.PapeC.et al. (2021). Profiling cellular diversity in sponges informs animal cell type and nervous system evolution. Science374, 717723. doi: 10.1126/science.abj2949

  • 133

    NabissiM.MorelliM. B.SantoniM.SantoniG. (2013). Triggering of the TRPV2 channel by cannabidiol sensitizes glioblastoma cells to cytotoxic chemotherapeutic agents. Carcinogenesis34, 4857. doi: 10.1093/carcin/bgs328

  • 134

    NadezhdinK. D.NeubergerA.NikolaevY. A.MurphyL. A.GrachevaE. O.BagriantsevS. N.et al. (2021). Extracellular cap domain is an essential component of the TRPV1 gating mechanism. Nat. Commun.12:2154. doi: 10.1038/s41467-021-22507-3

  • 135

    NavarroD.GasparyanA.NavarreteF.TorregrosaA. B.RubioG.Marín-MayorM.et al. (2022). Molecular alterations of the endocannabinoid system in psychiatric disorders. Int. J. Mol. Sci.23:4764. doi: 10.3390/ijms23094764

  • 136

    NeelandsT. R.JarvisM. F.HanP.FaltynekC. R.SurowyC. S. (2005). Acidification of rat TRPV1 alters the kinetics of capsaicin responses. Mol. Pain1:1744-8069-1-28. doi: 10.1186/1744-8069-1-28

  • 137

    NeubergerA.OdaM.NikolaevY. A.NadezhdinK. D.GrachevaE. O.BagriantsevS. N.et al. (2023). Human TRPV1 structure and inhibition by the analgesic SB-366791. Nat. Commun.14:2451. doi: 10.1038/s41467-023-38162-9

  • 138

    NewellK. A.DengC.HuangX.-F. (2006). Increased cannabinoid receptor density in the posterior cingulate cortex in schizophrenia. Exp. Brain Res.172, 556560. doi: 10.1007/s00221-006-0503-x

  • 139

    O’SullivanS. E. (2016). An update on PPAR activation by cannabinoids. Br. J. Pharmacol.173, 18991910. doi: 10.1111/bph.13497

  • 140

    Ortega-HernándezJ.Lerosey-AubrilR.LossoS. R.WeaverJ. C. (2022). Neuroanatomy in a middle Cambrian mollisoniid and the ancestral nervous system organization of chelicerates. Nat. Commun.13:410. doi: 10.1038/s41467-022-28054-9

  • 141

    OstlundI.Von GuntenM.SmithC.EdwardsJ. G. (2022). Chronic Δ9-tetrahydrocannabinol impact on plasticity, and differential activation requirement for CB1-dependent long-term depression in ventral tegmental area GABA neurons in adult versus young mice. Front. Neurosci.16:1067493. doi: 10.3389/fnins.2022.1067493

  • 142

    OvertonH. A.BabbsA. J.DoelS. M.FyfeM. C. T.GardnerL. S.GriffinG.et al. (2006). Deorphanization of a G protein-coupled receptor for oleoylethanolamide and its use in the discovery of small-molecule hypophagic agents. Cell Metab.3, 167175. doi: 10.1016/j.cmet.2006.02.004

  • 143

    ParryL.CaronJ.-B. (2019). Canadia spinosa and the early evolution of the annelid nervous system. Sci. Adv.5:eaax5858. doi: 10.1126/sciadv.aax5858

  • 144

    PasierskiM.SzulczykB. (2020). Capsaicin inhibits sodium currents and epileptiform activity in prefrontal cortex pyramidal neurons. Neurochem. Int.135:104709. doi: 10.1016/j.neuint.2020.104709

  • 145

    PatesS.BottingJ. P.MuirL. A.WolfeJ. M. (2022). Ordovician opabiniid-like animals and the role of the proboscis in euarthropod head evolution. Nat. Commun.13:6969. doi: 10.1038/s41467-022-34204-w

  • 146

    PecoitsE.KonhauserK. O.AubetN. R.HeamanL. M.VeroslavskyG.SternR. A.et al. (2012). Bilaterian burrows and grazing behavior at >585 million years ago. Science336, 16931696. doi: 10.1126/science.1216295

  • 147

    PezzoliM.ElhamdaniA.CamachoS.MeystreJ.GonzálezS. M.le CoutreJ.et al. (2014). Dampened neural activity and abolition of epileptic-like activity in cortical slices by active ingredients of spices. Sci. Rep.4:6825. doi: 10.1038/srep06825

  • 148

    ProudnikovD.KroslakT.SipeJ. C.RandesiM.LiD.HamonS.et al. (2010). Association of polymorphisms of the cannabinoid receptor (CNR1) and fatty acid amide hydrolase (FAAH) genes with heroin addiction: impact of long repeats of CNR1. Pharmacogenomics J.10, 232242. doi: 10.1038/tpj.2009.59

  • 149

    QinN.NeeperM. P.LiuY.HutchinsonT. L.LubinM. L.FloresC. M. (2008). TRPV2 is activated by Cannabidiol and mediates CGRP release in cultured rat dorsal root ganglion neurons. J. Neurosci.28, 62316238. doi: 10.1523/JNEUROSCI.0504-08.2008

  • 150

    R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. Available at: http://www.R-project.org/

  • 151

    RamakerR. C.BowlingK. M.LasseigneB. N.HagenauerM. H.HardiganA. A.DavisN. S.et al. (2017). Post-mortem molecular profiling of three psychiatric disorders. Genome Med.9:72. doi: 10.1186/s13073-017-0458-5

  • 152

    RanganathanM.Cortes-BrionesJ.RadhakrishnanR.ThurnauerH.PlanetaB.SkosnikP.et al. (2016). Reduced brain cannabinoid receptor availability in schizophrenia. Biol. Psychiatry79, 9971005. doi: 10.1016/j.biopsych.2015.08.021

  • 153

    RodríguezG.NeugebauerN. M.YaoK. L.MeltzerH. Y.CsernanskyJ. G.DongH. (2017). Δ9-tetrahydrocannabinol (Δ9-THC) administration after neonatal exposure to phencyclidine potentiates schizophrenia-related behavioral phenotypes in mice. Pharmacol. Biochem. Behav.159, 611. doi: 10.1016/j.pbb.2017.06.010

  • 154

    RStudio Team (2020). RStudio: Integrated development for R. Available at: http://www.rstudio.com/

  • 155

    RussoP.StrazzulloP.CappuccioF. P.TregouetD. A.LauriaF.LoguercioM.et al. (2007). Genetic variations at the endocannabinoid type 1 receptor gene (CNR1) are associated with obesity phenotypes in men. J. Clin. Endocrinol. Metab.92, 23822386. doi: 10.1210/jc.2006-2523

  • 156

    RybergE.LarssonN.SjögrenS.HjorthS.HermanssonN.-O.LeonovaJ.et al. (2007). The orphan receptor GPR55 is a novel cannabinoid receptor. Br. J. Pharmacol.152, 10921101. doi: 10.1038/sj.bjp.0707460

  • 157

    SaffarzadehF.EslamizadeM. J.GhadiriT.Modarres MousaviS. M.HadjighassemM.GorjiA. (2015). Effects of TRPV1 on the hippocampal synaptic plasticity in the epileptic rat brain. Synapse69, 375383. doi: 10.1002/syn.21825

  • 158

    SaitoS.FukutaN.ShingaiR.TominagaM. (2011). Evolution of vertebrate transient receptor potential vanilloid 3 channels: opposite temperature sensitivity between mammals and western clawed frogs. PLoS Genet.7:e1002041. doi: 10.1371/journal.pgen.1002041

  • 159

    SaitoS.NakatsukaK.TakahashiK.FukutaN.ImagawaT.OhtaT.et al. (2012). Analysis of transient receptor potential ankyrin 1 (TRPA1) in frogs and lizards illuminates both nociceptive heat and chemical sensitivities and coexpression with TRP vanilloid 1 (TRPV1) in ancestral vertebrates. J. Biol. Chem.287, 3074330754. doi: 10.1074/jbc.M112.362194

  • 160

    SaitoS.SaitoC. T.IgawaT.TakedaN.KomakiS.OhtaT.et al. (2022). Evolutionary tuning of transient receptor potential Ankyrin 1 underlies the variation in heat avoidance behaviors among frog species inhabiting diverse thermal niches. Mol. Biol. Evol.39:msac180. doi: 10.1093/molbev/msac180

  • 161

    SanchisP.LavignolleR.AbbateM.Lage-VickersS.VazquezE.CotignolaJ.et al. (2021). Analysis workflow of publicly available RNA-sequencing datasets. STAR Protoc.2:100478. doi: 10.1016/j.xpro.2021.100478

  • 162

    SbarskiB.AkiravI. (2018). Chronic exposure to cannabinoids before an emotional trauma may have negative effects on emotional function. Eur. Neuropsychopharmacol.28, 955969. doi: 10.1016/j.euroneuro.2018.05.008

  • 163

    SchrödingerL.DeLanoW. (2020). PyMol. Available at: http://www.pymol.org/pymol

  • 164

    SeillierA.MartinezA. A.GiuffridaA. (2020). Differential effects of Δ9-tetrahydrocannabinol dosing on correlates of schizophrenia in the sub-chronic PCP rat model. PLoS One15:e0230238. doi: 10.1371/journal.pone.0230238

  • 165

    ServadioM.MelanciaF.ManducaA.di MasiA.SchiaviS.CartocciV.et al. (2016). Targeting anandamide metabolism rescues core and associated autistic-like symptoms in rats prenatally exposed to valproic acid. Transl. Psychiatry6:e902. doi: 10.1038/tp.2016.182

  • 166

    ShaoZ.YinJ.ChapmanK.GrzemskaM.ClarkL.WangJ.et al. (2016). High-resolution crystal structure of the human CB1 cannabinoid receptor. Nature540, 602606. doi: 10.1038/nature20613

  • 167

    SiniscalcoD.SaponeA.GiordanoC.CirilloA.de MagistrisL.RossiF.et al. (2013). Cannabinoid receptor type 2, but not type 1, is up-regulated in peripheral blood mononuclear cells of children affected by autistic disorders. J. Autism Dev. Disord.43, 26862695. doi: 10.1007/s10803-013-1824-9

  • 168

    SrivastavaP. K.van EyllJ.GodardP.MazzuferiM.Delahaye-DuriezA.Van SteenwinckelJ.et al. (2018). A systems-level framework for drug discovery identifies Csf1R as an anti-epileptic drug target. Nat. Commun.9:3561. doi: 10.1038/s41467-018-06008-4

  • 169

    StarkT.IannottiF. A.Di MartinoS.Di BartolomeoM.Ruda-KucerovaJ.PiscitelliF.et al. (2022). Early blockade of CB1 receptors ameliorates schizophrenia-like alterations in the neurodevelopmental MAM model of schizophrenia. Biomol. Ther.12:108. doi: 10.3390/biom12010108

  • 170

    StarkusJ.JansenC.ShimodaL. M. N.StokesA. J.Small-HowardA. L.TurnerH. (2019). Diverse TRPV1 responses to cannabinoids. Channels (Austin)13, 172191. doi: 10.1080/19336950.2019.1619436

  • 171

    StefanoG. B.SalzetB.SalzetM. (1997). Identification and characterization of the leech CNS cannabinoid receptor: coupling to nitric oxide release. Brain Res.753, 219224. doi: 10.1016/s0006-8993(96)01484-9

  • 172

    StelzerG.RosenN.PlaschkesI.ZimmermanS.TwikM.FishilevichS.et al. (2016). The GeneCards suite: from gene data mining to disease genome sequence analyzes. Curr. Protoc. Bioinformatics54, 1.30.11.30.33. doi: 10.1002/cpbi.5

  • 173

    StraikerA.DvorakovaM.Bosquez-BergerT.BlahosJ.MackieK. (2023). A collection of cannabinoid-related negative findings from autaptic hippocampal neurons. Sci. Rep.13:9610. doi: 10.1038/s41598-023-36710-3

  • 174

    SuT.YanY.LiQ.YeJ.PeiL. (2021). Endocannabinoid system unlocks the puzzle of autism treatment via microglia. Front. Psych.12:734837. doi: 10.3389/fpsyt.2021.734837

  • 175

    SugayaY.KanoM. (2022). Endocannabinoid-mediated control of neural circuit excitability and epileptic seizures. Front. Neural Circuits15:781113. doi: 10.3389/fncir.2021.781113

  • 176

    SylantyevS.JensenT. P.RossR. A.RusakovD. A. (2013). Cannabinoid- and lysophosphatidylinositol-sensitive receptor GPR55 boosts neurotransmitter release at central synapses. Proc. Natl. Acad. Sci. USA110, 51935198. doi: 10.1073/pnas.1211204110

  • 177

    SzklarczykD.GableA. L.LyonD.JungeA.WyderS.Huerta-CepasJ.et al. (2019). STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res.47, D607D613. doi: 10.1093/nar/gky1131

  • 178

    SzűcsE.DvorácskóS.TömbölyC.BükiA.KékesiG.HorváthG.et al. (2016). Decreased CB receptor binding and cannabinoid signaling in three brain regions of a rat model of schizophrenia. Neurosci. Lett.633, 8793. doi: 10.1016/j.neulet.2016.09.020

  • 179

    TanakaM.SackettS.ZhangY. (2020). Endocannabinoid modulation of microglial phenotypes in neuropathology. Front. Neurol.11:87. doi: 10.3389/fneur.2020.00087

  • 180

    ThomasA.BaillieG. L.PhillipsA. M.RazdanR. K.RossR. A.PertweeR. G. (2007). Cannabidiol displays unexpectedly high potency as an antagonist of CB1 and CB2 receptor agonists in vitro. Br. J. Pharmacol.150, 613623. doi: 10.1038/sj.bjp.0707133

  • 181

    TominagaM.CaterinaM. J.MalmbergA. B.RosenT. A.GilbertH.SkinnerK.et al. (1998). The cloned capsaicin receptor integrates multiple pain-producing stimuli. Neuron21, 531543. doi: 10.1016/s0896-6273(00)80564-4

  • 182

    TrottO.OlsonA. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J. Comput. Chem.31, 455461. doi: 10.1002/jcc.21334

  • 183

    UrigüenL.García-FusterM. J.CalladoL. F.MorentinB.La HarpeR.CasadóV.et al. (2009). Immunodensity and mRNA expression of A2A adenosine, D2 dopamine, and CB1 cannabinoid receptors in postmortem frontal cortex of subjects with schizophrenia: effect of antipsychotic treatment. Psychopharmacology206, 313324. doi: 10.1007/s00213-009-1608-2

  • 184

    VasquezV.MugarraR. C.BellB.LeeJ.Cordero-MoralesJ. F. (2020). Inhibition of inositol Monophosphatase enhances TRPV1 function in vivo. Biophys. J.118:410a. doi: 10.1016/j.bpj.2019.11.2321

  • 185

    VetterI.WyseB. D.MonteithG. R.Roberts-ThomsonS. J.CabotP. J. (2006). The μ opioid agonist morphine modulates potentiation of capsaicin-evoked TRPV1 responses through a cyclic AMP-dependent protein kinase a pathway. Mol. Pain2:22. doi: 10.1186/1744-8069-2-22

  • 186

    VoetsT.DroogmansG.WissenbachU.JanssensA.FlockerziV.NiliusB. (2004). The principle of temperature-dependent gating in cold- and heat-sensitive TRP channels. Nature430, 748754. doi: 10.1038/nature02732

  • 187

    WangX.YangX.-L.KongW.-L.ZengM.-L.ShaoL.JiangG.-T.et al. (2019). TRPV1 translocated to astrocytic membrane to promote migration and inflammatory infiltration thus promotes epilepsy after hypoxic ischemia in immature brain. J. Neuroinflammation16:214. doi: 10.1186/s12974-019-1618-x

  • 188

    WatanabeH.VriensJ.PrenenJ.DroogmansG.VoetsT.NiliusB. (2003). Anandamide and arachidonic acid use epoxyeicosatrienoic acids to activate TRPV4 channels. Nature424, 434438. doi: 10.1038/nature01807

  • 189

    WeiD.DinhD.LeeD.LiD.AngurenA.Moreno-SanzG.et al. (2016). Enhancement of anandamide-mediated endocannabinoid signaling corrects autism-related social impairment. Cannabis Cannabinoid Res.1, 8189. doi: 10.1089/can.2015.0008

  • 190

    WilsonR. I.NicollR. A. (2001). Endogenous cannabinoids mediate retrograde signalling at hippocampal synapses. Nature410, 588592. doi: 10.1038/35069076

  • 191

    WongD. F.KuwabaraH.HortiA. G.RaymontV.BrasicJ.GuevaraM.et al. (2010). Quantification of cerebral cannabinoid receptors subtype 1 (CB1) in healthy subjects and schizophrenia by the novel PET radioligand [11C]OMAR. NeuroImage52, 15051513. doi: 10.1016/j.neuroimage.2010.04.034

  • 192

    WuY.-W.KimJ.-I.TawfikV. L.LalchandaniR. R.ScherrerG.DingJ. B. (2015). Input- and cell type-specific endocannabinoid-dependent LTD in the striatum. Cell Rep.10, 7587. doi: 10.1016/j.celrep.2014.12.005

  • 193

    WuH.-F.LuT.-Y.ChuM.-C.ChenP. S.LeeC.-W.LinH.-C. (2020). Targeting the inhibition of fatty acid amide hydrolase ameliorate the endocannabinoid-mediated synaptic dysfunction in a valproic acid-induced rat model of autism. Neuropharmacology162:107736. doi: 10.1016/j.neuropharm.2019.107736

  • 194

    Yañez-GuerraL. A.ThielD.JékelyG. (2022). Premetazoan origin of neuropeptide signaling. Mol. Biol. Evol.39:msac051. doi: 10.1093/molbev/msac051

  • 195

    YangJ.Ortega-HernándezJ.ButterfieldN. J.LiuY.BoyanG. S.HouJ.et al. (2016). Fuxianhuiid ventral nerve cord and early nervous system evolution in Panarthropoda. Proc. Natl. Acad. Sci.113, 29882993. doi: 10.1073/pnas.1522434113

  • 196

    YangJ.Ortega-HernándezJ.ButterfieldN. J.ZhangX. (2013). Specialized appendages in fuxianhuiids and the head organization of early euarthropods. Nature494, 468471. doi: 10.1038/nature11874

  • 197

    ZamberlettiE.GabaglioM.Woolley-RobertsM.BinghamS.RubinoT.ParolaroD. (2019). Cannabidivarin treatment ameliorates autism-like behaviors and restores hippocampal endocannabinoid system and glia alterations induced by prenatal Valproic acid exposure in rats. Front. Cell. Neurosci.13:367. doi: 10.3389/fncel.2019.00367

  • 198

    ZavitsanouK.GarrickT.HuangX. F. (2004). Selective antagonist [3H]SR141716A binding to cannabinoid CB1 receptors is increased in the anterior cingulate cortex in schizophrenia. Prog. Neuro-Psychopharmacol. Biol. Psychiatry28, 355360. doi: 10.1016/j.pnpbp.2003.11.005

  • 199

    ZhengW.WenH. (2019). Heat activation mechanism of TRPV1: new insights from molecular dynamics simulation. Temperature (Austin)6, 120131. doi: 10.1080/23328940.2019.1578634

  • 200

    ZuoY.IemoloA.Montilla-PerezP.LiH.-R.YangX.TeleseF. (2022). Chronic adolescent exposure to cannabis in mice leads to sex-biased changes in gene expression networks across brain regions. Neuropsychopharmacology47, 20712080. doi: 10.1038/s41386-022-01413-2

Summary

Keywords

CB1R, TRPV1, endocannabinoid system, evolution, protein structure, phylogenetic tree, transcriptomics, bioinformatics

Citation

Cherry AL, Wheeler MJ, Mathisova K and Di Miceli M (2024) In silico analyses of the involvement of GPR55, CB1R and TRPV1: response to THC, contribution to temporal lobe epilepsy, structural modeling and updated evolution. Front. Neuroinform. 18:1294939. doi: 10.3389/fninf.2024.1294939

Received

28 September 2023

Accepted

19 January 2024

Published

07 February 2024

Volume

18 - 2024

Edited by

Gabor Juhasz, Eötvös Loránd University, Hungary

Reviewed by

José Luis Castañeda-Cabral, University of Guadalajara, Mexico

Balapal Basavarajappa, New York University, United States

Rui Milton Patricio Silva Junior, University of São Paulo, Brazil

Updates

Copyright

*Correspondence: Mathieu Di Miceli,

†These authors have contributed equally to this work and share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics