A series of novel 2-pyridylbenzimidazole derivatives was rationally designed and synthesized based on our previous studies on benzimidazole 14, a CB1 agonist used as a template for optimization. In the present series, 21 compounds displayed high affinities with Ki values in the nanomolar range. JM-39 (compound 39) was the most active of the series (KiCB1 = 0.53 nM), while compounds 31 and 44 exhibited similar affinities to WIN 55212-2. CoMFA analysis was performed based on the biological data obtained and resulted in a statistically significant CoMFA model with high predictive value (q2 = 0.710, r2 = 0.998, r2pred = 0.823).
A series of novel 2-pyridylbenzimidazole derivatives was rationally designed and synthesized based on our previous studies on benzimidazole 14, a CB1 agonist used as a template for optimization. In the present series, 21 compounds displayed high affinities with Ki values in the nanomolar range. JM-39 (compound 39) was the most active of the series (KiCB1 = 0.53 nM), while compounds 31 and 44 exhibited similar affinities to WIN 55212-2. CoMFA analysis was performed based on the biological data obtained and resulted in a statistically significant CoMFA model with high predictive value (q2 = 0.710, r2 = 0.998, r2pred = 0.823).
The regulation of the endocannabinoid system (ECS) has been studied extensively because of its great potential in the treatment of many conditions. CB1 agonists has been associated to the treatment of pain, [1], cancer [2], and emesis [3] among others. Otherwise, CB1 antagonists have been mainly used in the treatment of obesity [4,5]. Its physiological effects are mediated through the so-called cannabinoid receptors, which are attractive targets for current drug development [6,7]. These proteins belong to the G protein-coupled receptor (GPCR) family and have been the focus of much recent work since the discovery of the CB1 and CB2 receptor subtypes [8,9]. Among the reported and non-classical cannabinergic ligands [10,11,12], several compounds contain a central heterocyclic scaffold such as indole or benzofuran, displaying a range of different pharmacological profiles [13,14]. For example, the aminoalkylindole derivative WIN55212-2 has been reported to be a potent agonist of the CB1 receptor (Ki = 9.4 nM) [15], AM679 is a non-selective CB1 and CB2 agonist (Ki = 13.5 nM CB1) [16], JWH-007 has been reported as a non selective CB1/CB2cannabinoid agonist (Ki = 9.5 nM) [17], while the benzofuran derivative LY320135 is a selective CB1 antagonist (Ki = 131 nM) [18,19] (Figure 1). In the field of cannabimimetic heterocycles, benzimidazoles have been recently reported as CB2 agonist with good CNS penetration by Pfizer [20]. By other hand, Abbott Laboratories have developed azaindoles with agonist activity over the CB1 and CB2 receptors [21], while imidazopyridines from Merck have been reported to be selective CB2 agonists [22] (Figure 1).
Figure 1
Previously and recently reported cannabinoid ligands, prototype 14 (JM-6), and 39 (JM-39), the most active compound of the present series. KiCB1 data for 14 (JM-6) and 39 (JM-39) reported herein.
Previously and recently reported cannabinoid ligands, prototype 14 (JM-6), and 39 (JM-39), the most active compound of the present series. KiCB1 data for 14 (JM-6) and 39 (JM-39) reported herein.For many years our research group has been interested in the synthesis of benzimidazole derivatives as ligands of the type 1 cannabinoid receptor. In recent screening tests using mice brain synaptosomes [23], compound 14 (JM-6) was identified as a high affinity CB1 agonist (Ki = 98.2 nM). Considering the structural resemblance of 14 to the reported ligands (Figure 1), we defined a strategy to develop new and more active CB1 ligand derivatives. Docking simulation studies of 14 with our reported CB1 model [24] allowed us to propose the existence of three major interaction regions: (a) a hydrogen bond between the pyridine ring and the quaternary amino group of the side chain of residue K3.28(192); (b) the naphthalene moiety establishing hydrophobic and aromatic interactions with the side-chain residues F2.57(170), F2.61(174), L3.32(196), I5.31(267), F7.35(379), A7.36(380); and (c) the benzene ring of the benzimidazolecore which displayed hydrophobic and aromatic interactions with the side chain residues F2.57(170), F3.36(200), W6.48(356) and L6.51(359).In order to explore SAR requirements of benzimidazoles and optimize the affinity with the CB1 receptor, compound 14 was subjected to the following changes (Figure 2A): (i) Region I: replacement of the 3-pyridyl ring by 2- and 4-pyridyl isomers to study the influence of the position of the nitrogen atom on hydrogen bond interactions. (ii) Region II and linker: the 1-naphthyl moiety was replaced by a 2-naphthyl group and the carbonyl function exchanged for a methylene group aimed at providing an improved orientation and a more favorable conformation in the receptor cavity. (iii) Region III. Finally, substitution on the benzimidazole ring (C-5 and C-6) with fluorine and methyl groups was thought to increase π-stacking interactions. Accordingly, we report the design, synthesis and biological evaluation of novel series of 1 and 2-naphthyl, and 1 and 2-naphthoyl -2-pyridyl-benzimidazole derivatives (Figure 2B) targeting the CB1 receptor. Structure-activity relationships were explored by molecular docking in an updated CB1 receptor model and a 3D-QSAR (CoMFA) was developed using the docking-based alignment.
Figure 2
A. Regions for modification on compound 14. B. General structure of synthesized compounds.
A. Regions for modification on compound 14. B. General structure of synthesized compounds.
2. Results and Discussion
2.1. Chemistry
The synthetic sequence for obtaining the target compounds is shown in Scheme 1. A series of substituted o-phenylenediamines (benzene-1,2-diamine; 4,5-dimethylbenzene-1,2-diamine and 4-fluorobenzene-1,2-diamine) were reacted with 2-, 3- and 4-pyridinecarbaldehyde in ethanol in the presence of cerium (IV) ammonium nitrate (CAN) and H2O2 as oxidizing agents [25], yielding benzimidazole derivatives 1–9. The obtained benzimidazoles 1–9 were finally acylated [26,27] with 1- or 2-naphthoyl chloride or alkylated with 1- or 2-(bromomethyl)naphthalene to yield compounds 10–56.
Scheme 1
Synthesis of Compounds 1–9 and 10–56.
Synthesis of Compounds 1–9 and 10–56.Reagents and Conditions: (i) Ce(NH4)2(NO3)6 (CAN), H2O2, 50 °C 1 h. (ii) THF, NaH, 1- or 2-naphthoyl chloride, or 1- or 2-(bromomethyl)naphthalene.Since the benzimidazole ring exists in tautomeric forms, the acylation or alkylation of 2-pyridin-5-fluoro-1H-benzo[d]imidazoles generate two regioisomers whose structures have been definitively assigned by HMBC experiments which showed an HMBCcorrelation between C-3a (Figure 3) and the hydrogen at C-7 in the 5-fluoro isomer, and with the hydrogens at C-5 and C-7 in the 6-fluoro isomer.
Figure 3
HMBC correlations in regioisomeric products.
HMBCcorrelations in regioisomeric products.
2.2. CB1 Binding Affinities, Docking and SAR
In our efforts to obtain novel CB1 ligands with improved affinities and to gain experimental evidence concerning the mode of interaction of cannabimimeticbenzimidazoles with this receptor, we have designed and synthesized a series of benzimidazoles 10–56. The CB1 receptor affinity values for these compounds are shown in Table 1.
Table 1
Radioligand displacement results for 2-arylbenzimidazoles 10–56
a.
Comp.
X
Y
Z
R
R1
R2
R3
CB1
(Ki, nM)
10
N
CH
CH
C=O
H
H
1-naphtyl
n.a
11
N
CH
CH
C=O
F
H
1-naphtyl
n.a
12
N
CH
CH
C=O
H
F
1-naphtyl
n.a
13
N
CH
CH
C=O
CH3
CH3
1-naphtyl
n.a
14
CH
N
CH
C=O
H
H
1-naphtyl
98.24 ± 20.0
15
CH
N
CH
C=O
F
H
1-naphtyl
534
16
CH
N
CH
C=O
H
F
1-naphtyl
n.a
17
CH
N
CH
C=O
CH3
CH3
1-naphtyl
n.a
18
CH
CH
N
C=O
H
H
1-naphtyl
37.81 ± 5.3
19
CH
CH
N
C=O
F
H
1-naphtyl
n.a
20
CH
CH
N
C=O
H
F
1-naphtyl
n.a
21
CH
CH
N
C=O
CH3
CH3
1-naphtyl
3.95 ± 0.8
22
N
CH
CH
C=O
H
H
2-naphtyl
n.a
23
N
CH
CH
C=O
H
F
2-naphtyl
11.41 ± 0.9
24
N
CH
CH
C=O
CH3
CH3
2-naphtyl
28.08 ± 2.3
25
CH
N
CH
C=O
H
H
2-naphtyl
n.a
26
CH
N
CH
C=O
F
H
2-naphtyl
5700 ± 200
27
CH
N
CH
C=O
H
F
2-naphtyl
963 ± 155
28
CH
N
CH
C=O
CH3
CH3
2-naphtyl
567 ± 27.8
29
CH
CH
N
C=O
H
H
2-naphtyl
4487 ± 133
30
CH
CH
N
C=O
F
H
2-naphtyl
1028 ± 57.9
31
CH
CH
N
C=O
H
F
2-naphtyl
5.55 ± 0.9
32
CH
CH
N
C=O
CH3
CH3
2-naphtyl
n.a
33
N
CH
CH
CH2
H
H
1-naphtyl
n.a
34
N
CH
CH
CH2
F
H
1-naphtyl
n.a
35
N
CH
CH
CH2
H
F
1-naphtyl
n.a
36
N
CH
CH
CH2
CH3
CH3
1-naphtyl
13.04 ± 1.1
37
CH
N
CH
CH2
H
H
1-naphtyl
75.87 ± 16.6
38
CH
N
CH
CH2
F
H
1-naphtyl
98.18 ± 33.5
39
CH
N
CH
CH2
H
F
1-naphtyl
0.53 ± 0.9
40
CH
N
CH
CH2
CH3
CH3
1-naphtyl
n.a
41
CH
CH
N
CH2
H
H
1-naphtyl
149 ± 1.9
42
CH
CH
N
CH2
F
H
1-naphtyl
n.a
43
CH
CH
N
CH2
H
F
1-naphtyl
614 ± 54.7
44
CH
CH
N
CH2
CH3
CH3
1-naphtyl
4.84 ± 2.0
45
N
CH
CH
CH2
H
H
2-naphtyl
31.82 ± 12.0
46
N
CH
CH
CH2
F
H
2-naphtyl
n.a
47
N
CH
CH
CH2
H
F
2-naphtyl
n.a
48
N
CH
CH
CH2
CH3
CH3
2-naphtyl
47.39 ± 11.1
49
CH
N
CH
CH2
H
H
2-naphtyl
1844 ± 79.8
50
CH
N
CH
CH2
F
H
2-naphtyl
327 ± 60.0
51
CH
N
CH
CH2
H
F
2-naphtyl
10.68 ± 2.8
52
CH
N
CH
CH2
CH3
CH3
2-naphtyl
140 ± 10.1
53
CH
CH
N
CH2
H
H
2-naphtyl
464 ± 30.3
54
CH
CH
N
CH2
F
H
2-naphtyl
44.42 ± 5.5
55
CH
CH
N
CH2
H
F
2-naphtyl
567 ± 7.0
56
CH
CH
N
CH2
CH3
CH3
2-naphtyl
241 ± 10.0
WIN55212
9.40 b
CP55940
0.60 b
Data are means ± SEM of n = 3 separate experiments and are expressed as Ki (nM). n.a. = compounds which did not displace [3H]CP-55,940 by more than 50% at 1 µM in the preliminary screen. b Ki values obtained from technical data sheet of Perkin Elmer.
Significant affinity (nanomolar range) was observed for 21 of the 47 synthesized compounds, while another eight compounds had affinities in the micromolar range. Benzimidazoles 23, 36 and 51 exhibited IC50 values around 20 nM and four molecules (21, 31, 39, 44) exhibited an exceptional affinity with IC50 values under 10 nM. Docking results showed that all the synthesized compounds adopt a Y-like conformation within the CB1 active site, with the pyridine and naphthalenyl or naphthoyl groups as the arms, and the benzimidazole scaffold as the body (Figure 4A,B). The binding site of the ligands is located amidst transmembrane helices 2, 3, 6 and 7. The binding mode of 39, the most active compound of the series, is represented in Figure 4C.
Figure 4
Superimposition of top-scored binding models for the ligands within the CB1 receptor binding site. (A) Side view, and (B) top view from the extracellular side; the loops have been omitted for greater clarity. (C) Schematic representation of the binding mode of the most active compound of the series (39). The 2D representation was obtained using LigandScout v3.03. The main interactions are: (1) Two hydrogen bonds, one between the pyridine nitrogen and K3.28(192) and the other between the benzimidazole N-3 and T3.33(197); (2) Hydrophobic and aromatic interactions of the naphthalene ring with the side-chains of F2.61(174), F7.35(379), A7.36 (380), I5.31(267), between the pyridine ring and L3.29(193), Y5.39(275), M6.55(363), F7.35(379) and of the benzimidazole with F2.57(170), F3.36(200), L6.51(359) and W6.48(356).
Radioligand displacement results for 2-arylbenzimidazoles 10–56
a.Data are means ± SEM of n = 3 separate experiments and are expressed as Ki (nM). n.a. = compounds which did not displace [3H]CP-55,940 by more than 50% at 1 µM in the preliminary screen. b Ki values obtained from technical data sheet of Perkin Elmer.In order to obtain useful structure-activity relationships concerning the mode of interaction of the benzimidazole series, we have focused our analysis on the molecular regions defined in Figure 2A.Superimposition of top-scored binding models for the ligands within the CB1 receptor binding site. (A) Side view, and (B) top view from the extracellular side; the loops have been omitted for greater clarity. (C) Schematic representation of the binding mode of the most active compound of the series (39). The 2D representation was obtained using LigandScout v3.03. The main interactions are: (1) Two hydrogen bonds, one between the pyridinenitrogen and K3.28(192) and the other between the benzimidazole N-3 and T3.33(197); (2) Hydrophobic and aromatic interactions of the naphthalene ring with the side-chains of F2.61(174), F7.35(379), A7.36 (380), I5.31(267), between the pyridine ring and L3.29(193), Y5.39(275), M6.55(363), F7.35(379) and of the benzimidazole with F2.57(170), F3.36(200), L6.51(359) and W6.48(356).
Region I
First, we assessed the effect on affinity of replacing the 3-pyridyl moiety. According to our docking results, this region is crucial for ligand stabilization through an H-bond with the amino group of Lys192. The remainder of the ring presents hydrophobic interactions with the side chains of residues L3.29(193), Y5.39(275) and L6.51(359). Sixteen benzimidazoles were synthesized with a 3-pyridyl group at C-2. Besides the reference compound 14, four derivatives displayed affinities in the nanomolar range (<200 nM, 39, 51, 37, 38, Table 1).As expected, variation of the N position in the pyridinyl substituent to generate the corresponding 2-pyridyl homologs (compounds 35, 47, 33, 34, 10) had a dramatic effect on CB1 affinities (Table 1). The inactivity of these compounds might be explained in terms of an unfavorable distance between Lys 192 and the nitrogen to establish a hydrogen bond.The 4-pyridyl analogs (compounds 43, 55, 41, 42, 56) showed modest affinities compared with the 3-pyridyl derivatives (Table 1). Formation of a weaker H-bond between the nitrogen atom and Lys192compared to the 3-pyridyl compounds seems to be a reasonable explanation. The increased distance between the pyridinenitrogen and Lys192 on going from 3- to 4-pyridyl derivatives observed in our docking models (2 Å and 2.5 Å respectively) corroborates this. However, the activity/inactivity seems actually to be related to the still unknown combined effect of more than one modification. That could explain some exemptions to the observed SAR trends as when homologous compounds 13, 17 and 21 are compared. In summary, 3- or 4-pyridyl seems to be a better substituent than 2-pyridyl in region I.
Region II
According to our docking studies, the naphthyl ring interacts with the CB1 receptor via hydrophobic and aromatic interactions with the side-chains of hydrophobic residues. A comparison of the most active 1-naphthyl derivatives (39, 44, 36, 37, 38 and 41) with their 2-naphthylcongeners (51, 56, 48, 49, 50 and 53) was carried out. In general, the 2-naphthyl derivatives displayed lower CB1 receptor affinities than their 1-naphthylcongeners (Table 1). For benzimidazoles bearing a 2-naphthyl moiety, the naphthalene ring is oriented in a less favourable spatial region in order to interact with the hydrophobic area delimited by Phe379, Phe174, Val196, Ala 380, Ile267 and Phe170.
Linker
The methylene bridge does not present any important interaction with the protein, but ligands containing a keto group as linker can establish, according to our docking analysis, hydrogen bonds with Ser383. However, when the carbonyl group of prototype 14 (Ki = 98.24 nM) was replaced by a methylene linker (compound 37, Ki = 75.87 nM) no significant variation of CB1 receptor affinity was observed. In fact, derivatives with both carbonyl and methylene are among the most active compounds of the series (Table 1).There appears to be no a priori explanation for this experimental evidence. In order to gain some insight into the origin of these unexpected differences in receptor affinity, a comparative molecular modeling and receptor docking study was carried out for the naphthylmethyl and the corresponding naphthoyl congeners. No significant differences were seen in the angular arrangement of the naphthalene ring and the pose of the benzimidazole was maintained. Discarding a different spatial orientation induced by the linker as the explanation for the loss in affinity, we focused our analysis on the penalty imposed by the ligand desolvation factor. As described in the literature [28], desolvation of the ligand plays a significant role in the binding affinity in the cases where the ligand must be entirely desolvated to perform its function. This is of special relevance in membrane proteins because the ligand has to be transferred from the aqueous environment to the transmembrane binding site. Replacement of the carbonyl bridge by a methylene group could enhance the CB1 receptor binding affinities due to a significant decrease in the desolvation penalty, facilitating the entrance of the ligands into the binding site and compensating for the weaker ligand-receptor interaction. This was confirmed by semiempirical AM1-RHF solvation energy calculations (VAMP software) which displayed an increase in the amount of heat required to raise the entropy of carbonyl compounds.
Region III
In order to study the effect of C5 and C6 substitution of the benzimidazole ring on CB1 affinity, the substitution pattern at these positions for the most potent compounds of the series was analyzed. As already mentioned, according to our docking studies the phenyl ring of the benzimidazolecore displayed hydrophobic and aromatic interactions with the side chain of residues Phe170, Trp356 and Leu359. Besides, the basicnitrogen of the benzimidazole nucleus forms a hydrogen bond with the side chain of Thr197 for all the compounds.In general, substitution at the C5 and C6 positions of the phenyl ring of the benzimidazoles with methyl or fluorine seems to be beneficial for CB1 receptor affinity. Indeed, Just two of the most potent derivatives displayed in Table 1 (compounds 45 and 18) are not substituted on this ring. The enhancement of π-stacking interactions between the benzimidazole ring and Phe170 and Trp356 is suggested from our docking analysis. The most potent derivatives are monofluorinated at C6 (39, 31, 51, 23). On the other hand, an important decrease in affinity is seen when the fluorine is at C5, keeping the rest unchanged (38, 30, 50). Our docking analysis revealed a possible solvation sphere of water mediating a H bond between the fluorine atom at C6 fluorine and the TM2 Phe170 backbone, which might be crucial to enhance affinity. The structure-activity relationships for the complete series are summarized in Figure 5.
Figure 5
Structure-activity relationships.
Structure-activity relationships.
2.3. 3D-QSAR (CoMFA)
Active conformation selection is a key step for CoMFA analysis [29]. On the other hand, structural alignment is viewed as one of the most sensitive parameters in CoMFA analysis [30]. The predictive accuracy of CoMFA models and the reliability of the contours are directly dependent on the structural alignment rule [30]. In this case this problem was solved using the conformations obtained by docking studies and active compounds were aligned according to the flexible-ligand docking solutions. As can be seen on Figure 6, there are some molecules that are not well overlayed with respect to many of the other derivatives. Specifically, compounds 21, 36 and 44 display a different docking conformation. Derivatives 21 and 44 are directing the pyridine ring to the same region as the rest of the series compounds. It is noteworthy that compounds 21, 36 and 44 are dimethylated at C5 and C6 positions of the benzimidazole and localize the dimethylated ring towards the same region where the other compounds of the series direct the naphthalene ring. Indeed, this region on the active site favors hydrophobic interactions with bulky substituents. Finally Gasteiger-Hückel charges [31,32] were assigned to all the molecules.
Figure 6
(A) Alignment of the 29 studied molecules. (B) Docking conformation of compounds 21, 36 and 44.
(A) Alignment of the 29 studied molecules. (B) Docking conformation of compounds 21, 36 and 44.A 3D-QSAR model was obtained from CoMFA analysis and its statistical parameters are listed in Table 2. For a reliable predictive model, the square q of the cross-validation coefficient should be greater than 0.5 [33]. This model has high r (0.998), F (944.13), and q (0.71), as well as a small SEE (0.057), suggesting that it is reliable and predictive. The steric and electrostaticcontributions were found to be 35% and 65% respectively.
Table 2
Statistical parameters of the CoMFA model .
Statistic Index Summary of CoMFA
N
q2
r2
r2pred
SEE
SD
Press
F
Contribution
Steric
Electrostatic
7
0.71
0.998
0.823
0.057
5.18
0.92
944.1
35
65
Number of Components Analysis
SEP
0.781
0.707
0.660
0.614
0.619
0.624
0.624
0.649
0.675
0.706
q2
0.376
0.511
0.594
0.667
0.679
0.692
0.710
0.706
0.703
0.698
N
1
2
3
4
5
6
7
8
9
10
N is the optimal number of components, q is the square of the LOO cross-validation (CV) coefficient, r is the square of non cross-validation coefficient, r is the predictive r based only on the test set molecules, SEE is the standard error of estimation of non CV analysis, SEP is the standard error of prediction of LOO analysis, SD is the sum of the squared deviations between the biological activities of molecules in the test set and the mean activity of the training-set molecules, PRESS is the sum of the squared deviations between predicted and actual biological activity values for every molecule in the test set, and F is the F-test value.
Statistical parameters of the CoMFA model .N is the optimal number of components, q is the square of the LOO cross-validation (CV) coefficient, r is the square of non cross-validation coefficient, r is the predictive r based only on the test set molecules, SEE is the standard error of estimation of non CV analysis, SEP is the standard error of prediction of LOO analysis, SD is the sum of the squared deviations between the biological activities of molecules in the test set and the mean activity of the training-set molecules, PRESS is the sum of the squared deviations between predicted and actual biological activity values for every molecule in the test set, and F is the F-test value.The activity values of compounds predicted by CoMFA are listed in Table 3. The plot of the predicted pIC50 values versus the experimental ones for CoMFA analysis is also shown in Figure 7, in which most points are evenly distributed along the line Y = X, suggesting that the quality of the 3D-QSAR model is good. The steric and electrostaticcontour maps of CoMFA are displayed in Figure 8.
Table 3
Experimental and predicted activity for training and test set.
Compd.
Experimental activity (pIC50, M)
Predicted activity (pIC50, M)
Residual
Training Set
39
9.009
9.015
−0.006
21
8.140
8.126
0.013
44
8.052
8.076
−0.024
51
7.708
7.767
−0.059
23
7.679
7.624
0.056
36
7.621
7.622
0.000
24
7.288
7.300
−0.011
45
7.234
7.237
−0.003
54
7.089
7.067
0.022
48
7.061
7.067
−0.006
37
6.857
6.954
−0.097
38
6.745
6.686
0.059
14
6.745
6.693
0.052
41
6.563
6.575
−0.012
56
6.355
6.351
0.004
50
6.222
6.208
0.014
53
6.070
6.009
0.061
15
6.009
6.000
0.008
28
5.983
6.005
−0.022
55
5.983
5.998
−0.015
27
5.753
5.692
0.061
30
5.725
5.726
−0.001
49
5.471
5.429
0.042
26
4.981
5.115
−0.134
Test Set
31
7.992
8.433
−0.441
18
7.159
7.032
0.128
52
6.592
6.777
−0.185
43
5.949
6.703
−0.755
29
5.085
5.409
−0.325
Figure 7
Experimental versus predicted activity for A. Training Set and B. Test Set.
Figure 8
Left: CoMFA steric contour map for the most active compound 39. The green contour indicates where a bulky group favors activity, while the yellow contour indicates where a small group favors activity. Right: CoMFA electrostatic contour map for the most active compound 39. The blue contour indicates where positive charge favors activity, while a red contour indicates where negative charge favors activity.
Experimental and predicted activity for training and test set.Experimental versus predicted activity for A. Training Set and B. Test Set.Left: CoMFA stericcontour map for the most active compound 39. The green contour indicates where a bulky group favors activity, while the yellow contour indicates where a small group favors activity. Right: CoMFA electrostaticcontour map for the most active compound 39. The blue contour indicates where positive charge favors activity, while a red contour indicates where negative charge favors activity.The CoMFA stericcontour map shows two yellow polyhedra near the pyridine ring suggesting that substitution or replacement of this ring with bulky groups is detrimental for activity. An optimal distance for the H-bond with Lys192 is a very important steric restriction in this region. On the other hand, the CoMFA electrostaticcontour map displays red areas above and below the pyridine ring, which means that electronegative atoms are required on this ring for optimal activity. This is consistent with the hydrogen bonding between the pyridinenitrogen and the Lys192 residue observed in our docking analysis. Moreover, next to the red areas two blue areas can be seen. A detailed analysis of the active site reveals the presence of amino acids like tyrosine and threonine, available for H-bonding interactions, beyond Lys192, suggesting that replacement of the pyridine with extended hydrogen-bonding donor groups able to reach these residues might also generate active compounds.
Linker and Region II
For region II the CoMFA stericcontour map shows large green polyhedra around the naphthyl group, suggesting that bulky groups at this position favor affinity. This result is coherent with our docking results, which displayed the naphthyl group oriented toward the lumen of the receptor without any spatial restrictions. With regard to the CoMFA electrostaticcontour map, high electron density rings would increase affinity in the region of the naphthyl ring. Not only surfaces provide valuable information, but empty areas too. That is clear analyzing the corresponding CoMFA electrostaticcontour map for the linker region, which from the standpoint of the CoMFA seems not to be involved in any H-bonding interaction.The CoMFA stericcontour map reveals a small green surface near the fluorine atom at C6, supporting the synthesis of new compounds substituted at that position, for example with chlorine or bromine. The other positions of the benzimidazole moiety (C4 and C5) do not display any steric restrictions. With respect to the CoMFA electrostaticcontour map, the C5 and C6 positions of the benzimidazole ring do not seen to be implicated in relevant electrostatic interactions.
3. Experimental
3.1. General
All organic solvents used for the syntheses were analytical grade. Melting points were determined on a Stuart ScientificSMP3 apparatus and are uncorrected. IR spectra were recorded on a Bruker Vector 22 spectrophotometer using KBr discs. 1H and 13C-NMR spectra were obtained on a Bruker AM-400 spectrometer. The chemical shifts are expressed in ppm (δ scale) downfield from TMS, J values are given in Hertz for solutions in CDCl3 unless otherwise indicated. Column chromatography was performed on Merck silica gel 60 (70–230 mesh). Thin layer chromatographic separations were performed on Merck silica gel 60 (70–230 mesh) chromatofoils. Purity of all final derivatives for biological testing was confirmed to be >95% as determined using elemental analysis carried out on a FISONS EA 1108 CHNS-O analyzer.
General Method for 2-Pyridinyl-1H-benzo[d]imidazoles 1–9
General Method for the Synthesis of (2-Pyridinyl-1H-benzo[d]imidazol-1-yl)(naphthalenyl)methanones 10–32 and 2-Pyridinyl-1-(naphthylmethyl)-1H-benzo[d]imidazoles 33–56
For the CB1 receptor binding assays the new compounds were first subjected to a preliminary screen carried out using three concentrations of test compounds (1, 0.1 and 0.01 µM), humancannabinoid CB1 receptor (PerkinElmer Inc., Waltham, MA, USA, Product No.: 6110129400UA, Lot No.: 509–845-A) and [3H]-(-)-cis-3-[2-hydroxy-4-(1,1-dimethylheptyl)phenyl]-trans-4-(3-hydroxy-propyl)cyclohexanol ([3H]CP-55,940; Kd = 500 pM) as the high affinity ligand as described by the manufacturer (Perkin-Elmer). Compounds which displaced [3H]CP-55,940 by more than 50% at 1 µM were further analyzed. Displacement curves were generated by incubating drugs with 0.5 nM of [3H]CP-55,940. In all cases, Ki values were calculated by applying the Cheng-Prusoff equation [34] to the IC50 values (obtained by GraphPad, GraphPad Software, Inc., La Jolla, CA, USA) for the displacement of the bound radioligand by increasing concentrations of the test compounds (after three separate experiments for each compound).
3.3. Construction of the hCB1 Receptor Homology Model
The sequence of the humanCB1 receptor was retrieved from UniProt database with accession number P21554. The crystal structure of humanadenosine A2A receptor (PDB ID: 3EML) [35], was identified from the Protein Data Bank (PDB) and selected as the template for homology modeling according to the results of a BLAST search. The CB1R model was generated by comparative modeling using MODELLER v9.10 [36,37]. The selection of the template was reinforced with evidence from molecular phylogenetic reports that indicate that this receptor is more evolutionarily related to CB1 than the other members of family A GPCRs for which a crystal structure exists. The Clustal obtained alignment was manually modified to preserve functionally important residues and the disulfide bridge between residues C257 and C264. The residues corresponding to segment 306–324 of intracellular loop 3 were not considered in the modeling. The final alignment obtained between both sequences presents 25.4% and 49.2% of sequence identity and similarity, respectively. Most conserved residues within each helix are denoted using the Ballesteros-Weinstein notation [38]. A set of 100 models was generated, and the best model according to the internal scoring function of the program (PDF score) was subjected to an energy minimization protocol in order to relax the structure and optimize bond geometry. Model validation was performed using the tools available at the NIH SAVES server (http://nihserver.mbi.ucla.edu/SAVES/).The model presents an alpha-carbon root-mean square deviation of 0.64 Å against the template structure (Figure 9A). The final model contains residues 113–420, with a disulfide bridge between C257 and C264, and presents more than 97% of the residues in allowed regions according to Ramachandran plot analysis [39] (Figure 9B).
Figure 9
(A) Superimposition of the obtained CB1R model (red) with the A2AR template (blue). The plasma membrane is delimited by the red and blue dotted lines (B) Ramachandran plot for the CB1R model.
3.4. Docking of the Library of Compounds
Molecular docking of the compounds was performed using FRED v2.2.5 (OpenEye Scientific Software, Santa Fe, NM, USA). The CB1 model binding site was prepared using the FRED receptor, considering the residues involved in the interaction of CB1 with WIN55212, as our group has reported previously [24]. A multiconformer library was prepared using OMEGA v2.4.3 (OpenEye Scientific Software). Candidate poses of the ligands within the receptor site (100) were obtained and optimized using the Chemgauss3 scoring function. Consensus structures of the poses returned from exhaustive docking and optimization were obtained by consensus scoring using the PLP, Chemscore and Chemgauss3 scoring functions. Finally, the top ranked binding modes for each compound were minimized using the CHARMM22 force-field in Discovery Studio v2.1 (Accelrys Inc., San Diego, CA, USA). The minimization protocol allows movement of the side chains of residues positioning within a radius of 6 Å from the mass centroid of all docked ligands, using the conjugate gradient algorithm to a convergence criterion of 0.001 kcal/mol Å for the RMS of energy gradient. Energy evaluation was performed for each obtained complex using the PLP, LigScore, PMF, and LUDI scoring functions. Additional binding energy calculations were performed using the PBSA (Poisson Boltzmann with non-polar Surface Area) implicit solvent model [40], with a dielectricconstant of 4 in order to simulate the membrane environment for the complexes.(A) Superimposition of the obtained CB1R model (red) with the A2AR template (blue). The plasma membrane is delimited by the red and blue dotted lines (B) Ramachandran plot for the CB1R model.
3.5. Generation of CoMFA and Partial Least Squares (PLS) Analysis
CoMFA studies were performed using SYBYL-X 1.2 molecular modeling software (Tripos Inc., St. Louis, MO, USA). Models of steric and electrostatic fields for CoMFA were based on both Lennard–Jones and Coulombic potentials [41]. Steric and electrostatic energies were calculated using an sp3 carbon atom with a van der Waals radius of 0.152 nm, +1 charge, and 0.2 nm grid spacing. The CoMFA cutoff values were set to 125.46 kJ/mol for both steric and electrostatic fields.PLS analysis was used to construct a linear correlation between the CoMFA descriptors (independent variables) and the affinity values (dependent variables) [42]. To select the best model, the cross-validation analysis was performed using the LOO method (and SAMPLS), which deduces the square of the cross-validation coefficient (q) and the optimal number of components N. The non-cross-validation was performed with a column filter value of 2.0 to speed up the analysis and reduce the noise. The predictive r (r) [43,44] was based only on molecules not included in the training set and was obtained according to the following equation: r = (SD − PRESS)/SD (Table 2).
4. Conclusions
In summary, we have discovered a series of novel and potent CB1 receptor ligands. The benzimidazole derivatives were rationally designed based on chemical modifications of the prototype compound 14. A major focus of the optimization effort was to increase CB1 affinity.Biological characterization of these compounds in radioligand binding assays established that the position of the nitrogen atom on the pyridine ring is crucial for activity. Both 3- and 4-pyridyl groups yield active compounds, but 3-pyridyl is a much better substituent than 4-pyridyl and the 2-pyridyl seems useless. These results suggest possible hydrogen bonding with K3.28(192). No significant differences in affinity were observed when the carbonyl linker was replaced by a methylene group.Changing the position of the naphthyl ring from α to β resulted in a decrease of CB1 receptor affinities. Only two of the most potent derivatives are not substituted on the benzimidazole ring, suggesting that substitution at the C5 and/or C6 positions is beneficial for CB1 affinity, probably by enhancing π-stacking interactions in which F2.57(170), F3.36(200) and W6.48(356) may be involved.The QSAR model gave good statistical results in terms of q and r values. Our CoMFA analysis based on binding affinity data of benzimidazole ligands for the CB1 receptor allowed us to approach optimal substitution in the regions labeled I, II and III. The analysis of the contour maps indicates that the activity of the compounds is due largely to their electronic properties rather than to steric or lipophiliccontributions, defining two key areas for further chemical modification: the naphthalene and the pyridine rings. The information gained in this work should lead to the development of additional high-affinity CB1 receptor ligands.
Authors: Antoni R Blaazer; Jos H M Lange; Martina A W van der Neut; Arie Mulder; Femke S den Boon; Taco R Werkman; Chris G Kruse; Wytse J Wadman Journal: Eur J Med Chem Date: 2011-08-23 Impact factor: 6.514
Authors: Christine Watson; Dafydd R Owen; Denise Harding; Kana Kon-I; Mark L Lewis; Helen J Mason; Miyako Matsumizu; Takasuke Mukaiyama; Margarita Rodriguez-Lens; Akiko Shima; Mifune Takeuchi; Isabelle Tran; Tim Young Journal: Bioorg Med Chem Lett Date: 2011-05-25 Impact factor: 2.823
Authors: C C Felder; K E Joyce; E M Briley; M Glass; K P Mackie; K J Fahey; G J Cullinan; D C Hunden; D W Johnson; M O Chaney; G A Koppel; M Brownstein Journal: J Pharmacol Exp Ther Date: 1998-01 Impact factor: 4.030
Authors: Veli-Pekka Jaakola; Mark T Griffith; Michael A Hanson; Vadim Cherezov; Ellen Y T Chien; J Robert Lane; Adriaan P Ijzerman; Raymond C Stevens Journal: Science Date: 2008-10-02 Impact factor: 47.728
Authors: Gastón Apablaza; Luisa Montoya; Cesar Morales-Verdejo; Marco Mellado; Mauricio Cuellar; Carlos F Lagos; Jorge Soto-Delgado; Hery Chung; Carlos David Pessoa-Mahana; Jaime Mella Journal: Molecules Date: 2017-03-05 Impact factor: 4.411