(1) Background: SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target to fight COVID-19, and many RdRp inhibitors nucleotide/nucleoside analogs, such as remdesivir, have been identified or are in clinical studies. However, the appearance of resistant mutations could reduce their efficacy. In the present work, we structurally evaluated the impact of RdRp mutations found at baseline in 39 patients treated with remdesivir and associated with a different degree of antiviral response in vivo. (2) Methods: A refined bioinformatics approach was applied to assign SARS-CoV-2 clade and lineage, and to define RdRp mutational profiles. In line with such a method, the same mutations were built and analyzed by combining docking and thermodynamics evaluations with both molecular dynamics and representative pharmacophore models. (3) Results: Clinical studies revealed that patients bearing the most prevalent triple mutant P323L+671S+M899I, which was present in 41% of patients, or the more complex mutational profile P323L+G671S+L838I+D738Y+K91E, which was found with a prevalence of 2.6%, showed a delayed reduced response to remdesivir, as confirmed by the increase in SARS-CoV-2 viral load and by a reduced theoretical binding affinity versus RdRp (ΔGbindWT = -122.70 kcal/mol; ΔGbindP323L+671S+M899I = -84.78 kcal/mol; ΔGbindP323L+G671S+L838I+D738Y+K91E = -96.74 kcal/mol). Combined computational approaches helped to rationalize such clinical observations, offering a mechanistic understanding of the allosteric effects of mutants on the global motions of the viral RNA synthesis machine and in the changes of the interactions patterns of remdesivir during its binding.
Keywords: PCA; RNA-dependent RNA polymerase; SARS-CoV-2; molecular dynamics; mutations; pharmacophore; remdesivir.