key: cord-0810124-aijgmkax
authors: Krishna Cheerala, Vijay Sai; Ghanta, Prasanth; Neelakantan, Sundaresan Chittor
title: Design, synthesis and in silico screening of benzoxazole–thiazolidinone hybrids as potential inhibitors of SARS-CoV-2 proteases
date: 2021-12-10
journal: RSC advances
DOI: 10.1039/d1ra07504g
sha: f0282e8b17447fceecab2707499289bfbd04aa19
doc_id: 810124
cord_uid: aijgmkax

Hybrid molecules in the recent years have gained significant importance in drug research as promising therapeutic agents. We report a novel combination of two such bioactive scaffolds (benzoxazole and 4-thiazolidinone B–T hybrids) as inhibitors of SARS-CoV-2. The study uses an in silico approach to identify the potential of B–T hybrids as possible inhibitors of the SARS-CoV-2 proteases. Molecular docking was employed to identify the interactions of B–T hybrids with the two proteases – 3CLp (the 3-chymotrypsin-like protease) and PLp (the papain-like protease). Docking results of the screened 81 hybrids indicated that BT10 and BT14 interacted with the catalytic dyad residue of 3CLp (Cys145) with the best binding energy. MD simulations revealed that BT10 formed stable interactions via 4 hydrogen bonds with the catalytic site residues of 3CLp. In the case of PLp, BT27 and MBT9 interacted with the catalytic triad residue of PLp (His272) with high binding energy. MD simulations demonstrated that the reference drug Tipranavir relocated to the thumb region of the protease whereas BT27 remained in the active site of PLp stabilized by 2 hydrogen bonds, while MBT9 relocated to the BL2 loop of the palm region. The MM-PBSA and interaction entropy (IE) analysis indicated that BT14 exhibited the best ΔG (of −6.83 kcal mol(−1)) with 3CLp, while BT27 exhibited the best ΔG (of −7.76 kcal mol(−1)) with PLp. A four-step synthetic procedure was employed to synthesize the B–T hybrids starting from ammonium thiocyanate. The short-listed compounds in the case of 3CLp were synthesized and characterized using IR, NMR, and HRMS spectroscopic techniques.

An unprecedented disease outbreak in recent times which has created havoc worldwide since December 2019, is still embroiled in controversy regarding its genesis. 1 The affected people exhibited severe symptoms of pneumonia, and a plethora of scientic investigations revealed that the cause was a virus. This new virus of zoonotic origin belonging to the Coronaviridae family infected millions of people globally as a pandemic and was found to share 79.5% viral genome sequence similarity with SARS-CoV 2,3 and hence was labelled as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). 4 As of November 1st 2021, there were 246 million cases of infection and 4.99 million fatalities across the globe (https:// covid19.who.int/) due to SARS-CoV-2. Previous to this, viruses belonging to the family of Coronaviridae have been observed to cause epidemic outbreaks. For example, the Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and the Middle East Respiratory Syndrome Coronavirus (MERS-CoV), have exacerbated as an epidemic, causing substantial morbidity and mortality in 2002 and 2012 respectively. 5 The viruses belonging to this family are single-stranded RNA (ssRNA) viruses with the largest RNA genome (in the range of 27 to 32 kb) and at least six open reading frames (ORFs). The genome of CoVs which is highly susceptible to mutation has caused the emergence of new CoV strains. 6 Natural selection also espouses the strains that have relatively higher infectivity rates thus causing severe outbreaks. During this pandemic, several strains of SARS-CoV-2 with an increased transmission rate have been identied viz. B.1.1.7, B.1.1.28, B.1.351, and B.1.617, which have 30 to 80% higher infectivity rates than the original Wuhan strain. 7, 8 With the mutations causing the newly discovered strains to offset the immunity provided by the vaccines as well as from previous infections, there is a concerted effort to forage molecules capable of inhibiting proteases which play a pivotal role in the entry and replication of the virus. 9 , 10 Based on the viral entry into the host cell, synthesis, and replication, several drug targets are known for inhibiting SARS-CoV-2. 11 Among these enzymes, two targets -3CLp (the 3-chymotrypsin-like protease) and PLp (the papain-like protease) have received a lot of attention, as they are involved in the posttranslational modication of the polyprotein. 12, 13 They achieve this by cleaving the polyprotein precursor into individual functional non-structural proteins (NSP's). 14 Polyproteins are single giant proteins that are transcribed by a single gene and are typically non-functional. 14 These proteases are essential for virus reproduction, and hence are focussed as prime targets for antiviral drug development.

The initial and perhaps one of the most signicant steps in the drug development process is the design and discovery of new drug candidates. The idea of hybrid molecules was conceived to overcome the several challenges encountered in the current drug development process that include low pharmacokinetic proles, multidrug resistance, drug-drug interactions, and less stability of single target-directed drugs. 15, 16 To combine the effects of each molecule into one single chemical entity, two or more different pharmacologically active molecules are fused into hybrid molecules with a suitable linker. The resultant molecules may have benecial synergistic effects in addition to the desired pharmacokinetic and pharmacodynamic proles. Recent studies have corroborated that the hybrid molecules are pharmacologically more efficacious than their predecessors. 16 For example, melatonin, a well-known antioxidant, when combined with an anti-Alzheimer drug "Tacrine" to create new hybrid molecules that target multiple enzymes, was found to be multi-enzyme inhibitor and emerged as a promising candidate for cognitive enhancement. 17 Benzoxazole derivatives have generated a lot of interest in recent years, as precursors of novel bioactive molecules. The biological activities exhibited by benzoxazoles are antiviral, antimicrobial, anti-inammatory, anti-cancer, antihistamine, melatonin receptor antagonism, 5-HT3 antagonist effect, inhibition of Rho-kinase, and anti-amyloid properties. 18 Further, benzoxazoles were explored as the inhibitors of the NTPase/ helicase of Hepatitis C Virus (HCV), 19 the reverse transcriptase of HIV, 20 and the NS3 helicase of Dengue Virus (DENV), 21 They are available as the central active moiety in a variety of marketed drug molecules, few examples of which are shown in Fig. 1 . Table 1 lists a few of the potent benzoxazole derivatives that target a wide range of enzymes associated with various disorders.

4-Thiazolidinone is one of the important bioactive vemembered rings and a major scaffold that is considered to be associated with many biological activities like antifungal, antibacterial, 22,23 anti-inammatory, 24,25 antiviral (anti-HIV), [26] [27] [28] anticancer, 29 anticonvulsant, 30 antidiabetic, 31,32 antitubercular activity, 33 etc. The presence of nitrogen atom and methylene groups in the thiazolidinone ring could be exploited for the development of new derivatives. The placement of arylidene groups at various positions on the thiazolidinone ring has been shown to improve the biological activity of the thiazolidinone template. 23 The antiviral activity of substituted-thiazolidin-4ones against HIV-1 RT was signicant, and emerged as potent non-nucleoside HIV reverse transcriptase inhibitors (NNRTIs). 34 They also emerged as strong hepatitis C virus RNA polymerase inhibitors exhibiting over 95% inhibition against HCV NS5B in vitro. 35 Table 2 lists some promising 4-thiazolidinone derivatives that target a variety of enzymes associated with various diseases.

These two scaffolds showed antiviral efficacy (both in vitro and in vivo) against a variety of viruses. These hybrids also share a similarity with some of the active compounds screened in vitro against SARS-CoV-2 (ESI, SD-1 †). This encouraged us to screen the benzoxazole-thiazolidinone hybrids (Fig. 2 ) as SARS-CoV-2 protease inhibitors. In order to achieve this, a combination of docking, dynamics, and free energy studies were employed in this study. "Tipranavir", a well-known antiviral protease inhibitor and a known inhibitor of 3CLp of SARS-CoV-2, was used as the reference drug. 12, [36] [37] [38] [39] Computational methods

X-ray crystal structure of 3CLp of SARS-CoV-2 (PDB ID 6LU7) was obtained from the protein data bank database. In case of PLp, the homology model of the protein (SMTL ID: 6w9c.1) was obtained from SWISS model database (ESI, SD-2 †). During the preparation of proteins, water molecules and other heteroatoms were excluded and the hydrogen atoms were added, Kollman charges were assigned, and Gasteiger charges were computed. The 2D structures of 81 B-T hybrids were sketched using Chemdraw soware ( Fig. 3 and Table 3 ). The 3D structures of the molecules were generated in their lower energy conformations using Gaussian 09 soware (semi-empirical method).

Later, the 3D structures of minimized ligands were loaded into Autodock Tools4 soware 48 and Gasteiger charges were computed for all the ligands.

AutoDock tools' molecular docking soware 48 was used for docking simulations of B-T hybrids with two SARS-CoV-2 protein targets. The grid box was created around the active site amino acid residues [(Cys145 & His41, in the case of 3CLp), (Cys111, His272 & Asp286, in the case of PLp)] with the size of 50 A Â 50 A Â 50 A. The Lamarckian genetic algorithm (LGA) was utilized to set up the docking with a maximum of 25 00 000 energy evaluations. The number of genetic algorithms was set to 250 runs and the population size was set to 300 for all the ligands. The program was executed retaining other parameters as default. Out of 250 poses observed in docking studies, the pose with best interactions with the active site residues and least binding energy was chosen as the best pose. The complexes of the protein with the best ligands were then subjected to MD simulations. The interactions of ligands and 

To study the stability aspects of complexes, MD simulations were performed using all-atom force eld (AMBER03) 50 in GROMACS2018 soware. 51 The protein was extracted from the complex and processed through H++ server 52 assuming a suitable pH (7.5 for 3CLp 53 and 8.0 for PLp 54 ). The ACPYPE server was employed to generate ligand topologies. 55 The complex topologies were created by appending the ligand topologies with their respective protein's topology. The complexes and proteins were then solvated using the TIP3P water model in a cubic box with a distance of 1.5 nm from the box boundaries and neutralized with 0.15 M NaCl. To hold all of the atoms within the simulation box, periodic boundary conditions were used. The systems were subjected to 50 000 steps of energy minimization using the steepest descent algorithm. The equilibration was performed in two steps employing NVT and NPT ensembles for 500 ps. Particle mesh Ewald was used to measure long-range interactions while a distance of 1.2 nm cut-off was used for short-range electrostatics and van der Waals. The LINCS (Linear Constraints Solver) algorithm was used to restrain all bonds. The production run was carried out for 50 ns using the NPT ensemble (at 1.0 bar and 300 K) and the trajectories were saved every 2 ps. Inbuilt modules of GROMACS were used to analyse root mean square deviation (RMSD), root mean square uctuations (RMSF), hydrogen bonds interactions, solvent accessible surface area (SASA) analysis, principal component analysis (PCA), and radius of gyration (R g ). 

Relative binding energies of the complexes were computed using MM-PBSA analysis on the 300 frames from the last 30 ns of the trajectory using g_mmpbsa module. 56 The binding free energies of complexes were calculated using the equations below.

where DG bind is the net binding free energy, DE MM is the molecular mechanics potential energy, DE vdw is the van der Waals energy, DE elec is the electrostatic energy, and ÀTDS accounts for the entropic energy. DE polar and DE SASA cannotes the polar contributions and the non-polar solvation energy respectively. The interaction entropy method was adopted from ref. [57] [58] [59] to calculate the contribution of entropy towards the binding energy of the complexes.

Molecular docking, in combination with molecular dynamics, is a well-established and effective method for drug design process. 60 It is a frequently used methodology for determining the stability of the molecular interactions between proteins and ligands. In the current study, the molecular interactions between B-T hybrids and selected SARS-CoV-2 proteases were analysed using molecular docking and the stability aspects of the complexes obtained from docking studies were studied through MD simulations.

A docking study was carried out to estimate the binding of B-T hybrids with the two proteases and gain insight into the interactions these hybrids make with proteases of SARS-CoV-2. In the case of 3CLp the active site residues Cys145 and His41, as well as neighbouring residues in the substrate-binding region such as Thr24, Thr25, Thr26, Thr45, Met49, Ser144, and Gly143, exhibited strong interactions by mediating hydrogen bonds with the B-T hybrids. These results indicated that the benzoxazole part of the hybrids consistently confront the catalytic dyad where it interacts especially with Cys145 in the active site, via hydrogen bonds with the thiol or the amide groups. In addition to this, it was observed that the aryl-thiazolidinone moiety of the molecule interacts with a group of threonine residues on the opposite side that belongs to S1 0 and S2 subunits of the binding pocket as shown in Fig. 4A . The binding mode of B-T hybrids on 3CLp demonstrated that the amide group present in the thiazolidinone unit of the molecule mediates two hydrogen bonds with Thr26 of the S2 subunit. The hybrids frequently interacted with Ser144 and Gly143 with a hydrogen bond involving oxygen and nitrogen atoms of the benzoxazole unit. Table 4 lists the top 15 B-T hybrids with binding energies less than À8.30 kcal mol À1 with their binding energy interactions with 3CLp. The docking results of 3CLp revealed that the compounds BT10 and BT14 were appropriately positioned within the binding pocket, which was surrounded by polar amino acid residues (Thr24, Thr25, Thr26, Thr45, Ser46, Asn142, Ser144), and hydrophobic residues (Leu27, Cys44, Met49, Phe140, Leu141, Met165) with binding energies À8.87 and À8.49 kcal mol À1 respectively. The binding energy of Tipranavir is À8.13 kcal mol À1 . The key interactions of Tipranavir, BT10 and BT14 at the binding site of 3CLp are shown in Fig. 4 .

The binding mode of B-T hybrids on PLp demonstrated various kinds of interactions that include hydrogen bonding, pi-pi stacking, pi-cation stacking, and hydrophobic interactions in the active site of PLp. The pi-pi stacking interactions were observed between benzoxazole moiety of a B-T hybrid and indole side chain of Trp106 while the pi-cation interactions 

were observed between aromatic moieties of the hybrids and histidine unit of His272. Surprisingly, the sulphur atom of thiazolidinone was found to act as an acceptor for a hydrogen bond interaction with the indole-ring nitrogen from Trp106 in MBT17 (ESI, Fig. 4 †) . The docking results of the hybrids PLp revealed that BT27 and MBT9 interacted and formed hydrogen bonds with the residue His272 of the catalytic triad, while Tipranavir, the reference drug formed H-bonds with Asp286 apart from other interactions with the binding pocket residues. MBT9 exhibited two additional hydrogen bonds with Trp106's indole-ring nitrogen, which is thought to play a critical role in the stability of tetrahedral intermediates (oxyanions) formed during the catalytic activity. 61 Table 5 lists the top 10 B-T hybrids with binding energies less than À7.50 kcal mol À1 with Fig. 5 . It was noticed that CBT and MBT series at times, had similar kinds of interactions with the proteases. For example, the pairs (CBT23, MBT23), (CBT26, MBT26) against 3CLp and the pairs (CBT13, MBT13), (CBT26, MBT26) against PLp had similar interactions with comparable binding energies, indicating that the groups at the 5th position on the benzoxazole ring (chloro and methyl groups) act similarly regardless of their chemical behaviour. These results indicate that B-T hybrids could be potential inhibitors of SARS-CoV-2 proteases and in order to further ascertain if these complexes are stable, the top B-T hybrids which complexed with 3CLp and PLp were subjected to MD.

The promising leads from B-T hybrids -BT10 and BT14, in complex with 3CLp, BT27 and MBT9, in complex with PLp, were chosen for the 50 ns MD simulations. In addition, Tipranavir in complex with both the proteases and apo proteins (proteins without any ligands in its binding pocket) were also subjected to simulations.

In general, RMSD determines the magnitude of a molecule's divergences from the initial reference structure (protein or ligand-protein complex). Large RMSD values are attributed to the changes in the conformation of the complex/protein during the simulation. From the protein backbone RMSD of the complexes (Fig. 6A) , it can be observed that the 3CLp trajectories attained stability within the rst 6 ns. However, the complexes exhibited relatively stable RMSD within the range of 1.5 AE 0.2 to 1.7 AE 0.3 A which was better than the average RMSD exhibited by the apo-protein (2.2 AE 0.4 A), suggesting that the interactions made by the ligands with the active site of 3CLp could improve the stability of the protease. The 3CLp-BT14 complex exhibited a slight deviation in the RMSD of the backbone from 25 to 30 ns and exhibited an average value of 1.6 A. As shown in Fig. 6B , Tipranavir and BT10 have a ligand backbone RMSD of 1 A in comparison to the protein backbone. The backbone of BT14 showed a deviation of around 4 A with respect to the protein backbone. Trajectory visualization of the 3CLp-BT14 indicated that the ligand BT14 rotated from the S2 subunit to the S3 subunit keeping the S1 subunit (catalytic dyad) as the point of rotation. The root-mean-square uctuation (RMSF) prole was calculated for each protein-ligand complex relative to the apo-state to gain additional insight into the stability of the complex. RMSF provides information on the protein residues' exibility. Fig. 6C depicts the overall residual uctuations of the 3CLp apostructure and complexes. All the 3CLp complexes exhibited a roughly similar pattern with the apo structure. In all the complexes, the loop region uctuated slightly, while active site residue regions remained relatively steady. The 3CLp-BT14 complex displayed a marginally higher uctuation rate, particularly for residues 100-120 and 125-137, signifying the greater conformational exibility than the other complexes.

The radius of gyration (R g ) is a measure of system density that signies the folding or unfolding of the protein/complex during the course of the simulation. The complexes and apo protein exhibited similar R g , and the average R g 's were in the range of 22.3 AE 0.1 A and 22.5 AE 0.1 A (Fig. 6D) .

Hydrogen bonds play an important role in ensuring the stability of a complex. The hydrogen-bonding interactions between binding site residues and B-T hybrids were analysed throughout the course of 50 ns MD simulation by using the g_hbond module of GROMACS. The average number of H bonds identied for Tipranavir (ESI, Fig. 1A †) is 1.8 with a maximum of 4H bonds, for BT10 (ESI, Fig. 1B †) it is 3.6 with a maximum of 7H bonds, and for BT14 (ESI, Fig. 1C †) it is 0.4 with a maximum of 4H bonds throughout the simulation. These results indicate that BT10 is capable of making strong interactions with key residues of the catalytic pocket which help stabilize the 3CLp-BT10 complex.

Results of RMSD analysis of the PLp complexes and protein revealed that all the trajectories achieved stability within 3 ns (Fig. 7A) . Thereaer, relatively negligible deviations were observed in the trajectories of the complexes indicating that the complexes had not undergone any signicant conformational changes. The RMSDs of the heavy atoms of ligands -Tipranavir and MBT9 with the backbone of PLp (Fig. 7B )suggested that they were unstable in the binding pocket while BT27 was relatively stable with an average RMSD of 3 A. The RMSF of the residues in the complexes followed a similar trend. As illustrated in Fig. 7C , the residues associated with the PLp- Tipranavir and the PLp-MBT9 complexes uctuated slightly more than the residues of apo structure and BT27 complex. Trajectory analysis of complexes revealed that both Tipranavir and MBT9, dried away from the catalytic triad, while BT27 retained the initial interactions in the active site.

PLp has two primary domainsthe catalytic or active site Cterminal domain and the ubiquitin-like (UBL) N-terminal domain. In terms of enzymatic functions and protein inhibition, the catalytic domain plays a crucial role. As shown in Fig. 9A , the catalytic domain can be further divided into thumb, palm, and ngers subdomains. 62 The active site which consists of three major residues (the catalytic triad), is situated between palm and thumb domains. Tipranavir began moving away from the catalytic site at 19 ns, stabilized at the thumb region of the target at 40 ns, and remained at the thumb region till the end of the simulation by interacting with the residues Asp40, Gly81, Arg82, Met84, Ser85, and Thr158. MBT9 shied from the catalytic site to the blocking loop 2 (BL2) of the palm region at 25 ns and remained there till the end of the 50 ns simulation by making interactions with the key residues of the BL2 in the palm subdomain, just outside the active site of PLp. Fig. 9B illustrates the interactions of MBT9 with BL2 loop of PLp.

The radius of gyration scores of PLp-apo and all complexes were calculated to be between 23.6 AE 0.1 A and 23.7 AE 0.1 A (Fig. 7D) , indicating that the ligand movement had no effect on the compactness of the target protein. The average number of H bonds identied for Tipranavir (ESI, Fig. 2A †) is 0.7 with a maximum of 4H bonds. In the case of BT27 (ESI, Fig. 2B †) it is 1.7 with a maximum of 5H bonds, while BT14 (ESI, Fig. 2C †) had 1.3 average hydrogen bonds with a maximum of 3H bonds during the simulation.

The relative binding free energy was estimated using MM-PBSA to deduce the forces responsible for the affinity of a ligand with its receptor. 300 equidistant frames were extracted from the trajectory of MD production run between 20 and 50 ns for the MM-PBSA analysis, while all the frames between 20 to 50 ns trajectory were used for computing the entropy contribution (by interaction entropy method) towards the binding energy. The resulting DG Bind of 3CLp-Tipranavir, 3CLp-BT10, and 3CLp-BT14 complexes was found to be À8.65 kcal mol À1 , À4.87 kcal mol À1 , and À6.83 kcal mol À1 , respectively (Fig. 8C ) whereas DG Bind of PLp-Tipranavir, PLp-BT27, and PLp-MBT9 was found to be À4.74 kcal mol À1 , À7.76 kcal mol À1 , and -4.68 kcal mol À1 , respectively (Fig. 8D) . Further, binding free energy was decomposed into the individual energy terms to understand the impact of each energy component on the binding process ( Fig. 8A and B) . The van der Waals energy component (DE vdw ) was the primary contributor to the binding free energy in each complex. Despite its substantial electrostatic component (DE elec ), the 3CLp-BT10 complex had a higher value of binding free energy due to the large unfavorable polar component of solvation (DE polar ), whereas the 3CLp-BT14 complex had a relatively lower binding free energy due to the lower polar component of solvation (DE polar ). The binding free energy of BT27 to PLp was compared to that of other ligands and it was observed that BT27 has a better affinity to PLp than others. The PLp-Tipranavir complex suffered from low contributions from van der Waals (DE vdw ) and electrostatic (DE elec ) energy terms.

To gain more insights into the individual energy contributions of the active site residues which are responsible for the stability of the ligand, per residue energy decomposition study was conducted. For the 3CLp-Tipranavir complex, the major contributor for the favourable energy contributions was His41 (À8.4 kcal mol À1 ) and Thr25 (À8.4 kcal mol À1 ). For the 3CLp-BT10 complex, the major contributions originated from Cys145 (À4.3 kcal mol À1 ) and Thr25 (À4.6 kcal mol À1 ). This shows that both Tipranavir and BT10 actively interacted with the catalytic dyad. As BT14 relocated from its original position, in the case of 3CLp-BT14 complex, relocated the favourable energy contributions are therefore derived from a different set of residues located in the S2 and S3 subunits of the binding pocket. The main contributors for the binding energy of the 3CLp-BT14 complex are Gln189 (À4.8 kcal mol À1 ), and Met165 (À5.1 kcal mol À1 ) whereas Arg188 contributed to the binding energy unfavourably by 3.0 kcal mol À1 . The PLp-Tipranavir complex has a much weaker binding free energy (À4.68 kcal mol À1 ), in comparison to the PLp-BT27 complex (À7.76 kcal mol À1 ). For the PLp-BT27 complex, Asn109 (À6.6 kcal mol À1 ) and Leu289 (À5.2 kcal mol À1 ) are the main contributors to the binding energy. The crucial BL2 active site residues -Pro248, Tyr264 and Thr265, made signicant contributions to the free energy of the PLp-MBT9 complex with each contributing À5.8 kcal mol À1 , À5.0 kcal mol À1 , and À5.7 kcal mol À1 , respectively. Overall, it was observed that the catalytic triad residues made repulsive interactions with Tipranavir. The active site of PLp is in a solvent-exposed region located on the protease's surface, causing the ligands with weak interactions to exit the catalytic pocket. This could explain the migration of Tipranavir and MBT9 from the active site. It is worth noting that, in spite of these conditions, BT27 was stable in the active site as it formed strong interactions with the binding site residues (other than the catalytic triad), preventing it from driing away from the vicinity of the catalytic triad.

The top two hits in case of 3CLp were synthesized in a simple four step process. The synthetic approach for preparing B-T hybrids is illustrated in Scheme 1. The synthesis began with the preparation of isoperthiocyanic acid. Kalson method was used to prepare isoperthiocyanic acid from ammonium thiocyanate. 63 The synthesis of 2-thioureidobenzoxazole was carried out in a single step from isoperthiocyanic acid and 2-aminophenol in good yield as previously reported. 64 Reaction of 2thioureidobenzoxazole with chloroacetic acid as cyclizing reagent, in presence of anhydrous sodium acetate in boiling ethanol, afforded 4-thiazolidinone linked with benzoxazole moiety with an imine linkage. It was clearly seen from the docking studies that the arylated hybrids exhibited strong Scheme 1 General synthetic strategy for the synthesis of B-T hybrids. interactions than BT1 with both the proteases of SARS-CoV-2. Hence, in the nal step, the hybrids were arylated with corresponding aldehydes using the methylene carbon of thiazolidinone moiety.

All reagents and solvents were purchased from Sigma-Aldrich and SD Fine Chemicals. The solvents were distilled prior to use. The conversions were monitored by thin layer chromatography (TLC) using Merck TLC silica gel 60 F 254 . Compounds on TLC plate were visualized by ultraviolet light at 254 nm and 365 nm. UV-Vis absorption spectra were obtained from Hitachi UH5300 spectrophotometer. IR spectra were recorded on a Per-kinElmer FT-IR C114728 spectrophotometer with ATR accessory. 1 H NMR and 13 C NMR spectra were recorded at 400.23 and 100.64 MHz on a Bruker 400 MHz spectrometer, respectively. Using tetramethyl silane (TMS) as the internal standard, and dimethyl sulfoxide (DMSO-d 6 ) as the solvent. HRMS were recorded on Agilent Technologies 6550 iFunnel LC-MS Q-TOF with ESI (electrospray ionization).

To a magnetically stirred mixture of 2-thioureidobenzoxazole (1.158 g, 0.006 mol) and chloroacetic acid (0.728 g, 0.007 mol) in 15 ml of ethanol, anhydrous NaOAc (0.631 g, 0.007 mol) was added, and heated at reux for 24 hours. Aer the reaction, the precipitate (1.291 g, 90%) was collected and recrystallized from ethyl acetate as slender off-white needles; mp 255-257 C; l max (THF): 325 nm; n max ( 

A mixture of BT1 (0.466 g, 0.002 mol), corresponding aldehyde (0.003 mol) and the base catalyst, anhydrous NaOAc (0.328 g, 0.004 mol) in glacial acetic acid (15 ml) was reuxed overnight. Aer the completion of the reaction, the precipitate was collected through ltration, dried and recrystallized from THF.

BT10 Conclusion COVID-19 is the largest global pandemic currently and it has emerged as the deadliest infection of the century. It's critical to nd novel compounds that can inhibit viral replication, in turn reducing the pace of infection. The current study used molecular docking and molecular dynamics to identify the potential of B-T hybrids as inhibitors of the SARS-CoV-2 proteases. Following the completion of molecular docking and molecular dynamics simulations, the ndings were analysed and the leads were identied. BT10 exhibited good interactions with the 3CLp by mediating the highest number of hydrogen bonds (7), and maintained the interactions throughout the MD simulation with binding free energy of À4.87 kcal mol À1 . BT14, with the best binding energy (À6.83 kcal mol À1 ), interacted strongly with the loop region residues. In the case of PLp complexesalthough the residues of the catalytic triad of PLp made repulsive interactions, BT27 could tightly bind to the binding site and sustained the interactions throughout the simulation with a good binding energy (À7.76 kcal mol À1 ). It was observed that the lead inhibitors of PLp bind to the BL2 loop rather than the catalytic site, preventing the substrate entry to the catalytic site. 65 MBT9 showed a tendency to bind to the blocking loop2 (BL2) with a binding free energy -4.68 kcal mol À1 . Hence, MBT9 could act as an allosteric inhibitor of PLp. Finally, a four-step synthetic procedure was employed to synthesize the B-T hybrids starting from ammonium thiocyanate. The promising compounds in the case of 3CLp were synthesized and characterized using IR, NMR and mass spectroscopic techniques. To further evaluate the therapeutic potential of the proposed ligands, the in vivo and in vitro studies could be carried out. The potential of the B-T hybrids as covalent binders to 3CLp can also be explored, as the conjugate amide group present in the thiazolidinone scaffold of the molecule can serve as a Michael acceptor (ESI, Fig. 6 †) .

Vijay Sai Krishna Cheerala: conceptualization, methodology, data curation, visualization, writingoriginal dra preparation. Prasanth Ghanta: methodology, validation, writingreviewing and editing. Sundaresan Chittor Neelakantan: conceptualization, validation, supervision, writingreviewing and editing.

The authors declare that there are no conicts of interest in this study.

Design of Hybrid Molecules for Drug Development

Heterocyclic Benzoxazole Compositions as Inhibitors of Hepatitis C Virus

3H-1,2,4-dithiazole-3-thione as a Synthon: new synthesis of 2-thioureidobenzheteroazoles