Molecular simulation of the interaction of kappa-conotoxin-PVIIA with the Shaker potassium channel pore

Pathophysiological Responses to Conotoxin Modulation of Voltage-Gated Ion Currents

Voltage-gated ion channels are plasma membrane proteins that generate electrical signals following a change in the membrane voltage. Since they are involved in several physiological processes, their dysfunction may be responsible for a series of diseases and pain states particularly related to neuronal and muscular systems. It is well established for decades that bioactive peptides isolated from venoms of marine mollusks belonging to the Conus genus, collectively known as conotoxins, can target different types and isoforms of these channels exerting therapeutic effects and pain relief. For this reason, conotoxins are widely used for either therapeutic purposes or studies on ion channel mechanisms of action disclosure. In addition their positive property, however, conotoxins may generate pathological states through similar ion channel modulation. In this narrative review, we provide pieces of evidence on the pathophysiological impacts that different members of conotoxin families exert by targeting the three most important voltage-gated channels, such as sodium, calcium, and potassium, involved in cellular processes.

A Transcriptomic Survey of Ion Channel-Based Conotoxins in the Chinese Tubular Cone Snail (Conus betulinus)

Conotoxins in the venom of cone snails (Conus spp.) are a mixture of active peptides that work as blockers, agonists, antagonists, or inactivators of various ion channels. Recently we reported a high-throughput method to identify 215 conotoxin transcripts from the Chinese tubular cone snail, C. betulinus. Here, based on the previous datasets of four transcriptomes from three venom ducts and one venom bulb, we explored ion channel-based conotoxins and predicted their related ion channel receptors. Homologous analysis was also performed for the most abundant ion channel protein, voltage-gated potassium (Kv; with Kv1.1 as the representative), and the most studied ion channel receptor, nicotinic acetylcholine receptor (nAChR; with α2-nAChR as the representative), in different animals. Our transcriptomic survey demonstrated that ion channel-based conotoxins and related ion channel proteins/receptors transcribe differentially between the venom duct and the venom bulb. In addition, we observed that putative κ-conotoxins were the most common conotoxins with the highest transcription levels in the examined C. betulinus. Furthermore, Kv1.1 and α2-nAChR were conserved in their functional domains of deduced protein sequences, suggesting similar effects of conotoxins via the ion channels in various species, including human beings. In a word, our present work suggests a high-throughput way to develop conotoxins as potential drugs for treatment of ion channel-associated human diseases.

Vstx1, a modifier of Kv channel gating, localizes to the interfacial region of lipid bilayers

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.

Ion conduction and selectivity in K(+) channels

Potassium (K(+)) channels are tetrameric membrane-spanning proteins that provide a selective pore for the conductance of K(+) across the cell membranes. These channels are most remarkable in their ability to discriminate K(+) from Na(+) by more than a thousandfold and conduct at a throughput rate near diffusion limit. The recent progress in the structural characterization of K(+) channel provides us with a unique opportunity to understand their function at the atomic level. With their ability to go beyond static structures, molecular dynamics simulations based on atomic models can play an important role in shaping our view of how ion channels carry out their function. The purpose of this review is to summarize the most important findings from experiments and computations and to highlight a number of fundamental mechanistic questions about ion conduction and selectivity that will require further work.

Electrostatic Recognition and Induced Fit in the κ-PVIIA Toxin Binding to Shaker Potassium Channel

Brownian dynamics (BD) and molecular dynamics (MD) simulations and electrostatic calculations were performed to study the binding process of kappa-PVIIA to the Shaker potassium channel and the structure of the resulting complex. BD simulations, guided by electrostatic interactions, led to an initial alignment between the toxin and the channel protein. MD simulations were then carried out to allow for rearrangements from this initial structure. After approximately 4 ns, a critical "induced fit" process was observed to last for approximately 2 ns. In this process, the interface was reorganized, and side chains were moved so that favorable atomic contacts were formed or strengthened, while unfavorable contacts were eliminated. The final complex structure was stabilized through electrostatic interactions with the positively charged side chain of Lys7 of kappa-PVIIA deeply inserted into the channel pore and other hydrogen bonds and by hydrophobic interactions involving Phe9 and Phe23 of the toxin. The validity of the predicted structure for the complex was assessed by calculating the effects of mutating charged and polar residues of both the toxin and the channel protein, with the calculated effects correlating reasonably well with experimental data. The present study suggests a general binding mechanism, whereby proteins are pre-aligned in their diffusional encounter by long-range electrostatic attraction, and nanosecond-scale rearrangements within the initial complex then lead to a specifically bound complex.

Binding of kappa-conotoxin PVIIA to Shaker K+ channels reveals different K+ and Rb+ occupancies within the ion channel pore

The x-ray structure of the KcsA channel at different [K⁺] and [Rb⁺] provided insight into how K⁺ channels might achieve high selectivity and high K⁺ transit rates and showed marked differences between the occupancies of the two ions within the ion channel pore. In this study, the binding of κ-conotoxin PVIIA (κ-PVIIA) to Shaker K⁺ channel in the presence of K⁺ and Rb⁺ was investigated. It is demonstrated that the complex results obtained were largely rationalized by differences in selectivity filter occupancy of this 6TM channels as predicted from the structural work on KcsA. κ-PVIIA inhibition of the Shaker K⁺ channel differs in the closed and open state. When K⁺ is the only permeant ion, increasing extracellular [K⁺] decreases κ-PVIIA affinity for closed channels by decreasing the “on” binding rate, but has no effect on the block of open channels, which is influenced only by the intracellular [K⁺]. In contrast, extracellular [Rb⁺] affects both closed- and open-channel binding. As extracellular [Rb⁺] increases, (a) binding to the closed channel is slightly destabilized and acquires faster kinetics, and (b) open channel block is also destabilized and the lowest block seems to occur when the pore is likely filled only by Rb⁺. These results suggest that the nature of the permeant ions determines both the occupancy and the location of the pore site from which they interact with κ-PVIIA binding. Thus, our results suggest that the permeant ion(s) within a channel pore can determine its functional and pharmacological properties.

The binding of kappa-Conotoxin PVIIA and fast C-type inactivation of Shaker K+ channels are mutually exclusive

Kappa-conotoxin PVIIA (kappa-PVIIA), a 27-amino acid peptide identified from the venom of Conus purpurascens, inhibits the Shaker K+ channel by blocking its outer pore. The toxin appears as a gating modifier because its binding affinity decreases with relatively fast kinetics upon channel opening, but there is no indication that it interferes with the gating transitions of the wild-type channels (WT), including the structural changes of the outer pore that underlie its slow C-type inactivation. In this report we demonstrate that in two outer pore mutants of Shaker-IR (M448K and T449S), that have high toxin sensitivity and fast C-type inactivation, the latter process is instead antagonized by and incompatible with kappa-PVIIA binding. Inactivation is slowed by the necessary preliminary unbinding of kappa-PVIIA, whereas toxin rebinding must await recovery from inactivation causing a double-exponential relaxation of the second response to double-pulse stimulations. Compared with the lack of similar effects in WT, these results demonstrate the ability of peptide toxins like kappa-PVIIA to reveal possibly subtle differences in structural changes of the outer pore of K+ channels; however, they also warn against a naive use of fast inactivating mutants as models for C-type inactivation. Unfolded from the antagonistic effect of inactivation, toxin binding to mutant noninactivated channels shows state- and voltage-dependencies similar to WT: slow and high affinity for closed channels; relatively fast dissociation from open channels at rate increasing with voltage. This supports the idea that these properties depend mainly on interactions with pore-permeation processes that are not affected by the mutations. In mutant channels the state-dependence also greatly enhances the protection of toxin binding against steady-state inactivation at low depolarizations while still allowing large responses to depolarizing pulses that relieve toxin block. Although not obviously applicable to any known combination of natural channel and outer-pore blocker, our biophysical characterization of such highly efficient mechanism of protection from steady-state outer-pore inactivation may be of general interest.

Molecular modeling of ion channels: structural predictions

Recent advances in membrane protein crystallography have greatly increased structural information of channels permeating metal ions. Structural bioinformatics techniques and molecular dynamics calculations are providing structural models of ion channels for which the three-dimensional structure is not known. Most of the reported structure prediction studies focus on K(+) channels and are based on the KcsA K(+) channel structure.

Molecular dynamics simulations of a K+ channel blocker: Tc1 toxin from Tityus cambridgei

Toxins that block voltage-gated potassium (Kv) channels provide a possible template for improved homology models of the Kv pore. In assessing the interactions of Kv channels and their toxins it is important to determine the dynamic flexibility of the toxins. Multiple 10 ns duration molecular dynamics simulations combined with essential dynamics analysis have been used to explore the flexibility of four different Kv channel-blocking toxins. Three toxins (Tc1, AgTx and ChTx) share a common fold. They also share a common pattern of conformational dynamics, as revealed by essential dynamics analysis of the simulation results. This suggests that some aspects of dynamic behaviour are conserved across a single protein fold class. In each of these three toxins, the residue exhibiting minimum flexibility corresponds to a conserved lysine residue that is suggested to interact with the filter domain of the channel. Thus, comparative simulations reveal functionally important conservation of molecular dynamics as well as protein fold across a family of related toxins.