A computational design approach for virtual screening of peptide interactions across K(+) channel families

Pore blocking mechanisms of centipede toxin SsTx-4 on the inwardly rectifying potassium channels

The peptide toxin SsTx-4 derived from venom of centipede Scolopendra subspinipes mutilans was characterized as a potent antagonist of the inwardly rectifying potassium (Kir) channel subtypes Kir1.1, Kir4.1, and Kir6.2 in our previous study. Alanine-scanning mutagenesis analysis identified key molecular determinants on the SsTx-4 toxin interacting with these Kir channels, as well as those on the Kir6.2 channel interacting with the toxin. However, the key residues on Kir1.1 and Kir4.1 channels responsible for binding SsTx-4 remain unclear. Here, using a combination of site-directed mutagenesis, patch-clamp analysis, molecular docking with AlphaFold 3, and molecular dynamic simulations, we revealed that SsTx-4 acted on the Kir channels as a pore blocker, with K13 on toxin serving as the functional pore-blocking residue and other residues on it contributing to stabilize the toxin-channel complex by binding to multiple residues on the wall of the channels' outer vestibule, involving E104 on Kir1.1; D100, L115, and F133 on Kir4.1; and E108, S113, H115, and M137 on Kir6.2. Collectively, these findings advanced our understanding on the interaction between Kir channels and this prototype Kir antagonist, providing insights that could inspire the development of more potent and specific Kir subtype blockers in the future.

Identification of Potential Modulators of a Pathogenic G Protein-Gated Inwardly Rectifying K+ Channel 4 Mutant: In Silico Investigation in the Context of Drug Discovery for Hypertension

Genetic abnormalities have been associated with primary aldosteronism, a major cause of secondary hypertension. This includes mutations in the KCNJ5 gene, which encodes G protein-gated inwardly rectifying K+ channel 4 (GIRK4). For example, the substitution of glycine with glutamic acid gives rise to the pathogenic GIRK4G151E mutation, which alters channel selectivity, making it more permeable to Na+ and Ca2+. While tertiapin and tertiapin-Q are well-known peptide inhibitors of the GIRK4WT channel, clinically, there is a need for the development of selective modulators of mutated channels, including GIRK4G151E. Using in silico methods, including homology modeling, protein-peptide docking, ligand-binding site prediction, and molecular docking, we aimed to explore potential modulators of GIRK4WT and GIRK4G151E. Firstly, protein-peptide docking was performed to characterize the binding site of tertiapin and its derivative to the GIRK4 channels. In accordance with previous studies, the peptide inhibitors preferentially bind to the GIRK4WT channel selectivity filter compared to GIRK4G151E. A ligand-binding site analysis was subsequently performed, resulting in the identification of two potential regions of interest: the central cavity and G-loop gate. Utilizing curated chemical libraries, we screened over 700 small molecules against the central cavity of the GIRK4 channels. Flavonoids, including luteolin-7-O-rutinoside and rutin, and the macrolides rapamycin and troleandomycin bound strongly to the GIRK4 channels. Similarly, xanthophylls, particularly luteoxanthin, bound to the central cavity with a strong preference towards the mutated GIRK4G151E channel compared to GIRK4WT. Overall, our findings suggest potential lead compounds for further investigation, particularly luteoxanthin, that may selectively modulate GIRK4 channels.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

A novel ion conducting route besides the central pore in an inherited mutant of G-protein-gated inwardly rectifying K+ channel

G-protein-gated inwardly rectifying K⁺ (GIRK; Kir3.x) channels play important physiological roles in various organs. Some of the disease-associated mutations of GIRK channels are known to induce loss of K⁺ selectivity but their structural changes remain unclear. In this study, we investigated the mechanisms underlying the abnormal ion selectivity of inherited GIRK mutants. By the two-electrode voltage-clamp analysis of GIRK mutants heterologously expressed in Xenopus oocytes, we observed that Kir3.2 G156S permeates Li⁺ better than Rb⁺ , while T154del or L173R of Kir3.2 and T158A of Kir3.4 permeate Rb⁺ better than Li⁺ , suggesting a unique conformational change in the G156S mutant. Applications of blockers of the selectivity filter (SF) pathway, Ba²⁺ or Tertiapin-Q (TPN-Q), remarkably increased the Li⁺ -selectivity of Kir3.2 G156S but did not alter those of the other mutants. In single-channel recordings of Kir3.2 G156S expressed in mouse fibroblasts, two types of events were observed, one attributable to a TPN-Q-sensitive K⁺ current and the second a TPN-Q-resistant Li⁺ current. The results show that a novel Li⁺ -permeable and blocker-resistant pathway exists in G156S in addition to the SF pathway. Mutations in the pore helix, S148F and T151A also induced high Li⁺ permeation. Our results demonstrate that the mechanism underlying the loss of K⁺ selectivity of Kir3.2 G156S involves formation of a novel ion permeation pathway besides the SF pathway, which allows permeation of various species of cations. KEY POINTS: G-protein-gated inwardly rectifying K⁺ (GIRK; Kir3.x) channels play important roles in controlling excitation of cells in various organs, such as the brain and the heart. Some of the disease-associated mutations of GIRK channels are known to induce loss of K⁺ selectivity but their structural changes remain unclear. In this study, we investigated the mechanisms underlying the abnormal ion selectivity of inherited mutants of Kir3.2 and Kir3.4. Here we show that a novel Na⁺ , Li⁺ -permeable and blocker-resistant pathway exists in an inherited mutant, Kir3.2 G156S, in addition to the conventional ion conducting pathway formed by the selectivity filter (SF). Our results demonstrate that the mechanism underlying the loss of K⁺ selectivity of Kir3.2 G156S involves formation of a novel ion permeation pathway besides the SF pathway, which allows permeation of various species of cations.

Conformational selection of vasopressin upon V1a receptor binding

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V_1aR, V_1bR and V₂R) are of high interest in a wide array of drug discovery programs. V_1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V_1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V_1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V_1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V_1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V_1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Structural Determinants Mediating Tertiapin Block of Neuronal Kir3.2 Channels

Tertiapin (TPN) is a 21 amino acid venom peptide from Apis mellifera that inhibits certain members of the inward rectifier potassium (Kir) channel family at a nanomolar affinity with limited specificity. Structure-based computational simulations predict that TPN behaves as a pore blocker; however, the molecular determinants mediating block of neuronal Kir3 channels have been inconclusive and unvalidated. Here, using molecular docking and molecular dynamics (MD) simulations with 'potential of mean force' (PMF) calculations, we investigated the energetically most favored interaction of TPN with several Kir3.x channel structures. The resulting binding model for Kir3.2-TPN complexes was then tested by targeted mutagenesis of the predicted contact sites, and their impact on the functional channel block was measured electrophysiologically. Together, our findings indicate that a high-affinity TPN block of Kir3.2 channels involves a pore-inserting lysine side chain requiring (1) hydrophobic interactions at a phenylalanine ring surrounding the channel pore and (2) electrostatic interactions with two adjacent Kir3.2 turret regions. Together, these interactions collectively stabilize high-affinity toxin binding to the Kir3.2 outer vestibule, which orients the ε-amino group of TPN-K21 to occupy the outermost K⁺ binding site of the selectivity filter. The structural determinants for the TPN block described here also revealed a favored subunit arrangement for assembled Kir3.x heteromeric channels, in addition to a multimodal binding capacity of TPN variants consistent with the functional dyad model for polybasic peptide pore blockers. These novel findings will aid efforts in re-engineering the TPN pharmacophore to develop peptide variants having unique and distinct Kir channel blocking properties.

Identification of Aethina tumida Kir Channels as Putative Targets of the Bee Venom Peptide Tertiapin Using Structure-Based Virtual Screening Methods

Venoms are comprised of diverse mixtures of proteins, peptides, and small molecules. Identifying individual venom components and their target(s) with mechanism of action is now attainable to understand comprehensively the effectiveness of venom cocktails and how they collectively function in the defense and predation of an organism. Here, structure-based computational methods were used with bioinformatics tools to screen and identify potential biological targets of tertiapin (TPN), a venom peptide from Apis mellifera (European honey bee). The small hive beetle (Aethina tumida (A. tumida)) is a natural predator of the honey bee colony and was found to possess multiple inwardly rectifying K⁺ (Kir) channel subunit genes from a genomic BLAST search analysis. Structure-based virtual screening of homology modelled A. tumida Kir (atKir) channels found TPN to interact with a docking profile and interface “footprint” equivalent to known TPN-sensitive mammalian Kir channels. The results support the hypothesis that atKir channels, and perhaps other insect Kir channels, are natural biological targets of TPN that help defend the bee colony from infestations by blocking K⁺ transport via atKir channels. From these in silico findings, this hypothesis can now be subsequently tested in vitro by validating atKir channel block as well as in vivo TPN toxicity towards A. tumida. This study highlights the utility and potential benefits of screening in virtual space for venom peptide interactions and their biological targets, which otherwise would not be feasible.

Combinatorial Assembly of Lumitoxins

Ion channels are among the most important proteins in neuroscience and serve as drug targets for many brain disorders. During development, learning, disease progression, and other processes, the activity levels of specific ion channels are tuned in a cell-type specific manner. However, it is difficult to assess how cell-specific changes in ion channel activity alter emergent brain functions. We have developed a protein architecture for fully genetically encoded light-activated modulation of endogenous ion channel activity. Fusing a genetically encoded photoswitch and an ion channel-modulating peptide toxin in a computationally designed fashion, this reagent, which we call Lumitoxins, can mediate light-modulation of specific endogenous ion channel activities in targeted cells. The modular lumitoxin architecture may be useful in a diversity of neuroscience tools. Here, we delineate how to construct lumitoxin genes from synthesized components, and provide a general outline for how to test their function in mammalian cell culture.

Discovery of peptide ligands through docking and virtual screening at nicotinic acetylcholine receptor homology models

Venom peptide toxins such as conotoxins play a critical role in the characterization of nicotinic acetylcholine receptor (nAChR) structure and function and have potential as nervous system therapeutics as well. However, the lack of solved structures of conotoxins bound to nAChRs and the large size of these peptides are barriers to their computational docking and design. We addressed these challenges in the context of the α4β2 nAChR, a widespread ligand-gated ion channel in the brain and a target for nicotine addiction therapy, and the 19-residue conotoxin α-GID that antagonizes it. We developed a docking algorithm, ToxDock, which used ensemble-docking and extensive conformational sampling to dock α-GID and its analogs to an α4β2 nAChR homology model. Experimental testing demonstrated that a virtual screen with ToxDock correctly identified three bioactive α-GID mutants (α-GID[A10V], α-GID[V13I], and α-GID[V13Y]) and one inactive variant (α-GID[A10Q]). Two mutants, α-GID[A10V] and α-GID[V13Y], had substantially reduced potency at the human α7 nAChR relative to α-GID, a desirable feature for α-GID analogs. The general usefulness of the docking algorithm was highlighted by redocking of peptide toxins to two ion channels and a binding protein in which the peptide toxins successfully reverted back to near-native crystallographic poses after being perturbed. Our results demonstrate that ToxDock can overcome two fundamental challenges of docking large toxin peptides to ion channel homology models, as exemplified by the α-GID:α4β2 nAChR complex, and is extendable to other toxin peptides and ion channels. ToxDock is freely available at rosie.rosettacommons.org/tox_dock.

Venom-derived peptides inhibiting Kir channels: Past, present, and future

Inwardly rectifying K⁺ (Kir) channels play a significant role in vertebrate and invertebrate biology by regulating the movement of K⁺ ions involved in membrane transport and excitability. Yet unlike other ion channels including their ancestral K⁺-selective homologs, there are very few venom toxins known to target and inhibit Kir channels with the potency and selectivity found for the Ca²⁺-activated and voltage-gated K⁺ channel families. It is unclear whether this is simply due to a lack of discovery, or instead a consequence of the evolutionary processes that drive the development of venom components towards their targets based on a collective efficacy to 1) elicit pain for defensive purposes, 2) promote paralysis for prey capture, or 3) facilitate delivery of venom components into the circulation. The past two decades of venom screening has yielded three venom peptides with inhibitory activity towards mammalian Kir channels, including the discovery of tertiapin, a high-affinity pore blocker from the venom of the European honey bee Apis mellifera. Venomics and structure-based computational approaches represent exciting new frontiers for venom peptide development, where re-engineering peptide 'scaffolds' such as tertiapin may aid in the quest to expand the palette of potent and selective Kir channel blockers for future research and potentially new therapeutics. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

Molecular dynamics of the honey bee toxin tertiapin binding to Kir3.2

Tertiapin (TPN), a short peptide isolated from the venom of the honey bee, is a potent and selective blocker of the inward rectifier K⁺ (Kir) channel Kir3.2. Here we examine in atomic detail the binding mode of TPN to Kir3.2 using molecular dynamics, and deduce the key residue in Kir3.2 responsible for TPN selectivity. The binding of TPN to Kir3.2 is stable when the side chain of either Lys16 (TPN^K16-Kir3.2) or Lys17 (TPN^K17-Kir3.2) of the toxin protrudes into the channel pore. However, the binding affinity calculated from only TPN^K17-Kir3.2 and not TPN^K16-Kir3.2 is consistent with experiment, suggesting that Lys17 is the most plausible pore-blocking residue. The alanine mutation of Kir3.2-Glu127, which is not present in TPN-resistant channels, reduces the inhibitory ability of TPN by over 50 fold in TPN^K17-Kir3.2, indicating that Kir3.2-Glu127 is important for the selectivity of TPN.

Computational Studies of Venom Peptides Targeting Potassium Channels

Small peptides isolated from the venom of animals are potential scaffolds for ion channel drug discovery. This review article mainly focuses on the computational studies that have advanced our understanding of how various toxins interfere with the function of K⁺ channels. We introduce the computational tools available for the study of toxin-channel interactions. We then discuss how these computational tools have been fruitfully applied to elucidate the mechanisms of action of a wide range of venom peptides from scorpions, spiders, and sea anemone.