Toward high-resolution modeling of small molecule-ion channel interactions

Ion channels are critical drug targets for a range of pathologies, such as epilepsy, pain, itch, autoimmunity, and cardiac arrhythmias. To develop effective and safe therapeutics, it is necessary to design small molecules with high potency and selectivity for specific ion channel subtypes. There has been increasing implementation of structure-guided drug design for the development of small molecules targeting ion channels. We evaluated the performance of two RosettaLigand docking methods, RosettaLigand and GALigandDock, on the structures of known ligand-cation channel complexes. Ligands were docked to voltage-gated sodium (Na_V), voltage-gated calcium (Ca_V), and transient receptor potential vanilloid (TRPV) channel families. For each test case, RosettaLigand and GALigandDock methods frequently sampled a ligand-binding pose within a root mean square deviation (RMSD) of 1-2 Å relative to the experimental ligand coordinates. However, RosettaLigand and GALigandDock scoring functions cannot consistently identify experimental ligand coordinates as top-scoring models. Our study reveals that the proper scoring criteria for RosettaLigand and GALigandDock modeling of ligand-ion channel complexes should be assessed on a case-by-case basis using sufficient ligand and receptor interface sampling, knowledge about state-specific interactions of the ion channel, and inherent receptor site flexibility that could influence ligand binding.