Abstract
The ionotropic γ-aminobutyric acid (GABA) receptor (GABAR) is a key target for the development of antiparasitic agents, particularly against ectoparasites, such as fleas and ticks. Binding stability and selectivity of sarolaner enantiomers for Ctenocephalides felis RDL receptors (RDLR) were investigated in the current study. Wild-type (WT) C. felis RDLR and its A285S mutant were constructed using homology-based, fragment-based threading and AI-driven approaches, of which, SWISS-MODEL generated the most reliable structures. Molecular docking showed that the sarolaner S-enantiomer had higher binding affinity for both receptors than the R-enantiomer, primarily due to hydrogen bonding with Ile256, π-π stacking with Phe326, and hydrophobic interactions with Ile267 and Ile268. Molecular dynamics simulations confirmed the binding stability of the S-enantiomer-receptor complex in which key residues maintained interactions throughout the trajectories. Binding free energy analysis supported these results and highlighted the role of nonpolar interactions in binding stability. The A285S mutation had minimal impact on the binding pocket, and the S-enantiomer remained selective for and bound to the mutant receptor. Insights into the insecticidal mechanism of sarolaner enantiomers are given, and the current findings may inform the development of veterinary drugs from novel isoxazoline-based NAMs targeting insect GABARs.