Abstract
INTRODUCTION: The rapid evolution of SARS-CoV-2 Omicron variants highlights the urgent need for therapeutic strategies that can target viral evolution and leverage host immune recognition mechanisms. This study uses molecular dynamics (MD) simulations to analyze the immune evasion mechanisms of class 1 nanobodies against emerging SARS-CoV-2 variants, and to develop an efficient in silico pipeline for rapid affinity optimization. METHODS: We employed MD simulations and binding free energy calculations to investigate the immune evasion mechanisms of four class 1 nanobodies (R14, DL4, V(H) ab6, and Nanosota9) against wild-type (WT) and Omicron variants, including BA.2, JN.1, and KP.3/XEC. Building on these findings, we established a streamlined nanobody optimization pipeline integrating high-throughput mutagenesis of complementarity-determining regions (CDRs) and hotspot residues, protein-protein docking, and MD simulations. RESULTS: MD analysis confirmed that the immune evasion mechanism of KP.3/XEC is significantly associated with the Q493E mutation, which weakens electrostatic interactions between the nanobodies and the receptor binding domain (RBD). Through our pipeline, we identified high-affinity mutants including 3 for R14, 3 for DL4, 11 for VH ab6, and 9 for Nanosota9. The optimized R14 variant L29W/S52C/A101V demonstrated exceptional performance, achieving a 62.6% binding energy improvement against JN.1 (-76.88 kcal/mol compared to -47.3 kcal/mol for original R14 nanobody) while maintaining < 15% affinity variation across variants (compared to > 40% for original R14 nanobody). DISCUSSION: This study demonstrates that in silico affinity enhancement is a rapid and resource-efficient approach to repurpose nanobodies against SARS-CoV-2 variants, significantly accelerating affinity optimization while reducing experimental demands. This computational approach expedites the optimization of nanobody binding affinities while minimizing experimental resource requirements. By enhancing nanobody efficacy, our method provides a viable framework for developing targeted countermeasures against evolving SARS-CoV-2 variants and other pathogens.