In a structure-activity relationship (SAR) study, 3-methoxy-1,4-phenanthrenequinones, calanquinone A (6a), denbinobin (6b), 5-OAc-calanquinone A (7a) and 5-OAc-denbinobin (7b), have significantly promising cytotoxicity against various human cancer cell lines (IC(50) 0.08-1.66 µg/mL). Moreover, we also established a superior pharmacophore model for cytotoxicity (r = 0.931) containing three hydrogen bond acceptors (HBA1, HBA2 and HBA3) and one hydrophobic feature (HYD) against MCF-7 breast cancer cell line. The pharmacophore model indicates that HBA3 is an essential feature for the oxygen atom of 5-OH in 6a-b and for the carbonyl group of 5-OCOCH(3) in 7a-b, important for their cytotoxic properties. The SAR for moderately active 5a-b (5-OCH(3)), and highly active 6a-b and 7a-b, are also elaborated in a spatial aspect model. Further rational design and synthesis of new cytotoxic phenanthrene analogs can be implemented via this model. Additionally, employing a ChemGPS-NP based model for cytotoxicity mode of action (MOA) provides support for a preliminary classification of compounds 6a-b as topoisomerase II inhibitors.
In a structure-activity relationship (SAR) study, 3-methoxy-1,4-phenanthrenequinones, <span class="Chemical">calanquinone A (6a), denbinobin (6b), 5-OAc-calanquinone A (7a) and 5-OAc-denbinobin (7b), have significantly promising cytotoxicity against various humancancer cell lines (IC(50) 0.08-1.66 µg/mL). Moreover, we also established a superior pharmacophore model for cytotoxicity (r = 0.931) containing three hydrogen bond acceptors (HBA1, HBA2 and HBA3) and one hydrophobic feature (HYD) against MCF-7 breast cancer cell line. The pharmacophore model indicates that HBA3 is an essential feature for the oxygen atom of 5-OH in 6a-b and for the carbonyl group of 5-OCOCH(3) in 7a-b, important for their cytotoxic properties. The SAR for moderately active 5a-b (5-OCH(3)), and highly active 6a-b and 7a-b, are also elaborated in a spatial aspect model. Further rational design and synthesis of new cytotoxicphenanthrene analogs can be implemented via this model. Additionally, employing a ChemGPS-NP based model for cytotoxicity mode of action (MOA) provides support for a preliminary classification of compounds 6a-b as topoisomerase II inhibitors.
Natural phenanthrenes are probably generated from photochemical cyclization of <span class="Chemical">stilbenes [1]. More than 270 phenanthrenes have been isolated from natural products, especially the family Orchidaceae, and some of them possess various biological activities, including cytotoxicity, antiplatelet aggregation, anti-inflammatory, antimicrobial, spasmolytic, anti-allergic activities and phytotoxicity [1]. In our previous study, calanquinone A [5-hydroxy-3,6,7-trimethoxy-1,4-phenanthrenequinone (5-hydroxy-3,6,7-trimethoxy-1,4-PQ)] (Figure 1), a new PQ isolated from Calanthe arisanensis in 2008, showed significant cytotoxic activity against human lung (A549), prostate (PC-3 and DU145), colon (HCT-8), breast (MCF-7), nasopharyngeal (KB), and vincristine-resistant nasopharyngeal (KBVIN) cancer cell lines with EC50 values of 0.03–0.45 µg/mL [2], [3]. Calanquinone A is related in structure to other naturally occurring cytotoxicPQs, including denbinobin (Figure 1) (5-hydroxy-3,7-dimethoxy-1,4-PQ), sphenone A (3,6,7-trimethoxy-1,4-PQ), densiflorol B (7-hydroxy-2-methoxy-1,4-PQ), and annoquinone A (3-methoxy-1,4-PQ) [1]. Denbinobin, first isolated from Dendrobium nobile in 1981, is the only 1,4-PQ that has been studied in terms of the cytotoxic mechanisms of humancolon (HCT-116 and COLO 205), lung adenocarcinoma (A549), myelogenous leukemia (K562), and pancreatic adenocarcinoma (BxPC-3) cancer cell lines [4]–[11]. The data implied that denbinobin could be a potential anticancer lead compound.
Figure 1
Structures of calanquinone A and denbinobin.
In our preliminary results of cytotoxic structure-activity relationship (SAR) studies, <span class="Chemical">calanquinone A (6a) displayed an up to 7-fold greater cytotoxic activity than denbinobin (6b), which is known as a potent cytotoxic agent [4]–[11], toward human liver (HepG2 and Hep3B), oral (Ca9-22), lung (A549) and breast (MEA-MB-231 and MCF7) cancer cell lines. Up to now, the SAR of PQs and phenanthrenes has only rarely been reported and is worthy of further study. In this research, calanquinone A (6a), denbinobin (6b) and their derivatives were synthesized [3], [12], [13] and evaluated for cytotoxic activity. In addition, employing a ChemGPS-NP based model provides the prediction for cytotoxicity mode of action (MOA) of calanquinone A (6a) and denbinobin (6b).
Results and Discussion
Chemistry
Eleven natural phenanthrene analogs (<span class="Gene">CA-1-11) (Figure 2) were isolated from C. arisanensis, and calanquinone A (CA-1) exhibited the highest potency against humancancer cell lines [2], [3]. According to the previous results [3], calanquinone A (CA-1) was selected as a lead compound and its derivatives were then synthesized for this study.
Figure 2
Structural sets used in the pharmacophore study.
CA-1∼CA-11 and 3a-9b are noted as natural and synthesized compounds, respectively.
Structural sets used in the pharmacophore study.
CA-1∼<span class="Gene">CA-11 and 3a-9b are noted as natural and synthesized compounds, respectively.
We modified the synthetic procedure of Dr. Kraus and his co-workers [12], [13] to synthesize all phenanthrene derivatives. <span class="Chemical">As shown in Figure 3, 2-aldehyde-1,4-quinone was prepared by DDQ oxidation of commercially available 2,5-dihydroxybenzaldehyde. The quinone was coupled with 3,4,5-trimethoxytoluene and 3,5-dimethoxytoluene in the presence of 1 equivalent of trifluoroacetic acid to produce 1a and 1b, respectively. Compounds 1a and 1b were methylated with Me2SO4 in the presence of K2CO3 (acetone, 60°C, 5 h) to give the desired 2a and 2b. Cyclization of 2a and 2b with P4-tBu (benzene, 110 or 140°C, 19–29 h) gave phenanthrenes 3a and 3b, which were oxidized with AgO (6 N HNO3, acetone, 50°C, 2–3 min) to phenanthrenequinone 4a, 4b and 4c. Addition of methanol to 4a, 4b and 4c catalyzed by ferric sulfate [14] gave 5a–f, respectively (Figure 4). Compounds 7a and 7b were obtained by treatment of 5a and 5b with TMSI (CH2Cl2, RT or 60°C, monitored by TLC) to give calanquinone A (6a; CA-1) and denbinobin (6b), followed by treatment with Ac2O (pyridine, RT, overnight) to selectively remove the methyl group and incorporate an acetyl group at C-5, respectively (Figure 5). It is a characteristic feature of the angular arrangement of 1,4-phenanthrenequinones which led to remarkable selectivity in the cleavage of sterically hindered methyl ether at C-5 even in preference of that at C-3. However, applying TMSI to remove the methyl groups in phenanthrenes 3a and 3b was unsuccessful. Finally, cleavage of the methyl ether groups in 3a and 3b with AlCl3 generated compounds 8a & 9a and 8b & 9b, respectively (Figure 6). The excess AlCl3 regioselectively cleaved the methyl ethers only at C-4 and C-5 or C-6 in order to release the steric strain.
Figure 3
Synthetic procedure of phenanthrene derivatives.
Reagents and conditions: (i) DDQ, benzene, RT. (ii) TFA, ether, RT. (iii) Me2SO4, K2CO3, acetone, reflux. (iv) P4-tBu, benzene, 140°C. (v) AgO, 6 N HNO3, acetone, 60°C.
Figure 4
Synthesis of methoxy-phenanthrenequinones.
Reagents and conditions: (i) MeOH, Fe2(SO4)3, PTSA, 70°C.
Figure 5
Selective demethylation and acetylation of phenanthrenequinones.
Reagents and conditions: (i) TMSI, CH2Cl2, RT or 60°C. (ii) Ac2O, pyridine, RT.
Figure 6
Demethylation of phenanthrenes.
Reagents and conditions: (i) AlCl3, benzene, 70°C.
Synthetic procedure of phenanthrene derivatives.
Reagents and conditions: (i) DDQ, <span class="Chemical">benzene, RT. (ii) TFA, ether, RT. (iii) Me2SO4, K2CO3, acetone, reflux. (iv) P4-tBu, benzene, 140°C. (v) AgO, 6 N HNO3, acetone, 60°C.
Synthesis of methoxy-phenanthrenequinones.
Reagents and conditions: (i) MeOH, <span class="Chemical">Fe2(SO4)3, PTSA, 70°C.
Selective demethylation and acetylation of phenanthrenequinones.
Reagents and conditions: (i) TMSI, <span class="Chemical">CH2Cl2, RT or 60°C. (ii) Ac2O, pyridine, RT.
Demethylation of phenanthrenes.
Reagents and conditions: (i) AlCl3, <span class="Chemical">benzene, 70°C.
Cytotoxicity
The cytotoxic <span class="Chemical">assay of 11 naturally occurring and 19 synthesized phenanthrenes was carried out on a diverse set of human liver (HepG2 and Hep3B), oral (Ca9-22), lung (A549) and breast (MEA-MB-231 and MCF7) cancer cell lines, and a human fetal lung fibroblast (MRC-5) cell line (Tables 1 and 2). Doxorubicin was used as a positive control and an IC50>4 µg/mL was considered inactive.
Table 1
Cytotoxicity data of natural phenanthrenes isolated from C. arisanensis.
IC50 (µg/mL)/Cell line
Compd
HepG2
Hep3B
Ca9-22
A549
MCF-7
MDA-MB-231
MRC-5
CA-1
0.21±0.01
0.22±0.00
0.17±0.01
0.11±0.00
0.09±0.00
0.64±0.06
0.65±0.00
CA-2
14.47±0.43
11.94±0.55
14.64±0.11
18.50±0.29
11.90±0.08
>20
>20
CA-3
19.76±0.16
6.02±0.12
12.55±0.76
>20
13.25±0.00
11.90±0.02
>20
CA-4
>20
11.05±0.73
15.12±0.14
>20
14.30±0.18
19.83±0.45
>20
CA-5
17.80±0.40
9.87±0.09
12.17±0.03
>20
12.70±0.07
11.66±0.13
>20
CA-6
>20
11.50±0.76
14.56±0.31
>20
>20
>20
>20
CA-7
>20
>20
>20
>20
14.62±0.34
>20
>20
CA-8
>20
>20
>20
>20
>20
>20
>20
CA-9
>20
12.89±0.32
12.99±0.42
18.11±0.26
10.25±0.27
19.48±0.58
>20
CA-10
7.52±0.00
7.24±0.18
6.11±0.51
7.46±0.32
6.77±0.12
7.56±0.03
16.39±0.36
CA-11
7.17±0.35
6.28±0.26
5.86±0.06
7.40±0.01
4.84±0.89
7.42±0.03
13.17±0.33
Doxoa
0.22±0.04
0.49±0.00
0.16±0.02
0.81±0.04
0.67±0.03
0.74±0.01
1.94±0.02
Doxorubicin (Doxo) was used as the positive control.
Table 2
Cytotoxicity data of synthesized phenanthrenes.
IC50 (µg/mL)/Cell line
Compd
HepG2
Hep3B
Ca9-22
A549
MCF-7
MDA-MB-231
MRC-5
3a
>20
15.55±0.08
>20
>20
>20
>20
>20
3b
>20
>20
18.44±0.11
>20
>20
>20
>20
4a
15.68±0.01
7.33±0.03
2.17±0.01
15.94±0.01
13.19±0.38
4.60±0.06
11.93±0.12
4b
18.64±0.14
12.10±0.33
3.45±0.03
>20
15.63±0.01
7.09±0.01
>20
4c
6.87±0.10
4.37±0.00
1.90±0.01
12.65±0.48
4.71±0.26
4.09±0.12
5.37±0.25
5a
4.65±0.07
4.67±0.81
5.99±0.24
11.04±0.44
14.29±0.50
>20
>20
5b
12.99±0.21
15.95±0.21
15.25±0.05
15.71±0.17
14.40±0.28
18.98±0.55
>20
5c
>20
13.78±0.38
9.07±0.09
11.99±2.18
11.10±0.02
19.11±0.43
>20
5d
1.49±0.01
6.66±0.13
4.36±0.03
9.30±0.48
13.78±0.06
14.76±0.25
>20
5e
1.24±0.08
7.38±0.02
6.21±0.12
7.49±0.07
9.40±0.00
6.74±0.05
8.11±0.25
5f
19.64±0.74
7.88±0.19
4.57±0.03
>20
14.77±0.21
4.74±0.02
>20
6a
0.08±0.00
0.19±0.01
0.59±0.01
0.14±0.02
0.20±0.00
0.89±0.01
0.60±0.05
6b
0.23±0.01
0.34±0.00
0.68±0.00
0.99±0.02
0.26±0.00
1.06±0.03
2.14±0.08
7a
11.30±0.14
0.36±0.14
0.84±0.05
0.60±0.01
0.16±0.00
1.13±0.05
1.14±0.01
7b
13.23±0.25
0.60±0.11
1.55±0.20
1.66±0.01
0.53±0.03
1.61±0.04
2.64±0.06
8a
19.38±0.56
5.77±0.63
3.91±0.07
>20
17.01±0.42
15.63±1.52
>20
9a
>20
>20
>20
>20
>20
>20
>20
8b
10.08±0.66
6.93±0.19
12.07±0.37
19.58±0.22
19.31±0.08
>20
>20
9b
>20
>20
>20
>20
>20
>20
>20
Doxoa
0.22±0.02
0.42±0.01
0.14±0.04
0.63±0.15
0.35±0.09
1.12±0.05
1.94±0.02
Doxorubicin (Doxo) was used as the positive control.
Doxorubicin (<span class="Chemical">Doxo) was used as the positive control.
Doxorubicin (<span class="Chemical">Doxo) was used as the positive control.
Among the 11 naturally occurring compounds [2], calanquinone A (<span class="Gene">CA-1) (5-OH, 6-OCH3) and calanquinone B (CA-2) (5-OCH3, 6-OH) simply have reversed placements of the OH and one OCH3 group, but CA-2 was much less potent than CA-1 (Figure 2 and Table 1). The SAR results of CA-1 and CA-2 could possibly be explained by intramolecular hydrogen bonding between C = O (C-4) and OH (C-5) groups in 3-methoxy-1,4-PQs that may be a necessary moiety for cytotoxicity. To set up SAR correlations and identify active phenanthrene analogs, calanquinone A (CA-1) was selected as a lead compound for further studies.
Accordingly, 19 analogues including calanquinone A (6a; <span class="Gene">CA-1) were synthesized and tested in cytotoxicityassays. As shown in Table 2, calanquinone A (6a) and denbinobin (6b) exhibited significant potency against all cancer cell lines (IC50 0.08–1.06 µg/mL). PQs 7a and 7b also showed very high potency against five cancer cell lines (IC50 0.16–1.66 µg/mL), not including HepG2. Conversely and interestingly, PQs 5d and 5e were active only against the HepG2cancer cell line with IC50 values of 1.49 and 1.24 µg/mL, respectively. PQs 4a, 4b, 4c and phenanthrene 8a displayed selective activity toward the Ca9-22cancer cell line with IC50 values of 2.17, 3.45, 1.90 and 3.91 µg/mL, respectively. The SAR study of cytotoxicity suggested that the skeleton of 1,4-PQ is preferable to that of phenanthrene. To evaluate a potential SAR of the intramolecular hydrogen bond between C-4 and C-5, 3-methoxy-1,4-PQs 5a–b, 6a–b and 7a–b were designed. Compounds 6a and 6b, with OH at C-5 and C = O at C-4, can form an intramolecular hydrogen bond. However, the hydrogendonors of 5a–b and 7a–b have been replaced with OCH3 and OAc groups, respectively. Among the six 3-methoxy-1,4-PQs, 6a and 6b exhibited the most significant potency, especially 6a (IC50 0.08–0.89 µg/mL). Compounds 5a and 5b showed marginal activities against all humancancer cell lines. Surprisingly, the new 7a and 7b, with OAc at C-5 and C = O at C-4, were active toward five humancancer cell lines (IC50 0.16–1.66 µg/mL), but not HepG2. These data also represent the first time we have found this phenomenon in a cytotoxicassay of PQ derivatives. To expand upon the SAR study, all natural and synthesized compounds were used for the 3D pharmacophore model building.
3D Pharmacophore Modeling
To further identify the critical structural features of the phenanthrene analogs, 29 compounds (Figure 2) were used for pharmacophore modeling with Catalyst HypoGen. In this spatial <span class="Chemical">aspect model, the phenanthrene structures and their cytotoxicity toward MCF-7 cancer cell line showed some interesting information.
The best pharmacophore model was established <span class="Chemical">as a result of thirty runs with various parameters and characterized by a best correlation coefficient (0.931), the lowest total cost value (109.366), the highest cost difference (42.417), and the lowest RMS (0.790) (for details see Tables S1, S2 and Figure S1; Text S1). As shown in Figure 7, four essential features, three hydrogen bond acceptors (HBA1, HBA2 and HBA3) and one hydrophobic feature (HYD) were defined. All mutual distances of the four features can be measured. The distances between HBA1 and HBA2 or HBA3 were found to be 5.13 and 7.40 Å, respectively. The distances between HBA2 and HBA3 or HYD were found to be 5.95 and 5.92 Å, respectively. The distances between HYD and HBA1 or HBA3 were found to be 3.83 and 4.78 Å, respectively. The distance between HBA2 and HBA3 is especially critical for the MCF-7cytotoxic effect in this model.
Figure 7
The mutual distances of the hydrogen bond triad and the distances between hydrogen bonds and hydrophobic group in the best HypoGen pharmacophore model.
The pharmacophore features are colored with green, as are the hydrogen-bond acceptors (HBA1, HBA2 and HBA3). The hydrophobic aromatic feature (HYD) is denoted in cyan. Distances between features are in Ångström units.
The mutual distances of the hydrogen bond triad and the distances between hydrogen bonds and hydrophobic group in the best HypoGen pharmacophore model.
The pharmacophore features are colored with green, <span class="Chemical">as are the hydrogen-bond acceptors (HBA1, HBA2 and HBA3). The hydrophobic aromatic feature (HYD) is denoted in cyan. Distances between features are in Ångström units.
The mappings of the best model with all compounds show fit values ranging from 6.18 to 8.57 (Table S3; Text S1). Calanquinone A (6a) mapped to the best hypothesis model with the fit value of 8.57 reveals significant features in Figure 8A. Obviously, the <span class="Gene">HBA1 links to the carbonyl group of quinone ring at C-1, HBA2 links to the oxygen atom of the methoxyl group at C-3, HBA3 links to the oxygen atom of the hydroxyl group at C-5, and HYD aligns to the aromatic ring (B-ring). Compound 6b has a similar alignment as 6a, with a high fit value of 7.95 (Figure 8B and Table S3). As shown in Figures 8C and 8D,
, which are highly toxic to MCF-7 cells, also match against all features of the best hypothesis in which 7a was originally designed to remove the intramolecular hydrogen bond and was previously speculated to be less cytotoxic. The main difference in structure between 7a–b and 6a–b is the carbonyl group of the acetoxyl substituent at C-5. Consequently, the conserved distance between HBA2 and HBA3 explains why 7a and 7b, with OCOCH3 at C-5 and C = O at C-4 but without the same intramolecular hydrogen bonds as 6a and 6b, can still possess significant cytotoxicity, in contrast to our previous speculation. For less MCF-7cytotoxic compounds, mismatching one hydrogen bond in the triad and thus disrupting the structures of 5a and 5b results in moderate inhibition, one order higher in µg/mL. The oxygen atom of the methoxyl group at C-6, the carbonyl group at the C-4 position, and the aromatic atom at C-9 of 5a fit into the HBA2, HBA3 and HYD, respectively. However, they do not fit into HBA1 (Figure 8E). Also, the carbonyl group at C-1, the oxygen atom of the methoxyl group at C-3, and the aromatic atom at C-9 of 5b match against HBA1, HBA2 and HYD features, but are not linked to HBA3 (Figure 8F).
Figure 8
The best HypoGen pharmacophore model mapping onto calanquinone A (6a) and its derivatives.
The light and dark green represent active and inactive features, respectively. The models mapped with the compounds 6a (A), 6b (B), 7a (C), 7b (D), 5a (E) and 5b (F) are shown here. Pharmacophore features are colored as in Figure 7.
The best HypoGen pharmacophore model mapping onto calanquinone A (6a) and its derivatives.
The light and dark green represent active and inactive features, respectively. The models mapped with the compounds 6a (A), 6b (B), 7a (C), 7b (D), 5a (E) and 5b (F) are shown here. Pharmacophore features are <span class="Species">colored as in Figure 7.
Thus, HBA1, 2 and 3 complete the triad of the <span class="Chemical">hydrogen acceptor feature and clearly explain the MCF-7cytotoxic variation of 5a–b, 6a–b and 7a–b. In addition, the hydrophobic feature HYD indicates a pharmacophore anchor for a three-ring core as in phenanthrene and PQ derivatives in the series. The hydrophobic feature links to all compounds and the HBA feature always links to the carbonyl group at C-1 or C-4 of the quinone ring in all PQs. As a whole, a pharmacophore and explicit SAR were established herein.
ChemGPS-NP Analysis of Calanquinone A (6a) and Denbinobin (6b)
ChemGPS-NP (chemical global positioning system for natural products) is a computational model based on principal component analysis of physical-chemical properties. Such properties can be estimated directly from structure data, and by performing score prediction in the ChemGPS-NP model, this provides a versatile tool for charting and navigating the biologically relevant chemical space [15]. In a previous study ChemGPS-NP h<span class="Chemical">as successfully been used to chart a set of known anticancer agents with different cytotoxic mechanisms. The resulting map has been used as a tool to predict the anticancer Mode of Action (MOA) for new and previously unstudied lead compound [16]. As shown in Figure 9, the two most potent cytotoxic compounds, calanquinone A (6a) and denbinobin (6b), were predicted in the model. Evaluating their resulting position on the chemical space map, it can be concluded that these phenanthrene derivatives do not unambiguously belong to any of the well defined groups representing alkylating agents, antimetabolites, proteasome inhibitions, tyrosine kinase inhibitors, topoisomerase I, and tubulin inhibitors except topoisomerase II inhibitors. The preliminary result of this ChemGPS-NP analysis indicates that calanquinone A (6a) and denbinobin (6b) might be members of a topoisomerase II inhibitor, which however, still remains to be further elucidated.
Figure 9
ChemGPS-NP analysis of calanquinone A (6a) and denbinobin (6b).
Score plot of the three dimensions (principal components 2–4) consisting of PC2 (yellow; aromaticity etc.), PC3 (green; lipophilicity etc.) and PC4 (orange; flexibility/rigidity), from analysis of most potent compounds 6a and 6b as medium seagreen cubes in the ChemGPS-NP model addressed by Rosén et al. in 2009 for prediction of MOA. A reference set of known anticancer agents includes alkylating agents (red), anti-metabolites (lime), proteasome inhibitions (cyan), tyrosine kinase inhibitors (orange), topoisomerase I (blue), topoisomerase II (magenta), and tubulin inhibitors (black).
ChemGPS-NP analysis of calanquinone A (6a) and denbinobin (6b).
Score plot of the three dimensions (principal components 2–4) consisting of PC2 (yellow; aromaticity etc.), <span class="Gene">PC3 (green; lipophilicity etc.) and PC4 (orange; flexibility/rigidity), from analysis of most potent compounds 6a and 6b as medium seagreen cubes in the ChemGPS-NP model addressed by Rosén et al. in 2009 for prediction of MOA. A reference set of known anticancer agents includes alkylating agents (red), anti-metabolites (lime), proteasome inhibitions (cyan), tyrosine kinase inhibitors (orange), topoisomerase I (blue), topoisomerase II (magenta), and tubulin inhibitors (black).
Topoisomerase II Assay
From the ChemGPS-NP analysis, it seems that the MOA for 3-methoxy-1,4-PQs might processes <span class="Disease">cytotoxicityas a topoisomerase II inhibitor. In the previous study [17], 1,4-benzoquinone has been found to poison human topoisomerase IIα. According to these results, we chose the most potential calanquinone A (6a) and its moderate compound 5a to test the DNA cleavage assay, in which known etoposide (VP-16) was used as the positive control. As shown in Figure 10, compound 6a showed the inhibition on Topo IIα in the result of the appearance of supercoided DNA instead of the relaxed one at the concentration of 100 µmol/L. Additionally, compound 5a also had the similar effect at higher concentration (200 µmol/L). Moreover, both compounds induced the formation of linear DNA, suggesting that they could possibly trap Topo IIα into DNA cleavage complex. Our data proved PQs inhibit hTopoII in vitro with inducing DNA strand breaks and protein covalently bound to DNA, ultimately leading to cell cycle arrest and death.
Figure 10
Topoisomerase II DNA cleavage assay.
An in vitro assay was used to assess the effect of compounds 6a and 5a on the DNA cleavage activity of human TopoII. Etoposide was the positive control. Control lane: TopoIIα + plasmid DNA. DMSO lane: TopoIIα + plasmid DNA + DMSO.
Topoisomerase II DNA cleavage assay.
An in vitro assay w<span class="Chemical">as used to assess the effect of compounds 6a and 5a on the DNA cleavage activity of human TopoII. Etoposide was the positive control. Control lane: TopoIIα + plasmid DNA. DMSO lane: TopoIIα + plasmid DNA + DMSO.
Conclusion
In summary, a series of phenanthrene derivatives, including the new derivatives (3a, 4c, 5c–f, <span class="Chemical">7a, 8a–b and 9a–b), were synthesized in this investigation. On the basis of our SAR studies, 3-methoxy-1,4-PQs 6a (calanquinone A), 6b (denbinobin), 7a (5-OAc-calanquinone A) and 7b (5-OAc-denbinobin) were identified as highly potent cytotoxic agents.
A best ligand-based pharmacophore model against the <span class="Disease">MCF-7 cancer cell line was successfully established. It explains the SAR of 3-methoxy-1,4-PQs 5a–b (5-OCH3), 6a–b (5-OH) and 7a–b (5-OCOCH3) in a spatial aspect model. Highly active 6a, 6b, 7a and 7b possess three hydrogen bond acceptors forming a hydrogen bond triad combined with one hydrophobic group as a pharmacophore that can interact with a potential target. The revealed pharmacophore model provides a bona fide basis for further design and synthesis of promising phenanthrene structures in vitro to study their anti-breast cancer properties. On the basis of ChemGPS-NP prediction and TopoII assay assessment, 1,4-PQs were suggested as the topoisomerase II inhibitors. This is the first time to apply ChemGPS-NP to previously untested cytotoxic compounds for MOA prediction. In the future, ChemGPS-NP could be used to effectively find the most possible MOA in the new drug discovery, as suggested by Rosén and co-workers [16].
Overall, our data demonstrate that PQs could be promising lead compounds for the further development of anti-<span class="Disease">cancer.
Materials and Methods
General
Melting points were determined on a Yanaco® digital micro melting point apparatus model MP-500D without correction. NMR spectra were recorded on Varian Unity-plus 400 MHz FT-NMR and Varian Mercury-plus 400 MHz FT-NMR instruments. Chemical shift (δ) values are in ppm (parts per million) with <span class="Chemical">CDCl3as the internal standard, and coupling constants (J) are in Hz. HRESI-MS and ESI-MS measurements were performed on a Bruker Daltonics APEX II 30e mass spectrometer. TLC was performed on Kieselgel 60, F 254 (0.25 nm, Merck), and spots were viewed under ultraviolet light at 254 and 356 nm. For column chromatography, silica gel (Kieselgel 60, 70–230, and 230–400 mesh, Merck) and a Biotage® SP system apparatus were used.
Cytotoxicity ASSAY
Compounds were tested against human liver (<span class="CellLine">HepG2 and Hep3B), oral (Ca9-22), lung (A549), breast (MEA-MB-231 and MCF7) cancer cell lines, and the human fetal lung fibroblast (MRC-5) cell line using an established colorimetric MTTassay protocol [18]. The absorbance was measured at 550 nm using a microplate reader. The IC50 is the concentration of agent that reduced cell growth by 50% under the experimental conditions.
3D Pharmacophore Model
The pharmacophore modeling with Catalyst HypoGen was performed via Discovery Studio 2.1 (Accelrys, San Diego, CA, USA) [19]. Twenty-nine <span class="Chemical">phenanthrene derivatives were collected from the natural plant, C. arisanensis, and from chemical synthesis (Figure 2). Cytotoxicity against MCF-7 cells was determined by the MTTassay and the concentration (µg/mL) of test compound which inhibited 50% of the cancer cells (IC50) was used in the generation of the pharmacophore model. An IC50 value of >20 µg/mL was defined as 20 µg/mL. All experimental IC50 values spanned about 2–3 orders of magnitude from 0.09 to 20 µg/mL. The 2D/3D structures of compounds were generated using ChemBioOffice 2008 (Cambridge Scientic Computing, Cambridge, Massachusetts, USA) and then optimized in a Dreiding force field. The conformational ensemble of each compound was generated using the best conformational analysis method based on a CHARMM force field with a 20 kcal/mol energy threshold above the global minimum. A maximum limit of 255 conformations was used to cover maximum conformational space. The best 3D arrangements of chemical functionalities should explain the activity variations among the 29 compounds. Thirty runs with different parameters were performed for the best pharmacophore hypothesis. Four chemical features, including hydrogen-bond acceptor (HBA), hydrogen-bond donor (HBD), hydrophobic (HYD), and aromatic ring (AR) features, were also tested during the building of pharmacophore hypotheses (Table S1; Text S1). The best hypotheses were selected via a correlation and a cost analysis in Catalyst HypoGen.
Three costs including the total cost (the sum of weight, error and configuration cost), the null cost and the fixed cost will be evaluated. A total cost that is similar to the fixed cost and far from the null cost indicates statistically significant pharmacophore hypotheses. A difference between the total cost and null cost ranging from 40 to 60 indicates a true correlation of the pharmacophore hypothesis with 75–90% high probability. The true correlation represents <50% probability when it is less than 40. Generally, the configuration cost should be smaller than 17 in a standard HypoGen model. According to the total cost (109.366), fixed cost (99.558), null cost (151.783), RMS value (0.790), and correlation coefficient (0.931) (Table S2), the best pharmacophore hypothesis, run 22, containing three <span class="Chemical">hydrogen-bond acceptors (HBA1, HBA2, HBA3) and one hydrophobic feature (HYD) was selected (Figure 7).
ChemGPS-NP
The PCA-based model ChemGPS-NP (http://chemgps.bmc.uu.se) is a tool for navigation in biologically relevant chemical space. It h<span class="Chemical">as eight principal components (PC; dimensions), derived from 35 molecular descriptors describing physical-chemical properties such as size, shape, polarizability, lipophilicity, polarity, flexibility, rigidity, and hydrogen bond capacity for a reference set of compounds. The ChemGPS-NP descriptors were calculated for compounds 6a and 6b on the basis of their structure information as simplified molecular input line entry specification (SMILES) using the software DRAGON Professional. Compounds 6a and 6b were then mapped onto ChemGPS-NP using interpolation in terms of PCA score together with a reference set of known anticancer agents with previously studied Mode of Action (MOA) (Anticancer Agent Mechanism Database; http://dtp.nci.nih.gov/docs/cancer/searches/standard_mechanism.html). Principal component and PCA score prediction were calculated employing the software SIMCA-P+, with the training set ChemGPS-NP. Prior to PCA determination, all data were centered and scaled to unit variance [16].
Topoisomeras II Assay
Topoisomerase II assay was performed by using a Topo II Drug Screening Kit (TopoGEN, Inc.). In brief, 0.1 µg of pHOT plasmid DNA was incubated with 2 units of topoisomerase IIα in 20 µL assay buffer at 37°C for 40 min in the presence of tested compounds (6a, 5a) and control drug, etoposide, respectively. 2 µL of 10% SDS and 2.5 µL of 10 mg/mL proteinase K were added into the reaction sample and then incubated for 30 min at 37°C to digest topoisomerase IIα. The samples were mixed with 2 µL of loading buffer and cleaned up by adding an equal volume of phenol:chloroform:isoamyl alcohol (25∶24∶1) according to the description. The sample was mixed by vortex and centrifuge for 10 sec. An aliquot (10 µL) of the upper aqueous part was analyzed by electrophoresis with 2% agarose gel containing 0.5 µg/mL of ethidium bromide [20].
2,5-dihydroxybenzaldehyde (3.01 g, 21.79 mmol) and <span class="Chemical">DDQ (9.18 g, 40.42 mmol) were dissolved in dry benzene. The mixture was stirred at RT overnight, and the benzene was then evaporated. The crude product reacted with equimolar amounts of 3,4,5-trimethoxytoluene and trifluoroacetic acid in Et2O at RT for 24 h or for an extended reaction time to obtain the best yield. The crude product was poured onto ice water and extracted three times with EtOAc. The EtOAc layer was washed once with brine and dried with Na2SO4. The product was chromatographed on silica gel and eluted with CH2Cl2 to give compound 1a (2.14 g, 31.16%). 1H NMR (CDCl3): δ 2.01 (s, 3 H), 3.65 (s, 3 H), 3.87 (s, 3 H), 3.91 (s, 3 H), 4.65 (s, OH), 6.68 (s, 1 H), 6.95 (d, 1 H, J = 8.8 Hz), 7.23 (d, 1 H, J = 8.8 Hz), 9.45 (s, 1 H), 11.38 (s, OH); 13C NMR (CDCl3): δ 19.9, 56.0, 61.0, 61.1, 109.7, 116.2, 117.8, 118.2, 125.2, 125.5, 133.9, 140.5, 145.6, 152.2, 154.4, 156.9, 196.7.
5-Hydroxy-3,6,7-trimethoxy-1,4-phenanthrenequinone; calanquinone A (6a)
Compound 5a (115.02 mg, 0.35 mmol) was dissolved in 5 mL <span class="Chemical">CH2Cl2 and iodotrimethylsilane (112.23 mg, 0.56 mmol) was added to the solution in portions. The mixture was stirred at 60°C overnight (TLC monitoring) and then MeOH was added to quench the reaction. The solvent was evaporated and the residue extracted with Et2O/H2O. The Et2O layer was dried over Na2SO4. The mixture was chromatographed on silica gel and eluted with EtOAc/n-hexane (1∶2) to yield 6a (15.72 mg, 14.28%). Black solid; mp 187°C; 1H NMR (CDCl3): δ 3.96 (s, 3 H), 4.01 (s, 3 H), 4.02 (s, 3 H), 6.14 (s, 1 H), 6.85 (s, 1 H), 8.03 (d, 1 H, J = 8.4 Hz), 8.08 (d, 1 H, J = 8.4 Hz), 10.73 (s, OH); 13C NMR (CDCl3): δ 55.9, 56.9, 60.8, 101.2, 107.2, 118.5, 121.7, 128.1, 132.8, 134.8, 136.9, 140.1, 148.1, 155.0, 161.5, 184.4, 186.0; HRESIMS m/z 313.0648 (calculated for C17H14O6 -H, 313.0712).
Authors: Jeong Eun Gong; Ji Eun Kim; Su Jin Lee; Yun Ju Choi; You Jeong Jin; Young Whan Choi; Sun Il Choi; Dae Youn Hwang Journal: Molecules Date: 2022-01-23 Impact factor: 4.411
Authors: Alexey A Lagunin; Varvara I Dubovskaja; Anastasia V Rudik; Pavel V Pogodin; Dmitry S Druzhilovskiy; Tatyana A Gloriozova; Dmitry A Filimonov; Narahari G Sastry; Vladimir V Poroikov Journal: PLoS One Date: 2018-01-25 Impact factor: 3.240