Next Article in Journal
NKG2D Polymorphism in Melanoma Patients from Southeastern Spain
Next Article in Special Issue
The d16HER2 Splice Variant: A Friend or Foe of HER2-Positive Cancers?
Previous Article in Journal
Pheochromocytomas and Paragangliomas: From Genetic Diversity to Targeted Therapies
Previous Article in Special Issue
Newer-Generation EGFR Inhibitors in Lung Cancer: How Are They Best Used?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 11
Issue 4
10.3390/cancers11040437
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in Non-Small Cell Lung Cancers: Combined Experimental and Theoretical Investigations

by
Panupong Mahalapbutr
1,
Piyanuch Wonganan
2,*,
Warinthorn Chavasiri
3 and
Thanyada Rungrotmongkol
1,4,*
1
Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok 10330, Thailand
2
Department of Pharmacology, Faculty of Medicine, Chulalongkorn University, Pathumwan, Bangkok 10330, Thailand
3
Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok 10330, Thailand
4
Ph.D. Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok 10330, Thailand
*
Authors to whom correspondence should be addressed.
Cancers 2019, 11(4), 437; https://doi.org/10.3390/cancers11040437
Submission received: 11 January 2019 / Revised: 22 March 2019 / Accepted: 25 March 2019 / Published: 28 March 2019
(This article belongs to the Special Issue Epidermal Growth Factor Receptor Signaling in Cancer)
Browse Figures
Versions Notes

Abstract

:
Epidermal growth factor receptor (EGFR) is the key molecular target for non-small cell lung cancer (NSCLC) due to its major contribution to complex signaling cascades modulating the survival of cancer cells. Targeting EGFR-mediated signaling pathways has been proved as a potential strategy for NSCLC treatment. In the present study, mansonone G (MG), a naturally occurring quinone-containing compound, and its semi-synthetic ether derivatives were subjected to investigate the anticancer effects on human NSCLC cell lines expressing wild-type EGFR (A549) and mutant EGFR (H1975). In vitro cytotoxicity screening results demonstrated that butoxy MG (MG3) exhibits the potent cytotoxic effect on both A549 (IC50 of 8.54 μM) and H1975 (IC50 of 4.21 μM) NSCLC cell lines with low toxicity against PCS201-010 normal fibroblast cells (IC50 of 21.16 μM). Western blotting and flow cytometric analyses revealed that MG3 induces a caspase-dependent apoptosis mechanism through: (i) inhibition of p-STAT3 and p-Akt without affecting upstream p-EGFR and (ii) activation of p-Erk. The 500-ns molecular dynamics simulations and the molecular mechanics combined with generalized Born surface area (MM/GBSA)-based binding free energy calculations suggested that MG3 could possibly interact with STAT3 SH2 domain and ATP-binding pocket of Akt. According to principal component analysis, the binding of MG3 toward STAT3 and Akt dramatically altered the conformation of proteins, especially the residues in the active site, stabilizing MG3 mainly through van der Waals interactions.

Graphical Abstract

1. Introduction

Lung cancer is the first leading cause of cancer-related mortality worldwide [1,2,3]. Approximately 85% of diagnosed lung cancer cases are classified as non-small cell lung cancer (NSCLC), for which the 5-year survival rate is only 17.8% [4]. The genetic alterations of epidermal growth factor receptor (EGFR) such as overexpression, exon 19 deletion, and exon 21 L858R substitution have been characterized as oncogenic drivers for NSCLC development [5,6,7]. EGFR (also known as HER1 or ErbB1), a receptor tyrosine kinase, is commonly overexpressed in several types of cancer, including lung carcinoma [8]. The activation of EGFR-mediated signaling pathways is triggered by the binding of growth factors, including EGF and transforming growth factor alpha (TGFα), to the extracellular portion of EGFR, which subsequently induces receptor dimerization and cross-phosphorylation of specific tyrosine residues located on the cytoplasmic tyrosine kinase (TK) domain. These structural modifications result in stimulating a vast array of downstream signaling cascades, e.g., mitogen-activated protein kinase (MAPK), phosphoinositide 3-kinase (PI3K)/Akt, and signal transducer and activator of transcription (STAT) pathways, leading to cell growth, proliferation, migration, and apoptosis evasion [8,9].
Cisplatin (CDDP), a platinum-based chemotherapeutic drug, has been commonly used to treat a wide range of solid malignancies, including lung cancer [10,11]. CDDP is the standard regimen in the first-line chemotherapy for patients with advanced stage NSCLC [10,12,13,14], especially patients carrying wild-type EGFR [15,16]. After cellular uptake, CDDP becomes a positively charged aquo complex that can interact with deoxyribonucleic acid (DNA) to form intra- and inter-strand cross-links, resulting in apoptosis induction [10,17]. However, the efficacy of CDDP-based chemotherapy is limited by numerous severe side effects [18], as well as acquired drug resistance [19,20,21,22].
The first-generation tyrosine kinase inhibitors (TKIs, e.g., erlotinib and gefitinib) have shown to significantly prolong the progression-free survival of NSCLC patients harboring EGFR mutations, primarily exon 19 deletion and exon 21 L858R substitution mutations [23,24]. TKIs compete with adenosine triphosphate (ATP) at the ATP-binding site of the receptor, inhibiting EGFR-mediated signal transduction [25]. However, acquired resistance caused by the secondary mutation T790M develops inevitably after a median response duration of 9 to 13 months [26,27,28]. The replacement of threonine (T) to methionine (M) causes a steric hindrance inside ATP-binding pocket and alters the conformation of TK domain, resulting in increasing its affinity for ATP substrate while decreasing the binding affinity for TKIs [28,29]. Given that NSCLC cells rapidly acquire resistance to both CDDP and TKIs; thus, there is an urgent need to search for a novel compound that can potentially overcome such problems by targeting alternative intracellular survival signaling pathways in NSCLC.
Mansonone G (MG, Figure 1A), a 1,2-naphthoquinone-containing compound, is the major product isolated from the heartwood extract of Mansonia gagei Drumm. from the Sterculiaceae family [30]. MG demonstrates various biological activities, including antitumor [31], antibacterial [32], antiestrogenic [33], anticholinesterase [34], and antifungal activities [32]. Recently, semi-synthetic ether derivatives of MG (Figure 1A) have been shown to exhibit higher antibacterial activity against Staphylococcus aureus [30] and inhibit adipocyte differentiation and lipid accumulation [35] more than the MG parent compound. Although several pharmacological effects of MGs have been reported, the anticancer activity of MG and its derivatives against human NSCLC remains largely unknown. Therefore, in the present study, we aimed to search for the most potent cytotoxic MG analog against human NSCLC cell lines expressing wild-type EGFR (A549) and L858R/T790M EGFR (H1975). Moreover, the mechanisms underlying cell death were experimentally and theoretically investigated.

2. Results

2.1. In Vitro Cytotoxicity Screening of MG Derivatives against NSCLC Cell Lines

Initially, we screened for the most potent cytotoxic MG analogs against NSCLC cell lines expressing wild-type EGFR (A549) and L858R/T790M EGFR (H1975) using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. MG derivatives displaying a percentage of cell viability at 10 μM (%CV10 µM) <50 were defined as potent compounds.
As shown in Table 1, MG4 exhibited the most potent cytotoxicity against A549 cells (%CV10 µM was 7.93 ± 0.43) followed by MG3 (%CV10 µM was 46.90 ± 1.21), indicating that these two ether analogs were the most potent compounds toward wild-type cells. In the H1975 mutant cell line, all semi-synthesized MGs (except MG9) exhibited high cytotoxic activity with %CV10 µM of <50, indicating that T790M-positve NSCLC cell line was more susceptible to MG derivatives than A549 cells. Remarkably, MG3 and MG4, exhibiting strong cytotoxicity toward A549 cells, also showed great cytotoxic effects on H1975 cells with the %CV10 µM of 17.04 ± 1.42 and 23.52 ± 1.44, respectively.
Taken together, ether analogs of MG were more cytotoxic to H1975 cells than A549 cells. The two semi-synthetic ether analogs MG3 and MG4, which exhibited high cytotoxicity toward both NSCLC cell lines, were selected for further studies.

2.2. Butoxy Mansonone G Exhibits a Potent Cytotoxicity against NSCLC Cells

The half maximal inhibitory concentration (IC50) of the two focused MG analogs, MG3 and MG4, and CDDP, the most commonly used chemotherapeutic agent for NSCLC patients, were evaluated using MTT assay on two NSCLC cell lines A549 and H1975. As shown in Figure 2A–C, all three tested compounds decreased the viability of both A549 and H1975 cells in a concentration-dependent manner. The IC50 of CDDP against two NSCLC cells obtained from this study correlated well with the previous reports [11,36]. Notably, the two selected MG derivatives were more toxic than CDDP toward cancer cells, as evidenced by the lower IC50 of (i) ~4-fold and ~9-fold toward A549 cells and (ii) ~8-fold and ~11-fold against H1975 cells for MG3 and MG4, respectively (Figure 2G).
We further evaluated the toxicity of our focused MGs toward PCS201-010 normal skin fibroblast cell line (Figure 2D–F). The obtained IC50 results indicated that MG4 and CDDP were toxic to normal cells, as clearly shown by a similarity of IC50 between cancer and normal fibroblast cells (selectivity index (SI) of ~1, Figure 2H). Intriguingly, MG3 was relatively safe to normal cells, in which the IC50 observed in PCS202-010 cells was higher than those of A549 (SI of 2.48) and H1975 (SI of 5.03). Since MG3, possessing potent cytotoxicity against cancer cells, was less toxic to normal cells than MG4, we then elucidated the mechanisms underlying MG3-induced cell death in NSCLC cell lines experimentally and theoretically.
According to the morphological alteration of A549 and H1975 NSCLC cells upon MG3 treatment for 48 h (Figure 2I), it can be clearly seen that MG3 dose-dependently induced cellular shrinking, a predominant characteristic of programmed cell death [37], suggesting that MG3 promoted cell death through an apoptosis-inducing effect.

2.3. Butoxy Mansonone G Induces Apoptosis in A549 and H1975 Cell Lines

To determine whether apoptotic mechanism is involved in MG3-induced cytotoxicity in NSCLC cells, flow cytometric analysis of Annexin V/PI stained cells was carried out. Following 24 h of treatment, MG3 at 16 µM significantly induced apoptotic cell death in A549 cells, whereas MG3 at 2 µM induced significant apoptosis in H1975 cells (Figure 3A,B), indicating that apoptosis-inducing effect of MG3 was more pronounced in NSCLC cells carrying mutant EGFR than in NSCLC cells carrying wild-type EGFR. Additionally, MG3 at 8 μM dramatically induced apoptosis in H1975 cells (~50%) which was higher than CDDP, reflecting a strong cytotoxic activity. It should also be noted that, when compared to 30 μM CDDP, the lower concentrations of MG3 at: (i) 16 μM (~2-fold lower) for A549 cells and (ii) 2 μM (15-fold lower) for H1975 cells can significantly trigger cell apoptosis.
To further confirm the apoptosis-inducing effect of MG3 on NSCLC cells, the cleavage of procaspase-3 and poly(ADP-ribose) polymerase (PARP), key hallmarks of apoptosis, was determined using western blotting. Note that for H1975 cell line, MG3 at 8 μM was highly toxic to the cells (as evidenced by flow cytometric analysis), leading to a low concentration of extracted proteins; and thus, this concentration was excluded from this study. As shown in Figure 3C,D, MG3 (16 μM for A549 and 2 μM for H1975) as well as 30 μM CDDP significantly induced the cleavage of procaspase-3 and PARP, which was in good agreement with a significant apoptotic cell death detected by flow cytometric analysis.
We next characterized whether caspase-3 activation (Figure 3C,D) is mandatory for MG3-induced apoptosis. NSCLC cells were pretreated with Z-Val-Ala-Asp-(OMe)-fluoromethylketone (Z-VAD(OMe)-FMK), an irreversible pan-caspase inhibitor, for 1 h prior to challenge with MG3 for 24 h. As shown in Figure 3E,F, both MG3 and CDDP decreased cell viability by ~40% in both A549 and H1975 cells, and Z-VAD(OMe)-FMK alone did not affect the cell viability of cancer cells (%CV of ~100). Intriguingly, blockage of caspase activation by Z-VAD(OMe)-FMK inhibitor significantly restored cell viability in A549 and H1975 cells for both MG3- and CDDP-treated groups. These findings clearly demonstrated that activation of caspase-3 enzyme plays a crucial role in MG3-induced cell apoptosis.

2.4. Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in NSCLC Cell Lines

To elucidate the effect of MG3 on EGFR-mediated survival signaling pathways, western blot analysis was performed. As shown in Figure 4A,B, MG3 and CDDP significantly inhibited the phosphorylation of STAT3 and Akt in a concentration-dependent manner in both A549 and H1975 cells. Conversely, the expression of p-Erk was significantly increased after treatment with such two compounds. Remarkably, MG3, although at lower concentrations (~2-fold lower for A549 and 15-fold lower for H1975), exhibited similar effects on EGFR-mediated survival signaling pathways as 30 µM CDDP.
We further investigated whether the downregulation of p-STAT3 and p-Akt caused by MG3 was mediated through the inhibition of p-EGFR. The data in Figure 4C,D demonstrate that the addition of EGF dramatically increased p-EGFR in A549 vehicle-treated cells. On the other hand, EGF has no effect on phosphorylation of EGFR in H1975 control cells, since T790M mutation enhances the ATP binding affinity without EGF binding [38]. The 10 µM erlotinib, an EGFR TK inhibitor, totally inhibited the phosphorylation of EGFR in A549 cells, which was in agreement with a previous study [39]. However, it should be noted that pre-treatment with MG3 and CDDP did not alter the expression levels of p-EGFR following EGF stimulation, indicating that EGFR was not the preferential binding site for both MG3 and CDDP.
Taken together, the proposed underlying mechanisms of MG3 against two NSCLC cell lines are illustrated in Figure 4E showing that MG3 inhibited downstream activity of STAT3 and Akt without interfering phosphorylation of EGFR, whilst the phosphorylation of Erk was significantly enhanced upon MG3 treatments. These signal transduction effects led to the activation of caspase-3, cleavage of PARP, and induction of apoptosis, respectively.
Since we found that MG3 inhibits p-STAT3 and p-Akt, we further elucidated the atomistic binding mechanisms of MG3 against such target proteins, which culminated in phosphorylation inhibition, using multiple computational modeling techniques. The structural and dynamics properties, ligand-protein interactions, and binding affinity of MG3 in complex with two signaling proteins were compared to the known inhibitors as well as the apo-proteins.

2.5. Predictive Binding Affinity of Butoxy Mansonone G against STAT3 and Akt Signaling Proteins

To estimate the binding affinity of MG3 against the focused proteins STAT3 and Akt in comparison with their known inhibitors, the molecular mechanics combined with generalized Born surface area (MM/GBSA) method was applied on the 200 molecular dynamics (MD) snapshots extracted from the last 200-ns simulations. The binding free energy (ΔGbind) together with its energy components are summarized in Table 2. In the case of STAT3, the calculated ΔGbind results were ranked in the order of MG3 (−8.54 kcal/mol) << cryptotanshinone (CTS, −5.09 kcal/mol) < S3I201 (−3.73 kcal/mol), suggesting that the susceptibility of MG3 was significantly higher than those of known inhibitors. By considering Akt models, the ΔGbind of MG3 (−9.19 kcal/mol) was in the range of uprosertib (−10.45 kcal/mol) and H8 (−9.68 kcal/mol) inhibitors. Due to the nonpolar structure of MG3 and CTS (Figure 1A,B), the molecular mechanics energy (ΔEMM) revealed that van der Waals interaction (ΔEvdW) was the main force driving protein-ligand complexation (ΔEvdW of −35.77, −39.14, and −35.61 kcal/mol for MG3/STAT3, MG3/Akt, and CTS/STAT3, respectively). In contrast, the electrostatic attraction (ΔEele) was found to mainly contribute toward S3I201/STAT3 (ΔEele of −110.36 kcal/mol), uprosertib/Akt (ΔEele of −146.76 kcal/mol), and H8/Akt (ΔEele of −168.13 kcal/mol) complexes, since these inhibitors contain the ionic charged moiety in the chemical structures (R-COO−, R-NH3+, and R-NH+-R groups, respectively, Figure 1B).
Notably, the obtained ΔGbind of all studied inhibitors agreed well with the trend of experimental ΔGbind (ΔGbind, exp), implying that our free energy calculations could successfully predict the binding affinity of protein/inhibitor complexes. Taken together, the experimental and theoretical results suggested that the anticancer activity of MG3 was due to binding to STAT3 and Akt signaling proteins.

2.6. Key Binding Residues

The per-residue decomposition free energy ( Δ G bind residue ) calculation based on the MM/GBSA method was used to investigate the crucial amino acid residues involved in ligand binding within the SH2 domain of STAT3 and the ATP-binding pocket of Akt. The total contributing amino acids of all complexes are shown in Figure 5, where the negative and positive Δ G bind residue values represent respectively the stabilization and destabilization energies of the considered residue.
In the case of STAT3, there are three subpockets in the SH2 domain, including (i) pY + 0 (residues 591 and 609–620), (ii) pY − X (residues 592–608), and (iii) pY + 1 (residues 621–639) pockets. The pY + 0 site contains several polar residues responsible for phosphotyrosine (pTyr) binding, while the two subsites pY − X and pY + 1 are the hydrophobic regions (Figure 5A). Note that among residues 458–722 of STAT3 model, only the contribution from the residues 540–660 is shown. The obtained results demonstrated that there were four and eight amino acids involved in the binding of the two STAT3 inhibitors CTS (e.g., I589, E594, L598, and I634) and S3I201 (e.g., K557, I589, E594, I597, L607, R609, I634, and Q635), respectively; whereas MG3 interacted with the residues L598, T632, and I634 inside the hydrophobic subsites pY − X and pY + 1. Notably, the contribution from residues L598 and I634 toward MG3/STAT3 complex matched to the known inhibitors. The lipophilic group of ligands, including the cyclohexane ring of CTS, the O-tosyl group of S3I201, and the alkyl side chain of MG3 were found to be encapsulated into the hydrophobic pY + 1/pY − X pocket of STAT3, while the polar moieties of S3I201 orientated in the pY + 0 site, forming hydrogen bond (H-bond) interactions (Figure S1).
For Akt signaling protein (Figure 5B), the amino acids establishing inhibitor binding were: (i) V164, E234, E278, M281, T291, F438, and F442 for uprosertib and (ii) V164, Y229, A230, M281, and T291 for H8. Notably, the key binding residues V164, E234, M281, and F438 involved in MG3 binding were identical to those of the Akt inhibitors. Interestingly, the ligand binding mode of Akt models shared a structurally-related characteristic, in which the aromatic moiety of all studied compounds approached the key residue M281.

2.7. In Silico Study on Conformational Change of STAT3 and Akt upon Butoxy Mansonone G Binding

The structurally relevant motions of STAT3 and Akt signaling proteins derived from MG3 recognition were investigated in comparison with the apo-protein using principal component analysis (PCA) on the 2000 MD snapshots taken from the last 200-ns simulations. The results are illustrated in Figure 6, where the arrow and its length indicate the direction and amplitude of motions, respectively. Note that among residues 458–722 of STAT3 model, only the protein motion from the residues 499–688 is shown.
The first 15 PC modes showed the % accumulated variance of (i) 72.99 and 93.68 for apo and holo forms of STAT3 and (ii) 55.21 and 57.95 for apo and holo forms of Akt, respectively. The percentage of variances for PC1 of all systems was much higher than that of PC2, indicating that this mode could represent the significant motions of proteins.
By considering STAT3 model (Figure 6A), the first principal component (PC1) showed that MG3 binding importantly converted the direction of motions of overall protein to approach the ligand in a different manner from that of the apo form. Remarkably, the residues 592–601 on pY − X (green) and 625–633 on pY + 1 (magenta) hydrophobic pockets of STAT3 SH2 domain displayed not only increased direction of motions, but also enhanced amplitude of motions upon MG3 binding mainly through vdW interaction (Table 2). For the Akt systems (Figure 6B), the residues 157–163 on glycine-rich loop (GRL, green) inside the ATP-binding pocket of apo form pointed outward from the binding site, making it an opened conformation. Intriguingly, a complexation with MG3 led to the conversion of direction of GRL motion to be located closer to the MG3 molecule, resulting in a closed conformation.

3. Discussion

Platinum-based chemotherapy and the first-generation TKIs have been used as the first-line treatment for NSCLC patients carrying wild-type and mutant EGFRs, respectively [15,16,23,24]. However, acquired drug resistance is inevitable after a progression-free period of approximately 9 to 13 months. Therefore, a novel anticancer compound that remains effective in both NSCLC cells expressing wild-type and T790M-positive EGFRs is critically needed. MG, a naphthoquinone-containing compound extracted from Mansonia gagei Drumm, was shown to exhibit anticancer activity toward the A2780 ovarian cancer cell line with an IC50 of 10.2 ± 0.9 μM [31]. Recently, the etherification of the hydroxyl group of MG has been reported to potentially result in antibacterial activity as well as to suppress adipocyte differentiation and lipid accumulation by more than natural MG [30,35]. In the present study, we experimentally and theoretically elucidated the cytotoxic activity of MG and its semi-synthetic derivatives against A549 (expressing wild-type EGFR) and H1975 (expressing L858R/T790M EGFR) NSCLC cell lines. We found that, among ten ether analogs, MG3 and MG4 displayed the most potent cytotoxicity toward both A549 and H1975 cells (Table 1). Notably, the IC50 values of such two MGs were much lower than IC50 of CDDP for both NSCLC cell lines (Figure 2G). Introduction of longer carbon side chain to MG makes its chemical structure more hydrophobic, which culminate in increasing cellular uptake [30,44]. Therefore, it is likely that chemical modification of MG via increasing number of carbon units of alkyl side chain enhances the cytotoxicity of MG1-MG4 against NSCLC cell lines. Conversely, MG5 containing twelve carbons showed a dramatic reduction of cytotoxic activity, which may be due to the cut-off effect [30,45].
Unfortunately, we found that MG4 was highly toxic toward the PCS201-010 normal skin fibroblast cell line (IC50 of 3.73 ± 0.23 μM), thus, this compound was excluded from our mechanistic studies. Our present findings are consistent with others showing that MG3 was less toxic to differentiated adipocytes than MG4 [35]. However, it should be noted that only one normal cell line was tested in this study; thus, further investigations on other normal cell lines as well as in vivo studies should be conducted. Interestingly, as compared to the commercial EGFR-targeted drugs, the SI of gefitinib (SI = 0.81) and osimertinib (SI = 1.15) toward A549 NSCLC and normal human bronchial epithelial (HBE) cell lines is much lower than MG3 (2.48, Figure 2H) [46]. In addition, the SI of chemotherapeutic drug methotrexate toward A549 cells and normal embryo fibroblast (NIH/3T3) cells was only 0.02, reflecting a strong toxicity [47,48]. Altogether, MG3, which showed a lower toxicity (IC50 of 21.16 ± 0.98 μM) against PCS201-010 cells, was then selected for further studying the mechanisms underlying cell death.
Activation of caspase enzymes leads to the generation of signaling cascades responsible for apoptotic events [49]. In the present study, we documented that MG3 induced caspase-dependent apoptosis in both A549 and H1975 cells, in which H1975 cells expressing mutant EGFR were more susceptible to apoptosis-inducing effect of MG3 rather than A549 cells expressing wild-type EGFR. The reason for this observation is that L858R/T790M EGFR stimulates cell growth through the dramatic enhancement of catalytic phosphorylating activity over wild-type EGFR [50], which can provide more preferential target signaling proteins for MG3, as evidenced by the significant downregulation of p-STAT3 and p-Akt at low concentrations of 1 μM and 2 μM, respectively, in H1975 cells (Figure 4B). Previous studies demonstrated that CDDP induced apoptosis through the activation of caspase-3 and PARP [51,52,53,54]. Similarly, this study revealed that CDDP promoted apoptotic cell death through the activation of caspase cascades (Figure 3C–F).
Many lines of evidence have shown that several naphthoquinone-containing compounds, such as shikonin, plumbagin, furano-1,2-naphthoquinone, ramentaceone, mansonone E, and CTS significantly inhibited STAT3 and Akt signaling pathways in various kinds of malignancies [40,55,56,57,58,59,60,61]. In agreement with these reports, data in Figure 4A,B show that our ortho-naphthoquinone MG3 concentration-dependently inhibited phosphorylation of STAT3 and Akt in both A549 and H1975 cell lines, which strongly correlated with MM/GBSA free energy calculations showing that MG3 could possibly interact with SH2 domain of STAT3 and ATP-binding pocket of Akt in a similar manner to that of known STAT3 and Akt inhibitors (Table 2). Remarkably, our structural analyses on both STAT3 and Akt inhibitors were consistent with other experimentally/theoretically derived data as follows. The binding orientations ((i) O-tosyl group of S3I201 occupying the pY − X pocket and (ii) cyclohexane group of CTS pointing toward the pY+1 hydrophobic region), key binding amino acid residues (e.g., I597, R609, and Q635 for S3I201 and I634 for CTS), and H-bond formation patterns of STAT3 inhibitors displayed a similar manner to several research works [40,62,63,64,65]. It was evidenced that M281 is the key binding residue for hydrophobic packing against indole/phenyl ring of Akt inhibitors through ε-CH···π interaction [66,67]. In correlation with this fact, our Δ G bind residue data demonstrated that the aromatic moiety of all studied ligands pointed toward M281 residue (Figure 5B). Since the α,β-unsaturated carbonyl (α,β-UC) compounds can covalently interact with biological thiols of several kinase proteins via hetero-Michael addition reaction [68], we then measured the distance between the center of mass of α,β-UC unit of MG3 and the thiol group (SH-) of cysteine residues using the final snapshot from the 500-ns MD simulation. The obtained results revealed that α,β-UC part of MG3 positioned far away (>15 Å) from the cysteine residues in Akt’s active site (Figure S2), implying that MG3 could not form the covalent adduct with Akt. Altogether, MG3 showed a somewhat similar binding pattern to STAT3 and Akt inhibitors upon molecular complexation, suggesting that MG3 could likely inhibit the phosphorylation of such proteins in a similar fashion to that of the known inhibitors. By conducting PCA, we discovered that MG3 importantly induced large conformational changes of STAT3 and Akt, especially in ligand-binding pocket, as strongly evidenced by the superimposed X-ray crystal structures between apo (PDB ID: 1GZN [69]) and holo forms (PDB ID: 4GV1 [70]) of Akt (Figure S3). This might explain the atomistic mechanisms underlying the inhibition of phosphorylation-induced activation of such signaling protein mediated by MG3, which culminate in cell apoptosis. It should be noted that in silico results from MD simulations were suggestive; thus, further experimental techniques (e.g., isothermal titration calorimetry (ITC), circular dichroism (CD), and surface plasmon resonance (SPR)) should be conducted to confirm our findings.
Although our results showed that MG3 inhibited STAT3 and Akt activities, MG3 did not interfere with the phosphorylation of EGFR, indicating that MG3 preferentially targeted EGFR’s downstream signaling molecules STAT3 and Akt rather than upstream EGFR. We also discovered that treatment with CDDP did not cause any significant changes in the expression of p-EGFR, which was in good agreement with the fact that DNA is a molecular target for CDDP [17,71], and the changes of signaling cascades are derived from platinum-DNA adduct [71]. Because PARP cleavage was more pronounced in cells treated with CDDP than cells treated with MG3, it is possible that NSCLC cells induced PARP in order to repair damaged DNA caused by CDDP-mediated platination. In contrast to CDDP, DNA was found to be the non-preferential binding site for MG3, as evidenced by: (i) the lower level of cleaved PARP than CDDP (Figure 3C,D), (ii) the higher CDOCKER interaction energy (−32.22 kcal/mol) than CDDP (−44.24 kcal/mol), and (iii) the high distance between the center of mass (Cm) of MG3 and DNA (d(m(MG3)-Cm(DNA))) obtained from three independent 100-ns MD simulations (Figure S4). Altogether, MG3 exhibited overall mechanisms of action similar to those of CDDP, but preferred targeting proteins rather than DNA, which can serve as a promising anticancer agent for NSCLC patients harboring CDDP resistance.
In contrast to STAT3 and Akt, we noticed that treatment with MG3 and CDDP induced phosphorylation of Erk in both A549 and H1975 cells. A significant upregulation of p-Erk was found to be correlated with the generation of cleaved caspase-3 and cleaved PARP (Figure 3C,D). These findings were consistent with previous studies demonstrating that high level of p-Erk can promote cell apoptosis through the activation of caspase-3 [72,73]. Erk plays a dual role in both cell proliferation and cell death [72], and activation of Erk is extremely important for CDDP-induced apoptosis [71,73,74]. Thus, it is possible that, in addition to inhibition of STAT3 and Akt, cytotoxicity of MG3 may be mediated through activation of Erk in NSCLC cells.
Altogether, our present study provided the first step of the underlying mechanisms of MG3 toward EGFR-mediated signaling pathways in NSCLC cell lines expressing wild-type and mutant EGFRs. However, further investigations on other signaling pathways as well as the kinase screening assays need to be performed in order to gain more insights into the signal transduction inhibitions caused by MG3.

4. Materials and Methods

4.1. Experimental Part

4.1.1. Chemical Reagents and Antibodies

MG was extracted from the heartwood of M. gagei, whereas MG ether derivatives were semi-synthesized according to the previous study [30]. Bovine serum albumin (BSA), dimethyl sulfoxide (DMSO), CDDP, MTT, and protease inhibitor were purchased from Sigma-Aldrich (St. Louis, MO, USA). RIPA lysis buffer was purchased from Thermo Fisher Scientific (Waltham, MA, USA). The protein assay reagents were purchased from Bio-Rad (Hercules, CA, USA). Human epidermal growth factor (EGF, 8916) and antibodies against phospho-EGFR (p-EGFR, 2234), total-EGFR (t-EGFR, 4267), phospho-STAT3 (p-STAT3, 9145), total-STAT3 (t-STAT3, 12640), phospho-Akt (p-Akt, 4060), total-Akt (t-Akt, 4691), phospho-Erk (p-Erk, 4377), total-Erk (t-Erk, 4695), caspase-3 (9662), PARP (9542), GAPDH (5174), and anti-rabbit IgG HRP-linked antibody (7074) were purchased from Cell Signaling Technology (Santa Cruz, CA, USA). The pan-caspase inhibitor Z-VAD(OMe)-FMK (ab120487) was purchased from Abcam (Cambridge, UK).

4.1.2. Cell Lines and Culture

Human NSCLC cell lines A549 and H1975 as well as human normal skin fibroblast cell line (PCS201-010) were purchased from American Type Culture Collection (ATCC, Manassas, VA, USA). A549 cells were grown in Dulbecco’s modified Eagle’s minimal essential medium (DMEM; Gibco, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS; Gibco), 100 U/mL penicillin, and 100 µg/mL streptomycin (Gibco). H1975 cells were cultured in RPMI-1640 medium containing 10% FBS, 100 U/mL penicillin, and 100 µg/mL streptomycin. The DMEM with high glucose (4500 mg/L) supplemented with 10% FBS, 100 U/mL penicillin, and 100 µg/mL streptomycin was used for culturing PCS201-010 cells. All cells were maintained at 37 °C in a humidified 5% CO2 atmosphere.

4.1.3. Cell Viability Assay

Cell viability was assessed using the MTT assay. Cells were seeded into 96-well plates at a density of 5 × 103 cells/well for H1975 and PCS201-010 as well as of 3 × 103 cells/well for A549. After overnight incubation, cells were treated with MGs at 10 and 100 μM for 48 h. Note that, due to the low solubility of MG4 and MG10, the highest prepared concentration was 50 μM. Subsequently, the MTT solution (5 mg/mL) was added and then incubated for 4 h. The medium was removed and 150 µL of DMSO was added to each well. Finally, the absorbance of formazan product was measured at a wavelength of 570 nm using a LabSystems Multiskan MS microplate reader (Thermo Scientific, Vantaa, Finland). The selectivity index (SI) was calculated according to the following equation: SI = IC50 for normal cells/IC50 for cancer cells.

4.1.4. Western Blotting

A549 and H1975 cells were seeded into a 6-well plate at a density of 2 × 105 cells/well and 3 × 105 cells/well, respectively. After overnight incubation, cells were treated with indicated compounds. Note that the concentration of MG3 was varied to two-fold, one-half, and one-fourth of its IC50, whereas the positive control CDDP at the IC50 of 30 μM was used. After 24 h of incubation, cells were rinsed twice with cold PBS, homogenized in RIPA buffer containing protease inhibitor, and incubated on ice for 45 min. Total protein (20 µg) were separated on 8% SDS-PAGE and subsequently transferred to a PVDF membrane. The membrane was blocked with 3% non-fat dry milk for 1 h and then incubated with primary antibody at 4 °C overnight. After incubation, the membrane was washed thrice with TBST buffer (5 min each) and incubated with HRP-linked secondary antibody for 2 h at room temperature. Immunoreactive bands were detected using HRP substrate (Millipore, Billerica, MA, USA) and quantitatively measured using Image Studio Lite software (LI-COR, Lincoln, NE, USA). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as internal control for protein normalization.
Note that for detecting p-EGFR, it was reported that the expression of p-EGFR could not be clearly detected at 24 h due to the short half-life of activated EGFR (~1.5–4 h) [75]; thus, we pre-incubated NSCLC cell lines with the indicated concentrations of MG3 and CDDP in serum free media for 1 h prior to stimulation of EGFR with EGF (50 ng/mL) for 10 min.

4.1.5. Flow Cytometric Evaluation of Apoptosis

NSCLC cells were plated on a 6-well plate at a density of 2 × 105 cells/well for A549 cells and 3 × 105 cells/well for H1975 cells. After overnight incubation, cells were treated with 1, 2, 4, 8, and 16 μM MG3 and 30 μM CDDP for 24 h. After treatment, cells were harvested by trypsinization and collected by centrifugation at 1500 rpm. Subsequently, cells were washed with cold PBS and stained with 3 μL of Annexin-V fluorescein dye and 1 μL of propidium iodide (PI) at room temperature in the dark for 20 min. After that, cells were resuspended in 400 μL of cold assay buffer containing 0.01 M HEPES, 2.8 mM CaCl2, and 125 mM NaCl. The percentage of apoptotic cells was quantitatively measured using BD FACSCalibur flow cytometer (BD Bioscience, Heidelberg, Germany).

4.1.6. Statistical Analysis

The quantitative data are expressed as mean ± standard error of mean (SEM) of triplicate experiments. Differences between groups were determined using one-way analysis of variance (ANOVA) followed by a Turkey post hoc test. Differences were considered to be significant at p ≤ 0.05.

4.2. Computational Part

4.2.1. Preparation of Initial Structures

The crystal structures of human STAT3 (PDB ID: 1BG1) [76] and Akt1 (PDB ID: 4GV1) [70] were obtained from Protein Data Bank (PDB). The missing amino acid residues were completed using SWISS-MODEL server [77]. The 3D structure of MG3 was obtained from a previous study [78], whereas the known inhibitors of STAT3 (CTS and S3I201) and Akt1 (uprosertib and H8) (Figure 1B) were built and subsequently optimized by the HF/6-31(d) level of theory using the Gaussian09 program [79]. The protein-ligand complexes were generated using CDOCKER module implemented in Accelrys Discovery Studio 2.5 (Accelrys Inc.) [80] with 100 docking runs. Note that for STAT3 SH2 domain, prior to perform docking, the protein was relaxed in aqueous solution by conducting a short MD simulation at 298.0 K for 100 ps (as detailed in the next section). The residues K591, R595, R609, E612, W623, and Q635 of STAT3 were defined as binding site with a docking sphere radius of 15 Å, whereas the co-crystalized inhibitor at ATP-binding pocket was used as the docking center for Akt. In addition, the MG3/DNA complex was also tested and simulated. Totally, there are nine simulated models in which the computational details of all system preparations are summarized in Table S1.
The protonation states of all ionizable amino acids were characterized using PROPKA 3.0 [81] at pH 7.0. The electrostatic potential (ESP) charges of ligand were computed at the HF/6-31(d) level of theory, whereas the restrained ESP (RESP) charges and corresponding parameters of ligands were generated respectively using antechamber and parmchk modules in AMBER16 according to previous studies [82,83,84]. The AMBER ff14SB force field [85] was applied for protein, whilst the ligand was treated using the general AMBER force field (GAFF) [86,87,88]. The missing hydrogen atoms were added using the LEaP module. The added hydrogen atoms were then minimized using 1000 steps of the steepest descents (SD) and 2500 steps of conjugated gradient (CG) approaches. Subsequently, each system was solvated using TIP3P water model [89] in truncated octahedron periodic box with the minimum distance of 10 Å from the system surface. The systems were neutralized using Cl− or Na+ counter ions. The minimization with the SD of 1000 steps and CG of 2500 steps was performed on the added water molecules and counter ions, and finally the entire system was wholly minimized using the same procedure.

4.2.2. Molecular Dynamics (MD) Simulations and Binding Free Energy Calculations

MD simulations of all studied complexes were performed under periodic boundary condition using AMBER 16 program. The short-range cutoff of 10 Å was employed for non-bonded interactions, whilst the long-range electrostatic interactions were treated using Particle Mesh Ewald (PME) summation method [90]. SHAKE algorithm [91] was applied to constrain all chemical bonds involving hydrogen. The prepared systems were heated up from 10.0 K to 298.0 K for 100 ps. Subsequently, the MD simulations with NPT ensemble were performed at this temperature until reaching 500 ns. From root-mean-square displacement (RMSD) analysis shown in Figure S5, the equilibrated MD trajectories in the last 200-ns simulations of all systems were extracted for further analysis. The cpptraj module was used to compute the structural and dynamics data, including RMSD, ligand-protein H-bond occupation, and protein motion via PCA. Binding free energy of all studied complexes was calculated using the MM/GBSA method [92]. Additionally, the per-residue decomposition free energy based on MM/GBSA approach was evaluated in order to identify crucial amino acids important for ligand recognition.

5. Conclusions

The experimental and theoretical results obtained in this study shed light on the anticancer activity and its underlying mechanisms of butoxy MG against human NSCLC cell lines expressing wild-type EGFR and mutant EGFR, which might be useful to develop this compound as a novel anticancer agent and/or can be used as a theoretical guidance for designing and developing a new compound targeting STAT3 and Akt signaling pathways.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/11/4/437/s1, Figure S1: The percentage of H-bond occupation of the amino acid residues contributing to all ligands during the last 200-ns simulations within (A) SH2 domain of STAT3 and (B) ATP-binding pocket of Akt; Figure S2: The binding orientation of MG3 against Akt signaling protein taken from the last snapshot of 500-ns MD simulation; Figure S3: (A) The PCA result of Akt1 model. (B) the superimposed crystal structures between apo (PDB ID: 1GZN, black) and holo forms (PDB ID: 4GV1, pink) of Akt; Figure S4: The distance between the Cm of MG3 and DNA (d(Cm(MG3)-Cm(DNA))) of three independent simulations (MD1-3); Figure S5: RMSD plots of (A) STAT3 and (B) Akt models; Figure S6: Morphology of PCS201-010 cells after treatment with MG3 and CDDP for 24 h; Table S1: The computational details of all initial structures used for MD simulations.

Author Contributions

T.R. and P.W. conceived and designed the experiments. P.M. conducted theoretical and experimental studies. W.C. purified, synthesized, and identified all compounds using NMR. T.R., P.W., and P.M. analyzed the data. P.M. wrote the original manuscript. All authors reviewed and edited the manuscript.

Funding

This work was supported by the Thailand Research Fund (grant number RSA5980069). P.M. thanks the Science Achievement Scholarship of Thailand for Ph.D. scholarship, the 90th Anniversary of Chulalongkorn University (CU) Fund (Ratchadaphiseksomphot Endowment Fund), and the Overseas Presentations of Graduate Level Academic Thesis from Graduate School of CU.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Ahmedin, J. Cancer statistics, 2018. CA A Cancer J. Clin. 2018, 68, 7–30. [Google Scholar] [CrossRef] [PubMed]
  2. Cheng, T.-Y.D.; Cramb, S.M.; Baade, P.D.; Youlden, D.R.; Nwogu, C.; Reid, M.E. The International Epidemiology of Lung Cancer: Latest Trends, Disparities, and Tumor Characteristics. J. Thorac. Oncol. 2016, 11, 1653–1671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Islami, F.; Torre, L.A.; Jemal, A. Global trends of lung cancer mortality and smoking prevalence. Transl. Lung Cancer Res. 2015, 4, 327–338. [Google Scholar]
  4. Wong, M.C.S.; Lao, X.Q.; Ho, K.-F.; Goggins, W.B.; Tse, S.L.A. Incidence and mortality of lung cancer: Global trends and association with socioeconomic status. Sci. Rep. 2017, 7, 14300. [Google Scholar] [CrossRef]
  5. Zhang, H. Osimertinib making a breakthrough in lung cancer targeted therapy. Onco Targets Ther. 2016, 9, 5489–5493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bethune, G.; Bethune, D.; Ridgway, N.; Xu, Z. Epidermal growth factor receptor (EGFR) in lung cancer: An overview and update. J. Thorac. Dis. 2010, 2, 48–51. [Google Scholar] [PubMed]
  7. Li, X.; Fan, X.X.; Jiang, Z.B.; Loo, W.T.; Yao, X.J.; Leung, E.L.; Chow, L.W.; Liu, L. Shikonin inhibits gefitinib-resistant non-small cell lung cancer by inhibiting TrxR and activating the EGFR proteasomal degradation pathway. Pharmacol. Res. 2017, 115, 45–55. [Google Scholar] [CrossRef]
  8. Seshacharyulu, P.; Ponnusamy, M.P.; Haridas, D.; Jain, M.; Ganti, A.K.; Batra, S.K. Targeting the EGFR signaling pathway in cancer therapy. Expert Opin. Ther. Targets 2012, 16, 15–31. [Google Scholar] [CrossRef] [Green Version]
  9. Yarden, Y.; Sliwkowski, M.X. Untangling the ErbB signalling network. Nat. Rev. Mol. Cell Biol. 2001, 2, 127–137. [Google Scholar] [CrossRef]
  10. Rose, M.C.; Kostyanovskaya, E.; Huang, R.S. Pharmacogenomics of Cisplatin Sensitivity in Non-small Cell Lung Cancer. Genom. Proteom. Bioinform. 2014, 12, 198–209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Su, T.-J.; Ku, W.-H.; Chen, H.-Y.; Hsu, Y.-C.; Hong, Q.-S.; Chang, G.-C.; Yu, S.-L.; Chen, J.J.W. Oncogenic miR-137 contributes to cisplatin resistance via repressing CASP3 in lung adenocarcinoma. Am. J. Cancer Res. 2016, 6, 1317–1330. [Google Scholar] [PubMed]
  12. Fennell, D.A.; Summers, Y.; Cadranel, J.; Benepal, T.; Christoph, D.C.; Lal, R.; Das, M.; Maxwell, F.; Visseren-Grul, C.; Ferry, D. Cisplatin in the modern era: The backbone of first-line chemotherapy for non-small cell lung cancer. Cancer Treat. Rev. 2016, 44, 42–50. [Google Scholar] [CrossRef]
  13. Smit, E.; Moro-Sibilot, D.; Carpeño, J.d.C.; Lesniewski-Kmak, K.; Aerts, J.; Villatoro, R.; Kraaij, K.; Nacerddine, K.; Dyachkova, Y.; Smith, K.T.; et al. Cisplatin and carboplatin-based chemotherapy in the first-line treatment of non-small cell lung cancer: Analysis from the European FRAME study. Lung Cancer (Amst. Neth.) 2016, 92, 35–40. [Google Scholar] [CrossRef]
  14. Ozkaya, S.; Findik, S.; Atici, A.G.; Dirica, A. Cisplatin-based chemotherapy in elderly patients with advanced stage (IIIB and IV) non-small cell lung cancer patients. Neoplasma 2011, 58, 348–351. [Google Scholar] [CrossRef] [Green Version]
  15. Raimbourg, J.; Joalland, M.P.; Cabart, M.; de Plater, L.; Bouquet, F.; Savina, A.; Decaudin, D.; Bennouna, J.; Vallette, F.M.; Lalier, L. Sensitization of EGFR Wild-Type Non-Small Cell Lung Cancer Cells to EGFR-Tyrosine Kinase Inhibitor Erlotinib. Mol. Cancer Ther. 2017, 16, 1634–1644. [Google Scholar] [CrossRef] [PubMed]
  16. Leighl, N.B. Treatment paradigms for patients with metastatic non-small-cell lung cancer: First-, second-, and third-line. Curr. Oncol. (Tor. Ont.) 2012, 19, S52–S58. [Google Scholar] [CrossRef] [PubMed]
  17. Dasari, S.; Tchounwou, P.B. Cisplatin in cancer therapy: Molecular mechanisms of action. Eur. J. Pharmacol. 2014, 740, 364–378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Sun, X.; Zheng, Y. Cisplatin or Carboplatin for Advanced Non–Small-Cell Lung Cancer? J. Thorac. Oncol. 2014, 9, e70. [Google Scholar] [CrossRef] [Green Version]
  19. Sarin, N.; Engel, F.; Kalayda, G.V.; Mannewitz, M.; Cinatl, J., Jr.; Rothweiler, F.; Michaelis, M.; Saafan, H.; Ritter, C.A.; Jaehde, U.; et al. Cisplatin resistance in non-small cell lung cancer cells is associated with an abrogation of cisplatin-induced G2/M cell cycle arrest. PLoS ONE 2017, 12, e0181081. [Google Scholar] [CrossRef]
  20. Rosell, R.; Taron, M.; Barnadas, A.; Scagliotti, G.; Sarries, C.; Roig, B. Nucleotide excision repair pathways involved in Cisplatin resistance in non-small-cell lung cancer. Cancer Control. J. Moffitt Cancer Cent. 2003, 10, 297–305. [Google Scholar] [CrossRef] [PubMed]
  21. Kuo, M.T.; Chen, H.H.; Song, I.S.; Savaraj, N.; Ishikawa, T. The roles of copper transporters in cisplatin resistance. Cancer Metastasis Rev. 2007, 26, 71–83. [Google Scholar] [CrossRef] [PubMed]
  22. Galluzzi, L.; Senovilla, L.; Vitale, I.; Michels, J.; Martins, I.; Kepp, O.; Castedo, M.; Kroemer, G. Molecular mechanisms of cisplatin resistance. Oncogene 2011, 31, 1869. [Google Scholar] [CrossRef]
  23. Weber, B.; Hager, H.; Sorensen, B.S.; McCulloch, T.; Mellemgaard, A.; Khalil, A.A.; Nexo, E.; Meldgaard, P. EGFR mutation frequency and effectiveness of erlotinib: A prospective observational study in Danish patients with non-small cell lung cancer. Lung Cancer (Amst. Neth.) 2014, 83, 224–230. [Google Scholar] [CrossRef] [PubMed]
  24. Matsuo, N.; Azuma, K.; Sakai, K.; Hattori, S.; Kawahara, A.; Ishii, H.; Tokito, T.; Kinoshita, T.; Yamada, K.; Nishio, K.; et al. Association of EGFR Exon 19 Deletion and EGFR-TKI Treatment Duration with Frequency of T790M Mutation in EGFR-Mutant Lung Cancer Patients. Sci. Rep. 2016, 6, 36458. [Google Scholar] [CrossRef] [Green Version]
  25. Wakeling, A.E.; Guy, S.P.; Woodburn, J.R.; Ashton, S.E.; Curry, B.J.; Barker, A.J.; Gibson, K.H. ZD1839 (Iressa). An Orally Active Inhibitor of Epidermal Growth Factor Signaling with Potential for Cancer Therapy. Cancer Res. 2002, 62, 5749–5754. [Google Scholar] [PubMed]
  26. Wang, S.; Cang, S.; Liu, D. Third-generation inhibitors targeting EGFR T790M mutation in advanced non-small cell lung cancer. J. Hematol. Oncol. 2016, 9, 34. [Google Scholar] [CrossRef] [Green Version]
  27. Mayor, S. Osimertinib effective in EGFR T790M-positive lung cancer. Lancet Oncol. 2017, 18, e9. [Google Scholar] [CrossRef]
  28. Ma, C.; Wei, S.; Song, Y. T790M and acquired resistance of EGFR TKI: A literature review of clinical reports. J. Thorac. Dis. 2011, 3, 10–18. [Google Scholar]
  29. Ko, B.; Paucar, D.; Halmos, B. EGFR T790M: Revealing the secrets of a gatekeeper. Lung Cancer Targets Ther. 2017, 8, 147–159. [Google Scholar] [CrossRef]
  30. Hairani, R.; Mongkol, R.; Chavasiri, W. Allyl and prenyl ethers of mansonone G, new potential semisynthetic antibacterial agents. Bioorg. Med. Chem. Lett. 2016, 26, 5300–5303. [Google Scholar] [CrossRef]
  31. Dai, Y.; Harinantenaina, L.; Brodie, P.J.; Callmander, M.W.; Randrianasolo, S.; Rakotobe, E.; Rasamison, V.E.; Kingston, D.G.I. Isolation and Synthesis of Two Antiproliferative Calamenene-type Sesquiterpenoids from Sterculia tavia from the Madagascar Rain Forest. Bioorg. Med. Chem. 2012, 20, 6940–6944. [Google Scholar] [CrossRef]
  32. Mongkol, R.; Chavasiri, W. Antimicrobial, herbicidal and antifeedant activities of mansonone E from the heartwoods of Mansonia gagei Drumm. J. Integr. Agric. 2016, 15, 2795–2802. [Google Scholar] [CrossRef]
  33. El-Halawany, A.M.; Salah El Dine, R.; Hattori, M. Anti-estrogenic activity of mansonone G and mansorin A derivatives. Pharm. Biol. 2013, 51, 948–954. [Google Scholar] [CrossRef] [Green Version]
  34. Changwong, N.; Sabphon, C.; Ingkaninan, K.; Sawasdee, P. Acetyl- and butyryl-cholinesterase inhibitory activities of mansorins and mansonones. Phytother. Res. PTR 2012, 26, 392–396. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, H.K.; Hairani, R.; Jeong, H.; Jeong, M.G.; Chavasiri, W.; Hwang, E.S. CBMG, a novel derivative of mansonone G suppresses adipocyte differentiation via suppression of PPARγ activity. Chem. Biol. Interact. 2017, 273, 160–170. [Google Scholar] [CrossRef]
  36. Cetintas, V.B.; Kucukaslan, A.S.; Kosova, B.; Tetik, A.; Selvi, N.; Cok, G.; Gunduz, C.; Eroglu, Z. Cisplatin resistance induced by decreased apoptotic activity in non-small-cell lung cancer cell lines. Cell Biol. Int. 2012, 36, 261–265. [Google Scholar] [CrossRef] [PubMed]
  37. Bortner, C.D.; Cidlowski, J.A. A necessary role for cell shrinkage in apoptosis. Biochem. Pharmacol. 1998, 56, 1549–1559. [Google Scholar] [CrossRef]
  38. Yen, H.-Y.; Liu, Y.-C.; Chen, N.-Y.; Tsai, C.-F.; Wang, Y.-T.; Chen, Y.-J.; Hsu, T.-L.; Yang, P.-C.; Wong, C.-H. Effect of sialylation on EGFR phosphorylation and resistance to tyrosine kinase inhibition. Proc. Natl. Acad. Sci. USA 2015, 112, 6955–6960. [Google Scholar] [CrossRef] [Green Version]
  39. Lelj-Garolla, B.; Kumano, M.; Beraldi, E.; Nappi, L.; Rocchi, P.; Ionescu, D.N.; Fazli, L.; Zoubeidi, A.; Gleave, M.E. Hsp27 Inhibition with OGX-427 Sensitizes Non–Small Cell Lung Cancer Cells to Erlotinib and Chemotherapy. Mol. Cancer Ther. 2015, 14, 1107–1116. [Google Scholar] [CrossRef]
  40. Shin, D.-S.; Kim, H.-N.; Shin, K.D.; Yoon, Y.J.; Kim, S.-J.; Han, D.C.; Kwon, B.-M. Cryptotanshinone Inhibits Constitutive Signal Transducer and Activator of Transcription 3 Function through Blocking the Dimerization in DU145 Prostate Cancer Cells. Cancer Res. 2009, 69, 193–202. [Google Scholar] [CrossRef] [Green Version]
  41. Siddiquee, K.; Zhang, S.; Guida, W.C.; Blaskovich, M.A.; Greedy, B.; Lawrence, H.R.; Yip, M.L.R.; Jove, R.; McLaughlin, M.M.; Lawrence, N.J.; et al. Selective chemical probe inhibitor of Stat3, identified through structure-based virtual screening, induces antitumor activity. Proc. Natl. Acad. Sci. USA 2007, 104, 7391–7396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Pachl, F.; Plattner, P.; Ruprecht, B.; Medard, G.; Sewald, N.; Kuster, B. Characterization of a chemical affinity probe targeting Akt kinases. J. Proteome Res. 2013, 12, 3792–3800. [Google Scholar] [CrossRef] [PubMed]
  43. Burns, S.; Travers, J.; Collins, I.; Rowlands, M.G.; Newbatt, Y.; Thompson, N.; Garrett, M.D.; Workman, P.; Aherne, W. Identification of small-molecule inhibitors of protein kinase B (PKB/AKT) in an AlphaScreenTM high-throughput screen. J. Biomol. Screen. 2006, 11, 822–827. [Google Scholar] [CrossRef]
  44. Samadi Moghaddam, M.; Heiny, M.; Shastri, V.P. Enhanced cellular uptake of nanoparticles by increasing the hydrophobicity of poly (lactic acid) through copolymerization with cell-membrane-lipid components. Chem. Commun. 2015, 51, 14605–14608. [Google Scholar] [CrossRef] [PubMed]
  45. Franks, N.P.; Lieb, W.R. Partitioning of long-chain alcohols into lipid bilayers: Implications for mechanisms of general anesthesia. Proc. Natl. Acad. Sci. USA 1986, 83, 5116–5120. [Google Scholar] [CrossRef]
  46. Li, Y.; Song, Z.; Jin, Y.; Tang, Z.; Kang, J.; Ma, X. Novel Selective and Potent EGFR Inhibitor that Overcomes T790M-Mediated Resistance in Non-Small Cell Lung Cancer. Molecules (Basel Switz.) 2016, 21, 1462. [Google Scholar] [CrossRef] [PubMed]
  47. Vazirian, M.; Naser Ostad, S.; Manayi, A.; Hadjiakhoondi, A.; Khanavi, M. Comparison of cytotoxic activity of some Iranian Stachys spp. extracts on different cancer cell lines. Res. J. Pharmacogn. 2014, 1, 23–28. [Google Scholar]
  48. Patel, R.M.; Patel, S.K. Cytotoxic activity of methanolic extract of artocarpus heterophyllus against a549, hela and mcf-7 cell lines. J. Appl. Pharm. 2011, 1, 167–171. [Google Scholar]
  49. McIlwain, D.R.; Berger, T.; Mak, T.W. Caspase functions in cell death and disease. Cold Spring Harb. Perspect. Biol. 2013, 5, a008656. [Google Scholar] [CrossRef]
  50. Suda, K.; Onozato, R.; Yatabe, Y.; Mitsudomi, T. EGFR T790M Mutation: A Double Role in Lung Cancer Cell Survival? J. Thorac. Oncol. 2009, 4, 1–4. [Google Scholar] [CrossRef] [Green Version]
  51. Henkels, K.M.; Turchi, J.J. Cisplatin-induced apoptosis proceeds by caspase-3-dependent and -independent pathways in cisplatin-resistant and -sensitive human ovarian cancer cell lines. Cancer Res. 1999, 59, 3077–3083. [Google Scholar] [PubMed]
  52. Cummings, B.S.; Schnellmann, R.G. Cisplatin-induced renal cell apoptosis: Caspase 3-dependent and -independent pathways. J. Pharmacol. Exp. Ther. 2002, 302, 8–17. [Google Scholar] [CrossRef]
  53. Kaushal, G.P.; Kaushal, V.; Hong, X.; Shah, S.V. Role and regulation of activation of caspases in cisplatin-induced injury to renal tubular epithelial cells. Kidney Int. 2001, 60, 1726–1736. [Google Scholar] [CrossRef] [Green Version]
  54. Dai, C.H.; Li, J.; Chen, P.; Jiang, H.G.; Wu, M.; Chen, Y.C. RNA interferences targeting the Fanconi anemia/BRCA pathway upstream genes reverse cisplatin resistance in drug-resistant lung cancer cells. J. Biomed. Sci. 2015, 22, 77. [Google Scholar] [CrossRef] [Green Version]
  55. Aziz, M.H.; Dreckschmidt, N.E.; Verma, A.K. Plumbagin, a Medicinal Plant–Derived Naphthoquinone, Is a Novel Inhibitor of the Growth and Invasion of Hormone-Refractory Prostate Cancer. Cancer Res. 2008, 68, 9024–9032. [Google Scholar] [CrossRef]
  56. Sandur, S.K.; Pandey, M.K.; Sung, B.; Aggarwal, B.B. 5-hydroxy-2-methyl-1, 4-naphthoquinone, a vitamin K3 analogue, suppresses STAT3 activation pathway through induction of protein tyrosine phosphatase, SHP-1: Potential role in chemosensitization. Mol. Cancer Res. MCR 2010, 8, 107–118. [Google Scholar] [CrossRef] [PubMed]
  57. Sand, J.M.; Bin Hafeez, B.; Jamal, M.S.; Witkowsky, O.; Siebers, E.M.; Fischer, J.; Verma, A.K. Plumbagin (5-hydroxy-2-methyl-1,4-naphthoquinone), isolated from Plumbago zeylanica, inhibits ultraviolet radiation-induced development of squamous cell carcinomas. Carcinogenesis 2012, 33, 184–190. [Google Scholar] [CrossRef] [PubMed]
  58. Qiu, H.Y.; Fu, J.Y.; Yang, M.K.; Han, H.W.; Wang, P.F.; Zhang, Y.H.; Lin, H.Y.; Tang, C.Y.; Qi, J.L.; Yang, R.W.; et al. Identification of new shikonin derivatives as STAT3 inhibitors. Biochem. Pharmacol. 2017, 146, 74–86. [Google Scholar] [CrossRef] [PubMed]
  59. Qiu, H.-Y.; Zhu, X.; Luo, Y.-L.; Lin, H.-Y.; Tang, C.-Y.; Qi, J.-L.; Pang, Y.-J.; Yang, R.-W.; Lu, G.-H.; Wang, X.-M.; et al. Identification of New Shikonin Derivatives as Antitumor Agents Targeting STAT3 SH2 Domain. Sci. Rep. 2017, 7, 2863. [Google Scholar] [CrossRef]
  60. Kawiak, A.; Lojkowska, E. Ramentaceone, a Naphthoquinone Derived from Drosera sp., Induces Apoptosis by Suppressing PI3K/Akt Signaling in Breast Cancer Cells. PLoS ONE 2016, 11, e0147718. [Google Scholar] [CrossRef]
  61. Nishina, A.; Miura, A.; Goto, M.; Terakado, K.; Sato, D.; Kimura, H.; Hirai, Y.; Sato, H.; Phay, N. Mansonone E from Mansonia gagei Inhibited alpha-MSH-Induced Melanogenesis in B16 Cells by Inhibiting CREB Expression and Phosphorylation in the PI3K/Akt Pathway. Biol. Pharm. Bull. 2018, 41, 770–776. [Google Scholar] [CrossRef] [PubMed]
  62. Fletcher, S.; Singh, J.; Zhang, X.; Yue, P.; Page, B.D.G.; Sharmeen, S.; Shahani, V.M.; Zhao, W.; Schimmer, A.D.; Turkson, J.; et al. Disruption of Transcriptionally Active Stat3 Dimers with Non-phosphorylated, Salicylic Acid-Based Small Molecules: Potent in vitro and Tumor Cell Activities. Chembiochem A Eur. J. Chem. Biol. 2009, 10, 1959–1964. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Liu, L.J.; Leung, K.H.; Chan, D.S.; Wang, Y.T.; Ma, D.L.; Leung, C.H. Identification of a natural product-like STAT3 dimerization inhibitor by structure-based virtual screening. Cell Death Dis. 2014, 5, e1293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Shao, S.; Yu, R.; Yu, Y.; Li, Y. Dual-inhibitors of STAT5 and STAT3: Studies from molecular docking and molecular dynamics simulations. J. Mol. Model. 2014, 20, 2399. [Google Scholar] [CrossRef] [PubMed]
  65. Sgrignani, J.; Garofalo, M.; Matkovic, M.; Merulla, J.; Catapano, C.V.; Cavalli, A. Structural Biology of STAT3 and Its Implications for Anticancer Therapies Development. Int. J. Mol. Sci. 2018, 19, 1591. [Google Scholar] [CrossRef] [PubMed]
  66. Davies, T.G.; Verdonk, M.L.; Graham, B.; Saalau-Bethell, S.; Hamlett, C.C.; McHardy, T.; Collins, I.; Garrett, M.D.; Workman, P.; Woodhead, S.J.; et al. A structural comparison of inhibitor binding to PKB, PKA and PKA-PKB chimera. J. Mol. Biol. 2007, 367, 882–894. [Google Scholar] [CrossRef]
  67. Allen, F.H. The Cambridge Structural Database: A quarter of a million crystal structures and rising. Acta Crystallogr. Sect. B Struct. Sci. 2002, 58, 380–388. [Google Scholar] [CrossRef]
  68. Jackson, P.A.; Widen, J.C.; Harki, D.A.; Brummond, K.M. Covalent Modifiers: A Chemical Perspective on the Reactivity of α,β-Unsaturated Carbonyls with Thiols via Hetero-Michael Addition Reactions. J. Med. Chem. 2017, 60, 839–885. [Google Scholar] [CrossRef]
  69. Yang, J.; Cron, P.; Thompson, V.; Good, V.M.; Hess, D.; Hemmings, B.A.; Barford, D. Molecular mechanism for the regulation of protein kinase B/Akt by hydrophobic motif phosphorylation. Mol. Cell 2002, 9, 1227–1240. [Google Scholar] [CrossRef]
  70. Addie, M.; Ballard, P.; Buttar, D.; Crafter, C.; Currie, G.; Davies, B.R.; Debreczeni, J.; Dry, H.; Dudley, P.; Greenwood, R.; et al. Discovery of 4-amino-N-[(1S)-1-(4-chlorophenyl)-3-hydroxypropyl]-1-(7H-pyrrolo[2,3-d]pyrimidin -4-yl)piperidine-4-carboxamide (AZD5363), an orally bioavailable, potent inhibitor of Akt kinases. J. Med. Chem. 2013, 56, 2059–2073. [Google Scholar] [CrossRef]
  71. Siddik, Z.H. Cisplatin: Mode of cytotoxic action and molecular basis of resistance. Oncogene 2002, 22, 7265. [Google Scholar] [CrossRef] [PubMed]
  72. Tan, B.J.; Chiu, G.N. Role of oxidative stress, endoplasmic reticulum stress and ERK activation in triptolide-induced apoptosis. Int. J. Oncol. 2013, 42, 1605–1612. [Google Scholar] [CrossRef] [PubMed]
  73. Zhuang, S.; Schnellmann, R.G. A death-promoting role for extracellular signal-regulated kinase. J. Pharmacol. Exp. Ther. 2006, 319, 991–997. [Google Scholar] [CrossRef] [PubMed]
  74. Persons, D.L.; Yazlovitskaya, E.M.; Pelling, J.C. Effect of extracellular signal-regulated kinase on p53 accumulation in response to cisplatin. J. Biol. Chem. 2000, 275, 35778–35785. [Google Scholar] [CrossRef]
  75. Greig, M.J.; Niessen, S.; Weinrich, S.L.; Feng, J.L.; Shi, M.; Johnson, T.O. Effects of Activating Mutations on EGFR Cellular Protein Turnover and Amino Acid Recycling Determined Using SILAC Mass Spectrometry. Int. J. Cell Biol. 2015, 2015, 798936. [Google Scholar] [CrossRef]
  76. Becker, S.; Groner, B.; Muller, C.W. Three-dimensional structure of the Stat3beta homodimer bound to DNA. Nature 1998, 394, 145–151. [Google Scholar] [CrossRef]
  77. Schwede, T.; Kopp, J.; Guex, N.; Peitsch, M.C. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res. 2003, 31, 3381–3385. [Google Scholar] [CrossRef]
  78. Mahalapbutr, P.; Chusuth, P.; Kungwan, N.; Chavasiri, W.; Wolschann, P.; Rungrotmongkol, T. Molecular recognition of naphthoquinone-containing compounds against human DNA topoisomerase IIα ATPase domain: A molecular modeling study. J. Mol. Liq. 2017, 247, 374–385. [Google Scholar] [CrossRef]
  79. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  80. Wu, G.; Robertson, D.H.; Brooks, C.L., 3rd; Vieth, M. Detailed analysis of grid-based molecular docking: A case study of CDOCKER-A CHARMm-based MD docking algorithm. J. Comput. Chem. 2003, 24, 1549–1562. [Google Scholar] [CrossRef]
  81. Olsson, M.H.; Sondergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J. Chem. Theory Comput. 2011, 7, 525–537. [Google Scholar] [CrossRef] [PubMed]
  82. Phanich, J.; Rungrotmongkol, T.; Kungwan, N.; Hannongbua, S. Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB (GB)SA study. J. Comput.-Aided Mol. Des. 2016, 30, 917–926. [Google Scholar] [CrossRef] [PubMed]
  83. Mahalapbutr, P.; Nutho, B.; Wolschann, P.; Chavasiri, W.; Kungwan, N.; Rungrotmongkol, T. Molecular insights into inclusion complexes of mansonone E and H enantiomers with various beta-cyclodextrins. J. Mol. Graph. Model. 2018, 79, 72–80. [Google Scholar] [CrossRef] [PubMed]
  84. Meeprasert, A.; Hannongbua, S.; Rungrotmongkol, T. Key Binding and Susceptibility of NS3/4A Serine Protease Inhibitors against Hepatitis C Virus. J. Chem. Inf. Model. 2014, 54, 1208–1217. [Google Scholar] [CrossRef]
  85. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [PubMed]
  86. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Sangpheak, W.; Khuntawee, W.; Wolschann, P.; Pongsawasdi, P.; Rungrotmongkol, T. Enhanced stability of a naringenin/2, 6-dimethyl β-cyclodextrin inclusion complex: Molecular dynamics and free energy calculations based on MM- and QM-PBSA/GBSA. J. Mol. Graph. Model. 2014, 50, 10–15. [Google Scholar] [CrossRef]
  88. Mahalapbutr, P.; Thitinanthavet, K.; Kedkham, T.; Nguyen, H.; Theu, L.T.H.; Dokmaisrijan, S.; Huynh, L.; Kungwan, N.; Rungrotmongkol, T. A theoretical study on the molecular encapsulation of luteolin and pinocembrin with various derivatized beta-cyclodextrins. J. Mol. Struct. 2019, 1180, 480–490. [Google Scholar] [CrossRef]
  89. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  90. York, D.M.; Darden, T.A.; Pedersen, L.G. The effect of long-range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods. J. Chem. Phys. 1993, 99, 8345–8348. [Google Scholar] [CrossRef]
  91. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H.J.C. Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327–341. [Google Scholar] [CrossRef]
  92. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 11 00437 g001 550
Figure 1. Two-dimensional (2D) chemical structures of (A) MG and its semi-synthetic ether derivatives MG1-MG10 [30] and (B) the known STAT3 (cryptotanshinone (CST) and S3I201) and Akt (uprosertib and H8) inhibitors. Three-dimensional (3D) structures of (C) STAT3 and (D) Akt1 signaling proteins. The SH2 domain of STAT3 and the ATP-binding pocket of Akt are shown by blue surface and black circle, respectively.
Figure 1. Two-dimensional (2D) chemical structures of (A) MG and its semi-synthetic ether derivatives MG1-MG10 [30] and (B) the known STAT3 (cryptotanshinone (CST) and S3I201) and Akt (uprosertib and H8) inhibitors. Three-dimensional (3D) structures of (C) STAT3 and (D) Akt1 signaling proteins. The SH2 domain of STAT3 and the ATP-binding pocket of Akt are shown by blue surface and black circle, respectively.
Cancers 11 00437 g001
Cancers 11 00437 g002 550
Figure 2. Cell viability of NSCLC (A–C) and PCS201-010 (D–F) cell lines after treatment with MG3, MG4, and CDDP for 48 h. The IC50 (μM) and SI of three focused compounds against all studied cell lines are shown in (G) and (H), respectively. (I) The morphological changes of two NSCLC cell lines treated with MG3 at various concentrations for 48 h.
Figure 2. Cell viability of NSCLC (A–C) and PCS201-010 (D–F) cell lines after treatment with MG3, MG4, and CDDP for 48 h. The IC50 (μM) and SI of three focused compounds against all studied cell lines are shown in (G) and (H), respectively. (I) The morphological changes of two NSCLC cell lines treated with MG3 at various concentrations for 48 h.
Cancers 11 00437 g002
Cancers 11 00437 g003 550
Figure 3. Flow cytometric analysis of Annexin V/PI stained cells after MG3 and CDDP treatments for 24 h on (A) A549 and (B) H1975 cells. Western blot analysis of apoptotic markers, caspase-3 and PARP, for (C) A549 and (D) H1975 cells. Inhibition of MG3-induced apoptosis by the pan caspase inhibitor Z-VAD(OMe)-FMK in (E) A549 and (F) H1975 cells. Data are expressed as mean ± SEM (n = 3). * p ≤ 0.05, ** p ≤ 0.01, and *** p ≤ 0.001 vs. control. @ p ≤ 0.05, @@ p ≤ 0.01, and @@@ p ≤ 0.001 vs. MG3. # p ≤ 0.05, ## p ≤ 0.01, and ### p ≤ 0.001 vs. CDDP.
Figure 3. Flow cytometric analysis of Annexin V/PI stained cells after MG3 and CDDP treatments for 24 h on (A) A549 and (B) H1975 cells. Western blot analysis of apoptotic markers, caspase-3 and PARP, for (C) A549 and (D) H1975 cells. Inhibition of MG3-induced apoptosis by the pan caspase inhibitor Z-VAD(OMe)-FMK in (E) A549 and (F) H1975 cells. Data are expressed as mean ± SEM (n = 3). * p ≤ 0.05, ** p ≤ 0.01, and *** p ≤ 0.001 vs. control. @ p ≤ 0.05, @@ p ≤ 0.01, and @@@ p ≤ 0.001 vs. MG3. # p ≤ 0.05, ## p ≤ 0.01, and ### p ≤ 0.001 vs. CDDP.
Cancers 11 00437 g003
Cancers 11 00437 g004 550
Figure 4. MG3 dose-dependently inhibits the phosphorylation of STAT3 (Tyr705) and Akt (Ser473) in (A) A549 and (B) H1975 cells at 24 h. The expression of p-Erk (Thr202/Tyr204) is found to be increased upon MG3 and CDDP treatments. The phosphorylation of EGFR (Y1068) for both (C) A549 and (D) H1975 cells is not significantly affected by MG3 and CDDP treatments. (E) Proposed mechanisms of MG3 against two studied NSCLC cell lines, in which MG3 promotes cell apoptosis through the inhibition of p-Akt and p-STAT3 as well as through the activation of MAPK signaling pathway. Data are expressed as mean ± SEM of three independent experiments. * p ≤ 0.05, ** p ≤ 0.01, and *** p ≤ 0.001 vs. control.
Figure 4. MG3 dose-dependently inhibits the phosphorylation of STAT3 (Tyr705) and Akt (Ser473) in (A) A549 and (B) H1975 cells at 24 h. The expression of p-Erk (Thr202/Tyr204) is found to be increased upon MG3 and CDDP treatments. The phosphorylation of EGFR (Y1068) for both (C) A549 and (D) H1975 cells is not significantly affected by MG3 and CDDP treatments. (E) Proposed mechanisms of MG3 against two studied NSCLC cell lines, in which MG3 promotes cell apoptosis through the inhibition of p-Akt and p-STAT3 as well as through the activation of MAPK signaling pathway. Data are expressed as mean ± SEM of three independent experiments. * p ≤ 0.05, ** p ≤ 0.01, and *** p ≤ 0.001 vs. control.
Cancers 11 00437 g004
Cancers 11 00437 g005 550
Figure 5. Δ G bind residue (kcal/mol) of (A) STAT3 and (B) Akt. The amino acids involved in ligand binding are shaded according to their Δ G bind residue , in which the highest and lowest energies are ranged from red to magenta, respectively.
Figure 5. Δ G bind residue (kcal/mol) of (A) STAT3 and (B) Akt. The amino acids involved in ligand binding are shaded according to their Δ G bind residue , in which the highest and lowest energies are ranged from red to magenta, respectively.
Cancers 11 00437 g005
Cancers 11 00437 g006 550
Figure 6. Porcupine plot of apo and holo forms of (A) STAT3 and (B) Akt across the PC1. The residues within 5 Å of ligand are presented in circle where the MG3 molecule is shown in vdW model. The PCA scree plot of quantitative characters for each protein is given below, in which the columns in black and grey represent the data for apo and holo forms, respectively.
Figure 6. Porcupine plot of apo and holo forms of (A) STAT3 and (B) Akt across the PC1. The residues within 5 Å of ligand are presented in circle where the MG3 molecule is shown in vdW model. The PCA scree plot of quantitative characters for each protein is given below, in which the columns in black and grey represent the data for apo and holo forms, respectively.
Cancers 11 00437 g006
Table
Table 1. In vitro cytotoxicity screening of all MG analogs against A549 and H1975 NSCLC cell lines. Cells were treated with indicated compounds at 10 μM and 100 μM for 48 h, and cell viability was determined using MTT assay. Data are expressed as mean ± SEM of two independent experiments. MG derivatives exhibiting a percentage of cell viability at 10 μM (%CV10 µM) <50 were defined as potent compounds.
Table 1. In vitro cytotoxicity screening of all MG analogs against A549 and H1975 NSCLC cell lines. Cells were treated with indicated compounds at 10 μM and 100 μM for 48 h, and cell viability was determined using MTT assay. Data are expressed as mean ± SEM of two independent experiments. MG derivatives exhibiting a percentage of cell viability at 10 μM (%CV10 µM) <50 were defined as potent compounds.
MG AnalogsA549H1975
%CV10 µM%CV100 µM%CV10 µM%CV100 µM
MG95.31 ± 2.3510.71 ± 0.9577.28 ± 5.808.26 ± 0.79
MG170.55 ± 0.279.59 ± 1.4840.46 ± 2.837.04 ± 0.22
MG263.57 ± 0.918.51 ± 0.6628.33 ± 0.827.29 ± 0.57
MG346.90 ± 1.218.94 ± 0.6717.04 ± 1.426.50 ± 0.28
MG47.93 ± 0.437.10 ± 0.22 a23.52 ± 1.447.42 ± 0.55 a
MG5121.54 ± 0.3612.12 ± 1.7144.35 ± 2.987.79 ± 0.12
MG679.56 ± 4.5312.54 ± 1.3320.64 ± 1.176.73 ± 0.15
MG771.84 ± 1.157.45 ± 0.0427.53 ± 1.546.57 ± 0.25
MG866.59 ± 5.987.93 ± 0.0620.24 ± 0.928.25 ± 0.32
MG9110.92 ± 1.7272.86 ± 3.7589.84 ± 2.3929.50 ± 2.80
MG1051.00 ± 0.998.74 ± 0.30 a27.85 ± 2.3523.19 ± 1.38 a
a The %CV was determined at 50 µM due to the low solubility of compound.
Table
Table 2. The MM/GBSA ΔGbind and its energy components (kcal/mol). The ΔGbind, exp was calculated using the equation of ΔGbind, exp = RTlnIC50, where R is the gas constant (1.985 × 10−3 kcal/mol/K), T is the experimental temperature (K), and IC50 is the half maximal inhibitory concentration (μM).
Table 2. The MM/GBSA ΔGbind and its energy components (kcal/mol). The ΔGbind, exp was calculated using the equation of ΔGbind, exp = RTlnIC50, where R is the gas constant (1.985 × 10−3 kcal/mol/K), T is the experimental temperature (K), and IC50 is the half maximal inhibitory concentration (μM).
STAT3Akt
CTSS3I201MG3UprosertibH8MG3
ΔEele−6.80 ± 0.44−110.36 ± 3.35−1.56 ± 0.30−146.76 ± 1.43−168.13 ± 1.16−9.43 ± 0.17
ΔEvdW−35.61 ± 0.21−37.47 ± 0.23−35.77 ±0.22−44.62 ± 0.20−33.92 ± 0.18−39.14 ± 0.17
ΔEMM−42.41 ± 0.48−147.83 ± 3.36−37.33 ±0.37−191.39 ± 1.48−202.05 ± 1.18−48.57 ± 0.24
ΔGsolv, non-polar−4.46 ± 0.02−5.71 ± 0.02−4.89 ± 0.02−6.04 ± 0.02−4.95 ± 0.01−4.85 ± 0.01
ΔGsolv, polar23.99 ± 0.39124.43 ± 3.1515.67 ± 0.27162.16 ± 1.22176.02 ± 1.1224.55 ± 0.15
ΔGsolv19.53 ± 0.39118.72 ± 3.1510.78 ± 0.27156.11 ± 1.22171.07 ± 1.1219.70 ± 0.15
ΔEele+ΔGsolv, polar17.19 ± 0.5814.07 ± 4.5914.11 ± 0.4015.40 ± 1.877.89 ± 1.6115.12 ± 0.22
ΔEvdW+ΔGsolv, non-polar−40.07 ± 0.21−43.18 ± 0.23−40.66 ±0.22−50.66 ± 0.20−38.87 ± 0.18−43.99 ± 0.17
−TΔS17.78 ± 1.6925.38 ± 2.1618.02 ± 1.9324.82 ± 0.9021.28 ± 0.6619.67 ± 1.54
ΔGbind−5.09 ± 0.42−3.73 ± 0.58−8.54 ± 0.48−10.45 ± 0.38−9.68 ± 0.29−9.19 ± 0.39
ΔGbind, exp−7.26−5.54n/a−9.19−6.46n/a
IC50 (μM)4.6 [40]86 [41]n/a0.18 [42]18 [43]n/a

Share and Cite

MDPI and ACS Style

Mahalapbutr, P.; Wonganan, P.; Chavasiri, W.; Rungrotmongkol, T. Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in Non-Small Cell Lung Cancers: Combined Experimental and Theoretical Investigations. Cancers 2019, 11, 437. https://doi.org/10.3390/cancers11040437

AMA Style

Mahalapbutr P, Wonganan P, Chavasiri W, Rungrotmongkol T. Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in Non-Small Cell Lung Cancers: Combined Experimental and Theoretical Investigations. Cancers. 2019; 11(4):437. https://doi.org/10.3390/cancers11040437

Chicago/Turabian Style

Mahalapbutr, Panupong, Piyanuch Wonganan, Warinthorn Chavasiri, and Thanyada Rungrotmongkol. 2019. "Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in Non-Small Cell Lung Cancers: Combined Experimental and Theoretical Investigations" Cancers 11, no. 4: 437. https://doi.org/10.3390/cancers11040437

APA Style

Mahalapbutr, P., Wonganan, P., Chavasiri, W., & Rungrotmongkol, T. (2019). Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in Non-Small Cell Lung Cancers: Combined Experimental and Theoretical Investigations. Cancers, 11(4), 437. https://doi.org/10.3390/cancers11040437

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop