key: cord-1023430-jn1vc1na authors: Mattoteia, Daiana; Schiano Moriello, Aniello; Taglialatela-Scafati, Orazio; Amodeo, Pietro; De Petrocellis, Luciano; Appendino, Giovanni; Vitale, Rosa Maria; Caprioglio, Diego title: The Combined Effect of Branching and Elongation on the Bioactivity Profile of Phytocannabinoids. Part I: Thermo-TRPs date: 2021-08-23 journal: Biomedicines DOI: 10.3390/biomedicines9081070 sha: 5be763216ee81580b3a9f663b28af7fc845e416f doc_id: 1023430 cord_uid: jn1vc1na The affinity of cannabinoids for their CB(1) and CB(2) metabotropic receptors is dramatically affected by a combination of α-branching and elongation of their alkyl substituent, a maneuver exemplified by the n-pentyl -> α,α-dimethylheptyl (DMH) swap. The effect of this change on other cannabinoid end-points is still unknown, an observation surprising since thermo-TRPs are targeted by phytocannabinoids with often sub-micromolar affinity. To fill this gap, the α,α-dimethylheptyl analogues of the five major phytocannabinoids [CBD (1a), Δ(8)-THC (6a), CBG (7a), CBC (8a) and CBN (9a)] were prepared by total synthesis, and their activity on thermo-TRPs (TRPV1-4, TRPM8, and TRPA1) was compared with that of one of their natural analogues. Surprisingly, the DMH chain promoted a shift in the selectivity toward TRPA1, a target involved in pain and inflammatory diseases, in all investigated compounds. A comparative study of the putative binding modes at TRPA1 between DMH-CBC (8b), the most active compound within the series, and CBC (8a) was carried out by molecular docking, allowing the rationalization of their activity in terms of structure–activity relationships. Taken together, these observations qualify DMH-CBC (8b) as a non-covalent TRPA1-selective cannabinoid lead that is worthy of additional investigation as an analgesic and anti-inflammatory agent. Phytocannabinoids from cannabis (Cannabis sativa L.) are characterized by a linear alkyl substituent spanning one to six carbons bound to their resorcinol core, with the pentyl substitution being by far the most common one [1] . The existence of branched chain cannabinoids of the isoand ante-iso series is, in principle, plausible from a biogenetic standpoint, since the alkylresorcinol moiety of phytocannabinoids is of ketide derivation, and branching would therefore simply require the replacement of the acetate-derived starter with one derived from a branched amino acid. However, branched phytocannabinoids of this type have only been tentatively detected as trace constituents of cannabis, largely remaining unconfirmed curiosities in its inventory of constituents [1] . Conversely, nonbiogenetic branching at the benzyl carbon of the alkyl residue has played a critical role in research on cannabis and cannabinoids ever since the early synthesis of these compounds by Adams in the early forties of the past century [2] . Adams was unsuccessful in the isolation of the narcotic principle of cannabis, but the generation of a mixture of narcotic tetrahydrocannabinols from the acidic treatment of cannabidiol (CBD, Figure 1, 1a) made him correctly postulate a structure of this type for this elusive principle, whose structure was eventually established as 2 by Mechoulam two decades later [3] . In the course of these studies, Adams discovered that ∆ 6a,10a -tetrahydrocannabinol (3a), an unnatural compound relatively easily available by total synthesis, could replicate most of the activity of the elusive narcotic principle of cannabis, and investigated the activity of analogues where the n-pentyl substituent of the model compound was replaced by longer, shorter, and branched alkyl groups [2] . During these studies, the remarkable potency of compounds bearing substituents on the benzyl and the homobenzyl positions of the alkyl group was discovered. One of these compounds, named pyrahexyl (3b) [4] , was next developed as an incapacitating non-lethal war weapon and as anti-riot agent, but these studies terminated when the US signed the convention on the ban of chemical weapons in the early sixties [2] . A decade later, capitalizing on these observations, Pfizer chemists developed a series of ultra-potent and less lipophilic analogues of ∆ 9 -tetrahydrocannabinol (2) that proved instrumental in identifying the two cannabinoid receptors, and that are still used today as a reference for cannabinoid activity, as exemplified by CP-55,940 (4). Unfortunately, these compounds are also popular designer drugs for herbal incenses (spices), as shown by cannabicyclohexanol (5), the first synthetic cannabinoid discovered in this type of product [4] . Phytocannabinoids are multi-target agents, capable of modulating a variety of endpoints that includes not only metabotropic receptors (CB 1 , CB 2 ), but also enzymes, nuclear receptors, and ionotropic receptors belonging to the family of thermo-TRPs (TRPV1-V4, TRPM8, TRPA1) [5] . In particular, thermo-TRPs represent relevant pharmacological targets in a wide array of pathological conditions ranging from pain to inflammatory and respiratory diseases [6] , including those linked to COVID-19 [7] . Surprisingly, no systematic attempt has been made so far to investigate the effect of branching and elongation on the overall biological profile of phytocannabinoids. To fill this gap, we have synthesized the α,α-dimethylheptyl analogues of the five major phytocannabinoids (cannabidiol (CBD, 1a), ∆ 8 -tetrahydrocannabinol (∆ 8 -THC, Figure 2 , 6a) ( This biologically similar ∆ 8 -analogue of 3 was used because of its major stability and ease of synthesis.), cannabigerol (CBG, 7a), cannabichromene (CBC, 8a), and cannabinol (CBN, 9a)) and compared their profile of activity against thermo-TRPs with the one of their natural analogues. IR spectra were recorded on an Avatar 370 FT-IR Techno-Nicolet apparatus. 1 H (400 MHz) and 13 C (100 MHz) NMR spectra were measured on a Bruker Avance 400 MHz spectrometer. Chemical shifts were referenced to the residual solvent signal (CDCl 3 : δ H = 7.21, δ C = 77.0). Homonuclear 1 H connectivities were determined by the correlation spectroscopy (COSY) experiment. One-bond heteronuclear 1 H-13 C connectivities were determined with the heteronuclear single quantum coherence (HSQC) spectroscopy experiment. Two-and three-bond 1 H-13 C connectivities were determined by gradient twodimensional (2D) heteronuclear multiple bond correlation (HMBC) experiments optimized for a 2, 3 J = 9 Hz. Low-and high-resolution electrospray ionization mass spectrometry (ESI-MS) data were determined on an LTQ OrbitrapXL (Thermo Scientific) mass spectrometer. Reactions were monitored by thin-layer chromatography (TLC) on Merck 60 F254 (0.25 mm) plates, visualized by staining with 5% H 2 SO 4 in EtOH and heating. Organic phases were dried with Na 2 SO 4 before evaporation. Chemical reagents and solvents were purchased from Sigma-Aldrich, TCI Europe, or Fluorchem, and were used without further purification unless stated otherwise. Petroleum ether with a boiling point of 40-60 • C was used. Silica gel 60 (70-230 mesh) was used for gravity column chromatography (GCC). All work-up solutions were dried with Na 2 SO 4 before evaporation. Under a nitrogen atmosphere, BF 3 OEt 2 (1.5 mL/g substrate, 150 µL) was added to a stirred suspension of Al 2 O 3 (10 g/g of substrate, 1 g) in dry DCM (10 mL). After stirring for 15 min at room temperature, the suspension was heated to 40 • C for 1 min, and then the resorcinol 10 (100 mg, 0.423 mmol) and (1S,4R)-p-mentha-2,8-dien-1-ol (11.52 mg, 0.34 mmol, 0.8 molar equiv.) were sequentially added. The reaction was stirred at 40 • C for 10 s. and then quenched with 5 mL of sat. Na 2 CO 3 . The resulting biphasic system was extracted with EtOAc, and the organic phase was washed with brine, dried, and concentrated under reduced pressure. Purification by flash column chromatography on silica gel (PE 100% to PE-EtOAc 95:5 as eluent) gave 1b (94 mg, 60%, Rf: 0.75 (PE-EtOAc 9:1)) as a colorless oil. 1 To a stirred suspension of 10 (450 mg, 1.69 mmol) in toluene (2.5 mL), p-toluensulfonic acid (PTSA, 60 mg, 0.34 mmol, 0.2 molar equiv.), and (1S,4R)-p-mentha-2,8-dien-1-ol (11,290 mg, 1.86 mmol, 1.1 molar equiv.) were sequentially added. The solution was heated at 120 • C for two hours and then cooled to room temperature and quenched with brine. The biphasic system was extracted with EtOAc, and the organic phase was dried and evaporated. Purification by flash column chromatography on silica gel (PE 100% to PE/DCM 8:2 as eluent) gave 6b (532 mg, 85%, Rf: 0.90 (PE-EtOAc 9:1)) as a brown oil. 1 Under a nitrogen atmosphere, BF 3 Et 2 O (1.5 mL/g substrate, 150 µL) was added to a stirred suspension of Al 2 O 3 (10 g/g of substrate, 1 g) in dry DCM (10 mL). After stirring 15 min at room temperature, the suspension was heated to 40 • C for 1 min, and then the resorcinol 10 (100 mg, 0.42 mmol) and geraniol (12, To a stirred suspension of 10 (350 mg, 1.48 mmol) in toluene (10 mL), n-butylamine (141 µL, 104 mg, 1.42 mmol, 1 molar equiv.) was added. The solution was heated to 60 • C for 10 min, then citral (13, 241 µL, 217 mg, 1.42 mmol, 1 molar equiv.) was added and the solution was refluxed overnight. The reaction was cooled to room temperature and quenched with 20 mL 2M H 2 SO 4 , and then extracted with EtOAc. The organic phase was washed with brine, dried, and concentrated under reduced pressure. Purification by flash column chromatography on silica gel (PE 100% to PE-EtOAc 95:5 as eluent) gave 210 mg 8b (210 mg, 40%, Rf: 0.85 (PE-EtOAc 9:1)) as a brown oil. 1 .2950). The same protocol was used for the preparation of 8c and 8d, starting, from 5-n-heptylresorcinol and 5-α,α-dimethylpentylresorcinol, respectively. Depentyl-5-n-heptylcannabichromene (8c): brown oil (42%, Rf: 0.87 (PE-EtOAc 9:1)). 1 To a stirred suspension of 8b (210 mg, 0.57 mmol) in toluene (40 mL), iodine (291 mg, 1.15 mmol, 2 eq) was added. The solution was refluxed for three hours then cooled to room temperature and quenched with sat. Na 2 SO 3 . The mixture was extracted with EtOAc, dried, and concentrated under reduced pressure. Purification by flash column chromatography on silica gel (PE 100% to PE-CH 2 Cl 2 9:1 as eluent) gave 9b (173 mg, 83%, Rf: 0.85 (PE-EtOAc 9:1)) as a brown oil. 1 Effects of compounds on intracellular Ca 2+ concentration ([Ca 2+ ] i ) were determined using Fluo-4, a selective intracellular fluorescent probe for Ca 2+ . Assays of rat TRPA1, TRPV2, TRPV3, TRPV4, and TRPM8 or human TRPV1 mediated elevation of intracellular Ca 2+ in transfected HEK-293 cells were performed with a continuous monitoring of [Ca 2+ ] i during the experiments [8] . Briefly, human embryonic kidney (HEK-293) cells, stably transfected with rat TRPA1, TRPV2, TRPV3, TRPV4, and TRPM8 or human TRPV1 (selected by geneticin 600 µg mL −1 ) or not transfected, were cultured in EMEM+ 2 mM glutamine + 1% nonessential amino acids +10% FBS and maintained at 37 • C with 5% CO 2 . TRP-HEK-293 cells express stably high levels of TRP transcripts, while these transcripts were virtually absent in wild-type HEK-293 cells as checked by real-time PCR. On the day of the experiment, the cells were loaded for 1 h in the dark at room temperature with 2+ ] i . In the TRPV3 assay, TRPV3-expressing HEK293 cells were first sensitized with the structurally unrelated agonist 2-aminoethoxydiphenyl borate (100 µM). Antagonist/desensitizing behavior is evaluated against the agonist of the TRP analyzed by adding the compounds directly in the quartz cuvette 5 min before stimulation of cells with the agonist. Allylisothiocyanate (100 µM) was used for TRPA1, capsaicin (0.1 µM) for TRPV1, cannabidiol (2 µM) for TRPV2, thymol (100 µM) for TRPV3, GSK1016790A (10 nM) for TRPV4, and icilin (0.25 µM) for TRPM8. The IC 50 value was expressed as the concentration exerting a half-maximal inhibition of agonist effect, taking as 100% the effect on [Ca 2+ ] i exerted by the agonist alone. Dose-response curve fitting (sigmoidal dose-response variable slope) and parameter estimation were performed with Graph-Pad Prism8 (GraphPad Software Inc., San Diego, CA, USA). All determinations were performed at least in triplicate. Docking studies were performed with AutoDock 4.21 [9] . The rTRPA1 model, along with the ligands, were processed with AutoDock Tools (ADT) package version 1.5.6rc1 to merge non polar hydrogens, calculate Gasteiger charges, and select rotatable sidechain bonds. Grid dimensions of 60 × 60 × 50, centered in the putative binding pocket, were generated with the program AutoGrid 4.2 included in Auto-dock 4.2 distribution, with a spacing of 0.375 Å. Docking runs were carried out by either keeping the whole protein fixed or allowing the rotation of selected residues (Leu873, Leu884, Phe912, Met915 and Met949). A total of 100 molecular AutoDock docking runs for each docking calculation were performed adopting a Lamarckian genetic algorithm (LGA) and the following associated parameters: 100 individuals in a population with a maximum of 15 million energy evaluations and a maximum of 37,000 generations, followed by 300 iterations of Solis and Wets local search. Flexibility was used for all rotatable bonds of both docked ligands. The representative poses for each receptor were selected for the subsequent energy minimization with Amber16 package [10] using ff14SB force field for the protein, and gaff parameters for the ligand [11, 12] . All α,α-DMH analogs were prepared from 5-(1,1-dimethylhepty)-resorcinol (10), in turn prepared from 3,5-dimethoxybenzoic acid [13] . Terpenylation was then carried out mutuating chemistry developed for the synthesis of the corresponding natural n-pentyl phytocannabinoids. C-Menthylation with (4R)-2,9-p-menthadien-1-ol (11) afforded either DMH-CBD (1b) or DMH-∆ 8 -THC (6b) depending on the reaction conditions [14] , while geranylation with BF 3 on alumina gave DMH-CBG (7b) [15] , and chromenyation with citral (13) under basic conditions generated DMH-CBC (8b) [16] , next aromatized to DMH-CBN (9b) by treatment with iodine (Scheme 1) [17] . All DMH-derivatives retained the full agonistic profile on TRPA1 typical of their naturally occurring n-pentyl analogues, along with a comparable efficacy and ability to desensitize this channel (Table 1) . Conversely, the replacement of the pentyl chain with the branched DMH group was detrimental for the activity toward TRPV1-4 and TRPM8 channels (Tables 1 and 2) . Thus, while CBD (1a) and CBG (7a) are TRPV1 antagonists in a sub/low micromolar range, their correspondent DMH derivatives (1b and 7b, respectively) were inactive against this receptor. Similarly, while the natural derivatives are antagonists in a low micromolar range at TRPV4, the activity was lost for all DMH-derivatives, that were only weak inhibitors of TRPV2, with IC 50 ranging from 9 to >50 µM. As for TRPV3, DMH-∆ 8 -THC (6b) was the only compound active as an inhibitor at a low micromolar range (IC 50 of~2 µM), with a slightly better potency than ∆ 8 -THC (6a). At TRPM8, al DMH analogues acted as pure antagonists albeit with a lower potency than the natural derivatives, except DMH-CBC (8b), which was totally inactive, and with only DMH-CBN (9b) retaining a certain affinity (IC 50~1 µM). To better rationalize the selectivity of the DMH-derivatives toward TRPA1 and dissect the relative contribution of chain elongation and branching, two point-mutated analogues of the most potent compound (DMH-CBC (8b)) were prepared according to the general synthetic sequence summarized in Scheme 1 starting from a different 5-alkylresorcinol. In the n-heptyl chromene 8c, the alkyl chain is linear and expanded by two carbon atoms compared to CBC (8a), while in the α,α-dimethylpentyl analogue 8d the benzyl branching is implanted on a pentyl chain. As shown in Table 1 , the n-heptyl derivative 8c was less potent than the dimethyl-pentyl analogue 8d (IC 50 = 3.8 and 0.40 µM, respectively). Hence, the α,α-dimethyl substitution is the major contributor to the activity of the DHM derivative (IC 50 = 0.32 µM). To shed light on the putative binding modes of DMH-CBC (8a) on TRPA1 channels, a molecular docking study was carried out using our previous homology model of rat-TRPA1 [20] . Since the compound is a racemate, both enantiomers were considered, using CBC (8a) as a reference compound. Different docking runs were carried out by either keeping the protein side chains rigid or allowing flexibility for selected residues, as described in detail in the Materials and Methods section. The best docking poses for both CBC and DMH-CBC, selected on the basis of binding energy value and visual inspection, were subjected to energy minimization ( Figure 3 ). Both enantiomers of CBC and DMH-CBC engage H-bonding with their phenolic hydroxyl and Thr877 (monB). (R)-CBC points its pentyl chain in a hydrophobic cleft formed by Ile953 (monA), Met956 (monB), and Phe880 (monB), while the methyl group of the pentyl chain forms a CH-π interaction with Phe880 and the chromene methyl is surrounded by the hydrophobic residues Ile881 (monB), Leu885 (monB), Phe912 (monB), and Ile945 (monA). Conversely, (S)-CBC adopts a flipped orientation with the pentyl chain now pointing toward Phe912 (monB), engaging a CH-π interaction with the aromatic ring of this amino acid. The same arrangement of (S)-CBC is adopted by both enantiomers of CBC-DMH (8b), with the gem-dimethyls sandwiched between Phe912 (monB) and Ile881 (mon B). In this arrangement, the alkyl chain stabilizes the complex by hydrophobic interactions, counteracting the higher flexibility of the heptyl chain in comparison with the pentyl one and in accordance with the lower activity of the n-heptyl chain analogue, with the length of the CBC pentyl chain being consistent with alternative CH-π interaction with Phe880/Phe912 within the binding pocket. The two DMH-CBC enantiomers share the same orientation of the DMH-chain but show a distinct orientation of the chromene methyl group: in the R enantiomer, the methyl group points toward Phe880, and the terpenoid chain is hosted in a hydrophobic cleft at the interfaces between two monomers. Conversely, in the S enantiomer the orientation of these group is just the opposite. The serendipitous discovery by Adams that a combination of branching and elongation of the alkyl chain boosts the narcotic potency of tetrahydrocannabinols provided the basis for the synthesis of highly potent molecular probes, eventually leading to the characterization of the two cannabinoid receptors (CB1 and CB2) and to the synthesis of nabilone (Cesamet), a clinically efficacious synthetic analogue of ∆ 9 -THC [2] . Despite the relevance of this maneuver for CB 1 and CB 2 affinity, and its potential to bias activity toward specific targets, its effect has not yet been investigated on other high-affinity targets of phytocannabinoids like the thermo-TRPs (TRPV1-V4, TRPM8 and TRPA1) [5] . This is surprising, since one major issue plaguing the bench-to-bed translation of the biological activity of cannabinoids is their promiscuous end-point profile, which spans metabotropic and ionotropic receptors as well as enzymes and transcription factors. For ∆ 9 -THC, the interaction with the two cannabinoid receptors dominates, albeit not fully recapping, the clinical profile, while the other phytocannabinoids show a constellation of end-points characterized by a generally modest affinity for macromolecular end-points, with sub-micromolar potency basically limited to some thermo-TRPs [5] . To systematically address this issue, we have synthesized the α,α-dimethylheptyl analogues of the five major phytocannabinoids and have comparatively evaluated their activity on thermo-TRPs. The results showed that, by quenching affinity for the other thermo-TRPs, the replacement of the linear pentyl substituent with an α,α-dimethylheptyl residue induces a selectivity switch toward TRPA1 in all major phytocannabinoid chemotypes. Molecular docking studies rationalized these results in terms of binding, capitalizing our previously developed homology model of ratTRPA1 [20] . Structural and mutagenic experiments have identified two binding sites for non-electrophilic agonists, centered, respectively, on hThr874 (rThr877, uniprot ID:Q6RI86, Site 1) and hTyr840 (Site 2). However, only Site1 was considered, since, unlike phytocannabinoids and their DMH analogues, the selective Site 2 agonist GNE551 lacks channel desensitizing properties. Moreover, Tyr840, critical for GNE551 binding, is not conserved between human and rat hortologues, being replaced by a Phe residue in the latter. The results of the docking studies suggest the gem-dimethyls on the benzyl carbon provide an additional anchoring site to engage the binding pocket into non-covalent interactions, rationalizing the retention of affinity for TRPA1. This maneuver is detrimental for the other thermos-TRPs presumably because of the incapacity to accommodate the steric and lipophilicity changes associated with benzyl branching into less broadly tuned binding sites [21] . Non-covalent modulators of TRPA1 are relatively few compared to the covalent ones [22] , and their biological profile is critically dependent on their binding mode. Site 2 ligands like GNE551 do not elicit channel desensitization and tachyphylaxis, thus inducing persistent pain insensitive to TRPA1 antagonists [23] . Conversely, Site 1 ligands induce only transient pain, followed by desensitization and analgesia [24] . The desensitization of TRPA1 has the potential to translate into beneficial effects for a wide range of inflammatory conditions [25] , traditionally focused in the realm of inflammation and pain [20] . However, growing evidence of involvement in the SARS-CoV-2 infection through the modulation of inflammation, pain, and fever, has provided an additional therapeutic area of investigation, with documented proof-of-concept clinical cases of the beneficial effects of sulforafane, a TRPA1 and Nrf2-interacting nutrient, on the course of the infection [7] . Replacement of the n-pentyl side chain with an α,α-dimethylheptyl group quenches the exuberant interaction of phytocannabinoids with thermo-TRPs, making them selective non-covalent modulators of TRPA1. Within the DMH-analogues of the most common phytocannabinoids, DMH-CBC (8b) emerged as the most potent compound in the series, and the molecular details of its interaction with this ion channel were investigated and compared to those of the parent compound CBC (8a). Both compounds are racemic, but only the enantiomers of the natural product CBC showed a distinct orientation at TRPA1 binding sites. These results provide a basis for a broader comparison that will also include additional end points of phytocannabinoids and their endogenous analogues, like anandamide and palmitoylethanolamide. Phytocannabinoids: A unified critical inventory The early history of cannabinoid research Isolation, Structure, and Partial Synthesis of an Active Constituent of Hashish Recreational drug discovery: Natural products as lead structures for the synthesis of smart drugs From phytocannabinoids to cannabinoid receptors and endocannabinoids: Pleiotropic physiological and pathological roles through complex pharmacology Natural compounds and synthetic drugs targeting the ion-otropic cannabinoid members of Transient Receptor Potential (TRP) channels Potential Interplay between Nrf2, TRPA1, and TRPV1 Assay of TRPV1 Receptor Signaling AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility Im-proving the accuracy of protein side chain and backbone parameters from ff99SB Development and testing of a general amber force field Synthesis of functionalized canna-binoids Synthetic pathways to tetrahydrocannabinol (THC): An overview Synthesis and antitumor activity of cannabigerol Iodinepromoted aromatization of p-menthane-type phytocannabinoid Effects of canna-binoids and cannabinoid-enriched cannabis extracts on TRP channels and endocannabinoid metabolic enzymes Cannabinoid actions at TRPV channels: Effects on TRPV3 and TRPV4 and their potential relevance to gastrointestinal inflammation Dimeric Phytocannabinoid from Hemp, Cannabis sativa L., Is a Selective Thermo-TRP Modulator Structural Pharmacology of TRP Channels The transient receptor potential channel TRPA1: From gene to patho-physiology A Non-covalent Ligand Reveals Biased Agonism of the TRPA1 Ion Channel N-Cinnamoylanthranilates as hu-man TRPA1 modulators: Structure-activity relationships and channel binding sites Spices: The Savory and Beneficial Science of Informed Consent Statement: Not applicable. The authors declare no conflict of interest. The funder had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.