Molecular modeling and structural analysis of nAChR variants uncovers the mechanism of resistance to snake toxins

Nicotinic acetylcholine receptors (nAChRs) are neuromuscular proteins responsible for muscle contraction upon binding with chemical stimulant acetylcholine (ACh). The α-neurotoxins of snake mimic the structure of ACh and attacks nAChRs, which block the flow of ACh and leads to numbness and paralysis. The toxin-binding site of alpha subunit in the nAChRs is highly conserved throughout chordate lineages with few exceptions in resistance organisms. In this study, we have analyzed the sequence and structures of toxin-binding/resistant nAChRs and their interaction stability with toxins through molecular docking and molecular dynamics simulation (MDS). We have reported the potential glycosylation residues within the toxin-binding cleft adding sugar moieties through N-linked glycosylation in resistant organisms. Residue variations at key positions alter the secondary structure of binding cleft, which might interfere with toxin binding and it could be one of the possible explanations for the resistance to snake venoms. Analysis of nAChR-α-neurotoxin complexes has confirmed the key interacting residues. In addition, drastic variation in the binding stability of Mongoose nAChR-α-Bungarotoxin (α-BTX) and human nAChR-α-BTX complexes were found at specific phase of MDS. Our findings suggest that specific mutations in the binding site of toxin are potentially preventing the formation of stable complex of receptor-toxin, which might lead to mechanism of resistance. This in silico study on the binding cleft of nAChR and the findings of interacting residues will assist in designing potential inhibitors as therapeutic targets.