This study demonstrates the inhibitory effect of 42 pyrimidonic pharmaceuticals (PPs) on the 3-chymotrypsin-like protease of SARS-CoV-2 (3CLpro) 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 3CLpro 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 3CLpro. The free energy of binding also demonstrates the stability of the PP-3CLpro 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.
Keywords: 3-chymotrypsin-like protease; COVID-19; binding free energy; coronavirus SARS-CoV-2; molecular dynamics simulations; pyrimidonic pharmaceuticals.