Structural conservation of the gabapentinoid binding site in human and Caenorhabditis elegans α2δ subunits: a docking and molecular dynamics perspective
Listed in
This article is not in any list yet, why not save it to one of your lists.Abstract
Pain is a global health burden, highlighting the need for effective therapeutic strategies. Gabapentin (GBP) and pregabalin (PGB), used for neuropathic pain, act primarily through α2δ auxiliary subunits of voltage-gated calcium channels. Caenorhabditis elegans expresses UNC-36, an ortholog of mammalian α2δ proteins involved in calcium-channel function and nocifensive behavior. However, whether UNC-36 preserves the molecular features required for gabapentinoid recognition remains unclear. We compared human α2δ-1 and UNC-36 using sequence and structural analyses, molecular docking, 500-ns molecular dynamics simulations, interaction profiling, principal component and free-energy landscape analyses, and MM/GBSA calculations. UNC-36 preserved the overall architecture of the mammalian gabapentinoid-binding region despite substantial sequence divergence, and both ligands remained associated with the modeled pockets. However, residue-level interaction networks differed between species. Human α2δ-1 showed greater contributions from aromatic interactions, whereas UNC-36 relied more prominently on cationic and hydrogen-bond donor interactions mediated by Arg501 and Arg503. The human α2δ-1–PGB complex maintained the most stable ligand pose, whereas α2δ-1–GBP showed greater positional variation. In UNC-36, PGB exhibited greater deviation from its initial binding pose than GBP. The first two principal components accounted for more conformational variance in UNC-36 complexes than in human complexes. MM/GBSA estimates showed that PGB was energetically favored over GBP in human α2δ-1, whereas GBP was favored over PGB in UNC-36. These findings show that conservation of the gabapentinoid-binding architecture is accompanied by species-specific differences in interaction chemistry, conformational dynamics, and estimated binding energetics, providing a molecular basis for interpreting C. elegans gabapentinoid responses in a translational context.