X.-Y. Li et al. / Bioorg. Med. Chem. Lett. 23 (2013) 4617–4621
4619
Leu346
Lue410
Ala350
Thr347
Met342
Met343
Phe404
Lue403
Lue402
Lue384
Lue525
Gly415
Met388
A
Ala350 Leu346
Leu346
Met342
Met342
Thr347
Ala350
Phe404
Thr347
Figure 1. Linear relation comparison between Gscores of QPLD simulation and
Phe404
Met343
Met343
log(IC50) of the seven compounds.
Lue403
Lue410
Lue403
Lue410
Lue402
Lue384
Lue525
Lue402
Lue384
Lue525
Gscore values is completely the same as that of the IC50 values of
the seven compounds against the ERa in vitro: compound 7 > 6 >
Gly415
Gly415C
B
5 > 4 > 3 > 2 >1 (Gscores are negative and a low value indicates
a good result).
Leu346
Leu346
Ala350
Phe404
Ala350
Phe404
Lue403
Met342
Met343
Met342
Met343
Thr347
Lue410
Thr347
Lue410
In order to reveal the relationship between docking calculation
and antiproliferative activity, we conducted a diagram based on
the data of Gscores of QPLD simulation32 and IC50 values (Fig. 1).
For the sake of reflecting the relationship more clearly in the
pattern, here we adopted the logarithms of the IC50 values
(log(IC50) = log10 (IC50)). This figure shows that the results of recep-
tor–ligand binding interactions and antiproliferative activity are
presented into a good linear relationship. The linear relationship
can be written as a mathematical equation IC50 = b0 + b1*Gscore,
where b0 = 6.73, b1 = 0.49 with correlation R2 = 0.95.
Lue402
Lue384
Lue525
Lue402
Lue384
Lue525
Gly415
Met522
Thr347
Leu346
E
D
Ala405
Glu353
Leu346
Ala350
Thr347
Ala350
Met342
Met343
Met342
Met343
Phe404
Phe404
Lue403
Lue410
Lue410
Lue402
Lue384
We further analyzed the interactions between the seven com-
Lue525
Gly415
Met522
Lue402
Lue384
Lue525
pounds and ER
pounds 1–7 docked into the binding sites to ER
a
. Figure 2 showed the simulation results of com-
. The main
Met388
a
Gly415 F
G
binding modes in these complexes can be described as following.
One major interaction is hydrogen bond. Compound 1 formed
hydrogen bonds with the amino acid residues 5-OH/Leu346, N–
Figure 2. Schematic 3-dimensional diagrams of compound 1–7 in ER
Hydrogen bonds are indicated by yellow dotted lines between ligand and protein.
a
active site.
Binding site of ligands to Er . (A) compound 1. (B) compound 2. (C) compound 3 (D)
a
H/Thr347 of ER
bond with 5-OH/Leu346 of ER
a
(Fig. 2A). Compounds 2 and 3 formed a hydrogen
(Fig. 2B and C), respectively. Com-
compound 4 (E) compound 5 (F) compound 6 (G) compound 7. (The yellow dashed
lines show the formation of the H-bonds. Active site residues are represented by
wires and colored by element. Atoms in the ligand are colored as follows: O, red; C,
brown; H, white; N, blue; F, bright green; Cl, dark green.)
a
pound 6 formed a hydrogen bond with 5-OH/Glu353 (Fig. 2F) and
compound 7 form a hydrogen bond with 5-OH/His524 (Fig. 2G).
Another major interaction is p–p interaction. We got the 2-dimen-
sional diagrams from the ligand interaction diagram module in
Schrodinger Suite 2010. From the 2-dimensional diagrams of li-
gand–receptor interactions (Fig. 3), we found that the benzene
rings of Phe 404, Trp383, His524 were able to coordinate the li-
der Waals interations (a*vdW = À2.320 kcal/mol) and hydrogen
bond interactions (Hbond = À1.297 kcal/mol) compared to the
other six compounds (Fig. 4), and thus contributed to its best QM
calculation score (Gscore = À10.138 kcal/mol). Though compound
gands position in the active site of ER
intercations. Compound 1 formed
tween the two rings of benzofuran skeleton and Phe404 (Fig. 3A).
Compounds 2 and 6 formed conjugate interactions between
a
throught Van der Waals
2 formed p–p conjugate interactions with Phe404 similar to com-
p–p conjugate interactions be-
pound 1, but lacked the hydrogen bond with N–H (Fig. 3A and B).
Therefore its hydrogen bond interaction (Hbond = À0.698 kcal/
mol) was much weaker than compound 1 (Hbond = À1.297 kcal/
mol). Similarly, its coulomb force with protein (b*Coul = À0.236
p–p
the two rings of benzofuran skeleton and Phe404, the ring of side
chain and Trp383, respectively (Fig. 3B and F). Compound 3 formed
kcal/mol) was also weaker than that of compound
1
p–p
conjugate interactions between the six-membered ring of
benzofuran skeleton and Phe404, the ring of side chain and
Trp383, respectively (Fig. 3C). Compounds 4 and 5 formed
(b*Coul = À0.538 kcal/mol). Compounds 2 and 3 both formed
hydrogen bonds with Leu346, but the hydrogen bond interaction
p–p
of compound
(Hbond = À0.480 kcal/mol). They also shaped
2
was stronger than that of the latter
conjugate inter-
conjugate interactions between the five-membered ring of benzo-
furan skeleton and His524, the ring of side chain and Phe404,
respectively (Fig. 3D and E).
p–p
actions with Phe404 and Trp383, but Phe404 was merely affected
the six-membered ring of benzofuran skeleton for compound 3
(Fig. 3C), thus brought about a weaker Van der Waals interaction
(a*vdW = À1.367 kcal/mol) compared to compound 2 (a*vdW =
À1.553 kcal/mol). This leads to a lower Gscore for compound 2
Figure 4 displayed the results of QPLD. (GScore = a*vdW +
b*Coul + Lipo + Hbond + Metal + Buryp + RotB + Site,
a = 0.065,
b = 0.130, here we only discussed the main items which are vdW,
Coul and Hbond). Compound 1 possesses the lowest score of Van
(Gscore = À9.512 kcal/mol) than that for compound
3