Ghadamali Khodarahmi*a,
Parvin Asadi*a,
Hossein Farrokhpour*b,
Farshid Hassanzadeha and
Mohammad Dinarib
aDepartment of Medicinal Chemistry, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, 81746-73461, I. R. Iran. E-mail: khodarahmi@pharm.mui.ac.ir; asadi@pharm.mui.ac.ir
bDepartment of Chemistry, Isfahan University of Technology, Isfahan, 84156-83111, I. R. Iran. E-mail: h-farrokh@cc.iut.ac.ir
First published on 22nd June 2015
Considering the potent cytotoxic activities of hybrid benzofuran–imidazolium and quinazolinone derivatives on the breast cancer cell line MCF-7, novel hybrid derivatives incorporating benzofuran, imidazole and quinazolinone pharmacophores were designed using a molecular hybridization approach. Since aromatase is highly expressed in the MCF-7 cell line, we tried to put these pharmacophores together in such a way that they would arrange themselves in a symmetrical shape, similar to aromatase inhibitors. Subsequently, the binding of these novel hybrid compounds to aromatase have been investigated using a docking procedure applying a combined quantum mechanical/molecular mechanical (QM/MM) method. The QM/MM calculation was performed on the reference structures to obtain atomic charges on the ligand atoms. The results indicated that the hybrid compounds were adopted properly within the aromatase binding site, suggesting that they could be potential inhibitors of aromatase. These novel designed compounds engage in hydrophobic and H-bond interactions with the aromatase binding site, which are in agreement with the basic physicochemical features of known aromatase inhibitors. To obtain more accurate results for the binding energies of the ligands, the structures of the ligands with the best interaction energies, obtained from the docking study, were re-optimized using a three-layer ONIOM method (QM:QM:MM) in which the binding pocket of the enzyme was considered as medium-level. The results demonstrated that when the optimized geometrical structures were subjected to re-docking, a better interaction energy was obtained, which strengthens the ability of these compounds to act as potential inhibitors of aromatase.
One of the more computationally inexpensive methods for predicting if and how a ligand will bind to a protein binding pocket, followed by an estimation of the strength of the ligand binding affinity, is molecular docking. This aims to achieve a relative orientation of the protein and ligand such that the free energy of the overall system is minimized. Given the biological and pharmaceutical significance of molecular docking, considerable efforts have been directed towards promoting use of this method.4–6 In the earliest docking approaches, both the ligand and the protein were treated as rigid bodies, while later semi-flexible docking was used, in which the ligand is treated flexibly by allowing bonds to rotate.7,8 Molecular docking simulation of proteins is a fast and inexpensive method for description of the ligand–protein interactions, but poses some difficulties.5 Two of the most important limitations of conventional docking are the assumption of no protein flexibility upon ligand binding and the use of a force field based fixed dielectric charges for both protein and ligand atoms and therefore having false positives and/or false negatives in the energetic quantification of protein–ligand binding.9 For the first problem, it should be mentioned that flexible protein docking methods, which treat the protein in a flexible manner, require high computational costs.10 About the latter, it has been shown that assuming fixed dielectric charges for both protein and ligand atoms leads to low accuracy in protein–ligand docking results.11 Therefore, to increase the accuracy in the docking result, it is reasonable to expend additional effort to improve the quality of the charge. The comprehensive study of polarization and charge transfer requires quantum mechanical methods but their use in biological models requires sophisticated computing systems.12,13 To limit the computational complexity, combined quantum mechanical and molecular mechanical (QM/MM) methods have been developed as an economical approach, in which a small portion of the protein–ligand system is treated in QM detail.14–18 It has been shown that by using an ab initio quantum chemical approach via QM/MM methods, a better assumption of the ligand charges, which take polarization into account, is obtained. Therefore this method could be used as a promising start towards the development of more accurate docking methods for lead optimization applications.11
Molecular hybridization as a rational approach for drug design has attracted much attention of researchers to discover new chemical entities with a potential to afford some promising drugs of the future.19,20 In this method, active compounds and/or pharmacophoric units which recognize and derive from known bioactive molecules are fused to each other directly or with spacer. This method has been employed to develop new anticancer, anti-alzheimer, and antimalarial agents.21–23
Imidazole, as one of the most important pharmacophores in medical chemistry, has a critical role especially in antifungal, anti-bacterial, anti-cancer and sedative drugs.24 Its biocompatibility provides a scaffold for the preparation of different derivatives to afford new bioactive compounds.24 It has also been reported that imidazolium salts show potent cytotoxic activities towards different cancerous cell lines.25–27 For example, benzofuran–imidazolium hybrid compounds (1, 2) have good cytotoxic activity towards the MCF-7 cell line (Scheme 1).26,27 Quinazolinone is another pharmacophore which has been explored for developing pharmaceutically important molecules. Its derivatives have drawn considerable attention due to their profound chemotherapeutic properties including anticancer, antiinflammatory, anticonvulsant, and antidiuretic activities.28,29
Scheme 1 Design of novel analogues based on hybridization of benzofuran, imidazole and quinazolinone pharmacophores. |
Cytochrome P450 enzymes, as a family of isozymes containing heme cofactors, are responsible for many critical enzymatic reactions in living systems. They are, in general, the terminal oxidase enzymes in electron transfer chains.30–33 Cytochrome P450 19A1, commonly known as aromatase, is an enzyme of the cytochrome P450 superfamily that catalyzes the final and rate-limiting step of the conversion of androgens, testosterone and androstenedione, into estrogens, estradiol and estrone, respectively. It is comprised of a polypeptide chain of 503 amino-acid residues and a prosthetic heme group at its active site. An androgen-specific cleft consisting of hydrophobic and polar residues is situated within the confines of the aromatase binding site. This cleft is specific for androstenedione binding to catalyze conversion of androgen to estrogen via a three-step process. Each step requires one mol of O2, one mol of NADPH and NADPH cytochrome reductase. This reaction converts androstenedione, testosterone and 16α-hydroxytestosterone to estrone, 17β-estradiol and 17β, 16α-estriol, respectively. The two initial steps are the typical C19-methyl hydroxylation, while aromatization of the steroid A-ring is catalyzed in the final step.34 To block estrogen production, it is necessary to inhibit the enzyme through the use of aromatase inhibitors.35
According to previous studies, benzofuran–imidazole analogs, as well as quinazolinone derivatives (2), have good cytotoxic effects on the MCF-7 cell line (Scheme 1).26,27,36,37 Therefore, combining quinazolinone and benzofuran–imidazole moieties together as a potential cytotoxic agent towards this cell line would be logical. On the other hand, benzofuran has been used in the structure of some potent aromatase inhibitors.38–40 For example, Whomsley et al.38 identified substituted l-[(benzofuran-2-yl)-phenylmethyl]-imidazoles (3) as a class of potent aromatase inhibitors with in vitro IC50 values <10 nM, which is 80–1000 times the inhibitory activity of aminoglutethimide (Scheme 1). Since aromatase is overexpressed in the MCF-7 cell line, we want to rationally design a novel structure that incorporates these moieties into a single molecular scaffold which could act as an aromatase inhibitor. According to a previous study,41 besides the physicochemical properties of the inhibitor, molecular shape is expected to be extremely important for the access and fit within the active site of the aromatase. Most aromatase inhibitors used to build the common features model have a similar shape that is expected to be complementary to the volume of the aromatase active site. For example letrozole (Scheme 1), a quite rigid third generation inhibitor of aromatase, has a high degree of symmetry. To obtain a final symmetrical shape, which consists of three pharmacophores, the best arrangement is to put each of the pharmacophores on the vertices of a triangle. As the basic physicochemical features of known aromatase inhibitors are a high degree of hydrophobicity and the potential to establish hydrogen bonds,42,43 it is expected that this structure with heterocyclic rings will provide the desired hydrophobicity. On the other hand, hetero atoms in the ring and substitutions located on the structure provide the hydrogen bond potential for this structure.
After designing them, the energetic and structural properties of these compounds were investigated as aromatase inhibitors using molecular docking and QM/MM calculations. Due to the positive charge on the imidazole ring of the designed ligands and also the presence of heme iron in the active site of the enzyme, QM/MM calculations were performed, in which the partial charges of the ligand were re-fitted according to the polarized active site environment of the enzyme to increase the accuracy of the docking results. Then, the structures with the best interaction energy obtained from the docking study were optimized using a three-layer ONIOM method and finally re-docked to the 3D structure of enzyme to obtain the interaction energy. The results demonstrated that these novel designed compounds have good interaction energy and could be used as potential aromatase inhibitors.
To further increase the accuracy of the docking results, three layer ONIOM (QM:QM:MM) calculations were performed on the structures with the best interaction energies obtained from the previous docking study. For this purpose, the geometry of the ligand along with the porphyrin structure of the heme-iron and several residues, including Met 374, Val 373, Val 370, Ile 305, Ala 306, Ile 133, Trp 224, Leu 372, Leu 477, Phe 134 and Thr 310, were optimized in the field of the protein. The density functional theory (DFT) method, employing the B3LYP/LANL2MB/6-31G(d) basis set, was used for the high-level part of system (ligand). The PM6 semiempirical method was used for the medium-level part of system (heme-iron and the residues) and UFF was employed for the rest of protein.
The prediction of the binding affinity in a molecular docking tool is estimated using a scoring function, which generally needs to be both fast and accurate. The electronic interaction is one of the important components of the energy model. So, assuming fixed dielectric charge for the protein and ligand atoms, when considered in the docking procedure, leads to low accuracy of the docking results. This problem is more important for proteins with a metal ion in their active site. In fact, the presence of a metal ion induces a higher polarization effect and enhances restriction of docking in the prediction of electronic interactions. In this study, owing to the positive charge on the imidazole ring of the designed compounds and also the presence of the heme iron in the active site, polarization effects are more important in the energy calculations. So, to enhance the accuracy of the results, improvement of the charge model should be considered. QM/MM methods can help to offer a superior estimate of the electronic interactions. Previous studies have shown that a docking program gives better results if the ligand partial charges are refitted with QM/MM.11 In the present study, we performed the QM/MM calculation on the reference structures to obtain atomic charges on the ligand atoms. These calculated charges are presumed to represent with reasonable accuracy the ligand atom charges that are polarized by the surrounding atoms of the receptor molecule within the binding sites. Then, the ligands with improved charges were employed for re-docking into the aromatase enzyme. Since the top-scoring pose in docking is dependent on the charges of the ligand atoms, these two steps were repeated until the change in the charge values became insignificant. AutoDock4 was used to dock the new hybrid compounds into the aromatase enzyme and record the top scoring structure. Next, QM/MM calculation was performed in which the ligand was treated with PM6 calculation as the QM level, and new atomic charges on the ligand atoms were obtained via Mulliken population analysis. Once the charge values in the ligand files were substituted with these new charge values, another AutoDock4 run was performed and the top scoring structure was recorded again. After repeating these operations three times, the charges were almost constant and were used for the final docking process. Subsequently, with the selection of the best complex between the ligand and protein according to its cluster and binding energy, important interactions were investigated. The free energies of binding (ΔGb) and inhibition constants (Ki) as calculated using AutoDock are summarized in Table 1.
No. | X | Y | Z | E | Cluster | Binding energya | Binding energyb | Binding energyc | Ki |
---|---|---|---|---|---|---|---|---|---|
a The first binding energy calculated with AutoDock.b The binding energy calculated after refitting charge with the values obtained from QM/MM calculation.c The binding energy calculated after fixed change values. | |||||||||
1 | H | H | Methyl | H | 61 | −7.77 | −7.83 | −7.81 | 2.00 μM |
2 | Cl | H | Methyl | H | 60 | −7.83 | −7.97 | −7.95 | 1.50 μM |
3 | Cl | Cl | Methyl | H | 49 | −8.21 | −8.35 | −8.34 | 632.00 nM |
4 | OH | H | Methyl | H | 52 | −7.77 | −8.80 | −8.82 | 409.00 nM |
5 | OH | Methyl | Methyl | H | 50 | −8.08 | −8.32 | −8.35 | 610.00 nM |
6 | OH | OH | Methyl | H | 73 | −9.14 | −9.32 | −9.34 | 68.00 nM |
7 | Methoxy | H | Methyl | H | 46 | −8.14 | −8.30 | −8.33 | 640.00 nM |
8 | Methoxy | Br | Methyl | H | 40 | −7.61 | −7.74 | −7.73 | 2.29 μM |
9 | Methoxy | H | Propyl | H | 32 | −6.27 | −6.41 | −6.43 | 17.50 μM |
10 | Br | H | Methyl | H | 34 | −7.29 | −7.39 | −7.41 | 4.25 μM |
11 | OH | H | Methyl | H | 74 | −8.18 | −8.38 | −8.37 | 1.07 μM |
12 | OH | Cl | Propyl | H | 20 | −7.77 | −8.01 | −8.07 | 1.20 μM |
13 | Br | OH | Propyl | H | 12 | −6.84 | −6.54 | −6.57 | 15.00 μM |
14 | Methoxy | OH | Propyl | H | 16 | −6.48 | −6.37 | −6.34 | 22.00 μM |
15 | OH | Methoxy | Methyl | H | 38 | −7.92 | −8.36 | −8.33 | 639.00 μM |
16 | OH | Cl | Methyl | H | 54 | −8.71 | −8.80 | −8.86 | 411.00 nM |
17 | OH | Br | Methyl | H | 48 | −8.76 | −8.86 | −8.85 | 423.00 nM |
18 | OH | Methoxy | Propyl | H | 12 | −6.80 | −6.63 | −6.66 | 17.03 μM |
19 | OH | OH | Propyl | H | 36 | −7.12 | −7.28 | −7.21 | 4.51 μM |
20 | OH | H | Propyl | H | 24 | −7.54 | −7.09 | −7.12 | 6.04 μM |
21 | OH | Methyl | Propyl | H | 8 | −6.19 | −6.54 | −6.56 | 15.35 μM |
22 | OH | Cl | Methyl | Thiophene | 20 | −6.13 | −5.80 | −5.83 | 48.81 μM |
23 | OH | Br | Methyl | Thiophene | 20 | −6.56 | −6.08 | −6.10 | 25.00 μM |
24 | OH | Methoxy | Methyl | Thiophene | 36 | −6.27 | −5.98 | −6.01 | 39.00 μM |
25 | OH | OH | Methyl | Thiophene | 42 | −6.34 | −7.01 | −7.08 | 5.72 μM |
26 | OH | H | Methyl | Thiophene | 66 | −5.65 | −5.40 | −5.44 | 106.00 μM |
27 | OH | Methyl | Methyl | Thiophene | 32 | −5.72 | −5.32 | −5.35 | 110.00 μM |
28 | OH | Cl | Propyl | Thiophene | 14 | −5.41 | −5.31 | −5.28 | 113.00 μM |
29 | OH | Br | Propyl | Thiophene | 10 | −5.74 | −5.31 | −5.38 | 112.00 μM |
30 | OH | Methoxy | Propyl | Thiophene | 8 | −5.62 | −5.46 | −5.50 | 213.00 μM |
31 | Methoxy | H | Methyl | H | 36 | −8.09 | −8.31 | −8.34 | 825.00 nM |
32 | Methoxy | Methyl | Methyl | H | 50 | −7.8 | −8.01 | −8.08 | 1.18 μM |
33 | Methoxy | OH | Methyl | H | 42 | −7.92 | −7.97 | −8.02 | 1.19 μM |
34 | Methoxy | Methoxy | Methyl | H | 40 | −8.15 | −8.55 | −8.56 | 530.00 nM |
35 | Methoxy | Cl | Methyl | H | 58 | −7.61 | −7.98 | −8.03 | 1.19 μM |
36 | Methoxy | Br | Methyl | H | 54 | −7.82 | −8.02 | −8.13 | 1.04 μM |
37 | Methoxy | Methoxy | Propyl | H | 30 | −6.68 | −6.34 | −6.33 | 22.00 μM |
38 | Methoxy | OH | Propyl | H | 79 | −6.74 | −6.8 | −6.79 | 11.40 μM |
39 | Methoxy | H | Propyl | H | 87 | −6.17 | −6.30 | −6.34 | 9.26 μM |
40 | Methoxy | Methyl | Propyl | H | 10 | −6.03 | −5.83 | −5.86 | 31.86 μM |
41 | Methoxy | Cl | Methyl | Thiophene | 41 | −6.01 | −5.97 | −5.96 | 55.74 nM |
42 | Methoxy | Br | Methyl | Thiophene | 20 | −5.78 | −5.60 | −5.59 | 70.00 μM |
43 | Methoxy | Methoxy | Methyl | Thiophene | 50 | −6.04 | −5.93 | −5.91 | 44.60 μM |
44 | Methoxy | OH | Methyl | Thiophene | 36 | −6.19 | −6.09 | −6.06 | 36.33 μM |
45 | Methoxy | H | Methyl | Thiophene | 50 | −6.01 | −5.97 | −5.98 | 44.45 μM |
46 | Cl | H | Methyl | H | 30 | −7.73 | −7.88 | −7.85 | 1.61 μM |
47 | Cl | Methyl | Methyl | H | 60 | −8.15 | −8.34 | −8.36 | 825.00 nM |
48 | Cl | OH | Methyl | H | 54 | −8.30 | −8.44 | −8.45 | 642.00 nM |
49 | Cl | Methoxy | Methyl | H | 54 | −8.00 | −8.65 | −8.59 | 531.00 nM |
50 | Cl | Cl | Methyl | H | 62 | −7.78 | −8.02 | −8.04 | 1.21 μM |
51 | Cl | Br | Methyl | H | 58 | −7.64 | −7.81 | −7.79 | 1.69 μM |
52 | Cl | Methoxy | Propyl | H | 14 | −5.59 | −6.01 | −5.93 | 44.58 μM |
53 | Cl | OH | Propyl | H | 12 | −6.02 | −6.21 | −6.17 | 30.89 μM |
54 | Cl | H | Propyl | H | 22 | −6.60 | −6.13 | −6.20 | 26.68 μM |
55 | Cl | Methyl | Propyl | H | 18 | −6.16 | −5.99 | −5.95 | 5.68 μM |
56 | Methyl | H | Methyl | H | 30 | −6.89 | −6.97 | −6.95 | 8.78 μM |
57 | Methyl | Methyl | Methyl | H | 40 | −7.02 | −7.62 | −7.65 | 14.74 μM |
58 | Methyl | OH | Methyl | H | 16 | −7.95 | −8.45 | −8.42 | 644.00 nM |
59 | Methyl | Methoxy | Methyl | H | 34 | −7.35 | −7.85 | −7.81 | 13.69 μM |
60 | Methyl | Cl | Methyl | H | 30 | −7.97 | −8.35 | −8.38 | 820.00 μM |
61 | Methyl | Methoxy | Propyl | H | 36 | −5.32 | −5.65 | −5.62 | 67.00 μM |
62 | Methyl | OH | Propyl | H | 24 | −6.24 | −6.01 | −6.05 | 38.12 μM |
63 | Methyl | H | Propyl | H | 52 | −6.18 | −6.45 | −6.44 | 18.74 μM |
64 | Methyl | Methyl | Propyl | H | 58 | −6.57 | −6.76 | −6.71 | 11.04 μM |
In addition to the physicochemical properties of the inhibitor, molecular shape is also expected to be extremely important for accessing and fitting within the active site of the aromatase. Docking analyses revealed that the novel hybrid compounds obtained from this investigation could fit well within the binding site cavity (Fig. 1) to form a three-branched shape similar to most third generation aromatase inhibitors.41 As shown previously, a high degree of hydrophobicity and the potential to establish hydrogen bonds with the aromatase enzyme are the basic physicochemical features of known aromatase inhibitors.42 These properties are related to the non-polar binding pocket of the enzyme which is dominated by aliphatic amino acid residues, such as Met 374, Val 373, Val 370, Ile 305, Ala 306, Ile 133, Trp 224, Leu 372, Leu 477, Phe 134 and Thr 310.42 The newly designed compounds could well make the required hydrophobic interactions through their quinazolinone, benzofuran and imidazole moieties. Examination of the best-ranked docking reveals that among the three heteroaromatic rings located in the substrate cavity, relatively, the quinazolinone nucleus was in close proximity with the heme iron. This feature reflects the binding mechanism found for this type of molecule, which explains binding through heterocyclic aromatic coordination to the heme iron of the P450 active site. However, no H-bond was predicted for N3 of quinazoline from the docking results. Additionally, π–π conjugate interactions are formed among Phe 221, Trp 224, Ile 133 and the phenyl ring of quinazolinone. The benzofuran ring of the hybrid compounds is found to bind through hydrophobic interactions with Phe 134, Val 370, Leu 372 and Met 374. In some cases the oxygen atom of the benzofuran ring formed a hydrogen bond with Asp 309, which is in close proximity to the benzofuran. The imidazole ring was oriented to make hydrophobic interactions with the Trp 224, Val 369 and Thr 310 residues. The hydrophobic interaction observed for the cationic imidazole ring is interesting. The reason that can be stated is that this cation isn’t a simple point charge that can only be involved in electrostatic interactions but is a delocalized charge on the imidazolium ring. This delocalization leads to distribution of a positive charge on the five atoms of the ring such that carbon and hydrogen are still able to participate in hydrophobic interactions, as seen in the figure drawn using Ligplot. It should be mentioned that with this scaffold, due to the presence of a positive charge on the nitrogen atom of imidazole, it seems logical for it to be accommodated away from the heme iron. Also, because of the resonance effect, none of the nitrogens of the imidazole could make hydrogen bonds with the residues in the aromatase active site.
Fig. 1 The binding mode of the new hybrid scaffold in the active site of aromatase, obtained from AutoDock4: (a and b) 3D structure, (c) 2D structure. |
In the next step, butyl, halogen, hydroxyl and methoxy groups were appropriately substituted on the heteroaromatic rings of the designed compounds and their binding modes were investigated through docking. The new analogs share the same binding mode, similar to the unsubstituted derivatives, with an extra anchoring point which might cause stronger biding to the aromatase active site. The main binding modes in these complexes can be described as follows: introduction of OH on the benzofuran ring caused formation of hydrogen bonds with the amino acid residue Ser 478 of aromatase (Fig. 2), while an OH group on quinazolinone ring exhibited a H-bond with Met 374 (Fig. 3). In the structure with two hydroxyl groups substituted on both quinazolinone and benzofuran, both of the mentioned H bonds are visible (Fig. 4).
Fig. 2 Binding mode and hydrogen bond interaction of hydroxyl group on the benzofuran ring with Ser 478 in the aromatase active site: (a) 3D structure, (b) 2D structure. |
Fig. 3 Binding mode and hydrogen bond interaction of hydroxyl group on the quinazolinone ring with Met 374 in the aromatase active site: (a) 3D structure, (b) 2D structure. |
Insertion of methoxy groups on the quinazolinone and/or benzofuran systems also caused the formation of H-bonds with hydrogen-bond donors in the active site of the enzyme. The methoxy group on the benzofuran ring and Ser 478 are in close proximity (Fig. 5) with favorable hydrogen bonding interactions, while Met 374 formed a hydrogen bond with the methoxy group on the quinazolinone ring (Fig. 6). Additionally, hydrophobic substituents, such as chloro and methyl groups, on the quinazolinone and benzofuran rings somewhat increase aromatase inhibition potency through increased hydrophobic interactions. Chloro and methyl substituents on the quinazolinone ring are stabilized by van der Waals interactions with the non-polar amino acids in the active site (Ala 306, Thr 310, Trp 224, Val 370, Ile 133, Phe 134, Leu 372 and Val 373), while chloro and methyl on the benzofuran ring are involved in hydrophobic interactions with the hydrophobic residues Ile 133, Phe 134, Leu 372 and Val 369. Finally, substitution of a small alkyl group such as methyl on the imidazolium ring make its access to and accommodation within the active site easier compared to the bulkier butyl group.
Fig. 5 Binding mode and hydrogen bond interaction of a methoxy group on the benzofuran ring with Ser 478 in the aromatase active site: (a) 3D structure, (b) 2D structure. |
Fig. 6 Binding mode and hydrogen bond interaction of a methoxy group on the quinazolinone ring with Met 374 in the aromatase active site: (a) 3D structure, (b) 2D structure. |
The structures with the best Gibbs interaction energies obtained through the previous docking experiments were re-optimized using the three-layer ONIOM method. In other words, the structure of the ligand, along with the interacting residues and heme-iron, was optimized in the electrostatic field of the rest of the protein, which is frozen. The advantage of this method is that the ligand is optimized in the presence of the protein and its structure is more reliable than before. In addition, the geometrical structures of the residues interacting with the ligand are improved by this optimization. For the high-level layer we used the designed inhibitors inside the aromatase active site, for the medium-level we used the heme-iron and several residues from the binding pocket, including as Met 374, Val 373, Val 370, Ile 305, Ala 306, Ile 133, Trp 224, Leu 372, Leu 477, Phe 134 and Thr 310 and the remainder of the protein is treated as the low-level layer, as depicted in Fig. 7. After the geometric parameters of the molecules are fully optimized, the ligands were extracted from the complexes and were subjected to rigid re-docking. This means that all rotatable bonds were to be held constant and the ligand charges were replaced with obtained values from the QM/MM calculation. The results showed that all runs extended to creation of one cluster and also the calculated ΔGs increases by as much as 1–2 kcal mol−1 for the ligands (Table 2). In the early QM/MM calculation we only used the obtained charges by this method for improving the docking results, but in the final QM/MM calculation the optimized ligand in the protein environment was used for AutoDock. After that the results showed that if the structure of ligand was optimized in the active sit of the enzyme and then the same optimized ligand was re-docked to protein; the lower binding energy was obtained.
No | ΔGb docking after refitting of ligand charge | Ki after refitting of ligand charge | RMSD | ΔGb in rigid docking | Ki in rigid docking | RMSD |
---|---|---|---|---|---|---|
4 | −8.82 | 409 nM | 0.31 | −10.11 | 37.00 nM | 0.14 |
6 | −9.34 | 68 nM | 0.27 | −11.47 | 4.39 nM | 0.16 |
16 | −8.86 | 411 nM | 0.25 | −9.94 | 65.00 nM | 0.14 |
17 | −8.85 | 423 nM | 0.36 | −9.75 | 66.00 nM | 0.11 |
34 | −8.56 | 530 nM | 0.41 | −9.95 | 65.00 nM | 0.21 |
48 | −8.45 | 642 nM | 0.28 | −10.81 | 9.69 nM | 0.11 |
49 | −8.59 | 531 nM | 0.25 | −10.41 | 12.44 nM | 0.16 |
58 | −8.42 | 644 nM | 0.38 | −10.23 | 17.00 nM | 0.12 |
According to the above results, the designed ligands formed three-branched structures in the protein environment, which were accommodated well into the active site. We found that in addition to the van der Waals interactions, hydrogen bonds are key factors for the ligand–receptor interactions. Compound 6, which yielded the highest ΔGb (−11.47) and the best performing Ki (48 nM) value, was assumed to be the best ligand. Since the assay method used for evaluating the anti-aromatase potencies of the compounds was theoretical, the precise assessment of aromatase inhibition as the mechanism of action for these compounds needs further studies. Based on the results of the present research, synthesis and cytotoxic evaluation of the designed derivatives are under way in our research group.
This journal is © The Royal Society of Chemistry 2015 |