“…75−81 A comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of S trimers for BA.1, BA.2, and BA.2.75 variants showed that Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots. 82 By combining atomistic simulations and a community-based network model of epistatic couplings, we found that convergent Omicron mutations such as G446S (BA.2.75, BA.2.75.2, XBB), F486V (BA.4, BA.5, BQ.1, BQ.1.1), F486S, F490S (XBB.1), F486P (XBB.1.5) can display epistatic relationships with the major stability and binding affinity hotspots which may allow for the observed broad antibody resistance. 83 In the current study, we perform multiple microsecond MD simulations and Markov state model (MSM) analysis to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 S RBD-ACE2 complexes for the recently emerged XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants.…”