key: cord-0833582-63bxiydw
authors: Abu-Melha, Sraa; Edrees, Mastoura M.; Riyadh, Sayed M.; Abdelaziz, Mohamad R.; Elfiky, Abdo A.; Gomha, Sobhi M.
title: Clean Grinding Technique: A Facile Synthesis and In Silico Antiviral Activity of Hydrazones, Pyrazoles, and Pyrazines Bearing Thiazole Moiety against SARS-CoV-2 Main Protease (M(pro))
date: 2020-10-06
journal: Molecules
DOI: 10.3390/molecules25194565
sha: cbc621ac6a986f534d937a279086bef1c2fa0210
doc_id: 833582
cord_uid: 63bxiydw

A novel series of some hydrazones bearing thiazole moiety were generated via solvent-drop grinding of thiazole carbohydrazide 2 with various carbonyl compounds. Also, dehydrative-cyclocondensation of 2 with active methylene compounds or anhydrides gave the respective pyarzole or pyrazine derivatives. The structures of the newly synthesized compounds were established based on spectroscopic evidences and their alternative syntheses. Additionally, the anti-viral activity of all the products was tested against SARS-CoV-2 main protease (M(pro)) using molecular docking combined with molecular dynamics simulation (MDS). The average binding affinities of the compounds 3a, 3b, and 3c (−8.1 ± 0.33 kcal/mol, −8.0 ± 0.35 kcal/mol, and −8.2 ± 0.21 kcal/mol, respectively) are better than that of the positive control Nelfinavir (−6.9 ± 0.51 kcal/mol). This shows the possibility of these three compounds to effectively bind to SARS-CoV-2 Mpro and hence, contradict the virus lifecycle.

Clean chemistry is a multidisciplinary concept for the judicious use of chemical reactions and prevents its pollution. It involves the designing of chemical products by using environmentally safe procedures that reduce the generation of hazardous substances, as well as the cost [1] . Grinding methods are appealing eco-friendly techniques for efficient organic synthesis with many advantages, such as: mild reaction conditions, easy separation and purifications [2, 3] , high efficiency and selectivity [4] , and environmental acceptability [5] . In the meantime, hydrazone is an effective nucleus because hydrazine-containing compounds (R 1 R 2 C=NNH-R 3 ) offer a wide range of biological activities, and selectivity [4] , and environmental acceptability [5] . In the meantime, hydrazone is an effective nucleus because hydrazine-containing compounds (R1R2C=NNH-R3) offer a wide range of biological activities, including anticholinesterase [6] , antimycobacterial [7] , anti-cancer [8] , antibacterial, anti-inflammatory, anti-viral, antifungal, anti-tubercular, analgesic, muscle relaxants, and antihistamines [9] [10] [11] [12] . Hydrazones possessing an azometine -NHN=CH-proton constitute an important class of compounds for new drug development. Both nitrogen atoms of the hydrazone group are nucleophilic and the carbon atom of the hydrazone group has both electrophilic and nucleophilic properties (Figure 1 ). These structural fragments are principally responsible for the physical and chemical properties of hydrazones. Therefore, many researchers have synthesized these compounds as target structures and evaluated their biological activities. These observations have been guiding the development of new hydrazones that possess varied biological activities [13] [14] [15] . Likewise, pyrazoline derivatives widely occur in nature in the form of alkaloids, vitamins, animal cells, and pigments. The unique structural array and pronounced pharmacological activities [16, 17] of pyrazolines have made them attractive synthetic targets for many researchers. In continuation of our interest in green synthesis of heterocycles [2, [18] [19] [20] [21] [22] [23] , herein, we propose an easy, efficient, solid-state grinding method for the preparation of hydrazones. This method can also be extended to the reactions of hyrazide with active methylene and anhydride compounds.

SARS-CoV-2 represents the century challenge to human beings. It represents the primary agent for the novel, deadly, pandemic of COVID-19 worldwide [24] . The main viral protease (M pro ) of SARS-CoV-2 sheds drug designers due to its vital role in keeping the virus alive [25] . M pro is a crucial enzyme for the viral polypeptide processing into the different enzymes that maintain the life cycle of the SARS-CoV-2. Hence, blocking M pro activity is vital in combating COVID-19. Different compounds (natural and synthetic) are used as possible potential therapeutics against SARS-CoV-2 M pro . In this study, we tested some hydrazones, pyarzoles, and pyrazines derivatives synthesized using a clean grinding technique.

Molecular modeling represents an excellent method used to judge the chemical reactivity in silico [26] . The fast improvement in the last two decades makes molecular modeling the second hand for drug designers. It helps reduce time, money, and effort in the journey of drug design [27] . It was successful in helping find potent anti-viral agents against many viruses, like Human Immunodeficiency Virus (HIV) and Hepatitis C Virus (HCV) [28] [29] [30] . Additionally, computer-aided drug design is used to find the potency of drugs and natural products against SARS-CoV-2 proteins during the last few months [31] [32] [33] [34] . In this study, molecular docking combined with molecular dynamics simulation is used to test the different compounds against SARS-CoV-2 M pro .

An initial endeavor was to investigate the condensation reactions under solid-state conditions. Thus, grinding of 3-acetyl-7-amino-1-aryl- [1, 2, 4] triazolo [4,3-a]pyrimidin-5(1H)-one (1) [35] and 4-methyl-2-phenylthiazole -5-carbohydrazide (2) [36] , in the presence of a few drops of acetic acid at Likewise, pyrazoline derivatives widely occur in nature in the form of alkaloids, vitamins, animal cells, and pigments. The unique structural array and pronounced pharmacological activities [16, 17] of pyrazolines have made them attractive synthetic targets for many researchers. In continuation of our interest in green synthesis of heterocycles [2, [18] [19] [20] [21] [22] [23] , herein, we propose an easy, efficient, solid-state grinding method for the preparation of hydrazones. This method can also be extended to the reactions of hyrazide with active methylene and anhydride compounds. SARS-CoV-2 represents the century challenge to human beings. It represents the primary agent for the novel, deadly, pandemic of COVID-19 worldwide [24] . The main viral protease (M pro ) of SARS-CoV-2 sheds drug designers due to its vital role in keeping the virus alive [25] . M pro is a crucial enzyme for the viral polypeptide processing into the different enzymes that maintain the life cycle of the SARS-CoV-2. Hence, blocking M pro activity is vital in combating COVID-19. Different compounds (natural and synthetic) are used as possible potential therapeutics against SARS-CoV-2 M pro . In this study, we tested some hydrazones, pyarzoles, and pyrazines derivatives synthesized using a clean grinding technique.

Molecular modeling represents an excellent method used to judge the chemical reactivity in silico [26] . The fast improvement in the last two decades makes molecular modeling the second hand for drug designers. It helps reduce time, money, and effort in the journey of drug design [27] . It was successful in helping find potent anti-viral agents against many viruses, like Human Immunodeficiency Virus (HIV) and Hepatitis C Virus (HCV) [28] [29] [30] . Additionally, computer-aided drug design is used to find the potency of drugs and natural products against SARS-CoV-2 proteins during the last few months [31] [32] [33] [34] . In this study, molecular docking combined with molecular dynamics simulation is used to test the different compounds against SARS-CoV-2 M pro .

An initial endeavor was to investigate the condensation reactions under solid-state conditions. Thus, grinding of 3-acetyl-7-amino-1-aryl- [1, 2, 4] at room temperature, gave the condensation products N'-[1-(7-amino-5-oxo-1-aryl-1,5-dihydro- [1, 2, 4] triazolo[4,3-a]pyrimidin-3-yl)ethylidene]-4-methyl-2-phenylthiazole-5-carbohydrazides (3a-c) (Scheme 1).

The spectral data are in accordance with the isolated products 3a-c (see the Supplementary Figure S1 ). 1 H-NMR revealed two singlet signals at δ = (2.37-2.45) and (2.58-2.64) ppm due to (CH 3 -C=N) and

(CH 3 -thiazole) [34] , respectively. Another three singlet signals at δ = (5.30-5.36), (6.92-6.98), and (9.42-9.59) ppm were revealed, attributed to (CH-pyrimidine), (NH 2 ) [32] , and (NH-CO) [34] , respectively. The molecular ion peaks and microanalysis results for the synthesized products were also in accordance with the theoretical calculation.

To verify the assigned structure of the products 3a-c, we have synthesized these compounds by alternative methods. Thus, treatment of aminothiouracil (4) with 2-[2-(4-methyl-2-phenylthiazole-5-carbonyl) hydrazono]-N-arylpropanehydrazonoyl chlorides (5a-c) [37] in dioxane, in the presence of a catalytic amount of triethylamine, under thermal conditions, gave the authentic products identical in all respects (mp, mixed mp, IR) with the isolated products 3a-c via Smiles rearrangement [38] and elimination of H 2 S (Scheme 2).

Next, we investigated the derivatization of acylhydrazide group with thiazole moiety to generate the novel hydrazone templets. Thus, the reaction of 4-methyl-2-phenylthiazole-5-carbohydrazide (2) Molecules 2020, 25, 4565 4 of 13 with 2-phenyl-4-methyl-5-acetylthiazole (8) [39] or cyclic ketones 10a-c under solid-state conditions furnished the corresponding acylhydrazones 9 or 11a-c, respectively (Scheme 3). Microanalyses and spectral data (IR, 1 H-NMR, 13 To generalize the grinding technique, the cyclo-condensation reactions were examined. As a consequence, treatment of acid hydrazide 2 with active methylene compounds, such as 2,4-pentanedione (12), ethyl 3-oxobutanoate (14) , diethylmalonate (16) , and ethyl cyanoacetate (18) , under grinding conditions led to the formation of (3,5- (17), and 3-amino-1-(4-methyl-2-phenylthiazole-5-carbonyl)-1,4-dihydropyrazol-5-one (19) , respectively (Scheme 4). 1 H-NMR spectrum of compound 13 showed that the characteristic singlet at δ = 6.05 ppm corresponds to pyrazole-CH [40] . Also, the characteristic (CH 2 ) protons of pyrazole ring appeared as a singlet at δ = 5.25 ppm [41] to confirm the structure of compound 15. Furthermore, the presence of a singlet at δ = 5.94 ppm corresponding to the (NH 2 ) protons proved the formation of the compound 19.

In the IR spectra, all pyrazole and pyrazoline compounds displayed a characteristic band near 1690 cm −1 attributed to C=O bridge stretching vibration (see the Supplementary Figure S3 ). Additionally, the reaction of acid hydrazide 2 with maleic anhydride (20) or phthalic anhydride (22) resulted in ring-opening and ring closure to give the corresponding 1-(4-methyl-2-phenylthiazole-5carbonyl)-1,2-dihydropyridazine-3,6-dione (21) or 2-(4-methyl-2-phenyl-thiazole-5-carbonyl)-2,3dihydro phthalazine-1,4-dione (23), respectively (Scheme 5). Spectral data (IR, 1 H-NMR, 13 C-NMR, and MS) and elemental analyses are in accordance with the structures of the isolated products (see the Supplementary Figure S4 ). 

The structures are tested for their anti-viral activity against SARS-CoV-2 main protease (M pro ) using molecular docking combined with molecular dynamics simulation (MDS) for up to 100 nanoseconds (see the Supplementary Figures S5 and S6 ). Figure 2 shows the average binding affinities of the compounds calculated using AutoDock Vina [42] . The anti-viral Nelfinavir (red column) is tested as a positive control. The error bars represent the standard deviation (SD) from all the docking trials using the different conformations of SARS-CoV-2 M pro . As reflected from Figure 2 , the average binding affinities of the compounds 3a, 3b, and 3c (green columns) (−8.1 ± 0.33 kcal/mol, −8.0 ± 0.35 kcal/mol, and −8.2 ± 0.21 kcal/mol, respectively) are better (less) than that of the positive control Nelfinavir (−6.9 ± 0.51 kcal/mol). This shows the possibility of these three compounds to effectively bind to SARS-CoV-2 M pro and hence, contradict the virus lifecycle. Additionally, the compounds 9, 11a, 11b, 11c, 13, 15, 17, 19, 21, and 23 have average binding affinities comparable to Nelfinavir (−6.4 ± 0.09 kcal/mol for compound 13 down to −7.4 ± 0.13 kcal/mol for compound 23); hence, it may bind to SARS-CoV-2 M pro and stop the viral infection as well. To further analyze the binding mode, the docking complexes are analyzed by the aid of the Protein-Ligand Interaction Profiler (PLIP) software [43] . Table 1 shows the interaction profiles of the compounds 3a, 3b, 3c, 9, 11a, 11b, 11c, 13, 15, 17, 19, 21, and 23 docked into the active site residues of SARS-CoV-2 PL pro . The most reported types of interaction were the H-bonds and hydrophobic contacts. Noticeably, the most reported interacting residue in almost all complexes is E166. For the best three compounds (3a, 3b, and 3c), the interaction is almost identical with 7 H-bonds (N142, G143, S144, E166, D187, and Q189) and one hydrophobic contact (E166). These three compounds have the lowest binding energy values (-8.1 to -8.3 kcal/mol). For the other compounds, the number of formed interactions is fewer. This is reflected in the binding energy values (-6.4 to -7.3 kcal/mol). The interaction pattern for the best three thiazole compounds, 3a, 3b, and 3c, are depicted in Figure 3 . These three compounds can interact with the same binding site pocket of the SARS-CoV-2 M pro as reflected in the superposition (Figure 4 ) of the compounds after docking. The compounds 3a, 3b, and 3c proved their ability to bind firmly to the active site of the protein, and hence could be successful candidates against COVID-19. Experimental validation of the binding potency of the three compounds against SARS-CoV-2 M pro is suggested as a future recommendation. Further in-depth MDS analysis of the best three compounds is suggested as future work. Additionally, we suggest testing these compounds against other enzymes for SARS-CoV-2, such as the papain-like protease (PL pro ) and RNA-dependent RNA polymerase (RdRp). It may bind firmly to these targets and may help inactivate SARS-CoV-2. 

An electrothermal Gallenkamp apparatus was operated to measure the melting points for the newly synthesized compounds (Bibby Sci. Lim. Stone, Staffordshire, UK). Pye-Unicam SP300 

An electrothermal Gallenkamp apparatus was operated to measure the melting points for the newly synthesized compounds (Bibby Sci. Lim. Stone, Staffordshire, UK). Pye-Unicam SP300 instrument in potassium bromide discs was used to measure Infra-Red, IR spectra (Shimadzu, Tokyo, Japan). A Varian Mercury VXR-300 spectrometer (300 MHz for 1 H-NMR and 75 MHz for 13 C-NMR) (Varian, Inc., Karlsruhe, Germany) was manipulated to measure the 1 H NMR, and 13 C-NMR spectra and the chemical shifts were related to that of the solvent. GCMS-Q1000-EX Shimadzu spectrometers (Tokyo, Japan) were conducted to record the mass spectra of the samples on the ionizing voltage at 70 eV. Elemental analyses were carried out at the Microanalytical Centre of Cairo University, Giza, Egypt.

General Procedure for the Synthesis of Compounds 3a-c, 9, 11a-c, 16-19, 22, and 23 A mixture of 4-methyl-2-phenylthiazole-5-carbohydrazide (2) (0.233 g, 1 mmol) and appropriate carbonyl compounds (1 mmol) was ground in a mortar at room temperature in the presence of a few drops of acetic acid for 10-20 min (monitored through Thin Layer Chromatography, TLC). The reaction mixture was poured into water, filtered, and washed with water and ethanol. The crude products were recrystallized from the appropriate solvent. Alternative Synthesis of Compounds 3a-c:

Equimolecular amounts of aminothiouracil (4) (0.143 g, 1 mmol) and 2-[2-(4-methyl-2-phenyl thiazole-5-carbonyl)hydrazono]-N-arylpropanehydrazonoyl chlorides (5a-c) (1 mmol) were dissolved in dioxane (15 mL), containing a few drops of triethylamine. The reaction mixture was refluxed for 2 h and poured into ice/HCl. The precipitate formed was collected by filtration, dried, and recrystallized from appropriate solvent to give the respective products identical in all respects (mp and mixed mp) with 3a-c. 3-Amino-1-(4-methyl-2-phenylthiazole-5-carbonyl)-1,4-dihydropyrazol-5-one (19 

Nanoscale Molecular Dynamics (NAMD) and Visualizing Molecular Dynamics (VMD) software were used in the molecular dynamics study (MDS) using the standard CHARMM 36 forcefield, and the proteins were solvated with the TIP3P water model [44] [45] [46] [47] . The temperature of 310 K and the 154 mM NaCl concentration were adjusted during the simulation period to simulate the protein at its physiological conditions. The MDS calculations were done on the Bibliotheca Alexandrina supercomputer, Alexandria, Egypt. Chimera software was used to cluster the trajectories after MDS using default parameters [48] . Five representative structures of SARS-CoV-2 M pro were used in the docking study after clustering the trajectories of the MDS of the M pro (PDB ID: 6Y84).

AutoDock Vina 1.1.2 software was used to test the thiazole compounds against SARS-CoV-2 main protease different conformations, while Nelfinavir was used as a positive control. The ligands were prepared using the PyMOL 2.4 and AutoDock Tools 1.5.6 software [49] . Any missed hydrogen atoms were added, followed by charges addition before input files generation. The docking utilizes flexible active site (H41 and C145) and flexible ligand protocol. The AutoDock Vina grid box size used in the docking experiments is 40 × 40 × 40 Å 3 , while the box centers were set to be at the dyads H41 and C145 (−3.9, 19. 

In this context, a series of novel hydrazones (3a-c, 9, 11a-c), pyarzoles (13, 15, 17, 19) , and pyrazines (21, 23) were synthesized using a solvent-drop, eco-friendly, environmental grinding method. The structure of the newly prepared compounds was established based on both elemental analysis and spectroscopic data and by an alternative method wherever possible. Also, the anti-viral activity of all the products was tested against SARS-CoV-2 main protease (M pro ) using molecular docking combined with molecular dynamics simulation (MDS). The average binding affinities of the compounds 3a, 3b, and 3c were better than that of the positive control Nelfinavir. So, these compounds can bind effectively to SARS-CoV-2 M pro and hence, contradict the virus lifecycle.

Supplementary Materials: The following are available online, Figure S1 : 1 H-NMR, 13 C-NMR and Mass spectra of compound 3a, Figure S2 : 1 H-NMR, 13 C-NMR and IR spectra of compound 9, Figure S3 : 1 H-NMR and IR spectra of compound 15, Figure S4 : 1 H-NMR and IR spectra of compound 23, Figure S5 : Molecular Simulations, Figure S6 : Molecular Simulations (RMSF).

Innovations and green chemistry

Solvent-drop grinding method: Efficient synthesis, dpph radical scavenging and anti-diabetic activities of chalcones, bis-chalcones, azolines, and bis-azolines

A facile solvent-free one-pot three-component method for the synthesis of 2-amino-4H-pyrans and tetrahydro-4H-chromenes at ambient temperature

Mechanochemical solvent-free and catalyst-free one-pot synthesis of pyrano[2,3-d]pyrimidine-2,4(1H,3H)-diones with quantitative yields

A simple and efficient method for solvent-free iodination of hydroxylated aromatic aldehydes and ketones using iodine and iodic acid by grinding method

The design of novel 4,6-dimethoxyindole based hydrazide-hydrazones: Molecular modeling, synthesis and anticholinesterase activity

Synthesis and in vitro antimycobacterial potential of novel hydrazones of eugenol

Microwave and conventional synthesis of ester based dicationic pyridinium ionic liquids carrying hydrazone linkage: DNA binding, anticancer and docking studies

Synthesis, characterization and antibacterial activity of (E)-4-((3-methyl-4-(methylsulfonyl)pyridin-2-yl)methoxy)-N -(substituted-benzylidene)benzohydrazide derivatives

Abdel-Shafy, A. Efficient Synthesis and In-Silico study of some novel pyrido

Green synthesis, molecular docking and anticancer activity of novel 1,4-dihydropyridine-3,5-dicarbohydrazones under grind-stone chemistry

Synthesis and properties of hydrazones bearing amide, thioamide and amidine functions

Synthesis and spectral characterisation of hydrazone based 14-membered octaaza macrocyclic Ni (II) complexes

X-ray crystal structures, DNA binding, oxidative cleavage and antimicrobial studies of two Cu (II) hydrazone complees

Synthesis of novel substituted pyrazole-5-carbohydrazide hydrazone derivatives and discovery of a potent apoptosis inducer in A549 lung cancer cells

Synthesis and studies on antidepressant and anticonvulsant activities of some 3-(2-furyl)-pyrazoline derivatives

Synthesis and antidepressant activity of some 1,3,5-triphenyl-2-pyrazolines and 3-(2"-hydroxy naphthalen-1"-yl)-1,5-diphenyl-2-pyrazolines

]thiadiazines using chitosan as ecofriendly catalyst under microwave irradiation

Synthesis and anticancer activity of arylazothiazoles and 1,3,4-thiadiazoles using chitosan-grafted-poly(4-vinylpyridine) as a novel copolymer basic catalyst

Green synthesis, molecular docking and pharmacological evaluation of new triazolo-thiadiazepinylcoumarine derivatives as sedative-hypnotic scaffold

DABCO-catalyzed green synthesis of thiazole and 1,3-thiazine derivatives linked to benzofuran

A Facile green synthesis and anti-cancer activity of bis-arylhydrazononitriles

Novel 4-heteroaryl-antipyrines as DPP-IV Inhibitors

World Health Organization Situation Report 72

Structure of M pro from SARS-CoV-2 and discovery of its inhibitors

Molecular Modelling: Principles and Applications

Reply to a letter to the editor

Applications of computer-aided approaches in the development of hepatitis c antiviral agents

Design and Development of Some Viral Protease Inhibitors by QSAR and Molecular Modeling Studies, in Viral Proteases and Their Inhibitors

Theoretical Study on Modified Boceprevir Compounds as NS3 Protease Inhibitors

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective

Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study

Anti-HCV, nucleotide inhibitors, repurposing against COVID-19

Natural products may interfere with SARS-CoV-2 attachment to the host cell

A convenient synthesis of novel derivatives of pyrido

New 5-thiazolylcarbohydrazon-n-allyl-thiazolines. Synthesis, characterization and antioxidant activity

A facile access and evaluation of some novel thiazole and 1,3,4-thiadiazole derivatives incorporating thiazole moiety as potent anticancer agents

pyrimidinones via tandem Japp-Klingemann, Smiles rearrangement, and cyclization reactions

General two-step preparation of chalcones containing thiazole

A Green protocol for catalyst-free synthesis of pyrazole in glycerol-water solution

Synthesis and characterization of some substituted pyrazoles as analgesics and anti-inflammatory agents

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

PLIP: Fully automated protein-ligand interaction profiler

Solvated docking: Introducing water into the modelling of biomolecular complexes

Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K

Scalable molecular dynamics with NAMD

VMD: Visual molecular dynamics

UCSF Chimera-a visualization system for exploratory research and analysis

The PyMOL Molecular Graphics System

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution

The authors declare that there is no conflict of interest regarding the publication of this paper.