Abstract
Molnupiravir, an FDA-approved nucleoside prodrug for treating COVID-19, converts into N4-hydroxycytidine triphosphate (NHC-TP), which integrates into SARS-CoV-2 RNA by its RNA-dependent RNA polymerase (RdRp) causing lethal mutations in viral proteins. Due to the risk of RdRp-mediated drug resistance and potential off-target effects on host polymerases (e.g., human polymerase II/HPolII), it is crucial to understand NHC-TP interactions at polymerase active sites for developing new, resistance-proof treatments. In this study, we used molecular dynamics (MD) simulations to probe key interactions between NHC-TP and SARS-CoV-2 RdRp and designed novel NHC-TP analogues with greater selectivity for SARS-CoV-2 RdRp over HPolII by a virtual screening workflow. We docked NHC-TP to a modified SARS-CoV-2 RdRp-Remdesivir triphosphate structure (PDB ID: 7BV2) and generated 71 NHC-TP analogues with bulky substituents to increase the interaction with RdRP and to reduce HPolII incorporation. MD simulations assessed the stability, binding affinity, and site interactions of these analogues. The top 7 candidates, with favorable ADMET properties, likely inhibit replication via potential dual mechanisms (the replicative stalling and the induction of lethal mutations) while maintaining selectivity for SARS-CoV-2 RdRp.