Abstract
RNA-dependent RNA polymerase (RdRp), is an enzyme essential component in the RNA replication within the life cycle of the severely acute respiratory coronavirus-2 (SARS-CoV-2), causing the deadly respiratory induced sickness COVID-19. Remdesivir is a prodrug that has seen some success in inhibiting this enzyme, however there is still the pressing need for effective alternatives. In this study, we present the discovery of four non-nucleoside small molecules that bind favorably to RdRp over adenosine-triphosphate (ATP) and active-form remdesivir-triphosphate (RTP) using high-throughput virtual screening (HTVS) coupled with extensive (total 4800 ns) molecular dynamics (MD) simulations with using the ZINC compounds database against SARS-CoV-2 RdRp (PDB: 7BV2). We found that the simulations with both ATP and RTP remained stable for the duration of their trajectories, and it was revealed that the phosphate tail of RTP was stabilized by a positive amino acid pocket near the entry channel of RTP and magnesium ions containing residues K551, R553, R555 and K621. It was also found that residues D623, D760, and N691 further stabilized the ribose portion of RTP with U10 on the template RNA strand forming hydrogen pairs with the adenosine motif. Using these models of RdRp, we employed them to screen the ZINC database of ~17 million molecules. Using docking and drug properties scoring, we narrowed down our selection to fourteen candidates. These were subjected to 200 ns simulations each underwent free energy calculations. We identified four hit compounds from the ZINC database that have similar binding poses to RTP while possessing lower overall binding free energies, with ZINC097971592 having a binding free energy two times lower than RTP.