Computational and Experimental Drug Repurposing of FDA-Approved Compounds Targeting the Cannabinoid Receptor CB1
Department of Experimental Medicine, Tor Vergata University of Rome, Via Montpellier 1, 00121 Rome, Italy; emanuele.criscuolo@alumni.uniroma2.eu (E.C.); catalin.nicoara@alumni.uniroma2.eu (C.N.)
Department of Medicine, Campus Bio-Medico University of Rome, Via Alvaro del Portillo 21, 00128 Rome, Italy; marialaura.desciscio@uniroma1.it (M.L.D.S.); angela.decristofaro@alumni.uniroma2.eu (A.D.C.)
Department of Biotechnological and Applied Clinical Sciences, University of L’Aquila, Via Vetoio, Coppito, 67100 L’Aquila, Italy
European Center for Brain Research/Santa Lucia Foundation IRCCS, Via Del Fosso di Fiorano 64, 00143 Rome, Italy
Abstract
The cannabinoid receptor 1 (CB1R) plays a pivotal role in regulating various physiopathological processes, thus positioning itself as a promising and sought-after therapeutic target. However, the search for specific and effective CB1R ligands has been challenging, prompting the exploration of drug repurposing (DR) strategies. In this study, we present an innovative DR approach that combines computational screening and experimental validation to identify potential Food and Drug Administration (FDA)-approved compounds that can interact with the CB1R. Initially, a large-scale virtual screening was conducted using molecular docking simulations, where a library of FDA-approved drugs was screened against the CB1R’s three-dimensional structures. This in silico analysis allowed us to prioritize compounds based on their binding affinity through two different filters. Subsequently, the shortlisted compounds were subjected to in vitro assays using cellular and biochemical models to validate their interaction with the CB1R and determine their functional impact. Our results reveal FDA-approved compounds that exhibit promising interactions with the CB1R. These findings open up exciting opportunities for DR in various disorders where CB1R signaling is implicated. In conclusion, our integrated computational and experimental approach demonstrates the feasibility of DR for discovering CB1R modulators from existing FDA-approved compounds. By leveraging the wealth of existing pharmacological data, this strategy accelerates the identification of potential therapeutics while reducing development costs and timelines. The findings from this study hold the potential to advance novel treatments for a range of CB1R -associated diseases, presenting a significant step forward in drug discovery research.
Untitled section
Keywords: cannabinoid receptor 1, drug repurposing, structure-based virtual screening
Article notes
Untitled section
Received 2023 Oct 31; Revised 2023 Nov 25; Accepted 2023 Nov 30; Collection date 2023 Dec.
1. Introduction
Structure-based virtual screening (SBVS) is a computational method for early stage drug discovery starting from novel bioactive molecules [1,2]. Its application appears more efficient than traditional drug discovery approaches, which are often rather complex, expensive and risky.
To be applied, SBVS needs available three-dimensional structures of the proteins of interest, which can be obtained from different experimental techniques such as X-ray diffraction and cryogenic electron microscopy (cryo-EM), and more recently Alphafold—an artificial intelligence that predicts the protein conformation from the amino acid sequence [3]. Then, SBVS can be used for a potent drug discovery technique named drug repurposing (DR), which starts from existing drugs that were commercialized for different therapeutic indications. Of note, DR increases the success rate of drug development, a process that takes more than 10 years and EUR 1 billion per drug, with a high (~90%) failure rate in clinical trials [4].
The combination of SBVS and DR leads to structure-based drug repurposing (SBDR), a technique that allows us to discover potentially new drugs from virtual libraries of approved compounds, by uncovering their interactions with selected proteins. The prediction of drug–protein interactions and the analysis of their binding affinity values can be exploited in different computational methods, such as molecular docking, alchemical binding free energy calculation, umbrella sampling and molecular mechanics/generalized Born and surface area solvation (MM-GBSA) or molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) [5].
DR also promises high approval rates, because the safety of any commercially available drug has been already assessed in preclinical and clinical trials; thus, the time (and money) needed for drug development can be reduced using this approach [6,7].
Aspirin is the oldest example of DR. Initially formulated by Bayer in 1899 as an analgesic, it underwent repurposing several decades later as an antiplatelet aggregation drug [8,9]. Then, aspirin was proposed to treat other conditions such as cancer [10,11] and cardiovascular diseases [12]. Thalidomide is another example of DR. It was first marketed in 1957 for the treatment of anxiety, sleeping trouble and morning sickness. Yet, it was soon withdrawn because of its teratogenicity [13,14,15], but later on, it was repurposed to treat cancers—particularly erythema nodosum leprosum, multiple myeloma and other hematological malignancies [16]. Maybe the most famous example of DR is sildenafil, that was developed for angina but then repurposed for erectile dysfunction [17]. Finally, among repurposed drugs worth mentioning, one can list baricitinib, remdesivir and tocilizumab, originally developed to treat Alopecia Areata, viral infections (SARS, MERS and AIDS) and rheumatoid arthritis, respectively [18,19,20], but then repurposed as potential COVID-19 treatments [21,22].
Here, we sought to use SBDR for repurposing Food and Drug Administration (FDA)-approved drugs on cannabinoid receptor 1 (CB1R), one of the most prominent G protein-coupled receptors in the human brain [23,24,25]. To this end, knowing the 3D structure details of CB1R is extremely important. This receptor shows the classical 7 transmembrane fold with a lid above the binding site that contains abundant acidic residues and a highly hydrophobic orthosteric binding pocket [26,27]. In addition, a ‘twin toggle switch’ of Phe200 and Trp356 appears essential for CB1R activation [28]. Overall, CB1R is considered a major pharmacological target, due to its many implications for diseases of the central nervous system (CNS) and peripheral organs [24,29]. Unsurprisingly, CB1R is widely distributed throughout the body and is highly expressed within the CNS (approximately 10 to 50 times more than receptors of classic neurotransmitters like opioid and dopamine receptors) in the basal ganglia, hippocampus, cerebellum, amygdala, cingulate cortex, medial hypothalamus and spinal cord [30]. Moreover, CB1R is present—though to a lesser extent—in the periphery, namely in adipose tissue, liver, skeletal muscles, kidney and pancreas [21,31,32,33,34].
CB1R was first discovered as the target of the main psychoactive ingredient of cannabis extracts, Δ9-tetrahydrocannabinol [24,35,36]. Later on, it was recognized that CB1R is a key element of a complex lipid signaling system, called the endocannabinoid (eCB) system [24,37]. The latter comprises eCBs, that are endogenous ligands of CB1R [24,37]. The most relevant eCBs derive from arachidonic acid, anandamide (N-arachidonoylethanolamine, AEA), which is an amide, and 2-arachidonoylglycerol (2-AG), which is an ester. In addition to CB1R, AEA and 2-AG have other receptor targets and a number of metabolic enzymes that have been recently reviewed in detail [24]. Here, only the main metabolic enzymes of 2-AG (biosynthesis: diacylglycerol lipases (DAGL) α and β; degradation: monoacylglycerol lipase (MAGL) and α/β hydrolase domain-containing (ABHD) proteins 6 and 12) [24,38] have been analyzed, along with fatty acid amide hydrolase (FAAH), the main catabolic enzyme of AEA [24,37].
During the last decade, hundreds of compounds able to bind to CB1R have been synthesized and tested, reporting diverse pharmacological effects [24,39,40,41]. Unfortunately, most of them have been withdrawn because of undesirable side effects. Noteworthy seems the case of SR141716A (SR1, also known as Rimonabant or Acomplia®). This selective antagonist/inverse agonist of CB1R was first approved in Europe in 2006 for the management of obesity, but it was then withdrawn two years later because of adverse effects like an increased incidence of depression and suicidal ideation [42,43]. Another potent synthetic CB1R agonist, called AMB-Fubinaca, has been named the “zombie drug” in 2016, because 33 people (25–59 year old) were adversely affected with a semicomatose state [44,45]. These examples simply remind us that there is still an urgent need to find effective drugs able to modulate CB1R with minimal (if at all) side effects.
Here, SBDR has been used on two different 3D structures of CB1R, in order to increase the accuracy of the calculation taking into full account the orthosteric site flexibility of the receptor. Indeed, CB1R shows different conformations upon binding lead by molecular properties. Significant structural modifications due to the ligand nature are evident in helices I and II [28]. In particular, helix I undergoes an inward bending of approximately 6.6 Å, while helix II rotates inwardly by approximately 6.8 Å compared to the agonist-bound and the antagonist-bound states [28]. Likewise, notable conformational changes are observed in the cytoplasmic section of the receptor, where helix VI moves outward by about 8 Å. As a result of the inward shifts of helices I and II, the ligand-binding pocket’s volume reduces from 822 Å3 in the antagonist-bound structure to 384 Å3 in the agonist-bound complex, signifying a substantial 53% reduction [28,46].
In pursuit of identifying commercial drugs capable of interacting with CB1R, an FDA-approved list of 1379 molecules was retrieved from a chemical database (Zinc15 database). To further validate the in silico results, in vitro competitive radioligand assays were also performed and the experimental binding affinity of the compounds were calculated. Radioligand binding assays are indeed a powerful tool in the early phase of the drug design and discovery process, because they allow us to study directly ligand–receptor interactions [47]. Then, a proteomic technique was used to evaluate the possible interaction of the most potent CB1R ligands towards other key elements of the eCB system, in order to predict possible side effects.
2. Results
2.1. Virtual Screening
To estimate the binding affinity of the investigated drugs towards CB1R, a molecular docking analysis was performed by means of the MOE (Molecular Operating Environment) software (Chemical Computing Group (CCG), Montreal, QC H3A 2R7, Canada). To this end, two different 3D structures of CB1R were used, with PDB codes 5XRA (at 2.80 Å resolution) and 5U09 (at 2.60 Å resolution) in the Protein Data Bank (PDB) (www.rcsb.org). In particular, the first structure is co-crystallized with the receptor agonist AM11542 (Ki = 0.29 nM), whereas the second is complexed with the inverse agonist Taranabant (Ki = 0.13 nM) [27,28]. The performance of the docking procedure was ascertained by redocking the co-crystallized ligands, AM11542 and Taranabant, that yielded ΔG values of −10.951 kcal/mol and −11.134 kcal/mol, respectively (Figure 1). All binding free variations were obtained through the use of Equation (1) computed using MOE.
In this context, it seems important to recall that CB1R’s conformation is affected by the nature of the ligands [26]. Therefore, cross-docking has also been performed by analyzing the binding affinity of AM11542 for 5U09 and Taranabant for 5XRA, that yielded ΔG = −9.6540 kcal/mol and ΔG = −8.2429 kcal/mol, respectively (Figure 1). The different ΔG values of each CB1R ligand in both structures, −10.951 kcal/mol and −9.6540 kcal/mol for AM11542 and −11.134 kcal/mol and −8.2429 kcal/mol for Taranabant, confirmed the different conformation of the agonist-bound state and the inverse agonist-bound state, with a decrease in the binding affinity (ΔG) of −1.480 and −2.708 kcal/mol, respectively (Figure 1).
Then, the ZINC15 database (https://zinc.docking.org/) of 1379 molecules was used [48], in order to perform virtual screening on both 3D structures of CB1R through a rigid-protein docking approach. All details are reported in the Methods section. In this first analysis, a threshold of the binding affinity of −8.5 kcal/mol has been established in order to obtain a restricted set of 200 molecules.
Next, to reveal only potential CB1R agonists among the selected molecules, all subsequent analyses were carried out on the 5XRA structure [2]. Afterwards, a blind induced-fit docking, localized in the channel that leads to the orthosteric site of the 5XRA structure, was performed using an induced-fit approach, which allows the binding site to move freely. For this purpose, a second filter stage, that selects only the molecules with several high affinity poses close to the binding site, was applied. Figure 2 shows a representative example of the comparison between AM11542 (a) and Raloxifene (b), with all poses into the channel, and a discarded compound (mupirocin) (c), with several scattered predicted poses.
These channel analyses may provide an explanation for the different potency and selectivity of the CB1R ligands [49,50]. With this filter, 10 drugs were selected: Aminopterin (APGA), Avanafil, Ceftriaxone, Methotrexate, Miltefosine, PGE-1, Raloxifene, Raltegravir, Riociguat and Valsartan, as shown in Table 1. Interestingly, among these compounds, Methotrexate and APGA show similar structures, also shown in Table 1 for comparison.
| Drugs | Chemical Structures | Indications |
|---|---|---|
| Aminopterin (APGA) | Pediatric leukemia | |
| Avanafil | Erectile dysfunction | |
| Ceftriaxone | Bacterial infections, such as endocarditis, meningitis, pneumonia, skin infections, urinary tract infections | |
| Methotrexate | Cancer, Autoimmune diseases, Ectopic pregnancies | |
| Miltefosine | Leishmaniasis, Breast cancer (topical treatment) | |
| PGE-1 | Erectile dysfunction, Neonatal congenital heart defects | |
| Raloxifene | Osteoporosis, Breast cancer prevention | |
| Raltegravir | HIV | |
| Riociguat | Pulmonary hypertension | |
| Valsartan | Hypertension, Heart failure, Diabetes, Kidney disease |
The selected compounds have different chemical structures, and expectedly, they have been marketed for different diseases such as cancer, hypertension, HIV, leukemia, erectile dysfunction and osteoporosis (Table 1). Of note, one of the selected drugs, Raloxifene, has been already reported to bind to CB1R with Ki = 210 nM [51].
2.2. Analysis of CB1R Binding
With the aim of validating in silico computational data, competitive radioligand binding assays were performed at 10 µM to estimate the drug potency. The potent CB1R synthetic cannabinoid [3H]CP55,940 was used as the ligand to be displayed in mouse brain membranes. In the same experiments, Raloxifene was used as a positive control [51].
Eight of the ten selected drugs significantly competed with [3H]CP55,940 (Figure 3).
In particular, the most potent appeared to be Raltegravir, that produced ~80% displacement of [3H]CP55,940. Interesting, Raloxifene was less potent than Raltegravir (Figure 3). Furthermore, Miltefosine and Methotrexate displaced ~60% and ~50% [3H]CP55,940, respectively. Instead, Valsartan, Avanafil, APGA and Riociguat were found to be weak competitors, that produced ~30% (the first three) and ~20% displacement of [3H]CP55,940, respectively (Figure 3). Finally, Ceftriaxone and PGE-1 failed to show any competition at CB1R (Figure 3), which in the case of PGE-1 could be due to its poor stability [52].
2.3. Activity-Based Protein Profiling
Competitive activity-based protein profiling (ABPP) was performed to assess potential interactions of the most potent drugs—Raltegravir, Methotrexate and Miltelfosine—excluding the positive control Raloxifene, with some of the key enzymes involved in eCB metabolism. ABPP is a powerful technique that takes advantage of chemical probes able to react with an amino acid located in the catalytic site of the target enzymes [53,54]. Here, membrane and cytosol fractions of mouse brain were used to analyze the interactions of FDA drugs with the following serine hydrolases: DAGLα/β, MAGL, ABHD6, ABHD12 and FAAH (Figure 4).
To this end, two ABPP probes (MB064 and FP-Bodipy) were used at 100 nM and 10 µM, respectively [55]. Moreover, inhibitors of DAGLα/β (DH376), MAGL and ABHD6 (ABX1431), ABHD12 (DO264) and FAAH (PF04458745) were used as positive controls. None of the tested drugs were able to interact with the selected enzymes, supporting the lack of off-targets (Figure 4).
2.4. Cell Viability
The mechanisms of action of FDA-approved drugs used in this study are known [56,57], except for Miltefosine, whose molecular targets have not been identified yet. However, this drug is known to inhibit the proliferation of human keratinocytes (HaCaT cells), at micromolar concentrations [58]. In line with this, and keeping in mind the role of CB1R in cancer [59], here, the anti-proliferative effects of Miltefosine were investigated in the micromolar range (Figure 5a).
Furthermore, the effects were analyzed at different incubation times (2, 4 and 6 h), in order to assess possible time-dependence (Figure 5b). Miltefosine failed to affect cell growth at lower concentrations and at any time of incubation. Instead, at 25 µM, it showed anti-proliferative effects that were time-dependent (Figure 5a,b). Remarkably, the CB1R antagonist SR1 attenuated the effect of Miltefosine, in a manner that was directly proportional to the preincubation time (Figure 5b). These results demonstrate that Miltefosine exerts an anti-proliferative effect via the engagement of CB1R.
3. Discussion
DR represents an efficient approach to drug discovery because it exploits approved drugs with a known safety profile. The main advantage of DR is the marked reduction in the costs of clinical trials, which account for ~60% of the total cost of drug development. Indeed, DR can lead to new therapeutic indications for existing drugs with extensive clinical history and toxicology.
Here, virtual screening performed on two different 3D structures of CB1R allowed us to select 200 molecules out of 1379 compounds reported in the ZINC15 database. Then, 10 compounds were isolated through the application of a second filter to find possible agonists for CB1R using the 5XRA structure (Table 1). To validate the screening method employed, one of the final compounds has previously been reported to bind to CB1R. This selective estrogen receptor modulator, Raloxifene, exhibited a Ki value of 210 nM on CB1R [51]. Additionally, the authors investigated various compounds belonging to five structurally distinct classes of selective estrogen receptor modulators [51]. Notably, all the tested molecules demonstrated weaker affinities for the receptor compared to Raloxifene, such as Bazedoxifene, Nafoxidine and Ospemifene. Importantly, the initial filtering stage selected only the most potent compound among the estrogen receptor modulators. Indeed, the chosen threshold of −8.5 kcal/mol discarded Ospemifene with a ΔG of −8.28 kcal/mol, as well as Nafoxidine and Bazedoxifene with a ΔG ~ −8 kcal/mol. Following the computational analysis, an experimental competitive binding assay was conducted to validate the in silico findings. Specifically, eight drugs displayed the ability to displace the potent CB1R agonist, CP55,940, with a moderate or high affinity (Figure 3) and achieving a prediction accuracy of 80%. Among them, Methotrexate, Miltefosine and Raltegravir as well as the positive control, Raloxifene, showed the most significant displacement of [3H]CP55,940, exceeding 50% (Figure 3). The doses and characteristics of these compounds are detailed in Table 2.
| Drugs | Commercial Name | Dose | Features |
|---|---|---|---|
| Methotrexate | Trexall, Otrexup, Rasuvo, Xatmep, RediTrex, Jylamvo | Dosage depends on pathology | Methotrexate is a folate derivative that inhibits several enzymes responsible for nucleotide synthesis. It is used to treat inflammation caused by arthritis or to control cell division in neoplastic diseases. |
| Miltefosine | Impavido | 50 mg BID/TID | Miltefosine is a broad spectrum antimicrobial, anti-leishmanial, phospholipid drug developed in the 1980s as an anti-cancer agent. |
| Raltegravir | Isentress | 600 mg BID | Raltegravir is an antiretroviral agent used for the treatment of HIV infections. It is the first of a new class of HIV drugs. |
Furthermore, to explore potential interactions between the selected drugs and representative eCB system enzymes and the related ones, a proteomic analysis of ABHD6/12, DAGL, FAAH and MAGL was performed. The results revealed that the tested compounds did not interact with any of these enzymes (Figure 4). Although the analyzed drugs are highly established in treating specific diseases [56,57], their interaction with CB1R may extend therapeutic actions or elucidate adverse effects through synergistic effects. For instance, the effects of Raloxifene on estrogen receptors and CB1R could account for both beneficial actions and potential adverse events [51]. The binding assay demonstrates the ability of these compounds to interact with CB1R; however, it is essential to consider the dosage for each drug. Raltegravir, the most potent compound, is used in combination with other medications to treat human immunodeficiency virus (HIV) infections (Table 1 and Table 2) [60]. It belongs to the class of medications called HIV integrase inhibitors, and the initial regimen involves 400 mg (0.9 mmol) tablets taken twice daily, with a relatively high concentration (Table 2). Some reported side effects include insomnia, abnormal dreams, loss of appetite, headache, nausea, vomiting, fast heartbeat and depression, which may partially be attributed to CB1R involvement. Methotrexate, on the other hand, belongs to the class of medications called antimetabolites, and its dosage depends on the specific pathology being treated. It is used to slow the growth of cancer cells in certain types of cancer and is also employed to treat severe psoriasis and rheumatoid arthritis by reducing the immune system activity (Table 1 and Table 2). Interestingly, THC has shown similar effects to Methotrexate in neutralizing the inflammatory process [61]. Daily oral administration of THC for 21 days in arthritic rats was well tolerated, without causing significant psychoactive side effects; simultaneously, it attenuated the severity of clinical manifestations [61]. Notably, patients with rheumatoid arthritis and osteoarthritis have high expression levels of CB1R [61]. The interaction between Methotrexate and CB1R highlighted in this study could be involved in the compound’s capacity to treat psoriasis and rheumatoid arthritis.
Among the tested drugs, Miltefosine is particularly interesting as a classic multi-target drug. Originally developed in the 1980s as an antineoplastic agent, it is now marketed for leishmaniasis and has been approved for human treatment since 2002 [62], due to repurposing studies. Another noteworthy feature of Miltefosine is its amphiphilic nature, bearing remarkable similarity to endocannabinoids, hence explaining its ability to bind to CB1R (Table 1). Additionally, Miltefosine has shown to inhibit the proliferation of HaCaT cells [58]. One of the goals achieved in this study lies in establishing a clear relationship between Miltefosine’s ability to bind to CB1R and its effect on HaCaT cells. Notably, this cell line expresses the full eCB system [63]. The data illustrate that the selective antagonist/inverse agonist SR1 inhibits the antiproliferation effect of Miltefosine in a time-dependent manner (Figure 5b), demonstrating that its mechanism of action is mediated by CB1R. Miltefosine is administered orally at a dosage of 50 mg (0.122 mmol) three times a day, with common side effects including decreased appetite, diarrhea and vomiting.
CB1R plays a critical role in several human physiological and pathological conditions. Therefore, extensive endeavors have been undertaken to develop ligands for treating various diseases; nevertheless, none of these ligands are currently utilized in clinical settings.
4. Materials and Methods
4.1. Materials
Chemicals were of the purest analytical grade (>95%). [3H]CP55,940 (126 Ci/mmol) was from Perkin Elmer Life Sciences, Inc. (Boston, MA, USA). Avanafil, Aminopterin, Ceftriaxone, CP55,940, Dofetilide, Methotrexate, Miltefosine, Prostaglandin E1, Raloxifene, Raltegravir, Riociguart and Valsartan were from Cayman Chemical Company, Ann Arbor, MI, USA. Human keratinocytes (HaCaT cells) were provided by CLS Cell Lines Service GmbH, (Eppelheim, Germany).
4.2. Virtual Screening
The crystal structure of CB1R, with a potent co-crystallized agonist AM11542 (PDB code: 5XRA), was retrieved from the PDB (www.rcsb.org), and the Experimental Data Snapshot and the PDB validation were analyzed: resolution (2.80 Å), RFree (0.254), Clashscore (6), Ramachandran outliers (0), sidechain outliers (0.6%) and RSRZ outliers (4.7%). The CB1R crystal structure with a potent inverse agonist Taranabant (PDB code: 5U09) was retrieved from the PDB (www.rcsb.org), and the Experimental Data Snapshot and the PDB validation were analyzed: resolution (2.60 Å), RFree (0.247), Clashscore (14), Ramachandran outliers (0), sidechain outliers (0) and RSRZ outliers (0.8%).
4.2.1. Proteins and Ligands Preparation
The simulation studies were carried out using the cutting-edge Molecular Operating Environment (MOE 2021.0102) by the Chemical Computing Group (2021). For the selected crystal structure, the “Structure Preparation” panel with the “Protonate 3D” function was employed to optimize the ionization states of the added hydrogen atoms.
Atoms further away than 8 Å from the co-crystalized ligand were fixed while constraints were applied to atoms within the active site before the final minimization step to refine the protein structure. The 3D structures of the molecules under investigation were sourced from the Zinc Database (https://zinc.docking.org/) [48]. Ligands in ZINC15 were prepared (protonation and tautomers generation at physiological pH) through ChemAxon’s Jchem while Omega was used to obtain the conformations for docking.
4.2.2. Re-Docking Validation
To validate the molecular docking algorithm, co-crystallized ligands (AM11542 and Taranabant) were re-docked into the orthosteric site. The validation was considered successful when the ligand conformation with the lowest energy score closely overlapped with the co-crystallized molecule in the protein X-ray structure, achieving an impressive RMSD < 2.0 Å [64,65].
4.2.3. Protein Rigid Docking
Initially, a protein rigid docking using FDA drugs was performed using 5XRA and 5U09 structures. The dock panel settings included the Triangle Matcher method for placement and the London dG (Equation (1)) as the scoring function [66]. This method takes into account several factors, including rotational and translational entropy (c), energy due to the loss of flexibility of the ligand (Eflex), geometric imperfections of hydrogen bonds (fHB), geometric imperfections of metal ligations (fM) and desolvation energy (Di).
All drugs demonstrating a binding affinity of at least −8.5 kcal/mol were identified in both structures. Subsequently, these molecules underwent further screening through an additional filter stage.
4.2.4. Blind Induced-Fit Docking
A blind docking approach was then conducted, specifically targeting the extracellular domain region up to the orthosteric site within the 5XRA structure. In this phase, an induced-fit docking strategy was employed, allowing for flexibility within the binding site.
An induced-fit refinement was executed, enabling both the ligand and the active site to move freely. The resulting poses were scored using the London dG scoring function.
The graphic phases of the virtual screening are shown in the following Scheme 1.
4.3. Binding Assay
Mouse brain was resuspended in 2 mM Tris–EDTA, 320 mM sucrose, 5 mM MgCl2 (pH 7.4), homogenized in a Potter homogenizer, centrifuged three times at 1000× g (10 min) and the pellet was discharged. The supernatant was further centrifuged at 18,000× g (30 min), and the resulting pellet was then resuspended in HBSS. For the rapid filtration assays, 50 μg of these membrane fractions were used for each test, along with the radiolabeled agonist [3H]CP55,940 at a concentration of 400 pM. The effect of different drugs on CB1R binding was tested by adding the substance, with 80 min of preincubation, directly to the incubation medium followed by incubation time of 40 min at 37 °C [67].
4.4. Activity-Based Protein Profiling (ABPP)
Mouse brain lysate (14.6 μL, 2.0 mg/mL lysate, cytosol or membrane fraction) was pre-incubated with vehicle or drugs (0.375 μL, 30 min, 37 °C) followed by an incubation with the activity-based probe (100 nM of MB064 and 100 nM of Bodipy, 10 min, RT). Reactions were quenched with 5 μL Laemmli buffer + β-mercaptoethanol. The reaction was resolved on a 30% acrylamide SDS-PAGE gel (180 V, 75 min). Gels were scanned using Cy 2 (80 s), Cy3 (120 s) and Cy5 (10 s) multichannel settings and subsequently stained with Coomassie after scanning. Fluorescence was normalized to Coomassie staining and quantified with Image Lab (Bio-Rad) [53,54].
4.5. Cell Viability
The HaCaT cells are an aneuploid immortalized keratinocyte cell line from adult human skin [68]. HaCaT cells were cultured at 37 °C, 5% CO2 in DMEM supplemented with 10% (v/v) fetal bovine serum (Gibco) and penicillin (100 unit/mL)-streptomycin (0.1 mg/mL). Miltefosine was used at 25 µM, while SR1 was at 1 µM. The incubation time of Miltefosine count was 30 min, and after a further 4 h, the percentage of viable cells was assessed. SR1 underwent a preincubation period of 30 min, 1 h, 1 h and 30 min, and 2 h, allowing for a comprehensive evaluation of its effects.
5. Conclusions
This study presents a compelling demonstration of how a combination of computational and experimental approaches can shed light on CB1R as a promising new target for marketed drugs. Moreover, the quest to design selective and safe drugs for CB1R can be addressed through a novel strategy. To this aim, compounds that act on specific targets and only partially interact with CB1R offer a potentially safer therapeutic approach, minimizing adverse effects compared to direct CB1R-targeting molecules. Raloxifene, Raltegravir and Miltefosine could represent an example of this strategy. Considering the widespread distribution of CB1R and its involvement in diverse signaling cascades, further investigations similar to ours become imperative to evaluate its role in the mechanisms of certain drugs.
Acknowledgments
The authors thank Mario van der Stelt and Annemarieke Postmus of the Leiden University for their support in the ABPP experiments.
Institutional Review Board Statement
We used mouse brain in stock. We did not use new animals.
Informed Consent Statement
Not applicable.
Data Availability Statement
Data are contained within the article.
Conflicts of Interest
The authors declare no conflict of interest.
Funding Statement
This investigation was partly supported by the Italian Ministry of University and Research (MUR) under the competitive PRIN 2022 grant n. 20224CPSYL to MM. EC was partly supported by a FEBS short term fellowship 2021.
Footnotes
Footnote Group
References
Untitled section
References
- 1.Li Q., Shah S. Structure-Based Virtual Screening. In: Wu C.H., Arighi C.N., Ross K.E., editors. Protein Bioinformatics. Volume 1558. Springer; New York, NY, USA: 2017. pp. 111–124. Methods in Molecular Biology.
- 2.Gambacorta N., Ciriaco F., Amoroso N., Altomare C.D., Bajorath J., Nicolotti O. CIRCE: Web-Based Platform for the Prediction of Cannabinoid Receptor Ligands Using Explainable Machine Learning. J. Chem. Inf. Model. 2023;63:5916–5926. doi: 10.1021/acs.jcim.3c00914.
- 3.Jumper J., Evans R., Pritzel A., Green T., Figurnov M., Ronneberger O., Tunyasuvunakool K., Bates R., Žídek A., Potapenko A., et al. Highly Accurate Protein Structure Prediction with AlphaFold. Nature. 2021;596:583–589. doi: 10.1038/s41586-021-03819-2.
- 4.Ou-Yang S., Lu J., Kong X., Liang Z., Luo C., Jiang H. Computational Drug Discovery. Acta Pharmacol. Sin. 2012;33:1131–1140. doi: 10.1038/aps.2012.109.
- 5.Genheden S., Ryde U. The MM/PBSA and MM/GBSA Methods to Estimate Ligand-Binding Affinities. Expert Opin. Drug Discov. 2015;10:449–461. doi: 10.1517/17460441.2015.1032936.
- 6.Pushpakom S., Iorio F., Eyers P.A., Escott K.J., Hopper S., Wells A., Doig A., Guilliams T., Latimer J., McNamee C., et al. Drug Repurposing: Progress, Challenges and Recommendations. Nat. Rev. Drug Discov. 2019;18:41–58. doi: 10.1038/nrd.2018.168.
- 7.Stasiulewicz A., Lesniak A., Setny P., Bujalska-Zadrożny M., Sulkowska J.I. Identification of CB1 Ligands among Drugs, Phytochemicals and Natural-Like Compounds: Virtual Screening and In Vitro Verification. ACS Chem. Neurosci. 2022;13:2991–3007. doi: 10.1021/acschemneuro.2c00502.
- 8.Catella-Lawson F., Reilly M.P., Kapoor S.C., Cucchiara A.J., DeMarco S., Tournier B., Vyas S.N., FitzGerald G.A. Cyclooxygenase Inhibitors and the Antiplatelet Effects of Aspirin. N. Engl. J. Med. 2001;345:1809–1817. doi: 10.1056/NEJMoa003199.
- 9.Gordon J.L., Macintyre D.E. Proceedings: Inhibition of Human Platelet Aggregation by Aspirin In Vitro and Ex Vivo. Br. J. Pharmacol. 1974;52:451P.
- 10.Menter D.G., Bresalier R.S. An Aspirin a Day: New Pharmacological Developments and Cancer Chemoprevention. Annu. Rev. Pharmacol. Toxicol. 2023;63:165–186. doi: 10.1146/annurev-pharmtox-052020-023107.
- 11.Joharatnam-Hogan N., Hatem D., Cafferty F.H., Petrucci G., Cameron D.A., Ring A., Kynaston H.G., Gilbert D.C., Wilson R.H., Hubner R.A., et al. Thromboxane Biosynthesis in Cancer Patients and Its Inhibition by Aspirin: A Sub-Study of the Add-Aspirin Trial. Br. J. Cancer. 2023;129:706–720. doi: 10.1038/s41416-023-02310-1.
- 12.Virk H.U.H., Escobar J., Rodriguez M., Bates E.R., Khalid U., Jneid H., Birnbaum Y., Levine G.N., Smith S.C., Krittanawong C. Dual Antiplatelet Therapy: A Concise Review for Clinicians. Life. 2023;13:1580. doi: 10.3390/life13071580.
- 13.Miller M.T. Thalidomide Embryopathy: A Model for the Study of Congenital Incomitant Horizontal Strabismus. Trans. Am. Ophthalmol. Soc. 1991;89:623–674.
- 14.Vargesson N. Thalidomide-Induced Teratogenesis: History and Mechanisms. Birth Defects Res. Part C Embryo Today Rev. 2015;105:140–156. doi: 10.1002/bdrc.21096.
- 15.Fabro S., Smith R.L., Williams R.T. Toxicity and Teratogenicity of Optical Isomers of Thalidomide. Nature. 1967;215:296. doi: 10.1038/215296a0.
- 16.Franks M.E., Macpherson G.R., Figg W.D. Thalidomide. Lancet. 2004;363:1802–1811. doi: 10.1016/S0140-6736(04)16308-3.
- 17.Langtry H.D., Markham A. Sildenafil: A Review of Its Use in Erectile Dysfunction. Drugs. 1999;57:967–989. doi: 10.2165/00003495-199957060-00015.
- 18.Biggioggero M., Crotti C., Becciolini A., Favalli E.G. Tocilizumab in the Treatment of Rheumatoid Arthritis: An Evidence-Based Review and Patient Selection. Drug Des. Devel. Ther. 2018;13:57–70. doi: 10.2147/DDDT.S150580.
- 19.Freitas E., Guttman-Yassky E., Torres T. Baricitinib for the Treatment of Alopecia Areata. Drugs. 2023;83:761–770. doi: 10.1007/s40265-023-01873-w.
- 20.von Hentig N. Repositioning HIV Protease Inhibitors and Nucleos(t)Ide RNA Polymerase Inhibitors for the Treatment of SARS-CoV-2 Infection and COVID-19. Eur. J. Clin. Pharmacol. 2021;77:1297–1307. doi: 10.1007/s00228-021-03108-x.
- 21.Beigel J.H., Tomashek K.M., Dodd L.E., Mehta A.K., Zingman B.S., Kalil A.C., Hohmann E., Chu H.Y., Luetkemeyer A., Kline S., et al. Remdesivir for the Treatment of COVID-19—Final Report. N. Engl. J. Med. 2020;383:1813–1826. doi: 10.1056/NEJMoa2007764.
- 22.Chakraborty C., Sharma A.R., Bhattacharya M., Agoramoorthy G., Lee S.-S. The Drug Repurposing for COVID-19 Clinical Trials Provide Very Effective Therapeutic Combinations: Lessons Learned from Major Clinical Studies. Front. Pharmacol. 2021;12:704205. doi: 10.3389/fphar.2021.704205.
- 23.Howlett A.C. International Union of Pharmacology. XXVII. Classification of Cannabinoid Receptors. Pharmacol. Rev. 2002;54:161–202. doi: 10.1124/pr.54.2.161.
- 24.Maccarrone M., Di Marzo V., Gertsch J., Grether U., Howlett A.C., Hua T., Makriyannis A., Piomelli D., Ueda N., Van Der Stelt M. Goods and Bads of the Endocannabinoid System as a Therapeutic Target: Lessons Learned after 30 Years. Pharmacol. Rev. 2023;75:885–958. doi: 10.1124/pharmrev.122.000600.
- 25.Yin H., Chu A., Li W., Wang B., Shelton F., Otero F., Nguyen D.G., Caldwell J.S., Chen Y.A. Lipid G Protein-Coupled Receptor Ligand Identification Using β-Arrestin PathHunterTM Assay. J. Biol. Chem. 2009;284:12328–12338. doi: 10.1074/jbc.M806516200.
- 26.Hua T., Vemuri K., Pu M., Qu L., Han G.W., Wu Y., Zhao S., Shui W., Li S., Korde A., et al. Crystal Structure of the Human Cannabinoid Receptor CB1. Cell. 2016;167:750–762.e14. doi: 10.1016/j.cell.2016.10.004.
- 27.Shao Z., Yin J., Chapman K., Grzemska M., Clark L., Wang J., Rosenbaum D.M. High-Resolution Crystal Structure of the Human CB1 Cannabinoid Receptor. Nature. 2016;540:602–606. doi: 10.1038/nature20613.
- 28.Hua T., Vemuri K., Nikas S.P., Laprairie R.B., Wu Y., Qu L., Pu M., Korde A., Jiang S., Ho J.-H., et al. Crystal Structures of Agonist-Bound Human Cannabinoid Receptor CB1. Nature. 2017;547:468–471. doi: 10.1038/nature23272.
- 29.Zou S., Kumar U. Cannabinoid Receptors and the Endocannabinoid System: Signaling and Function in the Central Nervous System. Int. J. Mol. Sci. 2018;19:833. doi: 10.3390/ijms19030833.
- 30.Herkenham M., Lynn A., Johnson M., Melvin L., De Costa B., Rice K. Characterization and Localization of Cannabinoid Receptors in Rat Brain: A Quantitative In Vitro Autoradiographic Study. J. Neurosci. 1991;11:563–583. doi: 10.1523/JNEUROSCI.11-02-00563.1991.
- 31.Cavuoto P., McAinch A.J., Hatzinikolas G., Janovská A., Game P., Wittert G.A. The Expression of Receptors for Endocannabinoids in Human and Rodent Skeletal Muscle. Biochem. Biophys. Res. Commun. 2007;364:105–110. doi: 10.1016/j.bbrc.2007.09.099.
- 32.DiPatrizio N.V. Endocannabinoids in the Gut. Cannabis Cannabinoid Res. 2016;1:67–77. doi: 10.1089/can.2016.0001.
- 33.Malenczyk K., Jazurek M., Keimpema E., Silvestri C., Janikiewicz J., Mackie K., Di Marzo V., Redowicz M.J., Harkany T., Dobrzyn A. CB1 Cannabinoid Receptors Couple to Focal Adhesion Kinase to Control Insulin Release. J. Biol. Chem. 2013;288:32685–32699. doi: 10.1074/jbc.M113.478354.
- 34.Osei-Hyiaman D., Liu J., Zhou L., Godlewski G., Harvey-White J., Jeong W., Bátkai S., Marsicano G., Lutz B., Buettner C., et al. Hepatic CB1 Receptor Is Required for Development of Diet-Induced Steatosis, Dyslipidemia, and Insulin and Leptin Resistance in Mice. J. Clin. Investig. 2008;118:3160–3169. doi: 10.1172/JCI34827.
- 35.Devane W.A., Dysarz F.A., Johnson M.R., Melvin L.S., Howlett A.C. Determination and Characterization of a Cannabinoid Receptor in Rat Brain. Mol. Pharmacol. 1988;34:605–613.
- 36.Pertwee R.G., Howlett A.C., Abood M.E., Alexander S.P.H., Di Marzo V., Elphick M.R., Greasley P.J., Hansen H.S., Kunos G., Mackie K., et al. International Union of Basic and Clinical Pharmacology. LXXIX. Cannabinoid Receptors and Their Ligands: Beyond CB1 and CB2. Pharmacol. Rev. 2010;62:588–631. doi: 10.1124/pr.110.003004.
- 37.Fezza F., Bari M., Florio R., Talamonti E., Feole M., Maccarrone M. Endocannabinoids, Related Compounds and Their Metabolic Routes. Molecules. 2014;19:17078–17106. doi: 10.3390/molecules191117078.
- 38.Baggelaar M.P., Maccarrone M., Van Der Stelt M. 2-Arachidonoylglycerol: A Signaling Lipid with Manifold Actions in the Brain. Prog. Lipid Res. 2018;71:1–17. doi: 10.1016/j.plipres.2018.05.002.
- 39.Götz M.R., Collado J.A., Fernández-Ruiz J., Fiebich B.L., García-Toscano L., Gómez-Cañas M., Koch O., Leha A., Muñoz E., Navarrete C., et al. Structure–Effect Relationships of Novel Semi-Synthetic Cannabinoid Derivatives. Front. Pharmacol. 2019;10:1284. doi: 10.3389/fphar.2019.01284.
- 40.Pertwee R.G. Cannabinoid Receptor Ligands: Clinical and Neuropharmacological Considerations, Relevant to Future Drug Discovery and Development. Expert Opin. Investig. Drugs. 2000;9:1553–1571. doi: 10.1517/13543784.9.7.1553.
- 41.Schoeder C.T., Hess C., Madea B., Meiler J., Müller C.E. Pharmacological Evaluation of New Constituents of “Spice”: Synthetic Cannabinoids Based on Indole, Indazole, Benzimidazole and Carbazole Scaffolds. Forensic Toxicol. 2018;36:385–403. doi: 10.1007/s11419-018-0415-z.
- 42.King A. Neuropsychiatric Adverse Effects Signal the End of the Line for Rimonabant. Nat. Rev. Cardiol. 2010;7:602. doi: 10.1038/nrcardio.2010.148.
- 43.Mitchell P.B., Morris M.J. Depression and Anxiety with Rimonabant. Lancet. 2007;370:1671–1672. doi: 10.1016/S0140-6736(07)61705-X.
- 44.Adams A.J., Banister S.D., Irizarry L., Trecki J., Schwartz M., Gerona R. “Zombie” Outbreak Caused by the Synthetic Cannabinoid AMB-FUBINACA in New York. N. Engl. J. Med. 2017;376:235–242. doi: 10.1056/NEJMoa1610300.
- 45.Kumar K.K., Shalev-Benami M., Robertson M.J., Hu H., Banister S.D., Hollingsworth S.A., Latorraca N.R., Kato H.E., Hilger D., Maeda S., et al. Structure of a Signaling Cannabinoid Receptor 1-G Protein Complex. Cell. 2019;176:448–458.e12. doi: 10.1016/j.cell.2018.11.040.
- 46.Manandhar A., Haron M.H., Klein M.L., Elokely K. Understanding the Dynamics of the Structural States of Cannabinoid Receptors and the Role of Different Modulators. Life. 2022;12:2137. doi: 10.3390/life12122137.
- 47.Maguire J.J., Kuc R.E., Davenport A.P. Radioligand Binding Assays and Their Analysis. In: Davenport A.P., editor. Receptor Binding Techniques. Volume 897. Humana Press; Totowa, NJ, USA: 2012. pp. 31–77. Methods in Molecular Biology.
- 48.Irwin J.J., Shoichet B.K. ZINC—A Free Database of Commercially Available Compounds for Virtual Screening. J. Chem. Inf. Model. 2005;45:177–182. doi: 10.1021/ci049714+.
- 49.Jakowiecki J., Orzeł U., Chawananon S., Miszta P., Filipek S. The Hydrophobic Ligands Entry and Exit from the GPCR Binding Site-SMD and SuMD Simulations. Molecules. 2020;25:1930. doi: 10.3390/molecules25081930.
- 50.Hua T., Li X., Wu L., Iliopoulos-Tsoutsouvas C., Wang Y., Wu M., Shen L., Brust C.A., Nikas S.P., Song F., et al. Activation and Signaling Mechanism Revealed by Cannabinoid Receptor-Gi Complex Structures. Cell. 2020;180:655–665.e18. doi: 10.1016/j.cell.2020.01.008.
- 51.Franks L.N., Ford B.M., Prather P.L. Selective Estrogen Receptor Modulators: Cannabinoid Receptor Inverse Agonists with Differential CB1 and CB2 Selectivity. Front. Pharmacol. 2016;7:503. doi: 10.3389/fphar.2016.00503.
- 52.Soli M., Bertaccini A., Carparelli F., Gotti R., Cavrini V., Andrisano V., Martorana G. Vasoactive cocktails for erectile dysfunction: Chemical stability of PGE1, paraverine and phentolamine. J. Urol. 1998;160:551–555. doi: 10.1016/S0022-5347(01)62957-4.
- 53.Barglow K.T., Cravatt B.F. Activity-Based Protein Profiling for the Functional Annotation of Enzymes. Nat. Methods. 2007;4:822–827. doi: 10.1038/nmeth1092.
- 54.Van Rooden E.J., Florea B.I., Deng H., Baggelaar M.P., Van Esbroeck A.C.M., Zhou J., Overkleeft H.S., Van Der Stelt M. Mapping In Vivo Target Interaction Profiles of Covalent Inhibitors Using Chemical Proteomics with Label-Free Quantification. Nat. Protoc. 2018;13:752–767. doi: 10.1038/nprot.2017.159.
- 55.Janssen A.P.A., Van Der Vliet D., Bakker A.T., Jiang M., Grimm S.H., Campiani G., Butini S., Van Der Stelt M. Development of a Multiplexed Activity-Based Protein Profiling Assay to Evaluate Activity of Endocannabinoid Hydrolase Inhibitors. ACS Chem. Biol. 2018;13:2406–2413. doi: 10.1021/acschembio.8b00534.
- 56.Mouscadet J.-F., Tchertanov L. Raltegravir: Molecular Basis of Its Mechanism of Action. Eur. J. Med. Res. 2009;14:5. doi: 10.1186/2047-783X-14-S3-5.
- 57.Scaranti M., Cojocaru E., Banerjee S., Banerji U. Exploiting the Folate Receptor α in Oncology. Nat. Rev. Clin. Oncol. 2020;17:349–359. doi: 10.1038/s41571-020-0339-5.
- 58.Wieder T., Orfanos C.E., Geilen C.C. Induction of Ceramide-Mediated Apoptosis by the Anticancer Phospholipid Analog, Hexadecylphosphocholine. J. Biol. Chem. 1998;273:11025–11031. doi: 10.1074/jbc.273.18.11025.
- 59.Falasca M., Maccarrone M. Cannabinoids and Cancer. Cancers. 2021;13:4458. doi: 10.3390/cancers13174458.
- 60.De Miguel R., Montejano R., Stella-Ascariz N., Arribas J.R. A Safety Evaluation of Raltegravir for the Treatment of HIV. Expert Opin. Drug Saf. 2018;17:217–223. doi: 10.1080/14740338.2018.1411903.
- 61.Ismail M., Hasan H., El-Orfali Y., Ismail H., Khawaja G. Anti-Inflammatory, Antioxidative, and Hepatoprotective Effects of Trans Δ 9-Tetrahydrocannabinol/Sesame Oil on Adjuvant-Induced Arthritis in Rats. Evid. Based Complement. Alternat. Med. 2018;2018:9365464. doi: 10.1155/2018/9365464.
- 62.Sundar S., Jha T.K., Thakur C.P., Engel J., Sindermann H., Fischer C., Junge K., Bryceson A., Berman J. Oral Miltefosine for Indian Visceral Leishmaniasis. N. Engl. J. Med. 2002;347:1739–1746. doi: 10.1056/NEJMoa021556.
- 63.Maccarrone M., Di Rienzo M., Battista N., Gasperi V., Guerrieri P., Rossi A., Finazzi-Agrò A. The Endocannabinoid System in Human Keratinocytes. J. Biol. Chem. 2003;278:33896–33903. doi: 10.1074/jbc.M303994200.
- 64.Davis I.W., Baker D. RosettaLigand Docking with Full Ligand and Receptor Flexibility. J. Mol. Biol. 2009;385:381–392. doi: 10.1016/j.jmb.2008.11.010.
- 65.McConkey B.J., Sobolev V., Edelman M. The Performance of Current Methods in Ligand–Protein Docking. Curr. Sci. 2002;83:845–856.
- 66.Kalinowsky L., Weber J., Balasupramaniam S., Baumann K., Proschak E. A Diverse Benchmark Based on 3D Matched Molecular Pairs for Validating Scoring Functions. ACS Omega. 2018;3:5704–5714. doi: 10.1021/acsomega.7b01194.
- 67.Fezza F., Oddi S., Di Tommaso M., De Simone C., Rapino C., Pasquariello N., Dainese E., Finazzi-Agrò A., Maccarrone M. Characterization of Biotin-Anandamide, a Novel Tool for the Visualization of Anandamide Accumulation. J. Lipid Res. 2008;49:1216–1223. doi: 10.1194/jlr.M700486-JLR200.
- 68.Boukamp P., Petrussevska R.T., Breitkreutz D., Hornung J., Markham A., Fusenig N.E. Normal Keratinization in a Spontaneously Immortalized Aneuploid Human Keratinocyte Cell Line. J. Cell Biol. 1988;106:761–771. doi: 10.1083/jcb.106.3.761.
Associated Data
Data Availability Statement
Data are contained within the article.