The α4β2 nAChR is implicated in a range of diseases and disorders including nicotine addiction, epilepsy and Parkinson's and Alzheimer's diseases. Designing α4β2 nAChR selective inhibitors could help define the role of the α4β2 nAChR in such disease states. In this study, we aimed to modify globular and ribbon α-conotoxin GID to selectively target the α4β2 nAChR through competitive inhibition of the α4(+)β2(-) or α4(+)α4(-) interfaces. The binding modes of the globular α-conotoxin [γ4E]GID with rat α3β2, α4β2 and α7 nAChRs were deduced using computational methods and were validated using published experimental data. The binding mode of globular [γ4E]GID at α4β2 nAChR can explain the experimental mutagenesis data, suggesting that it could be used to design GID variants. The predicted mutational energy results showed that globular [γ4E]GID is optimal for binding to α4β2 nAChR and its activity could not likely be further improved through amino-acid substitutions. The binding mode of ribbon GID with the (α4)3(β2)2 nAChR was deduced using the information from the cryo-electron structure of (α4)3(β2)2 nAChR and the binding mode of ribbon AuIB. The program FoldX predicted the mutational energies of ribbon [γ4E]GID at the α4(+)α4(-) interface, and several ribbon[γ4E]GID mutants were suggested to have desirable properties to inhibit (α4)3(β2)2 nAChR.
Keywords: FoldX; drug design; molecular dynamics simulation; α-conotoxin; α4β2 nAChR.