key: cord-1018987-kvuufoxv authors: DERAFA, Wassila; AGGOUN, Djouhra; MESSASMA, Zakia; HOUCHI, Selma; BOUACIDA, Sofiane; OURARI, Ali title: An Unexpected Single Crystal Structure of Nickel(II) Complex: Spectral, DFT, NLO, Magnetic and Molecular Docking Studies. date: 2022-04-30 journal: J Mol Struct DOI: 10.1016/j.molstruc.2022.133190 sha: bcdd437a06f31e91a88c140eab3f2609932f8180 doc_id: 1018987 cord_uid: kvuufoxv This work explores the study of a synthesized nickel complex as a possible inhibitor against the main protease (Mpro) of the recent emerging coronavirus disease (COVID-19). Overall, the template reaction of 3-acetyl-2-hydroxy-6-methyl-4H-pyran-4-one with nickel(II) chloride hexahydrate in N,N-dimethylformamide (DMF) medium leads to the formation of neutral nickel complex. This resulting complex is formulated as [Ni(DHA)(2)(DMF)(2)] on the basis of FT-IR, UV-Vis., single-crystal X-ray diffraction analysis, magnetic susceptibility and CV measurements as well as DFT quantum chemical calculations. Its single crystal suggests was found to be surrounded by the both pairs of molecules of DHA and DMF through six oxygen atoms with octahedral coordination sphere. The obtained magnetic susceptibilities are positive and agree with its paramagnetic state. In addition to the experimental investigations, optimized geometry, spectroscopic and electronic properties were also performed using DFT calculation with B3LYP/6-31G(d,p) level of theory. The nonlinear optical (NLO) properties of this complex are again examined. Some suitable quantum descriptors (E(HOMO), E(LUMO), Energy gap, Global hardness), Milliken atomic charge, Electrophilic potion and Molecular Electrostatic Potential) have been elegantly described. Molecular docking results demonstrated that the docked nickel complex displayed remarkable binding energy with Mpro. Besides, important molecular properties and ADME pharmacokinetic profiles of possible Mpro inhibitors were assessed by in silico prediction. Wuhan city, Hubei province of China, December 2019, a novel infectious Coronavirus disease (formally called COVID-19) caused by SARS-CoV-2 has been reported and raised as a global health anxiety. Two months later on March 11, 2020 , World Health Organization (WHO) declared it as a pandemic after its hurried propagation in Chinese territory and also spread through the world [1] . Coronavirus, positive and single-stranded RNA, is one of the most viruses that can infect both humans and animals [2] . Until now, the most identified promising targets are the spike protein, RNA-dependent RNA polymerase (RdRp), and the 3-chymotrypsin-like protease (3CLpro) otherwise known as main protease (Mpro) [3, 4] . It cleaves the polypeptide in to various non-structural proteins which are involved in the viral replication. Hence, it is essential for processing the polyprotein that led to the proteolytic activation of the viral functional proteins [5] . Consequently, now, the main promising target for researcher community is to obtain novel treatment of coronavirus disease based on synthetic compounds. In coordination chemistry, dehydroacetic acid or 3-acetyl-2-hydroxy-6-methyl-4H-pyran-4-one (Abbreviated as DHA) and its derivatives are excellent reactive heterocyclic compounds [6] . Against a lot of organic compounds, DHA can serve as an important pharmacologically active intermediate in organic synthesis because of its several reactive functional groups compared to electrophilic and nucleophilic reagents [7] . Furthermore, this property makes of them a very promising advantage for researchers due to their wide-ranging biological properties such as fungicidal [8] , bactericidal [9] , antibiotic [10] and as well for a lot of their use in food technology preservative [11] . Additionally, major advances were made recently in the development of anti-viral therapy of DHA derivatives. As example, these compounds act efficiently in inhibiting diverse targets and as potential HIV protease inhibitors due to the interaction of dehydroacetic acid derivatives against HIV-1 integrase [12] . Other numerous properties and applications have also been known for DHA with its derivatives. These compounds are harmless as cosmetic ingredients in the current practice [13] and they are commonly known for their use in the enhancement and stabilization of vitamin C [14] . Therefore, metal complexes based on DHA that have exceptional chelating properties were synthesized and fully discussed in the literature [15] . Remarkably, it is acknowledged that the nickel complexes with Ni 2+ are more active than either dehydroacetic acid ligand alone or the pure inorganic nickel salts [16] . As example, nickel(II) complexes of symmetrical and unsymmetrical Schiff bases containing dehydroacetic acid have been prepared and characterized by S.F. Tan and K.P. Ang [17] . Recently, the synthesis of α,β-unsaturated compounds derived from DHA and indol-3-carboxaldehyde has been particularly focused [18, 19] . Also, F.A. El-Saied et al. have synthesized nickel(II) complexes, derived from DHA, N4-dialkyl and 3azacyclothiosemicarbazones [20] . In the context of this current study and as a continuation of our previous works [21] , we report herein the synthesis and determination of the crystal structure of the nickel(II) complex [Ni(DHA) 2 (DMF) 2 ] as mononuclear complex. The structural and spectral data of the prepared complex are presented with its magnetic properties which are also well discussed. As for the red/ox behavior of this nickel complex, it was investigated by cyclic voltammetry. Furthermore, this study is supported by some appropriate quantum parameters quantified by using the density functional theory (DFT). Hirschfeld surface analysis and HOMO/LUMO energy were also used to forecast the structural and spectroscopic parameters of this synthesized compound. Finally, the ADMET/tox studies and the binding affinity, of the [Ni(DHA) 2 complex with the key target main protease (M pro ) of SARS-CoV-2 through molecular docking analysis. In this study, the chemicals and solvents were used as received without any additional purification and were bought from Aldrich and Fluka. FT-IR spectra were recorded on Perkin and Elmer 1000-FT-IR Spectrometer using KBr disks, while the electronic spectra (UV-Vis.) were obtained on a Unicam UV-300 Spectrophotometer having 1 cm as a path length cell. Electrochemical experiments were performed on a voltalab 40, Potentiostat galvanostat controlled by microcomputer. Thus, the cyclic voltammograms were recorded using an individual cell of 5 ml using a conventional three-electrode system in distilled H 2 O + NaCl (10 −1 M) and NaOH (2·10 −1 M) under nitrogen atmosphere. The working electrode was a disc of glassy carbon (diam. 3 mm) while the counter electrode was a platinum wire, and all potentials were expressed versus the saturated solutions of calomel electrode (SCE). The synthesis of nickel complex was carried out in the same experimental conditions, early applied for copper complexes, as described in the literature [21˗22] . The complex was obtained by mixing stoichiometric amounts of dehydroacetic acid 0.168 g (1mmol) with nickel(II) chloride hexahydrate (NiCl 2 ·6H 2 O) 0.238 g (1 mmol). To this mixture, 0.06 g (1 mmol) of ethylenediamine, dissolved in 10 ml of methanol, were slowly added. Then, an excess of pyridine 0.158 g (2 mmol) was dissolved in 5 ml of methanol. After two hours of reaction, a sky-blue precipitate was observed and recovered by filtration. It was then purified by successive washings using cold methanol and diethyl ether followed by drying over CaCl 2 to give a sky-blue powder. This blue powder was dissolved in hot N,N-dimethylformamide in which the suitable monocrystals of [Ni(DHA) 2 (DMF) 2 ] were formed by slow evaporation in 10 days. The blue prism crystals of the nickel complex were obtained after redissolution of the resulting sky-blue powder followed by slow evaporation in hot N,N-dimethylformamide. X-ray data were collected with a Bruker Apex II CCD area detector, equipped with diffractometer using a graphite monochromated Mo-Kα radiation source (λ = 0.71073 Å) at 150 K. The crystallographic data and experimental details for structural analysis are summarized in Table 1 . The described structure was resolved through direct methods with SIR2004 [23] to locate all the non-H atoms, which were anisotropically refined with SHELXL97 [24] and using full-matrix least squares on squared procedure from within the WINGX suite of software [25] used to prepare material for publication. All the H atoms were placed in the calculated positions and constrained to ride on their parent atoms. Drawings of molecules were produced with the program ORTEP-3 and Diamond [26, 27] . Computational study of [Ni(DHA) 2 (DMF) 2 ] complex were performed by Gaussian 09 along with Gauss View 5.0.8 for visualizations by using DFT method at the level of B3LYP/6-31G (d, p) basis set [28] . This basis set was one of the most acknowledged basis sets used in the study of medium and large sized hydrogen bonded systems, and also for yielding a very small Basis Set Superposition Error (BSSE) [29] . The geometries of this complex were fully optimized without any constraint on every bond length, bond angle and dihedral angle were compared with the crystalline structure. The optimized structures were confirmed to be local minima by performing harmonic vibration frequency analyzes. Further analyzes of electronic and optical properties were carried out by calculation of UV-Vis. absorption in the framework of Time-Dependent DFT (TD-DFT) [30, 31] . Under the same basic set, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) were as well achieved. The reactivity descriptors that include energy gap (ΔEgap), hardness (η), softness (S), global electronegativity (χ), electrophilicity (ω), molecular electrostatic potential (MEP), Mulliken atomic charge and NLO have also been computed. Molecular docking and ADME/Tox studies were performed using Molecular Operating Environment (MOE) version 2019.10 molecular modelling software and Discovery Studio (DS) 2.0 Software (Studio 2.5, Accelrys, Co. Ltd., San Diego, CA, USA), respectively. The protein data bank of the Research Collaboration for Structural Bioinformatics (RCSB), website (http://www.rcsb.org/pdb) was used to obtain the crystal structure the Main Protease (M pro ), also called 3Clike protease (3CLpro), which represents a pivotal role in the propagation of SARS-CoV-2. The crystal structure of this enzyme was downloaded in complex with its inhibitor N3, (PDB ID: 6LU7; 2.16 Å) [3] . The targets protein, Mpro (6LU7), used in this study, were protonated where hydrogen atoms were added with its 3D geometry, corrected for any found errors in the connection or type of different atoms, and then energy minimized at the end of the preparation steps. This step followed by site-finder to define and isolate the same binding pocket of the co-crystallized native inhibitor as dummy atoms over helix for the docking step. The 2D structure of the [Ni(DHA) 2 (DMF) 2 ] complex was drawn using Marvin software [32] , converted to 3D and optimized by HyperChem 8.03 software [33] . Then, the energy minimizing of this compound was done under the following experimental conditions like temperature = 300 K and pH = 7. Furthermore, the geometry was performed using the field strengths in the MMFF94x implanted in MOE and Austin model 1 (AM 1) with gradient value of 0.0001 kcal/mol. Then, the compound was saved in mdb format as a new database. After that the protein structures and ligands were minimized and the active binding site was located, the ligands were docked to the protein using the dock tool of Molecular Operating Environment (MOE) 2019.0102 software for fitting of the isolated [Ni(DHA) 2 (DMF) 2 ] complex database into the active site of SARS-CoV-2 Mpro enzyme. Docking site was selected as dummy atoms, alpha triangle as the placement methodology, and London dG as scoring methodology. The docking process was run followed by evaluation of poses. Poses with the highest energy scores and best ligand-enzyme interactions were selected and recorded. Filters us Molecular weight of the ligand (<500 Da), high lipophilicity (Log P<5), number of hydrogen bonds donors (<5), number of hydrogen bond acceptors (<10) and lip druglike were used to investigate the monocrystal's drug-likeliness using the Lipinski's Rule of Five [34] . Parameter details were calculated using Molecular Operating Environment (MOE) 2019.0102 software. The [Ni(DHA) 2 (DMF) 2 ] complex were screened using Discovery Studio (DS) 2.0 Software (Studio 2.5, Accelrys, Co. Ltd., San Diego, CA, USA) to predict their important pharmacokinetic properties. The prediction of pharmacokinetics ADME (Absorption, Distribution, Metabolism, Excretion and Hepatotoxicity) properties include eight mathematical models, such as solubility and its level, HIA (human intestinal absorption), ADME AlogP98, ADME PSA-2D, BBB penetration level PPB (plasma protein binding), ADME CYP2D6 and hepatotoxicity. To determine the therapeutic compatibility of the drug, another method was used, toxicity prediction, by using the USFDA (US FDA, United States Food and Drug Administration) standard toxicity risk predictor software TOPKAT in the Discovery Studio. In this module, Ames mutagenicity, TOPKAT mouse female NTP, TOPKAT mouse male NTP probability and Daphnia EC 50 were predicted via Discovery Studio (DS) 2.0 Software (Studio 2.5, Accelrys, Co. Ltd., San Diego, CA, USA). The complex was found air stable and insoluble in most of the common organic solvents but soluble in hot DMF and DMSO. The given formulation of this complex is based on crystal structure, infrared and UVvisible spectra, magnetic measurements, electrochemical analysis and theoretical studies. The structure of the synthesized complex is represented in the following Scheme 1. Ethylenediamine (1 mmol) and an excess of Ni Scheme 1. Synthetic pathway for preparation of the Ni(II) complex with 3-acetyl-2-hydroxy-6-methyl-4Hpyran-4-one. The good quality crystals of this nickel complex was crystallized by dissolving [Ni(DHA) 2 The crystallization time was 15 days at room temperature. This nickel complex crystallizes in a triclinic space group P-1 system with a = 7.7441 (10) Å, b = 8.2173 (10) Å, 9.5338 (12) Å, α = 85.152 (7) °, β = 85.953 (9) °, γ = 78.157 (9) ° of volume V = 590.77 (13) Å 3 and Z = 1. An ORTEP view of the complex with the selective atom numbering scheme is shown in Fig. 1 . A summary of the details about crystal data, collection parameters and refinement are documented in Table 1 . The structure of the nickel complex reveals that it crystallizes with two molecules of the crystallization solvent (DMF). As shown in Fig. 1 , each nickel(II) ion centre is six coordinated in a form of a distorted square-bipyramid environment. This octahedral coordination sphere is formed by the both deprotonated dehydroacetic acid ligands through their two ketone oxygen atoms and two phenoxy oxygen atoms and also coordinated by two DMF ligands. The bond distances and angles are summarized in Table 2 . The bond lengths involving the phenoxy oxygen (C4-O1), the phenyl carbon (C2-O2) and ketone carbon of DMF molecule (C9-O5) are 1.276 (4), 1.248 (4) and 1.325 (4) Å, respectively. These values were found to be in the range of those usually observed as standard of C(sp 2 )-OH and C(sp 2 )-O-bond length [35, 36] . The six Ni-O bond lengths are 1.974 (2) and 2.003 (2) Å for DHA molecules and 2.108 (2) Å ( Table 2) . These values are in good agreement with related reported complexes [37, 38] . The selected cis angles of the nickel environments deviate from their ideal values only for 90° angles and found to be in the ranges 87.65 (8)°-92.35 (9)° while, for trans angles, they were obtained with ideal values and equal to 180°. The crystal packing for [Ni(DHA) 2 (DMF) 2 ] complex, can be described as an alternating layers parallel to (011) plane (See Fig. 2 ) and molecules are also linked via C-H…O intermolecular hydrogen bonds forming an infinite three-dimensional network. Geometrical parameters of [Ni(DHA) 2 (DMF) 2 ] complex were determined with SC-XRD and B3LYP/6-31G(d,p) and collected in Table 2 been done in the gas phase and experimental measurement has been carried out in the solid phase. Hence, the assessment of these outcomes confirmed that most of the simulated values are barely larger than the experimental values. In present work, most of the computed bond angles are absolutely larger than the experimental values due to the low scattering factor of hydrogen atoms participating in the X-ray diffraction experiment producing a huge deviation from the theoretical C−H bond lengths [39] . In this case, it can concluded that there is an acceptable agreement between the DFT calculated geometries and corresponding The solid-state infrared spectrum has been recorded in the region of 4000˗400 cm −1 . The corresponding theoretical spectrum has been calculated by using B3LYP level and 6-31G(d,p) basis set. The experimental and calculated vibrational spectra are given in Fig. 4 . The experimental and theoretical vibrational wave numbers, FT-IR intensity and assignments of vibration modes have also been presented in Table 3 . As shown in Fig. 4 , the formation of this nickel complex has been confirmed by detecting the characteristic absorption band of lactone O-C=O group, observed at 1710 cm −1 [40] . The dehydroacetic acid containing as well C=O group would be presented as strong and sharp absorption band at 1690 cm −1 [41] while, those calculated were found at 1702 and 1734 cm −1 , respectively. In addition, the absorption band observed at 1352 cm −1 is assigned to the C-N stretching vibration, which is in good agreement with the calculated wave stretching vibration frequencies [43] while those calculated were detected in the region 2944-3275 cm −1 . As presented in Table 3 and and (B) Experimental spectrum in KBr. TD-DFT calculations of [Ni(DHA) 2 (DMF) 2 ] complex in gas phase were performed by using B3LYP-6-31G (d,p). The electronic absorption wavelengths, oscillator strengths (ƒ os ) and excitation energy transitions of this complex have been given in Table 3 . Experimental and calculated spectra of the complex are given in oscillator strength of 0.0021 at 366 nm, so this absorption can be assigned to n-π* corresponding to the transitions of C=O chromophore [44] . This complex, showed also a weak broad band at 420 nm. The energy absorption band at 420 nm is 0.5624 eV with oscillator strength of 0.0004 at 504 nm. This absorption can be corresponding to the ligand-metal charge transfer transition (LMCT) [45] . It was observed that the results for hybrid potential B3LYP are slightly augmented and the finest values were found for B3LYP functional. The highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) are frequently studied in order to impart key information regarding the electron-donor and electron-acceptor character of the complexes which shall lead to the interpretation of the charge transfer process [46] . The lower energy of the HOMO indicates the lower ability as an electrondonor, resulting from the higher energy of LUMO and higher resistance to accept an electron. This allows elucidation of chemical stability by observing the difference in energy between HOMO and LUMO (Egap) while the large E gap is preferred for the high stability of complex. The distribution of electron density of the [Ni(DHA) 2 (DMF) 2 ] complex is illustrated in the following Fig. 6 . In addition, this study can also explain the chemical concept of chemical softness and hardness. With small E gap values, the complex may be considered as a "soft" base due to its high energy of HOMO and therefore enhance the interaction with the LUMO of soft acids. Other than that, indices such as electron affinity and ionization potential are also commonly interconnected with the studies of HOMO and LUMO energies in pursuing a better grasp of how complex theoretically behave, chemical wise. The exchange-correlation functional: B3LYP was checked by comparing the results of the reactivity descriptors parameters data for [Ni(DHA) 2 (DMF) 2 ] complex as presented in Table 4 . The given low energy gap (2.43929 eV) thus indicates high reactivity of the complex due to ease of charge transfer process [47] . The high value of softness (1.85776) or the low value of hardness (1.21964) indicates lower energy is needed for electron transition from HOMO to LUMO which means that the complex is susceptible to deform and ready to interact with other nucleophilic active site such as amino acid. The Mulliken atomic charges of [Ni(DHA) 2 (DMF) 2 ] complex were calculated by DFT using B3LYP 6-31G (d,p) at a basis set, the data were tabulated in Table S2 . The Mulliken's atomic charges (Fig. 7) of the DFT calculation revealed a charge distribution in individual atoms. The charges on carbon atoms exhibited either positive or negative values. All hydrogen atoms displayed a neat positive charge and act as acceptor atoms. All oxygen atoms of the optimized compounds were shown to have a negative charge, which acts as donor atoms. The positively charged centers are the most susceptible sites for nucleophilic attacks i.e., electron donation. However, the most negatively charged centers are the most susceptible sites for electrophilic one [48] . The molecular electrostatic potential (MEP) provides a visual method to understand the relative polarity of a molecule, electrostatic potential (electron + nuclei) distribution [49] and to find the reactive sites of molecules [50] . The molecular electrostatic potential of the studied complex is calculated using the B3LYP/6-31G (d,p) level and seen in Fig. 8 . The total electron density mapped with electrostatic potential surface of the ligands are shown in Fig. 8A , The MEP mapped surface and electrostatic potential contour map for positive and negative potentials are also shown in Fig. 8B and 8C. As represented in this figure, the maximum negative region is the preferred sites for electrophilic attack, indicated as red color. So, an attacking electrophile will be attracted by the negatively charged sites, and the opposite situation for the blue regions. These negative and positive sites help to predict the regions where a compound is responsible of non-covalent interactions [51] . complex. NLO properties such as molecular static polarizability (α), anisotropy of polarizability (∆α) and first hyperpolarizability (ß) have been calculated by quantum chemical calculations with B3LYP/6-31G (d, p). The investigation of NLO character explores the non-centrosymmetric nature and architecture of the molecules. This can be achieved by positioning of appropriate donor-acceptor units and alteration of πconjugation length [52˗54]. The <α>, ∆α and ß tot parameters can be obtained by using following equations: The Fig. S2 shows in the forward scan two oxidation waves at Epa1 = +770 mV/ SCE and Epa 2 = +1256 mV/s. The first wave represents the oxidation of the nickel(II) species [55] according to the following equation (1): While the second wave can be attributed to the oxidation of the ligand. On the back sweep, one cathodic reduction wave is only observed and appears at Epc1 = -743 mV / SCE. This wave represents the reduction of the Ni(II) species into Ni(I) [56] as shown by the following electrochemical reaction: When the electrode potential is scanned from 0.0 to +0.80 V, a reduction couple is discerned in the reverse scan, the half-wave potential located at E 1/2 = +753.5 mV/SCE and peak-to-peak separation estimated to about ΔE =180 mV can be assigned to the Ni(III)/Ni(II) redox couple [57] . The effect of the scan rate on the voltammetric response of the synthesized complex was also analyzed. Fig. S2 shows that the Ni(III)/Ni(II) redox couple seems to be dependent on the scan rates, especially for the anodic wave. This anodic wave, at low rates (5-100 mV s −1 ), splits in two waves Epa(I) and Epa(II) becoming at higher scan rate 100 mV s −1 a unique wave. Here, it was observed that as the scan increases, both waves Epa(I) and Epc(II) shift to the more anodic and cathodic potentials, respectively. An increase in peak currents was as well noted [58] . However, it can also be observed that the current varies linearly with the root square of the scan rate indicating a diffusion redox process. The The magnetic behavior of the sample was studied in the temperature range varied from 5 to 300 K ( Table 5) . Molecular docking is a technique which explains how a molecule interacts with each other when it placed in a restricted system under a well-defined force field. It is widely used to predict and identify ligands that fit into the binding pocket of a protein of interest. The energy of interaction of the protein-drug complex, i.e., the binding energy is reflected in the docking score [62] . In this paper, the molecular docking is performed using the Molecular Operating Environment (MOE) 2019.0102 software to evaluate the binding mode of ligand and interactions in the active site. Our screening was performed against the key target of the main protease (M pro ) to identify a therapeutic agent of COVID-19. Table 6 demonstrates the binding energies and different bonds of our complex. As compared to the inhibitor N3 (-7.44 kcal/mol), the [Ni(DHA) 2 (DMF) 2 ] complex exhibited a promising energy scores (-7.13) kcal/mol. The active Mpro homodimer of SARS-CoV-2, according to the 3D structure, comprises two protomers, constituting three domains. Domain I contains 8-101 Amino acid residues, Domain II contains 102-184 Amino acid residues and amino acid residues 201-306 constitute Domain III. We have defined the Mpro binding site residues as Thr24, Thr25, Thr26, Leu27, His41, Met49, Tyr54, Phe140, Leu141, Asn142, Gly143, Ser144, Cys145, His163, His164, Met165, Glu166, Leu167, Pro168, His172, Asp187, Arg188, Gln189, Thr190, Ala191, and Gln192 residues. Thr24, Thr26 and Asn119 are predicted to play important roles in drug interactions [63, 64] . His163, His172 and Glu166 are believed to provide the opening gate for the substrate in the active state of the protomer [65] , His41 and Cys145 form the catalytic dyad in the active site [66] and His164 is essential for the enzyme activity. Table 6 During drug development, pharmacokinetics studies evaluate the concentration of drugs and their actions in the body. However, it demonstrates the movement of drug into the body after the administration of drug. The study Drug-like properties and pharmacokinetic properties is essential to assess the safety and efficacy of drugs and design appropriate regimens Lipinski's rule of 5, Rules for the drug-likeliness, has been studied to ensure the drug-likeliness of This Rule is able to predict absorption or permeation of a potential drug candidate combining specific parameters [4] . According to this, good oral bioavailability is more likely when there are Molecular weight (MW) not greater than 500 g/mol, a partition coefficient log P less than or equal five number of hydrogen bond donors (HBD) (NH and OH groups) not greater than five, and number of hydrogen bond acceptors (HBA) (O and N atoms) not exceeds ten (10) . Drugs having log P ranging from 0 to 5, have high possibility of oral absorption. In general, an orally active drug has no more than one violation of these criteria. Violation of these parameters will reduce the draggability of the molecule [67] . As reported in Table 7 , the tested complex followed Lipinski's rules for drug-likeness with no more than one violation. [Ni(DHA) 2 (DMF) 2 ] complex drug ligand exhibited two violation regarding hydrogen acceptor and molecular weight. Its molecular weight is slightly more than the recommended values but, within the permissible limit. [Ni(DHA) 2 (DMF) 2 ] complex showed the value of Log P less than five, which indicates that the molecules are more likely to be in the hydrophilic environment and are favorable for their drug-likeness. The most popular and authentic rule for confirming the drug-likeness of the ligand is Lipinski's rule. However, many drugs have enormous pharmacological properties but, do not pass through the Lipinski's filter. They have been agreed by the FDA as possible drug for clinical purposes [68] . Despite violation of some rules, approved anticancer and anti-infective drugs from natural products or their semi-synthetic derivatives such as taxol and amphotericin B have also some violations but, they are biologically effective as drugs. Therefore, these results do not interfere with the development of this complex SARS-CoV-2 therapeutic agent [69] . software to predict their important pharmacokinetic properties. ADME properties include absorption, distribution, metabolism, excretion and hepatotoxicity (ADME) are considered as important in order to determine their role in drug discovery and success of the development of compound for human therapeutic use [70] . The screened results of ADME were summarized in Table 8 Intestinal Absorption (HIA) of drug after oral administration was determined by obtained levels: 0 (good), 1 (moderate), 2 (low), 3 (very low). [Ni(DHA) 2 (DMF) 2 ] complex was found to has moderate HIA, therefore it has an adequate gastro intestinal absorption. Easy distribution of the drug through blood brain barrier (BBB) measured by its AlogP98 value which must be less than 5. The obtained absorption levels determine drug absorption, it decreases inversely with the level, i.e., level 0 denotes proper absorption, level 1 denotes moderate absorption and so on. The intestinal absorption and blood brain barrier penetration were predicted by developing an ADME model using descriptors 2D PSA and AlogP98 that include 95% and 99% confidence ellipses. These ellipses define regions where well-absorbed compounds are expected to be found. Compounds with more than 90 % absorbability would lie within the 99 % confidence ellipse as the green ellipse Compounds falling out of the ellipse are those with less than 30 % absorbability. As shown in Fig. 13 , the results of DS-ADME model screening showed that [Ni(DHA) 2 (DMF) 2 ] complex within the 99% adsorption ellipse possess 99% confidence levels for human intestinal absorption. Table 8 shows that this complex has undefined values for BBB penetration level. Hence, the complex may not be able to penetrate the blood brain barrier. Hence, neurotoxicity of CNS with the proposed complex is not expected. ADME plasma protein binding model predicts whether a compound is likely to be highly bound to carrier proteins in the blood. True symbolizes binding and false symbolizes non-binding. The ADME plasma protein binding property prediction clearly suggests that the complex is not likely to be highly bound to carrier proteins in the blood. Metabolism is predicted based on the Cytochrome P450 2D6 inhibition. This descriptor determines inhibitory effect by the predicted classes: non-inhibitor (false) and inhibitor (true). [Ni(DHA) 2 (DMF) 2 ] complex is predicted as non-inhibitors of CYP2D6 suggesting that this complex is well metabolized in Phase-I metabolism. Hence, the side effects (i.e., liver dysfunction) are not expected upon administration of this complex. The toxicity predictions of the compounds were also investigated with Discovery Studio using the TOPKAT protocol. TOPKAT is a computational tool for prediction of toxicity of compounds developed by Accelrys [71] . It employs robust and cross-validated Quantitative Structure Toxicity Relationship (QSTR) models to predict a probable value of toxicity. The toxicity profiles calculated for all the compounds are tabulated in Table 9 . The results indicate that [Ni(DHA) 2 (DMF) 2 ] complex is non-Mutagen, non-carcinogen with Mouse female and male NTP. Therefore, the toxicity of the ligands was found to be suitable for the development into a medical drug. Thus, they experiment significant first-pass effect. Since December 2019, the COVID-19 consequences have severely affected the social life of human and also impacted the global economy. The growth of safer and extensive available strategies to overcome the increasing number of cases is greatly predicted. In this paper, we note the following results:  A novel single crystal nickel(II)-complex has been synthesized and thoroughly characterized by conventional spectroscopic techniques, such as FT-IR, UV-Vis and cyclic voltammetry.  The structural characterization supported by XRD-diffraction analysis of a single crystal has been used to define the molecular structure of the synthesized nickel complex.  The electrochemical behavior of the nickel complex has also been investigated. The linear relationship observed between the anodic peak current (ip a ) and the square root of the scan rate (v 1/2 ) indicates that the electrochemical process is governed by diffusion.  The geometry was optimized using DFT approach based on B3LYP/6-31G(d) allowing to explain the frontier molecular orbitals in gas phase. Meanwhile, TD-DFT methods were explored to better understanding the electronic transitions in the UV-Vis. spectrum.  Obtained dipole moment, polarizability and hyperpolarizability values indicate that this nickel complex can be considered as an effective NLO materials.  Finally, the molecular surfaces showed the electrophilic and nucleophilic regions on the nickel complex. Subsequently, this complex has been screened for its inhibitory effect against the main protease targets of SARS-CoV-2. A remarkable binding energy value, non-carcinogic and non-mutagenic properties have been shown.  This study suggests that the screened complex could be further explored for structural modification and detailed investigations in order to synthesize as newer potent agents endowed with better anti SARS-CoV-2 therapeutic inhibition.  As last conclusion, the diversification of this kind of studies will be considered as an efficient tool to master this disease. Validation, Reviewing and Editing. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. WHO coronavirus disease (COVID-19) dashboard. WHO. table The molecular biology of coronaviruses Structure of M pro from SARS-CoV-2 and discovery of its inhibitors An overview of severe acute respiratory syndrome-coronavirus (SARS-CoV) 3CL protease inhibitors: peptidomimetics and small molecule chemotherapy Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors Zur Synthese der Dehydracetsäure aus Acetessigester (Synthesis of dehydracetic acid from acetoacetic ester) Oxovanadium(IV) and ruthenium(II) carbonyl complexes of ONS-donor ligands derived from dehydroacetic acid and dithiocarbazate: Synthesis, characterization, antioxidant activity, DNA binding and in vitro cytotoxicity Design, synthesis, characterization, antimicrobial evaluation and molecular modeling studies of some dehydroacetic acid-chalcone-1,2,3-triazole hybrids Synthesis and antibacterial activities of enamine derivatives of dehydroacetic acid Bicyclic 6-6 Systems: Two Heteroatoms 1:1, Comprehensive Heterocyclic Chemistry Biological Role of Chalcones in Medicinal Chemistry. Vector-Borne Diseases: Recent Developments in Epidemiology and Control Structures and biological activity of cinnamoyl derivatives of coumarins and dehydroacetic acid and their boron difluoride complexes Simultaneous analysis of dehydroacetic acid, benzoic acid, sorbic acid and salicylic acid in cosmetic products by solid-phase extraction and highperformance liquid chromatography New Developments and Novel Therapeutic Perspectives for Vitamin C Electronic structure and vibrational frequencies in dehydroacetic acid (DHA) transition-metal complexes: A DFT study Synthesis and structure of cadmium and zinc complexes of dehydroacetic acid Nickel(I1) and palladium(II) chelates of dehydroacetic acid Schiff bases derived from thiosemicarbazide and hydrazinecarbodithioate Nickel(II) and Zinc(II) complexes based on DHA: Synthesis, X-ray crystal structure, antibacterial activity and DFT computational studies Unsymmetrical bis-Schiff bases derived from ethylenediamine, dehydroacetic acid and another aldehyde/ketone and the preparation and characterisation of their copper(II) and nickel(II) complexes Synthesis, Xray crystal structures and catecholase activity investigation of new chalcone ligands Ni(II) complexes of chalcone derivatives: Synthesis, X-ray crystal structure, electrochemical properties and DFT computational studies nickel(II), and copper(II) complexes of dehydroacetic acid N4-dialkyl-and 3-azacyclothiosemicarbazones Acta Cryst A novel copper(ii) complex with an unsymmetrical tridentate-Schiff base: synthesis, crystal structure, electrochemical, morphological and electrocatalytic behaviors toward electroreduction of alkyl and aryl halides Bis[μ-(3-acetyl-2-hydroxy-6-methyl-4H-pyran-4-one-κ3O:O′,O′′)]diaquatetrakis(pyridine-κN)dicopper(II) diperchlorate SIR2004: an improved tool for crystal structure determination and refinement A short history of SHELX WinGX and ORTEP for windows: an update ORTEP-3 for windows: a version of ORTEP-III with a graphical user interface (GUI) Development of the Colic-Salvetti correlation-energy formula into a functional of the electron density On the Effect of the BSSE on Intermolecular Potential Energy Surfaces. Comparison of a priori and a posteriori BSSE Correction Schemes Time-dependent density-functional response theory for molecules Key concepts in time-dependent density-functional theory Molecular Modeling System) Hypercube, Inc., 1115 NW, 4th Street Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings New nickel(II) and palladium(II) complexes with unsymmetrical quadridentate Schiff bases derived from 8-amino-4-methyl-5-azaoct-3-en-2-one. Crystal and molecular structure of [4,10-dimethyl-3-(nitroso-χN)-5,9-(diaza-χ2N)-trideca-3,10-diene-2,12-(dionato-χO)(2-)]nickel(II) Synthesis, characterization and X-ray crystal structure of novel nickel Schiff base complexes and investigation of their catalytic activity in the electrocatalytic reduction of alkyl and aryl halides Tables of bond lengths determined by X-ray and neutron difraction. Part 1. Bond lengths in organic compounds Synthesis and characterization of nickel(II) and copper(II) complexes with tetradentate Schiff base ligands Quantum chemical calculation studies on 4-phenyl1-(propan-2-ylidene) thiosemicarbazide Synthesis, spectral, thermal and biological studies of transition metal complexes of 4-hydroxy-3-[3-(4-hydroxyphenyl)-acryloyl]-6-methyl-2H-pyran-2-one Bidentate coordinating behavior of chalcone based ligands towards oxocations: VO(IV) and Mo(V) Preparation, physical and chemical studies on metal complexes of Schif bases as a nucleus key to prepare nanometer oxides have catalytic applications: nickel(II) complexes derived from 4-aminoantipyrine derivatives Synthesis of metal complexes involving Schiff base ligand with methylenedioxy moiety: spectral, thermal, XRD and antimicrobial studies UV-vis, IR and 1H NMR spectroscopic studies of some Schiff bases derivatives of 4-aminoantipyrine Co(II) and Ni(II) complexes of 4-(2-thiazolylazo)resorcinol: syntheses, characterization, catalase-like activity, thermal and electrochemical behaviour Frontier Orbitals and Organic Chemical Reactions DFT studies of disubstituted diphenyldithiophosphates of nickel(II): structural and some spectral parameters Structural, optical, morphology characterization and DFT studies of nano sized Cu(II) complexes containing Schiff base using green synthesis Natural localized molecular orbitals The electrostatic molecular potential as a tool for the interpretation of molecular properties Synthesis, X-ray crystal structure, DFT calculation and catalytic activity of two new oxidovanadium(V) complexes containing ONO tridentate Schiff bases First Theoretical Framework of Triphenylamine-Dicyanovinylene-Based Nonlinear Optical Dyes: Structural Modification of π-Linkers Effect of π-conjugation spacer (CC) on the first hyperpolarizabilities of polymeric chain containing polyoxometalate cluster as a side-chain pendant: A DFT study A DFT Study on The Two-Dimensional Second-Order Nonlinear Optical (NLO) Response of Terpyridine-Substituted Hexamolybdates: Physical Insight on 2D Inorganic-Organic Hybrid Functional Materials, Eur Ni(II), and Zn(II) of multidentate ligands resulting from condensation of 2-pyridinecarboxaldehyde with alpha, omegatriamines: selective imidazolidine/hexahydropyrimidine ring opening revisited Schiff base complexes derived from meso-1,2-diphenyle-1,2-ethylenediamine: synthesis, characterization, crystal structures, electrochemistry, and catalytic studies Mn(II) and VO(IV) containing N,O-donor Schiff base ligand Theory and Application of Cyclic Voltammetry for Measurement of Electrode Reaction Kinetics A nickel(II) di-μ2-phenolato bridged dinuclear complex: Weak antiferromagnetic interactions in nickel(II) dimers Magnetic relationship in new Ni II -Cu II -Ni II trinuclear complexes according to nickel(II) geometry Transition metal complexes of ligands containing azomethine group Stereochemistry and reactivity of four-coordinate nickel(II) complexes derived from (S)-(+)-cysteine derivatives The Many Roles of Computation in Drug Discovery Potential inhibitors against 2019-nCoV coronavirus M protease from clinically approved medicines Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor Structural basis of SARS-CoV-2 3CL(pro) and anti-COVID-19 drug discovery from medicinal plants In search of SARS CoV-2 replication inhibitors: Virtual screening, molecular dynamics simulations and ADMET analysis Drugbank database site Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development Screening of cryptogamic secondary metabolites as putative inhibitors of SARS-CoV-2 main protease and ribosomal binding domain of spike glycoprotein by molecular docking and molecular dynamics approaches Current and future perspectives on the development, evaluation, and application of in silico approaches for predicting toxicity The authors thank the Algerian Ministry of Higher Education and Scientific Research (MESRS) and the Director General for Scientific Research and Technological Development (DGRSDT) for the financial support.