Synthesis, antimicrobial, molecular docking and molecular dynamics studies of lauroyl thymidine analogs against SARS-CoV-2: POM study and identification of the pharmacophore sites

Graphical abstract

Nucleoside precursors and nucleoside analogs occupy an important place in the treatment of viral respiratory pathologies, especially during the current COVID-19 pandemic. From this perspective, the present study has been designed to explore and evaluate the synthesis and spectral characterisation of 5-O-(lauroyl) thymidine analogs 2-6 with different aliphatic and aromatic groups through comprehensive in vitro antimicrobial screening, cytotoxicity assessment, physicochemical aspects, molecular docking and molecular dynamics analysis, along with pharmacokinetic prediction. A unimolar one-step lauroylation of thymidine under controlled conditions furnished the 5-O-(lauroyl) thymidine and indicated the selectivity at C-5́ position and the development of thymidine based potential antimicrobial analogs, which were further converted into four newer 3-O-(acyl)-5-O-(lauroyl) thymidine analogs in reasonably good yields. The chemical structures of the newly synthesised analogs were ascertained by analysing their physicochemical, elemental, and spectroscopic data. In vitro antimicrobial tests against five bacteria and two fungi, along with the prediction of activity spectra for substances (PASS), indicated promising antibacterial functionality for these thymidine analogs compared to antifungal activity. In support of this observation, molecular docking experiments have been performed against the main protease of SARS-CoV-2, and significant binding affinities and non-bonding interactions were observed against the main protease (6LU7, 6Y84 and 7BQY), considering hydroxychloroquine (HCQ) as standard. Moreover, the 100 ns molecular dynamics simulation process was performed to monitor the behaviour of the complex structure formed by the main protease under in silico physiological conditions to examine its stability over time, and this revealed a stable conformation and binding pattern in a stimulating environment of thymidine analogs. Cytotoxicity determination confirmed that compounds were found less toxic. Pharmacokinetic predictions were investigated to evaluate their absorption, distribution, metabolism and toxic properties, and the combination of pharmacokinetic and drug-likeness predictions has shown promising results in silico. The POM analysis shows the presence of an antiviral (O1 δ-, O2 δ-) pharmacophore site. Overall, the current study should be of great help in the development of thymidine-based, novel, multiple drug-resistant antimicrobial and COVID-19 drugs.

Introduction
Nucleoside analogs have, in recent times, been used as the backbone for many drugs in the treatment of infectious diseases caused by HIV, hepatitis B or C viruses, and the herpes viruses [1]. Nucleoside analogs are the subunits of DNA and RNA, and structurally they consist of a sugar moiety linked to a nitrogen base through an N-β-glycosidic bond. Interestingly, biological activity, as well as the chemical and physical properties of nucleoside analog, can be tuned remarkably by the simple alteration of the sugar moiety of a nucleoside with different substituents and heteroatoms, among others [2]. The nucleoside analogs usually show antiviral activity by inhibiting viral replication through the blockage of cellular division; impairment of DNA/RNA synthesis or inhibition of cellular or viral enzyme activity. The nucleoside analogs, Telbivudine (hepatitis B inhibitor), Entecavir (HIV/ AIDS and Hepatitis B inhibitor), Clevudine (Hepatitis B inhibitor), Zalcitabine (reversetranscriptase inhibitor), Taribavirin (prodrug of Ribavirin), Stavudine (HIV/AIDS inhibitors), Lamivudine (first-generation nucleoside reverse transcriptase inhibitor), Cordycepin (RNA synthesis inhibitor), and Cordycepin Triphosphate (polyadenylation inhibitors, antineoplastic, antioxidant, and anti-inflammatory agent), have been used in the treatment of many viral diseases. Moreover, nucleoside analogs have immense clinical importance as medicinal agents because of their antiviral and anticancer activities [3]; they have been the "drug of choice" for the treatment of various viral diseases, such as ebola [11], dengue, and Zika [4]. Additionally, 2ʹ-deoxynucleosides such as doxuridne, trifluridine, doxudine, vidarabine, and brivudine, have been used to treat herpes virus infections [5,6]. The modifications in the sugar moiety, namely ribofuranose or deoxyribofuranose of nucleosides, include changes in the sugar substituents, replacement of the oxygen with another atom, and the addition of heteroatoms in the sugar ring, ring size variations, and replacement with acyclic moiety [7][8][9][10][11][12]. These alterations can produce remarkable variations in biological activity and the degree of selective toxicity, according to their respective chemical and physical properties [13][14][15][16][17][18][19][20]. The modified compounds have exhibited broad-spectrum biological activity. Thymidine (1) derivatives such as telbivudine are antiviral drugs used in HBV treatment [21]. Another nucleoside, cytidine, is found in offal meat and pyrimidine-rich foods such as beer, tomatoes, broccoli and oats. Cytidine is one of the RNA components that translate instructions from DNA to proteins [22]. Some recent research revealed that thymidine analogs, as well as nucleobases, are a pharmacologically diverse family, which includes cytotoxic compounds, antimicrobial agents and immune-suppressive molecules [23][24][25][26]. However, one study indicates that selective acylation of thymidine (1) causes complexities due to the appearance of one primary OH group, one secondary OH group and an -NH group. Hence, various methods like direct acylation [8], the dibutyltin oxide method [27,28], coupling techniques, and the enzymatic method have been reported. Among these methods, the direct method generally furnished 5-O-substituted analogs [27,28], whereas the dibutyltin oxide method gave 3-O-substituted analogs [8,27,28]. However, most of the time, the direct method produced higher yields at a lower cost. The novel COVID-19 virus spread around the world from Wuhan, China in December 2019. The World Health Organization (WHO) declared the disease a public pandemic emergency on March 11th 2020, as the cause of Severe Acquired Respiratory Syndrome (SARS-CoV-2). As a result, the development of multiple drug-resistant solutions for pathogenic organisms and coronavirus disease 2019 (COVID-19) has become a current global concern for researchers. Most of the proposed COVID-19 drugs, such as remdesivir and favipiravir, are based on a nucleoside structure, and modifications of the hydroxyl (-OH) group of nucleosides, including thymidine, have revealed some potent SARS-CoV-2 candidates [29][30][31] and antimicrobial agents in the current study. However, there is an urgent need for the development of new chemotherapeutic agents with a novel mode of action to combat these multiple drug-resistant organisms, including SARS CoV-2. One of our recent works based on PASS prediction, In vitro antimicrobial, physicochemical and molecular docking studies revealed that thymidine derivatives are potential inhibitors of SARS CoV (1Q2W) [32].
Keeping in mind these features mentioned above, as well as the subsequent aim of finding novel drugs [32], we have made some advancements in our present study. Different acylating agents are inserted in the 5́ and 3́ thymidine structures to synthesize novel thymidine analogs. The in vitro antimicrobial evaluation of several thymidine-based analogs 2-6 with different aliphatic and aromatic chains against eight pathogens, molecular docking against the coronavirus (SARS CoV-2) main protease (6LU7, 6Y84 and 7BQY), along with the prediction of activity spectra for substances (PASS), are reported on herein. To confirm the stability of the docked complexes, molecular dynamics was performed for 100 ns. Cytotoxicity was determined to explore the safety profile synthesize thymidine analogs. Then we analyze the series of analogs 2-6 to identify their pharmacophore sites. Furthermore, steps have been taken to optimize the synthesized thymidine analogs to investigate their physicochemical behavior based on the DFT approach with pharmacokinetic and drug-likeness properties.

Material and equipment
The solvents used were of analytical grade and were purified using standard procedures. All reagents used are commercially available from Sigma-Aldrich (Germany) and were used as received unless otherwise specified. Melting points (mp) were determined using an electrothermal melting point apparatus. Evaporations were performed under diminished pressure on a Büchi rotary evaporator. Infrared spectral analyses were recorded using a Fourier-transform infrared (FTIR) spectrophotometer (IR Prestige-21, Shimadzu, Japan) at the Department of Chemistry, University of Chittagong. A Brucker advance DPX 400 MHz using tetramethylsilane as the internal standard was employed for recording the proton nuclear magnetic resonance ( 1 H NMR) spectra at WMSRC, JU, Bangladesh. Mass spectra of the synthesised compounds were obtained by liquid chromatography-electrospray ionisation tandem mass spectrometry in positive ionization mode. Thin-layer chromatography (TLC) was performed on Kieselgel GF 254 (Germany), and the chromatogram was visualised by spraying the plates with 1% H 2 SO 4 , followed by heating the plates at 150-200 • C until coloration appeared. Column chromatography was performed using silica gel G 60 . To identify the drug interactions with receptor proteins, molecular docking is the most suitable tool. The blind docking method employs a search throughout the whole surface of the protein molecule for binding sites. The following software was used in the present study: and SwissADME free web tools (https://www.swissadme.ch) were employed to calculate the pharmacokinetic properties, and (ix) YASARA for molecular dynamics.

5-O-(Lauroyl)thymidine (2)
A solution of thymidine (1) (200 mg, 0.82 mmol) in dry pyridine (3 ml) was cooled to − 5 • C, whereupon lauroyl chloride (0.205 ml, 1.1 M eq.) was added to it. The mixture was stirred at this temperature for 6 h and then stirred overnight at room temperature. T.L.C. examination (methanol-chloroform, 1:4) indicated full conversion of the starting material into a single product (R f = 0.55). The solution was poured into ice water with constant stirring. It was then extracted with chloroform (3 × 10 ml). The combined chloroform layer was washed successively with dilute hydrochloric acid (10%), saturated aqueous sodium hydrogen carbonate (NaHCO 3 ) solution and distilled water. The organic layer was dried with anhydrous magnesium sulphate (MgSO 4 ) and filtered, and the filtrate was concentrated under reduced pressure to leave a syrup. Purification using silica gel column chromatography (with methanol-chloroform 1:4 as eluant) afforded the title compound (2) (130 mg, 62.86%) as a white crystalline solid which was used in the next stage. Recrystallisation from chloroform-n-hexane provided the lauroyl derivatives (2)  General procedure for the preparation of lauroyl derivatives (3-6).

3-O-(4-t-Butylbenzoyl)-5-O-(lauroyl)thymidine (3)
A suspension of lauroyl derivative (2) (130 mg, 0.30 mmol) in dry pyridine (3 ml) was cooled to − 5 • C whereupon 4-t-butyl benzoyl chloride (129 mg, 2.2 M eq.) was added to it. The reaction mixture was stirred for 6 h at this temperature and then stirred overnight at room temperature. The progress of the reaction was monitored using T.L.C. (methanol-chloroform, 1:8), which showed the complete conversion of reactant into a single product (R f = 0.51). The excess reagent was destroyed by the addition of a few pieces of ice, and the reaction mixture was processed as described earlier. The resulting syrup was purified using column chromatography to afford compound (3) (3 ml) was cooled to − 5 • C whereupon myristoyl chloride (0.129 ml, 2.2 M eq.) was added to it. The mixture was stirred at this temperature for 6 h and then stirred overnight at room temperature. The progress of the reaction was monitored by T.L.C. (methanol-chloroform, 1:8) which exhibited complete conversion of the starting material into a single product (R f = 0.52). The excess reagent was destroyed by the addition of a few pieces of ice and the reaction mixture was processed as

In vitro antimicrobial activity test 2.3.1. Test organisms
The tested microorganisms (bacteria and fungi) were collected from the Microbiology Laboratory, Department of Microbiology, University of Chittagong, Chattogram, Bangladesh. As shown in Table S1, five human pathogenic bacteria and two plant pathogenic fungi were used for the antimicrobial evaluation.

Disc diffusion test to check antibacterial susceptibility
Antibacterial susceptibility testing of the compounds was performed using standard disk diffusion assay according to CLSI protocol [33,34]. The inoculum was prepared by picking the test bacteria from plates with a sterile inoculating loop and suspending it in sterile normal saline (0.85%). The density of the suspension to be inoculated was determined through a comparison with the McFarland 0.5 standard. After that, the test bacteria were swabbed uniformly over the Mueller-Hinton agar (MHA) plates, and filter paper discs (6 mm, Hi-media) incorporated into the test compounds were placed in the agar medium. After that, the plates were incubated for 18 h at 37 • C and the diameters of zones of inhibition around the discs were measured in millimeters (mm). The experiments were repeated three times.

Determination of MIC and MBC using the micro-broth dilution method
Microdilution methods are the most appropriate ones for the determination of MIC and MBC values since they offer the possibility of estimating the concentration of the antimicrobial agent being tested in the broth medium, and they may be used to quantitatively measure the in vitro antimicrobial activity against bacteria. A two-fold dilution of the test compound, 0.125, 0.25, 0.5, 1, 2, 4, 8, 16, 32, 64 and 128 µg/ mL) was prepared by dispensing approximately 100 μL of each 2x compound solution in column 1 of a 96-well microtiter plate containing 100 μL of MHB broth, and from columns 2 to 10, compound solutions were serially diluted. Each well, except the sterility control well (blank), was inoculated with a 5 μL test strain inoculum of 10 4 to 10 5 CFU/mL prepared after dilution of a standardised suspension adjusted to a 0.5 McFarland scale. The microplate was loosely wrapped with clingfilm wrap to prevent bacterial dehydration. Each plate included a growth control in column 11 (without compound) and sterility control (only media) in column 12. Then, the inoculated 96-well microtitre plates were incubated at 37 • C for 24 h. To each well, 10 µL of 2, 3 and 5-triphenyl tetrazolium chloride 0.5% (w/v) solution was added. The cultures were then incubated at 37 • C for 24 hrs. The appearance of a red color indicated the growth of bacteria, and MIC and MBC were interpreted visually.

Screening of mycelial growth
The poisoned food technique [6,12] was employed to screen the antifungal activity in which potato dextrose agar (PDA) was used as the culture medium. The test compounds were dissolved in dimethyl sulfoxide (DMSO) at a concentration of 1% (w/v), and 0.1 ml of the solution (containing 1 mg of the respective test compound) was transferred to a sterile petri dish using a sterilised pipette. After that, 20 ml of the medium was poured into the Petri dish and allowed to solidify. The inoculation was then performed in the centre of each petri dish with a 5 mm mycelium block of each fungus. The blocks were placed at the centre of each petri dish in an inverted position to maximise the contact between the mycelium and the culture medium. The inoculation plates were incubated at 27 ± 2 • C, and the experiment was conducted in triplicate. A control sample (i.e., PDA without test chemicals) was also maintained under the same conditions. After 5 days of incubation, the diameter of the fungal radial mycelial growth was measured. An average of three measurements was considered for the radial mycelia growth of the fungus in mm.

Structure-activity relationship (SAR)
Structure-activity relationship (SAR) analysis was employed to predict antimicrobial activity from the molecular structure of the pharmaceutical target. This well-known technology is often used in the drug designing process to guide the acquisition or synthesis of new compounds with desirable properties. In the present study, the SAR research was analyzed according to the Hunt [34] and Kim [35] membrane permeation concept.

PASS and bioactivity prediction
The online web application PASS (https://www.pharmaexpert.ru/pa ssonline/) has been employed to calculate the antimicrobial activity spectrum of the selected thymidine analogs [36]. Firstly, the thymidine analog structures were drawn, and then changed into their smiles formats by using the renowned SwissADME free online application (https ://www.swissadme.ch) to determine the antimicrobial spectrum using the PASS web tool. This server can surmise over 4000 types of antimicrobial functions, together with drug and non-drug activity, which helps to suggest the best potential objects with 90% validity. PASS outcomes are revealed by Pa (probability for active molecule) and Pi (probability for inactive molecule). Having potentialities, the Pa and Pi scores vary in the range of 0.00 to 1.00, and usually, Pa + Pi ∕ = 1, as these potentialities are predicted freely. The biological actions with Pa > Pi are only thought of as probable for a selected drug molecule. PASS calculation outcomes were explained and used flexibly, viz. (i) when Pa is greater than 0.7, the probability to identify the activity is analytically high, (ii) if 0.5 < Pa < 0.7, the probability to identify the activity is analytically low, again, the molecule is perhaps not very similar to well conversant pharmaceutically used drugs and (iii) if Pa < 0.5, the potentiality to identify the activity analytically is less. In the present study, the Molinspiration online server (https://www.molinspiration.com/cgi-bin/properties) was utilised to analyse the drug-like properties of lead compounds. The Molinspiration cheminformatics engine allows the fast prediction of biological activity and virtual screening of large collections of molecules, and the selection of molecules with the highest probability, to show biological activity. The screening is based on the identification of fragments or substructure features typical of the active molecules.

Details of geometry optimization
In computational chemistry, quantum mechanical methods are widely used to calculate thermal, molecular orbital and molecular electrostatic properties [37]. Geometry optimisation and further modification of all synthesised analogs were carried out using the Gaussian 09 program [38]. Density functional theory (DFT) with Beck's (B) [39] three-parameter hybrid model, Lee, Yang and Parr's (LYP) [40] correlation functional under a 3-21G basis set, was employed to optimize and predict their thermal and molecular orbital properties. Dipole moment, enthalpy, free energy, and electrical energy were calculated for all the analogs. Frontier molecular orbital features, HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) were counted at the same level of theory. For each of the thymidine analogs, the HOMO-LUMO energy gap, hardness (η), and softness (S) were calculated from the energies of the frontier HOMO and LUMO as reported, considering Parr and Pearson's interpretation of DFT and Koopman's theorem [40] on the correlation of chemical potential (µ), electronegativity (χ) and electrophilicity (ω) with HOMO and LUMO energy (ε). The following equations were used to calculate global chemical reactivity by analyzing molecular orbital features.

Protein selection and molecular docking
The crystal 3D structure of SARS-CoV-2 main protease proteins (pdb: 6LU7, 6Y84 and 7BQY) were recuperated in the pdb format from the protein data bank [41]. PyMol (version 1.3) software packages were employed to remove all heteroatoms and water molecules [42]. Energy minimization of the protein was performed by using a Swiss-PdbViewer (version 4.1.0) [43]. Furthermore, a molecular docking study against the SARS-CoV-2 main protease proteins 6LU7, 6Y84 and 7BQY (Fig. S1) was conducted on the optimized drugs. Finally, the PyRx application (version 0.8) was used to carry out molecular docking simulation [44], envisaging the target protein as a macromolecule and the thymidine analogs as a ligand. The protein and ligands were input by converting the pdb format to pdbqt, and the AutoDock Tools of the MGL software package were used to perform this job. In AutoDockVina, the size of the grid box was maintained at (53.3467, 65.9465, and 57.6051 Å for 6LU7), (36.0681, 63.8808, and 61.9745 Å for 6Y84) and (49.8375, 64.9838, and 59.1841 Å for 7BQY) along the X, Y, and Z axes. After docking, both the structures of macromolecule and ligand were saved in pdbqt format, and Accelrys Discovery Studio (version 4.1) was employed to explore the results of docking and to predict the non-bonding interactions among cytidine derivatives and amino acid chains of the receptor protein [45]. PDBsum online server was also used to check the validation of the main protease receptor with a Lig-plot and Ramachandran plot (Fig. S2), which revealed that >95.0% of residues were in the allowed region and no residues were missed.

Molecular dynamics simulation
The molecular dynamics simulation study was conducted in YASARA dynamics [46] with the aid of the AMBER14 force field [47]. The docked complexes were initially cleaned and optimised, and hydrogen bond networking was oriented. A cubic simulation cell was created whereby the TIP3P solvation model was used under periodic boundary conditions [48]. The simulation cell was extended 20 Å for each case from the docked complexes. The physiological condition of the simulation cell was set at 310 K, pH 7.4 and 0.9% NaCl. The initial energy minimisation was conducted using the steepest gradient algorithms through simulated annealing methods (5000 cycles) [49]. The long-range electrostatic interactions were calculated by the Particle Mesh Ewald (PME) using a cut off radius of 8.0 Å [50]. The time step of the simulation cell was set at 2.0 fs. The simulation's trajectories were saved after every 100 ps. By following the constant pressure and Berendsen thermostat, the simulation was extended for 100 ns. The simulation's trajectories were used to calculate the root means square deviation (RMSD), root mean square fluctuation (RMSF), solvent accessible surface area (SASA), and the radius of gyration (Rg), and hydrogen bond [51][52][53][54][55].
Therefore, the YASARA trajectories were utilised in the calculation of the binding free energy from MM-PBSA methods, where positive energy indicates better energy [56]. The equation used for the calculations for the binding free energy calculations was:

Cytotoxic activity evaluation
To determine the toxicity of the thymidine analogs the brine shrimp lethality assay (BSLA) technique was performed that stated by [57]. The tested thymidine analogs were dissolved in DMSO and prepared at four different concentrations of 20, 40, 80 and 160 µL by adding NaCl solution to each vial up to 5 ml volume. The final concentrations of the sample in the vials were 4, 8, 16 and 32 µL in vials labelled as A, B, C and D, respectively. The experiment was performed three times for each concentration and 10 brine shrimps nauplii were added to each vial. A control test was performed in a vial containing 10 nauplii in 5 ml seawater. Then the vials were kept for incubation at room temperature for 24 hr and 48 hr. After incubation, the vials were observed using a magnifying glass and the number of survivors in each vial was enumerated and noted. From the data, the average percentage of mortality of nauplii was calculated for each concentration. No deaths were found for controls.

Pharmacokinetic and drug-likeness prediction
The prediction of ADMET properties in drug development is important to research that helps to prevent drug failure during the clinical stages. For this reason, the best-identified esters were evaluated using pkCSM [58] for their in silico pharmacokinetic parameters to avoid the failure of the esters during clinical trials and to improve their chances of reaching the stage of potential candidate drugs. Its absorption in the human intestine; percolation of the blood-brain barrier and the central nervous system, and its metabolism indicate the chemical biotransformation of a drug by the body, total clearance of drugs and the toxicity levels of the molecules. The drug-likeliness of a molecule was expressed using Lipinski's rule of five parameters (molecular weight < 500 Da, no more than 5 hydrogen bond donors, no more than 10 hydrogen bond acceptors and logP should not be greater than 5). Lipinski's rule of five properties was obtained from the SwissADME server (https://www. swissadme.ch/index.php).

Statistical analysis
For each parameter investigated, experimental results were presented as mean ± standard error for three replicates. Two-tailed Student t-tests were used as appropriate for statistical analysis. Only ρ values that were<0.05 were considered to be statistically significant.

Characterization
The main objective of this research work (Fig. 1) reported in this paper was to carry out selective lauroylation (Scheme 1) of thymidine (1) with lauroyl chloride using the direct acylation method. The resulting lauroylation product was transformed into several analogs employing various aliphatic and aromatic agents (Table 1). These thymidine analogs were used in an attempt to perform an antimicrobial assessment, along with the prediction of PASS and bioactivity spectra, to realize the mode of their antimicrobial behavior. As a result of screening synthesized thymidine analogs for potential antimicrobial activity, it was observed that most of the analogs have potent bactericidal and fungicidal in vitro activity against tested bacterial and fungal pathogens. Furthermore, synthesized analogs were analyzed for their physicochemical properties using DFT-geometry optimization. The observed results were then rationalized by molecular docking and molecular dynamics simulation, along with the prediction of in silico pharmacokinetic and drug-likeness properties.

2d-NMR
The COSY spectrum of compound (2) showed that the starting point could well be the signal from Ar-NH proton which is the most downfield and therefore readily assigned. Thus the signal from Ar-NH at the bottom left of the diagonal has a cross-peak labelled as Ar-NH, H-5b connecting it to the signal from H-5b. Thus, the Ar-NH proton around δ 8.71 is coupled to the hydrogen whose signal appears around δ 4.31 (i.e. H-5b proton). Similarly, the signal from H-5b is further connected by a crosspeak to the signal from 2H, CH 3 (CH 2 ) 9 CH 2 CO-to show the coupling between H-5b and 2H, CH 3 (CH 2 ) 9 CH 2 CO-. The lauroyl derivative (2) was then allowed to react with a unimolecular amount of 4-t-butyl benzoyl chloride in dry C 6 H 5 N at freezing temperature. A usual work-up procedure provided the 4-t-butyl benzoyl derivative (3) in good yields. In its 1 H NMR spectrum, two two-proton multiplets at δ 8.04 (Ar-H) and δ 7.49 (Ar-H) and three singlets at δ 1.37, δ 1.29 and δ 1.22 {9H, (CH 3 ) 3 C-} corresponded to the presence of one 4-t-butyl benzoyl group in compound (3). Further support for the structure accorded to compound (2) was obtained through the preparation of its myristoyl derivative (4), and palmitoyl derivative (5). Thus, the treatment of compound (2) with myristoyl chloride and palmitoyl

Drugs
Chemical Structure  (5) with good yields. It was possible to propose a structure for compounds 4 and 5 by analysing their complete spectroscopic and elemental data. Finally, 4-bromobenzoyl chloride was used to derivatise compound 2 using the direct acylation method. After the usual workup and purification procedure, the 4-bromobenzoyl derivative (6) was obtained. The formation of a monosubstitution product was revealed by its 1 H NMR spectrum, which showed one two-proton multiplet at

Antibacterial potentiality
The use of nucleoside analogs can be considered a novel option as they are expected to act at the genomic level, thereby interfering with the transcription or replication processes required for microbial survival. As there are no alternative pathways in the pathogens for these basic metabolic processes, the nucleoside analogs, by inhibiting these basic pathways, are effective antimicrobial agents [59]. The synthesized thymidine analogs (2-6) were tested for their in vitro antimicrobial activity against a series of strains of Gram-positive bacteria (Bacillus subtilis and Bacillus cereus), and Gram-negative bacteria (Escherichia coli, Salmonella typhi and Pseudomonas aeruginosa). The methods included the disk diffusion method and broth microdilution method for MIC and MBC determination, in accordance with [33]. The results of the preliminary antibacterial activities are shown in Table 3. Analog 3 exhibited the highest zone of inhibition against Bacillus cereus (15.25 mm) and Bacillus subtilis (13.25 mm), and analog 4 showed a notable zone of inhibition against Bacillus subtilis (12.75 mm) and Bacillus cereus (12.75 mm). Furthermore, analog 2 also exhibited a promising zone of inhibition against both Gram-positive bacteria Bacillus subtilis (12.75 mm) and Bacillus cereus (10.75 mm).
Interestingly, analog 5 displayed a fluctuating zone of inhibition,   The data are presented as mean ± SD and the values are represented for triplicate experiments. Statistically significant inhibition (p < 0.05) is marked with an asterisk (*) for test compounds and a double asterisk (**) for the reference antibiotic azithromycin. NI = No inhibition.
where a 14.50 mm zone was found for Bacillus subtilis and a 9.00 mm inhibition zone was observed against Bacillus cereus. Analog 6 exhibited a moderate zone of inhibition against only one Gram-positive bacteria Bacillus cereus (9.50 mm). Based on the above observation, the Grampositive antibacterial activity of the thymidine analogs can be ordered as 3˃4˃2˃5˃6˃ (Fig. S5). Furthermore, based on the data presented, analogs 2 and 5 exhibited a zone of inhibition against all the three Gram-negative bacteria: Escherichia coli, Salmonella typhi and Pseudomonas aeruginosa, and analog 3 showed inhibition against only two bacteria: Salmonella typhi and Pseudomonas aeruginosa. Moreover, analog 4, showed inhibition against Escherichia coli, and analog 6 showed inhibition against Salmonella typhi. These derivatives provided moderate to high activity against all three Gram-negative bacteria: Escherichia coli, Salmonella typhi and Pseudomonas aeruginosa. Based on this activity, the compounds can be ordered as 2˃5˃3˃4˃6. Based on the results for zones of inhibition, the compounds showed better activity against Gram-positive bacteria than against Gram-negative bacteria (Fig. S6).
Moreover, to determine the antibacterial effects against the pathogenic bacteria (Fig. S7), the minimal inhibitory concentration (MIC) was measured, along with the minimal bactericidal concentration (MBC) values of the most active thymidine analogs. The results are listed in Table S2 and Figs. 3 and 4. The best antibacterial effects against the tested strains were found for thymidine analogs 2 and 3, which showed MIC values in the range of 0.125-8.0 µg/mL. Analog 2 was active against all of the test bacteria, and the best activity for this compound was recorded against B subtilis (0.125 µg/mL). Interestingly, analog 3 recorded good activity against most of the test pathogens in impressively low concentrations, and the best activity was recorded against S. typhi (0.5 µg/mL). The lowest value of MBC was found for both analogs 2 and 3 (8.00 µg/mL) against B. subtilis and S. typhi. Besides, the highest MBC value (16.00 µg/mL) was obtained for these derivatives against E. coli, B. subtilis and P. aeruginosa. The MBC values for these compounds killing the other tested organism, lie between these ranges (8.00-16.00 µg/mL).

Antifungal susceptibility
Most of the methyl thymidine (1) analogs were found to have outstanding inhibition to the mycelial growth of both A. niger and A. flavus (Table 4). Among the tested analogs, 2 inhibited 82.76 ± 1.0% against the A. flavus and 60.16 ± 1.2% against the A. niger in their antifungal potential assessment. Remarkable mycelial growth prevention also built up for analog 3 against the A. niger (65.25 ± 1.1%) and A. flavus (72.95 ± 1.0%) in their mycelial growth test. Moreover, promising mycelial growth prevention also built up for analog 5 against the A. niger (57.62 ± 1.1%) and A. flavus (81.55 ± 1.0%) in the mycelial growth test. It is reaveled that analogs 2, 3 and 5 were very much effective against A. niger and A. flavus, and their zone of inhibition was higher than standard antibiotic nystatin (Figs. S8 and S9). Furthermore, analog 6 exhibited high inhibition against A. flavus, while no inhibition was observed against A. niger. The only synthesized analog 4 was inactive against both the pathogens; likewise, the parent molecule thymidine 1. It was found that the acylation of thymidine improves antimicrobial activity. The observed results reveal that the insertion of different acyl moieties, including 4-tert-butylbenzoyl, lauroyl, palmitoyl and 4-bromobenzoyl groups, significantly enhanced the antimicrobial activity of thymidine analogs.

SAR study
The development of antimicrobial agents is an area of great urgency and activity due to the emergence of multidrug resistance in common pathogens, and the appearance of new infections. The majority of antimicrobial agents are diverse five-and six-membered heterocyclic molecules that play a crucial role in the metabolism of all living cells [59]. Moreover, a great deal of interest has been directed toward condensed ring systems due to their various types of physiological activities and the success in utilizing them as privileged medicinal scaffolds. Considering that nucleosides play an essential role in most of the fundamental cellular metabolic functions, it is not surprising that nucleoside analogs are capable of targeting a variety of enzymes, such as those involved in bacterial peptidoglycan biosynthesis, fungal chitin biosynthesis, and those involved in protein synthesis. SAR analysis is important in understanding the mechanisms of antibacterial activity for thymidine derivatives (Fig. 5). The structure-activity relationship (SAR) of thymidine analogs can be established from the results of the antimicrobial activities reported in Table 3 and 4. The thymidine itself showed no activity against pathogenic bacteria, so a change in the thymidine skeleton greatly affected the antibacterial activity. It is possible to see that there is a tendency for decreasing influence in the following order: fused butylbenzoyl > lauroyl > palmitoyl for Gram-positive bacteria and lauroyl > palmitoyl > butylbenzoyl for Gram-negative bacteria respectively (Fig. 6).
For most of the bacteria tested, fused butyl benzoyl moieties were more active than the palmitoyl, but myristoyl and bromobenzoyl containing derivatives were weaker inhibitors than lauroyl analogs. In the present study, higher concentrations of synthesised compounds were required to inhibit the growth of Gram-negative bacteria (MIC values in the 2 to 8 mg/L) in comparison to Gram-positive. This type of different behaviour is probably due to the different structures of the cell wall between Gram-positive and Gram-negative bacteria. An explanation for this could be the fact that the peptidoglycan of Gram-negative bacteria is surrounded by an outer membrane which restricts diffusion through its lipopolysaccharide (LPS) covering. The LPS layer plays an essential role in providing selective permeability [60,61] and serves as an efficient barrier against the rapid penetration of various compounds. Moreover, the periplasmic space of Gram-negative bacteria contains hydrolytic enzymes, which can inactivate foreign molecules introduced from the environment. On the other hand, Gram-positive bacteria lack the outer membrane but are characterised by a thick, hydrophilic, porous structure that makes them more permeable. Therefore, Gram-positive bacteria are expected to be more sensitive to synthesized derivatives compared to Gram-negative bacteria. Moreover, the introduction of the lauroyl and palmitoyl groups gradually increased the hydrophobicity of the thymidine derivatives. The hydrophobicity of the compound is essential for bioactivity, as it may affect the toxicity or alter the integrity of the membrane because it is directly related to membrane permeation [10]. Hunt [34] and Judge [62] suggest that the potency of aliphatic alcohols is directly related to their lipid solubility due to the hydrophobic interactions between the alkyl chains of the alcohol and the lipid membrane regions. Hydrophobic interaction might occur between acyl chains of thymidine accumulated in the lipid-like nature of bacteria membranes. As a consequence of their hydrophobic interaction, bacteria lose their membrane permeability, ultimately causing the death of the bacteria.

Assessment of antimicrobial activities (PASS) and bioactivity
The antimicrobial spectrum was also predicted by applying the web server PASS to all the thymidine analogs 2-6. The PASS results are expressed as Pa and Pi, and are presented in Table 5. It can be seen in Table 5 that the thymidine analogs 2-6 showed 0.30 < Pa < 0.41 for antibacterial, 0.27 < Pa < 0.32 for antifungal, 0.51 < Pa < 0.72 for antiviral and 0.69 < Pa < 0.83 for anti-carcinogenic. These results reveal that these molecules were more efficient against bacterial and virus pathogens in comparison with fungal pathogens. The attachment of additional aliphatic acyl chains (C12 to C16) increased antibacterial activity (Pa = 0.415) of thymidine (1, Pa = 0.332), whereas the insertion of Br-and -C(CH 3 ) 3 substituted aromatic groups decreased the activity somewhat. The same scenario was observed for antifungal and anticarcinogenic activity, where acyl chain analogs revealed improved values compared to the benzoyl analogs. However, analog 2, which contain the lauroyl group, exhibited the highest antifungal and anticarcinogenic activities. An attempt was also made to predict the antiviral parameters of these analogs. Therefore, PASS determination exhibited 0.51 < Pa < 0.72 for antiviral, which shows that the thymidine analogs have more potential as antiviral agents than other antimicrobial parameters.

Thermodynamic analysis
A simple alteration of the chemical structure significantly influences structural properties, including thermal and molecular orbital properties. The spontaneity of a reaction and the stability of a product can be calculated from the free energy and enthalpy values [63]. Highly negative values are more likely to gain thermal stability. In drug design, hydrogen bond formation and nonbonded interactions are also influenced by the dipole moment. Free energy (G) is a significant criterion to represent the interaction of binding partners, where negative values are favorable for spontaneous binding and interaction. In the current study, the thymidine analogs possessed greater negative values for E, H and G than the parent thymidine, thus indicating that the attachment of an acyl group could improve interaction and the binding of these molecules with different microbial enzymes. The highest free energy is (-4316.309 Hartree) observed for thymidine analog (6), which has also shown the highest enthalpy (-4316.195 Hartree) and highest electronic energy (-4316.194 Hartree). The dipole moment value of a molecule is very significant to describe its electronic properties, where a high dipole moment value of a molecule causes more intermolecular interactions. A high dipole moment value reveals a more polar nature [64]. As shown in  The data are presented as mean ± SD and the values are represented for triplicate experiments. Statistically significant inhibition (p < 0.05) is marked with an asterisk (*) for test compounds and a double asterisk (**) for the reference antibiotic azithromycin. NI = No inhibition. Table 6, two thymidine analogs (4 and 6) have an improved dipole moment, which enhances the polar nature of a molecule, and promotes binding affinity, hydrogen bonding, and nonbonding interaction with the receptor protein. The highest dipole moment (5.911 Debye) was observed for analog (6), whereas (3) showed the lowest value (3.229 Debye). The scores for all parameters gradually increased with the number of carbon atoms in the substituents (2-6). The bromobenzoyl analog had better scores for all parameters, as evidenced by analog (6) (4-Br.Bz), which had the highest free energy of the therapeutics under investigation and showed markedly improved dipole moment. Overall, it is reasonable to propose that the modification of hydroxyl (-OH) groups of thymidine significantly increase its thermodynamic properties.

Frontier molecular orbitals (FMO)
The frontier molecular orbitals are the principal orbitals in a molecule, and they are considered to study chemical reactivity and kinetic stability. These frontier molecular orbitals are known as the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). The electronic absorption relates to the transition from the ground to the first excited state, and is mainly represented by one electron excitation from HOMO to LUMO [65].
The extended energy gap of a molecule relates to high chemical stability and low chemical reactivity. A low energy gap is related to low chemical stability and high chemical reactivity because of the ease of the transition of electrons. As a result, the removal of electrons from ground state HOMO to excited state LUMO requires more energy.
The HOMO and LUMO energies, HOMO-LUMO gap (Δ), hardness (η), softness (S), chemical potential (µ) and electronegativity (χ) and electrophilicity (ω) index of all analogs are presented in Table S3. Table S3 and Figs. 6, S10 and S11 show the thymidine analog (3) had a slightly higher energy gap value (5.575 eV), and the thymidine analog (6) had a slightly lower energy gap value (5.122 eV) than that of the other analogs. The thymidine analog (3) had chemical hardness and softness values of 2.788 eV and 0.359 eV respectively, where the hardness value was highest from among all the analogs. On the other hand, the thymidine analog (6) had the lowest chemical hardness (2.561 eV) as well as the highest chemical softness (0.391 eV). Analog (6) had a Brfunctional group in the aromatic ring, which is the reason behind this high chemical reactivity.

Molecular electrostatic potential map analysis
In computer-aided drug design, atomic charges are employed to investigate the connectivity between the structure and biological activity of a drug. The molecular electrostatic potential (MEP) is used globally as a reactivity map displaying the most suitable region for the electrophilic and nucleophilic attack of charged point-like reagents on organic molecules [66]. It helps to illustrate the biological recognition process and hydrogen bonding interaction [67].
An MEP counter map is a simple way of predicting how different geometry could interact. The MEP of the title analog is obtained based on the B3LYP with the basis set at 3-21G for an optimised result, as shown in Fig. 7. The importance of MEP lies in the fact that it simultaneously shows the molecular size and shape, as well as positive, negative and neutral electrostatic potential regions in terms of colour grading, which is very useful in the research of molecular structure, along with the physicochemical properties relationship [68]. Molecular electrostatic potential (MEP) was calculated to forecast the reactive sites for the electrophilic and nucleophilic attack of the optimized structure of thymidine (1) and analogs (2)(3)(4)(5)(6). The different values of electrostatic potential are represented by different colours, with potential increases in the order red < orange < yellow < green < blue. The red colour displays the maximum negative area, which shows favorable sites for electrophilic attack; the blue colour indicates the maximum positive area favorable for the nucleophilic attack, and the green colour represents zero potential areas (Figs. 8 and 9).

Molecular docking simulation
Molecular docking is a type of bioinformatics modeling which deals with the interaction between two or more molecules to give a stable adduct. The key aim of molecular docking is to predict the potential binding geometries of a putative ligand of a known three-dimensional structure with a target protein. In this investigation, a series of thymidine analogs was studied in silico to highlight their possible binding energy and interaction modes with the active site of SARS-CoV-2 main protease (Tables 7, 8 and S4), using the AutoDock Vina software. The multiple sequence alignment of the closest homologs of the main protease 6LU7, 6Y84 and 7BQY structure is displayed in Fig. S12. According to the outcomes obtained from docking screening, three analogs (2, 3 and 6) with the strongest binding energies were selected to describe the binding mode of the potential thymidine analogs. As shown in Table 7, the aromatic analogs displayed better binding scores than the aliphatic analogs. The interactions between the inhibitor and bordering residues of the SARS-CoV-2 main protease are illustrated in the 2D schematics, which were attained by importing the docking results into the Discovery Studio Visualizer. Fig. 8 and S13 showed the amino acids that participated in the pattern of interactions between the ligand and enzyme with an important contribution to the total energy of the interaction. Most of these interactions included hydrophobic contacts, Van der Waals interactions, hydrogen bonds, electrostatic, carbonyl, and one specific atom-aromatic ring, as well as providing an insight into understanding molecular recognition. Fig. 9 and S14 depicts the docked conformation of the most active molecules (2, 3 and 6) based on docking studies. The results show that thymidine analogs (3 and 6) are the most promising ligand (-7.9, − 6.3, − 7.2, − 8.6, − 6.6 and − 7.5 kcal/mol), which bound with SARSCoV-2 main protease via many hydrophobic bonding and hydrogen interactions. The binding sites were mainly located in a hydrophobic cleft bordered by the amino acid residues Cys145, His41, Phe294, Met49, Met165, Val104, Pro293, Thr199 and Thr25. There are six hydrogen bond contacts with four different amino acids: Asp153, Cys145, Thr26, Ser46, His164, Leu287 and Gln110 respectively. The thymidine analogs (3 and 6) had an additional benzene ring in the thymidine, providing a high density of electrons in the molecule, indicated by the highest binding score. These results show that modification of -OH group along with an aromatic ring molecule increased the binding affinity, while the addition of hetero groups like -C(CH 3 ) 3 , -Br caused some fluctuations in binding affinities; however, modification with halogenated aromatic rings increased the binding affinity. The docked pose clearly shows that the drug molecules bind within the active site of the SARS-CoV-2 main protease macromolecular structure.
The parent molecule, thymidine (1), exhibited interactions with the key residues of the main protease Cys145 and HIS41 through hydrogen and hydrophobic bonding within a closer bond distance (2.7547 Å). Additionally, Met165, Glu166 and Ser144 interactions were found; whereas Glu166 revealed a shorter distance (2.7547 Å) because of the unique interaction of branched alkyl chain with cytosine base. Aliphatic chain substituted analogs (2 and 4-5) revealed a lower binding score with the main protease, which indicated the burying of the ligand in the receptor cavity. In spite of having lower binding affinity, they also interacted with the catalytic binding of the main protease, such as Cys145, His41, Arg22, His164, Phe294, Cys44, Asn23, Glu166, Met49, Arg298, Thr111, Asp153 and Gly143. These analogs exhibited diverse nonbonding interactions, such as pi-alkyl, pi-sigma, amide pi-stacked, and pi-donor hydrogen bonds with the active site of the main protease.
These outcomes clearly show that, due to having high electron density, aromatic substituents can easily increase the binding ability as well as the antimicrobial ability of the thymidine analogs. Along with Phe294, all the analogs displayed the maximum π-π interactions with the His41, denoting tight binding with the active site. Reports suggest that Phe294 is thought to be the principal component of the PPS, APS and PPT responsible for the accessibility of small molecules at the active site. Binding energy and binding mode were exalted in some of the analogs (2, 3, and 6) because of significant hydrogen bonding. It was observed that the alterations of the -OH group in thymidine increased the π-π interactions with the amino acid chain at the binding site, whereas their polarity improvement resulted in the formation of hydrogen bond interactions. The maximum numbers of H-bonds were observed for analog (2), with Cys145, Gly143, and Glu166 residues.
It has already been reported that ten commercial medicines may possibly form H-bonds with key residues of 2019-nCoV main protease [69]. H-bonds execute a vital function in shaping the specificity of ligand binding with the receptor; drug design in chemical and biological processes, and molecular recognition and biological activity. The H-bond surface and hydrophobic surfaces of analogs (2 and 5) with the main protease are therefore represented in Fig. 10. It was accomplished from the docking study of all thymidine analogs with the SARS-CoV-2 main protease that the molecules are generally surrounded by the abovementioned residues, like the standard drug, this suggests that this molecule may prevent the viral replication of SARS-CoV-2.
The distance of the ligands, along with the change in the accessible area of the two important catalytic residues (Cys145 and His41) within the active site of the protease, is shown in Table 8. It is noted that the analyzed thymidine analogs bind within the active site of the main protease of SARS-CoV-2, like the standard drug hydroxychloroquine (HCQ), which is inevitable to prevent the protein mutarotation of the virus by minimizing the viral replication. The binding affinities calculated varied in the range of (-4.0 to − 8.6 kcal/mol), suggesting that the synthesized molecules can spontaneously interact within the binding site of the SARS-CoV-2 main protease. Although the blind docking studies revealed that all the molecules can act as potential agents for the establishment of COVID-19 drugs, from the estimated free energy of binding values, it may be inferred that the analogs (2, 3 and 6) with the highest negative minimum binding energy value amongst all the studied analogs could be the best possible SARS-CoV-2 inhibitors.

Molecular dynamics
The molecular dynamics simulation study was conducted to understand the structural stability of the docked complexes. Therefore, the root means square deviations of C-alpha atoms from the simulation trajectories were explored to examine flexibility across the simulation periods. Fig. 11 (a-f) indicates that analogs 2, 3, and 4 had the upper trend from the very beginning, which indicates the flexible nature of the docked complexes. Therefore, the complexes reached a steady-state after 35 ns times and maintained a lower degree of deviations in RMSD. The complexes had an RMSD below 2.5 Å for the entire atomistic simulation times, which defines the structural integrity of the complexes (Fig. 11a). Also, the solvent-accessible surface area of the complexes was examined to understand the changes in the protein volume, where a higher SASA correlates with the expansion of the surface area, and a lower SASA defines the truncated nature of the protein. Fig. 11b indicates that all complexes extended their surface volumes upon ligand binding during the initial phases, but complexes 2 and 3 had relatively stable profiles and maintained a similar trend. Moreover, analog 6 exhibited lower SASA than the other complexes, which defines the condensed nature of the complexes. The radius of gyrations of the complexes defines their labile nature, where a higher Rg is related to greater mobility, and a lower Rg indicates the stable nature of the complexes. Fig. 11c indicates that all complexes had a similar Rg profile; also, lower degrees of the deviations were observed for the three simulated complexes. This trend in the radius of gyration profiling defines the strict and stable nature of the complexes. The hydrogen bond plays an important role in defining the stable nature of the complexes, where Fig. 11d indicates that the three complexes lowered fluctuations and maintained integrity in hydrogen bond profiling. The root means square fluctuations were also analyzed to understand the flexibility across the amino acid residues. Fig. 11e demonstrates that almost every residue had a lower RMSF than 2.5 Å, which indicates the stable nature of the three complexes. Furthermore, the binding free energy of the complexes was explored from the MM-PBSA methods, where a positive score defines more favorable bindings according to the algorithms of YASARA dynamics (Fig. 11f). The average MM-PBSA scores of analogs 2, 3 and 4 were 25.47, 45.488 and 3.054 KJ/mol respectively, which indicates the favorable bindings of the three complexes.

Cytotoxic activity of the thymidine analogs
The cytotoxic activity of the synthesized thymidine analogs (2)(3)(4)(5)(6) in the brine shrimp lethality bioassay method [57] is displayed in Fig. 12, which exhibited the percentage of mortality of shrimps at 24 h and 48 h. A long alkyl chain increased the cytotoxicity of the compounds and the introduction of a phenyl ring was suggested to increase hydrophobicity and cytotoxicity [70]. From the observed data, it was found that    (6) were found less toxic for brine shrimp (i.e. 41. 25% and 38.75% death). It was evident from the observation, that the benzoyl analogs had lower cytotoxic activity than the alkyl chain analogs. Moreover, the cytotoxic activity of the alkyl chain derivatives increases with the concentration.

The pharmacokinetic profile and drug-likeness analysis
In an attempt to predict pharmacokinetic properties such as absorption, distribution, metabolism, excretion (ADME) and toxicity of the compounds, the pkCSM ADMET descriptors algorithm protocol was used. Drug absorption depends on factors such as membrane permeability [indicated by the cell line of colon cancer (Caco-2)], intestinal absorption, skin permeability thresholds, and substrate or inhibitor of Pglycoprotein. A value of intestinal absorbance below 30% suggests poor absorbance. According to the results shown in Table 9, all the analogs showed excellent absorption as the absorption values are more than 30. Skin permeability is an important consideration for improving drug efficacy, and it is particularly of interest in the development of transdermal drug delivery. A molecule will barely penetrate the skin if log Kp is more than − 2.5 cm/h. From Table 9, it can be seen that the skin permeability (Kp) of the thymidine analogs is − 2.732 to − 2.829 cm/h (<-2.5). Therefore, it can be predicted that all analogs have good skin penetrability. For the pkCSM predictive model, high Caco-2 permeability is translated into predicted log Papp values greater than 0.90 cm/ s. As Table 9 shows, the value of the Caco-2 permeability (log Papp) of the thymidine analogs ranged from 0.47 to 0.78 cm/s, log Papp < 0.9 cm/s, so it is predicted that these have low Caco-2 permeability.
For the discovery of oral administrative drugs, solubility is one of the major descriptors. High water solubility is useful for delivering active ingredients in sufficient quantities in small volumes of such pharmaceutical dosage. The values for water solubility are given in log (mol/l) (Insoluble ≤ -10 < poorly soluble < -6 < Moderately < -4 < soluble < -2 < very soluble < 0 ≤ highly soluble). From the results shown in Table 9, it can be observed that the analogs tested are soluble.
Distribution volume (Vd) is a pharmacokinetic parameter that reflects the tendency of an individual substance to either linger in the plasma or redistribute to another tissue compartment. According to Pires et al. [58], VDss is considered low if it is below 0.71 L/kg (log VDss < -0.15), and high if it is above 2.81 L/kg (log VDss greater than 0.45. It can be seen in Table 10 that the value of thymidine analogs VDss ranged from − 0.551 to 0.433, with only one analog having a VDss value of < -0.15 (analog 3). The blood-brain partitioning and brain distribution are critical properties for drugs targeting the central nervous system. The analogs tested a logBB < -1 considered poorly distributed to the brain. From Table 10, it can be seen that the logPS (the central nervous system (CNS) permeability) value of thymidine analog range from -3.74 to − 3.25, logPS < -3 so it can be predicted that analogs (2)(3)(4)(5)(6) are unable to penetrate the CNS. The model provided by pkCSM pharmacokinetics predicts the total clearance log(CLtot) of a given compound in log(mL/ min/kg). The larger the CLtot value of the compound, the faster the excretion processes. It can be seen in Table 10 that the log CLtot value of thymidine analogs ranges from 0.02 to 1.65 ml/min/kg, and from those values, the rate of excretion of the analog can be predicted. Metabolism is predicted based on the CYP models for substrate or inhibition (CYP2D6, CYP1A2, CYP2C19 and CYP3A4). Table S5 shows that all the thymidine analogs do not affect or inhibit all the enzymes, except CYP3A4 for analog 6, so it can be predicted that analogs (2)(3)(4)(5) in the body tend to be metabolized by the P450 enzyme. Bioactivity radar Charts of the thymidine analogs (Fig. S15) reveal the promising pharmacokinetic profiles of all candidates. The toxicity results of the thymidine analogs are described in Table S6, and their high LD50 values (2.23 to 2.47) suggest that the analogs are lethal only at very high doses. A negative result in the AMES test suggests that an analog is not mutagenic. The results also suggest that all analogs tested may not inhibit the hERG channel and may not have skin sensitization. Molecular parameters like the membrane permeability of lead molecules depend on some fundamental features of molecules, such as partition coefficient (logP), molecular weight (MW), and several hydrogen bond acceptors/donors, which is associated with Lipinski's 'rule of five' [71]. Table S7 showed that analogs (3)(4)(5)(6) violate the rule of five in two cases, indicating the good bioavailability of thymidine analogs, whereas analog (2) followed Table 8 Nonbonding interaction data of thymidine analogs against main protease.  the rules of Lipinski. Low molecular weight containing derivatives are easily absorbed, diffused and transported, as compared to high molecular weight derivatives. As molecular weight increases past a certain limit, the bulkiness of the molecules also increases comparably. TPSA (Topological Polar Surface Area) is a very important physicochemical property of molecules that provides information on the polarity of analogs. This parameter was screened to analyse drug transport properties. Polar surface area is the sum of all polar atoms, mainly oxygen and nitrogen, and includes attached hydrogen.
3.14. POM analysis: Identification of the pharmacophore sites POM Theory (Petra/Osiris/Molinspiration) that is invented by the group of Taibi Ben Hadda in collaboration with the American NCI and TAACF led us to real success in pharmacology and drug design fields [72]. By using the POM (Petra/Osiris/Molinspiration) Theory, it becomes now easier, to identify and optimize most of the antibacterial [73][74][75], antifungal [76][77][78], antiviral [79][80][81], antiparasitic [82,83], and antitumor [85] pharmacophore sites, one by one, on the basis of their different physicochemical parameters and their different electronic charge repartition of corresponding heteroatoms. This young POM Theory was extended, with success, to other various and different bio targets [86]. Here we treat the series of compounds 2-6 to identify their pharmacophore sites. So the identification of the type of pharmacophore sites of these compounds was derived from the physical and chemical properties of the tested compounds by using Petra, Osiris and Molinspiration (POM) platform. Results of pharmacokinetic properties and bioactivity score analysis are shown in (Table 11).

Conclusion
In this research, a series of thymidine analogs were analysed in vitro and in silico for their antimicrobial, thermodynamic, molecular docking, molecular dynamics and drug-likeness properties. Most of the thymidine analogs have a HOMO-LUMO gap lower than thymidine, resulting in modified analogs being more reactive than the parent drug. Insertion of various aliphatic and aromatic groups into the thymidine structure can significantly improve their biological activity mode. The study revealed that 4-t-butyl and bromo substituted benzoyl analogs 3 and 6 had more potential against bacterial organisms with better pharmacokinetic and biological spectra. These observations were rationalised by conducting molecular docking, which revealed promising antiviral efficacy for thymidine analogs. Thymidine analogs (2, 3 and 6) showed promising binding interactions and binding energy with SARS-CoV-2 main protease which revealed analogs (2, 3 and 6) showed an in silico potent ability to fight SARS-CoV-2. The molecular electrostatic potential study also showed the most negative and positive surface area of the investigated ligand, and thus, anticipated the most suitable areas for hydrogen bonding sites. This result greatly supports molecular dynamic studies up to 100 ns, keeping in mind protein, which confirms the binding stability of the docked complex in the trajectory analysis. This also means that the protein-ligand complex is highly stable in biological systems. In addition, these analogs were analysed for their cytotoxicity and pharmacokinetic properties, which revealed that a combination of toxicity determination, in silico ADMET prediction, and drug-likeness has promising results because most of the designed molecules showed improved kinetic parameters, and maintained all drug-likeness rules as well as provided an interesting result in terms of biological activity. As this study has been carried out using synthetic, antimicrobial and in silico computational methods, these results require further wet-lab experiments to be carried out under in vivo as well as in vitro conditions for these analogs to be potential drug candidates to treat SARS-CoV-2. This hypothesis is also predicted by POM Theory which shows the presence of an antiviral (O1 δ-, O2 δ-) pharmacophore site.
Author Contributions. MAH and NSM synthetic experimental analysis; MAG and MB validation of article and improve the manuscript; FAA, TBH and HL performed the POM study; FA antimicrobial; SM and MAS performed the computational, dynamics and interpretation of the results; SMAK conceptualization, methodology, article writing, and supervision. All authors reviewed and approved the manuscript.