Molecular dynamics simulations of scorpion toxin recognition by the Ca(2+)-activated potassium channel KCa3.1

ER-to-lysosome Ca2+ refilling followed by K+ efflux-coupled store-operated Ca2+ entry in inflammasome activation and metabolic inflammation

We studied lysosomal Ca2+ in inflammasome. Lipopolysaccharide (LPS) + palmitic acid (PA) decreased lysosomal Ca2+ ([Ca2+]Lys) and increased [Ca2+]i through mitochondrial ROS, which was suppressed in Trpm2-KO macrophages. Inflammasome activation and metabolic inflammation in adipose tissue of high-fat diet (HFD)-fed mice were ameliorated by Trpm2 KO. ER→lysosome Ca2+ refilling occurred after lysosomal Ca2+ release whose blockade attenuated LPS + PA-induced inflammasome. Subsequently, store-operated Ca2+entry (SOCE) was activated whose inhibition suppressed inflammasome. SOCE was coupled with K+ efflux whose inhibition reduced ER Ca2+ content ([Ca2+]ER) and impaired [Ca2+]Lys recovery. LPS + PA activated KCa3.1 channel, a Ca2+-activated K+ channel. Inhibitors of KCa3.1 channel or Kcnn4 KO reduced [Ca2+]ER, attenuated increase of [Ca2+]i or inflammasome activation by LPS + PA, and ameliorated HFD-induced inflammasome or metabolic inflammation. Lysosomal Ca2+ release induced delayed JNK and ASC phosphorylation through CAMKII-ASK1. These results suggest a novel role of lysosomal Ca2+ release sustained by ER→lysosome Ca2+ refilling and K+ efflux through KCa3.1 channel in inflammasome activation and metabolic inflammation.

Insight on the interaction between the scorpion toxin blocker Discrepin on potassium voltage-gated channel Kv4.3 by molecular dynamics simulations

Discrepin is a 38-residue α-toxin extracted from the venom of the Venezuelan scorpion Tityus discrepans, which inhibits ionic transit in the voltage-dependent potassium channels (Kv) of A-type current. The effect of specific residues on the IC₅₀ between Discrepine and Kv4.3, the main component of A-type currents, is known; however, the molecular details of the toxin-channel interaction are not known. In this work, we present interaction models between Discrepin (wt) and two peptide variants (V6K/D20K and K13A) on the pore-forming domain of the Kv4.3 channel obtained from homology, docking, and molecular dynamics modeling techniques. The free energy calculations in these models correspond to the order of the experimentally determined IC₅₀ values. Our studies shed light on the role of the K13 residue as responsible for occluding the Kv4.3 selectivity filter and the importance of the V6K mutation in the approach and stabilization of toxin-channel complex interactions.Communicated by Ramaswamy H. Sarma.

Retroviral glycoprotein-mediated immune suppression via the potassium channel KCa3.1 - A new strategy for amelioration of inflammatory bowel diseases

Peptides derived from retroviral envelope proteins have been shown to possess a wide range of immunosuppressive and anti-inflammatory activities. We have previously reported identification of such a peptide derived from the envelope protein coded by a human endogenous retrovirus (HERV). In this study, we identified that in vitro the peptide inhibits the K_Ca3.1 potassium channel, a potential target for therapy of immune diseases. We describe in vitro ENV59-GP3 effects with respect to potency of inhibition on K_Ca3.1 channels and calcium influx. Furthermore, we asses in vivo the effect of blocking K_Ca3.1 with ENV59-GP3 peptide or K_Ca3.1-blocker NS6180 on protection against DSS-induced acute colitis. ENV59-GP3 peptide treatment showed reduction of the disease score in the DSS-induced acute colitis mice model, which was comparable to effects of the K_Ca3.1 channel blocker NS6180. Analysis of cytokine production from DSS-mice model treated animals revealed equipotent inhibitory effects of the ENV59-GP3 and NS6180 compounds on the production of IL-6, TNF-α, IL-1β. These findings altogether suggest that ENV59-GP3 functions as a K_Ca3.1 channel inhibitor and underline the implications of using virus derived channel blockers for treatment of autoimmune diseases. Additionally, they open the possibilities whether K_Ca3.1 inhibition is efficacious in patients with inflammatory bowel diseases.

Diverse Structural Features of Potassium Channels Characterized by Scorpion Toxins as Molecular Probes

Scorpion toxins are well-known as the largest potassium channel peptide blocker family. They have been successfully proven to be valuable molecular probes for structural research on diverse potassium channels. The potassium channel pore region, including the turret and filter regions, is the binding interface for scorpion toxins, and structural features from different potassium channels have been identified using different scorpion toxins. According to the spatial orientation of channel turrets with differential sequence lengths and identities, conformational changes and molecular surface properties, the potassium channel turrets can be divided into the following three states: open state with less hindering effects on toxin binding, half-open state or half-closed state with certain effects on toxin binding, and closed state with remarkable effects on toxin binding. In this review, we summarized the diverse structural features of potassium channels explored using scorpion toxin tools and discuss future work in the field of scorpion toxin-potassium channel interactions.

A Novel ShK-Like Toxic Peptide from the Transcriptome of the Cnidarian Palythoa caribaeorum Displays Neuroprotection and Cardioprotection in Zebrafish

Palythoa caribaeorum (class Anthozoa) is a zoantharian which, together with other cnidarians, like jellyfishes, hydra, and sea anemones, possesses specialized structures in its tissues, the cnidocytes, which deliver an array of toxins in order to capture prey and deter predators. The whole transcriptome of P. caribaeroum was deep sequenced, and a diversity of toxin-related peptide sequences were identified, and some retrieved for functional analysis. In this work, a peptide precursor containing a ShK domain, named PcShK3, was analyzed by means of computational processing, comprising structural phylogenetic analysis, model prediction, and dynamics simulation of peptide-receptor interaction. The combined data indicated that PcShK3 is a distinct peptide which is homologous to a cluster of peptides belonging to the ShK toxin family. In vivo, PcShK3 distributed across the vitelline membrane and accumulated in the yolk sac stripe of zebrafish larvae. Notably, it displayed a significant cardio-protective effect in zebrafish in concentrations inferior to the IC₅₀ (<43.53 ± 6.45 &micro;M), while in high concentrations (>IC₅₀), it accumulated in the blood and caused pericardial edema, being cardiotoxic to zebrafish larvae. Remarkably, PcShK3 suppressed the 6-OHDA-induced neurotoxicity on the locomotive behavior of zebrafish. The present results indicated that PcShK3 is a novel member of ShK toxin family, and has the intrinsic ability to induce neuro- and cardio-protective effects or cause cardiac toxicity, according to its effective concentration.

Development of constant-pH simulation methods in implicit solvent and applications in biomolecular systems

pH is a critical parameter for biological and technological systems directly related with electrical charges. It can give rise to peculiar electrostatic phenomena, which also makes them more challenging. Due to the quantum nature of the process, involving the forming and breaking of chemical bonds, quantum methods should ideally by employed. Nevertheless, due to the very large number of ionizable sites, different macromolecular conformations, salt conditions, and all other charged species, the CPU time cost simply becomes prohibitive for computer simulations, making this a quite complex problem. Simplified methods based on Monte Carlo sampling have been devised and will be reviewed here, highlighting the updated state-of-the-art of this field, advantages, and limitations of different theoretical protocols for biomolecular systems (proteins and nucleic acids). Following a historical perspective, the discussion will be associated with the applications to protein interactions with other proteins, polyelectrolytes, and nanoparticles.

Molecular Simulations of Disulfide-Rich Venom Peptides with Ion Channels and Membranes

Disulfide-rich peptides isolated from the venom of arthropods and marine animals are a rich source of potent and selective modulators of ion channels. This makes these peptides valuable lead molecules for the development of new drugs to treat neurological disorders. Consequently, much effort goes into understanding their mechanism of action. This paper presents an overview of how molecular simulations have been used to study the interactions of disulfide-rich venom peptides with ion channels and membranes. The review is focused on the use of docking, molecular dynamics simulations, and free energy calculations to (i) predict the structure of peptide-channel complexes; (ii) calculate binding free energies including the effect of peptide modifications; and (iii) study the membrane-binding properties of disulfide-rich venom peptides. The review concludes with a summary and outlook.

C-Terminal residues in small potassium channel blockers OdK1 and OSK3 from scorpion venom fine-tune the selectivity

We report isolation, sequencing, and electrophysiological characterization of OSK3 (α-KTx 8.8 in Kalium and Uniprot databases), a potassium channel blocker from the scorpion Orthochirus scrobiculosus venom. Using the voltage clamp technique, OSK3 was tested on a wide panel of 11 voltage-gated potassium channels expressed in Xenopus oocytes, and was found to potently inhibit Kv1.2 and Kv1.3 with IC₅₀ values of ~331nM and ~503nM, respectively. OdK1 produced by the scorpion Odontobuthus doriae differs by just two C-terminal residues from OSK3, but shows marked preference to Kv1.2. Based on the charybdotoxin-potassium channel complex crystal structure, a model was built to explain the role of the variable residues in OdK1 and OSK3 selectivity.

Modeling of the Binding of Peptide Blockers to Voltage-Gated Potassium Channels: Approaches and Evidence

Modeling of the structure of voltage-gated potassium (KV) channels bound to peptide blockers aims to identify the key amino acid residues dictating affinity and provide insights into the toxin-channel interface. Computational approaches open up possibilities for in silico rational design of selective blockers, new molecular tools to study the cellular distribution and functional roles of potassium channels. It is anticipated that optimized blockers will advance the development of drugs that reduce over activation of potassium channels and attenuate the associated malfunction. Starting with an overview of the recent advances in computational simulation strategies to predict the bound state orientations of peptide pore blockers relative to KV-channels, we go on to review algorithms for the analysis of intermolecular interactions, and then take a look at the results of their application.

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.

Scorpion toxins prefer salt solutions

There is a wide variety of ion channel types with various types of blockers, making research in this field very complicated. To reduce this complexity, it is essential to study ion channels and their blockers independently. Scorpion toxins, a major class of blockers, are charged short peptides with high affinities for potassium channels. Their high selectivity and inhibitory properties make them an important pharmacological tool for treating autoimmune or nervous system disorders. Scorpion toxins typically have highly charged surfaces and-like other proteins-an intrinsic ability to bind ions (Friedman J Phys Chem B 115(29):9213-9223, 1996; Baldwin Biophys J 71(4):2056-2063, 1996; Vrbka et al. Proc Natl Acad Sci USA 103(42):15440-15444, 2006a; Vrbka et al. J Phys Chem B 110(13):7036-43, 2006b). Thus, their effects on potassium channels are usually investigated in various ionic solutions. In this work, computer simulations of protein structures were performed to analyze the structural properties of the key residues (i.e., those that are presumably involved in contact with the surfaces of the ion channels) of 12 scorpion toxins. The presence of the two most physiologically abundant cations, Na(+) and K(+), was considered. The results indicated that the ion-binding properties of the toxin residues vary. Overall, all of the investigated toxins had more stable structures in ionic solutions than in water. We found that both the number and length of elements in the secondary structure varied depending on the ionic solution used (i.e., in the presence of NaCl or KCl). This study revealed that the ionic solution should be chosen carefully before performing experiments on these toxins. Similarly, the influence of these ions should be taken into consideration in the design of toxin-based pharmaceuticals.

Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data

In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca(2+)-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the "open" and "closed" conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM's inhibitors correlated well with available experimental data as the r(2) obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca(2+)-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca(2+)-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca(2+)-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.

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

Ion channels represent a large family of membrane proteins with many being well established targets in pharmacotherapy. The ‘druggability’ of heteromeric channels comprised of different subunits remains obscure, due largely to a lack of channel-specific probes necessary to delineate their therapeutic potential in vivo. Our initial studies reported here, investigated the family of inwardly rectifying potassium (Kir) channels given the availability of high resolution crystal structures for the eukaryotic constitutively active Kir2.2 channel. We describe a ‘limited’ homology modeling approach that can yield chimeric Kir channels having an outer vestibule structure representing nearly any known vertebrate or invertebrate channel. These computationally-derived channel structures were tested ""in silico for ‘docking’ to NMR structures of tertiapin (TPN), a 21 amino acid peptide found in bee venom. TPN is a highly selective and potent blocker for the epithelial rat Kir1.1 channel, but does not block human or zebrafish Kir1.1 channel isoforms. Our Kir1.1 channel-TPN docking experiments recapitulated published in vitro ""findings for TPN-sensitive and TPN-insensitive channels. Additionally, in silico site-directed mutagenesis identified ‘hot spots’ within the channel outer vestibule that mediate energetically favorable docking scores and correlate with sites previously identified with in vitro thermodynamic mutant-cycle analysis. These ‘proof-of-principle’ results establish a framework for virtual screening of re-engineered peptide toxins for interactions with computationally derived Kir channels that currently lack channel-specific blockers. When coupled with electrophysiological validation, this virtual screening approach may accelerate the drug discovery process, and can be readily applied to other ion channels families where high resolution structures are available.

Binding modes of two scorpion toxins to the voltage-gated potassium channel kv1.3 revealed from molecular dynamics

Molecular dynamics (MD) simulations are used to examine the binding modes of two scorpion toxins, margatoxin (MgTx) and hongotoxin (HgTx), to the voltage gated K+ channel, Kv1.3. Using steered MD simulations, we insert either Lys28 or Lys35 of the toxins into the selectivity filter of the channel. The MgTx-Kv1.3 complex is stable when the side chain of Lys35 from the toxin occludes the channel filter, suggesting that Lys35 is the pore-blocking residue for Kv1.3. In this complex, Lys28 of the toxin forms one additional salt bridge with Asp449 just outside the filter of the channel. On the other hand, HgTx forms a stable complex with Kv1.3 when the side chain of Lys28 but not Lys35 protrudes into the filter of the channel. A survey of all the possible favorable binding modes of HgTx-Kv1.3 is carried out by rotating the toxin at 3° intervals around the channel axis while the position of HgTx-Lys28 relative to the filter is maintained. We identify two possible favorable binding modes: HgTx-Arg24 can interact with either Asp433 or Glu420 on the vestibular wall of the channel. The dissociation constants calculated from the two binding modes of HgTx-Kv1.3 differ by approximately 20 fold, suggesting that the two modes are of similar energetics.

Theoretical frameworks for multiscale modeling and simulation

Biomolecular systems have been modeled at a variety of scales, ranging from explicit treatment of electrons and nuclei to continuum description of bulk deformation or velocity. Many challenges of interfacing between scales have been overcome. Multiple models at different scales have been used to study the same system or calculate the same property (e.g., channel conductance). Accurate modeling of biochemical processes under in vivo conditions and the bridging of molecular and subcellular scales will likely soon become reality.