Abstract

Four pyrazolone derivatives and their corresponding silver complexes were synthesized and characterized. Based on elemental analysis, 1 : 2 (M : L) molar ratio was suggested for all inspected complexes. 1H, 13C NMR, mass, UV-Vis, TGA, and IR were the spectral tools used for describing the formulae. Moreover, XRD patterns and SEM pictures were used to evaluate the particle sizes which appeared strongly in nanometer range. CT-DNA study is the major consideration in this study, to test the interacting ability of all synthesized cationic complexes towards cell DNA. Each binding constant was computed and correlated with the Hammett sigma constant. Antitumor activity was examined upon three carcinoma cell lines (MCF-7, HepG2, and HCT116). The high efficiency was recorded towards MCF-7 (breast carcinoma) cell line. Kinetic studies yield essential parameters to assert on the rule of metal atom on thermal feature of organic compounds. Molecular modeling was implemented to optimize the structures of compounds. Also, molecular docking was achieved to obtain a clear view about proposed drug behavior within the affected cells. This was achieved through comparing the calculated internal energy values of all docking complexes. All the tested compounds displayed a significant interaction with breast cancer protein (strong matching with practical result) followed by DNA polymerase protein.

1. Introduction

After discovering the cisplatin as an effective antitumor drug, the interest of researchers in the coordination chemistry increased to design highly effective complexes as antiantitumor operators [13]. But, the limitations of the clinical success of cisplatin due to critical side effects and resistance that cause retrogression [4] have stimulated searches for other metal complexes. The synthesis for series silver(I) complexes as antitumor agents conquers a considerable place in cancer chemotherapy, as evidenced by their remarkable therapeutic potential reflected in many reported research works [5]. The silver ion is of considerable solicitude in the design of coordination complexes due to the resilience of its coordination sphere and the ease with which it varies its coordination number (from 2 to 4, rarely 5 or 6). Coordination complexes of silver are known to be convertible through changes in ligand geometry, rigidity or functionality, and/or by modifications to the counter-ion or solvent system [6]. But, silver often tends to prepare coordination polymers [7, 8]. The main factor in the effectiveness of anti-cancer compounds is their ability to interact with the DNA of the cancer cell which prevents its division until the cell will die [9]. Ag+ is uniquely interesting because it binds exclusively to the bases rather than the backbone of DNA [10, 11]. On the other hand, the wide range of biological activities of pyrazolone as antitumor [12], anti-inflammatory [13] and antimicrobial [14] makes it a basic objective to synthesize biologically active transition metal complexes as it is known that the formation of the complexes enhances the biological efficiency of the ligand [15]. Based on all of the afore-mention, we focused in this work on the synthesis of nonpolymeric silver complexes with pyrazolone to improve its efficiency as anticancer drug for three carcinoma cell lines and to study their interaction with DNA. The structure of novel complexes will be investigated through all possible spectral and analytical tools. Computational studies (i.e., Gaussian 09 molecular modeling, QSAR, and Autodock tools) will be applied to prove and interpret the efficiency of novel complexes.

2. Experimental Work

2.1. Chemical Reagents

All chemicals used to prepare pyrazolone derivatives (ethyl cyanoacetate, aniline, 4-chloroaniline, 4-toluidine, 4-nitroaniline, hydrazine hydrate, ethanol, sodium acetate trihydrate, sodium nitrite, and HCl) were used as received without any pretreatment. They were purchased from Fluka and Sigma-Aldrich. Also, AgNO3 used for preparing the complexes was PDH.

2.2. Preparations
2.2.1. Synthesis of 3-Amino-4-(2-substituted phenylhydrazono)-1-H-pyrazol-5(4H)-one

Derivatives of pyrazolone (HL1–4) were synthesized and fully characterized in accordance with the previous studies [16, 17]. A fresh prepared solution of cyanoacetic hydrazide (0.99 g, 0.01 mol) in ethanol (50 ml) was coupled in presence of sodium acetate trihydrate (1.36 g, 0.01 mol). The coupling was carried out on 0.01 mol arenediazonium chloride (prepared by diazotizing corresponding aniline derivative) in 6 M hydrochloric acid (6 ml) with cold sodium nitrite solution (10 ml, 1 M; as known) in ice bath (at 0–5°C). The reaction mixture was set in a refrigerator for 3 hrs. The precipitating compound was filtered off and dried. Refluxing ethanolic arylhydrazone solution for 2 hr gives the aimed ligands, 3-amino-4-(2-substituted phenylhydrazono)-4,5-dihydropyrazol-5-ones (HL1–4) as appeared in Scheme 1. Whereas, their scientific names are HL1, 3-amino-4-(2-phenylhydrazono)-1H-pyrazol-5(4H)-one; HL2, 3-amino-4-(2-(4-chlorophenyl)hydrazono)-1H-pyrazol-5(4H)-one; HL3, 3-amino-4-(2-(4-mehylphenyl)hydrazono)-1H-pyrazol-5(4H)-one and HL4, 3-amino-4-(2-(4-nitrophenyl)hydrazono)-1H-pyrazol-5(4H)-one. The prepared pyrazolone compounds were structurally optimized to clarify the orientation of active groups (Figure 1) using Gaussian 09 software. 1H NMR (DMSO-d6) spectra of four ligands (Figure 1S) were assigned as follows: HL1, δ: 5.80 (s, 2H, NH2), 7.10–7.54 (m, 6H, Ar–H, and amidic-NH), 10.53 (s, 1H, NH). HL2, δ: 5.84 (s, 2H, NH2), 7.42–7.59 (m, 5H, Ar–H and amidic-NH), 10.53 (s, 1H, NH). HL3, δ: 2.23 (s, 3H, CH3), 5.71 (s, 2H, NH2), 7.01–7.60 (m, 5H, Ar–H and amidic-NH), 11.41 (s, 1H, NH), and HL4, δ: 6.02 (s, 2H, NH2), 7.73–8.25 (m, 5H, Ar–H and amidic-NH), 10.66 (s, 1H, NH). Also, mass spectral analysis of HL2 and HL3 ligands (Figure 2S) displays the following fragmentation peaks starting with molecular ion; m/z = 239 (M+ + 2, 21.76); 238 (M+ + 1, 8.61); 237 (M+, 64.69); 126 (100); 111 (26.80); 98 (36.92); 91 (9.09) and 218 (M+ + 1, 13); 217 (M+, 100); 126 (39.26); 111 (1.32); 106 (38.82); 77 (46.96), respectively.

2.2.2. Synthesis of Silver(I) Complexes

2 mmol AgNO3 was dissolved in a little amount of ethanol. Then, the solution was added gradually to 2 mmol from each ligand solution (in DMF). The mixture was stirred and refluxed at 50–60°C for 15–30 min. The complexes appeared in the reaction medium without external additions. The main reaction product was precipitated, then filtered off and washed with ethanol followed by diethyl ether, and after that, dried over CaCl2 desiccators.

2.3. DNA Binding Methodology

The binding efficiency of calf thymus-DNA (CT-DNA) has been deliberated by spectrophotometric titration method. 50 mg from CT-DNA were dissolved by continuous stirring overnight in deionized water (pH = 7.0) and preserved at 4°C. Tris(hydroxymethyl)-aminomethane was used to prepare 5.0 mM tris-HCl buffer (pH = 7.2) after adding 50 mM NaCl in deionized water. The prepared buffering solution is perfectly free from protein, based on the value measured for absorbance ratio (1.8–1.9; A260/A280) at 260 and 280 nm [18]. The concentration of prepared DNA (5.25 × 10−4 M) was computed by knowing molar absorptivity coefficient which is 6600 M−1·cm−1 at 260 nm. 2.5 × 10–5 M was the concentration used from each complex. At room temperature, the UV-Vis scanning process was carried out at the fixed range (200–900 nm) by 1 cm quartz cuvette. A gradual increment for CT-DNA concentration from 0.00 till 1.58 × 10−4 mol·L–1 was done, towards a fixed concentration from each tested compound. The same CT-DNA amount was added to a reference buffer solution with each measurement, to omit the absorbance of free DNA. Intrinsic binding constant (Kb) of the interaction between compound and CT-DNA was computed by; [DNA]/(єa − єf) = [DNA]/(єb − єf) + 1/Kb (єa − єf), equation [19]. Where, [DNA] is the molar concentration from DNA, єa is the extinction coefficient observed at each DNA concentration by; Aobs/[compound], єf is the extinction coefficient for each free tested compound in solution. Moreover, єb is the extinction coefficient for the fully bonded compound towards DNA. By drawing the relation of [DNA]/(єa − єf) versus [DNA], Kb is given by the following ratio: slope/intercept.

2.4. In Vitro Antitumor Study

The biological feature of synthesized pyrazolone derivatives (HL1–4) and their Ag(I) complexes were screened towards three carcinoma cell lines (MCF-7, HEPG-2, and HCT-116) by sulforhodamine B (SRB). Cells were kept in RPMI medium enriched with 100 µg/ml streptomycin, 100 units/ml penicillin, and 10% heat-inactivated genetic bovine serum inhibit, 5% (v/v) CO2 atmosphere at 37°C, and cells were subcultured two times/week. The growing cells were gathered using trypsin–EDTA (0.25%) and plated in 96-well plates at 1000 cells/well. Cells were treated to extraction for 72 h and subsequently fixed with TCA (10%) for 1 h at 4°C. After washing times, cells were exposed to SRB solution (0.4%) for 10 min in a dark place and subsequently washed with glacial acetic acid (1%). After keeping overnight, Tris–HCl was used to dissolve SRB-stained cells, and the intensity of color was measured to be 540 nm. The viability of cell was measured by Trypan blue (0.4) stain, which permits us to differentiate viable (nonstained) from dead (stained) cells. Viable and dead cells were counted with Nikon microscope implementing hemocytometer. Surviving cells percentage = (1.00 − (number of blue cells/number of total cells)) × 100. Statistical analysis was executed by using Sigma Plot version 12.0.

2.5. Equipments and Physical Analysis

Carbon, hydrogen, and nitrogen contents were performed at Perkin-Elmer 2400 CHN elemental analyzer. The silver content was evaluated gravimetrically using standard method [20]. The molar conductivities of metal complexes were obtained by JENWAY model 4070 conductance bridge (in DMSO solvent). KBr-IR and 13C, 1H NMR spectra were obtained by JASCO FT-IR-4100 spectrophotometer (400–4000 cm−1) and Burker 500 MHz, respectively. The purity of most compounds was checked from mass spectra at 70 eV by using AEIMS 30 mass spectrometer with heating rate 40°C/min and mass range (50–1000). Electronic spectra were recorded by UV2 Unicam UV/Vis spectrophotometer at room temperature using 1 cm quartz cell (in DMSO solvent). Thermal analysis was accomplished using Shimadzu thermogravimetric analyzer by heating rate of 10°C min−1 under nitrogen (20–900°C). X-ray diffraction patterns (XRD) were accomplished by Rigaku diffractometer using Cu/Kα radiation. Scanning electron microscopy (SEM) images were extracted by the Joel JSM-6390 equipment. Antitumor activity was tested in microbiology center.

2.6. Computational Implementation
2.6.1. Kinetics

Essential thermodynamic parameters were extracted over two degradation stages from TGA curves. All degradation processes were carried out at a satisfactory rate which introduces a good computation for kinetic parameters along most degradation stages. The kinetic parameters were calculated for the first and second degradation stages. The chosen steps exhibit the decomposition for great masses from the coordination sphere, and order (n) and energy of activation (E) were computed. Many researches were dealing with kinetic calculations [2129] and proposed equations stated for such purpose. Among these methods, Coats and Redfern [23] and Horowitz and Metzger [28]were used in this study.

2.6.2. Structural Optimization

Utilizing Gaussian 09 software [30], the structural optimization for pyrazolone derivatives (HL1–4) and their Ag(I) complexes were executed. DFT/B3LYP is a proper method achieving best configurations and computational files. The files extracted from the optimization process (log and chk files) were visualized by Gauss-View screen [31] to produce essential parameters. Also, HOMO and LUMO images were extracted from chk file visualization.

2.6.3. QSAR Calculation

QSAR parameters calculations were carried out over all synthesized compounds after optimization process by applying HyperChem software program (8.1). Firstly, preoptimization was accomplished by soft adjustment procedures, through molecular mechanics force field (MM+) followed by semiempirical (AM1). The program implemented over structural forms without fixing any parameter to execute the optimization at equilibrium. The minimization of energy was done by Polak–Ribière conjugated gradient algorithm. QSAR parameters were computed for all derivates and their Ag(I) complexes, for comparison.

2.6.4. Molecular Docking Methodology

The docking process was executed between protein receptors (3s7s, 3e1r, and 4dk7 assign for breast, colon, and liver cancers, respectively) and pyrazolone derivatives (HL1–4) (proposed inhibitors). Moreover, the docking process was done over CT-DNA polymerase (1bpb). The selected protein receptors were the same for the pathogen of tumor cells investigated in the experimental part to assert on the results. Autodock tools 4.2 is a software used for docking methodology. Gasteiger partial charges were added over elements of tested compounds (pyrazolones). The rotatable bonds were omitted and nonpolar hydrogen atoms were ruling out. The energy levels of docked complexes were computed after addition of fundamental hydrogen atoms, Kollman united atom type charges, and salvation parameters [32]. Affinity (grid) maps of ××Å grid points and 0.375 Å spacing were generated utilizing Autogrid program [33]. Also, Vander Waal forces and electrostatic terms were acquired, through the use of autodock parameter set-dependent and distance-dependent dielectric functions, respectively. The docking was accomplished using Solis and Wets local search method and Lamarckian Genetic algorithm (LGA) [34]. Initial position, orientation, and torsions of tested molecule (inhibitor) were set indiscriminately. All rotatable torsions were expelled during the docking process. Each experiment is the average value for 10 different runs that were set to close after maximum of 250000 energy assessments. 150 is the used population size. During the process, translational step of 0.2 Å, quaternion, and torsion steps of 5 were used.

3. Results and Discussion

3.1. Physicochemical Properties

Common specific features are aggregated in Table 1. All synthesized complexes have high melting points (˃300°C), while that of free organic derivatives (HL1–4) were 246, 300, 252, and 300°C. The analytical results are proper with 1 : 2 molar ratio between silver ion and chelating ligands. The molar conductivity (in DMSO) which measured for all complexes (1 mmol) is coinciding with conductivity feature. The conductivity values (Λm) are confined in the range of conducting feature for monoanionic conjugated nitrate with cationic coordination sphere [35]. This is proper feature with high acidic anions as nitrate.

3.2. IR and 1H, 13C NMR Assignments

IR spectral analysis of pyrazolone (HL1–4) derivatives (Table 1S) displays significant bands assigned for the following vibrations; νas(NH2), νs(NH2), ν(NHs), δ(NH2), δ(NH), ν(C=N), and ν(C=O). Observable lower shifted appearance for ν(NH2), δ(NH2), and ν(C=N) bands suggests the bidentate contribution mode of pyrazolones towards silver atom. However, more or less shifted appearance of NHs and C=O donors introduces their ruling out from coordination [36]. Each coordination sphere includes two coordinating ligands towards one central atom, in which the orientation of donors according to central atom is the cause for chelating mode (by NH2 and C=N)and also has the major impact for ruling out of other donors (NHs and C=O). The appearance of new band at ≈1383 cm−1 refers to ν(NO3) about the ionic attachment of nitrate. Moreover, new band appeared at the range 509–560 cm−1 in all complexes, refers to ν(M–N) vibration [3739]. 1H and 13C NMR (DMSO-d6) spectral data of investigated complexes are tabulated (Table 2). 1H NMR spectra (Figure 3S) display a considerable down field shift for NH2 group by an observable ppm value, whereas 13C NMR spectral data represent the absence of significant alteration in hybridization attitude upon carbon atoms inside the forms. This is omitting hydrogen relocation through tautomerism followed by the enolization process.

3.3. Electronic Spectra

The preparation of AgNPs has been assured through UV-Vis spectral analysis, which carried out for all complexes in the DMSO solvent (10−3 M). The measurements process was accomplished at 510, 440, 420, and 480 nm for four Ag(I)-HL (1–4) complexes, respectively. The excessive presence for AgNPs was expected from the deep color of Ag(I) complexes. The dark color of complex solution points to the formation of AgNPs. The electronic spectra of four complexes reflect a well-defined Plasmon band between 420 and 510 nm in the nanoscale range, according to silver nanoparticles association. According to theoretical and experimental studies [40], the optical features of AgNPs were reported based on electronic spectral band position. The band at 420 nm for Ag(I)-HL3 complex solution is referring to the presence of small spherical silver nanoparticles. The spherical particles display λmax band at shorter wave length, while the nonspherical shape of particles shifts the λmax to longer wave length [41]. The higher the intensity of the complex color, the higher the concentration of AgNPs [42]. The tetrahedral configuration of Ag(I) complexes (Figure 2) was proposed based on analytical, IR, 1H, 13C NMR, and mass spectral analysis, which is acceptable for d10 systems.

3.4. Microanalysis Mass Spectra

Two organic ligands were chosen for this study (HL2 and HL3), the molecular ion peaks recorded were m/z = 239 (M+ + 2, calcd. 237.63) and 218 (M+ + 1, calcd. 217.23), respectively. A high agreement between each recorded M+ and the calculated one supports high purity. Three silver complexes were studied to assert on the formulae proposed (Figure 4S). In the [Ag(HL2)2](NO3) complex spectrum, a well-defined parent peak was appeared at m/z = 647.66 (M+ + 2 (1.2), calcd. 645.15) is followed by a successive fragmentation corresponding to peaks at 344 (3.81), 243 (11.78), 238 (7.55), 220 (14.50), 197 (81.93), 140 (15.61), 127 (39.25), 112 (7.27), 104 (16.64), 95 (15.05), and 95 (65.29). Also, in the [Ag(HL3)2](NO3)2H2O complex spectrum, a significant ion peak was recorded at m/z = 605 (25.94, calcd. 640.36) which is attributing to M+ + 1 − (2H2O). This is not considered a strange behavior; the sudden expel for water molecules may happen during vaporization. The vaporization is a prefragmentation process before electron bombardment [43]. Essential fragmentation peaks recorded are 546 (54.26), 386 (5.28), 370 (100), 268 (12.45), 243 (19.73), and 202 (1.06). Moreover, in the [Ag(HL4)2](NO3) complex spectrum, the recorded molecular ion peak attributes to M+ − 1 (m/z = 665.22) (calcd. 666.27). The more proper fragmentation perception was displayed in Scheme 2.

3.5. Thermogravimetric Analysis

A proper degradation path for all pyrazolone derivatives and their corresponding silver complexes is tabulated (Table 2S). The TGA curves of pyrazolone ligands display two significant degradation stages starting from ≈150°C and ended to ≈620°C. The successive degradation process had three to four stages (ended at ≈800°C). All the coordination spheres around the silver atom were completely destroyed along the temperature range applied. Moreover, the residue includes AgN3 or silver atom meanwhile polluted by carbons. A satisfactory agreement between calculated and found weight loss corresponding to the expelled particles reflects the exact assignments for step borders, which suggests that well computed kinetics may be obtained.

3.6. XRD and SEM Analysis

XRD patterns of pyrazolone ligands and their silver complexes (solid samples) were recorded (Figures 3 and 5S). This analysis gives a good insight about the crystal lattice for tested compounds. Applying referenced methods, the purity of each compound can be asserted through comparative study with the reactants patterns [44]. Moreover, essential lattice parameters were deducted from high intense peak in patterns (FWHM). All tested compounds show nanocrystalline feature, which clearly impacted on the sharpness of patterns peaks. This may be due to the fixation of compounds’ nucleus, as the ligands family and silver central atom. Particle sizes, 2θ, relative intensity (%), d spacing, and FWHM are represented in Table 3. The sizes of crystals were estimated from the Deby–Scherrer equation: β = 0.94λ/(S cos θ), where S is the crystallite size, θ is the angle of diffraction, and β is FWHM, Cu/Kα (λ) = 1.5406 Å.

The d-spacing upon two following layers was extracted from the Bragg equation:  = 2d sin(θ) at . The crystallite sizes calculated for all tested compounds appeared excellently in nanometer range. This feature devotes the application towards antitumor activity as well as DNA binding study. Also, ε index (crystal strain) and δ index (dislocation density) were calculated by known relations [45] equation. The two indexes are reflecting the dislocation inside the compound network. The lower their values, the higher the quality of compounds. Also, Miller indexes (hkl) were calculated on the bases of cubic feature for investigated crystals. The calculated values display mainly the parallelism of planes towards y- and z-axes while, intercepts with x-axis, except for Ag(I)-HL2 and Ag(I)-HL4 complexes. This reflects the relatively perfect crystal building towards the axes with most compounds. SEM is another technique used to obtain a clear view about the surface topography for tested compounds (Figure 6S). The images of four organic ligands were found as needles or sticks by a regular way. Also, the high aggregation observed with HL4 ligand may reflect its distinguished nanometer size of separate crystals which coincide with XRD data. However, the images of complexes are completely different from that of corresponding ligands, which initially reflects the purity. Moreover, the spherical shape is the feature observed, which indicates the impact of metal central atom in the aggregation process of particles.

3.7. DNA Binding Studies

The degree of interaction of Ag(I) complexes against CT-DNA was examined by the spectrophotometric titration method. The freshly prepared solutions (as reported in experimental part) were scanned along the applied wave length range (200–900 nm) against referenced solutions (blank). The intrinsic binding constants (Kb) were calculated at λmax for each charge transfer band; 364.7, 366.4, 365.1, and 458 nm attributing to the binding yields of silver complexes (1–4) and DNA. The absorption relationships for all tested compounds are displayed in Figure 7S. An observable bathochromic effect for new CT band is conjugated with progressive increase DNA amount [46]. The red shift observed is found in the range 1–4 nm, with all investigated compounds. This indicates the formation of new compounds yielded from the binding with DNA helix. Also, a gradual altitude in the absorption is the other effective indicator. The higher the binding constant, the higher the stability of coupled DNA helix. The coupling process is mainly proper with occluding inside major and minor grooves of DNA, whereas the coupling process attributing to silver complexes may prolonged by electrostatic attraction for cationic coordination sphere as well as occlusion inside grooves. The coupling stabilization may be facilitated by the broadness of occluded compound surface as well as a charged coordination sphere (with silver complex). The occlusion process causes reorganization inside CT-DNA, which required partial disassembly or deterioration for the double helix at the exterior phosphate leading to form cavity suitable for entering [47]. The bathochromic shift is corelated with the substituent’s effect by a direct relation with electron withdrawing character (σR). Hammett’s equation established the significant effect of p-substituent’s on the binding constant value (σR and Kb) [48]. The constant of binding was arranged by 1.29 × 104, 1.82 × 104, 3.57 × 103, and 2.4 × 104, attributing to complexes for HL1–4, respectively, in agreement with Hammett’s relation. The binding constants reflect an acceptable interaction between the complexes and CT-DNA. The higher the inductive effect, the higher the interacting ability with DNA helix. Other essential plots against σR (Figure 4) were extracted to clarify the relationships as; ΔE, PKa, (docking parameter), and log P. A direct relation was shinned with ΔE (energy gaps) and PKa plots, while a reverse with log P. The effective impact for log P value is appearing with the minimization of its value.

3.8. Antitumor Study

All synthesized complexes were screened against three cell lines (HEPG2, MCF-7, and HCT116) to establish a view about their antitumor activity. The obtained results propose the cytotoxic feature of complexes towards the tested cell lines in comparison with congruent ligands. The IC50 values are represented in Table 4 and graphed in Figure 8S. The cytotoxic indicator IC50 is considered a good anticancer parameter on behalf of other significant indexes such as DNA fragmentation or morphology of the cell. Based on this vision, experiments were conducted on all compounds to identify their impact on reducing cells or not through morphological index. To achieve that, cells were treated with IC50 for selected high impacted compounds (HL3-4 and their complexes). MCF-7 breast cancer was the chosen cell line. After treatment with IC50 of selected compounds, MCF-7 cells were incubated for 24 hours. Cells were fixed using Carnoy’s fixative and stained by fluorescent stain [49]. The estimation was based on the nucleus shape, nucleus density, foci presence, and the number of cells. The outcome data clearly showed that nucleus remained in its regular shape and that there were no condensation and no degradation for detected nucleus (Figure 8S). Lack of micronucleus inside cytoplasm revealed that ligands and their Ag(I) complexes are not clastogenic agents, which may consider a very good feature. Also, cells preserved their capability to produce foci which means that the compounds may have no impact towards penetrating layers and have no ability to solve cell-to-cell connections. This may be considered useful for them because they do not help cells to spread, where they can be killed in situ (Figure 5). In terms of cell numbers, picture has assured the compatibility and regularity with the datum proposed by SRB assay. In the light of former results and discussion, ligands and complexes may be treated as in situ cytotoxic but not genototoxic.

3.9. Computational Utilizations
3.9.1. Kinetics

The kinetics calculations are a significant study that gives clear view about the impact of metal atom occluded with organic molecules upon the thermal degradation attitude over whole compound (Table 5). E, ∆H, ∆S, and ∆G are the parameters computed over two degradation stages in each investigating compound. These parameters were calculated based on evaluating fraction decompose (α) at variable temperatures inside step borders. Coats and Redfern [23] and Horowitz and Metzger equations [28] are the two comparative methods used. The relationship graphs are exhibited in Figures 9AS and 9BS. ∆H = E − RT and ∆G = ∆H− TS were the equations used. Entropy of activation (∆S) was extracted from the equation , where h, Kb, and Ts are Plank’s constant, Boltzmann constant, and DTG peak midpoint temperature, respectively. The following notes are grouped into calculated values: (i) The values of activation energy are suffering reduction with complexes in comparison with free ligands. This is a good evidence for the weakness appearing on the bonds of organic compound as a result of metal coordination. (ii) There is an observable increase in the activation energy from one step to another due to the reduction of decomposition rate (endothermic steps) from the original molecule. This will lead to increase of ∆G values in the following stages. (iii) The altitude of activation entropies (∆S, by negative sign) for complexes upon their corresponding ligands reflects the high-ordered complexes than the free ligands and/or slow reaction rate.

3.9.2. Structural Optimization Parameters

Applying Gaussian 09 software, the best orientation for atoms inside the structures was extracted. DFT/B3LYP is a proper method for this purpose. The optimization procedure exerted essential output files (chk and log). The visualization for chk file, HOMO, and LUMO images and their energy values were obtained. The frontier, HOMO, and LUMO images belong to free ligands and their complexes are displayed in Figure 6 and Figures 10AS and 10BS. The HOMO level for pyrazolone ligands were centered on pyrazol-5(4H)-one, moiety which includes the atoms interested in bonding. While, the LUMO level were spread over a whole molecule. However, the attitude of HOMO and LUMO levels in silver complexes were centered on the metal atom. This is normal with such electron-dense atom as silver. All effective parameters tabulated in Table 3S were computed concerning to frontier energy gaps (ELUMO – EHOMO, ΔE). The indexes χ (electronegativity), μ (chemical potential), η (global hardness), S (global softness), ω (electrophilicity), and ϭ (absolute softness) were calculated using referenced equations [50, 51]. Using electrophilicity index (ω), the toxicity and reactivity of tested compounds can be evaluated. This will give an acceptable expectation for the biological feature of compounds by which, the excellence for complexes upon their corresponding ligands will be expected with the progression for Ag(I)-HL4 complex. As well as, the constancy and electron affinity can be evaluated from two opposing indexes; η and ϭ [52, 53]. Also, the superiority was recorded with the complexes in comparison with their congruous ligands. Other two opposing indexes (S and η) attributing to the degree of softness displayed the soft character of all compounds with the priority of complexes. This indicates an expected distinct biological behavior for complexes. Furthermore, frontier energy gaps in complexes were reduced from original ligands. This is due to the effective impact of metal atom, which is favored in biological activity [54].

3.9.3. Considerable Log File Parameters and QSAR

Other additive parameters were deducted from log file data, from which important data are drawn (Table 4S) that give accurate specifications for structural form and confirmed it. Significant charges of coordinating atoms (N8, N15, N36, and N43) before (free ligand) and after complexation suffer reduction overall. Moreover, the net charge of silver ion is highly reduced due to charge transfer accompanied with coordination (L → M). Two bond lengths contributing with donor atoms (C11–N15 and C12–N8) were estimated for ligands only, while they could not be technically computed for complexes. So, effective comparison was absent. Oscillator strength (ʄ) was extracted for all studied compounds and compared. An overall lessening was recorded with the complexes values. This is a good indicator for facilitating electronic transition inside the whole complexes molecules [55]. The dipole moment values were extracted and display variable features in comparative view. Two of them suffer enhancement in complexes and the others significantly reduced. The enhancement may reflect increase in polarity over a molecule, whilst, the reduction purpose a reverse feature. QSAR parameters (Table 6) were obtained after optimizing the structures of all tested compounds using HyperChem (8.1) program. Surface area and the volume of tested compounds showed a significant increase in complexes compared to ligands. Also, reactivity and polarizability were generally enhanced in complexes with the priority of Ag(I)-HL4 complex in agreement with the previously displayed. Moreover, log P values were estimated and displayed excellence with Ag(I)-HL4 complex and its ligand. Partition coefficient (log P) points to the biological feature of tested compound by a reverse relation.

3.9.4. Docking Interpretation

In the last few decades, a great revolution is seen in drug design industry, which extended to preemptive theoretical studies to show the expected biological efficiency of proposed drug. Autodock 4.2. tools software is one of these tools used for this purpose. The program designed for simulation about the interaction between proposed drug (inhibitor) and protein receptors attributes to infected cells. 3s7s, 3e1r, 4dk7, and 1bpb were the receptors used for breast, colon, and liver cancers besides CT-DNA polymerase. The docking feature of all organic compounds was studied to evaluate the inhibition activity towards the selected proteins. The output internal energies computed over PDB files are tabulated (Table 5S). pKa (dissociation constant) value was calculated for each compound, which is considered a biopharmaceutical measure for drug-likeness feature. This constant is reflecting the degree of ionization for compound. The high pKa value means high degree of polarity around hydrogen atoms of function groups. This will lead to increase in the number of H-bonding between the proposed drug (ligand) and the protein receptor. Also, a direct relation between pKa value and σR (Figure 4) goes parallel with the former aspect. Computed different energies showed distinct correlations for docked drugs (ligands) with breast cancer followed by DNA polymerase proteins [56]. The best docking complexes (Figure 7),(3s7s-HL1–4) display high occlusion for drugs (ligands) which centered inside protein helix. Also, the binding towards DNA (1bpb) is highly satisfactory for the inhibition feature while the other docking complexes with 3e1r and 4dk7 proteins display binding in moderate level (Figure 11S). 2D plots (Figure 12S) verify the intense H-bonding with 3s7s protein in comparison with others. This study shows the promising inhibition activity towards breast cancer as well as CT-DNA which strongly and strangely agreeable with practical data that reported in our study.

4. Conclusion

A series of tetrahedral Ag(I)-pyrazolone complexes were synthesized and characterized. All compounds investigated were found in the nanometer range based on UV-Vis and XRD spectra. CT-DNA study yields intrinsic binding constants by the arrangement attributing to Hammett’s postulation. Also, antitumor activity offered a promising feature for the Ag(I)-HL4 complex. Structural formulae of all synthesized compounds were optimized theoretically. Also, simulation for the interaction between tested drugs and types from infected cell proteins was achieved.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

Dr. Ismail Althagafi and coauthors would like to thank Deanship of Scientific Research at the Umm Al-Qura University (project ID: 43405078) for the financial support to this research extracted from the project.

Supplementary Materials

Table 1S: characteristic IR spectral bands (cm−1) for Ag(I)-pyrazolone (HL1–4) complexes. Table 2S: plausible TGA data for Ag(I)-pyrazolone (HL1–4) complexes. Table 3S: energy parameters (eV) applying DFT method. Table 4S: significant log file data for pyrazolone ligands and their Ag(I) complexes. Table 5S: energy values (kcal/mol) for pyrazolone-protein docked complexes. Figure 1S: 1H NMR spectra of pyrazolone derivatives, HL(1–4). Figure 2S: mass spectra of two elected ligands. Figure 3S: 1H NMR spectra of three Ag(I)-pyrazolone derivatives. Figure 4S: mass spectra of selected ligands and two Ag(I) complexes. Figure 5S: X-ray patterns of three pyrazolone derivatives and their Ag(I) complexes. Figure 6S: SEM images of pyrazolone derivatives, HL(1–4) and their Ag(I) complexes. Figure 7S: spectrophotometric plots for DNA binding towards Ag(I) complexes (1–4). Figure 8S: dose response curves of HL (1–4) ligands and their Ag(I) complexes against MCF-7, HEPG-2 and HCT-116 cancer cells, cells were treated with various concentrations of ligands and complexes incubated for 72 hrs. %5 CO2. Figure 9AS: kinetic plots using Coats–Redfern and Horowitz–Metzger for pyrazolone ligands, HL (1–4). Figure 9BS: kinetic plots using Coats–Redfern and Horowitz–Metzger for pyrazolone-Ag(I) complexes. Figure 10AS: images of frontier molecular orbital’s (HOMO,1 and LUMO,2) for pyrazolone ligands (HL1–3). Figure 10BS: images of frontier molecular orbital’s (HOMO,1 and LUMO,2) for Ag(I)-HL1–3 complexes (A–C, respectively). Figure 11S: interacting complexes appeared with 3e1r, 4dk7, and 1bpb proteins and HL1–4 ligands. Figure 12S: 2D plots for interacting complexes appeared with 3e1r, 4dk7, 1bpb, and 3s7s proteins proteins and HL1–4 ligands. (Supplementary Materials)