Crystal Structures, Molecular Docking and In Vitro Investigations of Two 4-Substituted 2-(5,5-dimethyl-3-styrylcyclohex-2-enylidene)malononitrile Derivatives as Potential Topoisomerase II Inhibitors

: Type II topoisomerases (TOP2s) play a key role in altering the DNA topology by transiently cleaving both strands of a DNA duplex. Therefore, increased TOP2 activity is associated with many cancers. Herein, we present the synthesis, structural characterization, virtual screening, and structural exploration, as well as evaluation of the antiproliferative effects of two new 4-substituted 2-(5,5-dimethyl-3-styrylcyclohex-2-enylidene)malononitrile derivatives with potential application in the drug design of isoform-specific TOP2 inhibitors. Both compounds 1 and 2 were verified by ESI-TOF-MS, NMR, and single-crystal X-ray diffraction (SCXRD) analysis. Furthermore, we applied our recently proposed SCXRD/HYdrogen DEsolvation (HYDE) technology platform in order to perform molecular modeling, virtual screening, and structural exploration with 1 and 2 . For this purpose, we used the crystal structure of human TOP2 β complexed to DNA and the anticancer drug etoposide. Moreover, we further evaluated the antiproliferative activity of 1 and 2 on human hepatocarcinoma HepG2 cells and compared the observed effects with those of the reference hTOP2 β inhibitor etoposide. Based on the obtained results, compounds 1 and 2 showed a virtually higher binding affinity ( K i HYDE values) over etoposide towards hTOP2 β but lower antiproliferative activity compared to those of etoposide.


Introduction
Topoisomerases type II (TOP2s) are crucial enzymes in altering the DNA topology by cleaving both strands of a DNA duplex with the active participation of Mg 2+ and energy from ATP hydrolysis [1][2][3][4].Human TOP2s comprise two isoforms, topoisomerase IIA (TOP2α) and IIB (TOP2β) [1][2][3].TOP2s play a key role in the formation of a TOP2-DNA cleavage complex through an enzyme-bridged "DNA-gate", leading to enzyme-mediated DNA breakage and the ultimate death of cancer cells [1,2].TOP2α is involved in DNA replication, recombination, and chromosome condensation and segregation, while the TOP2β isoform regulates DNA transcription [1][2][3].Therefore, TOP2s are recognized as biological targets of several approved anticancer drugs [1,4,5].Etoposide (brand name Vepesid ® ), for example, is a TOP2 poison targeting both TOP2 isoforms [1,5,6].Etoposide leads to increased levels of TOP2-DNA covalent complexes via trapping and drug-induced stabilization, generating DNA strand breaks [4,5].It is widely used in the treatment of different types of cancer, including prostate cancer, lymphoma, myelocytic leukemia, lung Crystals 2024, 14, 496 2 of 13 cancer, neuroblastoma, and others [4][5][6][7][8][9].Other examples of chemotherapeutics able to fill up the spaces within the nucleotides in the double DNA helix include doxorubicin and idarubicin [10,11].Teniposide, an analog of etoposide, was discovered and clinically tested for its anticancer activity before etoposide [4,12,13].Teniposide has a similar mechanism of action to that of etoposide, causing DNA double-stranded breaks via the stabilization of the DNA-TOP2 complexes [4,14].Although teniposide was approved by the FDA in 1992, nine years after the FDA approval of etoposide, the adverse effects caused in patients were the reason for its slower development as a drug [4,15].
In the present study, we aim (i) to estimate the potential of two new malononitrile derivatives (designated as compounds 1 and 2) as inhibitors of human TOP2β and (ii) to evaluate their antiproliferative effects on the human hepatocarcinoma HepG2 cell line in comparison with the anticancer drug etoposide (for structures, see Figure 1).For this purpose, we applied our recently proposed single-crystal X-ray diffraction (SCXRD)/HYdrogen DEsolvation (HYDE) technology platform to design new isoform-specific hTOP2-targeting agents [16].Therefore, compounds 1 and 2 were synthesized and structurally characterized.The single-crystal X-ray structures of both compounds were further used to perform initial molecular docking and virtual screening of the generated conformers in order to predict their binding affinity (so called K i HYDE values) towards hTOP2β and the thermodynamic profile of binding.In addition, the structural exploration of the basic scaffold of 1 and 2 was conducted with the aim of virtually generating a set of optimized hTOP2 inhibitors.Finally, we investigated the antiproliferative effects of 1 and 2 on HepG2 cells and compared them with those observed for the reference anticancer agent etoposide.
Crystals 2024, 14, x FOR PEER REVIEW 2 of 14 Vepesid ® ), for example, is a TOP2 poison targeting both TOP2 isoforms [1,5,6].Etoposide leads to increased levels of TOP2-DNA covalent complexes via trapping and drug-induced stabilization, generating DNA strand breaks [4,5].It is widely used in the treatment of different types of cancer, including prostate cancer, lymphoma, myelocytic leukemia, lung cancer, neuroblastoma, and others [4][5][6][7][8][9].Other examples of chemotherapeutics able to fill up the spaces within the nucleotides in the double DNA helix include doxorubicin and idarubicin [10,11].Teniposide, an analog of etoposide, was discovered and clinically tested for its anticancer activity before etoposide [4,12,13].Teniposide has a similar mechanism of action to that of etoposide, causing DNA double-stranded breaks via the stabilization of the DNA-TOP2 complexes [4,14].Although teniposide was approved by the FDA in 1992, nine years after the FDA approval of etoposide, the adverse effects caused in patients were the reason for its slower development as a drug [4,15].
In the present study, we aim (i) to estimate the potential of two new malononitrile derivatives (designated as compounds 1 and 2) as inhibitors of human TOP2β and (ii) to evaluate their antiproliferative effects on the human hepatocarcinoma HepG2 cell line in comparison with the anticancer drug etoposide (for structures, see Figure 1).For this purpose, we applied our recently proposed single-crystal X-ray diffraction (SCXRD)/HYdrogen DEsolvation (HYDE) technology platform to design new isoform-specific hTOP2-targeting agents [16].Therefore, compounds 1 and 2 were synthesized and structurally characterized.The single-crystal X-ray structures of both compounds were further used to perform initial molecular docking and virtual screening of the generated conformers in order to predict their binding affinity (so called Ki HYDE values) towards hTOP2β and the thermodynamic profile of binding.In addition, the structural exploration of the basic scaffold of 1 and 2 was conducted with the aim of virtually generating a set of optimized hTOP2 inhibitors.Finally, we investigated the antiproliferative effects of 1 and 2 on HepG2 cells and compared them with those observed for the reference anticancer agent etoposide.
NMR measurements: 1 H and 13 C spectra were recorded on a Bruker Avance III HD Ascend 500 MHz spectrometer (Bruker Corp., Billerica, MA, USA).All NMR spectra were acquired at room temperature (298 K) in DMSO-d 6 as a solvent.Chemical shifts (δ) were reported in parts per million (ppm) relative to the residual solvent peaks in the corresponding spectra: DMSO-d 6 δ 2.50 (for 1 H) and 39.51 (for 13 C).The coupling constants J were given in Hertz (Hz), while the spin multiplicities were given as singlet (s), doublets (d), doublet of doublets (dd), triplet (t), doublet of triplets (dt), quartet (q), and multiplet (m).

Single-Crystal X-ray Diffraction (SCXRD)
A single crystal of compound 1 or 2 was examined on a Rigaku Supernova diffractometer using Cu Kα (for 1) or Mo Kα (for 2) radiation, and data were collected at 100 or 95 K, respectively (Table 1).Data were processed using the Olex2 suite [18], and the structures were solved with the ShelXT structure solution program [19] using intrinsic phasing.The structure of compound 1 was subsequently refined as a racemic twin with the ShelXL [20] refinement package using Least Squares minimization.The disorder of the S and R enantiomers of the carbon atom C14 had a ratio of 82:18.Suitable restraints and con-straints were applied.The structure of compound 2 was refined with the olex2.refine[21] refinement package using Gauss-Newton minimization.Hydrogen atoms were refined anisotropically.Subsequent refinement was performed using No-spherA2, which is the implementation of non-spherical atom-form factors in Olex2 [22].The CCDC deposition numbers 2350475 and 2350476 contain the supplementary crystallographic data for both compounds in this paper and can be obtained on request from the Cambridge Crystallographic Data Centre (CCDC).

Molecular Modeling 2.4.1. Protein Preparation and Ligand-Protein Docking
All molecular modeling studies were performed using the SeeSAR software tool [23].The preparation of the protein binding site and subsequent molecular docking were carried out by applying the respective integrated modules, as described previously [16].The X-ray structure of human topoisomerase II beta (hTOP2β) in complex with DNA and the anticancer drug etoposide (PDB accession code: 3QX3) [2] was used as an input 3D structure for molecular docking experiments.Based on an initial HYDE visual assessment for the two crystallographically recorded ligand positions (EVP_A_1 vs. EVP_D_1) in the etoposide-hTOP2β-DNA complex, we chose hTOP2β chain β, containing DNA chains D, E and F and ligand EVP_D_1 for all further computations.The ligand (e.g., EVP_D_1) possesses better estimated binding affinity compared to chain α (EVP_A_1 vs. EVP_D_1, mM vs. µM range).The single X-ray structures (CCDC dep.no.2350475 and 2350476) were used as 3D input ligands for docking experiments (e.g., the major occupied part of compound 1 and only one molecule of compound 2).The 3D structure of each ligand was processed using the integrated standard docking module in SeeSAR and the generated poses (max.10) for each ligand-protein complex were then post-processed with the HYdrogen DEsolvation (HYDE) algorithm in SeeSAR.After the visual inspection and re-scoring assessment, the best-scored docking poses of each ligand were then selected for further discussions, e.g., estimation of binding affinities, binding modes, and thermodynamic profiling [16].Further structural optimizations using the molecular editor module in SeeSAR and subsequent re-scoring were not required.The obtained ligand's best-scored poses showed high quality in terms of torsion/bindings and intra-and intermolecular clashes (e.g., HYDE estimated green points).

HYDE Scoring and Visual Assessment
The scoring function HYDE is used for the rapid prediction of binding affinities (∆G, in kJ/mol), as implemented in SeeSAR [23].Considering the hydrogen bonds/salt bridges (approximately interpreted as the hydrogen term ∆H) and dehydration terms (approximately interpreted as the desolvation term-T∆S), SeeSAR visualizes the (HYDE-)estimated free energy of binding ∆G, in kJ/mol) using translucent spheres (HYDE "coronas") that are colored from deep red (very unfavorable) to dark green (very favorable for affinity [16,[24][25][26][27].As a result, SeeSAR estimates the binding affinity (represented as K i HYDE ranges) and other relevant physicochemical parameters such as the approximate lipophilic ligand efficiency (LLE HYDE ), ligand efficiency (LE), and others [16].Therefore, based on the obtained HYDE score, the docked poses of ligands were ranked and visually inspected under consideration for their torsional quality and intra-and intermolecular clashes [25].

Evaluation of Cell Viability
The cytotoxicity and antiproliferation of the reference etoposide and the investigated compounds 1 and 2 (100 mM stock solutions in DMSO with final concentration ≤0.5% v/v) was evaluated on the human hepatocellular carcinoma (HepG2; ATCC, Manassas, VA, USA) cell line using a colorimetric MTT assay with slight modifications [28].Cells in their exponential phase of growth were seeded into 96-well flat-bottom plates (0.5 × 10 5 cells/mL) at a final volume of 100 µL/well and incubated overnight before treatment with the test substances.All compounds were tested in the concentration range of 1.0-500 µM and incubated for 24 or 72 h at 37 • C under a 5.0% CO 2 atmosphere.After the incubation period, the medium was replaced with a freshly prepared one, the MTT solution (0.5 mg/mL) was added, and the cells were incubated for 180 min.Then, the medium was removed, and DMSO was added to each well, and the plates were placed on a shaker at room temperature until complete dissolution of the purple formazan product [28].The measurements of formazans, which are produced from the biological sample by a reduction of MTT, were performed on a microplate ELISA reader Varioscan TM LUX (Thermo Fisher Scientific Inc., Waltham, MA, USA) at 550 and 570 nm.The cytotoxicity (for etoposide) and antiproliferation (all tested compounds) determined by the MTT assay were expressed as % cell viability using the following equation: where OD sample , OD blank and OD control are the measured absorption of the respective test, blank, and control sample.Statistical analysis and result graphs were performed and created using GraphPad Prism 9.0 (GraphPad Software, La Jolla, CA, USA).The IC 50 values of etoposide and compound 2 at 72 h were obtained by non-linear regression analysis using the equation log(inhibitor) vs. normalized response-variable slope.The results are presented as the mean % of the untreated controls ± standard deviation (SD) from three independent experiments (n = 3).The data from the MTT tests were analyzed by one-way ANOVA followed by Dunnett's multiple comparison tests vs. the untreated control group (Ctrl.), and p values ≤ 0.05 were considered statistically significant.

Synthesis
The target compounds 1 and 2 were synthesized following a standard procedure of Knoevenagel condensation reaction using commercially available 4-cyano-or 4-methylbenzaldehyde (4 or 5), respectively, and 2-(3,5,5-trimethylcyclohex-2-en-1-ylidene)malononitrile (3) as the starting materials [17] (Figure S1).The crude products were purified by column chromatography on silica gel, affording compounds 1 and 2 moderate overall yields as a red or yellow crystalline solid, respectively.The purity and structure of the compounds were verified by ESI-TOF-MS and NMR ( 1 H and 13 C) analysis (see Supporting Information, spectral data for 1 and 2).The chemical structures of 1 and 2 were later confirmed by full crystallographic characterization.
in Figure 3 (left).The intermolecular interactions in the crystal packing of compound 2 are mainly characterized by σ-π interactions and non-classical C-H…N hydrogen bonds with 2.45(1) Å.The former leads to the dimerization of the two molecules 1 and 2 that form an asymmetric unit (Figure S2B).The shortest distances between the hydrogen atoms H2B and H22A and the centers of the benzene rings C9-C14 (for molecule 1) and C29-C34 (for molecule 2) are 2.391 Å and 2.41(1) Å, respectively (Figure 3, right).The shortest C-H…N distances indicating non-classical hydrogen bonds from N2 to H5′ with 2.426(4) Å and N3 to H19″ with 2.658(4) Å for compound 1 are shown in Figure 4 (left), in which a rigid model for the hydrogen atoms was used.In the case of compound 2, the shortest hydrogen bonds below 3.5 Å are those from N1 to H26′ and H27′, both with 2.49(1) Å, and those from N4 to N28″ with 2.44(1) Å, as indicated in Figure 4 (right).The shortest C-H. ..N distances indicating non-classical hydrogen bonds from N2 to H5 ′ with 2.426(4) Å and N3 to H19 ′′ with 2.658(4) Å for compound 1 are shown in Figure 4 (left), in which a rigid model for the hydrogen atoms was used.In the case of compound 2, the shortest hydrogen bonds below 3.5 Å are those from N1 to H26 ′ and H27 ′ , both with 2.49(1) Å, and those from N4 to N28 ′′ with 2.44(1) Å, as indicated in Figure 4 (right).

Molecular Modeling Experiments
In order to investigate whether compounds 1 and 2 are able to inhibit the activity of human TOP2 enzymes, we performed initial molecular docking studies with both proposed molecules for further structural exploration.Following our strategy, compounds 1 and 2 were originally designed within our extended research program of novel chromophores, consisting of a phenol donor (D), π-conjugated bridge (π), and dicyanomethylydiene acceptor (A) moieties, namely D-π-A structures with nonlinear optical properties (NLO) [17].However, based on our experiences in rational drug design, we decided to evaluate the biological activity of both 1 and 2 due to their structural similarity to some drugs comprising a cyanomethyl functional group.Therefore, we used the single X-ray structures of 1 (the major occupied part of the disordered S and R enantiomers at C14, ratio 82:18) and 2 (molecule 1 of both enantiomeric molecules) without further 3Doptimization as input ligands to investigate their interactions within the binding pocket of the human TOP2 β isoform (cf.Section 2).The high-resolution X-ray structure of the cleavage core of hTOP2β (herein residues 452-1201, known as hTOP2β core ) in complex with DNA and the anticancer drug etoposide (PDB code: 3QX3, res.2.16 Å) [2] was used as a structural basis for the implementation of our SCXRD/HYDE platform for structurebased drug design [7] (Figures S3 and 5).

Molecular Modeling Experiments
In order to investigate whether compounds 1 and 2 are able to inhibit the activity of human TOP2 enzymes, we performed initial molecular docking studies with both proposed molecules for further structural exploration.Following our strategy, compounds 1 and 2 were originally designed within our extended research program of novel chromophores, consisting of a phenol donor (D), π-conjugated bridge (π), and dicyanomethylydiene acceptor (A) moieties, namely D-π-A structures with nonlinear optical properties (NLO) [17].However, based on our experiences in rational drug design, we decided to evaluate the biological activity of both 1 and 2 due to their structural similarity to some drugs comprising a cyanomethyl functional group.Therefore, we used the single X-ray structures of 1 (the major occupied part of the disordered S and R enantiomers at C14, ratio 82:18) and 2 (molecule 1 of both enantiomeric molecules) without further 3D-optimization as input ligands to investigate their interactions within the binding pocket of the human TOP2 β isoform (cf.Section 2).The high-resolution X-ray structure of the cleavage core of hTOP2β (herein residues 452-1201, known as hTOP2β core ) in complex with DNA and the anticancer drug etoposide (PDB code: 3QX3, res.2.16 Å) [2] was used as a structural basis for the implementation of our SCXRD/HYDE platform for structure-based drug design [7] (Figures S3 and 5).
optimization as input ligands to investigate their interactions within the binding pocket of the human TOP2 β isoform (cf.Section 2).The high-resolution X-ray structure of the cleavage core of hTOP2β (herein residues 452-1201, known as hTOP2β core ) in complex with DNA and the anticancer drug etoposide (PDB code: 3QX3, res.2.16 Å) [2] was used as a structural basis for the implementation of our SCXRD/HYDE platform for structurebased drug design [7] (Figures S3 and 5).The performed docking and re-docking assessment procedure with monomer B (chain β) of the dimeric hTOP2β complexed to a 20-base pair DNA duplex (with chain D-F, containing a 5 ′ -C-NNNNG-3' cleavage site, DNA-gate) [2] and stabilized by the etoposide, resulted in a good overlap for the obtained docking poses with the almost identical conformation of both compounds 1 and 2 (Figures S3 and 5).The proposed docking models of 1 and 2 suggest that both compounds occupy the same ligand-binding cavity of the active site of hTOP2β (composed of catalytic WHD and Mg 2+ -chelating TOPRIM domain) [2], sharing twelve amino acid residues and eight DNA bases (Table S1).Moreover, compounds 1 and 2 formed the same contacts within the binding site of hTOP2β, involving one hydrogen bond with residue Asp479 and hydrophobic interactions with five DNA bases and four amino acid residues (Table S1 and Figure S4).The binding modes of 1 and 2 did not show any additional stabilization by water molecules within the binding site of the hTOP2β-DNA complex (Figure 5).In contrast, the binding of the etoposide to the hTOP2β active site complexed to DNA strongly indicates the involvement of seven amino acid residues and the same DNA bases as 1 and 2 in the formation of hydrophobic contacts (Table S1 and Figure S4).The main differences between the binding modes of both compounds and etoposide are observed for the formed hydrogen bonds.While compounds 1 and 2 interact via one H-bond with Asp479, the etoposide establishes two with Asp479, three with water molecules (HOH241, 364, and 1355), and one H-bond with the DNA base G13 (cf.Table S1 and Figure 5).The binding of all ligands to the hTOP2β-DNA complex also suggests the formation of H-bonds between HOH187, Asp503, and HOH46 with the coordination of Mg 2+ (Figure 5).
Using the obtained best-scored docking poses, we further assessed the thermodynamic binding profiles of 1, 2, and the complexed ligand etoposide, as well as estimated their binding affinities (the so-called K i HYDE ranges) within the binding pocket of the hTOP2β-DNA complex (Figure 6).Based on the docking protocol and the (re-)scoring procedure described recently by us [16], we applied the implemented SeeSAR algorithm HYDE to compute the enthalpic (∆H) and entropic effects (−T∆S) of all non-hydrogen atoms of the ligands to the overall energy of binding (∆G, kJ/mol).The complete HYDE visual assessment of the ligands is shown in Figure S5.The results from the HYDE assessment proposed that the main binding activity of 1 and 2 was due to the enthalpic (∆H) H-bond interaction between the ligand and the amino acid residue Asp479 with −18.3 kJ/mol for both compounds.The contribution of the hydrophobic contacts (−T∆S) within the hTOP2β-DNA complex to ∆G was more than two-fold lower with −8.10 and −7.70 kJ/mol for 1 and 2, respectively (Figure 6, top).As obtained from its best-docking pose, the binding affinity of etoposide is mainly driven by hydrogen bonds within the hTOP2β-DNA complex, but it is compensated by non-favorable hydrophobic interactions, e.g., an entropic part of ∆G of −T∆S = +44.1 kJ/mol (Figure 6, top).Considering the HYDE computed energy of binding (∆G, kJ/mol), it can be summarized that there is a good agreement between ∆G and the predicted binding affinities (expressed as K i HYDE ranges, from low to high µM) (Figure 6, bottom).Compounds 1 and 2, for example, showed almost equal ∆G of −26.4 and −26.0 kJ/mol, respectively, with K i HYDE ranges of 2.26-225 (for 1) and 2.84-282 µM (for 2) (Figure 6).Based on its HYDE binding affinity, the etoposide seems to have lower predicted K i HYDE ranges of 9.10-904 µM than 1 and 2; however, they are still in the low µM range.Taking into account the obtained docking modes for 1 and 2 within the hTOP2β-DNA complex, we performed the virtual structural exploration of both molecules in order to generate optimized compounds with improved physicochemical properties.Using the Inspirator module in SeeSAR, a set of six (for 1) and seven (for 2) new structural analogs were designed, selected on the basis of their improved physchem characteristics, and virtually screened at the active site of the hTOP2β-DNA complex (Tables S2).According to the obtained results, the introduction of a NH2 function at position meta in the phenyl ring of 1 or a fluorine atom instead of a methyl group in the dicyanomethylene moiety of 2 led to an increase in the estimated biological activity (Ki HYDE values) and physicochemical properties in the proposed structures 1_6 and 2_1 (cf.Table S2).Thus, these results suggest that the best docking poses of 1 and 2 within the active site of the hTOP2β-DNA cleavage complex can be useful for the further design of isoform-specific hTOP2 inhibitors.

Evaluation of Cytotoxicity and Antiproliferation
In order to examine the anticancer potential of compounds 1 and 2, we performed further investigation of their antiproliferative effects on hepatocarcinoma HepG2 cells.This cell line is commonly used to test cytotoxicity and the antiproliferation of investigated drugs and, in particular, for comparative studies with reference anticancer agents [28].In addition, the HepG2 cell line is also used in the preclinical safety assessment in Taking into account the obtained docking modes for 1 and 2 within the hTOP2β-DNA complex, we performed the virtual structural exploration of both molecules in order to generate optimized compounds with improved physicochemical properties.Using the Inspirator module in SeeSAR, a set of six (for 1) and seven (for 2) new structural analogs were designed, selected on the basis of their improved physchem characteristics, and virtually screened at the active site of the hTOP2β-DNA complex (Table S2).According to the obtained results, the introduction of a NH 2 function at position meta in the phenyl ring of 1 or a fluorine atom instead of a methyl group in the dicyanomethylene moiety of 2 led to an increase in the estimated biological activity (K i HYDE values) and physicochemical properties in the proposed structures 1_6 and 2_1 (cf.Table S2).Thus, these results suggest that the best docking poses of 1 and 2 within the active site of the hTOP2β-DNA cleavage complex can be useful for the further design of isoform-specific hTOP2 inhibitors.

Evaluation of Cytotoxicity and Antiproliferation
In order to examine the anticancer potential of compounds 1 and 2, we performed further investigation of their antiproliferative effects on hepatocarcinoma HepG2 cells.This cell line is commonly used to test cytotoxicity and the antiproliferation of investigated drugs and, in particular, for comparative studies with reference anticancer agents [28].In addition, the HepG2 cell line is also used in the preclinical safety assessment in order to estimate the drug's potential to provoke hepatotoxicity [30].The experiments were carried out according to our previous protocols applying the MTT assay [28].For all tested substances, the cell viability of HepG2 cells was evaluated after a 72 h incubation period in a concentration range of 1.0-500 µM (Figure 7).In the case of etoposide, the cytotoxic effect was determined after 24 h to examine its efficacy in a standard incubation period with HepG2 cells.In the tested concentration range, etoposide inhibits the growth of HepG2 cells after 24 h (about 50% inhibition at the two highest test concentrations of 250 and 500 µM) and 72 h of incubation (Figure 7a,b).Compared to etoposide, only compound 2 showed an antiproliferative effect (after 72 h treatment) against HepG2 cells at the highest-tested concentrations of 250 and 500 µM (Figure 7c,d).Dose-dependent non-linear inhibitors curves were built to obtain the respective IC50 values for antiproliferative effects of etoposide and compound 2 on HepG2 cells (Figure 7g,h).The IC50 values for etoposide and compound 2 were determined to be 7.99 ± 2.43 and 387 ± 62.1 µM (50% inhibition of cell growth at 500 µM of 2), respectively.However, the obtained results can be considered as a good starting point for further experiments since starting from these malononitrile derivatives may be performed new optimized anticancer compounds with improved physicochemical parameters, such as an optimal hydrophilic-lipophilic balance.Dose-dependent non-linear inhibitors curves were built to obtain the respective IC 50 values for antiproliferative effects of etoposide and compound 2 on HepG2 cells (Figure 7g,h).The IC 50 values for etoposide and compound 2 were determined to be 7.99 ± 2.43 and 387 ± 62.1 µM (50% inhibition of cell growth at 500 µM of 2), respectively.

Figure 1 .
Figure 1.Chemical structures, IUPAC names, chemical formula, and molecular weight of etoposide (trade name Vepesid ® ) compounds 1 and 2. The images on the top show the colored crystals of both compounds.

Figure 1 .
Figure 1.Chemical structures, IUPAC names, chemical formula, and molecular weight of etoposide (trade name Vepesid ® ) compounds 1 and 2. The images on the top show the colored crystals of both compounds.

Crystals 2024 , 14 Figure 2 .
Figure 2. Molecular structure of compounds 1 and 2 in the solid state with atomic labeling.Anisotropic displacement ellipsoids are shown at the 50% level.The minor occupied part of compound 1 is translucent, and the hydrogen atoms are omitted for clarity.

Figure 2 .
Figure 2. Molecular structure of compounds 1 and 2 in the solid state with atomic labeling.Anisotropic displacement ellipsoids are shown at the 50% level.The minor occupied part of compound 1 is translucent, and the hydrogen atoms are omitted for clarity.

Figure 5 .
Figure 5. Docking model of compound 1 (magenta, left), compound 2 (cyan, middle), and etoposide (orange, right) obtained from their best-scored docking pose into the active site of hTOP2β (β chain in green ribbons).All ligands are represented as balls and sticks with HYDE and torsional analysis, as generated in SeeSAR.The amino acid residues and DNA bases involved in H-bonding (dotted lines) and lipophilic interactions are shown as a wireframe and labeled.Mg 2+ and some water molecules are represented as white and red spheres, as indicated, respectively.The water molecules forming H-bonds are shown as colored sticks.HYDE visual assessment scheme: green = favorable, red = unfavorable, and non-colored = not relevant for binding affinity.

Figure 5 .
Figure 5. Docking model of compound 1 (magenta, left), compound 2 (cyan, middle), and etoposide (orange, right) obtained from their best-scored docking pose into the active site of hTOP2β (β chain in green ribbons).All ligands are represented as balls and sticks with HYDE and torsional analysis, as generated in SeeSAR.The amino acid residues and DNA bases involved in H-bonding (dotted lines) and lipophilic interactions are shown as a wireframe and labeled.Mg 2+ and some water molecules are represented as white and red spheres, as indicated, respectively.The water molecules forming H-bonds are shown as colored sticks.HYDE visual assessment scheme: green = favorable, red = unfavorable, and non-colored = not relevant for binding affinity.

Crystals 2024 , 14 Figure 6 .
Figure 6.Approximate computed thermodynamic profiling of 1, 2, and the etoposide.(Top): Thermodynamic profile of compounds 1, 2 and the etoposide obtained from their best docking poses into hTOP2β (PDB: 3QX3), representing the semi-quantitative decomposition of the enthalpic (ΔH: sum of interactions; in grey bars) and entropic part (−TΔS: sum of dehydration terms; in red) for all non-hydrogen atoms of the Gibbs free energy (∆G = ∆H − T∆S, kJ/mol; in black).The respective values (in kJ/mol) and binding affinities (Ki HYDE ranges) are summarized below the diagram.(Bottom): Tabular representation of HYDE-estimated affinity ranges (white bars from µM to nM), lipophilic ligand efficiencies (LLE), torsions (Tor, green points), intra-and intermolecular clashes (green points), as well as molecular weight (MW), logP, and the topological polar surface area (TPSA) as the main physicochemical and drug-like characteristics.

Figure 6 .
Figure 6.Approximate computed thermodynamic profiling of 1, 2, and the etoposide.(Top): Thermodynamic profile of compounds 1, 2 and the etoposide obtained from their best docking poses into hTOP2β (PDB: 3QX3), representing the semi-quantitative decomposition of the enthalpic (∆H: sum of interactions; in grey bars) and entropic part (−T∆S: sum of dehydration terms; in red) for all nonhydrogen atoms of the Gibbs free energy (∆G = ∆H − T∆S, kJ/mol; in black).The respective values (in kJ/mol) and binding affinities (K i HYDE ranges) are summarized below the diagram.(Bottom): Tabular representation of HYDE-estimated affinity ranges (white bars from µM to nM), lipophilic ligand efficiencies (LLE), torsions (Tor, green points), intra-and intermolecular clashes (green points), as well as molecular weight (MW), logP, and the topological polar surface area (TPSA) as the main physicochemical and drug-like characteristics.

Figure 7 .
Figure 7. Cytotoxicity of etoposide (a) and antiproliferative effects of etoposide (b), compound 1 (c), and compound 2 (d) measured on human hepatocarcinoma HepG2 cells after 24 and 72 h of exposure to different concentrations of test substances (1.0 to 500 µM).Representative images of the HepG2 cell culture after 72 h of treatment with 100, 250, and 500 µM of etoposide (e) and compound 2 (f).The nonlinear regression curve of etoposide (g) and compound 2 (h) shows the inhibition of HepG2 cells after 72 h incubation with the respective test substance (the IC50 values are indicated).The results are expressed as the mean % of untreated controls ±SD (n = 3).Statistical analysis was performed by one-way ANOVA and Dunnett's multiple comparison test.*, p < 0.1; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001 vs. the untreated control group (Ctrl.).

Figure 7 .
Figure 7. Cytotoxicity of etoposide (a) and antiproliferative effects of etoposide (b), compound 1 (c), and compound 2 (d) measured on human hepatocarcinoma HepG2 cells after 24 and 72 h of exposure to different concentrations of test substances (1.0 to 500 µM).Representative images of the HepG2 cell culture after 72 h of treatment with 100, 250, and 500 µM of etoposide (e) and compound 2 (f).The nonlinear regression curve of etoposide (g) and compound 2 (h) shows the inhibition of HepG2 cells after 72 h incubation with the respective test substance (the IC 50 values are indicated).The results are expressed as the mean % of untreated controls ±SD (n = 3).Statistical analysis was performed by one-way ANOVA and Dunnett's multiple comparison test.*, p < 0.1; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001 vs. the untreated control group (Ctrl.).