Subtle changes at the RBD/hACE2 interface during SARS-CoV2 variant evolution: a molecular dynamics study

Aria Gheeraert,Vincent Leroux,Dominique Mias-Lucquin,Yasaman Karami,Laurent Vuillon,Isaure Chauvot de Bauchêne,Marie-Dominique Devignes,Ivan Rivalta,Bernard Maigret,Laurent Chaloin
DOI: https://doi.org/10.1101/2024.12.12.628120
2024-12-13
Abstract:The SARS-CoV-2 Omicron variants present a different behavior compared to the previous variants, all particularly in respect to the Delta variant, as it seems to promote a lower morbidity although being much more contagious. In this perspective, we performed new molecular dynamics (MD) simulations of the various spike RBD/hACE2 complexes corresponding to the WT, Delta and Omicron variants (BA.1 up to BA.4/5) over 1.5 μs timescale. Then, carrying out a comprehensive analysis of residue interactions within and between the two partners, allowed us to draw the profile of each variant by using complementary methods (PairInt, hydrophobic potential, contact PCA). Main results of PairInt calculations highlighted the most involved residues in electrostatic interactions that represent a strong contribution in the binding with highly stable contacts between spike RBD and hACE2 (importance of mutated residues at positions 417, 493 and 498). In addition to the swappable arginine residues (493/498), the apolar contacts made a substantial and complementary contribution in Omicron with the detection of two hydrophobic patches, one of which was correlated with energetic contribution calculations. This study brings new highlights on the global dynamics of spike RBD/hACE2 complexes resulting from the analysis of contact networks and cross-correlation matrices able to detect subtle changes at point mutations. The results of our study are also consistent with alternative approaches such as binding free energy calculations but are more informative and sensitive to transient or low-energy interactions. Nevertheless, the energetic contributions of residues at positions 501 and 505 were in good agreement with hydrophobic interactions measurements. The contact PCA networks could identify the intramolecular incidence of the S375F mutation occurring in all Omicron variants and likely conferring them an advantage in binding stability. Collectively, these data revealed the major differences observed between WT/Delta and Omicron variants at the RBD/hACE2 interface, which may explain the greater persistence of Omicron.
Bioinformatics
What problem does this paper attempt to address?