key: cord-0767060-5avhc3mb
authors: Topal, Tufan; Zorlu, Yunus; Karapınar, Nazan
title: Synthesis, X-ray crystal structure, IR and Raman spectroscopic analysis, quantum chemical computational and molecular docking studies on hydrazone-pyridine compound: as an insight into the inhibitor capacity of Main Protease of SARS-CoV2
date: 2021-04-21
journal: J Mol Struct
DOI: 10.1016/j.molstruc.2021.130514
sha: db98b03af6ffae93630932d68922e4eae416e19b
doc_id: 767060
cord_uid: 5avhc3mb

The characterization and synthesis of 3-chloro-2-{(2E)-2-[1-(4-chlorophenyl)ethylidene]hydrazinyl}pyridine (CCPEHP) was investigated in our study. Mass and UV-visible spectra were recorded in chloroform solvent. The CCPEHP molecule containing pyridine and chlorophenyl rings and hydrazone group crystallized in the triclinic system and P-1 space group. FTRaman and FTIR spectra were performed in the solid state. The optimized geometry of CCPEHP was computed by DFT/B3LYP method with 6-311G (d, p) and 6-31G (d, p) levels. The computed vibrational analysis, electronic absorption spectrum, electronic properties, molecular electrostatic potential, natural bond orbitals analysis and other calculated structural parameters were determined by using the DFT/B3LYP/6-31G (d, p) basis set. The correlation of fundamental modes of the compound and the complete vibrational assignments analysis were studied. The strong and weak contacts were identified by using Hirshfeld surface analysis. The molecular modeling results showed that CCPEHP structure strongly binds to COVID-19 main protease by relative binding affinity of -6.4 kcal/mol.

cancer drugs and cerivastatin as HMG-CoA enzyme inhibitor drugs [2, 3] . The importance of pyridine and phenol vibration spectra in biological and industrial systems attracts great attention from a spectroscopic point of view [4] . In addition, hydrazone-containing compounds show electrophilic or nucleophilic character thanks to the azomethine groups they contain, which makes them unique in drug development studies [5] . While forming the molecular system, atoms do not bond randomly, they interact using a certain order and rules according to the laws of nature. The best experimental technique used to analyze these seemingly complex structures is the X-ray diffraction method [6] . Ab inito methods are of great importance in examining the properties of molecules, which are very important for researchers [7] . Such calculations are based on quantum theory and all chemical properties of the material can be calculated. Again, thanks to this method, the properties of the drugs designed for the treatment of various fatal diseases are determined, especially in pharmacology, and the loss of money and time can be prevented by aiming to inhibit the target protein [8] . Computer-aided drug design has been a promising technique for designing new and powerful therapeutics. Many theoretical studies have been applied to identify inhibitors against COVID-19. , created a three-dimensional SARS-CoV2 M pro homology model based on SARS-CoV M pro and screened it through protein modeling and virtualizing against 1903 drug inhibitors, highlighting nelfinavir as a potential drug inhibitor against SARS-CoV2 M pro [9, 10] . There is currently no defined treatment or antiviral drug for COVID-19 treatment. However, the US Food and Drug Administration (FDA) has approved both Chloroquine (CQ) and Hydroxychloroquine (HCQ) for COVID-19 control and treatment for emergency purposes [11] . Chloroquine and Hydroxychloroquine are notorious for toxicokinetic properties and can reveal dose-dependent toxicity. Recently, multi-centered global clinical trials are being conducted to evaluate the therapeutic potential of Chloroquine and its hydroxyl form, Hydroxychloroquine as a treatment for this infection [12] [13] [14] . The synthesized CCPEHP compound has a similar molecular structure with (CQ) and (HCQ) drugs. In our study, the behaviors of CCPEHP compound and (CQ) and (HCQ) drugs in molecular docking studies have been tried to be correlated. CCPEHP molecule was designed by different substituting (pyridine, hydrazone and heterocyclic) groups. These groups are considered for their potential to interact with amino acids in 4 different types (hydrophobic, hydrogen bonds, ionic and water bridges) [15, 16] . More specifically, the characterization of the synthesized compound with X-ray diffraction, spectroscopic and packing analysis and with quantum chemistry calculations was evaluated in this paper. The crystalline reactive sites with MEP analyses, electronic properties with HOMO-LUMO gap energies, donor-acceptor stabilization energies E (2) with NBO approach, electronic transitions with UV-Vis spectrum analysis, the effects of the compound on crystal packaging and the definition of non-covalent forces with Hirshfeld surface analysis were investigated. Finally, the inhibitor-receptor 2D and 3D interactions of the CCPEHP compound were visualized. [17] [18] [19] .

Intensity data for CCPEHP were collected using the ω-θ2 scan technique on an Agilent SuperNova, (Single source at offset and Eos CCD detector) diffractometer with a SuperNova (Mo) X-ray Source (Mo-Kα, k = 0.71073 Å). The CRYSALISPRO software program was used for cell refinement, data collection and reduction [20] . The FTIR spectrum of the CCPEHP compound was obtained by the Perkin Elmer FTIR-Spectrometer Spectrum Two Model between 400-4000 cm -1 . Mass spectrum (ESI) was recorded on TSQ Fortis™ Triple Quadrupole Mass Spectrometer LC-MS/MS model device. The FTRaman spectrum was obtained in the region 100-4000 cm -1 with a Thermo Fisher Scientific model DXR™ 3xi Raman Imaging Microscope with 1064 nm Nd: YAG laser source with a resolution of 3 cm -1 . The melting point of the crystallized material was determined by the Stuart SMP10 device. The absorption spectrum was recorded using a Shimadzu UV-1800 spectrophotometer between 240-800 nm. 3-chloro-2-hydrazinopyridine, 4′chloroacetophenone, ethanol, acetonitrile were bought from Sigma Aldrich company.

3-chloro-2-{(2E)-2-[1-(4-chlorophenyl)ethylidene]hydrazinyl}pyridine (CCPEHP) was prepared according to the literature [21] . The 4′-chloroacetophenone (1 mmol, 0.1545 g) in 15 ml acetonitrile was mixed with 3-chloro-2-hydrazinopyridine (1 mmol, 0.1435 g) in 10 ml acetonitrile (Fig.1) . The mixture was slowly stirred for 24 hours at RT. After stirring, the solution was left at RT for 3 days and crem product was obtained and washed with water and recrystallized by acetonitrile, and dried over in a vacuum desicator. C 13 

The geometries of 3-chloro-2-{(2E)-2-[1-(4-chlorophenyl)ethylidene]hydrazinyl}pyridine was optimized with Gaussian 09W package program using Beck's triply parameter hybrid Exchange function B3 along with Lee-Yang-Parr (LYP) functions with 6-311G (d, p) and 6-31G(d, p) basis sets [22, 23] . Diffuse functions on all atoms support these basis sets. Vibrational frequencies and IR intensity values of CCPEHP compound optimized with DFT/B3LYP/6-31G (d, p) level were calculated. In addition, potential energy distributions of vibration modes were determined using the VEDA4 program [24] . UV spectrum of the CCPEHP compound in chloroform solvent was computed using the TD-DFT method via DFT/B3LYP/6-31G (d, p) level [25] . The LUMO and HOMO gap energies, molecular electrostatic potential (MEP) and natural bonding orbitals (NBO) were calculated using the B3LYP method and 6-31G (d, p) level [1, 26, 27] . The Gauss-view 6.0 visualization software and the Gaussian 09w package program were used for all the computerized calculations [28] [29] [30] . Gamess US software was used as an auxiliary theoretical calculation program [31] . CrystalExplorer package program generated Hirsfeld surfaces and 2D fingerprint plots which were used to determine the analysis of the intermolecular interactions in the crystal case [32] . 

Targeting the 3CL pro protease of the corona virus consisting of structural and nonstructural polyproteins could constitute a valid approach for the treatment of COVID-19 potential drug design. One of the structures of 3CL pro like protease protein was downloaded (PDBID: 6LU7, 2.16 Å) from Protein Data Bank (PDB) (http://www.rcsb.org), in 3D format [33] . Preparation for simulation and simulation processes of COVID-19 3CL pro /M pro structure and ligand were carried out by AutoDock Tools 1.5.6, MG Tools of AutoDock Vina program and Discovery Studio 2020 Client (Dassault Systemes BIOVIA) [34, 35] . In protein preparation, the missing residues were first checked and there were no missing residues of the viral protein. Chain A, one of two chains of main protein, was selected and the binding affinity of this chain was performed with our potential drug candidate molecule. Water molecules and hetero atoms were deleted from the COVID-19 main protease, then Kollman charges and polar hydrogen atoms were added to the target protein (chain A) and saved with AutoDock Tools in pdbqt format [10, 36] . Autodock vina algorithm with the most appropriate configurations (exhaustiveness= 8, binding modes= 9, energy difference= 3 kcal/mol, feasible with x, y, and z coordinate) parameters is used in simulation studies [34, 37] . Default settings were used in the other parameters. The grid box value and intermolecular interactions of the ligand-protein complex were visualized with the Discovery studio 2020 client program. Hydrophobic gap coordinates of the binding site were used in the docking of N3 (as control) and fungal metabolites. The active site of the amino acid was calculated as sized (30 Å x 30 Å x 30 Å) and centered (-10.71 Å x 12.41 Å x 68.83 Å). (Default grid spacing= 0.375 Å) [38] . In order to prepare the ligand for docking, the CIF data obtained from the X-ray crystal structure [39] was converted to mol2 format. The partial charge and torsion angels were changed and saved in pdbqt format with AutoDock Tools. Before the simulation, the separate pdbqt formats and docking parameters of the drug candidate compound and COVID-19 main protease were combined into a single .txt file and the complexation process was completed through the AutoDock Vina software.

3.1. X-ray data collection and structure refinement

The structure was solved using SHELXT [40] and then refined by full-matrix least-squares refinements on F 2 using the SHELXL [41] in Olex2 Software Package [42] . Crystallographic data and refinement details of the data collection for CCPEHP are tabulated in Table 1 . Geometrical calculations and crystal structure validations were performed using Platon software [43] . Mercury software [44] was used for visualization of the cif file. Additional crystallographic data with CCDC reference number 2068456 for CCPEHP has been deposited within the Cambridge Crystallographic Data Center via www.ccdc.cam.ac.uk/deposit. In order to get a better insight into the intermolecular interactions in the solid state, the clear colorless crystals of CCPEHP for single crystal X-ray diffraction were grown by slow evaporation, and structural analysis exhibits that CCPEHP crystallizes in the triclinic space group P-1, with the asymmetric unit containing two CCPEHP molecules. The structure exhibited chlorophenyl subunit with the chloropyridine moieties lying nearly in the same plane with the tilt angles 2. 

The optimized structure of the CCPEHP compound, namely C 13 H 11 Cl 2 N 3 , was calculated with the help of quantum theory [30] . The selected DFT optimization and metrical parameters from Xray diffraction are tabulated in Table 2 . It is understood from the results given in Table 2 that, experimental and theoretical data are in good agreement with each other. Although the results are compatible with each other, there are still slight differences. The reason for this difference should be kept in mind that the experimental results were obtained in the crystal phase, while the theoretical results obtained with 6-31G (d, p) and 6-311G (d, p) levels were obtained in the gas phase [46] . The bond length C(1)-N(8) in the molecule was observed as 1.323 Å and the same bond length was calculated as 1.335 and 1.331 Å for 6-31G (d, p) and 6-311G (d, p) basis sets [28] . Table 2 The single-crystal X-ray and optimized selected geometrical parameters of CCPEHP. 

Vibrational frequencies calculations were obtained using optimized structure parameters. 81 fundamental vibrational modes of CCPEHP were calculated based on quantum chemical calculations. DFT potential symmetrically makes vibration wavenumbers higher than experimental data. Inconsistencies in these wavenumbers are corrected by calculating non-harmonic corrections or corrected by presenting a scaled area or by an appropriate factor directly calculated. In order to eliminate these systematic errors, the vibration wavenumbers were scaled with a 0.961 scaling factor. [28] http://cccbdb.nist.gov/vsfx.asp. After scaling the wavenumbers, the observed and the theoretical vibration wavenumbers showed a very good agreement with each other. The basic vibration modes of CCPEHP compound which are characteristic N-H, C-Cl, C-H, C-C, C=N vibrations found in pyridine and chlorophenyl rings and hydrazone group were tried to be investigated. The scaled and unscaled calculated vibrational modes and intensities of FTIR with their PED and observed FTIR and FTRaman bands are reported in Table 3 . It is observed in Table 3 that the scaled FTIR bands of the CCPEHP have a close relationship between experimental frequencies [47] . The CCPEHP molecule has 29 atoms. All vibration modes are assigned thanks to the graphical interface GaussView 6.0 software program [29] . The pyridine ring of CCPEHP involves 1 N atom and 5 C atoms which generate 4 C-C and 2 C=N stretching vibrations [48] . Additionally, there are also C=N peaks belonging to the hydrazone group in this structure. The C=N vibration peaks are recorded at 1030-1590 cm -1 in FTRaman and 1030-1600 cm -1 in FTIR in the literature [49, 50] . The stretching vibration of C=N groups were observed at 1030, 1239, 1589 cm -1 in FTIR and 1031, 1257, 1590 cm -1 in FTRaman and scaled calculated bands at 1045, 1216, 1264, 1319, 1437, 1504, 1554, 1601 cm -1 for CCPEHP molecule. The C-H stretching vibrations in compounds containing aromatic groups are generally in the range 3000-3100 cm -1 [51, 52] . The absorption peaks at 3069 and 3025, 3055 cm -1 were assigned to the FTRaman and FTIR spectra. The bands attributed to 3039, 3046, 3079, 3086, 3088, 3096, 3100 and 3103 cm -1 were calculated as C-H stretching vibrations in chlorophenyl and pyridine rings. In literature, the C-H peaks in-plane bending vibrations of aromatic rings are observed at 1020-1300 cm -1 [28] . The C-H in-plane bending vibrations were calculated at 1057, 1069, 1096, 1107, 1165, 1216, 1278 cm -1 , respectively based on PED results. Generally, the stretching vibrations of C-C bands occur in the range 1430-1625 cm -1 in the aromatic ring [53] . This band was observed at 1506, 1557 cm -1 in FTIR and 1509, 1560, 1582 cm -1 in FTRaman and that it was well coherent with the simulated data at 1504, 1548, 1554, 1577 and 1580 cm -1 . The peaks calculated at 395, 392, 468, 620, 639, 711, 764, 886, 985, 995 cm -1 were attributed for C-C-C in-plane bending vibrations of CCPEHP and the peaks at 78, 211, 265, 349, 404, 488, 536, 697, 916, 955 cm -1 were attributed for C-C-C-C torsion vibrations [54] . Since compound CCPEHP has more than one chlorine atom, it can occur in symmetric and asymmetric modes [55] . The C-Cl peaks were observed in both FTRaman and FTIR at 427 and 440 cm -1 , respectively [56] . The stretching vibrations of C-Cl mode were allotted at 392, 711 and 995 cm -1 with significant PED (>10%) [1] . Symmetrical stretching vibrations of methyl (-CH 3 ) groups are generally found in between 2850-2940 cm -1 [57] . The vibrations of methyl group were observed at 2910 in FT-Raman and 2948 cm -1 in FTIR. Similarly, computed values of 2912 and 2961 cm -1 were attributed to C-H stretching vibrations in the CH 3 group. The N-H stretching vibrations of aliphatic amines are found between 3300-3500 cm -1 . In this study, this vibration was computed at 3431 cm -1 and observed at 3383 cm -1 in FTRaman and 3380 cm -1 in FTIR. PED column is contributing exactly 100% to this mode [3] . The combined spectra and correlation diagram of FTIR data between theoretical and experimental frequencies are given in Fig. 3a and 3b . The linear and the equality in the correlation diagram are y = 0.9922x + 6.9917 (R² = 0.9996). (-15) γ; out-of-plane bending, ν; stretching, δ; in-plane-bending, τ; torsional, PED; Potential Energy Distribution. Fig. 3 . The a) combined spectra, b) correlation diagram for the observed and calculated FTIR wavenumbers (with scaling factor) of CCPEHP compound.

The electronic absorption spectrum of CCPEHP compound was computed by TD-DFT with B3LYP/6-31G (d, p) level in chloroform solution. The oscillator strengths, excitation energies and wavelengths were estimated by using this method and basis set. The CCPEHP compound was studied at room temperature (RT) in 2 x 10 -5 mol L -1 chloroform solution and in the wavelength range of 250-450 nm. Ultraviolet-visible spectroscopy of the compound was compared with the electron transition possibilities. Theoretical calculations of the CCPEHP and corresponding experimental UV-Vis spectra are shown in Fig 4a and 4b . The observed spectrum demonstrated the maximum absorption (A max ) peaks at 292 and 331 nm which indicated that π _ π* and n _ π* transitions, respectively. According to the theoretical calculations, the bands at 296 and 334 nm were attributed to the π-π* and n-π* transitions. These transitions of oscillator strengths (f) 0.1551/0.9191 and excitation energies (eV) 4.1801/3.7101 were obtained from simulation results and these data are given in Table 4 . The π-π* transitions of the CCPEHP arise from the electrons of the pyridine ring, while the n-π* band is due to hydrazone groups of atoms [2, 58] . 

Hirshfeld surface analysis emerged from an attempt to define the area covered by the molecule forming the crystal lattice. It is an effective tool that provides more information about the crystal structure by color-coding the short or long contacts in the molecule system [59] . Crystal packing differences occur between the particular types of intermolecular interactions due to molecular abilities. The formation of intermolecular interactions between the donor and acceptor affects biological activity [60] . The Hirshfeld surface analysis and 2D fingerprint plots studies were carried out using CIF data [61] . CCPEHP compound mapped over d norm , shape index and curvedness forms and surfaces are shown in Fig 5. The volumetric and area surface properties of mapped CCPEHP were found as 325.08 Å 3 and 311.95 Å 2 , respectively [58] . Short and strong contacts on the d norm surface are attributed to red spot color, while long and weak contacts are identified with blue areas in Fig 5a. According to the d norm surface analysis, it is understood that weak intermolecular interactions are dominant in the CCPEHP molecule. The regions observed in shape index (the orange-red colored) and in curvedness (the large green straight surface, which is separated by a blue outline) represent π-π stacking interactions of the CCPEHP molecule (Fig 5b and 5c) . Furthermore, these defined regions indicated that intermolecular interactions with neighboring molecules [62] . The 2D fingerprints plots are associated with Hirshfeld Surface area and the interactions of H‧‧‧H, Cl‧‧‧H/H‧‧‧Cl and C‧‧‧H/H‧‧‧C were the most contributing close contact interactions to the total surface and are illustrated in Fig. 6a, 6b and 6c . The intermolecular interaction areas are shown in blue, and areas outside of this interaction are shown in grey; H‧‧‧H 31.2%, Cl‧‧‧H/H‧‧‧Cl 20.7%, C‧‧‧H/H‧‧‧C 20.3%, N‧‧‧H/H‧‧‧N 10.4%, Cl‧‧‧Cl 7.1%, C‧‧‧Cl/Cl‧‧‧C 4.0%, C‧‧‧C 3.8%, C‧‧‧N/N‧‧‧C 1.8%, N‧‧‧Cl/Cl‧‧‧N 0.6% were obtained from the 2D fingerprint plot study [17] . When the results of the surface analysis are evaluated, the stability of the CCPEHP compound is due to the presence of hydrogen bonds in the structure. Additionally, examining the physical and chemical properties of potential biologically active compounds in simulation studies allows us to predict their behavior in in-vitro studies. 

Natural bonding orbital gives very important theoretical information about conjugation and hyperconjugation and intra and intermolecular charge transfer of the molecular orbital system [2] . Donor-acceptor interaction indicates that delocalization of electron density between formally unoccupied non-Lewis type natural bonding orbitals and occupied Lewis type natural bonding orbitals. The more intense the interaction between the electron acceptor and donor, the higher the E (2) value [63] . Decolalizations effects of electrons have been obtained from the second order perturbation results performed by the NBO 3.1 package [27] . DFT/B3LYP method with 6-31G (d, p) basis set was used to understand and interpret different intermolecular interactions of atoms in molecules. For CCPEHP, the some strongest π to π* interactions were calculated as π(C1-N8)→π*(C4-C5) 28.20 kcal/mol, π(C2-C3)→π*(C1-N8) 24.70 kcal/mol, π(C2-C3)→π*(C4-C5) 13.10 kcal/mol, π(C4-C5)→π*(C2-C3) 27.10 kcal/mol, π(C17-C19)→π*(C18-C20) 18.80 kcal/mol, π(C17-C19)→π*(C22-C24) 21.90 kcal/mol, π(C18-C20)→π*(C17-C19) 18.60 kcal/mol, π(C18-C20)→π*(C22-C24) 21.00 kcal/mol, π(C22-C24)→π*(C17-C19) 18.10 kcal/mol, π(C22-C24)→π*(C18-C20) 17.80 kcal/mol. Regarding the lone pairs of electrons delocalization, LP 1 (N 9 ) to π*(C1-N8) and LP 1 (N 9 ) to π*(N11-C12) with the hyperconjugative stabilization energy data were obtained as 49.10 kcal/mol and 30.80 kcal/mol, respectively. The other strongest hyperconjugative interactions were calculated as LP 1 (N 8 )→π*(C1-C2) 11.15 kcal/mol and LP 1 (N 11 )→π*(C12-C13) 11.60 kcal/mol [64] . According to the results of NBO analysis, stabilization energy values were also obtained between the antibonding donor orbitals and the antibonding receptor orbitals.

Determining the electronic properties of the molecule is a very important tool in interpreting the results of biological activity [65] . In the molecular system, frontier molecular orbitals (FMOs) are of great importance for drug designing. The highest occupied molecular orbital (HOMO) corresponded to ionization potential, while the lowest unoccupied molecular orbital (LUMO) corresponded to electron affinity [66] . Thanks to the FMOs energy gap, important information about the chemical stability and reactivity of molecules can be obtained. The chemical stability of the compound plays a crucial role in drug design investigations. A large energy gap brings about to the molecule less reactive and more stable while smaller energy gab brings about to less stable and more reactive in the system [2] . The electronic properties of the CCPEHP and the theoretical energies of (E HOMO ) and (E LUMO ) orbitals of the CCPEHP compound were investigated (Fig 7) . The energy difference (ΔE) was found as 4.1720 eV in (LUMO-HOMO) and 6.0178 eV in (LUMO+1)-(HOMO-1). According to the result obtained from the energy gap values, the CCPEHP molecule has high kinetic stability. The reactivity descriptors parameters are calculated from HOMO and LUMO orbitals where electron affinity, chemical softness, chemical hardness, electronic chemical potential, ionization potential are tabulated in Table 5 [67, 68] . The chemical stability of the molecule was supported by chemical hardness and softness in addition to the energy gap [18] . In Fig. 7 , while the red color symbolizes the positive lobes, the green color symbolizes the negative lobes of the molecule. Generally, HOMO and LUMO charges are localized over the whole molecule [1] . 

Molecular electrostatic potential gives information about the partial charges and electronegativities of different atoms in the molecule [26] . The color scale for (MEP) map surface: red zones are for negative potential, electron rich, strong repulsion; blue zones for positive potential, electron lacking, strong attraction; green zones for neutral. The map of 2D contour and electrostatic potential for CCPEHP are shown in Fig. 8a and 8b . The MEP values of the CCPEHP molecule were predicted in the range of -6.880 e -2 and 6.880 e -2 . In CCPEHP crystal, the red surfaces of the negative potentials are spread over the pyridine ring and hydrazone group of nitrogen atoms. As a result of this, red zones (N8, N9 and N11) which have negative electrostatic potential leads to electrophilic attack and this indicates the strongest repulsion. In Fig. 8a , the 2D contour map of the molecule's reactants clearly confirms the proton transfer [18, 28, 69] . 

One of the best ways to understand ligand and protein binding interactions in the pharmaceutical design industry is molecular docking studies. Molecular docking is one of the best simulation methods applied to estimate drug-target interactions by calculating energy minimization and binding energy. It identifies the best position against the target protein. Autodock vina algorithm determines the most suitable ligand conformation with minimized energy with the help of some scoring functions [37, 70] . After the simulation process, The lowest binding energy and the most hydrogen binding pose were selected from the ViewDock window. Ligand_1-6lu7 was chosen as the best position is most associated with hydrogen binding. Hydrogen binding and other interactions between the ligand and the 6lu7 protein were examined in the Biovia Discovery Studio (DS) visualizer (Fig. 9 ). The highest binding affinity was determined with COVID-19 M pro , which exhibits covalent bond interaction thanks to molecular dynamic simulation [71] . The residues from the N3 binding site are Ser46, His41, Met49, Phe140, Tyr54, Leu141, Gly143, Asn142, Cys145, His164, His163, Met165, Leu167, Glu166, Pro168, Phe185, Asp187, His172, Gln189, Thr190, Ala191 and Gln192 [72] . When the molecular docking result applied in our study was analyzed, it was observed that the ligand binds to similar enzyme sites. CCPEHP molecule has shown five interactions with COVID-19 M pro [38] (Fig 10a) . The CCPEHP ligand was observed to have interacted with two conventional hydrogen bonds His164 (2.31 Å) and His41 (2.92 Å), a Pi-Pi T-Shaped type with His41 (4.65 Å), a Pi-Sulfur with Cys145 (5.76 Å) and a Pi-alkyl with Met165 (4.91 Å). Pi-Pi T-Shaped and Pi-alkyl type interactions are included in the hydrophobic category. The binding affinity of this complex was found to be -6.4 kcal/mol. When the simulated binding interactions of the molecule were evaluated, the binding generally occurred in regions suitable for electrophilic attacks. While Pi interactions occurred on aromatic pyridine and chlorophenyl rings, hydrogen bond interactions occurred in the hydrazone region. In biological processes, the presence of hydrogen bonds is of great importance in the interaction of the inhibitor with the receptor. As shown in Fig 10b, the green areas represent the electron acceptor, while the pink areas represent the electron donor. It was revealed that the theoretically calculated inhibitor-receptor complex interaction gave a very successful result and in this sense, the CCPEHP compound we designed has a good inhibitory capacity [18, 19] . The behavior of molecular docking studies of FDA-approved drugs in the literature is examined. In a study by Dossary et al. (2021), the docking pose analysis showed that the chloroquine (CQ) bonded to Leu141, Asn142, His163, Pro168, Cys145 and Glu166 binding residues in the 6LU7 protein [73] . In another study on the main SARS-CoV-2 protease carried out by Nobre Júnior et al. (2020), the CCPEHP compound was showing better docking scores (-6.4 kcal/mol) as well as better ligandreceptor interactions in comparison to hydroxychloroquine and chloroquine (-5.0, and -4.7 respectively) [66] . 

3-chloro-2-{(2E)-2-[1-(4-chlorophenyl)ethylidene]hydrazinyl}pyridine was synthesized and characterized the experimental technique by single-crystal X-ray crystallography and on the other hand, quantum theory to support the experimental data. The geometrical parameters, vibrational frequencies and UV-visible spectrum were obtained from the optimized structure of the CCPEHP compound and compared with the experimental results. The experimental values of vibrational frequencies were in good agreement with the scaled calculated values. Hirshfeld Surface analysis proved that important contacts were the H‧‧‧H 31.2%, Cl‧‧‧H/H‧‧‧Cl 20.7% and C‧‧‧H/H‧‧‧C 20.3% interactions. Stabilization energy values between donor-acceptor orbitals and hyperconjugation interactions were determined by NBO orbital analysis. The reactive zones of CCPEHP were determined by molecular electrostatic potential (MEP) results in the nitrogens N8, N9 and N11 were nucleophilic zones. Frontier molecular orbitals (FMOs) and properties related to molecular interaction ability were predicted. By examining almost all the properties of the CCPEHP molecule, its capacity in in-silico studies has been tried to be revealed. The results of molecular docking, spectroscopic and theoretical techniques provide better information about inhibitor-receptor interactions. The rational design of drugs developed against dangerous diseases and understanding the mechanisms of these drugs at the molecular level provide valuable information on drug-receptor interactions. From molecular docking analysis, CCPEHP compound was fitted snuggly into the sites of COVID-19 M pro and exhibited Pi-Pi T-Shaped, Pi-alkyl type and hydrogen bond interactions. CCPEHP successfully docked with the amino acid molecule of the main protease of SARS-CoV2 with a high negative binding affinity. In conclusion, we hope that the successful insilico results we obtained in the CCPEHP molecule will be confirmed by in vitro and in vivo studies.

CRediT authorship contribution statement T. Topal: designed, did the in silico and theoretical calculation studies, Writing, supervised the project and edited the manuscript. N. Karapınar: did the experimental work and spectroscopic analysis. Y. Zorlu: Methodology, Visualization, Writing, did the X-ray analysis and refinement.

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.

A complete synergy on the experimental and theoretical study of the pyridine derivatives -2-Hydroxy-5-Nitropyridine and 2-Chloro-5-Nitropyridine

Quantum chemical studies and spectroscopic investigations on 2-amino-3-methyl-5-nitropyridine by density functional theory

Synthesis, spectroscopic characterization and quantum chemical computational studies of (S)-N-benzyl-1-phenyl-5-(pyridin-2-yl)-pent-4-yn-2-amine

Vibrational spectroscopic ( FT-IR , FT-Raman and NMR ) and DFT analysis of 2-methoxy-3-( trifluoromethyl ) pyridine

Biological activities of hydrazone derivatives

tert-bütil) benziliden) amino)-1,5-dimetil-2-fenil-1,2-dihidro-3H-pirazol-3-on'un X-ışını ve Teorik Çalışmaları

NMR, spectra, normal co-ordinate analysis, NBO, NLO and DFT calculation of N, N-diethyl-4-methylpiperazine-1-carboxamide molecule

Bis[2-(metilamino)troponato]bakır(II) Molekülünün Yapısal ve Elektronik Özelliklerinin Deneysel ve Kuramsal Analizleri

Nelfinavir was predicted to be a potential inhibitor of 2019-nCov main protease by an integrative approach combining homology modelling, molecular docking and binding free energy calculation

Prediction of potential inhibitors of the dimeric SARS-CoV2 main proteinase through the MM/GBSA approach

Does zinc supplementation enhance the clinical efficacy of chloroquine/hydroxychloroquine to win todays battle against COVID-19?

Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro

Chloroquine is a potent inhibitor of SARS coronavirus infection and spread

Screening of Chloroquine, Hydroxychloroquine and its derivatives for their binding affinity to multiple SARS-CoV-2 protein drug targets

Side-chain flexibility in protein-ligand binding: The minimal rotation hypothesis

In-silico strategies for probing chloroquine based inhibitors against SARS-CoV-2

Synthesis, crystal structure, TD/DFT calculations and Hirshfeld surface analysis of 1-(4-((Benzo)dioxol-5-ylmethyleneamino)phenyl)ethanone oxime

Synthesis, experimental, theoretical study and molecular docking of 1-ethylpiperazine-1,4-diium bis(nitrate)

Azafluorene derivatives as inhibitors of SARS CoV-2 RdRp: Synthesis, physicochemical, quantum chemical, modeling and molecular docking analysis

CrysAlisPro Software System, Version 1.171, Rigaku Corporation

Hirshfeld Surface Analysis, Drug-likeness Properties and Molecular Docking Studies of New Oxime-pyridine Compounds

Self-Consistent Equations Including Exchange and Correlation Effects

Vibrational energy distribution analysis (VEDA): scopes and limitations

Excited-state potential energy curves from time-dependent densityfunctional theory: A cross section of Formaldehyde's 1A1 manifold

Molecular electrostatic potentials: Concepts and Applications

-phenyl) boronic acid pinacol ester by DFT method

General atomic and molecular electronic structure system

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Discovery Studio Modeling Environment

Arabinosyltransferase C enzyme of Mycobacterium tuberculosis, a potential drug target: An insight from molecular docking study

Neuroprotective immunity by essential nutrient "Choline" for the prevention of SARS CoV2 infections: An in silico study by molecular dynamics approach

Potential Inhibitor of COVID-19 Main Protease ( M pro ) from Several Medicinal Plant Compounds by Molecular Docking Study

A novel 3-((5-methylpyridin-2-yl)amino)isobenzofuran-1(3H)-one: Molecular structure describe, X-ray diffractions and DFT calculations, antioxidant activity, DNA binding and molecular docking studies

SHELXT -Integrated space-group and crystal-structure determination

Crystal structure refinement with SHELXL

OLEX2 : a complete structure solution, refinement and analysis program

Structure validation in chemical crystallography

Mercury: Visualization and analysis of crystal structures

Synthesis of Dicyanovinyl-Substituted 1-(2-Pyridyl)pyrazoles: Design of a Fluorescent Chemosensor for Selective Recognition of Cyanide

Synthesis, crystal structure

Density functional theory studies on the structure, spectra (FT-IR, FT-Raman, and UV) and first order molecular hyperpolarizability of 2-hydroxy-3-methoxy-N-(2-chloro-benzyl)-benzaldehyde-imine: Comparison to experimental data

Investigation of Structures, FTIR, FT-Raman, In Vivo Anti-Inflammatory, Molecular Docking and Molecular Characteristics of 2-Amino-3-Pyridine Carboxaldehyde and Its Copper(II) Complex Using Experimental and Theoretical Approach

Synthesis and structural study on (1E,2E,1′E,2′E)-3,3′-bis

The spectroscopic (FT-IR, FT-Raman, dispersive Raman and NMR) study of ethyl-6-chloronicotinate molecule by combined density functional theory

New Type Biomembrane: Transport and Biodegradation of Reactive Textile Dye

Molecular structure analysis and spectroscopic characterization of 9-methoxy-2H-furo[3,2-g]chromen-2-one with experimental (FT-IR and FT-Raman) techniques and quantum chemical calculations

Vibrational spectra, molecular structure, natural bond orbital, first order hyperpolarizability, TD-DFT and thermodynamic analysis of 4-amino-3-hydroxy-1-naphthalenesulfonic acid by DFT approach

The vibrational spectra of pyridine, pyridine-4-d, pyridine-2,6-d2, and pyridine-3,5-d2

Some considerations on the out-of-plane vibration bands of PhnX type molecules in the 800-670 cm−1 region in relation to the estimation of the twist angle θ of benzente rings from their intensities

Infrared and Raman characteristic group frequencies : tables and charts

UV and quantum chemical studies on monomeric and dimeric conformations of 3,5-dimethyl-4-methoxybenzoic acid

Synthesis and characterization of zinc(II) complexes with new pyridine-based ligands: crystal structure, Hirshfeld surface analysis, and molecular docking study of lung cancer cell

Bonded-atom fragments for describing molecular charge densities

Synthesis, X-ray crystal structure, Hirshfeld surface analysis, and molecular docking study of novel inhibitor of hepatitis B: methyl 4-fluoro-3-(morpholinosulfonyl)benzo[b]thiophene-2-carboxylate

Hirshfeld surface analysis

Synthesis, crystallographic structure, Hirshfeld surface analysis and DFT calculations of two salen-type halogenated Schiff-base Ni(II) complexes

The spectroscopic (FT-IR, FT-IR gas phase, FT-Raman and UV) and NBO analysis of 4-Hydroxypiperidine by density functional method

NMR) and computational analysis (DOS, NBO, Fukui) of 3,5-dimethylisoxazole and 4-(chloromethyl)-3,5-dimethylisoxazole: A DFT study

New thiazole, pyridine and pyrazole derivatives as antioxidant candidates: synthesis, DFT calculations and molecular docking study

Quantum mechanical, spectroscopic and docking studies of 2-Amino-3-bromo-5-nitropyridine by Density Functional Method

Applications of the conceptual density functional theory indices to organic chemistry reactivity

Absolute electronegativity and hardness: applications to organic chemistry

Structural characterization, quantum chemical calculations and Hirshfeld surface analysis of a new third order harmonic organic crystal: 2-Amino-4-methylpyridinium benzilate

Challenges and advances in computational docking: 2009 in review

Reckoning a fungal metabolite, Pyranonigrin A as a potential Main protease (Mpro) inhibitor of novel SARS-CoV-2 virus identified using docking and molecular dynamics simulation

Virtual screening based on molecular docking of possible inhibitors of Covid-19 main protease

DFT and molecular docking study of chloroquine derivatives as antiviral to coronavirus COVID-19

I thank lecturer Süleyman Çelik for his contributions. The calculations reported in this article were completely/partially carried out at TUBITAK ULAKBIM, Grid Computing Center and High Performance (TRUBA resources).