key: cord-0911936-o9wm5b0b
authors: Matin, Mohammed Mahbubul; Uzzaman, Monir; Chowdhury, Shagir Ahammad; Bhuiyan, Md. Mosharef Hossain
title: In vitro antimicrobial, physicochemical, pharmacokinetics and molecular docking studies of benzoyl uridine esters against SARS-CoV-2 main protease
date: 2020-12-10
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1850358
sha: 3d3eee1d25b2d99123c81ef993689c6fed6334c8
doc_id: 911936
cord_uid: o9wm5b0b

Different esters were found potential against microorganisms, and could be a better choice to solve the multidrug resistant (MDR) pathogenic global issue due to their improved biological and pharmacokinetic properties. In this view, several 4-t-butylbenzoyl uridine esters 4–15 with different aliphatic and aromatic groups were synthesized for antimicrobial, physicochemical and biological studies. In vitro antimicrobial tests against nine bacteria and three fungi along with prediction of activity spectra for substances (PASS) indicated promising antifungal functionality of these uridine esters compared to the antibacterial activities. In support of this observation their cytotoxicity and molecular docking studies have been performed against lanosterol 14α-demethylase (CYP51A1) and Aspergillus flavus (1R51). Significant binding affinities were observed against SARS-CoV-2 main protease (7BQY) considering hydroxychloroquine (HCQ) as standard. ADMET predictions were investigated to evaluate their absorption, metabolism and toxic properties. Most of the uridine esters showed better results than that of the HCQ. Overall, the present study might be useful for the development of uridine-based novel MDR antimicrobial and COVID-19 drugs.

Uridine (1), a natural and necessary nucleoside, is one of the four basic components of ribonucleic acid (RNA). It is a structural element of many bioactive compounds, and it has crucial role in the pyrimidine metabolism of the brain. Uridine and its nucleotide derivatives may also have an additional role in the function of the central nervous system as signaling molecules (Dobolyi et al., 2011) . It has anti depression activity, asthmatic airway inflammation, hepatocyte proliferation (Carlezon et al., 2005) . Uridine has sleep-promoting and anti-epileptic effects, might affect mood, improves memory function and influences neuronal plasticity (Siegel et al., 2017) . Uridine, sometimes, has been recognized as playing a role in preserving and enhancing hepatic mitochondrial function in the face of challenge by certain liver toxins (Banasch et al., 2006) . In addition, uridine (1) and some of its biological derivatives behave as antimetabolites, when their extracellular concentrations are higher than those physiological levels (Strasser et al., 2006; Yamamoto et al., 2011) .

Many advantageous findings were reported for the uridine derivatives especially its esters. For example, uridine 5 0 esters showed neuroprotective effect and is capable of compensating the harmful effect of neurotoxins with the increase the 5-HT and 5-HIAA levels (Susilo, 2008) . 3 0 ,4 0 -Uridine acetonide (2, Figure 1 ) showed a differential effect and a higher activity over cell viability on MCF-7 human cancer breast cell line than for the CHO-K1 cell line (Escober et al., 2016) . Some uridine esters were found to exhibit promising antimicrobial properties Kabir, Matin, Kawsar, & Anwar, 1998; Kabir, Matin, Mridha, et al., 1998; Kabir et al., 1998a Kabir et al., , 1998b such as 3 0 -O-decanoyl-2 0 ,5 0 -di-O-acetyluridine (3) was found very active against E. coli (Kabir et al., 2002) . Thus, it has been proven that the uridine derivatives especially acyl esters can be used as drugs in order to treat, for instance, Parkinson, Alzheimer, Chorea Huntington or depression with cellular regulatory functions (Escober et al., 2016) .

Although uridine analogues (e.g. doxifluridine) have shown higher activities in the treatment of viral and cancer treatment, their efficiency is sometimes reduced by the appearance of resistance mechanisms (De Clercq, 2011) . Hence, search for newer analogues of uridine (nucleoside) got more attention especially esterification of its sugar moiety (Du et al., 2018) . However, investigations indicated that selective esterification (acylation) uridine (1) creates difficulties due to the presence of one 1 OH group, two 2 OH groups and a NH group. Hence, various methods like direct acylation (Matin, Bhattacharjee, et al., 2019; Matin, Bhuiyan, et al., 2019; Matin, Chakraborty, et al., 2020; Matin, Hasan, et al., 2020; , dibutyltin oxide (DBTO) method Kabir, Matin, Kawsar, & Anwar, 1998; Kabir, Matin, Mridha, et al., 1998; Kabir et al., 1997 Kabir et al., , 1998a Kabir et al., , 1998b , coupling technique (Steglich esterification) (Clayden et al., 2012) , protection-deprotection (Matin et al., 2017) and enzymatic method (Du et al., 2018) were reported. Among these methods, direct method generally furnished 5 0 -O-substituted esters Kabir, Matin, Kawsar, & Anwar, 1998; Kabir, Matin, Mridha, et al., 1998; Kabir et al., 1997 Kabir et al., , 1998a Kabir et al., , 1998b whereas DBTOmethod gave 3 0 -O-substituted products Kabir, Matin, Kawsar, & Anwar, 1998; Kabir, Matin, Mridha, et al., 1998; Kabir et al., 1997 Kabir et al., , 1998a Kabir et al., , 1998b Kabir et al., , 2002 . However, in most of the cases direct method gave higher yield(s) with lower cost (Matin, 2014; Matin et al., 2017) . Chowdhury et al. (2016) applied direct method with hindered 4-t-butylbenzoyl chloride in pyridine at À5 C and afforded the 5 0 -O-(4-t-butylbenzoyl)uridine (4) in excellent yield. Now-a-days, emergence of multiple drug resistant (MDR) pathogenic organisms, and corona virus disease 2019 (COVID-19) are global concern (Ghosh et al., 2020; Lakshmi et al., 2020; Wahedi et al., 2020) in addition to numerous obstacles encountered on delivery to the site of microbial infection. Many researchers are engaged to overcome these situations (Awual, 2017; Laws et al., 2020; Matin, 2008; Sharavanan et al., 2020) . However, there is urgent need for establishment of new chemotherapeutic agents with novel mode of action to combat with these MDR organisms, and SARS CoV-2 (severe acute respiratory syndrome corona virus-2). In this regard as well as our continuous effort to find novel drugs (Matin, 2014; Matin & Ibrahim, 2010; Matin et al., 2015) we report herein the in vitro antimicrobial efficacy of several uridine-based esters 4-15 with different aliphatic and aromatic chains against eleven pathogens, molecular docking against lanosterol 14a-demethylase (3LD6), A. flavus (1R51) and corona virus (SARS CoV-2) main protease (Mpro; 7BQY) along with prediction of activity spectra for substances (PASS). Also, attempts have been taken to optimize some acylated uridine derivatives to investigate their physicochemical behavior on the basis of DFT approach with cytotoxicity, and ADMET properties.

In our present study, twelve uridine-based esters with different side chain length and different groups 4-15 were used as test chemicals. These compounds (mainly 5 0 -O-t-butylbenzoyl esters of uridine), as presented in Figure 2 , were synthesized from 1 by direct method and well characterized earlier (Chowdhury et al., 2016) . For comparative study, uridine (1), standard antibacterial antibiotic ampicillin and antifungal antibiotic nystatin were also used.

Web-based PASS (http://www.pharmaexpert.ru/ PASSonline/ index.php) was used for the prediction of biological spectrum of these uridine derivatives (Matin et al., 2016) . Initially, structures of the uridine esters were drawn, and then, converted into their SD file formats which were used to predict biological spectrum using PASS online software. This program is designed to anticipate a plethora of biological activities with 90% accuracy. PASS results are designated as Pa (probability for active compound) and Pi (probability for inactive compound). Being probabilities, the Pa and Pi values vary from 0.000 to 1.000 and in general, Pa þ Pi 6 ¼ 1, since these probabilities are calculated independently. Only the activities with Pa > Pi are considered as possible for a particular drug. The PASS prediction results were interpreted and used in a flexible manner, viz. (i) when Pa > 0.7, the chance to find the activity experimentally is high, (ii) if 0.5 < Pa < 0.7, the chance to find the activity experimentally is less, but the compound is probably not so similar to known pharmaceutical agents and (iii) if Pa < 0.5, the chance to find the activity experimentally is less, but the chance to find a structurally. So, the predicted activity of spectrum of a drug is known as its intrinsic property.

The tested microorganisms (bacteria and fungi) were collected from the Pharmacology laboratory, Department of Physiology, Biochemistry & Pharmacology, Faculty of Veterinary Medicine, Chittagong Veterinary and Animal Science University, Chattogram, Bangladesh. As shown in Table 1 , nine human pathogenic bacteria and three plant pathogenic fungi were used for the antimicrobial evaluation.

In vitro antibacterial activities of uridine 4-15 were conducted by paper disc diffusion method (Kabir et al., 2003; Matin et al., 2013) , which was determined by the presence or absence of diameter of zone of inhibition (in mm). In this method, initially a 2% solution of the uridine compound(s) in DMF was prepared. Mueller-Hinton (agar and broth) medium was used to culture the bacteria. Paper disc of 4 mm in a diameter and petri plate of 70 mm in a diameter was used throughout the experiment. These paper discs were sterilized in an autoclave and dried at 150 C in an electric oven. Then, the discs were socked with test compound(s) at the rate of 100 lg (dry weight) per disc. These paper discs were put in the petri-plates and were incubated at 37 C for 48 h for growth of test organisms followed by the cooling for 4 h at low temperature (4 C). The test chemicals diffused from disc to the surrounding medium by this time. The growth inhibition was calculated at 24 h intervals for two days. Each experiment was carried out three times. All the results were compared with the standard antibacterial antibiotic ampicillin (50 lg/disc, brand name Decapen, Beximco Pharmaceuticals Ltd., Bangladesh). Ampicillin was also used as a positive control and proper control was maintained with DMF without chemicals under identical conditions.

The in vitro antifungal screening of the synthesized uridine derivatives (4-15) were investigated employing food poisoning technique Kabir, Matin, Kawsar, & Anwar, 1998; Kabir, Matin, Mridha, et al., 1998; Kabir et al., 1998a Kabir et al., , 1998b Matin, 2014) . Sabouraud potato dextrose agar (PDA, agar and broth) medium was used for the culture of fungi. Necessary amount of medium was taken in a conical flask separately and was sterilized in an autoclave (at 121 C and 15 psi) for 15 min. After autoclaving, weighed amount of test chemical (2% in DMF) was added to the sterilized medium in conical flask at the point of pouring to obtain desired concentration. The flask was shaken thoroughly to mix the chemical with the medium homogeneously before pouring. The medium with definite concentration (2%) of chemical was poured at the rate of 100 lg/mL PDA in sterilized glass petri dishes individually and subjected for incubation. Finally, linear mycelial growth of fungus was measured after 2-5 days (24 h interval) of incubation. The percentage of radial mycelial growth inhibition of the test fungus was calculated using the following equation: 

Where, I ¼ percentage of inhibition, C ¼ diameter of the fungal colony in control (DMF), T ¼ diameter of the fungal colony in treatment. In all the cases, the results were compared with standard antifungal antibiotic nystatin (100 lg/mL medium, brand name Fungistin, Beximco Pharmaceuticals Ltd., Bangladesh). Proper control was maintained separately and with sterilized PDA medium without chemical and three replications were prepared for each treatment.

Quantum mechanical (QM) methods has gained attention on calculation of thermodynamic properties, molecular orbital features, dipole moment and as well as interpretation of different types of interactions. Molecular geometry optimization and further modification of all uridine derivatives carried out using Gaussian 16 program. All the structures were optimized using density functional theory (DFT) employing Becke's (B) three-parameter hybrid model, Lee, Yang and Parr's (LYP) correlation functional under 3-21G basis set (Kruse et al., 2012) . Initial optimization of all compounds was performed in the gas phase. Dipole moment, electronic energy, enthalpy and free energy are calculated for all the compounds.

Frontier molecular orbital features HOMO (highest occupied molecular orbital), LUMO (lowest unoccupied molecular orbital) were calculated at the same level of theory. For each of the structure HOMO-LUMO gap(De), hardness (g), softness (S) and chemical potential (m) were calculated from the energies of frontier HOMO and LUMO as reported considering Parr and Pearson interpretation of DFT and Koopmans theorem (Pearson, 1986 ) on the correlation of ionization potential and electron affinities with HOMO and LUMO energy (e). Hardness (g), softness (S) and chemical potential (m) are calculated by using the following equations:

The 3D crystal structure of main protease of COVID-19 (7BQY), lanosterol 14a-demethylase (3LD6/CYP51A1) and A. flavus (1R51) was obtained in pdb format from online protein data bank (PDB) database (Lucido et al., 2016) . All hetero atoms, inhibitors and water molecules were erased using PyMol (version 1.3) software packages. Energy minimization of the protein implemented by Swiss-Pdb viewer software (version 4.1.0) (Matin, Chakraborty, et al., 2020; Matin, Hasan, et al., 2020; . The optimized structures were subjected for molecular docking study against 7BQY, CYP51A1 and 1R51 considering the protein as macromolecule and uridine esters as ligand. (Wahedi et al., 2020) . Docking results are basically calculated as binding affinities of ligand-protease and reported in kcal/mol unit.

For the determination of toxicity of the uridine derivatives we employed brine shrimp lethality assay technique as described by McLaughlin (1991) . The tested compounds were dissolved in DMSO, and prepared 20, 40, 80 and 160 lL by adding NaCl solution to each vial up to 1 mL volume which were designated as type-A, B, C and D, respectively. Three sets of experiment were done for each concentration and 10 brine shrimps nauplii were placed in each vial. A control experiment was performed in a vial containing 10 nauplii in 1 mL seawater. After 24 h and 48 h of incubation at room temperature, the vials were observed using a magnifying glass and the number of survivors in each vial was counted and noted. From the data, the average percentage of mortality of nauplii was calculated for each concentration. No deaths were found for controls.

Online computational approaches were made to predict absorption, distribution, metabolism, excretion and toxicity (ADMET) for 5 0 -O-(4-t-butylbenzoyl)uridine (4) and its 2 0 ,3 0 -di-O-acyl derivatives 5-15 (Matin, Chakraborty, et al., 2020; Matin, Hasan, et al., 2020; . The computational prediction of drug could avoid the tremendous cost and time associated with the in vivo experiments, and has attracted more and more attention. In the present study, we utilized in silico screening approaches employing pkCSM server (http://biosig.unimelb.edu.au). The pkCSM signatures were successfully used across five main different pharmacokinetic property classes to develop predictive regression and classification models (Pires et al., 2015) . Firstly, all the structures of uridine esters were drawn in ChemDraw 16.0 to collect InChI Key, isomeric SMILES (simplified molecular-input line-entry system), and SD file format. ADMET of all the uridine esters were predicted by using pkCSM-pharmacokinetics and SwissADME free web tools (http://www.swissadme.ch).

In the present study, twelve uridine esters with different aliphatic and aromatic chains (4-15) were selected as test compounds, and screened for their in vitro antibacterial and antifungal activities against nine bacteria and, and three phytopathogenic fungi, respectively. Initially, we predicted these activities using PASS program followed by in vitro determination. The observed activities were then rationalized by calculating their physicochemical (DFT method), cytotoxicity (shrimp), molecular docking and ADMET properties. Finally, all these findings were used to summarize their structure-activity relationship (SAR).

Having all the uridine esters 4-15, we predicted their biological spectrum using online PASS software (http://www. pharmaexpert.ru/ PASSonline/index.php). The PASS results in its designated as Pa and Pi form is presented in Table 2 . It was evident from predication Table 2 for uridine esters 4-15 showed 0.16 < Pa < 0.35 for antibacterial and 0.38 < Pa < 0.53 for antifungal. This clearly indicated that these compounds were more prone against fungi as compared to that of bacterial pathogens. Also, incorporation of mesyl groups, as in compound 6, totally abolished antifungal potentiality of uridine. Addition of aliphatic acyl chains (C2 to C16, compound 5, 7-12) increased antifungal activity (Pa ¼ 0.526) of uridine (1, Pa ¼ 0.421), whereas attachment of chloro-substituted benzoyl groups did not enhance reasonably (13-15, Pa ¼ 0.389-0.436). We have also calculated anti-carcinogenic and antioxidant property of these compounds. Thus, PASS predication showed 0.46 < Pa < 0.67 for anti-carcinogenic and 0.13 < Pa < 0.43 for antioxidant, which indicated that the uridine esters were more potent as anti-carcinogenic agents than that of antioxidant properties. Interestingly, the results as shown in Table 2 , antifungal and anti-carcinogenic properties of uridine esters with saturated aliphatic chains (5, 7-12) were found more promising than the sulphonyl ester (6) or chloro-substituted benzoyl derivatives (13-15). These results are in agreement with our previous study for monosaccharide esters (Matin, Chakraborty, et al., 2020; Matin, Hasan, et al., 2020; .

For the in vitro antibacterial efficacy test we followed disc diffusion method (Kabir et al., 2003) . The results of diameter of inhibition zone (mm) of the uridine-based esters 4-15 against four Gram-positive, and five Gram-negative bacteria are presented in Table 3 ( Figure 3 ) and Table 4 (Figure 4) . It was evident that these esters are less susceptible towards both the Gram-positive and Gram-negative pathogens, and similar results were reported for other uridine esters Kabir, Matin, Kawsar, & Anwar, 1998; Kabir, Matin, Mridha, et al., 1998; Kabir et al., 1998a Kabir et al., , 1998b . Amongst the uridine derivatives, compound 8 (C6 chain) and 10 (C12 chain) showed moderate inhibition against B. megaterium, S. dysenteriae and V. cholerae. Also, mesyl derivative 6 and chlorobenzoyl derivatives 13-15 showed very little antibacterial functionality which is in agreement with the PASS results (Table 2 ). In general, compounds 4-15 were found more prone against Gram-negative pathogens than the Gram-positive organisms.

The susceptibility of uridine esters 4-15 against three pathogenic fungi are presented in Table 5 . The results, as presented in percentage inhibitions of mycelial growth, indicated that these compounds were more prone against tested fungal pathogens compared to bacterial organisms. Among the tested fungi, 4-15 were more susceptible against M. phaseolina. Aliphatic chain containing compound 5 (C2, 57.8%, Figure showed excellent inhibition against M. phaseolina and the results were comparable to nystatin ( Ã 71.8%) in contrast to their antibacterial and PASS results.

The molecular structure significantly influences the structural properties such as free energy (G), enthalpy (H), dipole moment and electrostatic potential. Spontaneity of a reaction and stability of a product can be predicted from Gibb's free energy, enthalpy and electronic energy. Free energy (G) is a significant criterion to represent the interaction of binding partners, where negative value favorable for spontaneous binding and interaction. Greater negative values predict better thermodynamic properties. In the present study, all the uridine esters possess greater negative value for E, H and G than the parent uridine, and hence, indicated that the attachment of ester group could improve interaction and binding of these molecules with different microbial enzymes. Due to these higher values (-2990 to -3903 Hartree) chlorobenzoyl esters 14-15 exhibited better inhibition against M. phaseolina. In addition, these increased negative values of compound 4-15 suggested that these are energetically and configurationally more stable. As shown in Table 6 , some of the uridine esters have improved dipole moment. Improved dipole moment enhances the polar nature of a molecule, and promotes the binding affinity, hydrogen bonding and nonbonding interaction with the receptor protein (Lien et al., 1982) . The dipole moment of compound 4 was 7.588 Debye and for compound 15 was 7.619 Debye which were slightly higher than the uridine (7.339 Debye) resulting their better binding affinity and interactions with the amino acid residues of receptor protein.

The HOMO and LUMO energies, HOMO-LUMO gap (D? ), hardness (g), softness (S) and chemical potential (m) index of all compounds are presented in Table 7 ( Figure 6 ). According to the frontier molecular orbital theory, energies of HOMO and LUMO play an important role in chemical reactivity. The HOMO-LUMO gap is also related to the chemical hardness, softness, chemical potential and electrophilic index of a molecule (Cohen & Benson, 1993; Lien et al., 1982) . Large HOMO-LUMO gap is responsible for high kinetic stability and low chemical reactivity. On the other hand, small HOMO-LUMO gap is important for low chemical stability, because addition of electrons to a high-lying LUMO and/or removal of electrons from a low-lying HOMO is energetically favorable in any potential reaction. In these esters 4-15, D? value was found lower than uridine (1), and thus, favorable for energetically favorable reaction. We noticed that, with the increase of the number of ester group(s) and chain length (4-15) hardness of these compounds gradually decreased while their softness gradually increased. All these properties may contribute to show higher chemical activity and polarizability in the drug related chemical and biochemical functionalities.

Estimation of the feasible binding geometries and interactions between drugs and active site of proteins (Matin, Bhattacharjee, et al., 2019; Matin, Bhuiyan, et al., 2019) is obtained from molecular docking (Ghosh et al., 2020) . Molecular docking score of esters 4-15 along with standard antibiotic ampicillin (AMP) and hydroxychloroquine (HCQ) are shown in Table 8 . We calculated docking score (binding affinity) of 13-14, as these chlorobenzoyl compounds showed promising antifungal activities, with lanosterol 14ademethylase (CYP51A1) and A. flavus (1R51) while all the uridine esters were considered for COVID-19 main protease (7BQY). In addition, other uridine esters 4-12 expressed, in average, somewhat lower docking score than 13-14 against CYP51A1 and 1R51. Ampicillin (AMP) and HCQ were used as standard drugs; Discovery Studio (version 4.1) was performed for analysis and view of docking results.

Docking score clearly indicated that chlorobenzoyl esters 13-15 had binding energy with both the CYP51A1 (-7.9 to -12.0 kcal/mol; Figure 7(a,b) ) and 1R51 (-7.5 to -9.3 kcal/mol; Figure 7 (c,d)) proteins which were higher than the standard antibiotic ampicillin (-7.9 and -7.5 kcal/mol). The results are in complete accord with the in vitro results ( Table 5 ) that these uridine esters had better antifungal activities.

Considering the current COVID-19 pandemic caused by the SARS-CoV-2 betacoronavirus (Ghosh et al., 2020; Lakshmi et al., 2020) we have extended our studies by docking with COVID-19 main protease (7BQY). As shown in Table 8 , most of the uridine esters showed higher binding affinity with 7BQY than HCQ (-6.1 Kcal/mol). It is worth to mention that saturated fatty acyl esters with 5, 7-8 smaller chain length exhibited better binding affinity than the higher chain length 9-12, sulphonyl ester 6 and chlorobenzoyl esters 13-15 ( Figure 8 ). The docking score of 5, 7 and 8 (-7.9 to -8.0 kcal/mol) were found to better than the previously published other compounds (Ghosh et al., 2020) . Thus, uridine aliphatic esters could be better choice than HCQ in treating COVID-19 infection.

The toxicity activity of the acylated uridine derivatives 4-15 in the brine shrimp lethality bioassay method (McLaughlin, 1991) is presented in Figure 9 , which showed the percentage of mortality of shrimps at 24 h and 48 h. From the data, it was evident that the test chemicals 5 0 ˗O˗(4˗t˗butylbenzoyl)˗ 2 0 ,3 0 ˗ di˗ O˗ pentanoyl-uridine (7) and 5 0 -O-(4-t-butylbenzoyl)-2 0 ,3 0 -di-O-lauroyluridine (10) showed highest levels of toxicity (i.e. $57% death) indicating its higher mortality. The rest of the compounds were found less toxic for brine shrimp (Figure 9 ). Table 7 . Energy (eV) of HOMO, LUMO, Gap (D? ), hardness (g), softness (S) and chemical potential (l) of uridine (1) and its esters 4-15. -9.0 -7.3 -7.9 6 -8.5 -7.2 -6.6 7 -7.5 -7.5 -8.0 8 -9.2 -6.4 -8.0 9 -9.0 -6.0 -6.4 10 -5.4 -5.8 -5.5 11 -6.6 -6.3 -4.9 12 -4.5 -5.8 -5.2 13 -12.0 -9.3 -7.7 14 -11.1 -7.7 -7.6 15 -8.3 -8.0 -7.8 AMP -7.9 -7.5 -HCQ ---6.1 Figure 6 . Frontier molecular orbital and related energy of uridine derivatives 5 and 15.

3.6. ADME/T properties Absorption, distribution, metabolism, excretion and toxicity (ADMET) related pharmacokinetic profile is mentioned in Table 9 . It is evident that these compounds possess better absorption value than standard antibiotics (ampicillin and nystatin) although they are P-glycoprotein inhibitor (except 15). Distribution like BBB (>-1) and CNS (>-2) permeability values are comparable to standard drugs. Compounds 4-14 are CYP3A4 enzyme substrate while 15 is non-substrate for this enzyme. Their excretion, presented as total clearance, is comparable to ampicillin. We have also calculated their toxicity values along with standard antibiotics. hERG (human ether-a-go-go-related gene I) showed negative values for these compounds like nystatin. Hence, these compounds could be safe and will not generate QT syndrome like fatal ventricular arrhythmia (Matin, Bhattacharjee, et al., 2019; Matin, Bhuiyan, et al., 2019; Matin, Chakraborty, et al., 2020) . In addition, it is important to consider the lethal dosage (LD 50 ) values and is used to assess the relative toxicity of different molecules. Here, most of the ester's LD 50 values (<3 mol/kg) are in good consistent with the value of nystatin (2.52 mol/kg).

Drug likeness calculations were calculated from SwissADME tool and are summarized in Table 10 . For comparison both the uridine esters and standard antibiotics' C2P: Caco-2 permeability (log Papp in 10 -6 cm/s, >0.90 indicates high permeability); HIA: Human intestinal absorption (% absorbed, >30% is better absorbed); P-gpI: P-glycoprotein inhibitor; BBB (blood brain barrier) is expressed in logBB (logBB >-1.0 is moderately cross blood brain barrier); CNS (central nervous system) is expressed as logPS (logPS > -2.0 can easily penetrate the CNS); Total clearance is expressed in log mL/min/kg; hERG: human ether-a-go-go-related gene I; Toxicity is calculated in oral rat acute toxicity (mol/kg); AMP: ampicillin; NYS: nystatin.

drug likeness values were calculated. It was found that all the esters (except mesylate 6) have good hydrogen bonds donor and acceptor which is in consistent with the Lipnski's rule of five. The more the value topological polar surface area (TPSA) more the polarity of the compounds and should be less than 140 Å 2 . The TPSA in these molecules was in excellent agreement ($143 Å 2 ) with the most important rules of drug likeness (Lipinski et al., 1997) .

In the present study, we have checked the CYP enzymes, particularly isoforms 1A2, 2C9, 2C19, 2D6 and 3A4, inhibitory activities which are responsible for about 90% oxidative metabolic reactions. Inhibition of CYP enzymes will lead to inductive or inhibitory failure of drug metabolism. These inhibitions of compound 4, 9-12 and 15 were found similar to that of the standard drugs (Table 10 ). In addition, Panassay interference compounds (PAINS) revealed no violation with these uridine esters. PAINS are chemical compounds that often give false positive results in high-throughput screens. PAINS tend to react nonspecifically with numerous biological targets rather than specifically affecting one desired target. Many studies indicated the acute toxic effects of medicines on non-targeted organisms, and hence, the MTT, MTS and Resazurin assay methods are used to check the cytotoxic concentration of any new molecule (synthetic or natural) which gives an idea about the maximum nontoxic dose of the drug molecule. We hope to conduct such studies for the more drug related validation of these promising uridine compounds.

It was observed that the most of the potential antifungal drugs (fluconazole; itraconazole; voriconazole; terbinafine) are designed to deter the activity of CYP51A1 enzyme to synthesize essential ergosterol (Matin, Chakraborty, et al., 2020; Matin, Hasan, et al., 2020; . On the basis of our study, it could be thought that the effective uridine esters 13-15 restricted CYP51 enzyme activity by binding with it, and hence, fungal ergosterol could not form. The unavailability of ergosterol ultimately causes disruption of the fungal cell wall causing organisms death.

On the other hand, saturated fatty acyl esters with 5, 7-8 smaller chain length (C2-C6) showed slightly better binding score (-7.9 to -8.0 kcal/mol) than chlorobenzoyl esters 13-15 (-7.6 to -7.8 kcal/mol) with 7BQY (COVID-19 main protease). Methanesuphonyl ester 6 (-6.6 kcal/mol) and HQC (-6.1 kcal/ mol) had lower binding energy. The present results, thus, might have pave the way for discovery of novel drugs related to COVID-19 treatment although related further study is needed.

Covid-19, toxicity and evolution of MDR organisms in the microbial world caused a huge threat, and hence, novel antiviral and antimicrobial drugs preferably with newer mode of action, and reduced side-effects are essential. In this context several uridine esters were analyzed in vitro and in silico for their antimicrobial, thermodynamic, molecular docking and drug likeness properties. The study revealed that chloro-substituted benzoyl esters 13-15 of uridine were more potential against fungal organisms with better pharmacokinetic and drug likeness spectra. In addition, some aliphatic acyl esters especially with lower carbon chains (C2-C6) showed promising binding score with COVID-19 main protease (7BQY). In this promising context, more drug likeness in vitro and in vivo study such as nontoxic concentration towards healthy cells may be conducted in the near future.

New type mesoporous conjugate material for selective optical copper(II) ions monitoring & removal from polluted waters

Uridine supplementation enhances hepatic mitochondrial function in thymidine-analogue treated HIV-infected patients

Antidepressant-like effects of uridine and omega-3 fatty acids are potentiated by combined treatment in rats

Simple and rapid synthesis of some nucleoside derivatives: Structural and spectral characterization

Organic chemistry

Estimation of heats of formation of organic compounds by additivity methods

A 40-year journey in search of selective antiviral chemotherapy

Uridine function in the central nervous system

Enzymatic synthesis of nucleoside analogues from uridines and vinyl esters in a continuous-flow microreactor

Synthesis and cytotoxic activity of tri-acyl ester derivatives of uridine in breast cancer cells

Evaluation of green tea polyphenols as novel corona virus (SARS CoV-2) main protease (Mpro) inhibitors -an in silico docking and molecular dynamics simulation study

Comparative studies on selective acylation and antimicrobial activities of some D-glucofuranose derivatives

Synthesis and antimicrobial activities of some acylated uridine derivatives

Selective acylation of uridine using the dibutyltin oxide and direct methods

Synthesis and antibacterial activities of some uridine derivatives

Antibacterial activities of some selectively acylated uridine derivatives

Antifungal activities of some methyl 6-O-trityl-a-D-mannopyranosides

Comparative studies on selective acylation of uridine using the dibutyltin oxide and direct methods

Selective myristoylation and palmitoylation of uridine

Why the standard B3LYP/6-31G Ã model chemistry should not be used in DFT calculations of molecular thermochemistry: Understanding and correcting the problem

Ethnomedicines of Indian origin for combating COVID-19 infection by hampering the viral replication: Using structure-based drug discovery approach

Current trends and future approaches in smallmolecule therapeutics for COVID-19

Use of dipole moment as a parameter in drug-receptor interaction and quantitative structureactivity relationship studies

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings

Crystal structure of aspirin-acetylated human cyclooxygenase-2: Insight into the formation of products with reversed stereochemistry

Synthesis of D-glucose derived oxetane: 1,2-O-isopropylidene-4-(S)-3-O,4-C-methylene-5-O-methanesulfonyl-b-L-threopento-1,4-furanose

Synthesis and antimicrobial study of some methyl 4-O-palmitoyl-a-L-rhamnopyranoside derivatives

Synthesis, PASS predication, in vitro antimicrobial evaluation and pharmacokinetic study of novel n-octyl glucopyranoside esters

Design and synthesis of benzyl 4-O-lauroyl-a-L-rhamnopyranoside derivatives as antimicrobial agents

Synthesis and comparative antimicrobial studies of some acylated Dglucofuranose and D-glucopyranose derivatives

Synthesis and comparative antibacterial studies of some benzylidene monosaccharide benzoates

Synthesis, characterization, ADMET, PASS predication, and antimicrobial study of 6-O-lauroyl mannopyranosides

Novel mannopyranoside esters as sterol 14a-demethylase inhibitors: Synthesis, PASS predication, molecular docking, and pharmacokinetic studies

Synthesis, spectroscopic characterization, molecular docking, and ADMET studies of mannopyranoside esters as antimicrobial agents

Synthesis of some methyl 4-O-octnoyl-a-L-rhamnopyranoside derivatives

PASS-predication and in vitro antimicrobial activity of benzyl 4-O-benzoyl-a-L-rhamnopyranoside derivatives

PASS Predication, antiviral, in vitro antimicrobial, and ADMET studies of rhamnopyranoside esters

Crown-gall tumours in potato discs and brine shrimp lethality: Two simple bioassays for higher plant screening and fractionation

Absolute electronegativity and hardness correlated with molecular orbital theory

pkCSM: Predicting small-molecule pharmacokinetic and toxicity properties using graphbased signatures

Pollutants inducing epigenetic changes and diseases

Discovery and synthesis of a phosphoramidate prodrug of a Pyrrolo[2,1-f][triazin-4-amino] adenine C-nucleoside (GS-5734) for the treatment of ebola and emerging viruses

5-FdUrd-araC heterodinucleoside re-establishes sensitivity in 5-FdUrd-and AraC-resistant MCF-7 breast cancer cells overexpressing ErbB2. Differentiation

Pharmaceutically active uridine esters

Stilbene-based natural compounds as promising drug candidates against COVID-19

Biochemistry of uridine in plasma