Orientation of μ-conotoxin PIIIA in a sodium channel vestibule, based on voltage dependence of its binding

In Silico Conotoxin Studies: Progress and Prospects

Cone snails of the genus Conus have evolved to produce structurally distinct and functionally diverse venom peptides for defensive and predatory purposes. This nature-devised delicacy enlightened drug discovery and for decades, the bioactive cone snail venom peptides, known as conotoxins, have been widely explored for their therapeutic potential, yet we know very little about them. With the augmentation of computational algorithms from the realms of bioinformatics and machine learning, in silico strategies have made substantial contributions to facilitate conotoxin studies although still with certain limitations. In this review, we made a bibliometric analysis of in silico conotoxin studies from 2004 to 2024 and then discussed in silico strategies to not only efficiently classify conotoxin superfamilies but also speed up drug discovery from conotoxins, reveal binding modes of known conotoxin-ion channel interactions at a microscopic level and relate the mechanisms of ion channel modulation to its underlying molecular structure. We summarized the current progress of studies in this field and gave an outlook on prospects.

Historical Perspective of the Characterization of Conotoxins Targeting Voltage-Gated Sodium Channels

Marine toxins have potent actions on diverse sodium ion channels regulated by transmembrane voltage (voltage-gated ion channels) or by neurotransmitters (nicotinic acetylcholine receptor channels). Studies of these toxins have focused on varied aspects of venom peptides ranging from evolutionary relationships of predator and prey, biological actions on excitable tissues, potential application as pharmacological intervention in disease therapy, and as part of multiple experimental approaches towards an understanding of the atomistic characterization of ion channel structure. This review examines the historical perspective of the study of conotoxin peptides active on sodium channels gated by transmembrane voltage, which has led to recent advances in ion channel research made possible with the exploitation of the diversity of these marine toxins.

Extremely Potent Block of Bacterial Voltage-Gated Sodium Channels by µ-Conotoxin PIIIA

µ-Conotoxin PIIIA, in the sub-picomolar, range inhibits the archetypal bacterial sodium channel NaChBac (NavBh) in a voltage- and use-dependent manner. Peptide µ-conotoxins were first recognized as potent components of the venoms of fish-hunting cone snails that selectively inhibit voltage-gated skeletal muscle sodium channels, thus preventing muscle contraction. Intriguingly, computer simulations predicted that PIIIA binds to prokaryotic channel NavAb with much higher affinity than to fish (and other vertebrates) skeletal muscle sodium channel (Nav 1.4). Here, using whole-cell voltage clamp, we demonstrate that PIIIA inhibits NavBac mediated currents even more potently than predicted. From concentration-response data, with [PIIIA] varying more than 6 orders of magnitude (10⁻¹² to 10⁻⁵ M), we estimated an IC₅₀ = ~5 pM, maximal block of 0.95 and a Hill coefficient of 0.81 for the inhibition of peak currents. Inhibition was stronger at depolarized holding potentials and was modulated by the frequency and duration of the stimulation pulses. An important feature of the PIIIA action was acceleration of macroscopic inactivation. Docking of PIIIA in a NaChBac (NavBh) model revealed two interconvertible binding modes. In one mode, PIIIA sterically and electrostatically blocks the permeation pathway. In a second mode, apparent stabilization of the inactivated state was achieved by PIIIA binding between P2 helices and trans-membrane S5s from adjacent channel subunits, partially occluding the outer pore. Together, our experimental and computational results suggest that, besides blocking the channel-mediated currents by directly occluding the conducting pathway, PIIIA may also change the relative populations of conducting (activated) and non-conducting (inactivated) states.

In Silico Analysis of the Subtype Selective Blockage of KCNA Ion Channels through the µ-Conotoxins PIIIA, SIIIA, and GIIIA

Understanding subtype specific ion channel pore blockage by natural peptide-based toxins is crucial for developing such compounds into promising drug candidates. Herein, docking and molecular dynamics simulations were employed in order to understand the dynamics and binding states of the µ-conotoxins, PIIIA, SIIIA, and GIIIA, at the voltage-gated potassium channels of the KV1 family, and they were correlated with their experimental activities recently reported by Leipold et al. Their different activities can only adequately be understood when dynamic information about the toxin-channel systems is available. For all of the channel-bound toxins investigated herein, a certain conformational flexibility was observed during the molecular dynamic simulations, which corresponds to their bioactivity. Our data suggest a similar binding mode of µ-PIIIA at KV1.6 and KV1.1, in which a plethora of hydrogen bonds are formed by the Arg and Lys residues within the α-helical core region of µ-PIIIA, with the central pore residues of the channel. Furthermore, the contribution of the K+ channel’s outer and inner pore loops with respect to the toxin binding. and how the subtype specificity is induced, were proposed.

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

Role of the disulfide bond on the structure and activity of μ-conotoxin PIIIA in the inhibition of NaV1.4

μ-Conotoxin PIIIA, a peptide toxin isolated from Conus purpurascens, blocks the skeletal muscle voltage-gated sodium channel NaV1.4 with significant potency. PIIIA has three disulfide bonds, which contribute largely to its highly constrained and stable structure. In this study, a combination of experimental studies and computational modeling were performed to assess the effects of deletion of the disulfide bonds on the structure and activity of PIIIA. The final results indicate that the three disulfide bonds of PIIIA are required to produce the effective inhibition of NaV1.4, and the removal of any one of the disulfide bonds significantly reduces its binding affinity owing to secondary structure variation, among which the Cys11-Cys22 is the most important for sustaining the structure and activity of PIIIA.

Predicting Structural Details of the Sodium Channel Pore Basing on Animal Toxin Studies

Eukaryotic voltage-gated sodium channels play key roles in physiology and are targets for many toxins and medically important drugs. Physiology, pharmacology, and general architecture of the channels has long been the subject of intensive research in academia and industry. In particular, animal toxins such as tetrodotoxin, saxitoxin, and conotoxins have been used as molecular probes of the channel structure. More recently, X-ray structures of potassium and prokaryotic sodium channels allowed elaborating models of the toxin-channel complexes that integrated data from biophysical, electrophysiological, and mutational studies. Atomic level cryo-EM structures of eukaryotic sodium channels, which became available in 2017, show that the selectivity filter structure and other important features of the pore domain have been correctly predicted. This validates further employments of toxins and other small molecules as sensitive probes of fine structural details of ion channels.

Determination of the μ-Conotoxin PIIIA Specificity Against Voltage-Gated Sodium Channels from Binding Energy Calculations

Voltage-gated sodium (NaV) channels generate and propagate action potentials in excitable cells, and several NaV subtypes have become important targets for pain management. The μ-conotoxins inhibit subtypes of the NaV with varied specificity but often lack of specificity to interested subtypes. Engineering the selectivity of the μ-conotoxins presents considerable complexity and challenge, as it involves the optimization of their binding affinities to multiple highly conserved NaV subtypes. In this study, a model of NaV1.4 bound with μ-conotoxin PIIIA complex was constructed using homology modeling, docking, molecular dynamic simulations and binding energy calculations. The accuracy of this model was confirmed based on the experimental mutagenesis data. The complex models of PIIIA bound with varied subtypes of NaV1.x (x = 1, 2, 3, 5, 6, 7, 8, or 9) were built using NaV1.4/PIIIA complex as a template, and refined using molecular dynamic simulations. The binding affinities of PIIIA to varied subtypes of NaV1.x (x = 1 to 9) were calculated using the Molecular Mechanics Generalized Born/Surface Area (MMGB/SA) and umbrella sampling, and were compared with the experimental values. The binding affinities calculated using MMGB/SA and umbrella sampling are correlated with the experimental values, with the former and the latter giving correlation coefficient of 0.41 (R²) and 0.68 (R²), respectively. Binding energy decomposition suggests that conserved and nonconserved residues among varied NaV subtypes have a synergistic effect on the selectivity of PIIIA.

Conotoxins That Could Provide Analgesia through Voltage Gated Sodium Channel Inhibition

Chronic pain creates a large socio-economic burden around the world. It is physically and mentally debilitating, and many sufferers are unresponsive to current therapeutics. Many drugs that provide pain relief have adverse side effects and addiction liabilities. Therefore, a great need has risen for alternative treatment strategies. One rich source of potential analgesic compounds that has emerged over the past few decades are conotoxins. These toxins are extremely diverse and display selective activity at ion channels. Voltage gated sodium (Na_V) channels are one such group of ion channels that play a significant role in multiple pain pathways. This review will explore the literature around conotoxins that bind Na_V channels and determine their analgesic potential.

Computational approaches for designing potent and selective analogs of peptide toxins as novel therapeutics

Peptide toxins provide valuable therapeutic leads for many diseases. As they bind to their targets with high affinity, potency is usually ensured. However, toxins also bind to off-target receptors, causing potential side effects. Thus, a major challenge in generating drugs from peptide toxins is ensuring their specificity for their intended targets. Computational methods can play an important role in solving such design problems through construction of accurate models of receptor-toxin complexes and calculation of binding free energies. Here we review the computational methods used for this purpose and their application to toxins targeting ion channels. We describe ShK and HsTX1 toxins, high-affinity blockers of the voltage-gated potassium channel Kv1.3, which could be developed as therapeutic agents for autoimmune diseases.

Folding similarity of the outer pore region in prokaryotic and eukaryotic sodium channels revealed by docking of conotoxins GIIIA, PIIIA, and KIIIA in a NavAb-based model of Nav1.4

Analyses of toxin binding to a homology model of Nav1.4 indicate similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.

Voltage-gated sodium channels are targets for many drugs and toxins. However, the rational design of medically relevant channel modulators is hampered by the lack of x-ray structures of eukaryotic channels. Here, we used a homology model based on the x-ray structure of the NavAb prokaryotic sodium channel together with published experimental data to analyze interactions of the μ-conotoxins GIIIA, PIIIA, and KIIIA with the Nav1.4 eukaryotic channel. Using Monte Carlo energy minimizations and published experimentally defined pairwise contacts as distance constraints, we developed a model in which specific contacts between GIIIA and Nav1.4 were readily reproduced without deformation of the channel or toxin backbones. Computed energies of specific interactions between individual residues of GIIIA and the channel correlated with experimental estimates. The predicted complexes of PIIIA and KIIIA with Nav1.4 are consistent with a large body of experimental data. In particular, a model of Nav1.4 interactions with KIIIA and tetrodotoxin (TTX) indicated that TTX can pass between Nav1.4 and channel-bound KIIIA to reach its binding site at the selectivity filter. Our models also allowed us to explain experimental data that currently lack structural interpretations. For instance, consistent with the incomplete block observed with KIIIA and some GIIIA and PIIIA mutants, our computations predict an uninterrupted pathway for sodium ions between the extracellular space and the selectivity filter if at least one of the four outer carboxylates is not bound to the toxin. We found a good correlation between computational and experimental data on complete and incomplete channel block by native and mutant toxins. Thus, our study suggests similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.

Mechanism of μ-conotoxin PIIIA binding to the voltage-gated Na+ channel NaV1.4

Several subtypes of voltage-gated Na+ (NaV) channels are important targets for pain management. μ-Conotoxins isolated from venoms of cone snails are potent and specific blockers of different NaV channel isoforms. The inhibitory effect of μ-conotoxins on NaV channels has been examined extensively, but the mechanism of toxin specificity has not been understood in detail. Here the known structure of μ-conotoxin PIIIA and a model of the skeletal muscle channel NaV1.4 are used to elucidate elements that contribute to the structural basis of μ-conotoxin binding and specificity. The model of NaV1.4 is constructed based on the crystal structure of the bacterial NaV channel, NaVAb. Six different binding modes, in which the side chain of each of the basic residues carried by the toxin protrudes into the selectivity filter of NaV1.4, are examined in atomic detail using molecular dynamics simulations with explicit solvent. The dissociation constants (Kd) computed for two selected binding modes in which Lys9 or Arg14 from the toxin protrudes into the filter of the channel are within 2 fold; both values in close proximity to those determined from dose response data for the block of NaV currents. To explore the mechanism of PIIIA specificity, a double mutant of NaV1.4 mimicking NaV channels resistant to μ-conotoxins and tetrodotoxin is constructed and the binding of PIIIA to this mutant channel examined. The double mutation causes the affinity of PIIIA to reduce by two orders of magnitude.

Computational methods of studying the binding of toxins from venomous animals to biological ion channels: theory and applications

The discovery of new drugs that selectively block or modulate ion channels has great potential to provide new treatments for a host of conditions. One promising avenue revolves around modifying or mimicking certain naturally occurring ion channel modulator toxins. This strategy appears to offer the prospect of designing drugs that are both potent and specific. The use of computational modeling is crucial to this endeavor, as it has the potential to provide lower cost alternatives for exploring the effects of new compounds on ion channels. In addition, computational modeling can provide structural information and theoretical understanding that is not easily derivable from experimental results. In this review, we look at the theory and computational methods that are applicable to the study of ion channel modulators. The first section provides an introduction to various theoretical concepts, including force-fields and the statistical mechanics of binding. We then look at various computational techniques available to the researcher, including molecular dynamics, brownian dynamics, and molecular docking systems. The latter section of the review explores applications of these techniques, concentrating on pore blocker and gating modifier toxins of potassium and sodium channels. After first discussing the structural features of these channels, and their modes of block, we provide an in-depth review of past computational work that has been carried out. Finally, we discuss prospects for future developments in the field.

Extension of Brownian dynamics for studying blockers of ion channels

We present new Brownian dynamics techniques for studying blockers of ion channels. By treating the channel as a fixed body, simulating the blocker molecules using rigid bodies, and using an implicit water force field with explicit ions, we are able to carry out fast simulations that can be used to investigate the dynamics of block and unblock, deduce binding modes, and calculate binding affinities. We test our program using the NavAb bacterial sodium channel, whose structure was recently solved (Payandeh et al. Nature, 2011, 475, 353-358) in conjunction with the μ-conotoxin PIIIA blocker. We derive an ohmic current-voltage relationship for channel permeation, calculate potentials of mean force for blocker unbinding, and deduce multiple binding modes for the blocker. Our results are shown to be compatible with other computational and experimental results. Finally, we discuss future improvements such as the inclusion of flexible side chains. After these improvements are carried out, we anticipate our program will be an extremely useful new tool that could be used to help develop new drugs to treat a range of ion-channel related diseases.

Conotoxins targeting neuronal voltage-gated sodium channel subtypes: potential analgesics?

Voltage-gated sodium channels (VGSC) are the primary mediators of electrical signal amplification and propagation in excitable cells. VGSC subtypes are diverse, with different biophysical and pharmacological properties, and varied tissue distribution. Altered VGSC expression and/or increased VGSC activity in sensory neurons is characteristic of inflammatory and neuropathic pain states. Therefore, VGSC modulators could be used in prospective analgesic compounds. VGSCs have specific binding sites for four conotoxin families: μ-, μO-, δ- and ί-conotoxins. Various studies have identified that the binding site of these peptide toxins is restricted to well-defined areas or domains. To date, only the μ- and μO-family exhibit analgesic properties in animal pain models. This review will focus on conotoxins from the μ- and μO-families that act on neuronal VGSCs. Examples of how these conotoxins target various pharmacologically important neuronal ion channels, as well as potential problems with the development of drugs from conotoxins, will be discussed.

Design of bioactive peptides from naturally occurring μ-conotoxin structures

To date, cone snail toxins ("conotoxins") are of great interest in the pursuit of novel subtype-selective modulators of voltage-gated sodium channels (Na(v)s). Na(v)s participate in a wide range of electrophysiological processes. Consequently, their malfunctioning has been associated with numerous diseases. The development of subtype-selective modulators of Na(v)s remains highly important in the treatment of such disorders. In current research, a series of novel, synthetic, and bioactive compounds were designed based on two naturally occurring μ-conotoxins that target Na(v)s. The initial designed peptide contains solely 13 amino acids and was therefore named "Mini peptide." It was derived from the μ-conotoxins KIIIA and BuIIIC. Based on this Mini peptide, 10 analogues were subsequently developed, comprising 12-16 amino acids with two disulfide bridges. Following appropriate folding and mass verification, blocking effects on Na(v)s were investigated. The most promising compound established an IC(50) of 34.1 ± 0.01 nM (R2-Midi on Na(v)1.2). An NMR structure of one of our most promising compounds was determined. Surprisingly, this structure does not reveal an α-helix. We prove that it is possible to design small peptides based on known pharmacophores of μ-conotoxins without losing their potency and selectivity. These data can provide crucial material for further development of conotoxin-based therapeutics.

Architecture and pore block of eukaryotic voltage-gated sodium channels in view of NavAb bacterial sodium channel structure

The X-ray structure of the bacterial sodium channel NavAb provides a new template for the study of sodium and calcium channels. Unlike potassium channels, NavAb contains P2 helices in the outer-pore region. Because the sequence similarity between eukaryotic and prokaryotic sodium channels in this region is poor, the structural similarity is unclear. We analyzed it by using experimental data on tetrodotoxin block of sodium channels. Key tetrodotoxin-binding residues are outer carboxylates in repeats I, II, and IV, three positions downstream from the selectivity-filter residues. In a NavAb-based model of Nav1 channels derived from the sequence alignment without insertions/deletions, the outer carboxylates did not face the pore and therefore did not interact with tetrodotoxin. The hypothesis that the evolutionary appearance of Nav1 channels involved point deletions in an ancestral channel between the selectivity filter and the outer carboxylates allowed building of a NavAb-based model with tetrodotoxin-channel contacts similar to those proposed previously. This hypothesis also allowed us to reproduce in Nav1 the folding-stabilizing contacts between long-side chain residues in P1 and P2, which are seen in NavAb. The NavAb-based inner-pore model of Nav1 preserved major features of our previous KcsA-based models, including the access pathway for ligands through the repeat III/IV interface and their interactions with specific residues. Thus, structural properties of eukaryotic voltage-gated sodium channels that are suggested by functional data were reproduced in the NavAb-based models built by using the unaltered template structure but with adjusted sequence alignment. Sequences of eukaryotic calcium channels aligned with NavAb without insertions/deletions, which suggests that NavAb is a promising basis for the modeling of calcium channels.

Binding modes of μ-conotoxin to the bacterial sodium channel (NaVAb)

Polypeptide toxins isolated from the venom of cone snails, known as μ-conotoxins, block voltage-gated sodium channels by physically occluding the ion-conducting pathway. Using molecular dynamics, we show that one subtype of μ-conotoxins, PIIIA, effectively blocks the bacterial voltage-gated sodium channel Na(V)Ab, whose crystal structure has recently been elucidated. The spherically shaped toxin, carrying a net charge of +6 e with six basic residues protruding from its surface, is attracted by the negatively charged residues on the vestibular wall and the selectivity filter of the channel. The side chain of each of these six arginine and lysine residues can wedge into the selectivity filter, whereas the side chains of other basic residues form electrostatic complexes with two acidic residues on the channel. We construct the profile of potential of mean force for the unbinding of PIIIA from the channel, and predict that PIIIA blocks the bacterial sodium channel with subnanomolar affinity.

Interactions of key charged residues contributing to selective block of neuronal sodium channels by μ-conotoxin KIIIA

Voltage-gated sodium channels are important in initiating and propagating nerve impulses in various tissues, including cardiac muscle, skeletal muscle, the brain, and the peripheral nerves. Hyperexcitability of these channels leads to such disorders as cardiac arrhythmias (Na(v)1.5), myotonias (Na(v)1.4), epilepsies (Na(v)1.2), and pain (Na(v)1.7). Thus, there is strong motivation to identify isoform-specific blockers and the molecular determinants underlying their selectivity among these channels. μ-Conotoxin KIIIA blocks rNa(v)1.2 (IC(50), 5 nM), rNa(v)1.4 (37 nM), and hNa(v)1.7 (97 nM), expressed in mammalian cells, with high affinity and a maximal block at saturating concentrations of 90 to 95%. Mutations of charged residues on both the toxin and channel modulate the maximal block and/or affinity of KIIIA. Two toxin substitutions, K7A and R10A, modulate the maximal block (52-70%). KIIIA-H12A and R14A were the only derivatives tested that altered Na(v) isoform specificity. KIIIA-R14A showed the highest affinity for Na(v)1.7, a channel involved in pain signaling. Wild-type KIIIA has a 2-fold higher affinity for Na(v)1.4 than for Na(v)1.7, which can be attributed to a missing outer vestibule charge in domain III of Na(v)1.7. Reciprocal mutations Na(v)1.4 D1241I and Na(v)1.7 I1410D remove the affinity differences between these two channels for wild-type KIIIA without affecting their affinities for KIIIA-R14A. KIIIA is the first μ-conotoxin to show enhanced activity as pH is lowered, apparently resulting from titration of the free N terminus. Removal of this free amino group reduced the pH sensitivity by 10-fold. Recognition of these molecular determinants of KIIIA block may facilitate further development of subtype-specific, sodium channel blockers to treat hyperexcitability disorders.