Abstract
The unusual and interesting architecture of the catalytic chamber of the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) was recently explored using Cryogenic Electron Microscopy (Cryo-EM), which revealed the presence of two distinctive binding cavities within the catalytic chamber. In this report, first, we mapped out and fully characterized the variations between the two binding sites, BS1 and BS2, for significant differences in their amino acid architecture, size, volume, and hydrophobicity. This was followed by investigating the preferential binding of eight antiviral agents to each of the two binding sites, BS1 and BS2, to understand the fundamental factors that govern the preferential binding of each drug to each binding site. Results showed that, in general, hydrophobic drugs, such as remdesivir and sofosbuvir, bind better to both binding sites than relatively less hydrophobic drugs, such as alovudine, molnupiravir, zidovudine, favilavir, and ribavirin. However, suramin, which is a highly hydrophobic drug, unexpectedly showed overall weaker binding affinities in both binding sites when compared to other drugs. This unexpected observation may be attributed to its high binding solvation energy, which disfavors overall binding of suramin in both binding sites. On the other hand, hydrophobic drugs displayed higher binding affinities towards BS1 due to its higher hydrophobic architecture when compared to BS2, while less hydrophobic drugs did not show a significant difference in binding affinities in both binding sites. Analysis of binding energy contributions revealed that the most favorable components are the ΔE
, ΔE
, and ΔG
, whereas ΔG
was unfavorable. The ΔE
and ΔG
for hydrophobic drugs were enough to balance the unfavorable ΔG
, leaving the ΔE
to be the most determining factor of the total binding energy. The information presented in this report will provide guidelines for tailoring SARS-CoV-2 inhibitors with enhanced binding profiles.