Discovery of Homobivalent Bitopic Ligands of the Cannabinoid CB2 Receptor**
Medicinal Chemistry Institute, Spanish Research Council, Madrid, Spain
Department of Biochemistry and Physiology, CIBERNED, Faculty of Pharmacy and Food Sciences, Universitat de Barcelona, Barcelona, Spain
Laboratory of Computational Medicine, Biostatistics Unit, Faculty of Medicine, Universitat Autónoma de Barcelona, Barcelona, Spain
Department of Biochemistry and Molecular Biology, CIBERNED, IRYCIS, Faculty of Medicine, Universidad Complutense de Madrid, Madrid, Spain
Department of Biochemistry and Molecular Biology, CIBERNED, School of Chemistry, Universitat de Barcelona, Barcelona, Spain
Present address: Computational Chemistry, Janssen Research & Development, Janssen Pharmaceutica N.V. Belgium
Abstract
Single chemical entities with potential to simultaneously interact with two binding sites are emerging strategies in medicinal chemistry. We have designed, synthesized and functionally characterized the first bitopic ligands for the CB2 receptor. These compounds selectively target CB2 versus CB1 receptors. Their binding mode was studied by molecular dynamic simulations and site‐directed mutagenesis.
Untitled section
Keywords: bitopic ligands, CB2 cannabinoid, G protein-coupled receptors, molecular dynamics, site-directed mutagenesis
Graphical
Single chemical entities with potential to simultaneously interact with two binding sites are emerging strategies in GPCR pharmacology. We present herein the design, synthesis and functional characterization of the first bitopic ligands for the CB2 receptor. Molecular understanding of the binding mode was guided by molecular dynamic simulations and site‐directed mutagenesis.
Boxed Text
Article notes
Untitled section
Received 2020 Jul 19; Revised 2020 Aug 12; Issue date 2020 Dec 4.
G protein‐coupled receptors (GPCRs) regulate a vast amount of cellular processes, [1] thus, they form one of the most important pharmaceutical drug‐target class (475 drugs in the market that represent ∼34 % of all drugs approved by the US Food and Drug Administration). [2] However, these drugs target only 108 unique GPCRs that are <15 % of the ∼800 genes (or ∼30 % of the ∼360 non‐olfactory GPCRs). [2] One of the reasons is that the orthosteric binding site for a particular endogenous ligand is often highly conserved across a GPCR subfamily, thus making it difficult to achieve high selectivity for specific receptor subtypes. [3] Novel approaches to overcome this problem include the discovery of bitopic ligands that bind the orthosteric site as well as a less conserved site within the same receptor unit.[ 4 , 5 , 6 ] This type of complementary cavity is often located at the entrance of the orthosteric binding site, as identified in ligand binding pathway simulations, which have been named extracellular vestibule [7] or entrance, [8] or secondary [9] or metastable [10] binding site, or exosite. [11] In this work, we will name this cavity as receptor vestibule or exosite. Bitopic ligands that target the orthosteric site and the receptor vestibule improve selectivity,[ 11 , 12 ] off‐rates and signaling bias,[ 4 , 13 , 14 ] maintaining bioavailability and brain penetration properties in mice. [15] Other type of comparable ligands are designed to simultaneously bind two orthosteric sites of a (homo/hetero) GPCR dimer. [16] These type of ligands have been recently reviewed. [17]
While most GPCRs recognize polar ligands, GPCRs for lipid mediators are activated by hormone‐like signaling molecules derived from lipid species, which possess long hydrophobic moieties. This subfamily is mostly composed of the sphingosine‐1‐phosphate (S1P), lysophosphatidic acid (LPA) and cannabinoid (CB1R and CB2R) receptors. [18] In the crystal structures of these receptors [19] the extracellular N‐terminus and extracellular loop 2 folds over the ligand binding pocket blocking the access to the orthosteric binding cavity from the extracellular environment. These structures together with binding pathway simulations suggest that binding of lipid‐like ligands to a lipid GPCR occurs through a narrow channel between transmembrane helices (TMs) 1 and 7 that connects the orthosteric binding site to the lipid bilayer. [20] Notably, the access to the ligand pocket of the MT1 melatonin receptor that binds polar ligands (serotonin‐derived compounds) is also via the lipid bilayer. [21] Thus, the design of bitopic ligands for lipid GPCRs is more challenging than for other GPCRs that fully expose the binding site to the extracellular environment, due to the narrow channel linking the binding site and the lipid bilayer.
In the present communication, we have designed bitopic ligands for the cannabinoid CB2R. We have selected CB2R, instead of CB1R, due to its lack of adverse psychotropic effects along with its wide therapeutic application in pathologies such as cancer, neuroinflammation and pain. [22] Bivalent ligands have already been reported for CB1R[ 23 , 24 , 25 ] and CB2R. [26] These ligands were reported before the release of crystal structures, thus, their binding characteristics remain unclear.[ 27 , 28 ] Here, we have used the recently released structure of CB2R in its active [29] Gi‐bound conformation to identify the binding mode of the designed ligands.
The design of bitopic ligands requires the selection of a moiety able to bind the orthosteric site (pharmacophore). In this regard, we have selected chromenopyrazole derivatives A and B (Figure 1), which have been previously identified as CB2R orthosteric agonists. [30] Then, a second pharmacophore unit needs to be developed for the vestibule or exosite. This is challenging because this additional cavity has not been properly characterized yet for most GPCRs. We have taken advantage of the simulated binding process of a lipid inhibitor to the S1P1 receptor. [20] The process consists of the diffusion of the ligand through the bilayer leaflet to contact the vestibule at the top of TM 7 (the rate‐limiting step), subsequently moving from this lipid‐facing vestibule to the orthosteric binding cavity through the channel between TMs 1 and 7. We propose that lipid GPCRs are capable to recognize orthosteric ligands at the vestibule of the receptor. Thus, we have also selected the chromenopyrazole moiety as the second pharmacophore so that the designed bitopic ligands are symmetrical (it contains two copies of the same pharmacophore). In addition, an appropriate length spacer to cover the distance between both pharmacophores is required. Importantly, this approach has been recently supported in the model of the bitopic ligand CTL01‐05‐B‐A05, a symmetrical agomelatine molecule linked by an ethoxyethane spacer, which binds both the orthosteric binding site and the exosite of the MT1 melatonin receptor. [21]
Chromenopyrazole A, in its most stable conformation (Figure S1), was docked into the orthosteric and vestibule sites of CB2R, and its stability was assessed by molecular dynamic (MD) simulations (Figure S2). Results showed that pharmacophore units remain highly stable at the orthosteric site and moderately stable at the exosite. Visual inspection of the models shows that the −CH3 group of the methoxy moieties of both pharmacophores are suitable attachment points to link the spacer moiety. Linker lengths from six to sixteen methylene units were chosen for the bivalent molecules.
Bivalent chromenopyrazoles and their monovalent analogues were synthesized, starting from chromenopyrazoles 4 and 5 (Scheme 1). [31] Preparation of 9‐alkoxychromenopyrazoles (6–17) was achieved in high yields by deprotonation of the hydroxyl group with sodium hydride, followed by rapid addition of an excess of the appropriate 1‐bromoalkane. Bivalent compounds 18–29 were achieved by alkylation of the corresponding chromenopyrazoles with 0.5 equivalents of the desired dibromoalkanes. Different bases were tested, and finally cesium carbonate was selected and used under inert atmosphere. Thus, the desired bivalent compounds (18–29) were achieved in low to moderate yields.
In vitro binding affinities of the alleged bitopic ligands 18–29 (Table 1) and the corresponding monovalent counterparts 6–17 (Table S1) were obtained from [3H]CP‐55,940 competition‐binding assays using membrane fractions of the human CB1R and CB2R, respectively, expressed in HEK‐293T cells. None of the monovalent chromenopyrazoles exhibits affinity towards any of the cannabinoid receptors. Addition of the second pharmacophore makes bivalent ligands capable of binding CB2R in a selective manner. Optimal spacer length for N1‐ and N2‐ethyl derivatives is from 10 (n=8) to 14 (n=12) methylene units.
| Compd | R1 | N [a] | CB1R Ki [μm][b] | CB2R Ki [μm][b] |
|---|---|---|---|---|
| 18 | N1‐Et | 4 | >40 | 28.1±1.6 |
| 19 | N2‐Et | 4 | >40 | 12.4±2.0 |
| 20 | N1‐Et | 6 | >40 | 2.2±0.7 |
| 21 | N2‐Et | 6 | nd | nd |
| 22 | N1‐Et | 8 | >40 | 0.9±0.2 |
| 23 | N2‐Et | 8 | >40 | 5.8±1.5 |
| 24 | N1‐Et | 10 | >40 | 0.4±0.14 |
| 25 | N2‐Et | 10 | >40 | 0.3±0.1 |
| 26 | N1‐Et | 12 | >40 | 0.8±0.1 |
| 27 | N2‐Et | 12 | >40 | 0.3±0.1 |
| 28 | N1‐Et | 14 | >40 | >40 |
| 29 | N2‐Et | 14 | >40 | 2.12±0.21 |
| A30 | N1‐Et | – | 5.0±0.7 | 0.16±0.03 |
| B30 | N2‐Et | – | 2.9±0.5 | 0.09±0.02 |
| WIN [c] | – | – | 0.04±0.01 | 0.003±0.002 |
Compounds with CB2R affinity constants in the low micromolar range (22, 24–27) were selected for functional evaluation by measuring their effect on forskolin‐induced cAMP levels in HEK‐293 cells expressing hCB2R (Figure 2 A). Dose‐response experiments demonstrated that these compounds are able to inhibit cAMP accumulation as efficiently as CP55,940 but with a slight drop of potency (27: pEC50=7.6 vs. CP55,940: pEC50=8.2) (Table 2), while their monovalent counterparts are inactive (Figure S3). None of the compounds displayed an effect in non‐transfected HEK‐293T cells, confirming that the results are fully mediated by CB2R (Figure S4).
| Wild type | Val361.35Met | Ala2827.36Met | ||||
|---|---|---|---|---|---|---|
| pEC50 | E max [a] | pEC50 | E max [a] | pEC50 | E max [a] | |
| 22 | 7.1±0.2 | 42±5.4 | 7.7±0.1 | 75±1.5 | 7.2±0.2 | 77±2.5 |
| 24 | 7.2±0.2 | 45±6.3 | nd | nd | ||
| 25 | 7.2±0.2 | 40±5.9 | nd | nd | ||
| 26 | 6.5±0.2 | 52±4.1 | nd | nd | ||
| 27 | 7.6±0.1 | 44±3.7 | 7.1±0.2 | 68±3.5 | 7.3±0.2 | 84±1.4 |
| CP [b] | 8.2±0.1 | 48±1.7 | 8.1±0.1 | 50±2.3 | 7.2±0.1 | 48±3.43 |
In order to assess the binding mode of homobivalent chromenopyrazoles 22, 25 and 27 (10, 12 and 14 methylene units), we docked their most favorable conformation (Figure S1) into CB2R in its active state in such a manner that both pharmacophoric units bind into the orthosteric and vestibule sites (Figure 3 A). Unbiased MD simulations show that these chain lengths can simultaneously bind both sites (Figure 3, S5). The methylene spacer expands toward the lipid‐facing vestibule interacting with the hydrophobic side chains of Phe912.61, Ala2827.36, Met2867.40, and Val361.35. In the vestibule, the chromenopyrazole moiety forms aromatic–aromatic interactions with Phe2837.37. In the simulations Phe2837.37 adopts the trans conformation that opens the channel between TMs 1 and 7 and permits the ligand to reach the membrane, in contrast to the gauche + conformation observed in the crystal structure of CB2R[ 29 , 32 ] that closes the channel. The heptyl chain is accommodated in a hydrophobic cavity of TMs 1 and 7, facing the membrane, which is formed by Leu391.38, Cys401.39, Met2867.40, Leu431.39, Ile2907.44, Met2937.47 and Leu461.45. In the case of the N1‐ethyl derivative the lone pair in the N2 atom forms a hydrogen bond with Gln321.31. A detailed description of these interactions is shown in Figure S6.
Figure 3 B shows the mesh surface formed by Val361.35 and Ala2827.36. Clearly, these TMs side chains, which are located midway between the orthosteric site and the receptor vestibule, delimit the channel between TMs 1 and 7. Thus, in order to validate the proposed binding mode of bitopic ligands, we mutated the side chains of Val361.35 (Figure 2 B) and Ala2827.36 (Figure 2 C) to the much larger Met side chain. As expected, these mutations do not influence the function of the orthosteric agonist CP55,490, but clearly impair signaling of bitopic ligands 22 and 27 (Table 2) by occupying the volume of the channel. Notably, these results contrast with other MD simulations, suggesting the entrance of the ligands from the extracellular environment.[ 33 , 34 ]
In summary, we present herein the discovery of CB2R homobivalent bitopic ligands. These compounds were designed as symmetrical bivalent compounds using the previously reported chromenopyrazole scaffold as potential pharmacophore for both units linked by a methylene spacer. Binding and functional cAMP studies revealed their ability to selectively activate CB2R versus CB1R. The longer Ile1.35 in CB1R than Val1.35 in CB2R seems responsible for the observed selectivity (Figure S7). MD simulations and site‐directed mutagenesis studies show that these bitopic ligands bind into the orthosteric site and in a vestibule/exosite located at the ligands entry/egress channel that connects the orthosteric site with the lipid bilayer membrane. Whether these bitopic ligands at CB2R show beneficial therapeutic application needs further investigation such as structure–activity relationship studies to improve potency.
Conflict of interest
The authors declare no conflict of interest.
Supporting information
Acknowledgements
This work was supported by the MICIU/FEDER (RTI2018‐095544‐B‐I00, PID2019‐109240RB‐I00, SAF2017‐84117‐R), and by the CSIC (PIE‐201580E033). P.M. is grateful to the CM program “Atracción de Talento Investigador” number 2018‐T2/BMD‐1081.
Untitled section
P. Morales, G. Navarro, M. Gómez-Autet, L. Redondo, J. Fernández-Ruiz, L. Pérez-Benito, A. Cordomí, L. Pardo, R. Franco, N. Jagerovic, Chem. Eur. J. 2020, 26, 15839.
Untitled section
Footnote Group
Contributor Information
Prof. Dr. Leonardo Pardo, Email: leonardo.pardo@uab.es.
Prof. Dr. Rafael Franco, Email: rfranco@ub.edu.
Dr. Nadine Jagerovic, Email: nadine@iqm.csic.es.
References
Untitled section
References
- 1. Wootten D., Christopoulos A., Marti-Solano M., Babu M. M., Sexton P. M., Nat. Rev. Mol. Cell Biol. 2018, 19, 638–653.
- 2. Hauser A. S., Attwood M. M., Rask-Andersen M., Schiöth H. B., Gloriam D. E., Nat. Rev. Drug Discovery 2017, 16, 829–842.
- 3. Foster D. J., Conn P. J., Neuron 2017, 94, 431–446.
- 4. Lane J. R., Sexton P. M., Christopoulos A., Trends Pharmacol. Sci. 2013, 34, 59–66.
- 5. Mohr K., Schmitz J., Schrage R., Tränkle C., Holzgrabe U., Angew. Chem. Int. Ed. 2013, 52, 508–516; Angew. Chem. 2013, 125, 530–538.
- 6. Bradley S. J., Tobin A. B., Annu. Rev. Pharmacol. Toxicol. 2016, 56, 535–559.
- 7. Dror R. O., Pan A. C., Arlow D. H., Borhani D. W., Maragakis P., Shan Y., Xu H., Shaw D. E., Proc. Natl. Acad. Sci. USA 2011, 108, 13118–13123.
- 8. Wang C., Jiang Y., Ma J., Wu H., Wacker D., Katritch V., Han G. W., Liu W., Huang X. P., Vardy E., McCorvy J. D., Gao X., Zhou X. E., Melcher K., Zhang C., Bai F., Yang H., Yang L., Jiang H., Roth B. L., Cherezov V., Stevens R. C., Xu H. E., Science 2013, 340, 610–614.
- 9. González A., Perez-Acle T., Pardo L., Deupi X., PLoS One 2011, 6, e23815.
- 10. Fronik P., Gaiser B. I., Sejer Pedersen D., J. Med. Chem. 2017, 60, 4126–4134.
- 11. Masureel M., Zou Y., Picard L. P., van der Westhuizen E., Mahoney J. P., Rodrigues J. P. G. L. M., Mildorf T. J., Dror R. O., Shaw D. E., Bouvier M., Pardon E., Steyaert J., Sunahara R. K., Weis W. I., Zhang C., Kobilka B. K., Nat. Chem. Biol. 2018, 14, 1059–1066.
- 12. Medina R. A., Vázquez-Villa H., Gómez-Tamayo J. C., Benhamú B., Martín-Fontecha M., De La Fuente T., Caltabiano G., Hedlund P. B., Pardo L., López-Rodríguez M. L., J. Med. Chem. 2014, 57, 6879–6884.
- 13. Valant C., Lane J. R., Sexton P. M., Christopoulos A., Annu. Rev. Pharmacol. Toxicol. 2012, 52, 153–178.
- 14. Gaiser B. I., Danielsen M., Marcher-Rørsted E., Røpke Jørgensen K., Wróbel T. M., Frykman M., Johansson H., Bräuner-Osborne H., Gloriam D. E., Mathiesen J. M., Sejer Pedersen D., J. Med. Chem. 2019, 62, 7806–7839.
- 15. Tan L., Zhou Q., Yan W., Sun J., Kozikowski A. P., Zhao S., Huang X., Cheng J., J. Med. Chem. 2020, 63, 4579–4602.
- 16. Pulido D., Casadó-Anguera V., Pérez-Benito L., Moreno E., Cordomí A., López L., Cortés A., Ferré S., Pardo L., Casadó V., Royo M., J. Med. Chem. 2018, 61, 9335–9346.
- 17. Newman A. H., Battiti F. O., Bonifazi A., J. Med. Chem. 2020, 63, 1770–1779.
- 18. Pándy-Szekeres G., Munk C., Tsonkov T. M., Mordalski S., Harpsøe K., Hauser A. S., Bojarski A. J., Gloriam D. E., Nucleic Acids Res. 2018, 46, D440–D446.
- 19. Audet M., Stevens R. C., Protein Sci. 2019, 28, 292–304.
- 20. Stanley N., Pardo L., De Fabritiis G., Sci. Rep. 2016, 6, 22639.
- 21. Stauch B., Johansson L. C., McCorvy J. D., Patel N., Han G. W., Huang X. P., Gati C., Batyuk A., Slocum S. T., Ishchenko A., Brehm W., White T. A., Michaelian N., Madsen C., Zhu L., Grant T. D., Grandner J. M., Shiriaeva A., Olsen R. H. J., Tribo A. R., Yous S., Stevens R. C., Weierstall U., Katritch V., Roth B. L., Liu W., Cherezov V., Nature 2019, 569, 284–288.
- 22. Dhopeshwarkar A., Mackie K., Mol. Pharmacol. 2014, 86, 430–437.
- 23. Zhang Y., Gilliam A., Maitra R., Damaj M. I., Tajuba J. M., Seltzman H. H., Thomas B. F., J. Med. Chem. 2010, 53, 7048–7060.
- 24. Fernández-Fernández C., Decara J., Bermúdez-Silva F. J., Sánchez E., Morales P., Gõmez-Cañas M., Gõmez-Ruíz M., Callado L. F., Goya P., Rodríguez De Fonseca F., Martín M. I., Fernández-Ruíz J., Meana J. J., Jagerovic N., Arch. Pharm. 2013, 346, 171–179.
- 25. Huang G., Pemp D., Stadtmüller P., Nimczick M., Heilmann J., Decker M., Bioorg. Med. Chem. Lett. 2014, 24, 4209–4214.
- 26. Nimczick M., Pemp D., Darras F. H., Chen X., Heilmann J., Decker M., Bioorg. Med. Chem. 2014, 22, 3938–3946.
- 27. Glass M., Govindpani K., Furkert D. P., Hurst D. P., Reggio P. H., Flanagan J. U., Trends Pharmacol. Sci. 2016, 37, 353–363.
- 28. Pérez-Benito L., Henry A., Matsoukas M. T., Lopez L., Pulido D., Royo M., Cordomí A., Tresadern G., Pardo L., Bioinformatics 2018, 34, 3857–3863.
- 29. Xing C., Zhuang Y., Xing C., Zhuang Y., Xu T., Feng Z., Zhou X. E., Chen M., Wang L., Cell 2020, 180, 645–654.
- 30. Morales P., Gómez-Cañas M., Navarro G., Hurst D. P., Carrillo-Salinas F. J., Lagartera L., Pazos R., Goya P., Reggio P. H., Guaza C., Franco R., Fernández-Ruiz J., Jagerovic N., J. Med. Chem. 2016, 59, 6753–6771.
- 31. Cumella J., Hernández-Folgado L., Girón R., Sánchez E., Morales P., Hurst D. P., Gómez-Cañas M., Gómez-Ruiz M., Pinto D. C. G. A., Goya P., Reggio P. H., Martin M. I., Fernández-Ruiz J., Silva A. M. S., Jagerovic N., ChemMedChem 2012, 7, 452–463.
- 32. Li X., Hua T., Vemuri K., Ho J. H., Wu Y., Wu L., Popov P., Benchama O., Zvonok N., Locke K., Qu L., Han G. W., Iyer M. R., Cinar R., Coffey N. J., Wang J., Wu M., Katritch V., Zhao S., Kunos G., Bohn L. M., Makriyannis A., Stevens R. C., Liu Z. J., Cell 2019, 176, 459–467.
- 33. Hua T., Vemuri K., Pu M., Qu L., Han G. W., Wu Y., Zhao S., Shui W., Li S., Korde A., Laprairie R. B., Stahl E. L., Ho J. H., Zvonok N., Zhou H., Kufareva I., Wu B., Zhao Q., Hanson M. A., Bohn L. M., Makriyannis A., Stevens R. C., Liu Z. J., Cell 2016, 167, 750–762.
- 34. Saleh N., Hucke O., Kramer G., Schmidt E., Montel F., Lipinski R., Ferger B., Clark T., Hildebrand P. W., Tautermann C. S., Angew. Chem. Int. Ed. 2018, 57, 2580–2585; Angew. Chem. 2018, 130, 2610–2615.