Green synthesis, structural analysis and anticancer activity of dihydropyrimidinone derivatives

In this study, for the first time, we have used Citrus macroptera juice to synthesize dihydropyrimidine (DHPM) derivatives via the Biginelli reaction, which showed better yield, shorter reaction time, and did not require an organic solvent for the reaction. A series of DHPM derivatives were synthesized, and characterized, and structural analysis was achieved through SCXRD & Hirshfeld surface analysis. We observed that these synthesized dihydropyrimidine (DHPM) derivatives showed C–H⋯π, C–H⋯O, C–H⋯N, C–H⋯C, lone pair⋯π, π⋯π, etc. interactions. We also performed in silico studies for their inhibitory activities against human kinesin Eg5 enzyme, and the cytotoxic activity of the synthesized compounds was carried out against A549 lung adenocarcinoma cells. In silico analysis demonstrated that compounds with a chloro-group at the 3- or 4-position in the substituted ring of DHPM showed higher binding affinity for the human kinesin Eg5 enzyme (−7.9 kcal mol−1) than the standard drug monastrol (−7.8 kcal mol−1). Furthermore, in vitro cellular studies revealed that compounds with a chloro-group at the 3- or 4-position in the substituted ring of DHPM induced significant cell death in human A549 lung adenocarcinoma cells. This result indicates that a deactivating group (chlorine) at the 3- or 4-position in the substituted ring of DHPM might be a promising anticancer drug candidate for treating different types of cancers, particularly cancer of the lung.


Introduction
Cancer accounts for the second leading cause of death worldwide despite many breakthroughs for cancer treatment. Almost 10.0 million cancer deaths occurred due to cancer in 2020 alone. 1 According to the WHO, if this trend continues, this is expected to rise to 13.1 million by the end of 2030. Due to drug resistance and toxicity associated with its treatment, cancer remains a focus of immense interest. 2,3 This warrants newer drug discovery and drug development efforts to mitigate these concerns in cancer therapy.
Heterocyclic molecules continue to be attractive targets for synthesis since they oen exhibit vital biological properties. In particular, dihydropyrimidinones (DHPMs) are well known for their wide range of bioactivities. DHPMs have vital applications in the eld of drug discovery. Three nucleic acid bases are pyrimidine derivatives, which comprise cytosine, found in DNA and RNA, uracil in RNA, and thymine in DNA (Fig. 1). Therefore, the synthesis and study of DHPM systems with methylene-linked ester makes it an exciting topic for synthetic organic chemistry.
DHPM classes of compounds are structural analogs of monastrol. These are important to medicinal chemistry because of many biological activities. 4 These compounds also possess interesting biological activities such as antiviral, 5 antitumor, 6 and antibacterial, 7 as well as calcium channel modulating activity. 8 Monastrol is the protagonist of the DHPMs class. Several studies revealed that DHPMs inhibitory action on human kinesin Eg5 leads to a mitotic arrest resulting in apoptosis. 9 This was the main action in apoptosis described for The synthesis of DHPM and evaluation of their anticancer properties is of ongoing interest to our laboratory. Synthesis of DHPM derivatives by "solvent-free and catalyst-free" green methodology has been done by others. 13 Pramanik et al. developed a safe, eco-friendly, and green method for synthesizing DHPM at room temperature using common fruit juice. 14 Microwave-assisted green synthesis of DHPM was developed by Bhatewara et al. and studied their antimicrobial properties. 15 Citrus macroptera (a fruit from citrus family, usually found in India, Bangladesh, Malaysia and Melanesia) juice medium never have been used for Biginelli reaction using electron-rich/ electron-decient aromatic aldehydes, but similar kind of citrus fruit like orange and lime has been reported earlier for similar reactions. 14 Here, we have reported a green synthesis for DHPM with electron-rich and electron-decient aromatic aldehydes in Citrus macroptera juice medium [Scheme 1]. The non-covalent interactions of these compounds were studied by single-crystal X-ray diffractometers (SCXRD). Hirshfeld surface analyses were performed to study the structural property of molecules. DFT calculations were carried out to know how the substituents impact the geometry of different synthesized DHPM derivatives. In silico analysis studies the molecules binding affinity in the active site of target proteins (Eg5). Subsequently, we also carried out cytotoxic studies of these compounds in the lung adenocarcinoma (A549) cells.

Experimental design
1 H NMR (300 MHz) and 13 C (75 MHz) NMR spectra were recorded on JEOL AL300 FTNMR spectrometer using TMS as an internal reference, and chemical shi values are expressed in d, ppm units. Melting points of all the compounds were recorded on the electrically heated instrument and were uncorrected. All the reactions were monitored by thin-layer chromatography (TLC) using an appropriate solvent system. The compound was puried from column chromatography (silica gel size 60-120 mesh) and ash chromatography (silica gel size 230-400 mesh) using the standard protocol.

Chemistry
Our main target was to perform green synthesis of DHPMs containing electron-rich and electron-decient phenyl rings from Citrus macroptera juice and study the non-covalent interactions of these compounds through X-ray crystallographic study. The required Citrus macroptera juices were extracted directly from the naturally obtained Citrus macroptera fruit, and it was used straightaway for the reaction. DHPMs synthesis via Biginelli reaction with electron-rich and electron-decient aromatic aldehydes in Citrus macroptera juice medium was done at room temperature (Scheme 1).
2.1.1. General procedure for the synthesis of DHPM derivatives (1.1-1.5). The equimolar quantities of ethyl acetoacetate (20 mmol, 2.6 g), thiourea (20 mmol, 1.52 g), aromatic aldehyde (20 mmol) were stirred in 4 mL of Citrus macroptera juice at room temperature for 12 h. Upon completion of the reaction, the solid product was precipitated out of the reaction medium. The crude solid product was ltered and recrystallized from hot ethanol to get the pure compounds (1.

X-ray crystallographic study
All crystal data has been used to solve and rene the crystal structure by SHELXL-97 16 in the program suite WinGX (version 1.80.05). 17 The molecular graphics were developed by using ORTEP-3 18 soware. For 3D structure visualization and packing diagrams, Mercury 3.3 soware has been used throughout the experiment. PARST-95 and PLATON 19 were used for geometrical calculations. Hydrogen bonds were represented by broken light green lines in the packing diagram. Carbon atoms were represented by grey color, hydrogen atoms white, oxygen atoms red, nitrogen atoms light blue, sulphur atoms yellow, chlorine atoms blue, etc., in the packing diagram. Single-crystal X-ray data for compounds 1.3, 1.4, and 1.5 were collected with an Oxford Diffraction Xcalibur CCD diffractometer.

Hirshfeld surface analysis
Hirshfeld surface analysis 20 is used to visualize intermolecular interactions in crystals using 3D molecular surface contours and 2D ngerprint plots. The Hirshfeld surface and associated 2D ngerprint plots for the compounds were calculated using Crystal Explorer 17. Differing van der Waals radii determined the distance between the points on the surface to a nucleus (atom) inside (di) and outside (de) the mean surfaces. The Hirshfeld surface, ngerprint plots, and interaction energies of compounds have been calculated using Crystal Explorer 17. 21 The C-H bond lengths converted to normalized values. This is based on neutron diffraction results. 22 The interaction energies of compounds have been calculated using the CE-B3LYP/6-31G(d,p) functional/basis set combination. Superposition energy (BSSE) was corrected for the basis set using the counterpoise (CP) method. 23

Density functional theory (DFT) calculations
Prior to investigating the effect of substituents on the synthesized DHPM derivatives on their non-covalent interactions and anticancer drug property, it is of importance to understand how the substituents impact the geometry of different synthesized DHPM derivatives (1.3, 1.4 and 1.5). We, therefore, carried out a series of density functional theory (DFT) calculations with the hybrid meta-GGA functional M06-2X 24 and Pople's basis set 6-31+G(d,p), as implemented in Gaussian 09 soware suite, 25 to evaluate the stability of different possible conformers formed due to twisting of the phenyl rings with respect to the nonaromatic dihydropyrimidinone ring in the synthesized DHPM derivatives. The M06-2X functional developed by Truhlar's group is reported to be a reliable functional for the study of relative stability of conformers, noncovalent interactions, rate constants and for thermochemistry of a wide variety of molecular systems. 26 2.5.1. Cell culture. Human adenocarcinoma A549 cell line was acquired from ATCC (Manassas, VA, USA) and grown in RPMI-1640 media complemented with 1% PSA and 10% FBS at 37 C and 5% CO 2 . The medium was changed every alternate day, and cells were trypsinized aer achieving the conuency.
2.5.2. MTT assay. The cytotoxicity or viability study of compounds 1.1, 1.2, 1.3, 1.4, and 1.5 was done with A549 cells. Approximately 8000 cells were plated per well in a 96-well plate, and aer attachment, cells were treated with different concentrations of compounds 1.1, 1.2, 1.3, 1.4, and 1.5. Media were removed from wells aer 24 and 48 h, followed by incubation with 50 ml of MTT solution in PBS (5 mg mL À1 ) for 4 h. Formazan crystals were dissolved in DMSO (150 ml), and absorbance was taken at 570 nm in a multimode plate reader (Synergy H1 Hybrid Reader, Biotek; Winooski, VT, USA).
2.5.3. Trypan blue assay. Trypan blue is an azo dye that stains dead cells with the compromised membrane in blue colour, and it does not stain live cells with the intact membrane, thus they appear clear. Nearly 60 000 cells were plated in 12-well plates, subjected to various doses of 1.2, 1.4, and 1.5 for 24 and 48 h. The cells were harvested and pelleted at 1200 rpm for 5 min. The cells were mixed in 200 ml 1Â PBS aer removing the supernatant, and 0.5% trypan blue dye was added. Cells were counted under a phase-contrast microscope (Zeiss, Germany).
2.5.4. Acridine orange-ethidium bromide staining. To assess apoptosis and necrosis in A549 cells, nearly 60 000 cells were seeded in a 12-well plate followed by treatment with selected compounds 1.2 and 1.5 (200 mM) for 24 and 48 h time points. Aerward, cells were harvested and pelleted at 1200 rpm for 5 min. The cell mix was prepared in 50 ml of 1Â PBS aer removing the supernatant, and to 10 ml of cell aliquot 1 ml acridine orange (1 mg mL À1 ) and 1 ml ethidium bromide (1 mg mL À1 ) were added on a slide and mounted with a cover slip. The slides were subjected to uorescence microscopy (Olympus, Japan), and photomicrographs were taken at 200Â.
2.5.5. Hoechst 33342 staining. Hoechst 33342 is a nucleic acid staining dye that is used to observe morphological changes in apoptotic cells. Hoechst 33342 assay was performed in a 12well plate in which tissue cultured cover slips were placed. 60 000 cells per well were seeded and incubated at 37 C in a humidied CO 2 incubator and treated with a 200 mM concentration of 1.2 and 1.5 compounds. Aer treatment, media were removed, and cells were rinsed with PBS. 1 mL Hoechst 33342 stain (1 mg mL À1 in 1Â PBS) was added to the wells and incubated in the dark for 15-20 min at 37 C. Hoechst 33342 was removed, and cells were washed thrice with chilled 1Â PBS. Cover slips were removed from the wells and placed upside down on the slide containing mounting media and subjected to uorescent microscopy (Olympus, Japan), and photomicrographs were taken at 200Â.

Molecular docking
Molecular docking analysis was performed using the AutoDock Vina. 28 The crystal structures of the human kinesin Eg5 and survivin protein were retrieved from the RSCB protein data bank (PDB id: 3UIH and 1X88, respectively). 29 The protein preparation was done in autodock tools 30 by removing co-crystallized ligands, embedded water molecules, and cofactors. It was further processed by adding polar hydrogens and saved in the pdbqt format. The crystal structure of ligands in .cif format were converted to .pdb format in Mercury, which is then imported to autodock tools to add gasteiger charges and saved in pdbqt format. The grid parameters for docking with the kinesin Eg5 motor domain were determined based on the native ligand monastrol. However, for survivin protein, the grid parameters were assigned carefully to encompass the binding sites at the dimerization interface and allosteric site located near it. The parameter was set to 9 modes for analyzing binding affinity, and subsequently, crystal compounds 1.1-1.5 were subjected to molecular docking with prepared proteins. The docking results were validated by re-docking. Further, validation of parameters for Eg5 was done by repeating the docking of the crystal structure of monastrol, which was extracted from Eg5 protein, then superimposing with native monastrol in Eg5 protein. The docked monastrol presented an RMSD value of 0.210Å

Results and discussion
A series of DHPM derivatives (1.1-1.5) were synthesized using urea, ethyl acetoacetate with electron-rich and electron-decient aromatic aldehydes through green catalyzed reaction given in Scheme 1. First-time Citrus macroptera juices were used as bio-catalyst for the synthesized DHPM derivatives, which showed better yield and shorter duration. This procedure is also an eco-friendly one, as it does not require any organic solvent for the reaction. Out of ve compounds, three compounds (1.3-1.5) crystallized as single crystals. Cytotoxic activity of synthesized compounds was carried out against A549 lung adenocarcinoma cells to assess their anticancer activity.
3.1. X-ray crystallographic studies of 1.3-1.5 All three prismatic transparent crystals (1.3, 1.4, and 1.5) were obtained using methanol solvent by slow evaporation, aer 15 days, at room temperature. The ORTEP diagrams of 1.3, 1.4, and 1.5 are shown in Fig. 3. In all three crystal structures, the nonaromatic dihydropyrimidinone ring and phenyl ring occupied different planes. Compounds 1.4 and 1.5 crystallized in triclinic space group P 1 with two independent molecules in the unit cell. Compound 1.3 crystallized in monoclinic space group C2/c, having eight molecules in the unit cell. Further details of crystal data and renement are presented in Table 1.
In the crystal structure of compound 1.3, the crystal packing is stabilized by dimers through N-H/S hydrogen bonds in a R 2 2 (8) ring motif (Fig. 4). One more intermolecular hydrogen bonding of N-H/O & C-H/O results in the generation of a dimeric structure, which can be described as R 1 2 (6) graph-set notation in which thioamide & ester oxygen are involved (Fig. 4B). The hydrogen-bonding network for compound 1.3 and crystal packing are given in Fig. 4 Fig. 4). respectively (Table 2 and Fig. 5). An intramolecular C-H/p interactions (H15/Cg ¼ 2.846Å) were also observed, which further stabilized interlayer connectivity of crystal structure.  (Table 2). Both rings are arranged in the AABB pattern in  (8), and R 4 4 (20) graph-set notation (Fig. 6). For crystal structure of 1.5, one p/p stacking interaction was also observed, and distance for p/p stacking interaction was 3.665Å (Fig. 6B). The stacking distances for C-H/p were 3.180, 3.611, 3.480Å (Table 2), which further assisted the interlayer connectivity and linear chain formation of crystal structure. The N-H/O and N-H/S bond distances for hydrogen bonding were 2.1876, 2.584Å, and angles on hydrogen atoms were 138.86 , 168.69 , respectively (Table 2 and Fig. 6). The C-H/O bond distance for R 1 2 (6) ring motif is 2.705 A and angles on hydrogen atoms was 128.86 , while the other C-H/O bond distance were 2.747, 3.456Å and angle on hydrogen atom were 159.68 and 119.28 , respectively ( Table 2).         Table 3. Spoke-like pattern in the ngerprint plots of 1.3, 1.4 and 1.5 represents C-H/O interactions in crystal lattice in region of di + de ¼ 2.00-3.0Å, with H/O/O/H contributions of 12.4%, 6.9% and 7.7%, respectively (Fig. 7A, 8A and 9A). The second spoke-like pattern in the ngerprint plots of 1.3, 1.4 and 1.5 represents C-H/S interactions in crystal lattice in region of di + de ¼ 2.30-3.4Å with H/S/S/H contributions of 15.1%, 15.5% and 16.7%, respectively (Fig. 7A, 8A and 9A). The C-H/p interactions in compounds 1.3, 1.4 and 1.5 can be seen as a pair of unique blue-colored wings in the region of di + de ¼ 3.2-3.6Å with H/C/C/H contribution of 15.2%, 14.2% and 12.9%, respectively (Fig. 7A, 8A and 9A). The C-H/N pair of contacts is also reected as two characteristic wings occupied in di + de ¼ 3.2-3.5Å in compounds 1.3, 1.4, and 1.5. The yellowish-red bin is absent on the ngerprint plots in compounds 1.3 and 1.4, which means the absence of weak p/p stacking in crystal packing ( Fig. 7A and 8A. However, there is a light yellowish-red bin (with C/C contact of contribution is 3.9%) present on the ngerprint plot indicates the existence of weak p/p stacking interactions between the phenyl rings in compound 1.5 (Fig. 9A). The Hirshfeld weak interactions calculation found a similar weak non-covalent intermolecular interactions pattern in crystal packing. C-H/p interactions, C-H/N, and C-H/O interactions, etc., of compound 1.3, 1.4 and 1.5 in crystal packing structure are in Fig. 7B, 8B, 9B) and Table 3. In all three crystal structure's, 1.3, 1.4, and 1.5, extensive hydrogen-bonding network of calculated interactions, terminal carbonyl oxygen, thio-pyrimidine nitrogen, and thio-pyrimidine sulfur are involved in weak interaction and forming eight membered R 2 2 (8) (Fig. 7B, 8B and 9B) rings. The role of p/p stacking is central importance in chemistry as well as bio medicine eld. They are key interactions inuencing the tertiary structure of proteins, the vertical base stacking in DNA, and the intercalation of different drugs into DNA, which helps in drug design. The role of p/p stacking also become prominent in drug-receptor interactions as most of the drugs are aromatic, and about 20% of amino acids are aromatic. The curvedness plots and shape index plots of 3D Hirshfeld is a very promising tool to visualize p/p stacking interactions in the crystals. 34 In curvedness plots, the yellow spots represent very weak intermolecular interactions, and red-yellow colored spots represent strong hydrogen-bonding interactions in the crystal structures of 1.3, 1.4, and 1.5 (Fig. 10). In curvedness plots of compounds 1.3 and 1.4, the absence of green-colored at regions bounded by red-colored rectangles on the phenyl rings further conrms the non-existence of p/p stacking effect in crystal packing (Fig. 10A & C). However, in the curvedness plots of compound 1.5, at green colored regions bounded by red-colored rectangles on phenyl rings, once again conrm the presence of p/p stacking effect in crystal packing of 1.5 (Fig. 10E). Red and blue regions represent the acceptor and the donor property, respectively, in the shape index of compound 1.3, 1.4, and 1.5 (Fig. 10B, D & F). Yellowish-red colored concave regions inside the counters indicate the presence of weak intermolecular interactions in the Shape index plots of compounds 1.3, 1.4, and 1.5 ( Fig. 10B and D).The absence of adjacent red and blue triangles in the shape-index map of 1.3 and 1.4 also proved the lack of p/p stacking effect in the crystal packing (Fig. 10B and D). There are adjacent red and blue triangles bounded by the black colored rectangles in the shape-index map of 1.5, once again conrms the p/p stacking effect in the crystal packing of 1.5 (Fig. 10F). We have also carried out the DFT calculations to know how the substituents impact the geometry of different synthesized DHPM derivatives and calculated the interaction energies of DHPM derivatives (1.3, 1.4 and 1.5) with their neighbouring molecules in the crystal structure at the M06-2X/6-31+G(d,p) level of theory. Torsional PES (potential energy surfaces) for twisting the phenyl rings in synthesized DHPM derivatives (1.3, 1.4 and 1.5) are shown in Fig. 11.
We note from Fig. 11 that each PES is characterized by only one minimum and the fully planar congurations (0 and 180 ) are unstable. The DHPM derivatives with torsion angle ($40 ) are found to be energetically most unstable. The global minimum for 1.3, 1.4 and 1.5 molecule are found to be at torsion angle of 155.5 , 146.8 and 152.3 , respectively (Fig. 11), which is in agreement with their crystal structure. The optimized geometries of 1.3, 1.4 and 1.5 are shown in Fig. S5 (ESI †). Thus, our DFT calculations show that substituents do not impact the geometry of synthesized DHPM derivatives appreciably.
We also calculated the interaction energies of DHPM derivatives (1.3, 1.4 and 1.5) with their neighbouring molecules in the crystal structure at the M06-2X/6-31+G(d,p) level of theory. We considered neighbouring molecules which are in direct contact with a central molecule for complex formation and thus complexes formed by 8, 5 and 7 molecules for 1.3, 1.4 and 1.5, respectively, were taken for calculation of interaction energy. These complexes are shown in Fig. S7 and S8 (ESI †). The interaction energy of a complex was calculated as follows: Interaction energy ¼ total energy of the complex À n Â total energy of one molecule where n is the number of molecules in the complex. The interaction energies of the complexes of 1.3, 1.4 and 1.5 are found to be À103.20, À24.61, À90.81 kcal mol À1 , respectively. The interaction energy per molecule for 1.3, 1.4 and 1.5 are À12.90, À6.15 and À12.97 kcal mol À1 , respectively. This indicates that 1.4 has very less tendency to associate with other molecule as compared to 1.3 and 1.5 molecule. In silico analysis was carried out for compounds 1.1, 1.2, 1.3, 1.4, and 1.5 to study the interactions and binding energies of these compounds with kinesin Eg5 protein (PDB id: 3UIH). Pymol and Discovery studio visualizer was for the analysis of docking results.
The common trend in the binding interactions of compounds 1.2, 1.5, and 1 (monastrol) in the cavity of Eg5 is observed. Ester group protrudes outside the cavity of Eg5. In silico analysis revealed that compounds 1.2 and 1.5 showed better binding energy (À7.9 kcal mol À1 ) than the standard drug monastrol (À7.8 kcal mol À1 ) ( Table 4). Chlorobenzene ring of compounds 1.2 and 1.5 is directed towards the hydrophobic region of the active site of Eg5 ( Fig. 12A and B). The chlorobenzene ring of 1.2 and 1.5 are orientated towards the hydrophobic region, so p-anion interactions with the negatively charged residue Glu116 will be favorable. Compound 1.5 also showed one halogen bond with the residue Trp127, which indicated that compound 1.5 ts nicely into the cavity of Eg5. Further, the residues Ala 133, Pro137, and Tyr211, Leu214, and Arg221 facilitate the hydrophobic p-alkyl and alkyl interactions with compound 1.5 (Fig. 12B) and (Table 4).
3.5. Cytotoxicity and anti-cancer activity of 1.1-1.    (Fig. 14B). The number of total cells decreased as the concentration of variants was increased compared to control groups at 24 h (Fig. 14A). At 200 mM concentration, 1.2 treatment caused the highest cell death at 48 h, and 1.2 and 1.5 treatments decreased total cell number to less than 50% at 48 h  (Fig. 14B). These results suggest that treatment with 1.2, 1.4, and 1.5 compounds induced signicant cell death in A549 cells. We chose 1.2 and 1.5 compounds for further study as the percentage of cell death was higher in these groups aer treatment.
3.5.3. Compounds 1.2 and 1.5 induced apoptosis in A549 cells. The compounds 1.2 and 1.5 were selected for studying apoptosis induction in lung adenocarcinoma cells because the percentage of cell death was higher with these two analogues. Acridine orange-ethidium bromide staining distinguishes live, apoptotic, and necrotic cells. Live cells appear green, and dead cells appear red. The control cells did not show any characteristics of cell death at any time point. Treatment with 1.2 and 1.5 compounds at 200 mM for 24 and 48 h showed typical apoptotic characteristics. Membrane blebbing, apoptotic body formation and nuclear fragmentation were observed at 24 and 48 h posttreatment. Treatment with 1.5 for 48 h caused necrosis in cells that appeared red (Fig. 15). These results suggest that 1.2 and 1.5 compound treatments induced cell death in lung adenocarcinoma cells.
3.5.4. Compounds 1.2 and 1.5 induced nuclear fragmentation A549 adenocarcinoma cells. Hoechst assay was performed to analyze the effect of 1.2 and 1.5 compounds on apoptosis of A549 cells. Control cells did not show any characteristics associated with cell death. In contrast, in 1.2 and 1.5 (200 mM) treated groups, DNA fragmentation, nuclear condensation, and apoptotic body formation were observed by bright blue uorescence aer 24 and 48 h of treatment (Fig. 16). These  results suggest that 1.2 and 1.5 treatments could induce apoptosis in A549 cells.
It has been observed that DHPM exists as a dimer as found in all 1.3, 1.4, and 1.5 crystal structures. Our results through the single crystal and Hirshfeld analysis found the existence of the dimeric structure of DHPM. The two such monomers lead to the dimer by forming an additional hydrogen bond yielding the bifurcated hydrogen bonds and N-H-S bridges. The presence of the thioamide group in the series exists in monomer-dimeric equilibrium. The two monomers lead to the dimer by forming by N-H-S hydrogen bond.
In addition to the N-H-S hydrogen bond (Figs. 4, 5, and 6), two more intermolecular hydrogen bonding interactions C-H/ O, N-H/O are observed between the oxygen atom of the ester group and hydrogen atom of the thioamide group. These additional H-bonds form tetramer of DHPM. The overall structure of the tetramer is close to planar. Hydrogen bonding represents one of the most versatile interactions that could be used for molecular recognition. It might be possible to interfere with the strength of the hydrogen bond effectively by linking the hydrogen-bonding site to a binding site of the molecule with the receptor. The presence of p/p stacking in compound 1.5 will also help to study the drug receptor interactions of this compound.
Both in silico and in vitro results also found the same order of anti-cancer activity. An in silico study showed that order of binding affinity is 1.2-1.5 > 1.1 > 1.4 > 1.3, whereas in vitro study showed that variants 1.2, 1.4, and 1.5 decreased the cell viability signicantly. Analogues 1.2, 1.4, and 1.5 also induced signicant cell death in A549, but 1.2 and 1.5 were found signicant in this group. Both 1.2 and 1.5 treatment caused apoptosisinduced cell death in A549 cells. Nuclear fragmentation in A549 adenocarcinoma cells by Hoechst staining showed that both 1.2 and 1.5 can induce apoptosis in cells.
As shown in tables 1.4, the values of inhibition zone for ligands are related to the nature of substituent as they increase according to the following order: p-Cl $ m-Cl > H > p-CH 3 > p-OCH 3 . This can be attributed to the fact that the effective charge experienced by the central ring increased due to de-activating substituent (p-Cl $ m-Cl). At the same time, it decreased by activating group H > p-CH 3 > p-OCH 3 . It is important to note that existence of a methyl and/or methoxy group enhances the electron density in the ring and simultaneously decreases inhibition zone values.
Electron withdrawing bridge would be expected to increase the acidity of proton donors and hence increase its binding ability as the electron-withdrawing character of a chloro group is relevant to the best analogue 1.2 and 1.5. In general, hydrogen  Inhibit proliferation of leukemia cell lines HL-60(TB) (IC 50 ¼ 0.05656 mM) and against MOLT-4 (IC50 ¼ 1.788 mM) Compound 3 Inhibit proliferation of breast cancer cell lines SK-BR-3 (GIC 50 ¼ AE0.4 mM)

1.2
Signicantly inhibit proliferation of lung adenocarcinoma cells A549 at 200 mM 1. 3 Signicantly inhibit proliferation of lung adenocarcinoma cells A549 at 200 mM bonds involving OH groups are proton donors, and their O atoms are proton acceptors. Both intra and intermolecular OH-N may form several structures in simultaneous equilibrium. The DHPM series of compounds can induce intracellular oxidation resulting in cytosolic toxicity and genome toxicity in cancer cells. These events can further cause changes at the molecular level, ultimately resulting in an altered cell cycle leading to apoptosis. These events could also reduce migration ability and impede tumor growth resulting in reduced tumor size and migration of cancer cells. These ndings suggest that 1.2 and 1.5 might be potential anticancer drug candidates. Thus, these analogues can be promising agents for treating various cancers, including lung cancer.

Conclusion
In summary, we elaborated green syntheses of dihydropyrimidinones using Citrus macroptera juice in a solvent-free reaction and deduced the primary SAR of anti-cancer. The advantages of this procedure include good yield, short reaction time, and eco-friendly procedure. The meaningful SARs can be drawn as follow: (i) deactivating group (Chlorine) at 3, 4 positions in the substituted ring of DHPM enhance activity; (ii) activating group in DHPM scaffold not suitable for active drug towards cancer cell. Moreover, we have studied the aromatic interactions using the single-crystal X-ray diffraction (SCXRD) method and the Hirshfeld surface analysis method. We found that C-H/p, C-H/O, C-H/N, C-H/C, lone pair/p, p/p interactions for synthesized compounds were benecial for drug design, stabilization of DNA structure, protein folding, crystal engineering, material science, etc. Hirshfeld's surface analysis showed various intermolecular interactions in the crystal structures supporting the single-crystal X-ray diffraction (SCXRD) study. DFT calculations show that substituent does not impact the geometry of synthesized DHPM derivatives and conformational study of molecule reveals that it can change its conformation to get best interaction with active site. By these studies, exposure of different groups in analogues in the solidstate provides a clue for further study for biological activity concerning different receptors. As far as cytotoxicity and anticancer activity of these compounds are concerned, we found that they showed good anticancer activity. The anticancer activity was carried out in Human lung adenocarcinoma (A549) cells, where two derivative compounds (1.2 and 1.5) from Scheme 1 showed anticancer activity. The anticancer activity of these compounds (1.2 and 1.5) was well supported by in silico analysis, where they showed better binding energy than standard drug monastrol. Compound 1.4 showed poor anticancer activity compare to 1.2 and 1.5 molecules due to weak tendency of weak interactions as per nding in DFT calculation. Thus, anticancer activities suggest that our synthesized compounds could be potential candidates to develop anticancer drugs.

Conflicts of interest
There are no conicts to declare.