Multiscale Modeling and Dynamic Mutational Profiling of Binding Energetics and Immune Escape for Class I Antibodies with SARS-CoV-2 Spike Protein: Dissecting Mechanisms of High Resistance to Viral Escape Against Emerging Variants

The rapid evolution of SARS-CoV-2 has underscored the need for a detailed understanding of antibody binding mechanisms to combat immune evasion by emerging variants. In this study, we investigated the interactions between Class I neutralizing antibodies—BD55-1205, BD-604, OMI-42, P5S-1H1, and P5S-2B10—and the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein using multiscale modeling, which combined molecular simulations with the ensemble-based mutational scanning of the binding interfaces and binding free energy computations. A central theme emerging from this work is that the unique binding strength and resilience to immune escape of the BD55-1205 antibody are determined by leveraging a broad epitope footprint and distributed hotspot architecture, additionally supported by backbone-mediated specific interactions, which are less sensitive to amino acid substitutions and together enable exceptional tolerance to mutational escape. In contrast, BD-604 and OMI-42 exhibit localized binding modes with strong dependence on side-chain interactions, rendering them particularly vulnerable to escape mutations at K417N, L455M, F456L and A475V. Similarly, P5S-1H1 and P5S-2B10 display intermediate behavior—effective in some contexts but increasingly susceptible to antigenic drift due to narrower epitope coverage and concentrated hotspots. Our computational predictions show strong agreement with experimental deep mutational scanning data, validating the accuracy of the models and reinforcing the value of binding hotspot mapping in predicting antibody vulnerability. This work highlights that neutralization breadth and durability are not solely dictated by epitope location, but also by how binding energy is distributed across the interface. The results provide atomistic insight into mechanisms driving resilience to immune escape for broadly neutralizing antibodies targeting the ACE2 binding interface—which stems from cumulative effects of structural diversity in binding contacts, redundancy in interaction patterns and reduced vulnerability to mutation-prone positions.