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.

The continued evolution of SARS-CoV-2 has enabled an escape from most monoclonal antibodies, yet a subset of broadly neutralizing antibodies targeting three newly identified super-conserved RBD epitopes—SCORE-A, SCORE-B, and SCORE-C—retains remarkable activity against even the most recent JN.1-derived sublineages. Here, we employed an integrated computational framework combining conformational dynamics, mutational scanning, MM-GBSA binding energetics, and frustration profiling to dissect the molecular mechanisms by which XGI antibodies achieve broad neutralization and resistance to immune escape. Structural analysis revealed that all three SCORE epitopes share a common architecture: a highly conserved, minimally frustrated core that provides stable anchoring, flanked by peripheral regions that accommodate antibody-specific variations. Conformational dynamics showed that SCORE-A antibodies (XGI-183) rigidify the lateral epitope while leaving the RBM partially mobile; SCORE-B antibodies (XGI-198, XGI-203) clamp the RBM apex, directly blocking ACE2; and SCORE-C antibodies (XGI-171) allosterically loosen the RBM loop, impairing receptor engagement indirectly. Mutational scanning identified a hierarchical hotspot organization where primary hotspots (e.g., K356, T500, Y380, T385) are evolutionarily constrained and minimally frustrated, while secondary hotspots (e.g., V503, Y508, S383) are neutrally frustrated and represent the principal sites of immune-driven mutations. MM-GBSA decomposition revealed that van der Waals-driven hydrophobic packing dominates binding, with electrostatic interactions providing auxiliary stabilization. Critically, frustration analysis demonstrated that immune escape hotspots reside precisely in zones of neutral frustration—“energetic playgrounds” that permit mutational exploration without destabilizing the RBD—while minimally frustrated cores are evolutionarily locked. The comparative analysis of conformational versus mutational frustration distributions revealed a unifying principle: aligned neutral frustration yields permissive, escape-prone interfaces; decoupling enables the targeting of constrained cores; and the convergence of minimal frustration in both distributions creates invulnerable interfaces. These findings establish that broad neutralization arises not from ultra-high-affinity anchors but from strategic energy distribution across rigid, evolutionarily informed interfaces, providing a roadmap for designing next-generation therapeutics that target the invulnerable cores of viral surface proteins.

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 coarse-grained simulations and atomistic reconstruction of conformational landscapes together with mutational scanning of the binding interfaces, dynamic profiling of binding and immune escape using molecular mechanics generalized Born surface area (MM-GBSA) analysis. A central theme emerging from this work is the critical role of epitope breadth and interaction diversity in determining an antibody’s resilience to mutations. BD55-1205 antibody exemplifies the advantages of broad epitope coverage and distributed hotspot mechanisms. By engaging an extensive network of residues across the RBD, BD55-1205 minimizes its dependence on individual side-chain conformations, allowing it to maintain robust binding even when key residues are mutated. This adaptability is particularly evident in its tolerance to mutations at positions such as L455 and F456 , which severely compromise other antibodies. The ability of BD55-1205 to sustain cumulative interactions underscores the importance of targeting diverse epitopes through multiple interaction mechanisms, a strategy that enhances resistance to immune evasion while maintaining functional integrity. In contrast, BD-604 and OMI-42, with localized binding mechanisms, are more vulnerable to escape mutations at critical positions such as L455, F456, and A475. P5S-1H1 and P5S-2B10 exhibit intermediate behavior, balancing specificity and adaptability but lacking the robustness of BD55-1205. Mutational scanning identified key residues like Y421, Y489, and F456 as critical hotspots for RBD stability and antibody binding, highlighting their dual role in viral fitness and immune evasion. The computational predictions generated through mutational scanning and MM-GBSA analysis demonstrate excellent agreement with experimental data on average antibody escape scores. This study underscores the diversity of binding mechanisms employed by different antibodies and molecular basis for high affinity and excellent neutralization activity of the latest generation of antibodies.

Evolution of SARS-CoV-2 has led to the emergence of variants with increased immune evasion capabilities, posing significant challenges to antibody-based therapeutics and vaccines. The cross-neutralization activity of antibodies against Omicron variants is governed by a complex and delicate interplay of multiple energetic factors and interaction contributions. In this study, we conducted a comprehensive analysis of the interactions between the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein and four neutralizing antibodies—S309, S304, CYFN1006, and VIR-7229. Using integrative computational modeling that combined all-atom molecular dy-namics (MD) simulations, mutational scanning, and MM-GBSA binding free energy calculations, we elucidated the structural, energetic, and dynamic determinants of antibody binding. Our findings reveal distinct dynamic binding mechanisms and evolutionary adaptation driving broad neutralization effect of these antibodies. We show that S309 targets conserved residues near the ACE2 interface, leveraging synergistic van der Waals and electrostatic interactions, while S304 focuses on fewer but sensitive residues, making it more susceptible to escape mutations. The analysis of CYFN-1006.1 and CYFN-1006.2 antibody binding highlights broad epitope coverage with critical anchors at T345, K440, and T346, enhancing its efficacy against variants carrying the K356T mutation which caused escape from S309 binding. Our analysis of broadly potent VIR-7229 antibody binding to XBB.1.5 and EG.5 Omicron variants emphasized a large and structurally complex epitope, demonstrating certain adaptability and compensatory effects to F456L and L455S mutations. Mutational profiling identified key residues crucial for antibody binding, including T345, P337, and R346 for S309, and T385 and K386 for S304, underscoring their roles as evolutionary "weak spots" that balance viral fitness and immune evasion. The results of this energetic analysis demonstrate a good agreement between the predicted binding hotspots and critical mutations with respect to the latest experiments on average antibody escape scores. The results of this study dissect distinct energetic mechanisms of binding and importance of targeting conserved residues and diverse epitopes to counteract viral resistance. Broad-spectrum antibodies CYFN1006 and VIR-7229, which maintain efficacy across multiple variants achieve neutralization by targeting convergent evolution hotspots while enabling tolerance to muta-tions in these positions through structural adaptability and compensatory interactions at the binding interface. The results of this study underscore the diversity of binding mechanisms employed by different antibodies and molecular basis for high affinity and excellent neutraliza-tion activity of the latest generation of antibodies.

Understanding the atomistic basis of multi-layer mechanisms employed by broadly reactive neutralizing antibodies of the SARS-CoV-2 spike protein without directly blocking receptor engagement remains an important challenge in coronavirus immunology. Class 4 antibodies represent an intriguing case: they target a deeply conserved, cryptic epitope on the receptor-binding domain yet exhibit variable neutralization potency across subgroups F1 (CR3022, EY6A, COVA1-16), F2 (DH1047), and F3 (S2X259). The molecular basis for this variability is not fully understood. Here, we employed a multi-modal computational approach integrating atomistic and coarse-grained molecular dynamics simulations, binding free energy calculations, mutational scanning, and dynamic network analysis to elucidate how these antibodies engage the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein and influence its function. Our results reveal that neutralization efficacy arises from the interplay of direct interfacial interactions and allosteric effects. Group F1 antibodies (CR3022, EY6A, COVA1-16) primarily operate via classic allostery, modulating flexibility in RBD loop regions to indirectly interfere with the ACE2 receptor binding through long-range effects. Group F2 antibody DH1047 represents an intermediate mechanism, combining partial steric hindrance—through engagement of ACE2-critical residues T376, R408, V503, and Y508—with significant allosteric influence, facilitated by localized communication pathways linking the epitope to the receptor interface. Group F3 antibody S2X259 achieves potent neutralization through a synergistic mechanism involving direct competition with ACE2 and localized allosteric stabilization, albeit with potentially increased escape vulnerability. Dynamic network analysis identified a conserved “allosteric ring” within the RBD core that serves as a structural scaffold for long-range signal propagation, with antibody-specific extensions modulating communication to the ACE2 interface. These findings support a model where Class 4 neutralization strategies evolve through the refinement of peripheral allosteric connections rather than epitope redesign. This study establishes a robust computational framework for understanding the atomistic basis of neutralization activity and immune escape for Class 4 antibodies, highlighting how the interplay of binding energetics, conformational dynamics, and allosteric modulation governs their effectiveness against SARS-CoV-2.

KRAS is a pivotal oncoprotein that regulates cell proliferation and survival through interactions with downstream effectors such as RAF1. Despite significant advances in understanding KRAS biology, the structural and dynamic mechanisms of KRAS allostery remain poorly understood. In this study, we employ microsecond molecular dynamics simulations, mutational scanning, and binding free energy calculations together with dynamic network modeling to dissect how engineered DARPin proteins K27, K55, K13, and K19 engage KRAS through diverse molecular mechanisms ranging from effector mimicry to conformational restriction and allosteric modulation. Mutational scanning across all four DARPin systems identifies a core set of evolutionarily constrained residues that function as universal hotspots in KRAS recognition. KRAS residues I36, Y40, M67, and H95 consistently emerge as critical contributors to binding stability. Binding free energy computations show that, despite similar binding modes, K27 relies heavily on electrostatic contributions from major binding hotspots while K55 exploits a dense hydrophobic cluster enhancing its effector-mimetic signature. The allosteric binders K13 and K19, by contrast, stabilize a KRAS-specific pocket in the α3–loop–α4 motif, introducing new hinges and bottlenecks that rewire the communication architecture of KRAS without full immobilization. Network-based analysis reveals a strikingly consistent theme: despite their distinct mechanisms of recognition, all systems engage a unifying allosteric architecture that spans multiple functional motifs. This architecture is not only preserved across complexes but also mirrors the intrinsic communication framework of KRAS itself, where specific residues function as central hubs transmitting conformational changes across the protein. By integrating dynamic profiling, energetic mapping, and network modeling, our study provides a multi-scale mechanistic roadmap for targeting KRAS, revealing how engineered proteins can exploit both conserved motifs and isoform-specific features to enable precision modulation of KRAS signaling in oncogenic contexts.