Artemisinin Derivatives Target Topoisomerase 1 and Cause DNA Damage in Silico and in Vitro

DNA topoisomerases 1 and 2 are enzymes that maintain DNA topology and play important essential genome functions, including DNA replication and transcription. Aberrant topoisomerases cause genome instability and a wide range of diseases, cancer in particular. Both Topo 1 and 2 are the targets of valuable anticancer drugs, such as camptothecin. It has been previously shown that artemisinin, a sesquiterpene lactone from Artemisia annua L. also known as qinghaosu, possesses anti-cancer effects and one of its derivatives, artesunate inhibits Topo 2. In this study, we evaluated artemisinin and 40 derivatives as potential Topo 1 inhibitors at first by in silico molecular docking analyses. Five compounds that showed comparable binding energies and similar docking poses were selected for in vitro cytotoxicity test and Comet assay for DNA damage. WWLL-013, WWLL-022, and WWLL-1098 showed the lowest binding energy also induced DNA damage in the Comet assay. CMK-0298 and CMK-0398 intercalated into DNA and induced mild DNA damage. All selected compounds, WWLL-013 in particular, were more cytotoxic toward the rat tumor cells than to the normal cells. In conclusion, the artemisinin derivatives such as CMK-0298, CMK-0398, WWLL-013, WWLL-022, and WWLL-1098 can be further developed as Topo 1 inhibitors.

Introduction

DNA topoisomerases are classified as type 1 and type 2 (Wang, 2002; Forterre et al., 2007). Type 1 topoisomerases (Topo 1) change the degree of supercoiling of DNA by single-strand breaks and religation, while type 2 topoisomerases (Topo 2) cause double-strand breaks. The enzymes relax supercoiling and torsional tension arisen during transcription and DNA replication, then permit DNA strand passage during early mitosis (Wang, 1985; Brill and Sternglanz, 1988; Kim and Wang, 1989). Six topoisomerase genes were discovered: two Topo 1 genes (TOP1 and TOP1mt), two Topo 2 genes (TOP2α and β), and two topoisomerases 3 (TOP3α and β) (Wang, 2002; Pommier, 2009).

Eukaryotic Topo 1 inserts transient breaks into duplex DNA. It breaks and rejoins one strand of duplex DNA resulting in the relaxation of both positively and negatively supercoiled DNA (Camilloni et al., 1989), whereas Topo 1 from E. coli can only relax negatively supercoiled DNA (Tse-Dinh et al., 1983). The relaxed DNA migrates as a series of bands of different sizes rather than a single band, known as topoisomers, which are DNAs of different linking number. This is a measure of a fundamental property of a closed-circular DNA molecule (Shimada and Yamakawa, 1984). Topo 1 catalyzes the equilibrium between these isomers. The energy difference between molecules of similar linking number around the average linking number of relaxed DNA is less than the thermal energy available at normal temperatures. Hence, the equilibrated sample consists of a mixture of topoisomers (Bates and Maxwell, 2005). Migration in an electric field through an agarose gel for circular molecules of the same molecular weight, is essentially a measure of writhe, which describes how the helix axis coils in space. Although adjacent bands differ in linking number, associated differences in writhe, which can cause their different mobility in the gel (Bates and Maxwell, 1997; Bates et al., 2013).

In general, whether a topoisomer of a given linking number is positively or negatively supercoiled or relaxed depends on the conditions, i.e., on the number of helical DNA repeats (Bates and Maxwell, 2005). An important factor influencing the twist and helical repeats of DNA is the presence of intercalating molecules. An intercalator such as ethidium bromide (Neidle and Abraham, 1984) contains a planar, usually polycyclic, aromatic structure, which inserts itself between two base pairs of double-stranded DNA. This causes local unwinding of the DNA helix, resulting in overall increase in the helical repeat that is a decrease of the DNA twist. A decreased twist of closed circular DNA results in increased writhe, causing increased DNA mobility. However, some classes of DNA-binding molecules such as netropsin have opposite effects on the DNA helix, i.e., they increase DNA twist.

Topoisomerases play important roles in DNA replication, transcription, mitotic chromosome formation, DNA recombination, and pre-mRNA splicing (Wang, 2002; Soret et al., 2003). Topo 1 is elevated in malignant tumors, including carcinomas of the colon, prostate, ovary and lung, (Husain et al., 1994; Karachaliou et al., 2013; Romer et al., 2013). High expression of Topo 1 and 2 correlates with poor response to treatment and survival of patients with small cell lung cancer (Karachaliou et al., 2013). Topo 1 expression is also closely associated with tumor metastasis and poor prognosis of patients with colorectal cancer (Romer et al., 2013; Silvestris et al., 2014).

Dysfunction of topoisomerases during cell division causes strand breaks, DNA damage and a collapse of the DNA replication complex, which ultimately leads to apoptosis of growing cells. Previously, Topo 1 stimulated research as the molecular target for cancer therapy owing to the discovery of the chemotherapeutic drug, camptothecin (Hsiang et al., 1985). This is a known Topo 1 inhibitor, and Topo 1 mutations can cause resistance to camptothecin (Chrencik et al., 2004). Therefore, the development of new compounds better targeting this enzyme has high priority.

Artemisinin, a sesquiterpene lactone from Artemisia annua L., reveals promising activity toward cancer cells in vitro and in vivo (Efferth, 2006; Lai et al., 2013; Takatani-Nakase, 2014; Tilaoui et al., 2014). Artesunate, a semisynthetic derivative of artemisinin was previously identified as Topo 2 inhibitor (Youns et al., 2009). In the present study, we investigated the activity of artemisinin derivatives to inhibit Topo 1 by molecular docking. Selected compounds were further evaluated by in vitro experiments.

Molecular docking can be applied as a virtual drug screening method based on 3D crystallographic structures of target proteins to predict the protein-ligand interactions. The bioinformatical calculations are based on parameters, which influence attraction and repulsion between target and compound, such as van der Waals interactions, hydrogen bonds, electrostatic interactions and hydrophobic interactions (Morris et al., 2008, 2009). The calculated binding energies and interacting amino acid residues characterize the affinity and the activity of a test compound toward its target. The lower the binding energy, the better the binding is of the compound to the target (Schneidman-Duhovny et al., 2004; Mihasan, 2012).

Materials and Methods

Cell Lines and Culture Conditions

All experiments were carried out on normal and transformed Rat 6 cells. The normal Rat 6 (R6) is an embryonal fibroblast cell line subcloned from an immortalized F2408 rat embryo cell line (Freeman et al., 1970). It has been a useful immortalized normal cell line in cancer research. Transformed Rat 6 (R6T24) cells are derived from the normal Rat 6 cell line transfected with the plasmid pT24, which carries the c-H-ras oncogene. Rat 6 and R6T24 cell lines were employed in this study to test the differential responses of normal and transformed cells to the test compounds. The culture medium used was Dulbecco's Modified Eagle's Medium (low glucose formulation) (DMEM; Gibco BRL, USA) supplemented with Basal Medium Eagle amino acids (BME; Sigma-Aldrich, Taufkirchen, Germany), 10% bovine calf serum (CS; Hyclone, USA) plus 5 mg/ml of each streptomycin and penicillin. Cell cultures were maintained in a humidified incubator at 37°C with 5% CO₂. Cells were subcultured by trypsinizing the cells with PBS containing 0.05% trypsin and 1 mM EDTA and replated at the appropriate seeding cell density per tissue culture (TC) plate (Corning, USA). Culture medium was changed twice a week.

Chemicals and Reagents

Agarose (normal melting and low melting) was purchased from FMC (Houston, USA). Boric acid was purchased from BDH Chemicals (Dawsonville, GA, USA). Dimethyl sulfoxide, ethylenediaminetetraacetic acid-disodium salt, ethidium bromide, trizma base, triton-X-100, cis-platinum, propidium iodide, ribonuclease A, and paraformaldehyde were purchased from Sigma. Ethanol, hydrochloric acid, sodium hydroxide are purchased from Honeywell-Riedel-de Haën (Seelze, Germany). Sodium chloride, topoisomerase and trypsin were purchased from Gibco BRL (USA). Sodium sarcosinate was purchased from Fluka Biochemika (Germany).

Molecular Docking

The protocol for molecular docking was previously reported by us (Zeino et al., 2014). An X-ray crystallography-based 3D structure of wild-type human Topo 1 (PDB ID: 1K4T) and its N722S mutant (PDB ID: 1RRJ) were obtained from Protein Data Bank (http://www.rcsb.org/pdb). This mutant confers camptothecin resistance (Chrencik et al., 2004). Both structures possess a phosphorylated tyrosine at position 723. The selected residues for defined docking were obtained from the docking pose of topotecan (a camptothecin derivative) on human Topo 1 and its interacting residues (Staker et al., 2002), i.e., DA5, DG6, DA7, DC8, DT9, DC117, DT118, DT119, Arg364, Lys532, Asp533, Asn722, Tyr723.

Chemical structures of artemisinin derivatives were retrieved from Pubchem database and the 3D structures were created by Corina Online. Chemical properties and classification of the derivatives were described in the Table S1. The known Topo 1 inhibitors; camptothecin (Baikar and Malpathak, 2010) and topotecan (Kollmannsberger et al., 1999) as well as hydrolyzed topotecan were used as control drugs. A grid box was then constructed to define docking spaces according to the selected residues for defined docking, whereas the whole protein surface was covered for blind docking. Docking parameters were set to 250 runs and 2,500,000 energy evaluations for defined docking, 25,000,000 for blind docking. Docking was performed three times independently by Autodock4 implemented on AutodockTools-1.5.7rc1 (Morris et al., 2009) using the Lamarckian Genetic Algorithm. Corresponding lowest binding energies and predicted inhibition constants were represented as mean ± SD. Visual Molecular Dynamics (VMD) was used to depict the docking poses of artemisinin derivatives and control drugs.

DNA Decatenation Analysis

DNA topoisomerases are enzymes that unwind supercoiled DNA. Topo 1 was used to determine, whether the test compounds intercalate into DNA. The test compounds were incubated with closed circular pUC18 DNA in the presence of Topo 1 and electrophoresis was performed in an agarose gel. A competition in the direction of the writhe may result from the simultaneous reaction of the DNA with both Topo 1 and test compound. Two opposite forces appear on DNA during this reaction. Topo 1 relaxes the closed circular DNA decreasing the writhe to zero. On the other hand, if the test compound intercalates into DNA, the writhe increases due to decreased DNA twist. Depending on the strength of the pull from each side, DNA has different conformational states and its mobility shifts.

The plasmid pUC18 (2.69 Kb) was extracted from DH5α, an E. coli bacterial strain using Qiagen Spin column and was used in all experiments. The reaction mix contained 50 ng of pUC18 plasmid DNA (1 μl), 2 μl of 5 × reaction buffer (50 mM Tris HCl pH 7.5, 50 mM KCl, 10 mM MgCl₂, 0.5 mM DTT, 0.1 mM EDTA and 30 μg/ml BSA), two units of Topo 1 enzyme, 1 μl of the test chemical and dH₂O to make up a final volume of 10 μl. The reaction mix was incubated at 37°C water bath for 2 h and 2 μl of 6 × loading dye was added. The entire volume of 12 μl was then loaded onto a 1% agarose gel without ethidium bromide (EtBr) and electrophoresed in TAE buffer (40 mM Tris- acetate, 1 mM EDTA). The gel was then stained with 50 μg/ml of EtBr for 20 min and viewed with Digital Imaging System (BIO-RAD, USA).

Comet Assay

Cells were treated as described below under “Cytotoxicity Test.” After 6 days treatment, cells were counted and resuspended in an appropriate volume of PBS to reach a final density of ~3,000–10,000 cells/μl. Cell suspensions were kept on ice prior to the comet assay. We applied the same treatment duration to the cell cultures as the duration used for MTT assay.

Fisher's fully frosted microscope slides (Fisher Scientific, UK) were each covered with 100 μl of 0.5% normal melting agarose (NMA) in PBS at 45°C, immediately covered with a coverslip and then kept at 4°C for 10 min to allow the agarose to solidify. Treated or untreated control cells suspended in 10 μl PBS were mixed with 75 μl 0.5% low melting point agarose (LMA) at 37°C. After gently removing the coverslip, the cell suspension was rapidly pipetted onto the first layer of agarose, spread using a coverslip and maintained at 4°C for 10 min. After removal of the coverslip, a third layer of 0.5% LMA (80 μl) at 37°C was added and spread using a coverslip and allowed to solidify at 4°C (this layer was optional). After the removal of the coverslip, the slides were immersed in freshly prepared cold lysing solution (2.5 M NaCl, 100 mM Na₂EDTA, 10 mM Tris; pH10, 1% sodium sarcosinate with 1% Triton X-100 and 10% DMSO) for at least 1 h at 4°C. The slides were then removed from the lysing solution, drained and placed in a horizontal gel electrophoresis tank side by side. The tank was filled with fresh electrophoresis buffer (1 mM Na₂EDTA and 300 mM NaOH) to a level slightly above the slides. The slides were left in this buffer for 20 min to allow unwinding of the DNA before electrophoresis. Electrophoresis was conducted at room temperature for 20 min at 25 V, adjusted to 300 mA by raising or lowering the buffer level in the tank. After electrophoresis, the slides were washed gently to remove alkali and detergents which would interfere with ethidium bromide staining, then placing on a tray and flooding them slowly with three changes of 0.4 M Tris pH7.5, each for 5 min. After neutralization, the slides were stained with 50–100 μl of 2 μg/ml ethidium bromide in distilled water, covered with a coverslip, and stored in the dark. Observations were made using Zeiss Axiophot fluorescence microscope (Zeiss, Jena, Germany). Photographs of single cells were taken at 400 × magnification using P1600 Ektachrome film (Kodak, Stuttgart, Germany).

MTT Cytotoxicity Assay

The lead compound artemisinin and its derivatives were initially tested for their cytotoxicity effects on Rat 6 and R6T24 cell lines. The toxicity effects of the chemicals were assessed using MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) assay at the end of the treatment periods. For the assay, 800 cells were seeded in each well of a 96-well plate in a final volume of 100 μl of DMEM medium supplemented with 10% calf serum. Six duplicate wells were used for each sample. On the next day, cells were fed with drug-containing medium and kept incubated for another 3 or 6 days. Drugs were freshly added to the cultures when the cells were fed every other day with fresh medium. The final concentration of each compound ranged from 0 to 100 μM. DMSO or tetrahydrofuran was used as solvent. Depending on the solubility of the compounds, 1,000 × stock (i.e., 100 mM) of the compounds was prepared. Due to the limitation of the solubility of the test compounds, the maximum test concentrations varied from 25 to 100 μM as designated in the experiment. At the end of the drug treatment, 50 μl of the MTT labeling reagent was added at a final concentration of 0.5 mg/ml to each well. The plate was re-incubated for 4 h in a humidified atmosphere prior to the addition of 100 μl of the solubilizing solution (10% SDS in 0.01 M HCl) into each well. The plate was allowed to stand overnight in the incubator to enhance the complete solubilization of the purple formazan crystals. The plate was scanned for measuring the absorbance using a Dynatech MR5000 (Dynatach, USA) at wavelengths between 550 and 600 nm. The data were expressed as the mean ± SEM (n = 6) and were analyzed by the Student's t-test with a p < 0.05 considered to be significant.

Real-Time PCR-Based Melting DNA Shifting Assay in the Presence of Cisplatin and Artemisinin

RNA was isolated from HeLa cells using Qiagen RNasy mini kit (Qiagen). Reverse transcription was performed in 20 μl reaction mixture [Transcriptor Universal cDNA Master (Roche)] with 2.5 μg total RNA. The synthesized cDNA was diluted 1:10 before used as the template. β-Actin primers (5″- CTC TTC CAG CCT TCC TTC CT-3″ –F & 5″ -AGC ACT GTG TTG GCG TAC AG- 3″ -R) were used in real-time PCR (qPCR) assay. The 20 μl qPCR reaction mixture contained 2 μl template; 0.3 μl each of 1 mM forward and reverse primers; 10 μl PowerUp SYBR™ Green Master Mix (Thermo Fisher); 2 μl 0.1 mM compound (artemisinin or cisplatin) or DMSO; 5.4 μl ddH2O. The qPCR was conducted as following: an initial 2 min at 50°C, 10 min at 95°C; followed by 40 cycles of 15 s at 95°C and 1 min at 60°C. The last melting curve stage was conducted at 15 s at 95°C, 1 min at 60°C, 15 s at 95°C. The dissociation curve was analyzed using ViiA 7 RUO software (Life technology, Carlsbad, USA) to demonstrate the shift of melting curve due to the binding of the test compound to DNA molecules.

Results

Evaluation of 40 Artemisinin Derivatives as Potential Topo 1 Inhibitors using Molecular Docking Analysis

Artemisinin derivatives were firstly evaluated by in silico molecular docking on human Topo 1 to identify potential inhibitors. Two protein structures were used: wildtype and N722S mutant Topo 1. Docking poses and binding energies were compared with known inhibitors, i.e., topotecan, hydrolyzed topotecan and camptothecin. The docking results were summarized in Table 1 for docking to the wildtype protein and in Table 2 to the mutant protein. Artemisinin derivatives that showed comparable binding energies and similar docking poses with the control inhibitors were selected for further analyses (they were labeled bold in Tables 1, 2). Chemical structures of the selected artemisinin derivatives (CMK0298, CMK0398, WWLL-013, WWLL-022, and WWLL-1098) are depicted in Figure 1. These selected compounds possessed similar docking poses as the control inhibitors both on wildtype (Figure 2) and mutant Topo 1 (Figure 3).

TABLE 1

Table 1. Molecular docking analyses of 40 artemisinin derivatives on wild-type human Topo 1 (PDB ID:1K4T).

TABLE 2

Table 2. Molecular docking analyses of 40 artemisinin derivatives on N722S mutant human Topo 1 (PDB ID:1RRJ).

FIGURE 1

Figure 1. Chemical structures of selected artemisinin derivatives (CMK0298, CMK0398, WWLL-013, WWLL-022, and WWLL-1098).

FIGURE 2

Figure 2. CMK0298 (cyan), CMK0398 (black), WWLL-1098 (violet), WWLL-022 (yellow), WWLL-013 (blue), topotecan (green), hydrolyzed topotecan (red), camptothecin (silver) docking on human Topo 1 (magenta. PDB ID:1K4T).

FIGURE 3

Figure 3. CMK0298 (cyan), CMK0398 (black), WWLL-1098 (violet), WWLL-022 (yellow), WWLL-013 (blue), topotecan (green), hydrolyzed topotecan (red), camptothecin (silver) docking on N722S mutant human Topo1 (orange. PDB ID 1RRJ).

The Artemisinin Derivatives Interfered with Topo 1 Activity as Demonstrated by DNA Decatenation Assay

To test whether the selected artemisinin derivatives interfere Topo 1 activities, CMK0298 and CMK0398 were incubated with supercoiled pUC18 DNA in the presence of Topo 1 and subjected to agarose electrophoresis. Bleomycin, a cancer therapeutic agent causing DNA strand breaks, was used as positive control for the assay. The force to increase the writhe by intercalating into DNA overrides the force of Topo 1 activity, if the concentration was above 1 μg/ml. Therefore, the DNA isomers were converted to a more compact conformation causing a faster migration in the gel. At saturating amounts of bleomycin (>2 μg/ml), the DNA migrated at similar rates as untreated closed circular DNA (refer to lane 2). At around the equilibrium concentration (point of balance between the opposing forces) or below (0.75 μg/ml), the DNA band pattern was similar to that of treatment only with Topo 1 (refer to lane 3). This indicated that Topo 1 activity was superior to DNA intercalation (if any) in the reaction. The conformation of the topoisomers at these low concentrations appeared not to be affected by treatment of with test compounds. Similar findings were observed with bleomycin. Closed circular pUC18 DNA that was incubated with Topo 1 together with the compounds, CMK-0298 and CMK-0398, which showed similar results as bleomycin (Figure 4). However, a noticeable difference was observed for the equilibrium concentrations of these test compounds, which were higher than those of the classical intercalators. One possible explanation may be the smaller size of these test compounds. Hence, more molecules were required to achieve equivalent unwinding effects compared to those of classical DNA intercalators, such as bleomycin.

FIGURE 4

Figure 4. DNA decatenation assay for CMK-0298 and CMK-0398. Bleomycin was used as control compounds. Lane 1, 1 kb DNA marker; Lane 2, supercoiled pUC18 plasmid DNA; Lane 3, Topo 1-treated pUC18.The other lanes show pUC18 DNA treated with Topo 1 plus test compounds at dosages indicated below: CMK-0298: lanes 4–8: 1, 10, 100, 200, 500 μM; CMK-0398: lanes 4–7: 0.01, 0.1, 1, 2.5 mM; Bleomycin: lanes 4–7: 0.1 mg/mL. 10 ng/mL, 1 ng/mL, 0.1 ng/mL; Ethidium bromide: lanes 4–9: 2, 1.5, 1, 0.75, 0.5, 0.2 μg/ml.

Artemisinin Derivatives Showed Various Degree of DNA Damage as Tested by Comet Assay

The dysfunction of Topo I induced by inhibitors is often associated with DNA damage. We then tested the effect of the selected artemisinin derivatives on DNA damage based on the formation of comet tail upon drug treatment. The data showed that CMK-0298, CMK-0398, WWLL-013, WWLL-022, and WWLL-1098 all showed various degree of positivity on Comet assay (Figure 4). The induced dosages ranged from 0.5 to 1 μM, which were in general higher than the IC₅₀ values of the tested derivative (Figure 5, Table 3). Bleomycin caused DNA damage at relative concentrations (Figure 6). Since CMK-0298 and CMK-0398 also showed positive results in the DNA decatenation assay, these compounds might indeed induce DNA damage.

FIGURE 5

Figure 5. MTT assays of artemisinin derivatives on the normal R6 and transformed R6 (R6T24) cells treated with CMK-0298 (A), CMK-0398 (B), WWLL-013 (C), WWLL-022 (D) and WWLL-1098 (E) respectively.

TABLE 3

Table 3. IC₅₀ values of selected artemisinin derivatives on the normal R6 and transformed R6 (R6T24) cells assessed by MTT assay.

FIGURE 6

Figure 6. Comet assays of R6 cells treated with CMK-0298, CMK-0398, WWLL-013, and WWLL-022. Bleomycin was used as the control compound.

As can be seen in Table 3, CMK-0298, CMK-0398, WWLL-013, WWLL-022, and WWLL-1098 were cytotoxic toward both cell lines in the low micromolar range. It is of interest to know that WWLL-013 showed a big discrepancy in IC₅₀ between normal and the transformed R6 cells. This implies that this compound possess differential cytotoxicity against transformed cells, suggesting that there might be a leverage to use this compound to treat cancer with less harmful effect to the normal cells.

Artemisinin does not Directly Intercalate into DNA to Interfere with Topo 1 Activity

In order to discriminate between DNA intercalation and Topo 1 inhibition, we conducted real time PCR assays for the house-keeping gene, actin and monitored its DNA dissociation curve and melting temperature in the presence of artemisinin or the known DNA intercalator and adduct-forming drug, cisplatin (Sangeetha Gowda et al., 2014). The melting temperature shifted from 80.05 to 84.42°C in the presence of cisplatin as determined by its DNA dissociation curve (Figure 7). On the other hand, the DNA dissociation curve in the presence of artemisinin was closely overlapping with that of the DMSO solvent control, suggesting that artemisinin did not exert DNA binding activity. This finding provides evidence that artemisinin acts directly on Topo 1 as inhibitor.

FIGURE 7

Figure 7. DNA melting curves of PCR products incubated in DMSO or test chemicals. The qPCR was conducted according to the “Materials and Methods.” DMSO, 10 μM cisplatin, or 10 μM artemisinin were added to the qPCR reaction.

Discussion

Molecular docking studies pointed to artemisinin derivatives as Topo 1 inhibitors. Interestingly, the selected artemisinin derivatives (CMK-0298, CMK-0398, WWLL 1098, WWLL-022, and WWLL-013) might be indeed Topo 1 inhibitors, since comparable binding energies and similar docking poses were observed with the known inhibitors. Further in vitro analyses supported their potential to be used as Topo 1 inhibitors.

WWLL-013 (−12.42 kcal/mol) and WWLL-022 (−12.21 kcal/mol) showed stronger binding than the known inhibitors on the wild-type Topo 1, whereas WWLL-1098 (−10.96 kcal/mol) and CMK-0398 (−10.72 kcal/mol) showed stronger binding than hydrolyzed topotecan (−10.53 kcal/mol) and camptothecin (−10.29 kcal/mol). CMK-0298 showed comparable binding energies (−9.97 kcal/mol) with the inhibitors. The selected compounds shared common interacting residues with the known inhibitors (Kollmannsberger et al., 1999; Baikar and Malpathak, 2010). WWLL-013, WWLL-022, and WWLL-1098 revealed six common residues, whereas CMK-0298 and CMK-0398 revealed two common residues. Molecular docking analysis on the N722S mutant Topo 1 revealed that camptothecin (−9.94 kcal/mol) docked to a different binding site, which is in corroboration with the hypothesis that the indicated mutation causes camptothecin resistance. The missing −NH₂ group at serine 722 might be important for its binding. WWLL-013 (−12.58 kcal/mol) and WWLL-022 (−11.72 kcal/mol) showed comparable binding energies with topotecan (−12.66 kcal/mol) and four common interacting residues were observed with topotecan. CMK-0398 (−10.63 kcal/mol) and WWLL-1098 (–10.56 kcal/mol) showed stronger binding than hydrolyzed topotecan (−9.95 kcal/mol) and camptothecin (−9.94 kcal/mol), CMK-0298 showed comparable binding energy (–9.37 kcal/mol). The selected artemisinin derivatives showed strong interaction on wild-type and N722S mutant Topo 1 and shared common interacting residues, implying that they might be used as potent Topo 1 inhibitors. In order to prove the in silico findings, we performed experiments in vitro.

DNA decatenation assay (Sahai and Kaplan, 1986; Boos and Stopper, 2001) is a well-established method to detect DNA topoisomers, if an intercalator molecule is present. Under normal conditions, Topo 1 relaxes the closed circular DNA, and, hence, decreases the writhe to zero. In presence of an intercalator compound, the writhe will increase due to the decrease in the DNA twist leading to conformational changes and, thus, different bands at agarose gel. CMK-0298 and CMK-0398 indeed showed intercalator properties similar to bleomycin and ethidium bromide, implying that they possess DNA intercalator ability and thus Topo 1 inhibitory effect. Further evidence was provided by comet assay in order to evaluate DNA damage by the selected compounds. CMK-0298, CMK-0398, WWLL-013, and WLL-022 indeed induced DNA break. The fact that CMK-0298 and CMK-0398 also showed activity in decatenation assay, implies that they might be promising topoisomerase 1 inhibitors, as they induce DNA damage and possess intercalation activity.

Selected artemisinin derivatives showed cytotoxicity toward transformed rat fibroblast cells transfected with c-H-Ras oncogene in a low micromolar range. Moreover, the transformed cells were more sensitive than the normal cells to the all test derivatives, WWLL-013 in particular. This points out that they might possess anti-cancer activity as they are able to kill tumor model cells. Based on the DNA disassociation curve experiment (Figure 7), we showed that artemisinin does not directly intercalate into DNA to interfere with Topo 1 activity.

It has been recently reported that artemisinin and its derivatives such as dihydroartemisinin, artesunate, artemether, arteether possess a multi-target effect toward cancer as it influences various signaling pathways playing important roles in cancer progression (Li et al., 2016; Efferth, 2017). In the present study we identified artemisinin derivatives that could target Topo 1.

In conclusion, CMK-0298, CMK-0398, WWLL-013, WWLL-022, and WWLL-1098 possess anti-cancer activity and potential Topo 1 inhibitor properties as shown in silico and in vitro. Further preclinical and clinical studies are warranted in order to evaluate their therapeutic potential for clinical use.

Author Contributions

TE and WH conceived the study. OK performed the in silico experiments. VW, HF contributed to the in silico experiments. OK and TE wrote the manuscript. AC, ACLQ, VW performed the in vitro experiments. All the authors read the manuscript.

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.

Acknowledgments

We thank Prof. Richard K. Haynes, Centre of Excellence for Pharmaceutical Sciences, Faculty of Health Sciences G2, North-West University Potchefstroom Campus, Potchefstroom, South Africa) for the provision of the artemisinin derivatives. This study was partially supported by Macau FDCT Grants 015/2014/A1 and 103/2016/A3 to WH.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2017.00711/full#supplementary-material

References

Baikar, S., and Malpathak, N. (2010). Secondary metabolites as DNA topoisomerase inhibitors: a new era towards designing of anticancer drugs. Pharmacogn. Rev. 4, 12–26. doi: 10.4103/0973-7847.65320

PubMed Abstract | CrossRef Full Text | Google Scholar

Bates, A. D., and Maxwell, A. (1997). DNA topology: Topoisomerases keep it simple. Curr. Biol. 7, R778–R781. doi: 10.1016/S0960-9822(06)00403-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Bates, A. D., and Maxwell, A. (2005). DNA Topology. New York, NY: Oxford University Press.

Google Scholar

Bates, A. D., Noy, A., Piperakis, M. M., Harris, S. A., and Maxwell, A. (2013). Small DNA circles as probes of DNA topology. Biochem. Soc. Trans. 41, 565–570. doi: 10.1042/BST20120320

PubMed Abstract | CrossRef Full Text | Google Scholar

Boos, G., and Stopper, H. (2001). DNA methylation influences the decatenation activity of topoisomerase II. Int. J. Biol. Macromol. 28, 103–106. doi: 10.1016/S0141-8130(00)00157-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Brill, S. J., and Sternglanz, R. (1988). Transcription-dependent DNA supercoiling in yeast DNA topoisomerase mutants. Cell 54, 403–411. doi: 10.1016/0092-8674(88)90203-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Camilloni, G., Di Martino, E., Di Mauro, E., and Caserta, M. (1989). Regulation of the function of eukaryotic DNA topoisomerase I: topological conditions for inactivity. Proc. Natl. Acad. Sci. U.S.A. 86, 3080–3084. doi: 10.1073/pnas.86.9.3080

PubMed Abstract | CrossRef Full Text | Google Scholar

Chrencik, J. E., Staker, B. L., Burgin, A. B., Pourquier, P., Pommier, Y., Stewart, L., et al. (2004). Mechanisms of camptothecin resistance by human topoisomerase I mutations. J. Mol. Biol. 339, 773–784. doi: 10.1016/j.jmb.2004.03.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Efferth, T. (2006). Molecular pharmacology and pharmacogenomics of artemisinin and its derivatives in cancer cells. Curr. Drug Targets 7, 407–421. doi: 10.2174/138945006776359412

PubMed Abstract | CrossRef Full Text | Google Scholar

Efferth, T. (2017). From ancient herb to modern drug: Artemisia annua and artemisinin for cancer therapy. Semin. Cancer Biol. 46, 65–83. doi: 10.1016/j.semcancer.2017.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Forterre, P., Gribaldo, S., Gadelle, D., and Serre, M. C. (2007). Origin and evolution of DNA topoisomerases. Biochimie 89, 427–446. doi: 10.1016/j.biochi.2006.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Freeman, A. E., Price, P. J., Igel, H. J., Young, J. C., Maryak, J. M., and Huebner, R. J. (1970). Morphological transformation of rat embryo cells induced by diethylnitrosamine and murine leukemia viruses. J. Natl. Cancer Instit. 44, 65–78.

PubMed Abstract | Google Scholar

Hsiang, Y. H., Hertzberg, R., Hecht, S., and Liu, L. F. (1985). Camptothecin induces protein-linked DNA breaks via mammalian DNA topoisomerase I. J. Biol. Chem. 260, 14873–14878.

PubMed Abstract | Google Scholar

Husain, I., Mohler, J. L., Seigler, H. F., and Besterman, J. M. (1994). Elevation of topoisomerase I messenger RNA, protein, and catalytic activity in human tumors: demonstration of tumor-type specificity and implications for cancer chemotherapy. Cancer Res. 54, 539–546.

PubMed Abstract | Google Scholar

Karachaliou, N., Papadaki, C., Lagoudaki, E., Trypaki, M., Sfakianaki, M., Koutsopoulos, A., et al. (2013). Predictive value of BRCA1, ERCC1, ATP7B, PKM2, TOPOI, TOPOmicron-IIA, TOPOIIB and C-MYC genes in patients with small cell lung cancer (SCLC) who received first line therapy with cisplatin and etoposide. PLoS ONE 8:e74611. doi: 10.1371/journal.pone.0074611

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, R. A., and Wang, J. C. (1989). Function of DNA topoisomerases as replication swivels in Saccharomyces-cerevisiae. J. Mol. Biol. 208, 257–267. doi: 10.1016/0022-2836(89)90387-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kollmannsberger, C., Mross, K., Jakob, A., Kanz, L., and Bokemeyer, C. (1999). Topotecan - A novel topoisomerase I inhibitor: pharmacology and clinical experience. Oncology 56, 1–12. doi: 10.1159/000011923

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai, H. C., Singh, N. P., and Sasaki, T. (2013). Development of artemisinin compounds for cancer treatment. Invest. New Drugs 31, 230–246. doi: 10.1007/s10637-012-9873-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Li, Q., Wu, J., Wang, M. Y., and Yu, J. X. (2016). Artemisinin and its derivatives as a repurposing anticancer agent: what else do we need to do? Molecules 7, 21. doi: 10.3390/molecules21101331

CrossRef Full Text | Google Scholar

Mihasan, M. (2012). What in silico molecular docking can do for the bench-working biologists. J. Biosci. 37, 1089–1095. doi: 10.1007/s12038-012-9273-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, G. M., Huey, R., and Olson, A. J. (2008). Using Auto Dock for ligand-receptor docking. Curr. Protoc. Bioinform. 8, 14. doi: 10.1002/0471250953.bi0814s24

CrossRef Full Text | Google Scholar

Morris, G. M., Huey, R., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). AutoDock4 and Auto Dock tools4: automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791. doi: 10.1002/jcc.21256

CrossRef Full Text | Google Scholar

Neidle, S., and Abraham, Z. (1984). Structural and sequence-dependent aspects of drug intercalation into nucleic acids. CRC Crit. Rev. Biochem. 17, 73–121. doi: 10.3109/10409238409110270

PubMed Abstract | CrossRef Full Text | Google Scholar

Pommier, Y. (2009). DNA topoisomerase I inhibitors: chemistry, biology, and interfacial inhibition. Chem. Rev. 109, 2894–2902. doi: 10.1021/cr900097c

PubMed Abstract | CrossRef Full Text | Google Scholar

Romer, M. U., Nygard, S. B., Christensen, I. J., Nielsen, S. L., Nielsen, K. V., Muller, S., et al. (2013). Topoisomerase 1(TOP1) gene copy number in stage III colorectal cancer patients and its relation to prognosis. Mol. Oncol. 7, 101–111. doi: 10.1016/j.molonc.2012.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahai, B. M., and Kaplan, J. G. (1986). A quantitative decatenation assay for type II topoisomerases. Anal. Biochem. 156, 364–379. doi: 10.1016/0003-2697(86)90267-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneidman-Duhovny, D., Nussinov, R., and Wolfson, H. J. (2004). Predicting molecular interactions in silico: II. Protein-protein and protein-drug docking. Curr. Med. Chem. 11, 91–107. doi: 10.2174/0929867043456223

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimada, J., and Yamakawa, H. (1984). DNA-topoisomer analysis on the basis of the helical wormlike chain. Biopolymers 23, 853–857. doi: 10.1002/bip.360230503

PubMed Abstract | CrossRef Full Text | Google Scholar

Silvestris, N., Simone, G., Partipilo, G., Scarpi, E., Lorusso, V., Brunetti, A. E., et al. (2014). CES2, ABCG2, TS and Topo-I primary and synchronous metastasis expression and clinical outcome in metastatic colorectal cancer patients treated with first-line FOLFIRI regimen. Int. J. Mol. Sci. 15, 15767–15777. doi: 10.3390/ijms150915767

PubMed Abstract | CrossRef Full Text | Google Scholar

Sangeetha Gowda, K. R., Mathew, B. B., Sudhamani, C. N., and Naik, H. S. B. (2014). Mechanism of DNA binding and cleavage. Biomed. Biotechnol. 2, 1–9. doi: 10.12691/bb-2-1-1

CrossRef Full Text | Google Scholar

Soret, J., Gabut, M., Dupon, C., Kohlhagen, G., Stevenin, J., Pommier, Y., et al. (2003). Altered serine/arginine-rich protein phosphorylation and exonic enhancer-dependent splicing in Mammalian cells lacking topoisomerase I. Cancer Res. 63, 8203–8211.

PubMed Abstract | Google Scholar

Staker, B. L., Hjerrild, K., Feese, M. D., Behnke, C. A., Burgin, A. B. Jr., and Stewart, L. (2002). The mechanism of topoisomerase I poisoning by a camptothecin analog. Proc. Natl. Acad. Sci. U.S.A. 99, 15387–15392. doi: 10.1073/pnas.242259599

PubMed Abstract | CrossRef Full Text | Google Scholar

Takatani-Nakase, T. (2014). [Artemisinin: a natural product for fighting against cancer]. Nippon. Yakurigaku Zasshi 143, 61–64. doi: 10.1254/fpj.143.61

PubMed Abstract | CrossRef Full Text | Google Scholar

Tilaoui, M., Mouse, H. A., Jaafari, A., and Zyad, A. (2014). Differential effect of artemisinin against cancer cell lines. Nat. Prod. Bioprospect. 4, 189–196. doi: 10.1007/s13659-014-0024-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Tse-Dinh, Y. C., Mccarron, B. G., Arentzen, R., and Chowdhry, V. (1983). Mechanistic study of E. coli DNA topoisomerase I: cleavage of oligonucleotides. Nucleic Acids Res. 11, 8691–8701. doi: 10.1093/nar/11.24.8691

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. C. (1985). DNA topoisomerases. Annu. Rev. Biochem. 54, 665–697. doi: 10.1146/annurev.bi.54.070185.003313

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. C. (2002). Cellular roles of DNA topoisomerases: a molecular perspective. Nat. Rev. Mol. Cell Biol. 3, 430–440. doi: 10.1038/nrm831

PubMed Abstract | CrossRef Full Text | Google Scholar

Youns, M., Efferth, T., Reichling, J., Fellenberg, K., Bauer, A., and Hoheisel, J. D. (2009). Gene expression profiling identifies novel key players involved in the cytotoxic effect of Artesunate on pancreatic cancer cells. Biochem. Pharmacol. 78, 273–283. doi: 10.1016/j.bcp.2009.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeino, M., Saeed, M. E., Kadioglu, O., and Efferth, T. (2014). The ability of molecular docking to unravel the controversy and challenges related to P-glycoprotein–a well-known, yet poorly understood drug transporter. Invest. New Drugs 32, 618–625. doi: 10.1007/s10637-014-0098-1

PubMed Abstract | CrossRef Full Text | Google Scholar