Characterization and Fluctuations of an Ivermectin Binding Site at the Lipid Raft Interface of the N-Terminal Domain (NTD) of the Spike Protein of SARS-CoV-2 Variants

Most studies on the docking of ivermectin on the spike protein of SARS-CoV-2 concern the receptor binding domain (RBD) and, more precisely, the RBD interface recognized by the ACE2 receptor. The N-terminal domain (NTD), which controls the initial attachment of the virus to lipid raft gangliosides, has not received the attention it deserves. In this study, we combined molecular modeling and physicochemical approaches to analyze the mode of interaction of ivermectin with the interface of the NTD-facing lipid rafts of the host cell membrane. We characterize a binding area that presents point mutations and deletions in successive SARS-CoV-2 variants from the initial strain to omicron KP.3 circulating in many countries in 2024. We show that ivermectin has exceptional flexibility, allowing the drug to bind to the spike protein of all variants tested. The energy of interaction is specific to each variant, allowing a classification according to their affinity for ivermectin in the following ascending order: Omicron KP.3 < Delta < Omicron BA.5 < Alpha < Wuhan (B.1) < Omicron BA.1. The binding site of ivermectin is subject to important variations of the NTD, including the Y144 deletion. It overlaps with the ganglioside binding domain of the NTD, as demonstrated by docking and physicochemical studies. These results suggest a new mechanism of antiviral action for ivermectin based on competitive inhibition for initial virus attachment to lipid rafts. The current KP.3 variant is still recognized by ivermectin, although with an affinity slightly lower than the Wuhan strain.