A therapeutic approach to target mitochondrial dysfunction using molecular docking studies: Screening of natural drugs for oral carcinoma
Manish Singh1, Manish Kumar Tripathi2, Alok Kumar Singh3, Chandra Shekhar Azad3, Indrajeet Singh Gambhir3, Brijesh Kumar1, Suresh Purohit1
1 Department of Pharmacology, Institute of Medical Sciences, Banaras Hindu University, Varanasi, Uttar Pradesh, India
2 Department of Biological Sciences and Engineering, Maulana Azad National Institute of Technology, Bhopal, Madhya Pradesh, India
3 Department of Medicine, Institute of Medical Sciences, Banaras Hindu University, Varanasi, Uttar Pradesh, India
|Date of Submission||13-Oct-2017|
|Date of Acceptance||06-Dec-2017|
|Date of Web Publication||28-Jun-2018|
Alok Kumar Singh
Department of Medicine, Institute of Medical Sciences, Banaras Hindu University, Varanasi - 221 005, Uttar Pradesh
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: Mitochondrial dysfunction is the major cause of various types of cancer, leading to death worldwide. The present study investigated the in silico binding potential of natural flavonoids and essential oils with human cyclophilin D (CyPD) protein. CyPD protein is a major molecular marker for apoptosis and has been reported to be elevated in oral carcinoma. Methods: PubChem database was used to check the efficacy of different active phytoconstituents (kaempferol, quercetin, eugenol, oxyresveratrol, tanshinone 2a, catechin, epicatechin, cinnamaldehyde, and emodin). These compounds were used as ligands to check their potential as anticancer agents against the inner mitochondrial membrane protein, CyPD. Docking studies were performed with the help of Discovery Studio 2.5 and Autodock. Emodin was used as a reference inhibitor to compare the results. Results: The binding energy (B.E.) of the reference inhibitor (known/established drug) emodin was observed −28.9 kcal/mol while novel inhibitors (catechin, cinnamaldehyde, epicatechin, eugenol, kaempferol, oxyresveratrol, quercetin, and tanshinone 2a) exhibited a range from −51.51 to −5.89 kcal/mol. Quercetin, kaempferol, and epicatechin (B.E.: −51.51, −34.79, and −30.62 kcal/mol, respectively) showed strong affinity as compared to reference inhibitor (B.E.: −28.9 kcal/mol). Conclusion: Quercetin, kaempferol, and epicatechin can be used as lead inhibitors against targeting CyPD.
Abbreviations used: CyPD: Cyclophilin D, BE: Binding Energy, PTPC: Permeability transition pore complex, mPTP: Mitochondrial permeability transition pore.
Keywords: Apoptosis, cyclophilin D pathway, mitochondrial dysfunction, molecular docking, oral cancer
|How to cite this article:|
Singh M, Tripathi MK, Singh AK, Azad CS, Gambhir IS, Kumar B, Purohit S. A therapeutic approach to target mitochondrial dysfunction using molecular docking studies: Screening of natural drugs for oral carcinoma. Phcog Mag 2018;14, Suppl S1:192-6
|How to cite this URL:|
Singh M, Tripathi MK, Singh AK, Azad CS, Gambhir IS, Kumar B, Purohit S. A therapeutic approach to target mitochondrial dysfunction using molecular docking studies: Screening of natural drugs for oral carcinoma. Phcog Mag [serial online] 2018 [cited 2022 Jan 17];14, Suppl S1:192-6. Available from: http://www.phcog.com/text.asp?2018/14/55/192/235280
- The present study investigates the in-silico binding potential of natural flavonoids and essential oils with human cyclophilin D (CyPD) protein
- It was observed that quercetin, kaempferol, and epicatechin showed strong affinity as compared to reference inhibitor emodin.
| Introduction|| |
Oral cancer is the major cause of cancer morbidity and mortality worldwide. The prevalence of oral cancer in Southeast Asia and in India is approximately 40% of total malignancies. The most affected oral subsite is buccal mucosa and alveolus in the developing countries.
Most of the cancer cells survive and proliferate by manufacturing ATP via glycolysis more willingly than oxidative phosphorylation even in the presence of sufficient oxygen as it was postulated by Warburg., Corresponding with this, increase in glycolysis is a suppression of mitochondrial activity in cancer. Under physiological conditions, mitochondria harbor a robust mitochondrial transmembrane potential and a low conductance state of the permeability transition pore complex (PTPC) might contribute to the exchange of small metabolites between the cytosol and mitochondrial matrix, a process that is mainly controlled by mitochondrial solute carriers. PTPC would be composed of voltage-dependent anion channel present in the outer membrane, adenine nucleotide translocase in the inner membrane, and cyclophilin D (CyPD) in the mitochondrial matrix.
Cell apoptosis or survival depends on the balance of pro-apoptotic or anti-apoptotic BCL-2 proteins primarily at the mitochondrial membrane. BCL-2 anti-apoptotic protein is considered as a primary target for anticancer drug development. There are reports for upregulated CyPD in many human cancers that may represent a suppression of apoptosis in cancer cells. It was reported that BCL-2 interacts with CyPD and may be significant for tumor development. A novel function of CyPD is that it serves as the main regulator of the mitochondrial permeability transition pore (mPTP). CyPD is proposed to exert an anti-apoptotic effect by binding to BCL-2, a key regulator of apoptosis. A CyPD inhibitor, cyclosporine-A which is a well-known drug, disrupts the CyPD–BCL-2 interaction and enhances the release of cytochrome c from mitochondria, which ultimately leads to apoptosis.
Naturally occurring active compounds have been used for the prevention and treatment of various diseases including cancer. These natural compounds are more beneficial than synthetic compounds due to less toxicity, more availability, and less expensive. Flavonoids and essential oils are a group of natural compounds having antioxidative, anti-inflammatory, and antitumor activity.
Computational and bioinformatics tools became a very useful resource for screening potential therapeutic agent for a particular protein receptor. It has been found that only limited reports of CyPD inhibition with natural active compound are available., Therefore, the present study aims to determine the possible interaction and binding-free energy of natural active compounds with CyPD protein using molecular docking approach.
| Materials and Methods|| |
Preparation of protein and active site identification
The protein selected for the present study is human CyPD (protein data bank [PDB] ID: 4ZSC) whose three-dimensional (3D) structure was obtained from PDB with crystal resolution 1.5 Ǻ. Water molecules in the crystal were removed and the crystal structure was prepared for docking with the help of prepare protein protocol of Discovery Studio 2.5. Active sites were identified with the help of PDB of Japan. The amino acid residues, namely ARG-97, GLN-105, GLY-114, ARG-124, ALA-143, ASN-144, THR-149, GLN-153, PHE-155, and HIS-168, were found to be present in the active sites of the selected protein. CASTp an online tool was also used to measure active site pockets on 3D protein structures, and the area and volume of active sites were found to be 267 Ǻ2 and 330.1 Ǻ3 respectively.
Preparation of ligand
Literature suggested that some of the herbal compounds have anticancer properties; therefore, the present study plans to use these compounds as ligands for the molecular interaction study. Ligands selected for the present study, namely catechin, cinnamaldehyde, epicatechin, eugenol, oxyresveratrol, quercetin, tanshinone 2a, kaempferol, and emodin, were retrieved from PubChem database (www.pubmed.com) as SDF format. Further, compounds were prepared for docking study with the help of “prepare ligand” protocol of Discovery Studio. The 2D structure of the selected ligands is shown in [Table 1].
Molecular docking protocol with Discovery Studio
Molecular docking was performed with the help of CDOCKER program of Discovery Studio. It allowed full ligand flexibility to the ligands which help in the generation of ligand conformation. The crystal structure was prepared from the “prepare protein” protocol of Discovery Studio with a default parameter. The selected phytoconstituents used as ligands were prepared in Discovery Studio with the help of “prepare ligand” protocol. Further, docking was performed using the CDOCKER protocol which is a grid-based molecular docking method that uses CHARMm force filed. The default parameters used for the docking study were top hits, 10; random conformations, 10; orientations to refine, 10; force field, CHARMm; and use full potential, false. The obtained results are evaluated on the basis of total docking energy, i.e., CDOCKER energy.
Molecular docking protocol Autodock 4.0
Further, docking study was also performed with the help of Autodock to validate our results obtained from Discovery Studio 2.5. The Lamarckian genetic algorithm was used to identify the binding modes and confirmation of ligands in the active site of the protein. Autodock tools were used to prepare the crystal structure of protein. Heteroatoms and water molecules were removed from the crystal structure of the protein, and polar hydrogen and charges were added to the macromolecule during the protein preparation. Ligands torsion was made rotatable to perform flexible docking. Grid maps were set around the active site residues of the receptor using AutoGrid with a grid-point spacing of 0.375 Ǻ and grid box measuring 51 × 51 × 55 Ǻ. Docking studies were performed using a Lamarckian genetic algorithm with a population size of 150 in combination with grid-based energy evaluation method for the calculation of grid maps. Docking parameters such as mutation rate, crossover rate, and population size were varied to perform docking studies. The best docked conformations were obtained with a mutation rate of 0.02, population size of 150, and crossover rate of 0.8. The results obtained were clustered with a root mean square tolerance of 1.0 Ǻ to obtain final docking results.
| Results and Discussion|| |
Cancer possesses a noteworthy health concern worldwide. Standard treatment strategy for cancer depends on the stage of the disease; more advanced carcinomas often require multimodality therapy with surgery, radiation, and chemotherapy, which can result in very high morbidity. Treatment modalities of cancer have numerous side effects; therefore, development of new treatment modalities such as herbal medicine is crucial for prevention and to reduce mortality.
Docking studies reveal the binding mode of the selected phytocompounds with the human CyPD protein and give an insight to the amino acid residues involved in the binding process. The details of the binding residues involved in the binding process, hydrophobic interaction pattern, their bond energy and their 2D binding forces are shown in [Figure 1](a-i) respectively are shown in [Table 2]. Emodin is docked and used as a standard/reference drug in this study and its binding energy (B.E.) was found to be −28.96 kcal/mol. Recently, Zhang et al. suggested that emodin is a potent inhibitor of CyPD through their in vitro study on HepG2 cells (liver cancer cell line) and also revealed that three hydrogen bonds exist in between CyPD–emodin complex through molecule docking. The docking results of all the selected phytoconstituents were compared with the emodin. Compounds namely epicatechin, quercetin, and kaempferol were found to show higher B.E. −30.62, −51.51, and −34.79 kcal/mol, respectively, as compared to emodin. The selected phytoconstituents were also docked using Autodock to perform cross-validation of the results. Docking results with Autodock also show that epicatechin, quercetin, and kaempferol show higher B.E. as compared to the other phytoconstituents [Table 2].
|Figure 1: T wo-dimensional interaction map of selected phyoconstituents with human cyclophilin D protein: (a) interaction of catechin with target protein, (b) interaction of cinnamaldehyde with target protein, (c) interaction of epicatechin with target protein, (d) interaction of eugenol with target protein, (e) interaction of kaempferol with target protein, (f) interaction of oxyresveratrol with target protein, (g) interaction of quercetin with target protein, (h) interaction of tanshinone 2a with target protein, (i) interaction of emodin with target protein|
Click here to view
|Table 2: Interaction of human cyclophilin-D proteins with selected phytoconstituents|
Click here to view
Docking study with epicatechin
Docking results obtained from the interaction of human CyPD protein with epicatechin reveal that it shows good binding affinity ΔG −30.62 kcal/mol with the receptor. The docking complex of epicatechin and human CyPD and their corresponding interacting residues are shown in [Figure 1]c. The amino acid residues GLY-114, THR-115, ARG-124, and ALA-143 showed that interaction with the phenolic oxygen of epicatechin and the bond length of this interaction was found to be 1.96, 2.40, 2.47, and 2.25 Ǻ, respectively. From the docking complex, it showed that amino acid residues, namely HIS-96, GLY-116, THR-149, SER-52, GLY-151, ASN-150, MET-142, and ASN-144, were found to be involved with Van der Waals interaction [Figure 1]c. Docking results from Autodock show that epicatechin interacts with the active site residue THR-149, and this binding and interaction pattern support the result of Discovery Studio docking results. B.E. of this interaction was found to be 5.40 kcal/mol which is greater as compared to the emodin ΔG −5.31 kcal/mol.
Docking study with quercetin
Docking result of quercetin with human CyPD protein reveals that it also shows good binding affinity ΔG −51.51 kcal/mol with the receptor. The dock complex and the residues involved in the binding process are shown in [Figure 1]g. The amino acid residues GLN-153 and ASN-144 showed interaction with the phenolic oxygen of quercetin; bond lengths of this interaction were found to be 2.61 and 3.06 Ǻ, respectively. Docking of quercetin with human CyPD showed that amino acid residues, namely GLN-105, SER-152, GLY-151, GLY-116, THR-115, and GLY-114, were found to be involved with Van der Waals interaction [Figure 1]g. For comparative study, we also performed docking from Autodock, the results obtained shows that quercetin interacts with the active site residue GLN-153, and B.E. of this interaction was found to be 6.4 kcal/mol.
Docking study with kaempferol
The dock complex of kaempferol with human CyPD protein and their interacting residues during the binding process is shown in [Figure 1]e. The B.E. obtained after the binding process ΔG −34.79 kcal/mol reveals that kaempferol also shows a good binding affinity for the receptor as compared to standard. In the binding of kaempferol with human CyPD amino acid residues, GLY-116 shows interaction with the phenolic oxygen of kaempferol with the bond length of 2.29 Ǻ. In the docking process, Van der Waals interactions are also formed and the amino acid residues involved in this interaction are GLY-114, THR-115, ARG-124, GLN-153, and GLN-105 [Figure 1]e. Docking study of kaempferol with Autodock also shows that it interacts with the GLY-114 residue of the active site and B.E. of this interaction was found to be Δ −5.83 kcal/mol.
Mitochondrial CyPD has appeared as a vital target in cancer chemotherapy because it plays a key role in mPTP, being an integral component of the mPTP. The immunosuppressant CsA is the CyPD inhibitor that has effectively been used for studying the role of CyPD in mitochondria-mediated cell death. CyP-D is proposed to exert an anti-apoptotic effect by binding to BCL-2. The present study revealed that the inhibitors quercetin, kaempferol, and epicatechin have higher binding potency than emodin; therefore, these drugs may disrupt the CyPD–BCL-2 interaction which may ultimately enhance the apoptotic process.
| Conclusion|| |
Our finding showed that quercetin, epicatechin, and kaempferol have potential to bind with human CyPD protein, which is a mitochondrial matrix protein and plays a very important role in mitochondrial membrane pore opening. The obtained results reveal that quercetin, epicatechin, and kaempferol have better binding ability than the standard drug emodin against human CypD. Therefore, it is worth mentioning that these phytoconstituents could serve as potential anticancer drugs through mitochondrial dysfunction.
This study was financially supported by DST-SERB Young Scientist Scheme [No. YSS/2014/000297] from the DST. We thank the Institute of Medical Sciences, Banaras Hindu University, Varanasi, India, for providing the infrastructure, Ph. D. fellowship and other research facilities.
Financial support and sponsorship
This study was financially supported by DST-SERB Young Scientist Scheme [No. YSS/2014/000297] from the DST.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Jain V, Dharkar D, Nandini H, Jain SM. Overall treatment duration with multimodality treatment approach and outcome results in oral cavity carcinoma: A study from an institute from central India. J Cancer Res Ther 2017;13:351-5.
Pathak KA, Das AK, Agarwal R, Talole S, Deshpande MS, Chaturvedi P, et al.
Selective neck dissection (I-III) for node negative and node positive necks. Oral Oncol 2006;42:837-41.
Giang AH, Raymond T, Brookes P, de Mesy Bentley K, Schwarz E, O'Keefe R, et al.
Mitochondrial dysfunction and permeability transition in osteosarcoma cells showing the Warburg effect. J Biol Chem 2013;288:33303-11.
Warburg O. On the origin of cancer cells. Science 1956;123:309-14.
Fulda S, Galluzzi L, Kroemer G. Targeting mitochondria for cancer therapy. Nat Rev Drug Discov 2010;9:447-64.
Wenner CE. Targeting mitochondria as a therapeutic target in cancer. J Cell Physiol 2012;227:450-6.
Muhammad SA, Fatima N. In silico
analysis and molecular docking studies of potential angiotensin-converting enzyme inhibitor using quercetin glycosides. Pharmacogn Mag 2015;11:S123-6.
Chen G, Izzo J, Demizu Y, Wang F, Guha S, Wu X, et al.
Different redox states in malignant and nonmalignant esophageal epithelial cells and differential cytotoxic responses to bile acid and honokiol. Antioxid Redox Signal 2009;11:1083-95.
Zhang L, He D, Li K, Liu H, Wang B, Zheng L, et al.
Emodin targets mitochondrial cyclophilin D to induce apoptosis in HepG2 cells. Biomed Pharmacother 2017;90:222-8.
Kinjo AR, Suzuki H, Yamashita R, Ikegawa Y, Kudou T, Igarashi R, et al.
Protein Data Bank Japan (PDBj): Maintaining a structural data archive and resource description framework format. Nucleic Acids Res 2012;40:D453-60.
Dundas J, Ouyang Z, Tseng J, Binkowski A, Turpaz Y, Liang J, et al.
CASTp: Computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues. Nucleic Acids Res 2006;34:W116-8.
Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, et al
. AutoDock. Ver. 4.0.1. La Jolla: The Scripps Research Institute; 2007.
Aghbali A, Hosseini SV, Delazar A, Gharavi NK, Shahneh FZ, Orangi M, et al.
Induction of apoptosis by grape seed extract (Vitis vinifera
) in oral squamous cell carcinoma. Bosn J Basic Med Sci 2013;13:186-91.
[Table 1], [Table 2]