Dimerization and Ligand Binding Rewire Allosteric Networks in SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (MPro) is an essential enzyme for viral replication and a primary target for antiviral drug development. Despite extensive structural and biochemical characterization, the allosteric mechanisms by which dimerization informs conformational changes at active site lack an explicit comparison across the different states that identify key residues that connect substrate binding, dimerization, and catalytic activation. Here, we integrate microsecond time scale all-atom molecular dynamics (MD) simulations with dynamical network analysis to characterize how ligand binding and dimerization modulate the allosteric communication landscape of MPro. We performed triplicate 1-μs simulations of MPro in the monomer and dimer states. For each of these states, we simulated MPro in the apo state, as well as bound to a natural peptide substrate (nsp 15/16), the covalent inhibitor nirmatrelvir (Paxlovid) and the noncovalent inhibitor ensitrelvir (Xocova). Dynamical network analyses from the resulting simulations reveal that dimerization redirects the highest correlated motions from the interdomain loop towards the domain II and III interface. At the dimer interface, we identified N-terminal and domain II β-hairpin residues that act as central communication hubs creating networks that connect both chains in the dimer. Small molecule binding to the active site further modulates these networks in distinct ways: nirmatrelvir and peptide substrate binding results in the formation of allosteric networks within the oxyanion loop, while ensitrelvir-bound monomeric MPro results in a dimer-like network, suggesting an inhibitory "allosteric switch" mechanism that may hinder dimerization upon binding. Across all systems, domain III emerges as an allosteric "pivot", providing a platform that allows the most relevant networks to connect inter-chain communication to the active site upon dimerization. Together, these findings define how correlated motion networks couple active-site dynamics to dimerization and ligand binding, providing molecular insight into the principles governing allosteric regulation in MPro. This framework highlights potential avenues for developing antivirals that target not only the catalytic site but also the communication pathways sustaining dimer stability and enzymatic function.