Abstract
The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019 has posed a serious threat to global health and the economy for over two years, prompting the need for development of antiviral inhibitors. Due to its vital role in viral replication, RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target. Herein, we analyzed amino acid sequence conservation of RdRp across coronaviruses. The conserved amino acids at the catalytic binding site served as the ligand-contacting residues for in silico screening to elucidate possible resistant mutation. Molecular docking was employed to screen inhibitors of SARS-CoV-2 from the ZINC ChemDiv database. The top-ranked compounds selected from GOLD docking were further investigated for binding modes at the conserved residues of RdRp, and ten compounds were selected for experimental validation. Of which, three compounds exhibited promising antiviral activity. The most promising candidate showed a half-maximal effective concentration (EC50) of 5.04 µM. Molecular dynamics simulations, binding free-energy calculation and hydrogen bond analysis were performed to elucidate the critical interactions providing a foundation for developing lead compounds effective against SARS-CoV-2.