Abstract
The clinical durability of SARS-CoV-2 main protease (M(pro)) inhibitors depends on their resilience to emerging resistance mutations. Recent genomic surveillance and functional reports have highlighted substitutions at positions 49, 165, and 301, raising questions about the robustness of the noncovalent inhibitor WU-04 in variant backgrounds. Here, we combined μs-scale, triplicate molecular dynamics simulations with end-state binding free energy estimates and a network-rewiring inference (NRI) framework that maps long-range dynamical communication across the full protease dimer. We evaluated wild type (WT), single mutants M49K, M165V, S301P, and selected double mutants (M49K & M165V, M49K & S301P). Relative to WT, single substitutions produced reductions in computed binding affinity of up to ~12kcal/mol, accompanied by loss or reshaping of the S2 subsite and altered ligand burial. Notably, the M49K/S301P double mutant partially restored WU-04 engagement, narrowing the ΔΔG(restore) gap to within ΔΔG(restore) of WT and re-establishing key hydrophobic and hydrogen-bond contacts. NRI analysis revealed that distal residue 301 participates in a communication corridor linking the C-terminal helical domain to the active-site cleft; its substitution rewires inter-domain coupling that can compensate for local disruptions at residue 49. Together, these results identify structural hotspots and network pathways that may inform the design of next-generation M(pro) inhibitors with improved mutation tolerance-specifically by strengthening interactions that do not rely solely on the mutable S2 pocket and by engaging conserved backbone features near the 165-166 region.