Mechanism of molecular recognition revealed through dynamic drug binding pathways to SARS-CoV-2 main protease

Characterization of drug-binding pathways remains experimentally limited by transient intermediates and computationally challenging due to long timescales intractable for conventional molecular dynamics. To address these challenges, we combined solution NMR titrations with weighted ensemble (WE) enhanced sampling simulations to resolve atomistic pathways of nirmatrelvir binding to the SARS-CoV-2 main protease. NMR titration revealed residue-dependent heterogeneity spanning fast, intermediate, and slow exchange regimes. WE simulations complement the NMR by providing insights into unassigned residues and adding time-resolved and three-dimensional structural context. We map key interactions along two distinct binding pathways, provide dynamic explanations for residues involved in resistance, and capture unique backbone conformations compared to those sampled in unbound or bound states. Our comprehensive binding model is consistent with a combined conformational selection and induced fit mechanism in which early transient contacts are made with residues E47 and L50 and allosteric motions are centered around residue V204 of the distal domain. This synergistic application of WE and titration NMR enables a more comprehensive characterization of drug binding than either method alone, providing an integrated framework that may have broader applicability to defining structure-kinetic relationships and guiding design of next-generation inhibitors.