Abstract
Annual mortality of 8.2 million could be attributable to cancer globally, posing a serious health issue; particularly, the high number of nonsmall cell lung cancer (NSCLC) diagnosed cases in recent years highlight the need for development in anticancer agents. In NSCLC, a number of specific inhibitors of phosphatidylinositol-3-kinase (PI3K), Protein kinase B (AKT), and mammalian target of rapamycin are currently under development; however, the early evidence has yielded disappointing results. Ent-kaurane diterpenoid compounds from Cronton tonkinensis have been investigated for several bioactivities such as antibacterial, cytotoxic activity, and so on;; however, lung cancer is not yet studied. In this study, we conducted a molecular docking study of 7 ent-kaurane diterpenoids from C tonkinensis against PI3K targeted anticancer therapies; furthermore, their cytotoxicity effects against A549 lung cancer cells were also evaluated. Obtained results indicated that compounds
Introduction
Approximately 1.79 million deaths could be attributed to lung cancer annually, making it one of the most prevalent cancers worldwide, particularly in Eastern Asia.1,2 Lung cancer often manifests as heterogeneous tumors, in which nonsmall cell lung cancer (NSCLC) occurs in ∼70% to 85% of patient cases. 3 Despite many efforts that have been made to improve the systematic therapies of NSCLC patients, the rate of survival for the patients is only around 15%.4,5 Currently, combinational treatment of NSCLC involving surgery, chemotherapy, and radiotherapy is suggested to be inadequate to improve the prognosis for the patient. 4 Alternatively, patients seek possible treatment in alternative medicines, most notably traditional medicine. In Vietnam, traditional medicine involving the use of herbals and therapies has been widely adopted since a thousand years; it is well established and grounded with various scientific evidences regarding the treatment of various diseases including cancer.6,7
Croton tonkinensis Gagnep. (Euphorbiaceae), or commonly referred to as “Kho sam cho la” in Vietnamese, is a native plant in the north of Vietnam. 8 It is reported that ent-kaurane diterpenoids from C tonkinensis possess a myriad of biological effects including inhibition of silent information regulator 2 ortholog 1, stimulation of osteoblast differentiation, and anti-inflammatory, antibacterial, or antitumor activity.9-12 Of the various ent-kaurane diterpenoids from C tonkinensis described so far in the literature, only a few compounds have been studied against the NSCLC model; therefore, the mechanism by which they affect this cell line is still scarcely explored.10,13-15
In NSCLC, the phosphatidylinositol-3-kinase (PI3K) pathway has been heavily implicated in both tumorigenesis and the progression of disease. 16 There are various genomic episodes that could distort the normal regulation of the pathway such as mutation, amplification, and displacement. In addition, the pathway could be affected by growth factors and regulators to give enhanced metastatic suitability and improved resistance to therapies. 17 The PI3K target proteins are divided into classes and subclasses based on their functions and structures in which class 1A is highly accountable for human cancer development. It has been reported that the activation of this class causes a signaling protein phosphoinositide-dependent kinase-1 (PDK1) to promoted the protein Protein Kinase B (AKT);18-20 thus, mouse double minute 2 homolog (MDM2) movement to the nucleus could also be facilitated ultimately restricting MDM2 to the cytoplasm by undermining p52 functioning (Supplemental Figure S1).21-23 This suggests that the lack of feedback inhibition causes activation of upstream proteins, including AKT, solely through inhibition of the mammalian target of rapamycin (mTOR). These notions raise the need for multiple site targeting when it comes to intervening in the feedback loops, minimizing cancer cell survival and proliferation in the PI3K pathway.24-28 However, the PI3K pathway could also be altered by the inflammatory pathway, which involves arachidonic acid and its metabolism, operating through transformations from cyclooxygenases (COXs) to prostaglandins. 29 This process improves AKT and the aforementioned downstream signals; in addition, previous studies have shown that the tumor multiplicity is inversely correlated with inhibition of 5-lipoxygenase (5-LOX) and that greater effectiveness could be achieved with simultaneous inhibition of both 5-LOX and COX-2.30,31 Currently, there have been continuous efforts conducted by scientists in developing novel agents that specifically inhibit PI3K pathway targets such as gefitinib, perifosine, sirolimus, and so on. Several of them have processed to the clinical trials stage; however, the early evidences showed disappointing results.32,33
In this study, 7 ent-kaurane diterpenoid compounds originated from C tonkinensis were selected for molecular docking simulation to investigate potential inhibition activities with different targets of the PI3K pathway; in addition, their in vitro cytotoxicity effects against human lung cancer cell line A549 were then assessed. The simulated and experimental results are expected to contribute to enhance cancer therapy and guide further in vitro and in vivo studies.
Results and Discussion
Molecular Docking Study
Molegro Virtual Docker (MVD) 34 and Ligplot+ v2.1 35 were utilized to give insight into the orientation of the molecule at the active site of a targeted protein and identify residues participating with a ligand. In the docking results, detailed analysis including binding affinity calculation and hydrogen bond interaction were made for the top-ranked dock pose of each compound.
Docking score analysis revealed that among the studied molecules, compound

Hydrogen bonding interaction: (A) Compound
Docking Scores of Ent-Kaurane Diterpenoids and Oxaliplatin Against AKT, mToR, COX-2, MDM2, and PDK1 Protein Targets.
Abbreviations: AKT, Protein Kinase B; mToR, mammalian target of rapamycin; COX-2, cyclooxygenase-2; MDM2, mouse double minute 2 homolog; PDK1, phosphoinositide dependent kinase-1.
Comparing results with the reference compound, we found that for the AKT model, 4 c/mpounds recorded o have higher docking scores than oxaliplatin including compounds
Hydrophobic interaction analysis showed that all the compounds formed interactions within the active site residues. The residues involved for hydrophobic bonding along with the contributing ligands are described in Supplemental Table S1 and Supplemental Figures S6 to S10. As compound

A similar hydrophobic binding pattern of compound
In general, the obtained results suggested that compound
Cell Viability Study
The A549 cell line, human alveolar basal epithelium cells from an explanted culture of lung carcinomatous tissue, was selected as the target for evaluation of the cytotoxicity of the compounds using dimethylthiazol-diphenyltetrazolium bromide (MTT) assay. The control cells showed high proliferation, which could be considered as 100%. Cells were treated with all 7 compounds with various concentrations for 48 h.
From Table 2, the IC50 of compounds
Cytotoxic Activity (IC50, µM) of Ent-Kaurane Diterpenoids (Compounds
Materials and Methods
Protein and Ligand Preparation
Three-dimensional (3D) coordinates for PI3K proteins AKT (PDB id: 1UNQ), mToR (PDB id: 4JSV), COX-2 (PDB id: 3NT1), MDM2 (PDB id: 4JRG), and PDK1 (PDB id: 5LOV) were collected from PDB. The MVD program (6.0, Molegro-a CLC Bio Company) was selected to prepare proteins by assignment of bonds, bond order, charges, and hybridization state. This process removed water molecules and co-factor metal ions followed by adding hydrogen bonds and Tripos atom types in the next step. Piecewise-linear potential (PLP) was adopted for side-chain minimization for hydrogen and steric interactions. Coulomb potential was used for electrostatic forces with constant backbone and varying torsional angles. Cavities were predicted with a probe size of 1.2 Å and a grid size of 0.8 Å.
The 3D structures of small molecules were prepared using MarvinSketch version 19.27.0 and PyMOL version 2.2.2. The energy minimization was carried out using Gabedit 2.5.0.
Docking of Ligands at Active Sites of the Targets and Binding Affinity Calculations
To score each compound, a grid-based scoring system namely MolDock score (GRID) of MVD was adopted. MolDock SE (simplex evolution) was used as search algorithms. The score includes more interaction types like buried, repulsive, nonpolar, hydrogen bonding, and metal, which are scored as follows:
MolDock score:
The re-rank score for each compound was then used to sort compound poses; the highest scored was considered for evaluating its binding affinity. This process was continued for all the ligands and for all the targets. Binding affinity for the best-ranked pose of each ligand was evaluated in Molrego Data Modeler 3.0 (MDM). Oxaliplatin was considered as a reference for comparing the results.
Hydrophobic Interaction Calculation
Residues involved in hydrophobic interactions were calculated using Ligplot+ v2.1 for all the compounds inside the binding pocket along with the reference compound.
Plant Material Preparation
The ent-kaurane diterpenoids were previously isolated from the leaves of the plant C tonkinensis.10-15 The purification and structure identification of the compounds were performed at the Institute of Natural Products Chemistry-Vietnam Academy of Science and Technology. The purity of the compounds was determined as ≥97% by high-performance liquid chromatography. The ent-kaurane diterpenoids and oxaliplatin structures are described in Figure 3.

(A) The common chemical structure of ent-kaurane-type diterpenoid compounds
Cell Line and Cell Culture
The human lung cancer cell line, A549, was obtained from the Center for Nutraceutical and Pharmaceutical Materials, Myongji University. The culture of cells was performed in Dulbecco’s Modified Eagle Medium. The supplementation for the medium consisted of L-alanyl-L-glutamine (2 mM), penicillin (100 g/ml), streptomycin (100 U/ml), and 10% heat-inactivated fetal bovine serum. The atmosphere in which cells were cultured was humid and had a temperature of 37°C and 5% CO2.
Cell Viability Assay
Upon reaching ∼80% confluency, the detachment of cells was performed using 0.05% trypsin/ethylenediamine tetraacetic acid. Cells were then counted, re-suspended, and added onto a 96-well plate, followed by 24 h of attachment. Every tumor cell line was subjected to 48 h of exposure to 7 diterpenoid compounds at 4 concentrations, respectively (100, 10, 5, and 1 μg/ml). The control was selected as oxaliplatin. An MTT assay was adopted to assess cell viability.
Statistical Analysis
At least 3 experimental attempts were performed for all experiments. Results were expressed as mean ± standard deviation. Statistically significant differences were realized at P < .05 via the Student's t-test.
Conclusions
A total of 7 ent-kaurane diterpenoids from C tonkinensis were screened in silico for their cancer inhibitory properties through the PI3K signaling pathway targeting AKT, mToR, COX-2, MDM2, and PDK1. Furthermore, all these compounds were tested for their cytotoxic effects against the NSCLC cell line (A549). Compound
Supplemental Material
sj-docx-1-npx-10.1177_1934578X211033211 - Supplemental material for Identification of Ent-Kaurane Diterpenoid Compounds as Potential Inhibitors of the PI3K Pathway in Nonsmall Cell Lung Cancer Through Molecular Docking Simulations
Supplemental material, sj-docx-1-npx-10.1177_1934578X211033211 for Identification of Ent-Kaurane Diterpenoid Compounds as Potential Inhibitors of the PI3K Pathway in Nonsmall Cell Lung Cancer Through Molecular Docking Simulations by Pham T. Hong Minh, Tran T. Hoai Van, Tran Q. Toan, Le M. Bui, Nguyen H. Thuan Anh and Pham M. Quan in Natural Product Communications
Footnotes
Acknowledgments
The authors thank Dr Surovi Saikia and Dr Manobjyoti Bordoloi of the Natural Products Chemistry Division, CSIR North East Institute of Science & Technology, India, for their scientific assistance.
Ethical Approval
Ethical Approval is not applicable for this article.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Research projects granted by Ministry of Agriculture & Rural Development (Code: 04/HĐ-KHCN).
Trial Registration
Not applicable, because this article does not contain any clinical trials.
Informed Consent
There are no human subjects in this article and informed consent is not applicable.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
