Original Research ARTICLE
The Combined Use of in Silico, in Vitro, and in Vivo Analyses to Assess Anti-cancerous Potential of a Bioactive Compound from Cyanobacterium Nostoc sp. MGL001
- 1Laboratory of Microbial Genetics, Department of Botany, Banaras Hindu University, Varanasi, India
- 2Biochemistry and Molecular Biology Laboratory, Department of Zoology, Banaras Hindu University, Varanasi, India
Escalating incidences of cancer, especially in developed and developing countries, demand evaluation of potential unexplored natural drug resources. Here, anticancer potential of 9-Ethyliminomethyl-12-(morpholin-4-ylmethoxy)-5,8,13,16-tetraaza -hexacene-2,3-dicarboxylic acid (EMTAHDCA) isolated from fresh water cyanobacterium Nostoc sp. MGL001 was screened through in silico, in vitro, and in vivo studies. For in silico analysis, EMTAHDCA was selected as ligand and 11 cancer related proteins (Protein Data Bank ID: 1BIX, 1NOW, 1TE6, 2RCW, 2UVL, 2VCJ, 3CRY, 3HQU, 3NMQ, 5P21, and 4B7P) which are common targets of various anticancer drugs were selected as receptors. The results obtained from in silico analysis showed that EMTAHDCA has strong binding affinity for all the 11 target protein receptors. The ability of EMTAHDCA to bind active sites of cancer protein targets indicated that it is functionally similar to commercially available anticancer drugs. For assessing cellular metabolic activities, in vitro studies were performed by using calorimetric assay viz. 3-(4,5-dimethylthiazol-2-yl)-2,5 diphenyltetrazolium bromide (MTT). Results showed that EMTAHDCA induced significant cytotoxic response against Dalton's lymphoma ascites (DLA) cells in a dose and time dependent manner with an inhibitory concentration (IC50) value of 372.4 ng/mL after 24 h of incubation. However, in case of normal bone marrow cells, the EMTAHDCA did not induce cytotoxicity as the IC50 value was not obtained even with higher dose of 1,000 ng/mL EMTAHDCA. Further, in vivo studies revealed that the median life span/survival days of tumor bearing mice treated with EMTAHDCA increased significantly with a fold change of ~1.9 and 1.81 corresponding to doses of 5 and 10 mg/kg body weight (B.W.) of EMTAHDCA respectively, as compared to the DL group. Our results suggest that 5 mg/kg B.W. is effective since the dose of 10 mg/kg B.W. did not show any significant difference as compared to 5 mg/kg B.W. Taken together, our findings based on in silico, in vitro, and in vivo analyses suggest that EMTAHDCA has potential anticancer effects, and thus, can be considered for cancer treatment.
Cancer is one of the world's deadliest diseases and is becoming a matter of worldwide concern. According to the World Health Organization (WHO), 7.9 million deaths (around 13% of all deaths) have been accounted in 2007 mainly due to the cancer; of which, 38% cancer incidences have been observed in developed countries, and 68% cases have been notified in developing countries (http://www.who.int/cancer/en/). A project of the International Agency for Research on Cancer (IARC), known as GLOBOCAN 2012, has estimated that more than 13.2 million deaths and nearly 21.4 million new cancer incidences will appear by 2030 (Bray et al., 2013; Ferlay et al., 2013). Escalating cancer evidences, ineffectiveness of anticancer drugs, and adverse side-effects of drug over-dose demand more investigations on identifying newer, broad-spectrum anticancer molecules from less explored, potential natural products. Development of modern technologies offers unique opportunities for finding and discovering new drug compounds from natural resources (Cragg et al., 1997). Over the past few years, natural products and their derivatives have been recognized as important sources of new medicines and therapeutic agents (Cragg and Newman, 2013). da Rocha et al. (2001) have mentioned about the importance of natural products derived from microbes in anticancer therapy. Various microorganisms viz. actinomycetes, fungi, and bacteria including cyanobacteria are known for producing bioactive metabolites (Bérdy, 2005; El-Elimat et al., 2012). Of these organisms, cyanobacterial cultivation without organic substrates is considered as a cost effective approach for scientific studies as compared to other microorganisms (Bullerjahn and Post, 2014; Dias et al., 2015).
Cyanobacteria are the most diverse group of oxygenic photosynthetic prokaryotic organisms distributed across a broad range of ecosystem (Singh et al., 2005; Sivonen and Börner, 2008). Nowadays, cyanobacteria have become a potential microorganism for biotechnological applications for producing wide array of unprecedented biologically active natural products with anticancer, antifungal, antimicrobial, antiviral, antioxidants, and anti-inflammatory effects (Gul and Hamam, 2005; Mayer and Hamann, 2005; Singh et al., 2005; Prasanna et al., 2010; Harnedy and FitzGerald, 2011). Several scientific studies have demonstrated that cyanobacteria are emerging source for drug discovery (Patterson et al., 1994; Khairy and El-Kassas, 2010; Kumar et al., 2010). Recently, Lee et al. (2014) have reported that multifunctional photo-luminescent green carbon nanotags (G-Tags) synthesized from harmful cyanobacteria has been used as a drug delivery system in cancer therapy. Some fresh water filamentous cyanobacterial strains are also considered as important sources of anticancer biomolecules (Srivastava et al., 2015).
Among different Nostocales species, Nostoc sp. is known to produce more than 100 different natural products with antifungal, antimicrobial, antimalarial, and anticancer activities (Shanab et al., 2012). Cryptophycin1 isolated from cyanobacterium Nostoc sp. GSV224 has shown to exhibit potent anticancer activity against Adriamycin resistant M17 breast cancer and DMS 273 lung cancer cell lines, KB human nasopharyngeal cancer cells, and LoVo human colorectal cancer cells (Patterson et al., 1991). Several chemically synthesized and naturally occurring analogs of cryptophycin are under second generation clinical trials as potential future anticancer drugs (Liang et al., 2005). Cyanobacterium Nostoc sp. is also known to produce bioactive secondary metabolites having antiviral and antitumor activities (Dembitsky and Rezanka, 2005). Nostodione A is a potent antimitotic compound first time isolated from Nostoc commune (Kobayashi et al., 1994). Some of the marine cyanobacteria belonging to this particular Nostoc genera are Nostoc linckia, N spongiaeforme var. tenue, and Nostoc sp. ATCC 53789 and GSV 224 are well-recognized sources of potent cytotoxic compounds working against human tumor cell lines (Burja et al., 2001).
In our previous study, we have isolated and characterized novel bioactive compound i.e. Ethyliminomethyl-12-(morpholin-4-ylmethoxy)-5,8,13,16-tetraaza-hexacene-2,3-dicarboxylic acid (EMTAHDCA) from cyanobacterium Nostoc sp. MGL001 and proved its application as an antibacterial agent against multidrug resistant gram negative bacteria. The compound has been observed to produce maximum zone of inhibition at 150 μg/mL concentration (Niveshika et al., 2016). In this study, we intended to evaluate the potential application of EMTAHDCA as an anticancer drug resource apart from antibacterial activity through in silico, in vivo, and in vitro studies.
Materials and Methods
In Silico Study
Structural elucidation of EMTAHDCA was performed, as described previously (Niveshika et al., 2016). Three dimensional crystal structures of 11 cancer target receptors with PDBID of 1BIX (Gorman et al., 1997), 1NOW (Mark et al., 2003), 1TE6 (Chai et al., 2004), 2RCW (Madej et al., 2014), 2UVL (Herman et al., 2009), 2VCJ (Brough et al., 2008), 3CRY (Oakley et al., 2010), 3HQU (Chowdhury et al., 2009), 3NMQ (Yun et al., 2011), 5P21 (Pai et al., 1990), and 4B7P (Fogliatto et al., 2013) were retrieved from the protein databank (PDB) (www.rcsb.org/pdb; Table S1). Discovery Studio 3.1 tool was used to optimize 3D models of ligand (EMTAHDCA) as well as 11 cancer related proteins (Gao and Huang, 2011). Preparation of terminus of each proteins were not require to perform. Again discovery studio version 3.1 software was used to remove heteroatom, water molecules, and ligands before docking (Gao and Huang, 2011). MetaPocket 2.0 was used for predicting the active site of a target protein (http://projects.biotec.tu-dresden.de/metapocket/; Huang, 2009; Zhang et al., 2011). An advanced docking tool, Yet Another Scientific Artificial Reality Application (YASARA), was used for docking calculation between receptors and ligand (Krieger and Vriend, 2014). The target protein was set in YASARA to run macro file (dock_run.mcr). YASARA structure provides Autodock and VINA tools to dock ligands with proteins (http://www.yasara.org/docking.htm) using a very simple approach (Morris et al., 1998; Trott and Olson, 2010). AMBER force field was used by docking software YASARA for calculation of energy. AMBER force-field parameters were adjusted iteratively so that the damage done to the structures was minimal (Krieger et al., 2004). Energy minimization was done using default parameter of YASARA Autodock Vina tool (http://www.yasara.org/yamber.htm).
In Vitro Study
For all experiments, 15–20 weeks old healthy adult mice (AKR strain) were used because these mice are highly susceptible to develop tumor and have a shorter life-span. Under aseptic condition, about 1 × 106 DLA cells were transplanted intraperitoneally in mice. The life-span of DL bearing mice was 20 ± 2 days, however the life-span of healthy mice is about 18 months, as reported by Bharti and Mishra (2008). The success rate of developing cancer using DLA cells is found to be 100% (Goldie and Felix, 1951). Experimental mice were maintained with standard feeding and drinking at the temperature of 25 ± 2°C for 12/12 h light/dark cycle throughout the study period in the animal house facility of the Zoology Department, Institute of Science, Banaras Hindu University, Varanasi, India as per guidelines of the Institutional Animal Ethical Committee (IAEC), Banaras Hindu University (Bharti and Mishra, 2011).
After successful development of the tumor, mice were anesthetized with chloroform for collecting ascites fluids from the peritoneal cavity. The fluid was then centrifuged at 3,000 rpm for 6 min, and cell pellet was taken and re-suspended in 1 mL of 1X Phosphate Buffered Saline (PBS). Next, 100 μL of cell suspension was mixed with 100 μL of 2X trypan blue (Sigma) to check viable cell density using hemocytometer. The formula used to count viable cells is as follows: No. of cells × dilution factor/ volume (0.4 mm3 i.e., 4 squares were counted, each with area of 1 mm2, and depth of counting chamber of 0.1 mm).
A preliminary study of DLA cell viability under different concentrations of EMTAHDCA (100, 250, 500, 750, and 1,000 ng/mL) were performed. The 1 mg/mL concentration of EMTAHDCA was prepared by dissolving 1 mg of purified EMTAHDCA in 1 mL of 1% DMSO. Different concentrations of EMTAHDCA were prepared separately using 1 mg/mL concentration as stock solution by formula C1V1 = C2V2 where C1 represent the stock concentration i.e., 1 mg, C2 represent the working concentrations i.e., 100, 250, 500, 750, and 1,000 ng, V1 represent the required volume of stock concentration to make working concentration in desired volume (V2).
Approximately, 1 × 106 cells were seeded per well in 6-well plate using complete growth medium (RPMI+10% FBS), then cells were treated with different concentrations of EMTAHDCA and incubated for 24 h inside the CO2 incubator with 5% CO2. Cells treated with only medium served as a control group. Cell imaging was done under EVOS FLoid Cell Imaging Station.
Cytotoxicity of EMTAHDCA at various concentrations (100, 250, 500, 750, and 1,000 ng/mL) were assessed by using MTT assay in DLA cells and normal mouse bone marrow cells (BMC). BMC were isolated from femur bones of 15–20 weeks old adult mice. The mice was anesthetized using chloroform for collecting bone marrow. After that with the help of 24-gauge needle, bone marrow was flushed with pre-warmed PBS and single cells suspension was prepared. The cells were suspended in RPMI-1640 medium supplemented with antibiotic solution and 10% FBS. Cells were seeded in culture plates and maintained in 5% CO2 at 37°C for 24 h (Prasad and Koch, 2014).
For MTT assay, 1 × 106 cells/mL DLA and BMC (positive control) were seeded per well in 96-well plate and treated with different concentrations (100–1,000 ng/mL) of EMTAHDCA along with the control (only medium) and vehicle control (1% DMSO). After 24 h, supernatant was removed from each well, washed twice with PBS and MTT solution (5 mg ml−1 in PBS) with 100 μL of medium was added. After incubation for another 4 h, resultant formazan crystals were dissolved in dimethyl sulfoxide (100 μL), and the absorbance was measured by a microplate reader (Bio-RAD 680, USA) at 590 nm with a reference wavelength of 620 nm (Prasad and Koch, 2014) and IC50 value of EMTAHDCA was calculated using GraphPad Prism 7 software. All experiments were performed in quadruplicate, and the relative cell viability (%) was expressed as a percentage relative to control group.
In case of time dependent study, cells were exposed to 250, 500, and 750 ng/mL of EMTAHDCA along with control and vehicle control and incubated for 12 and 24 h, respectively, and the cell viability test was performed as described above.
In Vivo Study
Overall survival time and body weight of DL bearing mice after treatment with EMTAHDCA were evaluated to check anti-cancer and cytotoxic effect of EMTAHDCA in DL induced mice.
DL transplanted mice were randomly divided into two groups, DL group (DL) and DL+ EMTAHDCA treated group. EMTAHDCA treated groups were further divided into two subgroups, 5 and 10 mg/kg B.W. groups, and were administered with respective doses at 6th, 12th, and 18th day of DL transplantation. For the survival time analysis, the percent of mean survival time of EMTAHDCA treated mice (T) over mean survival time of the control group (C) mice were calculated followed by the evaluation of antitumor activities of EMTAHDCA by computing the T/C value of respective groups. The DL group mice were taken as control group for survival time analysis. The body weight was measured at day 0, 3, 6, 9, 12, 15, and 18 because an increase in body weight in ascites tumor model represents tumor mass as explained by Sharma and Koch (2016). Briefly, male mice of 22 g B.W. were taken for the experiment and considered as blank. The graph was plotted after subtracting the body weight of DL mice and 5 and 10 mg/kg B.W. of EMTAHDCA treated DL mice with that of the healthy mice. Kaplan–Meier survival curve was plotted and analyzed by GraphPad Prism7 software using log-rank analysis to examine the level of significance, and p < 0.05 was obtained by comparing between normal mice represented as control, DL induced mice and EMTAHDCA treated DL mice groups.
Statistical analysis were performed using software SPSS (version 16.0). For each experiment, six adult male mice were taken in each group (n = 6). Each experiments were done 4 times independently. Results were expressed as mean ± standard error of mean (SEM). The IC50 value of compound was calculated using GraphPad Prism 7 software. Data were analyzed by chi-square test and One-way ANOVA followed by Tukey post-hoc test. A p < 0.05 was considered as statistically significant.
Molecular docking was successfully performed between the selected ligand (EMTAHDCA) and 11 cancer protein targets (PDB ID: 1BIX, 1NOW, 1TE6, 2RCW, 2UVL, 2VCJ, 3CRY, 3HQU, 3NMQ, 5P21, and 4B7P; Table S1) using YASARA software. The active site of cancer protein targets was determined using Metapocket. Summary of molecular docking results of selected ligand with cancer target protein is represented in Table S2. EMTAHDCA interacted efficiently with the cancer target 3CRY with binding energy of 9.2 and dissociation constant of 169507.6. Contracting receptor residues involved in their interaction through the formation of electrostatic bond, van der Waals force, covalent bond, and hydrogen bond were Gly23, Ser24, Asn25, Arg30, Arg34, Gly58, Lys59, Thr60, Ser61, Gln62, Thr63, Trp64, Ile68, Gln97, Gln98, Met104, Tyr105, and Tyr139 (Figure 1).
Figure 1. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (3CRY), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) The 2D level interation clearly mentioned that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
Moreover, EMTAHDCA efficiently interacted with the cancer target 2VCJ with a good binding energy of 8.5 and less dissociation constant of 567576.2. In this case, active site 1 residues viz. Leu48, Asn51, Ser52, Ala55, Lys58, Tyr61, Glu62, Asp93, Ile96, Gly97, Met98, Asp102, Leu107, Gly108, Phe138, Val150, Thr184, and Val186 were found to be involved in the interaction. Residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding (Figure 2).
Figure 2. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (2VCJ), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) The 2D level interation depicted that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
Ligand EMTAHDCA also interacted very efficiently with the cancer protein target 1TE6 with a good binding score of 8.4 and lowest dissociation constant [pM] of 647439.1. Active site 1 residues viz. Lys27, Gly28, Leu29, Phe30, Arg31, Ala32, Lys119, Ala122, Ala123, Glu126, Leu127, Pro128, Glu376, Asp377, Thr378, and Phe379 were involved in their interaction (Figure 3).
Figure 3. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (1TE6) (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) The 2D level interation indicated that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
In case of interaction between EMTAHDCA and 2RCW, the binding score was found to be 8.4, whereas the dissociation constant was 689160.3. Amino acid residues involved in this interaction were Lys42, Val112, Ser115, Leu116, Gly119, Gly120, Ser121, Val131, Glu134, Lys135, Lys137, Arg180, Glu181, Gly182, Glu183, Arg186, Ser213, Gln214, and Ile333. In addition, 2D view of ligand and protein interaction depicted very clearly that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding (Figure 4).
Figure 4. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (2RCW), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) Residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding depicted through 2D level interation.
Interaction between EMTAHDCA and 1BIX secured binding score of 8.2 and dissociation constant of 1017760.4. Amino acid residues involved in their interaction were Glu96, Tyr128, Arg156, Tyr171, Asn174, Gly176, Arg177, Gly178, Arg181, Asn222, Asn226, Asn229, Phe266, Trp267, Thr268, Met271, Ala273, Val278, Gly279, and Trp280. These residues were interacted through electrostatic, van der Waals, covalent, and hydrogen bonding (Figure 5).
Figure 5. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (1BIX), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) Amino acid residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding clearly indicated by their 2D level interation.
EMTAHDCA interacted with 1NOW at its active site with a high binding score of 8.2 and low dissociation constant of 1024654.8. Active site residues were Ile143, Ser144, Ser145, Asp146, Glu147, Ser148, Ser197, Pro198, Arg199, Phe200, Ser201, Arg203, Pro442, Leu469, Gln475, Lys478, Gln479, Leu480, Phe481, Ile482, and Lys517 (Figure 6).
Figure 6. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (1NOW), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) The 2D level interation of ligand and protein target clearly mentioned that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
In case of EMTAHDCA and 5P21 interaction, binding score was 8.2 and dissociation constant was 1608019.0. Active site 1 residues were Gly12, Gly13, Val14, Gly15, Lys16, Ser17, Ala18, Val29, Asp30, Glu31, Tyr32, Asp33, Pro34, Gln61, Asn85, Asn86, Thr87, Lys88, Lys117, and Thr124. The 2D view clearly mentioned that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding (Figure 7).
Figure 7. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (5P21), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) Residues were interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
In case of 3NMQ, EMTAHDCA interacted with a binding score of 7.8 and dissociation constant of 1955781.4. Active site residues viz. Asn51, Ala55, Lys58, Ile96, Gly97, Met98,Asp102, Leu107, Ile110, Ser113, Gly114, Phe134, Gly135, Val136, Phe138, Val150, His154, Thr184, and Val186 interacted through electrostatic, van der Waals, covalent, and hydrogen bonding (Figure 8).
Figure 8. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (3NMQ), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) 2D level interaction indicated that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
Cancer target protein 4B7P interacted with EMTAHDCA with positive energy of 7.8 and low dissociation constant of 1975688.0. Residues such as Arg46, Glu47, Ser50, Asn51, Ser53, Asp54, Asp57, Arg60, Lys112, Gly132, Gln133, Gly135, Val136, Gly137, Phe138, Phe213, Ile214, Gly215, and Tyr216 were involved in this interaction (Figure 9).
Figure 9. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (4B7P), (B) Poses of docked complexes, green color in sphere indicates prominent active site where the ligand interacted, (C) Residues were interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
Interaction between EMTAHDCA and 2UVL secured a binding score of 7.8 and dissociation constant of 2081812.8. Active site amino acid residues such as Arg244, Thr246, Val247, Tyr282, Tyr283, Val284, Gly285, Asn286, Asp289, Lys291, Gly298, Leu299, and Arg300 were involved in their interaction. Electrostatic, van der Waals, covalent, and hydrogen bond were involved in this interaction as clearly mentioned in the 2D view (Figure 10).
Figure 10. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (2UVL), (B) Poses of docked complexes, green color in sphere indicates prominent active site where the ligand interacted, (C) Amino acid residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding clearly indicated by their 2D level interation.
Interaction between EMTAHDCA and 3HQU occurred with positive energy of 7.6 and low dissociation constant of 2645629.8. Amino acid residues such as Leu188, Pro189, Ala190, Leu193, Tyr197, Val210, Asp211, Asp212, Phe213, Leu214, Gly215, Lys216, Glu217, Thr218, Ser275, Asp278, and Leu279 were involved in their interaction using electrostatic, van der Waals, covalent, and hydrogen bonding (Figure 11).
Figure 11. Poses of docked complexes (A) In silico interaction between selected ligand (EMTAHDCA) and cancer target protein (3HQU), (B) Green color in sphere indicates prominent active site where the ligand interacted, (C) The 2D level interation depicted that residues interacted through electrostatic, van der Waals, covalent, and hydrogen bonding.
In Vitro and in Vivo Studies
In order to test the anticancer activity of EMTAHDCA, DLA cells were initially cultured with different concentrations of EMTAHDCA i.e., 100, 250, 500, 750, and 1,000 ng/mL for 24 h (Figures 12A–F) and observed to have cytotoxic effect. Cytotoxic effects of different concentrations of EMTAHDCA on DLA cells and normal BMC were also examined through MTT assay. The percent (%) viability of DLA cells were observed to be 100, 99.8, 97.7, 92.5, 51.2, 46.6, 40.2 at control, vehicle control, 100, 250, 500, 750, and 1,000 ng/mL, respectively whereas percent (%) viability in case of normal BMC were found to be 100, 99.7, 98.9, 98.9, 98.8, 98.7, and 98.7 at control, vehicle control, 100, 250, 500, 750, and 1,000 ng/mL, respectively (Figure 13A, Table S3). The calculated IC50 value was found to be 372.4 ng/mL EMTAHDCA against DLA cells. Further, in order to determine time dependent cytotoxicity of EMTAHDCA, MTT assay was performed using 250, 500, and 750 ng/mL concentration of EMTAHDCA at 12 and 24 h time points. The cell viability of DLA cells treated with 500 ng EMTAHDCA was observed 85% after 12 h and 52.5% after 24 h of incubation as compared to the control (100%). However, treatment with 250 ng/mL EMTAHDCA showed 95 and 94% cell viability after 12 and 24 h, respectively. Cell viability of DLA cells treated with 750 ng/mL EMTAHDCA decreased from 89% after 12 h to 45.6% after 24 h (Figure 13B, Table S4).
Figure 12. (A) Control slide representing DLA cells cultured for 24 h at 37°C with 5% CO2 in CO2 incubator. (B–F) Effect of different concentrations of EMTAHDCA (Control, 100, 250, 500, 750, and 1,000 ng/mL, respectively) on the viability of DLA cells after incubation of 24 h at 37°C with 5% CO2 in CO2 incubator. DC, Dead cells; VC, Viable cells.
Figure 13. Cell viability was determined by MTT assay, (A) Graph represents the cytotoxicity profile of EMTAHDCA against DLA cells and normal mice BMC at different concentrations (100, 250, 500, 750, and 1,000 ng/mL) along with control and vehicle control after 24 h of incubation, (B) Effect of EMTAHDCA on the DLA cell viability post 12 and 24 h treatment, respectively.
In vivo results showed a significant inhibition of tumor growth in mice of EMTAHDCA treated group as compared to the DL group. More importantly, median life span/survival days of tumor bearing mice treated with 5 and 10 mg/kg B.W. of EMTAHDCA significantly increased to 190 and 181%, respectively as compared to DL group (Figure 14A). Increased longevity of EMTAHDCA treated mice as compared to the DL mice was further verified by Kaplan's Meier survival curve using log-rank statistics which showed a significant increase in survival, ~1.9- and 1.81-fold with doses of 5 and 10 mg/kg B.W. of EMTAHDCA respectively as compared to DL mice (Figure 14B). An increase in body weight of DL mice was observed 6 days after the intraperitoneal injection of DLA cells until the death of mice. In contrast, a significant decrease in body weight was observed in DL mice treated with 5 and 10 mg/kg B.W. of the compound on 6th, 12th, and 18th days as compared to DL group (Figure 15). On 9th day, the body weight decreased to 5.8 and 5.7 g after the treatment with 5 and 10 mg/kg B.W. of the compound, respectively, as compared to DL group (7.2 g body weight). On 18th day, the weight decreased to 15 and 17 g after the treatment with 5 and 10 mg/kg B.W. of the compound, respectively, as compared to DL group (21 g body weight). We also found that the administration of EMTAHDCA in healthy mice at a same dose did not cause any major side effects.
Figure 14. (A) Graph represents mean survival/lifespan of DL mice and DL mice treated with EMTAHDCA at two different doses i.e., 5 and 10 mg/kg B.W. The treated/control (T/C) values are shown in percentage (%), (B) Kaplan–Meier survival curve for survival of healthy mice (control), DL mice, and DL mice treated with EMTAHDCA at two different doses i.e., 5 and 10 mg/kg B.W. Graph was plotted and analyzed by GraphPad Prism 7 software with log-rank analysis to examine the level of significance and P < 0.05 was obtained in between DL and DL treated with EMTAHDCA groups.
Figure 15. Graph represents effect of two different doses of EMTAHDCA i.e., 5 and 10 mg/kg B.W. on increase in body weight of DL mice from day 0 to 18 after DL transplantation. Values are expressed as mean ± S.E. obtained from 4 independent sets of experiments.
Cyanobacteria are known to produce a broad spectrum of metabolites with diverse activities. Since a significant portion of important cyanobacterial drug candidates is still underexplored, it is establishing a new research area of the pharmaceutical science. In our previous study, we reported for the first time that the bioactive compound, EMTAHDCA, from fresh water cyanobacterium Nostoc sp. MGL001 has antibacterial activity (Niveshika et al., 2016). Now, this study was designed using in silico, in vivo, and in vitro approaches to investigate whether this particular natural product also has anticancer effects.
In silico docking approach is considered as an emerging field for the rational design of drugs, and therefore, gaining significance in the field of pharmaceutical science (Terstappen and Reggiani, 2001). Here, in silico approaches were adopted in order to predict whether this compound has anticancer potential or not. Eleven cancer target proteins were selected (PDBID: 1BIX, 1NOW, 1TE6, 2RCW, 2UVL, 2VCJ, 3CRY, 3HQU, 3NMQ, 5P21, and 4B7P) based on their affiliation with different cancer diseases such as breast cancer, brain tumor, ovarian cancer, colorectal cancer, skin cancer, liver cancer, lung cancer, colon cancer, and bladder cancer (Table S1). Prominent active sites of selected cancer protein targets, which are generally used by commercially available anticancer drugs for binding, have already been reported (Gorman et al., 1997; Mark et al., 2003; Brough et al., 2008; Chowdhury et al., 2009; Herman et al., 2009; Oakley et al., 2010; Yun et al., 2011; Fogliatto et al., 2013; Madej et al., 2014).
In this study, comparative molecular docking calculation indicate that selected ligand (EMTAHDCA) interacts with good binding affinities and best positive energies with all the 11 cancer target proteins at their prominent active sites (predicted by Metapocket) that are generally used by commercially available anticancer drugs for binding. More positive energies indicate stronger binding, and negative energies mean no binding (Morris et al., 1998). Therefore, with the help of advanced molecular docking tools, it was predicted that EMTAHDCA has the ability to work as an anticancer drug agent.
Prediction of anticancer drug potential of EMTAHDCA through in silico approaches was further validated using in vitro and in vivo studies. For the in vitro study, DLA cells were selected because of their similarities with some lymphoma and leukemia of human origin (Kumar and Singh, 1996). DLA cells can be maintained by serial transplantation of ascites cells in the intraperitoneal cavity of one mouse to other mouse. To check cytotoxic effect of EMTAHDCA on DLA cells, MTT assay was performed. MTT is one of the reliable and sensitive colorimetric techniques, which is generally used for measuring in vitro cytotoxic effects of drugs on cancer cell lines and assessing the viability and proliferation of cancer cells (Meerloo et al., 2011). Five different doses of EMTAHDCA were selected to find out the most effective dose at a minimum concentration. Therefore, we started with the concentration ranging from 100 to 1,000 ng/mL for in vitro studies. Cytotoxic effects of EMTAHDCA observed using MTT assay on DLA cells indicate dose and time dependent activities. Our results on in vitro studies indicate that treatment with EMTAHDCA cause death of cancerous DLA cells and could act as cytotoxic agent in concentration and time dependent manner with an IC50 value at 372.4 ng/mL concentration against DLA cells without affecting normal bone marrow cells even with higher dose of 1,000 ng/mL EMTAHDCA. Since, we do not observed any cytotoxic effect of EMTAHDCA on bone marrow cells in concentration dependent study for 24 h, therefore, we have not taken positive control i.e., bone marrow cells for time dependent study.
In vivo studies showed significant increase in median lifespan/survival days of 5 mg/kg B.W. EMTAHDCA-treated mice and decrease in body weight as compared to 10 mg/kg B.W. EMTAHDCA treated and DL mice. However, the drug did not show any significant adverse side-effect on normal mice. The decrease in body weight and increased survival of EMTAHDCA treated mice are in agreement with earlier reports (Ajith and Janardhanan, 2003; Koiri et al., 2009; Singh and Trigun, 2013; Sharma and Koch, 2016) and could be due to regression of tumor growth in the presence of EMTAHDCA as compared to DL mice.
On the whole, our findings suggest that the selected compound EMTAHDCA has potential anticancer properties.
Our findings based on in silico, in vitro, and in vivo studies prove that bioactive compound EMTAHDCA isolated from fresh water cyanobacterium Nostoc sp. MGL001 has anticancer potential. Therefore, it is worth mentioning that EMTAHDCA can serve as an anticancer drug resource.
N and AM: designed the experiments; N and SM: performed the experiments; N, EV, AM, SM, and RM: analyzed the data; N and SM: wrote the manuscript and AM critically reviewed the paper.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Head, Department of Botany, Banaras Hindu University, Varanasi, India is gratefully acknowledged for providing laboratory facilities. We are also thankful to the Prof. Sukh Mahendra Singh, School of Biotechnology, Banaras Hindu University for providing YASARA software facility. Two of us (N and EV) are thankful to the UGC, New Delhi for financial support in the form of SRF.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2017.00873/full#supplementary-material
Ajith, T. A., and Janardhanan, K. K. (2003). Cytotoxic and antitumor activities of a polypore macrofungus, Phellinus rimosus (Berk) Pilat. J. Ethnopharmacol. 84, 157–162. doi: 10.1016/S0378-8741(02)00292-1
Bharti, B., and Mishra, R. (2011). Isoforms of Pax5 and co-regulation of T- and B-cells associated genes influence phenotypic traits of ascetic cells causing Dalton's lymphoma. Biochim. Biophys. Acta 1813, 2071–2078. doi: 10.1016/j.bbamcr.2011.08.003
Brough, P. A., Aherne, W., Barril, X., Borgognoni, J., Boxall, K., Cansfield, J. E., et al. (2008). 4, 5-Diarylisoxazole Hsp90 chaperone inhibitors: potential therapeutic agents for the treatment of cancer. J. Med. Chem. 51, 196–218. doi: 10.1021/jm701018h
Burja, A. M., Banaigs, B., Abou-Mansour, E., Burgess, J. G., and Wright, P. C. (2001). Marine cyanobacteria—a prolific source of natural products. Tetrahedron 57, 9347–9377. doi: 10.1016/S0040-4020(01)00931-0
Chai, G., Brewer, J. M., Lovelace, L. L., Aoki, T., Minor, W., and Lebioda, L. (2004). Expression, purification and the 1.8 Å resolution crystal structure of human neuron specific enolase. J. Mol. Biol. 341, 1015–1021. doi: 10.1016/j.jmb.2004.05.068
Chowdhury, R., McDonough, M. A., Mecinovic, J., Loenarz, C., Flashman, E., Hewitson, K. S., et al. (2009). Structural basis for binding of hypoxia-inducible factor to the oxygen-sensing prolyl hydroxylases. Structure 17, 981–989. doi: 10.1016/j.str.2009.06.002
Dias, E., Oliveira, M., Jones-Dias, D., Vasconcelos, V., Ferreira, E., Manageiro, V., et al. (2015). Assessing the antibiotic susceptibility of freshwater Cyanobacteria spp. Front. Microbiol. 6:799. doi: 10.3389/fmicb.2015.00799
El-Elimat, T., Zhang, X., Jarjoura, D., Moy, F. J., Orjala, J., Kinghorn, A. D., et al. (2012). Chemical diversity of metabolites from fungi, Cyanobacteria, and plants relative to FDA-approved anticancer agents. ACS Med. Chem. Lett. 3, 645–649. doi: 10.1021/ml300105s
Ferlay, J., Soerjomataram, I., Ervik, M., Dikshit, R., Eser, S., Mathers, C., et al. (2013). GLOBOCAN 2012 v1.0. Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. Lyon: International Agency for Research on Cancer. Available online at: http://globocan.iarc.fr
Fogliatto, G., Gianellini, L., Brasca, M. G., Casale, E., Ballinari, D., Ciomei, M., et al. (2013). NMS-E973, a novel synthetic inhibitor of Hsp90 with activity against multiple models of drug resistance to targeted agents, including intracranial metastases. Clin. Cancer Res. 19, 3520–3532. doi: 10.1158/1078-0432.CCR-12-3512
Gao, Y. D., and Huang, J. F. (2011). An extension strategy of Discovery Studio 2.0 for non-bonded interaction energy automatic calculation at the residue level (in Chinese). Zool. Res. 32, 262–266. doi: 10.3724/SP.J.1141.2011.03262
Gorman, M. A., Morera, S., Rothwell, D. G., de La Fortelle, E., Mol, C. D., Tainer, J. A., et al. (1997). The crystal structure of the human DNA repair endonuclease HAP1 suggests the recognition of extra-helical deoxyribose at DNA abasic sites. EMBO J. 16, 6548–6558. doi: 10.1093/emboj/16.21.6548
Gul, W., and Hamam, M. J. (2005). Indole alkaloid marine natural products: an established source of cancer drug leads with considerable promise for the control of parasitic, neurological and other diseases. Life Sci. 78, 442–453. doi: 10.1016/j.lfs.2005.09.007
Herman, M. D., Moche, M., Flodin, S., Welin, M., Tre'saugues, L., Johansson, I., et al. (2009). Structures of BIR domains from human NAIP and cIAP2. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 65, 1091–1096. doi: 10.1107/S1744309109038597
Koiri, R. K., Trigun, S. K., Mishra, L., Pandey, K., Dixit, D., and Dubey, S. K. (2009). Regression of Dalton's lymphoma in vivo via decline in lactate dehydrogenase and induction of apoptosis by a ruthenium(II)-complex containing 4-carboxy N-ethylbenzamide as ligand. Invest. New Drugs 27, 503–516. doi: 10.1007/s10637-008-9202-8
Krieger, E., Darden, T., Nabuurs, S. B., Finkelstein, A., and Vriend, G. (2004). Making optimal use of empirical energy functions: force-field parameterization in crystal space. Proteins 57, 678–683. doi: 10.1002/prot.20251
Kumar, A., and Singh, S. M. (1996). Effect of tumor growth on the blastogenic response of splenocytes: a role of macrophage derived nitric oxide. Immunol. Invest. 25, 413–423. doi: 10.3109/08820139609055731
Lee, H. U., Park, S. Y., Park, E. S., Son, B., Lee, S. C., Lee, J. W., et al. (2014). Photoluminescent carbon nanotags from harmful cyanobacteria for drug delivery and imaging in cancer cells. Sci. Rep. 4:4665. doi: 10.1038/srep04665
Liang, J., Moore, R. E., Moher, E. D., Munroe, J. E., Al-awar, R. S., Hay, D. A., et al. (2005). Cryptophycin-309249 and other cryptophycins analogs: preclinical efficacy studies with mouse and human tumors. Invest. New Drugs 23, 213–224. doi: 10.1007/s10637-005-6729-9
Madej, T., Lanczycki, C. J., Zhang, D., Thiessen, P. A., Geer, R. C., Marchler-Bauer, A., et al. (2014). MMDB and VAST+: tracking structural similarities between macromolecular complexes. Nucleic. Acids Res. 42, D297–D303. doi: 10.1093/nar/gkt1208
Mark, B. L., Mahuran, D. J., Cherney, M. M., Zhao, D., Knapp, S., and James, M. N. G. (2003). Crystal structure of human β-hexosaminidase B: understanding the molecular basis of Sandhoff and Tay–Sachs disease. J. Mol. Biol. 327, 1093–1109. doi: 10.1016/S0022-2836(03)00216-X
Mayer, A. M. S., and Hamann, M. J. (2005). Marine Pharmacology in 2001-2002, Marine compounds with anthelmintic, antibacterial, anticoagulant, antidiabetic, antifungal, anti-inflammatory, antimalarial, antiplatelet, antiprotozoal, antituberculosis and antiviral activities, affecting the cardiovascular, immune and nervous systems and other miscellaneous mechanisms of action. Comp. Biochem. Physiol. C Toxicol. Pharmacol. 140, 265–286. doi: 10.1016/j.cca.2005.04.004
Meerloo, J. V., Kaspers, G. J. L., and Cloos, J. (2011). “Cell sensitivity assays: The MTT assay,” in Cancer Cell Culture: Methods and Protocols, 2nd ed. Methods in Molecular Biology, ed I. A. Cree (New York, NY: Springer Science; Humana Press), 237–245. doi: 10.1007/978-1-61779-080-5_20
Morris, G. M., Goodsell, D. S., Halliday, R. S., Huey, R., Hart, W. E., Belew, R. K., et al. (1998). Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 19, 1639–1662. doi: 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-B
Oakley, A. J., Coggan, M., and Board, P. G. (2010). Identification and characterization of -glutamylamine cyclotransferase, an enzyme responsible for -glutamyl -lysine catabolism. J. Biol. Chem. 285, 9642–9648. doi: 10.1074/jbc.M109.082099
Pai, E. F., Krengel, U., Petsko, G. A., Goody, R. S., Kabsch, W., and Wittinghofer, A. (1990). Refined crystal structure of the triphosphate conformation of H-ras p21 at 1.35 A resolution: implications for the mechanism of GTP hydrolysis. EMBO J. 9, 2351–2359.
Patterson, G. M. L., Baldwin, C. L., Bolis, C. M., Caplan, F. R., Karuso, H., Larsen, L. K., et al. (1991). Antineoplastic activity of cultured blue-green algae (Cyanophyta). J. Phycol. 27, 530–536. doi: 10.1111/j.0022-3646.1991.00530.x
Prasad, R., and Koch, B. (2014). Antitumor Activity of ethanolic extract of Dendrobium formosum in T-Cell lymphoma: an in vitro and in vivo study. Biomed. Res. Int. 2014:753451. doi: 10.1155/2014/753451
Prasanna, R., Sood, A., Jaiswal, P., Nayak, S., Gupta, V., Chaudhary, V., et al. (2010). Rediscovering cyanobacteria as valuable sources of bioactive compounds. Appl. Biochem. Microbiol. 46, 119–134. doi: 10.1134/S0003683810020018
Shanab, S. M. M., Mostafa, S. S. M., Shalaby, E. A., and Mahmoud, G. I. (2012). Aqueous extracts of microalgae exhibit antioxidant and anticancer activities. Asian Pac. J. Trop. Biomed. 2, 608–615. doi: 10.1016/S2221-1691(12)60106-3
Sharma, G., and Koch, B. (2016). In vivo therapeutic efficacy of two ruthenium complexes and tumor regression in dalton's lymphoma ascites tumor model. Int. J. Pharmacol. Biol. Sci. 6, 121–127. doi: 10.21276/ijpbs.2016.6.3.14
Singh, K. B., and Trigun, S. K. (2013). Apoptosis of Dalton's lymphoma due to in vivo treatment with emodin is associated with modulations of hydrogen peroxide metabolizing antioxidant enzymes. Cell Biochem. Biophys. 67, 439–449. doi: 10.1007/s12013-011-9305-2
Sivonen, K., and Börner, T. (2008). “Bioactive compounds produced by cyanobacteria,” in The Cyanobacteria; Molecular Biology, Genomics and Evolution, ed A. Herrero and E. Flores (Norflok: Caster Academic Press), 159–197.
Srivastava, A., Tiwari, R., Srivastava, V., Singh, T. B., and Asthana, R. K. (2015). Fresh water cyanobacteria Geitlerinema sp. CCC728 and Arthrospira sp. CCC729 as an anticancer drug resource. PLoS ONE 10:e0136838. doi: 10.1371/journal.pone.0136838
Trott, O., and Olson, A. J. (2010). AutoDock VINA: improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J. Comput. Chem. 31, 455–461. doi: 10.1002/jcc.21334
Niveshika Verma, E., Mishra, A. K., Singh, A. K., and Singh, V. K. (2016). Structural elucidation and molecular docking of a novel antibiotic compound from cyanobacterium Nostoc sp. MGL001. Front. Microbiol. 7:1899. doi: 10.3389/fmicb.2016.01899
Yun, T. J., Harning, E. K., Giza, K., Rabah, D., Li, P., Arndt, J. W., et al. (2011). Ec144, a synthetic inhibitor of heat shock protein 90, blocks innate and adaptive immune responses in models of inflammation and autoimmunity. J. Immunol. 186, 563–575. doi: 10.4049/jimmunol.1000222
Zhang, Z., Li, Y., Lin, B., Schroeder, M., and Huang, B. (2011). Identification of cavities on protein surface using multiple computational approaches for drug binding site prediction. Bioinformatics 27, 2083–2088. doi: 10.1093/bioinformatics/btr331
Keywords: EMTAHDCA, anticancer drug potential, in silico, in vivo, in vitro analysis
Citation: Niveshika, Verma E, Maurya SK, Mishra R and Mishra AK (2017) The Combined Use of in Silico, in Vitro, and in Vivo Analyses to Assess Anti-cancerous Potential of a Bioactive Compound from Cyanobacterium Nostoc sp. MGL001. Front. Pharmacol. 8:873. doi: 10.3389/fphar.2017.00873
Received: 06 June 2017; Accepted: 13 November 2017;
Published: 27 November 2017.
Edited by:Amit K. Tiwari, University of Toledo, United States
Reviewed by:Amarjit Luniwal, North American Science Associates Inc., United States
Nikhil S. Khadtare, St. John's University, United States
Copyright © 2017 Niveshika, Verma, Maurya, Mishra and Mishra. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.