Modeling the structure of agitoxin in complex with the Shaker K+ channel: a computational approach based on experimental distance restraints extracted from thermodynamic mutant cycles

Deciphering the shape selective conformational equilibrium of E- and Z-locked azobenzene-tetraethylammonium ion in regulating photo-switchable K+-ion channel blocking

The search for photo-switchable optopharmacological agents that can block ion channels has been a prevalent area owing to its prime advantages of reversibility and specificity over the traditional blockers. However, the quest for a higher blocking ability shown by a less stable photo-isomer to perfectly suit the requirement of the optopharmacological agents is still ongoing. To date, only a marginal improvement in terms of blocking ability is observed by the less stable E-isomer of para-substituted locked azobenzene with TEA (LAB-TEA) for the K+-ion channel. Thus, rationalization of the limitation for achieving high activity by the E-isomer is rather essential to aid the improvement of the efficiency of photoswitchable blocker drugs. Herein, we report a molecular-level analysis on the mechanism of blocking by E- and Z-LAB-TEA with the bacterial KcsA K+-ion channel using Molecular Dynamics (MD) simulation and Quantum Mechanical (QM) calculations. The positively charged TEA fragment engages in stronger electrostatic interactions, while the neutral LAB fragment engages in weaker dispersive interactions. The binding free energy calculated by Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) for E-LAB-TEA (-22.3 kcal mol-1) shows less thermodynamic preference for binding with K+-ion channels than Z-LAB-TEA (-21.6 kcal mol-1) corroborating the experimental observation. The correlation between the structure and the binding ability of E- and Z-isomers of LAB-TEA indicates that the channel gate is narrow and acts as a bottleneck for the entry of the binder molecule inside the large cavity. Upon irradiation, the Z-isomer converts into a less stable but long and planar E-isomer (ΔE of photoisomerism = 7.0 kcal mol-1, at SA2-CASPT2(6,4)/6-31+G(d)//CASSCF(6,4)/6-31+G(d)), which is structurally more suitable to fit into the narrow channel gate rather than the curved and non-planar Z-LAB-TEA. Thus, a reduction in the ionic current is observed owing to the preferential entry and subsequent blocking by E-LAB-TEA. Discontinuing the irradiation leads to conversion to the Z-isomer, the curved nature of which hinders its spontaneous release outside the cavity, thereby contributing only a small increase in the ionic current.

SKCa- and Kv1-type potassium channels and cancer: Promising therapeutic targets?

Ion channels are transmembrane structures that allow the passage of ions across cell membranes such as the plasma membrane or the membranes of various organelles like the nucleus, endoplasmic reticulum, Golgi apparatus or mitochondria. Aberrant expression of various ion channels has been demonstrated in several tumor cells, leading to the promotion of key functions in tumor development, such as cell proliferation, resistance to apoptosis, angiogenesis, invasion and metastasis. The link between ion channels and these key biological functions that promote tumor development has led to the classification of cancers as oncochannelopathies. Among all ion channels, the most varied and numerous, forming the largest family, are the potassium channels, with over 70 genes encoding them in humans. In this context, this review will provide a non-exhaustive overview of the role of plasma membrane potassium channels in cancer, describing 1) the nomenclature and structure of potassium channels, 2) the role of these channels in the control of biological functions that promotes tumor development such as proliferation, migration and cell death, and 3) the role of two particular classes of potassium channels, the SKCa- and Kv1- type potassium channels in cancer progression.

Computational Assessment of Protein-Protein Binding Specificity within a Family of Synaptic Surface Receptors

Atomic-level information is essential to explain the formation of specific protein complexes in terms of structure and dynamics. The set of Dpr and DIP proteins, which play a key role in the neuromorphogenesis in the nervous system of Drosophila melanogaster, offer a rich paradigm to learn about protein-protein recognition. Many members of the DIP subfamily cross-react with several members of the Dpr family and vice versa. While there exists a total of 231 possible Dpr-DIP heterodimer complexes from the 21 Dpr and 11 DIP proteins, only 57 "cognate" pairs have been detected by surface plasmon resonance (SPR) experiments, suggesting that the remaining 174 pairs have low or unreliable binding affinity. Our goal is to assess the performance of computational approaches to characterize the global set of interactions between Dpr and DIP proteins and identify the specificity of binding between each DIP with their corresponding Dpr binding partners. In addition, we aim to characterize how mutations influence the specificity of the binding interaction. In this work, a wide range of knowledge-based and physics-based approaches are utilized, including mutual information, linear discriminant analysis, homology modeling, molecular dynamics simulations, Poisson-Boltzmann continuum electrostatics calculations, and alchemical free energy perturbation to decipher the origin of binding specificity of the Dpr-DIP complexes examined. Ultimately, the results show that those two broad strategies are complementary, with different strengths and limitations. Biological inter-relations are more clearly revealed through knowledge-based approaches combining evolutionary and structural features, the molecular determinants controlling binding specificity can be predicted accurately with physics-based approaches based on atomic models.

Comparison of K+ Channel Families

K⁺ channels enable potassium to flow across the membrane with great selectivity. There are four K⁺ channel families: voltage-gated K (K_v), calcium-activated (K_Ca), inwardly rectifying K (K_ir), and two-pore domain potassium (K_2P) channels. All four K⁺ channels are formed by subunits assembling into a classic tetrameric (4x1P = 4P for the K_v, K_Ca, and K_ir channels) or tetramer-like (2x2P = 4P for the K_2P channels) architecture. These subunits can either be the same (homomers) or different (heteromers), conferring great diversity to these channels. They share a highly conserved selectivity filter within the pore but show different gating mechanisms adapted for their function. K⁺ channels play essential roles in controlling neuronal excitability by shaping action potentials, influencing the resting membrane potential, and responding to diverse physicochemical stimuli, such as a voltage change (K_v), intracellular calcium oscillations (K_Ca), cellular mediators (K_ir), or temperature (K_2P).

Dual-Site Binding of Quaternary Ammonium Ions as Internal K+-Ion Channel Blockers: Nonclassical (C-H···O) H Bonding vs Dispersive (C-H···H-C) Interaction

A molecular-level study of the influence of the alkyl chain length of quaternary ammonium ions (QAs) on the blocking action and the mode of binding with the bacterial KcsA K⁺-ion channel is carried out by molecular dynamics (MD) simulations as well as quantum mechanics/molecular mechanics (QM/MM) methods. The present work unveils distinct modes of binding for different QAs, due to differences in size and hydrophobicity. The QAs bind near the channel gate as well as at the central cavity, leading to a possible dual-site blocking action. Small-sized tetraethylammonium (TEA) and tetrabutylammonium (TBA) ions enter inside the channel cavity in the open state of KcsA but bind strongly in the closed state. TEA binds to the polar hydroxyl group of threonine residues situated at the channel gate via nonclassical H-bonding interaction (C-H···O), while TBA binds to a second binding site, the central cavity, with hydrophobic benzyl and sec-butyl side chains of phenylalanine and isoleucine residues via alkyl-π and hydrophobic interactions (C-H···H-C). On the contrary, large tetrahexylammonium (THA) and tetraoctylammonium (TOA) ions bind the hydrophobic side-chain methyl and isopropyl of alanine and valine at the channel gate both in the open and closed states, thereby restricting the free movement of large QAs toward the center of the cavity. However, the binding to the hydrophobic benzyl and sec-butyl side chains of phenylalanine and isoleucine residues in the closed state is thermodynamically preferable. Also, the binding energy is found to increase with an increase in the alkyl chain length from ethyl (-16.4 kcal/mol) to octyl (-65.5 kcal/mol), due to an almost linear increase in dispersive interaction.

Hydrophobic Drug/Toxin Binding Sites in Voltage-Dependent K+ and Na+ Channels

In the Na_v channel family the lipophilic drugs/toxins binding sites and the presence of fenestrations in the channel pore wall are well defined and categorized. No such classification exists in the much larger K_v channel family, although certain lipophilic compounds seem to deviate from binding to well-known hydrophilic binding sites. By mapping different compound binding sites onto 3D structures of Kv channels, there appear to be three distinct lipid-exposed binding sites preserved in K_v channels: the front and back side of the pore domain, and S2-S3/S3-S4 clefts. One or a combination of these sites is most likely the orthologous equivalent of neurotoxin site 5 in Na_v channels. This review describes the different lipophilic binding sites and location of pore wall fenestrations within the K_v channel family and compares it to the knowledge of Na_v channels.

The binding process of BmKTX and BmKTX-D33H toward to Kv1.3 channel: a molecular dynamics simulation study

The potassium channel Kv1.3 is an important pharmacological target and the Kaliotoxin-type toxins (α-KTX-3 family) are its specific blockers. Here, we study the binding process of two kinds of Kaliotoxin-type toxins:BmKTX and its mutant (BmKTX-D33H) toward to Kv1.3 channel using MD simulation and umbrella sampling simulation, respectively. The calculated binding free energies are -27 kcal/mol and -34 kcal/mol for BmKTX and BmKTX-D33H, respectively, which are consistent with experimental results. The further analysis indicate that the characteristic of electrostatic potential of the α-KTX-3 have important effect on their binding modes with Kv1.3 channel; the residue 33 in BmKTX or BmKTX-D33H plays a key role in determine their binding orientations toward to Kv1.3 channel; when residue 33 (or 34) has negative electrostatic potential, the anti-parallel β-sheet domain of α-KTX-3 toxin peptide will keep away from the filter region of Kv1.3 channel, as BmKTX; when residue 33(or 34) has positive electrostatic potential, the anti-parallel β-sheet domain of α-KTX-3 toxin peptide will interact with the filter region of Kv1.3 channel, as BmKTX-D33H. Above all, electrostatic potential differences on toxin surfaces and correlations motions within the toxins will determine the toxin-potassium channel interaction model. In addition, the hydrogen bond interaction is the pivotal factor for the Kv1.3-Kaliotoxin association. Understanding the binding mechanism of toxin-potassium channel will facilitate the rational development of new toxin analogue.Communicated by Ramaswamy H. Sarma.

Marine Toxins Targeting Kv1 Channels: Pharmacological Tools and Therapeutic Scaffolds

Toxins from marine animals provide molecular tools for the study of many ion channels, including mammalian voltage-gated potassium channels of the Kv1 family. Selectivity profiling and molecular investigation of these toxins have contributed to the development of novel drug leads with therapeutic potential for the treatment of ion channel-related diseases or channelopathies. Here, we review specific peptide and small-molecule marine toxins modulating Kv1 channels and thus cover recent findings of bioactives found in the venoms of marine Gastropod (cone snails), Cnidarian (sea anemones), and small compounds from cyanobacteria. Furthermore, we discuss pivotal advancements at exploiting the interaction of κM-conotoxin RIIIJ and heteromeric Kv1.1/1.2 channels as prevalent neuronal Kv complex. RIIIJ's exquisite Kv1 subtype selectivity underpins a novel and facile functional classification of large-diameter dorsal root ganglion neurons. The vast potential of marine toxins warrants further collaborative efforts and high-throughput approaches aimed at the discovery and profiling of Kv1-targeted bioactives, which will greatly accelerate the development of a thorough molecular toolbox and much-needed therapeutics.

String Method for Protein-Protein Binding Free-Energy Calculations

A powerful computational strategy to determine the equilibrium association constant of two macromolecules with explicit-solvent molecular dynamics (MD) simulations is the "geometric route", which considers the reversible physical separation of the bound complex in solution. Nonetheless, multiple challenges remain to render this type of methodology reliable and computationally efficient in practice. In particular, in one, formulation of the geometric route relies on the potential of mean force (PMF) for physically separating the two binding partners restrained along a straight axis, which must be selected prior to the calculation. However, practical applications indicate that the calculation of the separation PMF along the predefined rectilinear pathway may be suboptimal and slowly convergent. Recognizing that a rectilinear straight separation pathway is generally not representative of how the protein complex physically separates in solution, we put forth a novel theoretical framework for binding free-energy calculations, leaning on the optimal curvilinear minimum free-energy path (MFEP) determined from the string method. The proposed formalism is validated by comparing the results obtained using both rectilinear and curvilinear pathways for a prototypical host-guest complex formed by cucurbit[7]uril (CB[7]) binding benzene, and for the barnase-barstar protein complex. On the basis of multi-microsecond MD calculations, we find that the calculations following the traditional rectilinear pathway and the string-based curvilinear pathway agree quantitatively, but convergence is faster with the latter.

Trans-toxin ion-sensitivity of charybdotoxin-blocked potassium-channels reveals unbinding transitional states

In silico and in vitro studies have made progress in understanding protein–protein complex formation; however, the molecular mechanisms for their dissociation are unclear. Protein–protein complexes, lasting from microseconds to years, often involve induced-fit, challenging computational or kinetic analysis. Charybdotoxin (CTX), a peptide from the Leiurus scorpion venom, blocks voltage-gated K⁺-channels in a unique example of binding/unbinding simplicity. CTX plugs the external mouth of K⁺-channels pore, stopping K⁺-ion conduction, without inducing conformational changes. Conflicting with a tight binding, we show that external permeant ions enhance CTX-dissociation, implying a path connecting the pore, in the toxin-bound channel, with the external solution. This sensitivity is explained if CTX wobbles between several bound conformations, producing transient events that restore the electrical and ionic trans-pore gradients. Wobbling may originate from a network of contacts in the interaction interface that are in dynamic stochastic equilibria. These partially-bound intermediates could lead to distinct, and potentially manipulable, dissociation pathways.

High-speed AFM reveals accelerated binding of agitoxin-2 to a K+ channel by induced fit

Agitoxin-2 (AgTx2) from scorpion venom is a potent blocker of K⁺ channels. The docking model has been elucidated, but it remains unclear whether binding dynamics are described by a two-state model (AgTx2-bound and AgTx2-unbound) or a more complicated mechanism, such as induced fit or conformational selection. Here, we observed the binding dynamics of AgTx2 to the KcsA channel using high-speed atomic force microscopy. From images of repeated binding and dissociation of AgTx2 to the channel, single-molecule kinetic analyses revealed that the affinity of the channel for AgTx2 increased during persistent binding and decreased during persistent dissociation. We propose a four-state model, including high- and low-affinity states of the channel, with relevant rate constants. An induced-fit pathway was dominant and accelerated binding by 400 times. This is the first analytical imaging of scorpion toxin binding in real time, which is applicable to various biological dynamics including channel ligands, DNA-modifier proteins, and antigen-antibody complexes.

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.

Combining the polarizable Drude force field with a continuum electrostatic Poisson-Boltzmann implicit solvation model

In this work, we have combined the polarizable force field based on the classical Drude oscillator with a continuum Poisson-Boltzmann/solvent-accessible surface area (PB/SASA) model. In practice, the positions of the Drude particles experiencing the solvent reaction field arising from the fixed charges and induced polarization of the solute must be optimized in a self-consistent manner. Here, we parameterized the model to reproduce experimental solvation free energies of a set of small molecules. The model reproduces well-experimental solvation free energies of 70 molecules, yielding a root mean square difference of 0.8 kcal/mol versus 2.5 kcal/mol for the CHARMM36 additive force field. The polarization work associated with the solute transfer from the gas-phase to the polar solvent, a term neglected in the framework of additive force fields, was found to make a large contribution to the total solvation free energy, comparable to the polar solute-solvent solvation contribution. The Drude PB/SASA also reproduces well the electronic polarization from the explicit solvent simulations of a small protein, BPTI. Model validation was based on comparisons with the experimental relative binding free energies of 371 single alanine mutations. With the Drude PB/SASA model the root mean square deviation between the predicted and experimental relative binding free energies is 3.35 kcal/mol, lower than 5.11 kcal/mol computed with the CHARMM36 additive force field. Overall, the results indicate that the main limitation of the Drude PB/SASA model is the inability of the SASA term to accurately capture non-polar solvation effects. © 2018 Wiley Periodicals, Inc.

Molecular Dynamics Simulation Reveals Specific Interaction Sites between Scorpion Toxins and Kv1.2 Channel: Implications for Design of Highly Selective Drugs

The K_v1.2 channel plays an important role in the maintenance of resting membrane potential and the regulation of the cellular excitability of neurons, whose silencing or mutations can elicit neuropathic pain or neurological diseases (e.g., epilepsy and ataxia). Scorpion venom contains a variety of peptide toxins targeting the pore region of this channel. Despite a large amount of structural and functional data currently available, their detailed interaction modes are poorly understood. In this work, we choose four K_v1.2-targeted scorpion toxins (Margatoxin, Agitoxin-2, OsK-1, and Mesomartoxin) to construct their complexes with K_v1.2 based on the experimental structure of ChTx-K_v1.2. Molecular dynamics simulation of these complexes lead to the identification of hydrophobic patches, hydrogen-bonds, and salt bridges as three essential forces mediating the interactions between this channel and the toxins, in which four K_v1.2-specific interacting amino acids (D353, Q358, V381, and T383) are identified for the first time. This discovery might help design highly selective K_v1.2-channel inhibitors by altering amino acids of these toxins binding to the four channel residues. Finally, our results provide new evidence in favor of an induced fit model between scorpion toxins and K⁺ channel interactions.

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.

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.

Complexes of Peptide Blockers with Kv1.6 Pore Domain: Molecular Modeling and Studies with KcsA-Kv1.6 Channel

Potassium voltage-gated Kv1.6 channel, which is distributed primarily in neurons of central and peripheral nervous systems, is of significant physiological importance. To date, several high-affinity Kv1.6-channel blockers are known, but the lack of selective ones among them hampers the studies of tissue localization and functioning of Kv1.6 channels. Here we present an approach to advanced understanding of interactions of peptide toxin blockers with a Kv1.6 pore. It combines molecular modeling studies and an application of a new bioengineering system based on a KcsA-Kv1.6 hybrid channel for the quantitative fluorescent analysis of blocker-channel interactions. Using this system we demonstrate that peptide toxins agitoxin 2, kaliotoxin1 and OSK1 have similar high affinity to the extracellular vestibule of the K⁺-conducting pore of Kv1.6, hetlaxin is a low-affinity ligand, whereas margatoxin and scyllatoxin do not bind to Kv1.6 pore. Binding of toxins to Kv1.6 pore has considerable inverse dependence on the ionic strength. Model structures of KcsA-Kv1.6 and Kv1.6 complexes with agitoxin 2, kaliotoxin 1 and OSK1 were obtained using homology modeling and molecular dynamics simulation. Interaction interfaces, which are formed by 15-19 toxin residues and 10 channel residues, are described and compared. Specific sites of Kv1.6 pore recognition are identified for targeting of peptide blockers. Analysis of interactions between agitoxin 2 derivatives with point mutations (S7K, S11G, L19S, R31G) and KcsA-Kv1.6 confirms reliability of the calculated complex structure.

Structural correlates of affinity in fetal versus adult endplate nicotinic receptors

Adult-type nicotinic acetylcholine receptors (AChRs) mediate signalling at mature neuromuscular junctions and fetal-type AChRs are necessary for proper synapse development. Each AChR has two neurotransmitter binding sites located at the interface of a principal and a complementary subunit. Although all agonist binding sites have the same core of five aromatic amino acids, the fetal site has ∼30-fold higher affinity for the neurotransmitter ACh. Here we use molecular dynamics simulations of adult versus fetal homology models to identify complementary-subunit residues near the core that influence affinity, and use single-channel electrophysiology to corroborate the results. Four residues in combination determine adult versus fetal affinity. Simulations suggest that at lower-affinity sites, one of these unsettles the core directly and the others (in loop E) increase backbone flexibility to unlock a key, complementary tryptophan from the core. Swapping only four amino acids is necessary and sufficient to exchange function between adult and fetal AChRs.

Adult and fetal nicotinic acetylcholine receptors (AChRs) have different functional requirements and affinity for ACh. Here, the authors use molecular dynamics and electrophysiology to investigate this affinity, and identify four amino acids that when swapped exchange function between adult and fetal AChRs.

The Scorpion Toxin Analogue BmKTX-D33H as a Potential Kv1.3 Channel-Selective Immunomodulator for Autoimmune Diseases

The Kv1.3 channel-acting scorpion toxins usually adopt the conserved anti-parallel β-sheet domain as the binding interface, but it remains challenging to discover some highly selective Kv1.3 channel-acting toxins. In this work, we investigated the pharmacological profile of the Kv1.3 channel-acting BmKTX-D33H, a structural analogue of the BmKTX scorpion toxin. Interestingly, BmKTX-D33H, with its conserved anti-parallel β-sheet domain as a Kv1.3 channel-interacting interface, exhibited more than 1000-fold selectivity towards the Kv1.3 channel as compared to other K⁺ channels (including Kv1.1, Kv1.2, Kv1.7, Kv11.1, KCa2.2, KCa2.3, and KCa3.1). As expected, BmKTX-D33H was found to inhibit the cytokine production and proliferation of both Jurkat cells and human T cells in vitro. It also significantly improved the delayed-type hypersensitivity (DTH) responses, an autoreactive T cell-mediated inflammation in rats. Amino acid sequence alignment and structural analysis strongly suggest that the “evolutionary” Gly11 residue of BmKTX-D33H interacts with the turret domain of Kv1 channels; it appears to be a pivotal amino acid residue with regard to the selectivity of BmKTX-D33H towards the Kv1.3 channel (in comparison with the highly homologous scorpion toxins). Together, our data indicate that BmKTX-D33H is a Kv1.3 channel–specific blocker. Finally, the remarkable selectivity of BmKTX-D33H highlights the great potential of evolutionary-guided peptide drug design in future studies.

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.

Interaction of peptides with cell membranes: insights from molecular modeling

The investigation of the interaction of peptides with cell membranes is the focus of active research. It can enhance the understanding of basic membrane functions such as membrane transport, fusion, and signaling processes, and it may shed light on potential applications of peptides in biomedicine. In this review, we will present current advances in computational studies on the interaction of different types of peptides with the cell membrane. Depending on the properties of the peptide, membrane, and external environment, the peptide-membrane interaction shows a variety of different forms. Here, on the basis of recent computational progress, we will discuss how different peptides could initiate membrane pores, translocate across the membrane, induce membrane endocytosis, produce membrane curvature, form fibrils on the membrane surface, as well as interact with functional membrane proteins. Finally, we will present a conclusion summarizing recent progress and providing some specific insights into future developments in this field.

Combination of Ambiguous and Unambiguous Data in the Restraint-driven Docking of Flexible Peptides with HADDOCK: The Binding of the Spider Toxin PcTx1 to the Acid Sensing Ion Channel (ASIC) 1a

Peptides that bind to ion channels have attracted much interest as potential lead molecules for the development of new drugs and insecticides. However, the structure determination of large peptide-channel complexes using experimental methods is challenging. Thus structural models are often derived from combining experimental information with restraint-driven docking approaches. Using the complex formed by the venom peptide PcTx1 and the acid sensing ion channel (ASIC) 1a as a case study, we have examined the effect of different combinations of restraints and input structures on the statistical likelihood of (a) correctly predicting the structure of the binding interface and (b) the ability to predict which residues are involved in specific pairwise peptide-channel interactions. For this, we have analyzed over 200,000 water-refined docked structures obtained with various amounts and types of restraints of the peptide-channel complex predicted using the docking program HADDOCK. We found that increasing the number of restraints or even the use of pairwise interaction data resulted in only a modest improvement in the likelihood of finding a structure within a given accuracy. This suggests that shape complementarity and the force field make a large contribution to the accuracy of the predicted structure. The results also showed that there are large variations in the accuracy of the predicted structure depending on the precise combination of residues used as restraints. Finally, we reflect on the limitations of relying on geometric criteria such as root-mean square deviations to assess the accuracy of docking procedures. We propose that in addition to currently used measures, the likelihood of finding a structure within a given level of accuracy should be also used to evaluate docking methods.

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.

Structural basis for the inhibition of voltage-dependent K+ channel by gating modifier toxin

Voltage-dependent K⁺ (K_v) channels play crucial roles in nerve and muscle action potentials. Voltage-sensing domains (VSDs) of K_v channels sense changes in the transmembrane potential, regulating the K⁺-permeability across the membrane. Gating modifier toxins, which have been used for the functional analyses of K_v channels, inhibit K_v channels by binding to VSD. However, the structural basis for the inhibition remains elusive. Here, fluorescence and NMR analyses of the interaction between VSD derived from K_vAP channel and its gating modifier toxin, VSTx1, indicate that VSTx1 recognizes VSD under depolarized condition. We identified the VSD-binding residues of VSTx1 and their proximal residues of VSD by the cross-saturation (CS) and amino acid selective CS experiments, which enabled to build a docking model of the complex. These results provide structural basis for the specific binding and inhibition of K_v channels by gating modifier toxins.

Unraveling the mechanism of selective ion transport in hydrophobic subnanometer channels

Recently reported synthetic organic nanopore (SONP) can mimic a key feature of natural ion channels, i.e., selective ion transport. However, the physical mechanism underlying the K(+)/Na(+) selectivity for the SONPs is dramatically different from that of natural ion channels. To achieve a better understanding of the selective ion transport in hydrophobic subnanometer channels in general and SONPs in particular, we perform a series of ab initio molecular dynamics simulations to investigate the diffusivity of aqua Na(+) and K(+) ions in two prototype hydrophobic nanochannels: (i) an SONP with radius of 3.2 Å, and (ii) single-walled carbon nanotubes (CNTs) with radii of 3-5 Å (these radii are comparable to those of the biological potassium K(+) channels). We find that the hydration shell of aqua Na(+) ion is smaller than that of aqua K(+) ion but notably more structured and less yielding. The aqua ions do not lower the diffusivity of water molecules in CNTs, but in SONP the diffusivity of aqua ions (Na(+) in particular) is strongly suppressed due to the rugged inner surface. Moreover, the aqua Na(+) ion requires higher formation energy than aqua K(+) ion in the hydrophobic nanochannels. As such, we find that the ion (K(+) vs. Na(+)) selectivity of the (8, 8) CNT is ∼20× higher than that of SONP. Hence, the (8, 8) CNT is likely the most efficient artificial K(+) channel due in part to its special interior environment in which Na(+) can be fully solvated, whereas K(+) cannot. This work provides deeper insights into the physical chemistry behind selective ion transport in nanochannels.

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.

Decrypting the structural, dynamic, and energetic basis of a monomeric kinesin interacting with a tubulin dimer in three ATPase states by all-atom molecular dynamics simulation

We have employed molecular dynamics (MD) simulation to investigate, with atomic details, the structural dynamics and energetics of three major ATPase states (ADP, APO, and ATP state) of a human kinesin-1 monomer in complex with a tubulin dimer. Starting from a recently solved crystal structure of ATP-like kinesin-tubulin complex by the Knossow lab, we have used flexible fitting of cryo-electron-microscopy maps to construct new structural models of the kinesin-tubulin complex in APO and ATP state, and then conducted extensive MD simulations (total 400 ns for each state), followed by flexibility analysis, principal component analysis, hydrogen bond analysis, and binding free energy analysis. Our modeling and simulation have revealed key nucleotide-dependent changes in the structure and flexibility of the nucleotide-binding pocket (featuring a highly flexible and open switch I in APO state) and the tubulin-binding site, and allosterically coupled motions driving the APO to ATP transition. In addition, our binding free energy analysis has identified a set of key residues involved in kinesin-tubulin binding. On the basis of our simulation, we have attempted to address several outstanding issues in kinesin study, including the possible roles of β-sheet twist and neck linker docking in regulating nucleotide release and binding, the structural mechanism of ADP release, and possible extension and shortening of α4 helix during the ATPase cycle. This study has provided a comprehensive structural and dynamic picture of kinesin's major ATPase states, and offered promising targets for future mutational and functional studies to investigate the molecular mechanism of kinesin motors.

Normal mode dynamics of voltage-gated K(+) channels: gating principle, opening mechanism, and inhibition

Voltage-dependent potassium channels open in response to changes in membrane potential and become partially inactivated upon binding of inhibitors. Here we calculate normal mode motion of two voltage-dependent K(+) channels, KvAP and Shaker, and their complexes with inhibitors and address the gating principle, opening mechanism, and inhibition. The normal modes indicate that pore expansion and channel opening is correlated with a displacement of the arginine gating charges and a tilting of the voltage-sensor paddles. Normal modes of Shaker in complex with agitoxin, which blocks the central pore, do not display significantly altered paddle tilting and pore expansion. In contrast, normal modes of Shaker in complex with hanatoxin, which binds to the voltage sensor paddle, display decreased paddle tilting and pore expansion. This study presents a unified motion for the gating principle and channel opening, and offers insight into the voltage sensor paddle motion and its inhibition.

Effect of confinement on DNA, solvent and counterion dynamics in a model biological nanopore

The application of recent advances in nanopore technology to high-throughput DNA sequencing requires a more detailed understanding of solvent, ion and DNA interactions occurring within these pores. Here we present a combination of atomistic and coarse-grained modeling studies of the dynamics of short single-stranded DNA (ssDNA) homopolymers within the alpha-hemolysin pore, for the two single-stranded homopolymers poly(dA)40 and poly(dC)40. Analysis of atomistic simulations along with the per-residue decomposition of protein-DNA interactions in these simulations gives new insight into the very complex issues that have yet to be fully addressed with detailed MD simulations. We discuss a modification of the solvent properties and ion distribution around DNA within nanopore confinement and put it into the general framework of counterion condensation theory. There is a reasonable agreement in computed properties from our all-atom simulations and the resulting predictions from analytical theories with experimental data, and our equilibrium results here support the conclusions from our previous non-equilibrium Brownian dynamics studies with a recently developed BROMOC protocol that cations are the primary charge carriers through alpha-hemolysin nanopores under an applied voltage in the presence of ssDNA. Clustering analysis led to an identification of distinct conformational states of captured polymer and depth of the current blockade. Therefore, our data suggest that confined polymer may act as a flickering gate, thus contributing to excess noise phenomena. We also discuss the extent of water structuring due to nanopore confinement and the relationship between the conformational dynamics of a captured polymer and the distribution of blocked current.

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.

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

The Ca(2+)-activated channel of intermediate-conductance (KCa3.1) is a target for antisickling and immunosuppressant agents. Many small peptides isolated from animal venoms inhibit KCa3.1 with nanomolar affinities and are promising drug scaffolds. Although the inhibitory effect of peptide toxins on KCa3.1 has been examined extensively, the structural basis of toxin-channel recognition has not been understood in detail. Here, the binding modes of two selected scorpion toxins, charybdotoxin (ChTx) and OSK1, to human KCa3.1 are examined in atomic detail using molecular dynamics (MD) simulations. Employing a homology model of KCa3.1, we first determine conduction properties of the channel using Brownian dynamics and ascertain that the simulated results are in accord with experiment. The model structures of ChTx-KCa3.1 and OSK1-KCa3.1 complexes are then constructed using MD simulations biased with distance restraints. The ChTx-KCa3.1 complex predicted from biased MD is consistent with the crystal structure of ChTx bound to a voltage-gated K(+) channel. The dissociation constants (Kd) for the binding of both ChTx and OSK1 to KCa3.1 determined experimentally are reproduced within fivefold using potential of mean force calculations. Making use of the knowledge we gained by studying the ChTx-KCa3.1 complex, we attempt to enhance the binding affinity of the toxin by carrying out a theoretical mutagenesis. A mutant toxin, in which the positions of two amino acid residues are interchanged, exhibits a 35-fold lower Kd value for KCa3.1 than that of the wild-type. This study provides insight into the key molecular determinants for the high-affinity binding of peptide toxins to KCa3.1, and demonstrates the power of computational methods in the design of novel toxins.

Probing the flexibility of tropomyosin and its binding to filamentous actin using molecular dynamics simulations

Tropomyosin (Tm) is a coiled-coil protein that binds to filamentous actin (F-actin) and regulates its interactions with actin-binding proteins like myosin by moving between three positions on F-actin (the blocked, closed, and open positions). To elucidate the molecular details of Tm flexibility in relation to its binding to F-actin, we conducted extensive molecular dynamics simulations for both Tm alone and Tm-F-actin complex in the presence of explicit solvent (total simulation time >400 ns). Based on the simulations, we systematically analyzed the local flexibility of the Tm coiled coil using multiple parameters. We found a good correlation between the regions with high local flexibility and a number of destabilizing regions in Tm, including six clusters of core alanines. Despite the stabilization by F-actin binding, the distribution of local flexibility in Tm is largely unchanged in the absence and presence of F-actin. Our simulations showed variable fluctuations of individual Tm periods from the closed position toward the open position. In addition, we performed Tm-F-actin binding calculations based on the simulation trajectories, which support the importance of Tm flexibility to Tm-F-actin binding. We identified key residues of Tm involved in its dynamic interactions with F-actin, many of which have been found in recent mutational studies to be functionally important, and the rest of which will make promising targets for future mutational experiments.

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.

Design of a truncated cardiotoxin-I analogue with potent insulinotropic activity

Insulin secretion by pancreatic β-cells in response to glucose or other secretagogues is tightly coupled to membrane potential. Various studies have highlighted the prospect of enhancing insulin secretion in a glucose-dependent manner by blocking voltage-gated potassium channels (K(v)) and calcium-activated potassium channels (K(Ca)). Such strategy is expected to present a lower risk for hypoglycemic events compared to KATP channel blockers. Our group recently reported the discovery of a new insulinotropic agent, cardiotoxin-I (CTX-I), from the Naja kaouthia snake venom. In the present study, we report the design and synthesis of [Lys(52)]CTX-I(41-60) via structure-guided modification, a truncated, equipotent analogue of CTX-I, and demonstrate, using various pharmacological inhibitors, that this derivative probably exerts its action through Kv channels. This new analogue could represent a useful pharmacological tool to study β-cell physiology or even open a new therapeutic avenue for the treatment of type 2 diabetes.

Hydrophobic plug functions as a gate in voltage-gated proton channels

Voltage-gated proton (Hv1) channels play important roles in the respiratory burst, in pH regulation, in spermatozoa, in apoptosis, and in cancer metastasis. Unlike other voltage-gated cation channels, the Hv1 channel lacks a centrally located pore formed by the assembly of subunits. Instead, the proton permeation pathway in the Hv1 channel is within the voltage-sensing domain of each subunit. The gating mechanism of this pathway is still unclear. Mutagenic and fluorescence studies suggest that the fourth transmembrane (TM) segment (S4) functions as a voltage sensor and that there is an outward movement of S4 during channel activation. Using thermodynamic mutant cycle analysis, we find that the conserved positively charged residues in S4 are stabilized by countercharges in the other TM segments both in the closed and open states. We constructed models of both the closed and open states of Hv1 channels that are consistent with the mutant cycle analysis. These structural models suggest that electrostatic interactions between TM segments in the closed state pull hydrophobic residues together to form a hydrophobic plug in the center of the voltage-sensing domain. Outward S4 movement during channel activation induces conformational changes that remove this hydrophobic plug and instead insert protonatable residues in the center of the channel that, together with water molecules, can form a hydrogen bond chain across the channel for proton permeation. This suggests that salt bridge networks and the hydrophobic plug function as the gate in Hv1 channels and that outward movement of S4 leads to the opening of this gate.

Hydrogen bonds as molecular timers for slow inactivation in voltage-gated potassium channels

Voltage-gated potassium (Kv) channels enable potassium efflux and membrane repolarization in excitable tissues. Many Kv channels undergo a progressive loss of ion conductance in the presence of a prolonged voltage stimulus, termed slow inactivation, but the atomic determinants that regulate the kinetics of this process remain obscure. Using a combination of synthetic amino acid analogs and concatenated channel subunits we establish two H-bonds near the extracellular surface of the channel that endow Kv channels with a mechanism to time the entry into slow inactivation: an intra-subunit H-bond between Asp447 and Trp434 and an inter-subunit H-bond connecting Tyr445 to Thr439. Breaking of either interaction triggers slow inactivation by means of a local disruption in the selectivity filter, while severing the Tyr445–Thr439 H-bond is likely to communicate this conformational change to the adjacent subunit(s).

DOI:http://dx.doi.org/10.7554/eLife.01289.001

Proteins are made from long chains of smaller molecules, called amino acids. These chains twist and bend into complex three-dimensional shapes, and sometimes two or more chains, or ‘subunits’, are packed into a protein. These shapes are often held together by hydrogen bonds between some of the amino acids. Moreover, since the shape of a protein defines its function, some proteins must be able to switch between different shapes to function properly.

Ion channels are proteins that form pores through cell membranes, allowing ions to flow in and out of the cell. Potassium ion channels, which are found in neurons and heart muscle cells, have four subunits that move to open or close the central pore in response to various signals.

The closing of the channels can be ‘fast’ or ‘slow’. When the channels are closed quickly (called fast inactivation), a small part of the protein ‘plugs’ the pore from the inside of the cell. However, the mechanism behind slow inactivation remained obscure. It was thought to involve hydrogen bonds between some of the bulky amino acids that are found at the edge the pore. However, testing this hypothesis—by replacing these amino acids with alternatives that cannot form hydrogen bonds—was tricky because none of the 20 naturally occurring amino acids were alike enough to be suitable replacements.

Now, Pless et al. have overcome this limitation by using synthetic amino acids that form hydrogen bonds that are stronger or weaker than those formed by the amino acids they are replacing. The results suggest that two types of hydrogen bond keep the pore open: one is a bond between two amino acids in the same subunit, and the other is an inter-subunit bond between amino acids in neighbouring subunits. Pless et al. suggest that opening the channel causes small movements that gradually weaken, and eventually break, these bonds in one of the four subunits. Specific amino acids within the pore are then free to twist and—via a cascade of similar movements in the other three subunits—block the pore and halt the flow of ions. As such, these networks of hydrogen bonds act as pre-set breaking points allowing channels to close, even in response to continued stimulation.

Since regulated potassium channel activity underpins healthy neurons and heart muscles; understanding what controls their inactivation rate may lead to new approaches to tune their activity and treatments for important diseases.

DOI:http://dx.doi.org/10.7554/eLife.01289.002

All-atom molecular dynamics simulations of actin-myosin interactions: a comparative study of cardiac α myosin, β myosin, and fast skeletal muscle myosin

Myosins are a superfamily of actin-binding motor proteins with significant variations in kinetic properties (such as actin binding affinity) between different isoforms. It remains unknown how such kinetic variations arise from the structural and dynamic tuning of the actin-myosin interface at the amino acid residue level. To address this key issue, we have employed molecular modeling and simulations to investigate, with atomistic details, the isoform dependence of actin-myosin interactions in the rigor state. By combining electron microscopy-based docking with homology modeling, we have constructed three all-atom models for human cardiac α and β and rabbit fast skeletal muscle myosin in complex with three actin subunits in the rigor state. Starting from these models, we have performed extensive all-atom molecular dynamics (MD) simulations (total of 100 ns per system) and then used the MD trajectories to calculate actin-myosin binding free energies with contributions from both electrostatic and nonpolar forces. Our binding calculations are in good agreement with the experimental finding of isoform-dependent differences in actin binding affinity between these myosin isoforms. Such differences are traced to changes in actin-myosin electrostatic interactions (i.e., hydrogen bonds and salt bridges) that are highly dynamic and involve several flexible actin-binding loops. By partitioning the actin-myosin binding free energy to individual myosin residues, we have also identified key myosin residues involved in the actin-myosin interactions, some of which were previously validated experimentally or implicated in cardiomyopathy mutations, and the rest make promising targets for future mutational experiments.

The insecticidal potential of venom peptides

Pest insect species are a burden to humans as they destroy crops and serve as vectors for a wide range of diseases including malaria and dengue. Chemical insecticides are currently the dominant approach for combating these pests. However, the de-registration of key classes of chemical insecticides due to their perceived ecological and human health risks in combination with the development of insecticide resistance in many pest insect populations has created an urgent need for improved methods of insect pest control. The venoms of arthropod predators such as spiders and scorpions are a promising source of novel insecticidal peptides that often have different modes of action to extant chemical insecticides. These peptides have been optimized via a prey-predator arms race spanning hundreds of millions of years to target specific types of insect ion channels and receptors. Here we review the current literature on insecticidal venom peptides, with a particular focus on their structural and pharmacological diversity, and discuss their potential for deployment as insecticides.

Unique mechanism of the interaction between honey bee toxin TPNQ and rKir1.1 potassium channel explored by computational simulations: insights into the relative insensitivity of channel towards animal toxins

BACKGROUND

The 21-residue compact tertiapin-Q (TPNQ) toxin, a derivative of honey bee toxin tertiapin (TPN), is a potent blocker of inward-rectifier K(+) channel subtype, rat Kir1.1 (rKir1.1) channel, and their interaction mechanism remains unclear.

PRINCIPAL FINDINGS

Based on the flexible feature of potassium channel turrets, a good starting rKir1.1 channel structure was modeled for the accessibility of rKir1.1 channel turrets to TPNQ toxin. In combination with experimental alanine scanning mutagenesis data, computational approaches were further used to obtain a reasonable TPNQ toxin-rKir1.1 channel complex structure, which was completely different from the known binding modes between animal toxins and potassium channels. TPNQ toxin mainly adopted its helical domain as the channel-interacting surface together with His12 as the pore-blocking residue. The important Gln13 residue mainly contacted channel residues near the selectivity filter, and Lys20 residue was surrounded by a polar "groove" formed by Arg118, Thr119, Glu123, and Asn124 in the channel turret. On the other hand, four turrets of rKir1.1 channel gathered to form a narrow pore entryway for TPNQ toxin recognition. The Phe146 and Phe148 residues in the channel pore region formed strong hydrophobic protrusions, and produced dominant nonpolar interactions with toxin residues. These specific structure features of rKir1.1 channel vestibule well matched the binding of potent TPNQ toxin, and likely restricted the binding of the classical animal toxins.

CONCLUSIONS/SIGNIFICANCE

The TPNQ toxin-rKir1.1 channel complex structure not only revealed their unique interaction mechanism, but also would highlight the diverse animal toxin-potassium channel interactions, and elucidate the relative insensitivity of rKir1.1 channel towards animal toxins.

Complex structures between the N-type calcium channel (CaV2.2) and ω-conotoxin GVIA predicted via molecular dynamics

The N-type voltage-gated Ca(2+) channel CaV2.2 is one of the important targets for pain management. ω-Conotoxins isolated from venoms of cone snails, which specifically inhibit CaV2.2, are promising scaffolds for novel analgesics. The inhibitory action of ω-conotoxins on CaV2.2 has been examined experimentally, but the modes of binding of the toxins to this and other related subfamilies of Ca(2+) channels are not understood in detail. Here molecular dynamics simulations are used to construct models of ω-conotoxin GVIA in complex with a homology model of the pore domain of CaV2.2. Three different binding modes in which the side chain of Lys2, Arg17, or Lys24 from the toxin protrudes into the selectivity filter of CaV2.2 are considered. In all the modes, the toxin forms a salt bridge with an aspartate residue of subunit II just above the EEEE ring of the selectivity filter. Using the umbrella sampling technique and potential of mean force calculations, the half-maximal inhibitory concentration (IC50) values are calculated to be 1.5 and 0.7 nM for the modes in which Lys2 and Arg17 occlude the ion conduction pathway, respectively. Both IC50 values compare favorably with the values of 0.04-1.0 nM determined experimentally. The similar IC50 values calculated for the different binding modes demonstrate that GVIA can inhibit CaV2.2 with alternative binding modes. Such a multiple-binding mode mechanism may be common for ω-conotoxins.

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.

Nano-positioning system for structural analysis of functional homomeric proteins in multiple conformations

Proteins may undergo multiple conformational changes required for their function. One strategy used to estimate target-site positions in unknown structural conformations involves single-pair resonance energy transfer (RET) distance measurements. However, interpretation of inter-residue distances is difficult when applied to three-dimensional structural rearrangements, especially in homomeric systems. We developed a positioning method using inverse trilateration/triangulation to map target sites within a homomeric protein in all defined states, with simultaneous functional recordings. The procedure accounts for probe diffusion to accurately determine the three-dimensional position and confidence region of lanthanide LRET donors attached to a target site (one per subunit), relative to a single fluorescent acceptor placed in a static site. As first application, the method is used to determine the position of a functional voltage-gated potassium channel's voltage sensor. Our results verify the crystal structure relaxed conformation and report on the resting and active conformations for which crystal structures are not available.

Structural basis of the selective block of Kv1.2 by maurotoxin from computer simulations

The 34-residue polypeptide maurotoxin (MTx) isolated from scorpion venoms selectively inhibits the current of the voltage-gated potassium channel Kv1.2 by occluding the ion conduction pathway. Here using molecular dynamics simulation as a docking method, the binding modes of MTx to three closely related channels (Kv1.1, Kv1.2 and Kv1.3) are examined. We show that MTx forms more favorable electrostatic interactions with the outer vestibule of Kv1.2 compared to Kv1.1 and Kv1.3, consistent with the selectivity of MTx for Kv1.2 over Kv1.1 and Kv1.3 observed experimentally. One salt bridge in the bound complex of MTx-Kv1.2 forms and breaks in a simulation period of 20ns, suggesting the dynamic nature of toxin-channel interactions. The toxin selectivity likely arises from the differences in the shape of the channel outer vestibule, giving rise to distinct orientations of MTx on block. Potential of mean force calculations show that MTx blocks Kv1.1, Kv1.2 and Kv1.3 with an IC₅₀ value of 6 µM, 0.6nM and 18 µM, respectively.