key: cord-1039891-4caax43q authors: Alesawy, Mohamed S.; Elkaeed, Eslam B.; Alsfouk, Aisha A.; Metwaly, Ahmed M.; Eissa, Ibrahim H. title: In Silico Screening of Semi-Synthesized Compounds as Potential Inhibitors for SARS-CoV-2 Papain-like Protease: Pharmacophoric Features, Molecular Docking, ADMET, Toxicity and DFT Studies date: 2021-10-30 journal: Molecules DOI: 10.3390/molecules26216593 sha: 13148b2327a862461c262acf4de304fc03b03328 doc_id: 1039891 cord_uid: 4caax43q Papain-like protease is an essential enzyme in the proteolytic processing required for the replication of SARS-CoV-2. Accordingly, such an enzyme is an important target for the development of anti-SARS-CoV-2 agents which may reduce the mortality associated with outbreaks of SARS-CoV-2. A set of 69 semi-synthesized molecules that exhibited the structural features of SARS-CoV-2 papain-like protease inhibitors (PLPI) were docked against the coronavirus papain-like protease (PLpro) enzyme (PDB ID: (4OW0). Docking studies showed that derivatives 34 and 58 were better than the co-crystallized ligand while derivatives 17, 28, 31, 40, 41, 43, 47, 54, and 65 exhibited good binding modes and binding free energies. The pharmacokinetic profiling study was conducted according to the four principles of the Lipinski rules and excluded derivative 31. Furthermore, ADMET and toxicity studies showed that derivatives 28, 34, and 47 have the potential to be drugs and have been demonstrated as safe when assessed via seven toxicity models. Finally, comparing the molecular orbital energies and the molecular electrostatic potential maps of 28, 34, and 47 against the co-crystallized ligand in a DFT study indicated that 28 is the most promising candidate to interact with the target receptor (PLpro). As of 10 September 2021, the uncontrolled SARS-CoV-2 had infected 222,406,582 and killed 4,592,934 people all over the world, according to the WHO [1]. The alarming spread necessitates contentious work until a cure is discovered. Computational (in silico, computer-aided, or cheminformatics) approaches have been used in various fields related to drug discovery [2] , such as molecular docking [3] [4] [5] , pharmacophore studies [6] , drug molecular design [7, 8] , QSAR [9] , toxicity prediction [10] [11] [12] , ADMET assessment [13] [14] [15] , and DFT calculations [16] . These approaches have been practiced fruitfully and frequently in several scientific studies to find a treatment against COVID-19 with the advantage of taking less effort, time, and cost [17] [18] [19] [20] [21] . MOE software was used to conduct docking studies (Supplementary Materials) on the investigated derivatives, with co-crystallized ligand S88 as a reference. The study aimed at getting a deeper insight into the binding modes of the examined semi-synthesized molecules in the active site of the coronavirus papain-like protease (PLpro) enzyme (PDB ID: (4OW0)). The docking method was validated through redocking of the co-crystallized ligand in the enzyme active site. The protocol's applicability was confirmed through the demonstration of small RMSD (0.54 Å) between the co-crystallized pose and the re-docked one ( Figure 4 ). In this study, we relied on the corrected mode of binding of the examined semisynthesized molecules and S88 as well as the values of the binding free energy (∆G) between them. Table 1 illustrates the calculated ∆G of the tested semi-synthesized molecules and (S88) against the coronavirus papain-like protease enzyme. The semi-synthesized molecules 34 and 58 showed affinity values of −8.97 and −8.65, respectively, that were higher than that of the redocked ligand S88 (−8.59 kcal/mol). Moreover, the semisynthesized molecules 17, 28, 31, 40, 41, 43, 47, 54 and 65 revealed binding energy scores ranging from −8.33 to −8.57 kcal/mol, which were highly close to the redocked ligand S88. On the other hand, the other semi-synthesized molecules demonstrated affinity values lower than S88. The proposed binding mode of the redocked ligand S88 revealed an affinity value of −8.59 kcal/mol. This high binding affinity is probably attributed to the formation of two hydrogen-bonding interactions. One was formed between the N-H group of the amide moiety and Tyr269 while the other was formed between the nitrogen atom of the pyridine ring and Asp165. Additionally, the naphthyl moiety formed four hydrophobic interactions with Tyr269 and Pro249. These results were found to be consistent with the reported data [45] ( Figure 5 ). The docking simulation of compound 28 revealed that it has a good fitting into the enzyme active site with a docking score of −8.48 kcal/mol. The oxygen of the carbonyl group formed one hydrogen bond with the essential amino acid Tyr269. Additionally, the NH group of the pyrrole ring formed one hydrogen bond with Ala247. The 2,3,4,5tetrahydro-1H-pyrido [4,3-b] indole moiety formed Tyr265, Asp165, and Pro248 ( Figure 6 ). As illustrated in Figure 7 , compound 34 possessed a significant potential binding affinity (∆G = −8.97 kcal/mol) into the papain-like protease active site. This high binding affinity, which is higher than ligand S88, presumably attributed to the formation of one hydrogen bond interaction with Arg167. In addition, the 1H-indole formed two hydrophobic interactions with Lys158 and Leu163. Investigation of the top docking poses of compounds 47 and 54 (affinity values of −8.57 and 8.33 kcal/mol) respectively, demonstrated that compound 47 formed one hydrogen bond interaction with Tyr269. In addition, it formed two hydrophobic interactions with Pro249 and Lys158. Compound 54 formed two hydrogen bonds with the essential amino acid Tyr269 and Tyr265. In addition, it formed two hydrophobic interactions with Lys158 and Arg167 in the active site of PLpro. Figures 8 and 9 . The binding mode of compound 58 (affinity value of −8.65 kcal/mol) was better than ligand S88. In detail, the amide moiety formed one hydrogen bond with fundamental amino acid Asp165, and NH of the pyrrole ring formed another hydrogen bond with Glu168. Furthermore, it formed aromatic stacking interactions (4 pi-cation bonds) with Tyr269, Tyr265, and Arg167 ( Figure 10 ). An in silico computational evaluation of the physicochemical properties and profiling pharmacokinetics for the most active eleven semi-synthesized molecules, with ligand S88 as a reference compound, were conducted. The oral absorption of a drug is more likely to be better if the molecule fulfills at least three of the four principles of the Lipinski rules, listed below: (1) H bond donors (OH, NH, and SH) ≤ 5; (2) H bond acceptors (N, O, and S atoms) ≤ 10; (3) molecular weight < 500; (4) log P < 5. The bioavailability of compounds that violate more than one of these requirements is unlikely to be high. Moreover, reduced molecular flexibility, as measured by the number of rotatable bonds, and low polar surface area are found to be important predictors of good oral bioavailability [50] . Compounds with 10 or fewer rotatable bonds and a polar surface area of 140 Å or less have a high probability of good oral bioavailability [50, 51] . The results, given in Table 2 , revealed that all tested semi-synthesized molecules and reference ligand S88 showed no violation of Lipinski's rule except compound 31 (the Log P of compounds 31 was anticipated to be more than 5). S88 and favipiravir were used as reference drugs in ADMET studies for the most active eleven semi-synthesized molecules using Discovery studio 4.0 software. ADMET studies include many descriptors. The predicted descriptors are listed in Table 3 . All tested semi-synthesized molecules and favipiravir showed BBB penetration levels ranging from medium to low except compound 31, which displayed a very low BBB penetration level, and ligand S88 showed a high BBB penetration level. All semi-synthesized molecules, favipiravir, and ligand S88 have good absorption behavior except compound 31, which is expected to have a moderate absorption level. Moreover, the solubility level of the semi-synthesized molecules is projected to be better than or even comparable to that of the S88, which showed a low solubility level, except compound 31 that showed a very low solubility level. On the other hand, favipiravir demonstrated an optimal solubility level. All examined semi-synthesized molecules and favipiravir were predicted to be noninhibitors of CYP2D6 except compounds 31, 34, 47, and S88. Hepatotoxicity predictions found that all of the tested compounds and ligand S88 are predicted to be non-toxic except compounds 17, 31, 41, 43 and favipiravir, which have unfavorable hepatotoxic effects. All tested semi-synthesized molecules and S88 were expected to bind to plasma proteins more than 90% except compounds 28, 40, 43, 54, and favipiravir ( Figure 11 ). Discovery Studio 4.0 software was used to generate toxicity predictions for the most active eleven semi-synthesized molecules, which were based on validated and assembled models as follows: FDA rat carcinogenicity [52, 53] , carcinogenic potency TD 50 [54] , rat maximum tolerated dose (MTD) [55, 56] , rat oral LD 50 [57] , rat chronic LOAEL [58, 59] , ocular irritancy [60] and skin [19, 60, 61] . As shown in Table 4 , most of the examined semi-synthesized molecules have low toxicity. All the tested semi-synthesized molecules are non-carcinogens except 54, 58, and S88, which were predicted to be carcinogens. All tested semi-synthesized molecules showed TD 50 values ranging from 0.31 to 1.86 mg/kg body weight/day, which were lower than S88 (2.54 mg/kg body weight/day), except compounds 17, 31, 40, and 47 that showed TD 50 values of 5.57, 5.32, 10.31 and 14.54 mg/kg body weight/day, respectively, which were higher than S88. All the investigated semi-synthesized molecules revealed a maximum tolerated dose with a range of 0.045 to 0.122 g/kg body weight that was lower than S88 (0.124 g/kg body weight) except compounds 16, 28, 31, and 34, which demonstrated maximum tolerated doses of 0.142, 0.144, 0.369 and 0.328 g/kg body weight, respectively, which were higher than S88. All tested semi-synthesized molecules showed oral LD 50 values ranging from 2.97 to 32.39 mg/kg body weight/day that were higher than the LD 50 value of S88 (1.229 mg/kg body weight/day), except compound 31, which revealed an oral LD 50 value of 0.251 mg/kg body weight/day, which was lower than S88. Semi-synthesized molecules 54, 58, and 65 showed LOAEL values of 0.015, 0.001, and 0.018 g/kg body weight, respectively. These values were lower than S88 (0.035 g/kg body weight), while other semi-synthesized molecules revealed rat chronic LOAEL values ranging from 0.039 to 0.539 g/kg body weight, which were higher than S88. Moreover, all the tested semi-synthesized molecules and S88 were predicted to be mild irritants against the ocular irritancy model. On the other hand, the examined semi-synthesized molecules and S88 were expected to be non-irritant against the skin irritancy model except compound 58, which was anticipated to be a mild skin irritant. Figure 11 . The expected ADMET study of the most potent semi-synthesized molecules. DFT parameters, including total energy [62] , HOMO [63] , LUMO [63] , gap energy [64] , and dipole moment [65, 66] , were studied for the most semi-synthesized molecules 28, 34, and 47 using Discovery studio 4.0 software. The co-crystallized S88 was used as a reference molecule. The total energies of 28, 34, 47, and S88 were −1422.912, −1285.184, −1252.334, and 1242.947kcal/mol, respectively. These results indicated that compound 28 has higher total energy than 34 and 47 and is expected to have a more efficient interaction with PLpro. Accordingly, compound 28 can bind more efficiently with PLpro than 34 and 47. Furthermore, the dipole moment values of 28, 34, 47, and S88 were 2.790, 1.558, 2.249, and 3.542, respectively (Table 5) . These values indicate that 28 has a higher dipole moment than compounds 34 and 47. Based on these findings, it was expected that compound 28 can easily form hydrogen bonds and non-bonded interactions with PLpro, which, consequently, leads to an increased binding affinity with the target receptor during SARS-CoV-2 inhibition. Therefore, compound 28 is considered the most promising candidate to interact with the target receptor. As reported, HOMO and LUMO have a key role in chemical stability and reactivity [67] . Compound 28 had a gap energy value of 0.134 kcal/mol, which is higher than that of compounds 34 (0.099 kcal/mol) and 47 (0.097kcal/mol). The increased gap energy of compound 28 indicates the higher stability of this compound. Figure 12 showed the spatial distribution of molecular orbitals for the tested compounds. MEP demonstrates the total electrostatic potential of a molecule in three dimensions depending on its partial charges, electronegativity, and chemical reactivity [68] . Identifying the electrostatic potential will help in the understanding of the drug's binding mode against a PLpro [69] . MEP displays the electronegative atoms (negative values) in red. Electronegative atoms act as hydrogen bonding acceptors. On the other hand, it displays electron-poor atoms (positive value) in blue. Electron-poor atoms act as hydrogen bonding donors. It displays the neutral atoms (zero values) in a green to yellow color. Neutral atoms can form πand other types of hydrophobic interactions. Such information facilitates the prediction of the chemical reaction and the binding mode with the biological target [70] . Compound 28 showed five red patches and two blue patches, which can form hydrogen bond acceptors and hydrogen bond donors, respectively. The aromatic moieties showed yellow patches, which can form hydrophobic interactions with hydrophobic amino acid residues (Figures 12 and 13) . Compounds 34 and 47 showed four red patches, which can form hydrogen bond acceptors. Compound 34 showed three red patches and two blue patches. The aromatic moieties of these compounds showed yellow patches which can form hydrophobic interactions with hydrophobic amino acid residues (Figures 12 and 13) . As compound 28 showed five red patches, this explains its high binding energy (−8.48 kcal/mol) and ability to form two hydrogen bonds in the docking procedure. The yellow patches on aromatic moieties facilitated the hydrophobic interaction with the target receptor. Compounds 34 and 47 showed four red patches in MEP, which clarified the formation of two and three hydrogen bonds, respectively. In addition, these compounds showed high binding energies of −8.97 and −8.57 kcal/mol, respectively. MEP demonstrates the total electrostatic potential of a molecule in three dimensions depending on its partial charges, electronegativity, and chemical reactivity [68] . Identifying the electrostatic potential will help in the understanding of the drug's binding mode against a PLpro [69] . MEP displays the electronegative atoms (negative values) in red. Electronegative atoms act as hydrogen bonding acceptors. On the other hand, it displays electron-poor atoms (positive value) in blue. Electron-poor atoms act as hydrogen bonding donors. It displays the neutral atoms (zero values) in a green to yellow color. Neutral atoms can form π-and other types of hydrophobic interactions. Such information facilitates the prediction of the chemical reaction and the binding mode with the biological target [70] . Compound 28 showed five red patches and two blue patches, which can form hydrogen bond acceptors and hydrogen bond donors, respectively. The aromatic moieties showed yellow patches, which can form hydrophobic interactions with hydrophobic As compound 28 showed five red patches, this explains its high binding energy (−8.48 kcal/mol) and ability to form two hydrogen bonds in the docking procedure. The yellow patches on aromatic moieties facilitated the hydrophobic interaction with the target receptor. Compounds 34 and 47 showed four red patches in MEP, which clarified the formation of two and three hydrogen bonds, respectively. In addition, these compounds showed high binding energies of −8.97 and −8.57 kcal/mol, respectively. A set of 69 semi-synthesized molecules that exhibited the structural features of PLpro inhibitors (PLPI) were screened in silico to select the most potent inhibitor of PLpro enzyme (PDB ID: (4OW0). Docking studies showed that 11 molecules exhibited good binding modes and binding free energies. The pharmacokinetic profiling study excluded an unsuitable compound. Furthermore, ADMET and toxicity studies favored three molecules. Finally, a DFT study has been carried out and indicated that N-(3,4-dimethoxyphenethyl)-4-oxo-4-(1,3,4,5-tetrahydro-2H-pyrido [4,3-b] indol-2-yl)butanamide (28) is the most promising PLpro inhibitor. Further work has to be conducted to build on the presented results in the hopes of finding a cure. A set of 69 semi-synthesized molecules that exhibited the structural features of PLpro inhibitors (PLPI) were screened in silico to select the most potent inhibitor of PLpro enzyme (PDB ID: (4OW0). Docking studies showed that 11 molecules exhibited good binding modes and binding free energies. The pharmacokinetic profiling study excluded an unsuitable compound. Furthermore, ADMET and toxicity studies favored three molecules. Finally, a DFT study has been carried out and indicated that N-(3,4-dimethoxyphenethyl)-4-oxo-4-(1,3,4,5-tetrahydro-2H-pyrido [4,3-b] indol-2-yl)butanamide (28) is the most promising PLpro inhibitor. Further work has to be conducted to build on the presented results in the hopes of finding a cure. Crystal structure of human coronavirus papain-like protease inhibitor [PDB ID: 4OW0, resolution: 2.10 Å] was obtained from Protein Data Bank. The docking investigation was accomplished using MOE2014 software. At first, the crystal structure of SARS-CoV-2 helicase was prepared by removing water molecules. Only one chain was retained beside the co-crystallized ligand (S88). Then, the selected chain was protonated and subjected to minimization of energy process. Next, the active site of the target protein was defined. Structures of the tested compounds and the co-crystallized ligand were drawn using ChemBioDraw Ultra 14.0 and saved as MDL-SD format. Such file was opened using MOE to display the 3D structures which were protonated and subjected to energy minimization. Formerly, validation of the docking process was performed by docking the co-crystallized ligand against the isolated pocket of active site. The produced RMSD value indicated the validity of process. Finally, docking of the tested compounds was done through the dock option inserted in computer window. For each docked molecule, 30 docked poses were produced using ASE for scoring function and force field for refinement. The retain was kept at 30. The crystal parameters were adjusted at default values (Coordinates: Normal, Lattice Style: the same, Lattice: (1)). The results of the docking process were then visualized using Discovery Studio 4.0 software [71] . Pharmacokinetic profile of the compounds was determined using Discovery studio 4.0. [72] . ADMET descriptors (absorption, distribution, metabolism, excretion and toxicity) of the compounds were determined using Discovery studio 4.0. At first, the CHARMM force field was applied, then the tested compounds were prepared and minimized according to the preparation of small molecule protocol. Then ADMET descriptors protocol was applied to carry out these studies [73] . The toxicity parameters of the tested compounds were calculated using Discovery studio 4.0. Indinavir was used as a reference drug. At first, the CHARMM force field was applied then the compounds were prepared and minimized according to the preparation of small molecule protocol. Then different parameters were calculated from the toxicity prediction (extensible) protocol [73] [74] [75] . The DFT parameters (total energy, binding energy, HOMO, LUMO, gap energy, dipole moment, and electrostatic potential) were calculated using Discovery studio software. The tested compounds were prepared using prepare ligand protocol. Then, the prepared compounds were subjected to DFT calculation protocol using the default option [76] . Supplementary Materials: The following are available online, Table S1 : Detailed toxicity report, in addition to the method (Docking studies, ADMET studies, Toxicity studies and DFT studies). WHO. WHO Coronavirus (COVID-19) Dashboard. Available online In silico pharmacology for drug discovery: Applications to targets and beyond Structure-based drug design: Computational advances Screening of Some Sulfonamide and Sulfonylurea Derivatives as Anti-Alzheimer's Agents Targeting BACE1 and PPARγ Synthesis, biological evaluation and 2D-QSAR study of halophenyl bis-hydrazones as antimicrobial and antitubercular agents Pharmacophore modeling, molecular docking, QSAR, and in silico ADMET studies of gallic acid derivatives for immunomodulatory activity Computational multitarget drug design 3-dihydroxybenzoate alleviates oligodendrocyte damage following chronic cerebral hypoperfusion through IGF-1 receptor 2D and 3D-QSAR, molecular docking and ADMET properties in silico studies of azaaurones as antimalarial agents Discovery of thieno [2, 3-d] pyrimidine-based derivatives as potent VEGFR-2 kinase inhibitors and anti-cancer agents Design, synthesis, docking, ADMET studies, and anticancer evaluation of new 3-methylquinoxaline derivatives as VEGFR-2 inhibitors and apoptosis inducers Discovery of new 3-methylquinoxalines as potential anti-cancer agents and apoptosis inducers targeting VEGFR-2: Design, synthesis, and in silico studies c] quinoxaline derivatives as anticancer agents and apoptosis inducers New quinoxaline derivatives as VEGFR-2 inhibitors with anticancer and apoptotic activity: Design, molecular modeling, and synthesis 3-a: 4 , 3 -c] quinazoline derived DNA intercalators: Design, synthesis, in silico ADMET profile, molecular docking and anti-proliferative evaluation studies Discovery of new anticancer thiourea-azetidine hybrids: Design, synthesis, in vitro antiproliferative, SAR, in silico molecular docking against VEGFR-2, ADMET, toxicity, and DFT studies Isolation and In Silico Anti-COVID-19 Main Protease (Mpro) Activities of Flavonoids and a Sesquiterpene Lactone from Artemisia sublessingiana Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of ACEIs Against SARS-CoV-2 Targeting the hACE2 Receptor In Silico Studies of Some Isoflavonoids as Potential Candidates against COVID-19 Targeting Human ACE2 (hACE2) and Viral Main Protease (Mpro) Comprehensive virtual screening of the antiviral potentialities of marine polycyclic guanidine alkaloids against SARS-CoV-2 (COVID-19) Comprehensive In Silico Screening of the Antiviral Potentialities of a New Humulene Glucoside from Asteriscus hierochunticus against SARS-CoV-2 Traditional ancient Egyptian medicine: A review The Chinese herbal formulae (Yitangkang) exerts an antidiabetic effect through the regulation of substance metabolism and energy metabolism in type 2 diabetic rats Black ginseng and its saponins: Preparation, phytochemistry and pharmacological effects Chemical constituents of stems and leaves of Tagetespatula L. and its fingerprint Comparative biological evaluation of four endophytic fungi isolated from nigella sativa seeds New α-Pyrone derivatives from the endophytic fungus Embellisia sp Two new antileishmanial diketopiperazine alkaloids from the endophytic fungus Trichosporum sp Biological evaluation and molecular docking study of metabolites from Salvadora Persica L. Growing in Egypt Chemical Constituents of Tagetes patula and Their Neuroprotecting Action Induction of apoptosis in human cancer cells through extrinsic and intrinsic pathways by Balanites aegyptiaca furostanol saponins and saponin-coated silvernanoparticles New combination approaches to combat methicillin-resistant Staphylococcus aureus (MRSA) Isolation, Crystal Structure, and In Silico Aromatase Inhibition Activity of Ergosta-5, 22-dien-3β-ol from the Fungus Gyromitra esculenta Nigrosphaerin A a new isochromene derivative from the endophytic fungus Nigrospora sphaerica Antileukemic α-pyrone derivatives from the endophytic fungus Alternaria phragmospora In Vitro and In Silico Cytotoxic and Antibacterial Activities of a Diterpene from Cousinia alata Schrenk Antileishmanial Derivatives of Humulene from Asteriscus hierochunticus with in silico Synthesis and Molecular Docking of Some Grossgemin Amino Derivatives as Tubulin Inhibitors Targeting Colchicine Binding Site Natural products as sources of new drugs from 1981 to 2014 Semisynthetic and newly designed derivatives based on natural chemical scaffolds: Moving beyond natural products to fight Trypanosoma cruzi ex vivo evaluation, and SAR studies of coumarin derivatives as potential antiasthmatic drugs Papain-like protease regulates SARS-CoV-2 viral spread and innate immunity The SARS-coronavirus papain-like protease: Structure, function and inhibition by designed antiviral compounds A noncovalent class of papain-like protease/deubiquitinase inhibitors blocks SARS virus replication X-ray structural and biological evaluation of a series of potent and highly selective inhibitors of human coronavirus papain-like proteases Structure of papain-like protease from SARS-CoV-2 and its complexes with non-covalent inhibitors Design and synthesis of new 4-(2-nitrophenoxy) benzamide derivatives as potential antiviral agents: Molecular modeling and in vitro antiviral screening Activity profiling and crystal structures of inhibitor-bound SARS-CoV-2 papain-like protease: A framework for anti-COVID-19 drug design Natural-Product-Based Library Molecular properties that influence the oral bioavailability of drug candidates Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings Classification of kinase inhibitors using a Bayesian model Structure-Activity Relationships for Carcinogenic Potential Comparison of the dependence of the TD50 on maximum tolerated dose for mutagens and nonmutagens Correlation between carcinogenic potency and the maximum tolerated dose: Implications for risk assessment Comparison of in silico tools for evaluating rat oral acute toxicity In silico models for repeated-dose toxicity (RDT): Prediction of the no observed adverse effect level (NOAEL) and lowest observed adverse effect level (LOAEL) for drugs Assessment of the oral rat chronic lowest observed adverse effect level model in TOPKAT, a QSAR software package for toxicity prediction The Draize eye test Review on skin and mucous-membrane irritation tests and their application FT-IR spectra and total energy distribution of 3-pentyl-2, 6-diphenylpiperidin-4-one: DFT method Visible light driven photo-degradation of Congo red by TiO2ZnO/Ag: DFT approach on synergetic effect on band gap energy Synthesis, Characterization, DFT calculation and Antimicrobial Activity of Co (II) and Cu (II) complexes with azo dye Frontier Orbitals and Organic Chemical Reactions Ab initio HF, DFT and experimental (FT-IR) investigation of vibrational spectroscopy of PN, N-dimethylaminobenzylidenemalononitrile (DBM) Correlation found between the HOMO-LUMO energy separation and the chemical reactivity at the most reactive site for isolated-pentagon isomers of fullerenes Vibrational and electronic investigations, thermodynamic parameters, HOMO and LUMO analysis on Lornoxicam by density functional theory Electrostatic Potential Energy in Protein-Drug Complexes Synthesis, spectroscopic characterization, molecular docking, and ADMET studies of mannopyranoside esters as antimicrobial agents Structural and Biological Evaluation of a Series of Potent and Highly Selective Inhibitors of Human Coronavirus Papain-like Proteases Synthesis, docking, QSAR, ADMET and antimicrobial evaluation of new quinoline-3-carbonitrile derivatives as potential DNA-gyrase inhibitors New quinoxaline-2 (1H)-ones as potential VEGFR-2 inhibitors: Design, synthesis, molecular docking, ADMET profile and anti-proliferative evaluations Anticancer activity, spectroscopic and molecular docking of some new synthesized sugar hydrazones, Arylidene and α-Aminophosphonate derivatives Design and discovery of new 1, 2, 4-triazolo [4, 3-c] quinazolines as potential DNA intercalators and topoisomerase II inhibitors Silico Exploration of Potential Natural Inhibitors against SARS-CoV-2 nsp10 Institutional Review Board Statement: Not applicable. Data Availability Statement: Data is contained within the article. The authors declare no conflict of interest.Sample Availability: Samples of the compounds are not available from the authors.