Investigations on Anticancer Potentials by DNA Binding and Cytotoxicity Studies for Newly Synthesized and Characterized Imidazolidine and Thiazolidine-Based Isatin Derivatives

Imidazolidine and thiazolidine-based isatin derivatives (IST-01–04) were synthesized, characterized, and tested for their interactions with ds-DNA. Theoretical and experimental findings showed good compatibility and indicated compound–DNA binding by mixed mode of interactions. The evaluated binding parameters, i.e., binding constant (Kb), free energy change (ΔG), and binding site sizes (n), inferred comparatively greater and more spontaneous binding interactions of IST-02 and then IST-04 with the DNA, among all compounds tested under physiological pH and temperature (7.4, 37 °C). The cytotoxic activity of all compounds was assessed against HeLa (cervical carcinoma), MCF-7 (breast carcinoma), and HuH-7 (liver carcinoma), as well as normal HEK-293 (human embryonic kidney) cell lines. Among all compounds, IST-02 and 04 were found to be cytotoxic against HuH-7 cell lines with percentage cell toxicity of 75% and 66%, respectively, at 500 ng/µL dosage. Moreover, HEK-293 cells exhibit tolerance to the increasing drug concentration, suggesting these two compounds are less cytotoxic against normal cell lines compared to cancer cell lines. Hence, both DNA binding and cytotoxicity studies proved imidazolidine (IST-02) and thiazolidine (IST-04)-based isatin derivatives as potent anticancer drug candidates among which imidazolidine (IST-02) is comparatively the more promising.


Chemistry of Synthesized Compounds
Four new hydrazonoindolin-2-one compounds (IST-01-04) were synthesized with a general formula of [3-(hydrazono)1-benzyl-5-chloro-indolin-2-one]. Three different strategies were adopted to produce such derivatives ( Figure 3). In the first strategy ( Figure 3A), we focused on big lipophilic groups (coumarins and benzocoumarins) in the thiazole ring in the lead structure to achieve the target compounds IST-01 and IST-03. In the second strategy ( Figure 3B), a hydrophilic moiety was utilized to replace the lipophilic groups to obtain the target structure IST-04. In the third strategy ( Figure 3C), a bio-isosteric tactic was utilized to replace the thiazole ring by 2-thioxo-1,3-imidazolidine-4-one to obtain the target compound IST-02.
The final products were obtained by performing the reactions as shown in Scheme 1. The reagent 5-chloroisatin was treated with benzyl bromide to obtain product (3), which was further reacted with thiosemicarbazide to obtain the Schiff base product. In the last step, the Schiff base product (5) was reacted with alpha-halogenated compounds to yield desired compounds 1-benzyl-5-chloro-3-(2-heteroarylhydrazono) indoline-2-ones (IST-01-04). FT-IR analysis revealed that all synthesized compounds expressed characteristic absorption band for N-H at 3233 cm −1 , C-H aromatic at 3119 cm −1 , C=O of carbonyl at 1738 cm −1 , and C=N of imine at 1549 cm −1 . The 1 H-NMR spectra of IST-01 showed NH proton at 13.34 ppm, coumarin proton at 8.63 ppm, thiazoline proton at 8.10 ppm. For IST-02, the 1 H-NMR spectra showed isatin proton at 8.03 ppm, phenyl proton as a multiplet at 7.34-7.28 ppm, thiohydantoin proton at 4.97 ppm. For IST-03, 1 H-NMR spectra showed N-H proton at 12.88 ppm, thiazolidine proton as a multiplet at 7.58-7.25 ppm. For IST-04, 1 H-NMR spectra showed N-H proton at 13.24 ppm, thiazoline proton at 7.79 ppm. 13  The density functional theory (DFT) is an important theoretical method that could help in a better understanding of chemical and optical properties of compounds and their complexes. The electron density distributions and energy levels of isatin complexes with HOMO and LUMO are given in Figure 4. Theoretically, the orbital energy difference between HOMO and LUMO is named the HOMO-LUMO gap (DH-L). The optical band gap (E g ) is the difference between the ground (S 0 ) and excited state (S 1 ) rather than the difference between orbital energies of HOMO and LUMO. Experimentally, the optical band gap (E g ) obtained from the spectra is excitation or the lowest transition energy from the ground state to the excited state (first dipole state), considering only an electron, which is stimulated from the HOMO to the LUMO. The HOMO-LUMO energy gap decreases in the order as: IST-01 = IST-03 > IST-04 > IST-02.
The optimization energy and DH-L energy for IST-02 is the least, which demonstrates that it is the most reactive complex amongst other. The energy values for the S 0 and S 1 are shown in Table 1. The attached electron-donating groups, thiazole and 4-oxo-2thioxoimidazolidin, withdraw the electron density toward the middle of the substituents. Similarly, for the complex IST-01 and IST-04, electron-donating group thiazole pulls the electron density towards the substituent. From the DFT studies, it can be concluded that the IST-02 conduction of electrons from the ground state to the excited state is comparatively easier due to the least energy gap. Moreover, the increase in electron-withdrawing capability results in lowering of the HOMO-LUMO gap.

Molecular Docking Studies
To grasp the idea of modes of DNA binding and to identify the sites where a compound could bind in major or minor grooves of DNA, molecular docking is implemented as a valuable technique in theoretical drug design [33,34]. Molecular docking was performed by subjecting DNA to isatin compounds, with a number of docking poses per compound where each pose was assigned a score by MOE. The conformational search has been performed to identify the most favorable binding conformation. Best poses having most favorable interactions for all the derivatives were selected and shown in Figure 5. IST-01 furnished no hydrogen bonding interaction but hydrophobic interactions with the core of DNA, whereas the other three compounds developed H-bonding with the DNA base pairs. IST-02 developed a single H-bond with its receptor. However, IST-03 and IST-04 exhibited 2 H-bonds with the DNA base pairs. From the interaction pattern and binding constants (K b ) calculated through molecular docking, it can be established that hydrogen bond interactions along with hydrophobic interaction are responsible for the highest K b value of the IST-02. When the MOE predicted conformation of IST-02, IST-03, and IST-04 was matched with the conformations obtained from an X-ray diffractometer, the hydrogen bonding mode with DNA was found to be the same.  UV-visible spectrophotometry is a prime procedure to investigate compound-DNA binding. The variations in the UV-visible spectrum of a compound upon DNA addition could help identify the interactions involved between a compound and DNA.
The four isatin derivatives (IST-01-IST-04) have been evaluated spectrophotometrically to obtain an understanding of their interactions with ds-DNA. The spectra of all the derivatives were recorded, separately, at different concentrations and shown in Supplementary Materials as Figure S13, while graphs were drawn for the absorbance at various concentration and provided in Supplementary Materials as Figure S14. The molar extinction coefficients (ε; M −1 cm −1 ) were evaluated from the slopes and the values obtained for IST-01-IST-04 were found to be 38,933, 20,750, 42,767, and 14,533, respectively. These values predicted π-π* transitions in all compounds in the λ max range of 380-590 nm.
For DNA binding interactional studies, DNA titrations were performed, separately, into the fixed concentrations of each compound under physiological pH (7.4) and temperature (37 • C). The spectral responses before and after DNA titrations are shown in Figure 6. The peak intensity was observed as the drop in absorbance and the trend was found to be similar for all compounds after DNA titration. The percent drop in the absorption peak intensity was calculated by the following equation.
The resultant hypochromic effect was found to be 84.94% for IST-01, 65.88% for IST-02, 83.22% for IST-03, and 75.50% for IST-04 followed by a blue shift of magnitude 4.5 nm, 1 nm, and 5.5 nm for IST-01, IST-03, and IST-04, respectively, and a red shift of magnitude 10 nm for IST-02. Generally, hypochromism and a significant blue/red shift are the parameters that arise due to the change in the structure of DNA, with it reported for the intercalation mode of compound-DNA binding [3,[35][36][37][38][39]. The spectral observation and related data showed that all four compounds interacted with DNA via intercalative binding [3,36,37].
Generally, the Benesi-Hildebrand and van't Hoff equations were used to evaluate the binding constant and Gibbs free energy change, respectively.
Here, A o and A are the absorbance of the compound in the absence and presence of DNA, respectively, and ε G and ε H-G are the molar extinction coefficients of free compound and the compound-DNA complex, respectively. The binding constant (K b ) was calculated from the intercept-to-slope ratio from the plot of  Table 3. The order of magnitude for K b was evaluated to be 10 3 which reflects spontaneous and substantial binding of the compounds with DNA, with the same order of binding constant reported for some intercalating molecules [3,38,39].    Similar DNA titrations were carried out under identical conditions (pH 7.4, 37 • C) as those set for UV-visible experiments. All compounds were found to be luminescent; hence, their fluorescence spectra were recorded separately and then DNA titrations were carried out by adding different DNA concentrations directly into the fixed concentration of each compound. The fluorescence responses are shown in Figure 8. A remarkable rise in peak intensity was observed in the compounds' spectra when DNA was added up to 90 µM. This enhancement in the emission intensity was found to be 8.88, 4.18, 3.67, and 2.54 times higher for IST-01-04, respectively, after the addition of maximum DNA concentration and revealed interaction between the tested compounds and DNA. The literature has reported such spectral changes in compound-DNA spectra as interaction via intercalative binding [3,40].
Fluorescent enhancement also indicated that the molecule is strictly shackled at the binding site and there is no hydrophobic interaction between the nitrogenous bases of the DNA and water molecules. So, after excitation, the masking of the fluorophore and the interaction between stacked nitrogenous bases of the DNA decreases the vibrational modes that contribute towards the increase in fluorescence peak intensity [40].
The fluorescence intensity data of each compound before and after DNA titration was used in the following equation to evaluate the binding constant and binding site size [41].
where 'K b ' and 'n' are the binding constant and binding site size, respectively. From the plots of log (Figure 9), 'K b 'and 'n' of each compound-DNA complex were evaluated as the antilog of the intercept and slope values, respectively, while ∆G was evaluated by using K b value in van't Hoff equation.
The values of all binding parameters are provided in Table 3. The 'K b ' values were calculated in the same order as done in UV-visible spectroscopy, i.e., IST-02 > IST-04 > IST-03 > IST-01. This confirmed the relatively more stable and stronger interaction of IST-02 with DNA. As the binding site size for all compound-DNA complexes was greater than 1 (n > 1), it showed that along with intercalation, other types of reversible interactions such as groove binding or electrostatic interaction may be present due to more site availability. Gibbs free energy changes (Table 3) were also in agreement with the values calculated in UV-visible spectroscopy.

DNA Binding Studies by Cyclic Voltammetry
DNA binding predictions by spectroscopies (UV-visible, Flu-) could further be complemented by cyclic voltammetry that describes compound-DNA interaction in terms of changes in fundamental electrochemical parameters (current/potential). The variations in the cyclic voltammograms of all compounds after DNA titrations, at physiological conditions (pH 7.4, 37 • C), are shown in Figure 10. At first, the individual scanning was performed for IST-01-04 (−1.5 to +1.5V) and it was found that each compound exhibited irreversible cyclic process with a reduction peak within the scan range of −0.5 to −1V. The system irreversibility remained unchanged after DNA titrations for each compound. However, the reduction peak current significantly dropped down to 40.61%, 24.21%, 23.47%, and 20.74%, respectively, along with a positive shift in peak potential. Such voltametric observations showed that all tested compounds interacted with the DNA. The reduction peak current changes were further implied in the following equation for the calculation of binding constants of all tested compounds [41].
where I po and I p represents the peak current in the absence and presence of DNA, respectively. The binding constant (K b ) was determined by using the above equation in a graph of I p 2 vs. I po 2 − I p 2 /[DNA] (Figure 11). The evaluated values of K b , along with ∆G, are displayed in Table 3.  (Figure 12). The evaluated values are given in Table 3 and are consistent with those obtained from fluorescence data. All binding parameters (K b , ∆G, n) obtained spectroscopically and electrochemically complimented each other and indicated that all compounds have a significant attraction for DNA via mixed mode, with preferably greater binding interactions of compound IST-02 and then IST-04 with the DNA. To further assure the compound-DNA complex formation, the irreversible redox activity of all compounds were scanned at different scan rates (0.1-0.13 V/s), with and without DNA. The overlay of cyclic voltammograms at various scan rates are provided in Supplementary Materials as Figure S15. The current variations with the square root of scan rates were plotted for all compounds and their DNA-bound complexes. The graphs are shown in Supplementary Materials as Figure S16. The data was further implemented in Randles-Sevcik Equation (7) The calculated values of the diffusion coefficient are presented in Table 4. They clearly indicate slow diffusion, with lower diffusion coefficient values for compounds with DNA due to the formation of bulky compound-DNA complexes [42].

DNA Binding Studies by Viscometry
Another method to investigate compound-DNA interaction is to determine the viscosity of DNA in the presence of a compound's concentration. Viscosity of DNA increases when the planar portion of a drug/drug-like molecule is inserted into the DNA binding pockets and stacked between the base pairs resulting in the widening of the DNA helical structure. This insertion of the planar part is referred as intercalation and the increase in DNA viscosity arises due to the interaction of electronic density of the intercalating part of a compound with the stacked DNA base pairs. DNA is a natural biopolymer, hence the changes in DNA viscosity will be in minute fractions due to the smaller size of the inserted molecule in comparison to DNA. A mixed binding mode is predicted if a constant behavior is observed after a linear rise in DNA viscosity [36]. On the other hand, if DNA becomes denatured, then a drop in viscosity is observed [43]. Denaturation can occur by a variation in pH, rise in temperature or by using organic solvents such as DMF and DMSO [43].
The viscosity of a DNA solution was first measured at its fixed concentration of 10 µM and then after the addition of 10 to 90 µM of each compound, separately. These changes in DNA viscosity were shown in (η/η o ) 1/3 vs.
[Compound]/[DNA] plots for all compounds ( Figure 13). They indicated that, initially, the relative specific viscosity increased linearly before becoming constant after certain concentrations of the compound. The graphical trend in viscosity verified that all compounds exhibit mixed binding modes for compound-DNA interaction, further supporting the docking investigations for binding interactions via partial intercalation and groove binding.

Cytotoxicity Analysis by Cell Line Studies
The in vitro anticancer potential of the tested compounds IST-01-04 was carried out on three cancerous cell lines, namely HeLa (cervical carcinoma), MCF-7 (breast carcinoma), and HuH-7 (liver carcinoma). None of the compounds showed significant cytotoxic activity against cervical (HeLa) and breast cancer cells (MCF-7), suggesting that these drugs might not be potent against these specific cancers (Figures S17 and S18 in Supplementary Materials). However, among all compounds, IST-02 and IST-04 have shown significant cytotoxic activity against the liver cancer cell line (HuH-7), shown in Figure 14. IST-02 has shown more than 50% cytotoxicity against HuH-7 at 50 ng/µL to 500 ng/µL concentrations, while IST-04 has shown more than 50 percent cytotoxicity at all tested concentrations, i.e., 10 to 500 ng/µL (Table 5).  The cytotoxic effects of IST-02 and IST-04 were also tested on normal HEK-293 (human embryonic kidney) cell line, with both compounds also showing dose-dependent cytotoxic activity against normal HEK-293 cells ( Figure 15). However, the percentage cytotoxicity remained less than 50% at concentration range of 10 to100 ng/µL (p < 0.001) for both the compounds on HEK-293 cell line. At the same concentrations, the percentage cytotoxicity on the HuH-7 cancer cell lines was greater than 50%, indicating the potential of these drugs being selectively cytotoxic against hepatocellular carcinoma cells (Table 5). The half maximal inhibitory concentration (IC 50 ) values, as given in Table 5, also indicated that both IST-02 and IST-04 have greater anticancer potency against HuH-7 hepatoma cell line. IC 50 values indicate that a specific dose of a compound is essential to kill 50% of the cancerous cells. Most of the time, IC 50 values for a cancerous cell line are less, indicating that a lower concentration of drug is needed to kill the cancerous cells, whereas larger values of IC 50 on a healthy cell line indicate that a higher dose of drug is needed for the death of healthy cells. This type of phenomenon corresponds to the least side effects of a drug. In the present research, IST-02 and IST-04 had IC 50 value of 3.07 ± 9.47 and 14.60 ± 2.49 µM, respectively, which suggests that these drugs are sufficiently powerful enough to kill cancerous cells at lower concentrations. For healthy cells, the values of IC 50 for these drugs were 488.30 ± 2.44 and 197.80 ± 8.75 µM, respectively, which indicates lesser cytotoxicity for healthy cell.
The IC 50 values of isatin-dehydroepiandrosterone conjugates (4a, 4d, and 4e) were reported as 28.19 ± 3.59, 13.90 ± 3.91, and 22.26 ± 6.21µM, respectively, against HuH-7 [44]. In another recent study on triazole-linked 1,4-dihydropyridine-isatin scaffolds (N1, N2, N13), the IC 50 values were evaluated as 13.66, 6.73, and 19.36 µM, respectively, for the same cell line [45]. The extract of Isatis indigotica leaves and indole alkaloids (Ih, Ij) were tested on the same cell line and the IC 50 values were reported as 58 µM [46] and 4.88 ± 0.78 and 6.60 ± 1.20 [47], respectively. The IC 50 values evaluated for isatin derivatives in current studies have been low compared to the literature cited above, especially for the compound IST-02. These results indicate that our compounds are more potent than previously reported derivatives/compounds against liver cancer. To further validate these in vitro results and to authenticate the anticancer potential of tested compounds, in vivo experiments on animal models of respective cancers are required in future.

Materials and Methods
Analytical grade reagents and chemicals were used throughout the experimental procedures. Solvent purification and drying were performed using standard methods. Tetramethyl silane (TMS) was used as an internal reference and each compound and TMS were dissolved, separately, in deuterated solvents to record NMR spectra. For the expression of chemical shift values, the parts per million (ppm) unit was used. For binding studies, salmon sperm DNA was used as received from Sigma-Aldrich. To prepare the stock solution of DNA, deionized water was used. The DNA solution was further diluted to record its absorbance at λ max of 260 nm on a UV-visible spectrophotometer. The Beer-Lambert law was used (C = A 260 /ε × l, where ε is 6600 cm −1 M −1 and l is 1 cm), and the concentration of DNA was determined to be 5 × 10 −5 M. The absorbance ratio A 260 /A 280 was calculated as 1.6, which ensured that the DNA solution was free of other cell organelles including proteins and hence was pure for DNA binding studies [48].

Instrumentations
Gallenkamp melting point apparatus was used to determine the melting point via open-end capillary technique. IR analysis was performed on a Bruker FT-IR Bio-Rad-Excalibur spectrometer bearing series mode no. FTS 300 MX. To acquire 1 H-NMR of the synthesized compounds, a Bruker 300 MHz NMR spectrometer was used. A 75MHz NMR channel was used to obtain 13 C-NMR decoupled spectra of compounds using proton decoupling frequency. For DNA binding titration experiments, a UV-visible spectrophotometer (Shimadzu-1800, Series mode no. TCC-240 manufactured in Japan), fluorescence spectrophotometer F-7000 model FL2133-007, and AUTOLAB PGSTAT-302 (Metrohm, manufactured in Netherlands) version 4.9 of GPES software were used. The cell temperature of each instrument was controlled at 37 • C with the help of a temperature monitoring device. Wavelength scan, emission scan mode, 2400 nm/s scan speed, 0.0 s time delay, 400 V PMT voltage, 0.5 s response time, 5.0 nm each EX-slit and EM slit, 200.0 nm EM start wavelength, and 900 nm EM end wavelength were the instrumental parameters used in the fluorescence spectrophotometer. For cyclic voltammetric investigations, a conventional double-walled cell was used containing 99.9% pure platinum wire with 0.5 mm diameter as the counter electrode, glassy carbon with diameter 0.3 mm as the working electrode, and a saturated calomel (SCE) filled with 3.5 M potassium chloride (KCl, saturated) as the reference electrode. To obtain a shiny surface of glassy carbon electrode, soft rubbing of the electrode was performed on alumina slurry. It was later washed with distilled water and further sonicated for 30-45 s to remove any type of contamination. A digital viscometer bearing model no. AVS 310 (Schott Gerate automated) was used for viscosity measurements.

Synthesis of Isatin Derivatives
The synthesis protocol for the preparation of 5-chloroisatin-based heterocycles is shown in Scheme 1. In the first step, 5-chloroisatin (1) was treated with benzyl bromide (2) to synthesize the product (3) which was further reacted with thiosemicarbazide (4) to obtain Schiff base product (5). In the last step, Schiff base product was reacted with alpha-halogenated compounds to yield desired compounds 1-benzyl-5-chloro-3-(2heteroarylhydrazono) indoline-2-ones (IST-01, IST-02, IST-03, IST-04). The first reaction was carried out in anhydrous DMF under an inert atmosphere with a constant reflux of 12 h, while the second step was accomplished in the presence of dry ethanol and under controlled temperature and reflux of 3 h. The problem of side products formation in the final step was minimal as the reaction was carried out in dry distilled ethanol. The reactions were supervised by aluminum TLC plates having silica gel coated with F 254 indicator using mobile phase ethyl acetate:n-hexane (4:6). All synthesized derivatives were analyzed by FT-IR, 1 H NMR, and 13 C NMR spectroscopic techniques.

Synthesis of Isatin Derivatives
The synthesis protocol for the preparation of 5-chloroisatin-based heterocycles is shown in Scheme 1. In the first step, 5-chloroisatin (1) was treated with benzyl bromide (2) to synthesize the product (3) which was further reacted with thiosemicarbazide (4) to obtain Schiff base product (5). In the last step, Schiff base product was reacted with alphahalogenated compounds to yield desired compounds 1-benzyl-5-chloro-3-(2-heteroarylhydrazono) indoline-2-ones (IST-01, IST-02, IST-03, IST-04). The first reaction was carried out in anhydrous DMF under an inert atmosphere with a constant reflux of 12 h, while the second step was accomplished in the presence of dry ethanol and under controlled temperature and reflux of 3 h. The problem of side products formation in the final step was minimal as the reaction was carried out in dry distilled ethanol. The reactions were supervised by aluminum TLC plates having silica gel coated with F254 indicator using mobile phase ethyl acetate: n-hexane (4: 6). All synthesized derivatives were analyzed by FT-IR, 1 H NMR, and 13 C NMR spectroscopic techniques. In this study, all calculations were performed with the Gaussian 09 program using its default criteria [49]. The equilibrium structures were fully optimized without any symmetry constraints by DFT methods at B3LYP level of the theory using 6-311++G(d,p) basis sets. A conformational analysis was performed to obtain the lowest energy structures of the isatin derivatives. The electronic parameters including HOMO, LUMO, and HOMO-LUMO band gaps were determined as a difference of highest occupied molecular orbitals (HOMO) and lowest unoccupied molecular orbitals (LUMO) energies. Output structures of isatin derivatives were visualized using GaussView.

Molecular Docking
Docking studies were performed on the crystal structure of berenil-dodecanucleotide complex at a resolution of 2.5 Å, as obtained from Protein Data Bank (PDB) code 2DBE. Structure energy was minimized by using Amber10 forcefield. The active site of DNA was selected using MOE 2015. The selected binding region contained all the residue of the active site reported in the up-to-date literature. Moreover, the compounds were also prepared by refining the chemical correctness (3D protonation), ionization, and stereo-chemical variation. The energy minimization was achieved at physiological pH of 7.4. Molecular docking of inhibitors within the active site generated maximum 5 poses per ligand where each pose was assigned a score by MOE. The pose having highest score has the most favorable conformation for the ligand. The best pose based on the score was selected.
3.6. Experimental Procedures for Compound-DNA Binding 3.6.1. UV-Visible and Fluorescence Spectrometric Experiments UV-visible and fluorescence spectrophotometric investigations were performed for each compound (IST-01-04), individually, and then in presence of DNA concentrations at physiological pH (7.4) and temperature (37 • C) [36,37,[40][41][42]. The individual spectrum of each compound was recorded at optimized concentration (5 × 10 −6 M), then DNA (10-80 µM) titrations were carried out with each compound, separately. Prior to each spectral run, the establishment of an equilibrium for compound-DNA adduct was assured by keeping the reaction mixture in a cuvette to rest for at least 4-5 min at physiological temperature (37 • C).

Cyclic Voltammetric Experiments
Cyclic voltametric (CV) experiments were run under similar conditions as set in photophysical titration experiments. CV responses of the tested compounds before and after DNA titrations were recorded within the potential range of −1.5 to +1.5 V at 100 mV/s scan rate. For the diffusion coefficient calculation of each compound and its compound-DNA complex, the scanning of the CV responses was performed in the scan range of 30-150 mV/s, with a difference of 20 mV/s. Oxygen was completely flushed out from the system by bubbling argon gas for 5-8 min, while to maintain the cell at 37 • C, a continuous water circulation within the double wall of the cell was provided by a water-circulating bath at this temperature.

Viscosity Experiments
DNA viscosity (η o ) was measured at a fixed concentration of 10 µM. In contrast to UV-visible, Flu-and CV experimental work, the titration was carried out by using varying concentrations of each compound, separately, in the DNA solution. Changes in DNA viscosity were measured as (η). The relative viscosities were calculated from η o and η values and their cube roots were plotted against compound and DNA concentration ratios.

Cell Line Assays
Cell line studies of tested compounds (IST-01-04) were carried out on three cancerous cell lines, namely cervical cancer (HeLa), breast cancer (MCF-7), and liver cancer (HuH-7), and one normal human embryonic kidney cell line (HEK-293). The percentage cytotoxicity of all compounds for cancerous as well as healthy cell lines were evaluated by dosedependent MTT analysis [4,5,32]. The Huh-7 and HeLa cells were grown in DMEM (high-glucose medium) supplemented with 10% fetal bovine serum (FBS) and 1% penicillinstreptomycin (pen-strep), whereas the HEK-293 and MCF-7 cells were cultured in RPMI (high-glucose medium) supplemented with 10% fetal bovine serum (FBS) and 1% penicillinstreptomycin (pen-strep). Exponentially growing cells were counted and 10,000 cells per well were plated, in triplicate, in flat-bottomed 96-well plates (Nunc, Roskilde, Denmark). The volume of the cells was kept at 100 µL per well. Each tested compound was dissolved in supplemented media and 10% DMSO, separately, to obtain different concentrations (1.0 ng/µL, 10 µg/µL, 50 ng/µL, 100 ng/µL, and 500 ng/µL). Each concentration of the drugs was added to the 96-well plate, to obtain a final volume of~200 µL/well. In addition, each concentration was tested in triplicate on all HeLa, MCF-7, HuH-7, and HEK-293 cell lines. Control wells contained solvent control (without drug) and blank media (without cells). Subsequently, 5 mg/mL of MTT was dissolved in 1 mL PBS. Accordingly, 15 µL of prepared MTT solution was added to each well and incubated for 3 h at 37 • C, so that intracellular purple formazan crystals became visible under microscope. Following the formation of formazan crystals, all the solution from each well was removed. Then solubilizing solution, i.e., 150 µL DMSO, was added in each well. The plates were left at room temperature for a few minutes while DMSO solution was mixed thoroughly by pipetting up and down to dissolve the formazan crystals. Finally, the absorbance of the cells was measured by spectrophotometer at 550 nm. From the obtained data of absorbance, % cell viability was calculated using the following equation.

Conclusions
Four new isatin derivatives, IST-01-04, have been synthesized and characterized for their structures by various spectroscopic techniques and further examined for DNA binding and cytotoxicity activities for their anticancer potential candidacy. Results obtained from molecular docking and DFT exhibited the receptive nature of IST-02 and IST-04 and their relatively greater binding nature than IST-01 and IST-03. The computed electronic parameters predicted higher charge transfer and hence reactivity for IST-02 and IST-04 due to smaller HOMO-LUMO gap. Theoretical findings were further confirmed by experimental DNA binding investigations by spectroscopies (UV-visible, Flu-) and cyclic voltammetry (CV) under physiological pH of 7.4 and at 37 • C. The trends in spectral and voltametric responses, as well as evaluated binding parameters (K b , ∆G), indicated that among all compounds, IST-02 and IST-04 have greater and more spontaneous binding affinity for DNA. Binding site size (n) calculations from Flu-and CV data depicted more sites' availability for interactions and confirmed a mixed type of interactions, as predicted in docking analysis. DNA viscosity measurements in the presence of tested compounds further authorized our findings of mixed binding mode, as assured from theoretical and experimental DNA binding studies. In vitro cell line studies of the compounds revealed that IST-02 and IST-04 have better and selective anticancer potential on hepatocellular carcinoma cells (HuH-7) compared to normal cell lines. Furthermore, these compounds have lower IC 50 values compared to previously reported derivatives of the same class, indicating their high potency as potential anticancer drug candidates. All findings were in good agreement with each other and thus merit further in vivo testing of these compounds. The overall DNA binding and cytotoxicity results were found to be more significant for the compound IST-02.