key: cord-0807493-9nqxcj0q
authors: Elzupir, Amin Osman
title: Molecular Docking and Dynamics Investigations for Identifying Potential Inhibitors of the 3-Chymotrypsin-like Protease of SARS-CoV-2: Repurposing of Approved Pyrimidonic Pharmaceuticals for COVID-19 Treatment
date: 2021-12-09
journal: Molecules
DOI: 10.3390/molecules26247458
sha: 90cf9f0050f56c506dbb07cb3bed993498e3df63
doc_id: 807493
cord_uid: 9nqxcj0q

This study demonstrates the inhibitory effect of 42 pyrimidonic pharmaceuticals (PPs) on the 3-chymotrypsin-like protease of SARS-CoV-2 (3CL(pro)) through molecular docking, molecular dynamics simulations, and free binding energies by means of molecular mechanics–Poisson Boltzmann surface area (MM-PBSA) and molecular mechanics–generalized Born surface area (MM-GBSA). Of these tested PPs, 11 drugs approved by the US Food and Drug Administration showed an excellent binding affinity to the catalytic residues of 3CL(pro) of His41 and Cys145: uracil mustard, cytarabine, floxuridine, trifluridine, stavudine, lamivudine, zalcitabine, telbivudine, tipiracil, citicoline, and uridine triacetate. Their percentage of residues involved in binding at the active sites ranged from 56 to 100, and their binding affinities were in the range from −4.6 ± 0.14 to −7.0 ± 0.19 kcal/mol. The molecular dynamics as determined by a 200 ns simulation run of solvated docked complexes confirmed the stability of PP conformations that bound to the catalytic dyad and the active sites of 3CL(pro). The free energy of binding also demonstrates the stability of the PP–3CL(pro) complexes. Citicoline and uridine triacetate showed free binding energies of −25.53 and −7.07 kcal/mol, respectively. Therefore, I recommend that they be repurposed for the fight against COVID-19, following proper experimental and clinical validation.

Over a year has passed since the COVID-19 pandemic began. Some vaccines, such as those by Pfizer and Moderna, and some drugs, such as remdesivir, have been approved for use in therapy. The efforts by governments, health organizations, and other sectors to stem the alarmingly increasing numbers of deaths and cases were unprecedented [1] [2] [3] [4] [5] [6] . However, SARS-CoV-2 continues to threaten the world, with over four million deaths and 227 million cases as of 16 September 2021 (https://www.worldometers.info/coronavirus/ accessed on 27 November 2021). COVID-19 was declared a pandemic by the World Health Organization on 11 March 2020. Today, the new SARS-CoV-2 virus, the causative agent of COVID-19, has been detected in almost every country on the planet [5, [7] [8] [9] .

Coronaviruses are positive-stranded RNA viruses with the largest viral genomes ever known, ranging from 16 to 32 kb. The 3-chymotrypsin-like protease (3CL pro ) produced by SARS-CoV-2 is a cysteine protease encoded as nonstructural protein 3 in the polyprotein. 3CL pro is responsible for the cleavage of 11 specific sites of polyproteins (pp1a, pp1ab) produced by the 229E gene. These polyproteins are involved in the production of a functional polypeptide essential for viral replication and transcription. Further, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . ther, the specificity of 3CL pro is dissimilar to that of human host-cell protease. Thus, 3CL pro has become the focus of drug repurposing and development programs to combat the COVID-19 pandemic [10] [11] [12] [13] . Recent and ongoing research has reported that some pharmaceutical, synthetic, and natural products can act as 3CL pro inhibitors or against SARS-CoV-2 in general. These include selenium-containing heterocyclic compounds, chloroquine phosphate, indinavir, darunavir, lopinavir, eravacycline, naproxen, salix cortex, antioxidants, chiral phytochemicals from Opuntia ficus-indica, elbasvir, valrubicin, favipiravir isoflavone, and myricitrin [6, [14] [15] [16] [17] [18] [19] [20] [21] [22] . Although some of these have entered human clinical trials or were even approved, more studies are still needed. The importance of the pyridone ring was highlighted in synthetic materials and drugs containing pyridone [11, 23] . The pyrimidone ring has the exact shape of pyridine but is more functionalized and electrondeficient. Herein, we screened the inhibitory activity of 42 approved pyrimidonic pharmaceuticals (PPs) against 3CL pro using a combination of molecular docking analyses, molecular dynamics simulations, and calculations of the MM-PBSA and MM-GBSA binding free energies. The sites of action of active inhibitors were investigated, discussed, and explored.

The PPs were selected using the search engine of the drug bank database. The search uncovered 46 PPs; the pharmaceuticals containing caffeine were entirely excluded as all except enprofylline have previously been studied. In addition, the macropolymeric drug mipomersen, drugs composed of a mixture of drugs, and withdrawn drugs were not included in this study. The chosen drugs were classified into four categories according to their structures. 1PPs have only one heterocycle, 2aPPs have two, 2bPPs have two heterocycles with a pyrimidone ring having an extra carbonyl group, and 3PPs have three or more heterocycles (Table 1) . 

The 3D structures of the selected PPs were downloaded from the PubChem website as SDF files; their energy was minimized for 10,000 steepest descent steps at 5000 conjugate gradient steps using antechamber plugin UCSF Chimera [24, 25] . For alogliptin, the 3D structure was obtained by utilizing OpenBabel converter tools and ChemSkech [26] . The crystal structure of SARS-CoV-2 3CL pro was obtained from the Protein Data Bank database website (PDB ID: 6Y2E). For analysis, water was removed from the 3CL pro structure, and the energy was then minimized for 1000 steepest descent steps at 20 conjugate gradient steps.

Blind molecular docking experiments were performed using the AutoDock Vina tool implemented with the interactive visualization and analysis program UCSF Chimera. The default parameter values were adopted with a grid box (−15 × −25 × 15) Å, centered at (35, 65 , 65) Å. The predicted affinity values of the score were observed using the View Dock tool. The binding between ligands and 3CL pro active sites and the images were processed and visualized using UCSF Chimera [24] [25] [26] [27] [28] .

MD simulations were performed as previously described [29] . The PP ligands were separated from the docked complexes using UCSF Chimera. The missed hydrogens were added and saved as PDB files using AMBER's large-structure serial numbering. Topology files and parameters of the receptor and the ligands were made using leap and antechamber of Amber Tools 21 [30, 30] , utilizing Amber force fields of GAFF2 [31] and ff14SB [32] to assign inhibitors and 3CL pro structure, respectively. The systems were solvated with TIP3P water molecules [33] and were neutralized via sodium ions. Subsequently, molecular dynamics (MD) simulations were performed by means of the Nanoscale Molecular Dynamics (NAMD) Simulation 2.6 program [34] . Each system was 

The 3D structures of the selected PPs were downloaded from the PubChem website as SDF files; their energy was minimized for 10,000 steepest descent steps at 5000 conjugate gradient steps using antechamber plugin UCSF Chimera [24, 25] . For alogliptin, the 3D structure was obtained by utilizing OpenBabel converter tools and ChemSkech [26] . The crystal structure of SARS-CoV-2 3CL pro was obtained from the Protein Data Bank database website (PDB ID: 6Y2E). For analysis, water was removed from the 3CL pro structure, and the energy was then minimized for 1000 steepest descent steps at 20 conjugate gradient steps.

Blind molecular docking experiments were performed using the AutoDock Vina tool implemented with the interactive visualization and analysis program UCSF Chimera. The default parameter values were adopted with a grid box (−15 × −25 × 15) Å, centered at (35, 65 , 65) Å. The predicted affinity values of the score were observed using the View Dock tool. The binding between ligands and 3CL pro active sites and the images were processed and visualized using UCSF Chimera [24] [25] [26] [27] [28] .

MD simulations were performed as previously described [29] . The PP ligands were separated from the docked complexes using UCSF Chimera. The missed hydrogens were added and saved as PDB files using AMBER's large-structure serial numbering. Topology files and parameters of the receptor and the ligands were made using leap and antechamber of Amber Tools 21 [30, 30] , utilizing Amber force fields of GAFF2 [31] and ff14SB [32] to assign inhibitors and 3CL pro structure, respectively. The systems were solvated with TIP3P water molecules [33] and were neutralized via sodium ions. Subsequently, molecular dynamics (MD) simulations were performed by means of the Nanoscale Molecular Dynamics (NAMD) Simulation 2.6 program [34] . Each system was 

The 3D structures of the selected PPs were downloaded from the PubChem website as SDF files; their energy was minimized for 10,000 steepest descent steps at 5000 conjugate gradient steps using antechamber plugin UCSF Chimera [24, 25] . For alogliptin, the 3D structure was obtained by utilizing OpenBabel converter tools and ChemSkech [26] . The crystal structure of SARS-CoV-2 3CL pro was obtained from the Protein Data Bank database website (PDB ID: 6Y2E). For analysis, water was removed from the 3CL pro structure, and the energy was then minimized for 1000 steepest descent steps at 20 conjugate gradient steps.

Blind molecular docking experiments were performed using the AutoDock Vina tool implemented with the interactive visualization and analysis program UCSF Chimera. The default parameter values were adopted with a grid box (−15 × −25 × 15) Å, centered at (35, 65 , 65) Å. The predicted affinity values of the score were observed using the View Dock tool. The binding between ligands and 3CL pro active sites and the images were processed and visualized using UCSF Chimera [24] [25] [26] [27] [28] .

MD simulations were performed as previously described [29] . The PP ligands were separated from the docked complexes using UCSF Chimera. The missed hydrogens were added and saved as PDB files using AMBER's large-structure serial numbering. Topology files and parameters of the receptor and the ligands were made using leap and antechamber of Amber Tools 21 [30, 30] , utilizing Amber force fields of GAFF2 [31] and ff14SB [32] to assign inhibitors and 3CL pro structure, respectively. The systems were solvated with TIP3P water molecules [33] and were neutralized via sodium ions. Subsequently, molecular dynamics (MD) simulations were performed by means of the Nanoscale Molecular Dynamics (NAMD) Simulation 2.6 program [34] . Each system was 

The 3D structures of the selected PPs were downloaded from the PubChem website as SDF files; their energy was minimized for 10,000 steepest descent steps at 5000 conjugate gradient steps using antechamber plugin UCSF Chimera [24, 25] . For alogliptin, the 3D structure was obtained by utilizing OpenBabel converter tools and ChemSkech [26] . The crystal structure of SARS-CoV-2 3CL pro was obtained from the Protein Data Bank database website (PDB ID: 6Y2E). For analysis, water was removed from the 3CL pro structure, and the energy was then minimized for 1000 steepest descent steps at 20 conjugate gradient steps.

Blind molecular docking experiments were performed using the AutoDock Vina tool implemented with the interactive visualization and analysis program UCSF Chimera. The default parameter values were adopted with a grid box (−15 × −25 × 15) Å, centered at (35, 65, 65) Å. The predicted affinity values of the score were observed using the View Dock tool. The binding between ligands and 3CL pro active sites and the images were processed and visualized using UCSF Chimera [24] [25] [26] [27] [28] .

MD simulations were performed as previously described [29] . The PP ligands were separated from the docked complexes using UCSF Chimera. The missed hydrogens were added and saved as PDB files using AMBER's large-structure serial numbering. Topology files and parameters of the receptor and the ligands were made using leap and antechamber of Amber Tools 21 [30, 30] , utilizing Amber force fields of GAFF2 [31] and ff14SB [32] to assign inhibitors and 3CL pro structure, respectively. The systems were solvated with TIP3P water molecules [33] and were neutralized via sodium ions. Subsequently, molecular dynamics (MD) simulations were performed by means of the Nanoscale Molecular Dynamics (NAMD) Simulation 2.6 program [34] . Each system was -

The 3D structures of the selected PPs were downloaded from the PubChem website as SDF files; their energy was minimized for 10,000 steepest descent steps at 5000 conjugate gradient steps using antechamber plugin UCSF Chimera [24, 25] . For alogliptin, the 3D structure was obtained by utilizing OpenBabel converter tools and ChemSkech [26] . The crystal structure of SARS-CoV-2 3CL pro was obtained from the Protein Data Bank database website (PDB ID: 6Y2E). For analysis, water was removed from the 3CL pro structure, and the energy was then minimized for 1000 steepest descent steps at 20 conjugate gradient steps.

Blind molecular docking experiments were performed using the AutoDock Vina tool implemented with the interactive visualization and analysis program UCSF Chimera. The default parameter values were adopted with a grid box (−15 × −25 × 15) Å, centered at (35, 65, 65) Å. The predicted affinity values of the score were observed using the View Dock tool. The binding between ligands and 3CL pro active sites and the images were processed and visualized using UCSF Chimera [24] [25] [26] [27] [28] .

MD simulations were performed as previously described [29] . The PP ligands were separated from the docked complexes using UCSF Chimera. The missed hydrogens were added and saved as PDB files using AMBER's large-structure serial numbering. Topology files and parameters of the receptor and the ligands were made using leap and antechamber of Amber Tools 21 [30, 30] , utilizing Amber force fields of GAFF2 [31] and ff14SB [32] to assign inhibitors and 3CL pro structure, respectively. The systems were solvated with TIP3P water molecules [33] and were neutralized via sodium ions. Subsequently, molecular dynamics (MD) simulations were performed by means of the Nanoscale Molecular Dynamics (NAMD) Simulation 2.6 program [34] . Each system was minimized for 1 ps at 273.15 K using the NVE ensemble. The temperature was gradually increased to 310 K using the NVT ensemble in a protocol consisting of 1600 minimization steps. Then, each system was minimized for 10 ps at 310 K followed by 200 ns of MD simulation control using the NVT ensemble at 310 K and a time step of 2 fs. In order to calculate electrostatic interactions, the particle mesh Ewald process and periodic boundary conditions were applied [35, 36] . The root-mean-square fluctuation (RMSF) and the root-mean-square deviation (RMSD) for each system were obtained by analyzing the trajectory using the VMD 1.8 program [37] .

The binding free energies of the PP-3CL pro complexes were calculated by means of molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) and molecular mechanics-generalized Born surface area (MM-GBSA) using the MMPBSA.py module of Amber Tools 21 [38] . The MD simulation over 200 ns provided several conformations sampled after equilibrium, using the last frames to lessen the computational cost. CPPTRAJ was used to obtain the snapshots [39] . The conformational changes were evaluated through quasi-harmonic entropy approximation [40] . The free energy of the binding interaction between inhibitors and 3CL pro complexes can be obtained via the following equations:

∆G gas = E vdw + E elec (3)

where ∆H represents enthalpy change, T∆S represents the entropic contribution, E vdw represents the van der Waals interaction energy, E ele represents the electrostatic interaction energy, ∆G sol represents the polar solvation energy, and E np represents the nonpolar solvation energy.

The results of the molecular docking are tabulated in Tables S1-S4. Figure 1 shows the catalytic dyad and the active sites of 3CL pro . The crucial residues HIS 41, GLY 143, SER 144, and CYS 145 forming the S1 site are shown in black. Then, PHE 140, LUE 141, ASN 142, HIS 163, GLU 166 (magenta), and the N-terminal amino acid residues (blue) are involved in the formation of the S1 subsite of the substrate-binding pocket. The MET 49, TYR 54, HIS 164, ASP 187, and ARG R188 residues form the S2 site (green). MET 165, LEU 167, GLN 189, THR 190, and GLN 192 comprise the S4 site (cyan). The SER 284, ALA 285, and LEU 286 residues (yellow) are a result of genetic mutation leading to an increase in the SASR-CoV-2 3CL pro activity of 3.6 fold over that of the 3CL pro predecessor of SARS-CoV [12, 41] . the SASR-CoV-2 3CL pro activity of 3.6 fold over that of the 3CL pro predecessor of SARS-CoV [12, 42] . 

The docked complexes of the top 11 candidates are depicted in Figure 2 . Their binding affinities to the active sites of 3CL pro are shown in Table 2 . The 3PPs showed significant interactions with the residues LEU 286, SER 284, and ALA 285, and a relatively lower interaction ratio to the catalytic dyad, in contrast to the other groups. Of the 3PPs, alogliptin and flavin mononucleotide were found to have the highest binding percentage with the catalytic dyad and to form hydrogen bonds with the S1 and S'1 sites. These were followed by riboflavin and sofosbuvir with an advantage in binding to the LEU 286 residue (Table S1 ). Flavin adenine dinucleotide showed excellent binding affinity to LEU 286 but not with the catalytic dyad. Zidovudine and gemcitabine demonstrated similar activity to alogliptin and flavin mononucleotide (Tables S2 and S3) .

Among the 2bPPs, anti-hepatitis B infection telbivudine, anti-orotic aciduria uridine triacetate, and anticancer tipiracil were found to have the highest binding to 3CL pro active sites, followed by antimetabolite floxuridine, anti-herpesvirus trifluridine, and anti-HIV stavudine. Here, it is worth noting the importance of the molecular structure, as this set differed from the previous one by its increased ability to bind to the 3CL pro catalytic dyad. The 2aPPs showed similar activity to that of the 2bPPs. Anticancer cytarabine, anti- Color indicates the residues involved in the formation of the S1 site (shown in magenta), S1 site from the other promotor (blue), S2 site (green), S4 site (cyan), and S1' site (black), in addition to SER 284, ALA 285, and LEU 286 (yellow). (b) Only the catalytic dyad and GLU 166 residues.

The docked complexes of the top 11 candidates are depicted in Figure 2 . Their binding affinities to the active sites of 3CL pro are shown in Table 2 . The 3PPs showed significant interactions with the residues LEU 286, SER 284, and ALA 285, and a relatively lower interaction ratio to the catalytic dyad, in contrast to the other groups. Of the 3PPs, alogliptin and flavin mononucleotide were found to have the highest binding percentage with the catalytic dyad and to form hydrogen bonds with the S1 and S'1 sites. These were followed by riboflavin and sofosbuvir with an advantage in binding to the LEU 286 residue (Table S1 ). Flavin adenine dinucleotide showed excellent binding affinity to LEU 286 but not with the catalytic dyad. Zidovudine and gemcitabine demonstrated similar activity to alogliptin and flavin mononucleotide (Tables S2 and S3) .

Among the 2bPPs, anti-hepatitis B infection telbivudine, anti-orotic aciduria uridine triacetate, and anticancer tipiracil were found to have the highest binding to 3CL pro active sites, followed by antimetabolite floxuridine, anti-herpesvirus trifluridine, and anti-HIV stavudine. Here, it is worth noting the importance of the molecular structure, as this set differed from the previous one by its increased ability to bind to the 3CL pro catalytic dyad. The 2aPPs showed similar activity to that of the 2bPPs. Anticancer cytarabine, antiglaucoma citicoline, and anti-HIV drugs lamivudine and zalcitabine showed promising inhibitory activity (Table S3) . Finally, but very importantly, of the 1PPs, the chemotherapy drug uracil mustard showed binding to the catalytic dyad with all of its simulated conformations, followed by anti-cytomegalovirus cidofovir (Table S4) . glaucoma citicoline, and anti-HIV drugs lamivudine and zalcitabine showed promising inhibitory activity (Table S3) .

Finally, but very importantly, of the 1PPs, the chemotherapy drug uracil mustard showed binding to the catalytic dyad with all of its simulated conformations, followed by anti-cytomegalovirus cidofovir (Table S4 ). (a) 

MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. Table 3 . The binding interactions of the potential pyrimidonic pharmaceuticals/3-chymotrypsin-like protease 3CL pro complexes at different times throughout the production runs.

Name 100 ns 150 ns 200 ns

Molecules 2021, 26, x FOR PEER REVIEW 10 of 18 MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. MD was performed on the hole complexes of the top 11 PPs candidates. Based on the conformer score energy from docking, the complex with the conformer with the lowest value and interacting with the 3CL pro active site was selected. The RMSDs were computed along the trajectories using the initial structure as a reference. Figure 3 shows that the binding of PPs significantly affected the equilibration states of 3CL pro , as the majority of the tested systems reached their equilibrium at around 100 ns. The PP-3CLpro complexes revealed relatively lower average values for the RMSDs, between 0.41 and 0.52 Å, throughout the simulation, clarifying their good behavior in forming stable complexes. Moreover, the fluctuations in the 3CL pro backbone residues were analyzed by means of the RMSF (Figure 4) . The 3CL pro /PPs exhibited lower fluctuations, particularly at the active site. The fluctuations at the catalytic dyad and GLU 166 were minor, demonstrating the loss of flexibility at these regions upon binding to the PPs. Table 3 shows the superior stability of the PP-3CL pro complexes formed throughout the production runs; these results support the use of these PPs as 3CL pro inhibitors. 

The data on the binding free energies of the PP-3CL pro complexes are tabulated in Table 4 . In both MM-GBSA and MM-PBSA, van der Waals and electrostatic interactions 

The data on the binding free energies of the PP-3CL pro complexes are tabulated in Table 4 . In both MM-GBSA and MM-PBSA, van der Waals and electrostatic interactions acted as driving forces for the PP ligands to bind to 3CL pro , contrasting the solvation energies. The MM-GBSA and MM-PBSA results suggest that the PPs have an excellent ability to inhibit 3CL pro . Of these PPs, citicoline revealed the most promising inhibitory activity, followed by uridine triacetate (Table 4 ). Among the challenges of discovering 3CL pro inhibitors for COVID-19 treatment, these inhibitors must be highly bioavailable inside the cytosol [13] . The 2aPP and 2bPP structures contain a primidone heterocycle and ribose ring. These heterocycles increase their hydrophilicity and solubility in plasma. Thus, they can satisfy the requirement of bioavailability. Citicoline, the most promising inhibitor among the PPs investigated, has high hydrophilicity and good ADME properties [42, 43] . The 2aPPs and 2bPPs also have an intermediate structure size among the PP groups. This sheds light on the importance of the size and general structural features of PPs acting as 3CL pro inhibitors.

Recent reports suggest a general hypothesis that 3CL pro inhibitors comprise electrophilic sites such as Michael acceptors [12] . That the pyrimidone ring is highly electrondeficient clarifies and confirms this hypothesis. The pyrimidone ring plays an essential role in PPs' inhibitory activity and has a high tendency to form hydrogen bonds, particularly with the GLU 166 residue. This may preclude the formation of the S1 pocket. Contacts between the pyrimidone ring and HIS 41 were observed in floxuridine, stavudine, and telbivudine. In all these cases, HIS 41 interacts with the oxygen of the pyrimidone group. Further, the electrophilic carbon, nitrogen, and oxygen in the pyrimidone ring were attracted to bind with the sulfur of the CYS 145 residue.

Interestingly, most PPs investigated were previously studied against SARS-CoV-2. For example, flavin mononucleotide and flavin adenine dinucleotide have been suggested as good 3CL pro and RNA-dependent RNA polymerase inhibitors, respectively [22, 44] . Riboflavin and sofosbuvir were shown to be suitable inhibitors of the spike protein S1 domain/ACE2 and RNA-dependent RNA polymerase [45] [46] [47] . Alogliptin was also suggested as a 3CL pro inhibitor; however, the enzymatic assay demonstrated its inactivity against 3CL pro [48] . Compelling clues have been found regarding the use of zidovudine and gemcitabine against spike protein/human ACE2 and in the inhibition SARS-CoV-2 in cell culture [49] [50] [51] [52] [53] . Gemcitabine and cidofovir were reported to inhibit SARS-CoV and SARS-CoV-2 proteins with IC 50 values of 4.95 µM and 36 µM [54, 55] . Further, telbivudine, tipiracil, cytarabine, and citicoline were recommended as 3CL pro inhibitors [56] [57] [58] [59] [60] [61] . This transitory literature scanning confirms these pharmaceuticals' activity against 3CL pro of SARS-CoV-2, as demonstrated in the present study.

To conclude, the inhibitory effect on 3CL pro by PPs was investigated based on their ability to form hydrogen bonds and van der Waals interactions with the 3CL pro active side through molecular docking, MD simulations, and the calculation of binding free energy. The overall analysis revealed 11 candidates from the initial set of 42 investigated PPs are promising 3CL pro inhibitors. These include citicoline and uridine triacetate as the best choices, followed by telbivudine, trifluridine, lamivudine, cytarabine, stavudine, zalcitabine, tipiracil, floxuridine, and flavin mononucleotide. The interactions of PPs with the catalytic dyad and the active sides of 3CL pro of SARS-CoV-2 were comprehensively and thoroughly investigated. The pyrimidone ring was found to play an essential role in the PPs' inhibitory activity.

Supplementary Materials: The following are available online, Table S1 : The binding affinities of the pyrimidonic pharmaceuticals (group 3PPs) with 3-chymotrypsin-like protease (3CLpro), Table S2 : The binding affinities of the pyrimidonic pharmaceuticals (group 2bPPs) with 3-chymotrypsin-like protease (3CLpro), Table S3 : The binding affinities of the pyrimidonic pharmaceuticals (group 2aPPs) with 3-chymotrypsin-like protease (3CLpro), Table S4 : The binding affinity of the pyrimidone containing-pharmaceuticals (group 1PCPs) with 3-chymotrypsinlike protease (3CLpro). 

The UK has approved a COVID vaccine-here's what scientists now want to know

Covid-19: UK approves Moderna vaccine to be given as two doses 28 days apart

Covid-19: Moderna applies for US and EU approval as vaccine trial reports 94.1% efficacy

Remdesivir: First approval

FDA approval of remdesivir-A step in the right direction

Silico Screening of Natural Products Isolated from Mexican Herbal Medicines against COVID-19. Biomolecules 2021

World Health Organization

A pneumonia outbreak associated with a new coronavirus of probable bat origin

World Health Organization. COVID-19 Weekly Epidemiological Update

The SARS-CoV-2 main protease as drug target

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors

Protease targeted COVID-19 drug discovery: What we have learned from the past SARS-CoV inhibitors?

Challenges for targeting SARS-CoV-2 proteases as a therapeutic strategy for COVID-19

Predicting commercially available antiviral drugs that may act on the novel coronavirus (SARS-CoV-2) through a drug-target interaction deep learning model

Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants

Therapeutic intervention of COVID-19 by natural products: A population-specific survey directed approach

Fast Identification of Possible Drug Treatment of Coronavirus Disease-19 (COVID-19) through Computational Drug Repurposing Study

Potential therapeutic agents for COVID-19 based on the analysis of protease and RNA polymerase docking

Antiviral Properties of the NSAID Drug Naproxen Targeting the Nucleoprotein of SARS-CoV-2 Coronavirus

Comparative Anti-Inflammatory Effects of Salix Cortex Extracts and Acetylsalicylic Acid in SARS-CoV-2 Peptide and LPS-Activated Human In Vitro Systems

Effectiveness of Natural Antioxidants against SARS-CoV-2?

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

Inhibition of SARS-CoV-2 main protease 3CLpro by means of α-ketoamide and pyridone-containing pharmaceuticals using in silico molecular docking

UCSF Chimera-A visualization system for exploratory research and analysis

Automatic atom type and bond type perception in molecular mechanical calculations

Open Babel: An open chemical toolbox

A smoothed backbone-dependent rotamer library for proteins derived from adaptive kernel density estimates and regressions

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

Caffeine and caffeine-containing pharmaceuticals as promising inhibitors for 3-chymotrypsin-like protease of SARS-CoV-2

ACPYPE-Antechamber python parser interface

Development and testing of a general amber force field

Improving the accuracy of protein side chain and backbone parameters from ff99SB

Comparison of simple potential functions for simulating liquid water

NAMD: A parallel, object-oriented molecular dynamics program

Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants

A smooth particle mesh Ewald method

VMD: Visual molecular dynamics

py: An efficient program for end-state free energy calculations

Software for processing and analysis of molecular dynamics trajectory data

Conformational entropy of biomolecules: Beyond the quasi-harmonic approximation

Dynamically-driven enhancement of the catalytic machinery of the SARS 3C-like protease by the S284-T285-I286/A mutations on the extra domain

Update on a promising and widely available agent for neuroprotection and neurorepair

Pharmacokinetics of 14C CDP-choline

Combined Deep Learning and Molecular Docking Simulations Approach Identifies Potentially Effective FDA Approved Drugs for Repurposing Against SARS-CoV-2

Virtual screening and molecular dynamics study of approved drugs as inhibitors of spike protein S1 domain and ACE2 interaction in SARS-CoV-2

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

Nucleotide analogues as inhibitors of SARS-CoV-2 polymerase, a key drug target for COVID-19

Action of dipeptidyl peptidase-4 inhibitors on SARS-CoV-2 main protease

Probing antiviral drugs against SARS-CoV-2 through virus-drug association prediction based on the KATZ method

lycorine and oxysophoridine inhibit novel coronavirus (SARS-CoV-2) in cell culture

Influence of 4 -Substitution on the Activity of Gemcitabine and Its ProTide against VZV and SARS-CoV-2

Repurposing potential of FDA-approved and investigational drugs for COVID-19 targeting SARS-CoV-2 spike and main protease and validation by machine learning algorithm

Nucleoside analogs and nucleoside precursors as drugs in the fight against SARS-CoV-2 and other coronaviruses

A Study of 3CLpros as Promising Targets against SARS-CoV and SARS-CoV-2. Microorganisms

Recent discovery and development of inhibitors targeting coronaviruses

Virtual screening, ADME/T, and binding free energy analysis of anti-viral, anti-protease, and anti-infectious compounds against NSP10/NSP16 methyltransferase and main protease of SARS CoV-2

Virtual screening and repurposing of FDA approved drugs against COVID-19 main protease

Highthroughput screening of the ReFRAME, Pandemic Box, and COVID Box drug repurposing libraries against SARS-CoV-2 nsp15 endoribonuclease to identify small-molecule inhibitors of viral activity

Drug repurposing to identify therapeutics against COVID 19 with SARS-Cov-2 spike glycoprotein and main protease as targets: An in silico study

COVID19 Approved Drug Repurposing: Pocket Similarity Approach

Virtual screening of natural products against type II transmembrane serine protease (TMPRSS2), the priming agent of coronavirus 2 (SARS-CoV-2)

Institutional Review Board Statement: Not applicable.

Data Availability Statement: The data will be available upon request.

There is no conflict of interest to be reported.Sample Availability: Not applicable.