Synthesis, Characterization, and Biological Evaluation of Meldrum’s Acid Derivatives: Dual Activity and Molecular Docking Study

In the presented study, eight novel Meldrum’s acid derivatives containing various vanillic groups were synthesized. Vanillidene Meldrum’s acid compounds were tested against different cancer cell lines and microbes. Out of nine, three showed very good biological activity against E. coli, and HeLa and A549 cell lines. It is shown that the O-alkyl substituted derivatives possessed better antimicrobial and anticancer activities in comparison with the O-acyl ones. The decyl substituted molecule (3i) has the highest activity against E. coli (MIC = 12.4 μM) and cancer cell lines (HeLa, A549, and LS174 = 15.7, 21.8, and 30.5 μM, respectively). The selectivity index of 3i is 4.8 (HeLa). The molecular docking study indicates that compound 3i showed good binding affinity to DNA, E. coli Gyrase B, and topoisomerase II beta. The covalent docking showed that 3i was a Michael acceptor for the nucleophiles Lys and Ser. The best Eb was noted for the topoisomerase II beta-LYS482-3i cluster.


Introduction
Meldrum's acid presents a compound that has an active methylene group, which was discovered by Meldrum [1]. Therefore, Meldrum proposed the wrong chemical structure. Four decades later, the structure of Meldrum's acid was solved by Davidson [2]. Its carbonyl function easily can be attacked predominately by oxygen or nitrogen nucleophiles. The C5 position is involved in electrophilic substitution. At high-temperature regimes, Meldrum's acid and its derivatives could suffer from a fragmentation process that leads to reactive ketene formation [3]. Although it was discovered more than a century ago, its use continues to attract the attention of synthetic chemists. Electrophilicity, dienophilicity, and its relatively low acidity (pKa = 4.9) make Meldrum's acid attractive and quite useful for a wide number of synthetic manipulations such as heterocycles synthesis, Friedel-Crafts acylation, total synthesis, alkylidene and asymmetric synthesis, domino reactions, and

Results
Our initial idea was to synthesize a set of benzylidene Meldrum's acid derivatives starting from different vanillic aldehydes that have long-chain fragments. The first reaction was tested with 5 -nitrovanilline (2a) and Meldrum's acid (1). Applying similar catalytic conditions such as those described in our published paper [9], we isolated product 3a (79%). In the presence of the complex salt (PhNH 3 ) 2 CuCl 4 as the catalyst in a concentration of 1 mg/mL (0.85 mol%), the reaction was very fast (5 min). Therefore, PTSA (p-toluene sulfonic acid) as catalyst (5 and 10 mol%) produced 3a in lower yield (62%) during 12 h. The same experimental conditions were tested in reactions 2b and 2f with 1. In both cases PTSA gave lower yield in comparison with (PhNH 3 ) 2 CuCl 4 . In some cases, yields were significantly lower. For example, employing PTSA the yield of 3f was 51%, and 75% when (PhNH 3 ) 2 CuCl 4 was loaded. The results prompted us to apply (PhNH 3 ) 2 CuCl 4 (1 mg/mL) as the catalyst in reaction 1 and various O-acyl or -alkyl vanillic aldehydes. Applying the catalytic system, good to excellent yields (54-92%) of vanillidene Meldrum's acid derivatives were obtained (Scheme 1). Almost all products were isolated very easily by filtration. However, there were some problems during the workup of compound 3e. The increased number of oxygen atoms in the side chain made 3e very soluble in ethanol. The reaction mixture was evaporated, and upon the addition of the DCM product 3e was precipitated, filtered, and washed with the DCM. Acyl groups (propionyl, cyclopropanecarbonyl, and p-metoxybenzoyl) as structural motifs in 3b-d (61-69%) negatively affected the yield in comparison with long-chain alkyl groups in 3f-i (75-92%).
Proton and carbon NMR spectra of 3a-i samples (cca. 40 mg/mL) were carried out in CDCl 3 or DMSO-d 6 (Figures S1-S18). In every sample, the methyliden proton was the best measure of condensation efficiency. It can be found at approximately 8.35 ppm in proton NMR. In the 1 H NMR spectrum of 3a ( Figure S1), 1.74 (s), 3.92 (s), 8.14 (s), 8.33 (s), and 8.54 (s) ppm are assigned to 2,2-dimethyl, methoxy, aromatic, and methyliden protons, respectively. The phenolic proton in 3a is invisible in its NMR, which is presumably caused by the strong electron-withdrawing effect of the nitro group attached to the ortho position. The specific carbon peak located at 155 ppm ( Figure S2) was assigned to phenolic Pharmaceuticals 2023, 16, 281 3 of 13 carbon, fragment C-OH. Moreover, in all samples, the high-intensity peak at 1.7-1.8 ppm is assigned to the 2,2-dimethyl backbone of Meldrum's acid fragment. As compared to all carbon NMR, the 13  Proton and carbon NMR spectra of 3a-i samples (cca. 40 mg/mL) were carried out in CDCl3 or DMSO-d6 (Figures S1-S18). In every sample, the methyliden proton was the best measure of condensation efficiency. It can be found at approximately 8.35 ppm in proton NMR. In the 1 H NMR spectrum of 3a ( Figure S1), 1.74 (s), 3.92 (s), 8.14 (s), 8.33 (s), and 8.54 (s) ppm are assigned to 2,2-dimethyl, methoxy, aromatic, and methyliden protons, respectively. The phenolic proton in 3a is invisible in its NMR, which is presumably caused by the strong electron-withdrawing effect of the nitro group attached to the ortho position. The specific carbon peak located at 155 ppm ( Figure S2) was assigned to phenolic carbon, fragment C-OH. Moreover, in all samples, the high-intensity peak at 1.7-1.8 ppm is assigned to the 2,2-dimethyl backbone of Meldrum's acid fragment. As compared to all carbon NMR, the 13 C NMR spectrum of 3b-d ( Figures S4, S6 and S8) showed new peaks at around 171 ppm that originated from the acyl carbonyl function.
The strong and broad FT-IR peak at around 3204 cm −1 is assigned as an overlapped peak of phenolic OH and methyliden function. The band at 2949 cm −1 is related to C-H stretching of a chemically non-equivalent double bond functional within 3a. Two strong bands at 1751 and 1703 cm −1 denote the vibration of a carbonyl group from Meldrum's acid skeleton.
Vanillidene Meldrum's acid derivatives 3a-i were subjected to biological evaluation. For these experiments, we treated various microbes (Escherichia coli, Bacillus subtilis, Staphylococcus aureus, and Bacillus cereus) and cancer cell lines (HeLa, K562, A549, LS174, and PaCa-2). Selectivity was checked towards normal human fibroblast (MRC-5). Streptomycin and cis-platinum (cis-DDP) were chosen as reference standards. The data presented in Tables 1 and 2 are average antimicrobial and anticancer activities of 3a-i delivered after three independent experiments. The strong and broad FT-IR peak at around 3204 cm −1 is assigned as an overlapped peak of phenolic OH and methyliden function. The band at 2949 cm −1 is related to C-H stretching of a chemically non-equivalent double bond functional within 3a. Two strong bands at 1751 and 1703 cm −1 denote the vibration of a carbonyl group from Meldrum's acid skeleton.
Vanillidene Meldrum's acid derivatives 3a-i were subjected to biological evaluation. For these experiments, we treated various microbes (Escherichia coli, Bacillus subtilis, Staphylococcus aureus, and Bacillus cereus) and cancer cell lines (HeLa, K562, A549, LS174, and PaCa-2). Selectivity was checked towards normal human fibroblast (MRC-5). Streptomycin and cis-platinum (cis-DDP) were chosen as reference standards. The data presented in Tables 1 and 2 are average antimicrobial and anticancer activities of 3a-i delivered after three independent experiments. Going deeper inside interactions that should be responsible for antimicrobial or anticancer activity, molecular docking was applied to the most active compounds (3e, 3h, and 3i). The mentioned compounds were docked to DNA, Escherichia coli Gyrase B, and topoisomerase II beta. The energies of binding (E b ) of 3e, 3h, and 3i are presented in Table 3. The E b of reference compounds ellipticine (for DNA), P3C (for Gyrase B), and etoposide (for topoisomerase II) are −8.98, −7.9, and −14.04 kcal mol −1 , respectively. Considering these results, the most promising molecule that showed good binding affinity towards reference compounds is 3i.

Discussion
The results of antimicrobial testing of 3a-i against selected microbes showed significant activity (Table 1) against Gram-positive (B. subtilis, S. aureus, and B. cereus) and Gramnegative bacteria (E. coli). All the tested compounds showed significantly better activity against Gram-negative in comparison with Gram-positive strains. From a structural point of view, vanillidene with the O-alkyl fragment (3e-i) showed to be much more active in contrast to ones with O-acyl compounds (3b-d) against all treated bacteria. The principal reason could be found in the cell wall structural difference. The Gram-positive bacteria cell wall contains peptidoglycans and teichoic acids, while Gram-negative bacteria have lipopolysaccharides and lipopoliproteins as building units of the cell wall [25]. Hence, the presence of a long alkyl chain in structures 3e-i make them more lipophilic, and consequently enables easier passage of these molecules through the cell wall. Out of all nine, the compound with the decyl fragment (3i) has the most promising activity against E. coli (12.4 µM). Compound 3a has the lowest activity against all bacterial strains.
The cytotoxic potential of 3a-i was studied against various cancer and normal cell lines. The results are presented as IC 50 values in Table 2. Several compounds have good activity, while others possessed a moderate to weak cytotoxic effect. Generally, all nine molecules demonstrated the best activity against A549 cell lines in the range of 21.8-45.7 µM. The best activity was shown by 3i. In particular, this molecule stands out against HeLa, A549, and LS174 cell lines, with values of 15.7, 21.8, and 30.5, respectively. The achieved IC 50 values are similar to the ones obtained after treatment of HeLa, A549, and LS174 with cis-DDP. The investigation of MRC-5 cell lines suggests no significant toxicity of 3a-g, while 3h and 3i showed moderate toxicity. The three most active compounds possess good selectivity indexes (SI) for HeLa cells (3h, 3e, and 3i, with values of 2.9, 4.0, and 4.8, respectively). Significant SIs were obtained for A549 cells at 4.3, 2.6, and 3.6 for 3h, 3e, and 3i, respectively. Acyl fragments in structures of 3b-d influenced the obtained IC 50 value being at least double compared to ones delivered to alkyl substituted vanillidene Meldrum's acid (3e-i). The compounds of interest are 3e, 3h, and especially 3i (Figure 1), which has the highest cytotoxicity and antimicrobial activity against all treated cell lines and bacteria. Considering all results noted during antimicrobial and cytotoxicity testing, we employed molecular docking to investigate the affinity of 3e, 3h, and 3i to DNA, E. coli Gyrase B, and topoisomerase II beta at the molecular level.
The selected compounds (3e, 3h, and 3i,) were computationally evaluated for intercalating ability by docking to a six base pair DNA structure (PDB: 1z3f) against ellipticine, a co-crystalized substrate from the used file [26]. The results are presented in Table 3. Testing of each compound produced results of higher energy when compared against the reference substrate, indicating weaker binding. In each case, the pseudo planar part of the resulting lowest energy structure was packed into an intercalating cavity with slightly different positioning, indicating that each compound may bind in this way. The compound with the longest aliphatic chain (3i) showed the most promising results and is displayed in Figure 2.
Pharmaceuticals 2023, 16, x FOR PEER REVIEW 5 of 13 and 4.8, respectively). Significant SIs were obtained for A549 cells at 4.3, 2.6, and 3.6 for 3h, 3e, and 3i, respectively. Acyl fragments in structures of 3b-d influenced the obtained IC50 value being at least double compared to ones delivered to alkyl substituted vanillidene Meldrum's acid (3e-i). The compounds of interest are 3e, 3h, and especially 3i (Figure 1), which has the highest cytotoxicity and antimicrobial activity against all treated cell lines and bacteria. Considering all results noted during antimicrobial and cytotoxicity testing, we employed molecular docking to investigate the affinity of 3e, 3h, and 3i to DNA, E. coli Gyrase B, and topoisomerase II beta at the molecular level. The selected compounds (3e, 3h, and 3i,) were computationally evaluated for intercalating ability by docking to a six base pair DNA structure (PDB: 1z3f) against ellipticine, a co-crystalized substrate from the used file [26]. The results are presented in Table 3. Testing of each compound produced results of higher energy when compared against the reference substrate, indicating weaker binding. In each case, the pseudo planar part of the resulting lowest energy structure was packed into an intercalating cavity with slightly different positioning, indicating that each compound may bind in this way. The compound with the longest aliphatic chain (3i) showed the most promising results and is displayed in Figure 2. It was decided based on our biological tests to computationally examine new compounds for their binding affinity towards enzyme Gyrase B of E. Coli. Molecular docking of our compounds was performed on the crystal structure of Gyrase B (PDB:6f86) where the target cavity was the binding place of 4-(4-bromo-1H-pyrazol-1-yl)-6-[(ethylcarbamoyl)amino]-N-(pyridine-3-yl)pyridine-3-carboxamide (P3C). Since the inhibitory activity of P3C is known (IC50 = 0.037 μM [27]), it serves as a good reference for comparison of obtained docking results for new inhibitors. Binding energies are presented in Table 3. Based on the obtained binding energies, we can conclude that most of our compounds would generally show weaker inhibitory activity when compared to P3C, but the result obtained from docking 3i is comparable to the result displayed by the reference compound. The position of the docked structure for 3i is similar to the position of PC3, and E. coli Gyrase B, and topoisomerase II beta at the molecular level. The selected compounds (3e, 3h, and 3i,) were computationally evaluated for intercalating ability by docking to a six base pair DNA structure (PDB: 1z3f) against ellipticine, a co-crystalized substrate from the used file [26]. The results are presented in Table 3. Testing of each compound produced results of higher energy when compared against the reference substrate, indicating weaker binding. In each case, the pseudo planar part of the resulting lowest energy structure was packed into an intercalating cavity with slightly different positioning, indicating that each compound may bind in this way. The compound with the longest aliphatic chain (3i) showed the most promising results and is displayed in Figure 2. It was decided based on our biological tests to computationally examine new compounds for their binding affinity towards enzyme Gyrase B of E. Coli. Molecular docking of our compounds was performed on the crystal structure of Gyrase B (PDB:6f86) where the target cavity was the binding place of 4-(4-bromo-1H-pyrazol-1-yl)-6-[(ethylcarbamoyl)amino]-N-(pyridine-3-yl)pyridine-3-carboxamide (P3C). Since the inhibitory activity of P3C is known (IC50 = 0.037 μM [27]), it serves as a good reference for comparison of obtained docking results for new inhibitors. Binding energies are presented in Table 3. Based on the obtained binding energies, we can conclude that most of our compounds would generally show weaker inhibitory activity when compared to P3C, but the result obtained from docking 3i is comparable to the result displayed by the reference compound. The position of the docked structure for 3i is similar to the position of PC3, and  Table 3. Based on the obtained binding energies, we can conclude that most of our compounds would generally show weaker inhibitory activity when compared to P3C, but the result obtained from docking 3i is comparable to the result displayed by the reference compound. The position of the docked structure for 3i is similar to the position of PC3, and this similarity is evident from Figure 3 and Figure S19. The difference is, for example, that 3i does not possess a side group bound to a central aromatic ring like the bromopyrazolyl moiety of P3C, which allows interactions with and in the vicinity of ILE94, while 3i has a long aliphatic chain that extends much deeper into the cavity of the enzyme than the ethylcarbamoyl group of P3C.
To further examine the potential anticancer application of examined compounds, molecular docking experiments were performed on the crystal structure of type II topoisomerase beta (TOP2β) in a complex with DNA and the anticancer drug etoposide (PDB: 3qx3) [28]. Benzylidene Meldrum's acid was noted as the key precursor to novel topoisomerase II inhibitors [29]. One of them is α-lapachone, which is approved as a topoisomerase II inhibitor [30,31]. Redocking of etoposide was used to obtain the reference value of binding energy. The main reason for choosing the selected binding site was comparison with the selected substrate (etoposide) that binds in that exact site. Topoisomerase II is symmetric, and there are two identical sides. Docking experiments were conducted on both sides and identical results (geometry and energy of binding) were obtained. For better clarity we present only one of the identical binding modes. this similarity is evident from Figures 3 and S19. The difference is, for example, that 3i does not possess a side group bound to a central aromatic ring like the bromopyrazolyl moiety of P3C, which allows interactions with and in the vicinity of ILE94, while 3i has a long aliphatic chain that extends much deeper into the cavity of the enzyme than the ethylcarbamoyl group of P3C. To further examine the potential anticancer application of examined compounds, molecular docking experiments were performed on the crystal structure of type II topoisomerase beta (TOP2β) in a complex with DNA and the anticancer drug etoposide (PDB: 3qx3) [28]. Benzylidene Meldrum's acid was noted as the key precursor to novel topoisomerase II inhibitors [29]. One of them is α-lapachone, which is approved as a topоisomerase II inhibitor [30,31]. Redocking of etoposide was used to obtain the reference value of binding energy. The main reason for choosing the selected binding site was comparison with the selected substrate (etoposide) that binds in that exact site. Topoisomerase II is symmetric, and there are two identical sides. Docking experiments were conducted on both sides and identical results (geometry and energy of binding) were obtained. For better clarity we present only one of the identical binding modes.
The results derived from docking of the examined compounds are compared against the reference value in Table 3. Significantly higher binding energies of tested compounds when compared to etoposide suggest much weaker binding in the same pocket. Upon investigation of structures obtained from the best docking results, it was found that contrary to etoposide, which intercalates with DNA and interacts with the amino acid residues of the enzyme that can be found on both sides of the DNA chain, our compounds position themselves in such a way that only allows interactions with residues from a single side of the chain (Figures 4 and S20). The results derived from docking of the examined compounds are compared against the reference value in Table 3. Significantly higher binding energies of tested compounds when compared to etoposide suggest much weaker binding in the same pocket. Upon investigation of structures obtained from the best docking results, it was found that contrary to etoposide, which intercalates with DNA and interacts with the amino acid residues of the enzyme that can be found on both sides of the DNA chain, our compounds position themselves in such a way that only allows interactions with residues from a single side of the chain (Figure 4 and Figure S20).  Docking of covalently bound ligands is important for gaining insight into enzymatic processes and designing superior covalent ligands. AutoDock 4 was used to implement two methods for docking covalently bound ligands: a grid-based method and a modification of the flexible side chain approach [32]. Covalent docking has also been applied to oligopeptidase [33] and proteasome [34] inhibitors and glycoside hydrolyses ligands [35]. Benzylidene Meldrum's acid was used as a valuable Michael acceptor [36,37]. Consequently, in this paper, covalent docking was employed to simulate the formation of the covalent adduct using the flexible side chain technique proposed by Bianco et al. [32]. This protocol needs to adapt the residue taking part in the covalent bond by attaching the ligand to its side chain; this modified residue is then considered flexible during the docking calculation. Docking was attained for the ligand, keeping the remodeled cysteine and serine residues flexible. This permitted us to sample the torsional flexibility of the ligand within the Gyrase B and topoisomerase II beta. It was shown that the best conformational compatibility between the subjected ligand Michael acceptor (3i) and substrate Gyrase B or topoisomerase II beta had residues LYS139 (Gyrase B), LYS482, and SER480 (topoisomerase II beta). The molecular interactions between 3i and the amino acids lysine and serine in Gyrase B and topoisomerase II beta were estimated by using covalent docking Docking of covalently bound ligands is important for gaining insight into enzymatic processes and designing superior covalent ligands. AutoDock 4 was used to implement two methods for docking covalently bound ligands: a grid-based method and a modification of the flexible side chain approach [32]. Covalent docking has also been applied to oligopeptidase [33] and proteasome [34] inhibitors and glycoside hydrolyses ligands [35]. Benzylidene Meldrum's acid was used as a valuable Michael acceptor [36,37]. Consequently, in this paper, covalent docking was employed to simulate the formation of the covalent adduct using the flexible side chain technique proposed by Bianco et al. [32]. This protocol needs to adapt the residue taking part in the covalent bond by attaching the ligand to its side chain; this modified residue is then considered flexible during the docking calculation. Docking was attained for the ligand, keeping the remodeled cysteine and serine residues flexible. This permitted us to sample the torsional flexibility of the ligand within the Gyrase B and topoisomerase II beta. It was shown that the best conformational compatibility between the subjected ligand Michael acceptor (3i) and substrate Gyrase B or topoisomerase II beta had residues LYS139 (Gyrase B), LYS482, and SER480 (topoisomerase II beta). The molecular interactions between 3i and the amino acids lysine and serine in Gyrase B and topoisomerase II beta were estimated by using covalent docking calculations. The obtained results for E b are presented in Table 4. The more negative E b values indicate that the investigated ligand 3i inhibits the receptor better. As shown in Table 4, the investigated ligand strongly binds to the target receptors. The docking analyses revealed that covalent interactions existed between the investigated molecule and target receptors. These interactions occur between 3i and the amino acids LYS and SER in positions 482 and 480 in the primary structure of the topoisomerase II beta-receptor. On the other hand, 3i forms the covalent bond with the amino acid LYS139 from Gyrase B ( Figure 5). Additionally, the docking results show that several non-covalent interactions occurred between the investigated molecule and target receptors. The important interactions are hydrogen bonds, carbon-hydrogen bonds, alkyl-alkyl, and alkyl-π interactions ( Figure 5). The obtained results indicate that the 3i ligand could act as a potential covalent inhibitor. The α-and β-lapachones were described as covalent inhibitors of topoisomerase II [38]. The mentioned compounds have also been investigated as irreversible catalytic inhibitors of topoisomerase II [30]. The high affinity of 3i to form a covalent bond with topoisomerase II beta makes it similar to the aforementioned Meldrum's acid analogs.

Materials and Methods
The IR spectra were recorded by a Perkin-Elmer Spectrum One FT-IR spectrometer on a KBr pellet. The NMR spectra of compounds 3a-i were performed in DMSO-d6 and The αand β-lapachones were described as covalent inhibitors of topoisomerase II [38]. The mentioned compounds have also been investigated as irreversible catalytic inhibitors of topoisomerase II [30]. The high affinity of 3i to form a covalent bond with topoisomerase II beta makes it similar to the aforementioned Meldrum's acid analogs.

Materials and Methods
The IR spectra were recorded by a Perkin-Elmer Spectrum One FT-IR spectrometer on a KBr pellet. The NMR spectra of compounds 3a-i were performed in DMSO-d 6 and CDCl 3 with TMS as the internal standard on a Varian Gemini 200 MHz NMR spectrometer ( 1 H at 200 and 13 C at 50 MHz). Abbreviations for the NMR signals that were used are as follows: s = singlet, d = doublet, dd = double of double, t = triplet, and m = multiplet. The 1 H and 13 C spectra are given in Supplementary Information (Figures S1-S18).

Synthesis of Benzylidene Meldrum's Acid Derivatives (3)
In the 10 mL round-bottomed flask, Meldrum's acid (1.5 mmol) was dissolved in 5 mL of absolute ethanol. Selected aldehyde was added at room temperature. Immediately after, (PhNH 3 ) 2 CuCl 4 (1 mg/mL, 0.85 mol%) was loaded as catalyst. Precipitation of products occurred in the interval of 5 to 10 min. One hour later, amorphous powder was filtrated and washed out with 96% ethanol and water. Knoevenagel adducts were isolated in good purity grade for NMR measurements. NMR data for 3a-i are given below.

Determination of Antimicrobial Activity
The following strains of bacteria were used as test organisms in this study: Staphylococcus aureus (ATCC 25923), Bacillus subtilis (ATCC 6633), Bacillus cereus (ATCC 10987), and Escherichia coli (ATCC 25922). All the bacteria used were obtained from the American Type Culture Collection (ATCC). Cultures were stored at 4 • C and subcultured every 15 days. The minimal inhibitory concentration (MIC) was determined by the broth microdilution method using 96-well micro-titer plates [39]. A series of dilutions of the tested compounds was used in the experiment against every microorganism. The starting solutions of tested compounds were obtained by measuring off a certain quantity of the compounds and dissolving it in 5% DMSO. The MIC was determined with resazurin. The inoculated plates were incubated at 37 • C for 24 h. Resazurin is an oxidation-reduction indicator used for the evaluation of microbial growth. It is a blue non-fluorescent dye that becomes pink and fluorescent when reduced to resorufin by oxidoreductases within viable cells. The boundary dilution without any changing in color of resazurin was defined as the MIC for the tested microorganism at a given concentration. As a positive control of growth inhibition, streptomycin was used. A 5% DMSO solution was used as a negative control for the influence of the solvents.

Evaluation of Cytotoxicity
Cell line cells were obtained from the American Type Culture Collection. Human cervical adenocarcinoma HeLa, human chronic myelogenous leukemia K562, non-small cell lung carcinoma A549, human colon carcinoma LS174, pancreatic carcinoma PaCa-2, and normal human lung fibroblast MRC-5 cell lines were grown in RPMI-1640 medium with 10% fetal bovine serum, L-glutamine, and penicillin-streptomycin solution. Cells were plated into the 96-well cell culture plates. Cells were incubated at 37 • C in CO 2 incubator. Adherent cell lines were incubated for 20 h before addition of tested compounds. Two hours before addition of compounds, K562 cells were seeded into 96-well plates. The tested compounds were applied at five different concentrations ranging from 10 to 300 µM for 72 h. Thereafter, the survival of cells was determined by MTT assay in accordance with the protocol established by Mosmann [40] and Ohno and Abe [41]. The solution of MTT was added to blank and cell samples. The plates were incubated for 5 h at 37 • C, and then 10% solution of sodium dodecyl sulfate was added to the wells. The absorbance was measured at 570 nm using Thermo Scientific Multiskan EX plate reader. Three independent experiments were performed.

Molecular Docking
Molecular docking experiments were prepared using AutoDockTools and were performed using AutoDock 4 [42]. Structures of target molecules were obtained from rscb.org, and were prepared by first removing co-crystalized substrates, ions, and water molecules, followed by calculating Gasteiger charges and removing non-polar hydrogens. Structures of investigated molecules were optimized using PM7 method of MOPAC2016 [43], but Gasteiger charges were calculated for use in docking experiments. Each co-crystalized substrate was re-docked to its original parent molecule to obtain reference value of binding energy and to validate the applied method. Grid maps were calculated using cube-shaped grid boxes 60 pts wide for protein targets and 40 pts wide for DNA target (1 pt = 0.375 A). Grid boxes were centered using coordinates of co-crystalized substrates. Each experiment consisted of 10 hybrid genetic algorithm-local search runs, with 2.5 × 10 7 energy evaluations per run.
The two-point attractor and flexible side chain approaches [32] are two types of covalent docking methods that can be used with AutoDock 4. Both methods use precalculated grid maps and atom probes to speed up the scoring process. However, they simulate ligand conformations in different ways during the scoring process. In this paper, we used the flexible side chain method. This approach uses tethered docking to mimic the way covalent ligands are bound in the pocket. For this plan to work, the electrophilic center of the ligand must bond to the two nucleophilic atoms at the ends of the protein. By figuring out the right SMARTS pattern, these two atoms are put on top of the right residue atoms in the protein to make the structure that is wanted. The bound ligand is then treated as a flexible residue, and the standard AutoDock 4 method for flexible residues is used to look at its different shapes in the pocket. With the aid of the scripts provided by the AutoDock 4 website [44], the ei ligand and cysteine and serine residues were overlapped. Subsequently, the receptor grid maps were calculated with the AutoGrid4 software, mapping the receptor interaction energies using the ligand atom types as probes. The grid of 60 Å × 60 Å × 60 Å with 0.375 Å spacing was centered on the coordinates of the ligand originally present in the Gyrase B and topoisomerase II beta. Finally, the actual docking was attained for the ei ligand, keeping the cysteine and serine residues as flexible. This permitted the sampling of the torsional flexibility within the receptors.

Conclusions
In this paper, eight novel Meldrum's acid derivatives containing vanillic fragments were synthesized under soft reaction conditions in good to excellent yields (54-92%). All compounds showed good to moderate anticancer and antimicrobial activities. Vanillidene Meldrum's acid derivatives with an O-acyl group attached on the vanillic motif have lower activity against cancer cell lines and bacteria. However, the presence of the O-alkyl fragment significantly influenced the biological activity, both anticancer and antimicrobial. The most active compound that showed dual activity was 3i. This molecule contains a decyl chain that is probably responsible for the good activity against E. coli (12.4 µM) and cancer cell lines (IC 50 for HeLa, LS174, and A549: 15.7, 20.8, and 21.8 µM, respectively). Significant dual activity was displayed by compounds 3e and 3h. Very good selectivity indices were accomplished with 3e, 3h and 3i. The highest selectivity index was for 3i against HeLa (4.8), while the next highest was 3h against A549 cell lines (4.3). In further investigation, the most active compounds (3e, 3h, and 3i) were subjected to molecular docking. The binding affinities of mentioned compounds to DNA, E. coli Gyrase B, and topoisomerase II beta were tested. Compounds 3e, 3h, and 3i showed good non-covalent interactions and consequently good energy of binding. Therefore, 3i possessed the highest affinity and was subjected to a covalent docking study. Under this investigation, 3i was played as the Michael acceptor in the reaction with the amino or hydroxy group from Lys and Ser residues, respectively. The lowest binding energy was realized for the topoisomerase II beta-LYS482-3i cluster (E b = −7.40 kcal mol −1 ). From molecular docking and experimental outputs, it was shown that 3i is a molecule of interest. Going forward, we believe that the presented concept and dual-active compounds described herein have bright futures.