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

Peptide Toxins Targeting KV Channels

A number of peptide toxins isolated from animals target potassium ion (K⁺) channels. Many of them are particularly known to inhibit voltage-gated K⁺ (K_V) channels and are mainly classified into pore-blocking toxins or gating-modifier toxins. Pore-blocking toxins directly bind to the ion permeation pores of K_V channels, thereby physically occluding them. In contrast, gating-modifier toxins bind to the voltage-sensor domains of K_V channels, modulating their voltage-dependent conformational changes. These peptide toxins are useful molecular tools in revealing the structure-function relationship of K_V channels and have potential for novel treatments for diseases related to K_V channels. This review focuses on the inhibition mechanism of pore-blocking and gating-modifier toxins that target K_V channels.

DelPhi Suite: New Developments and Review of Functionalities

Electrostatic potential, energies, and forces affect virtually any process in molecular biology, however, computing these quantities is a difficult task due to irregularly shaped macromolecules and the presence of water. Here, we report a new edition of the popular software package DelPhi along with describing its functionalities. The new DelPhi is a C++ object-oriented package supporting various levels of multiprocessing and memory distribution. It is demonstrated that multiprocessing results in significant improvement of computational time. Furthermore, for computations requiring large grid size (large macromolecular assemblages), the approach of memory distribution is shown to reduce the requirement of RAM and thus permitting large-scale modeling to be done on Linux clusters with moderate architecture. The new release comes with new features, whose functionalities and applications are described as well. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Structural and Functional Analyses of Cone Snail Toxins

Cone snails are marine gastropod mollusks with one of the most powerful venoms in nature. The toxins, named conotoxins, must act quickly on the cone snails´ prey due to the fact that snails are extremely slow, reducing their hunting capability. Therefore, the characteristics of conotoxins have become the object of investigation, and as a result medicines have been developed or are in the trialing process. Conotoxins interact with transmembrane proteins, showing specificity and potency. They target ion channels and ionotropic receptors with greater regularity, and when interaction occurs, there is immediate physiological decompensation. In this review we aimed to evaluate the structural features of conotoxins and the relationship with their target types.

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

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

Defensins, a novel type of animal toxin-like potassium channel inhibitor

The classical potassium channel inhibitors are toxin peptides from venomous animals, and whether there are peptide inhibitors from other species is an open question. Due to both the independent and interdependent relationships between the spear (peptide inhibitors) and the shield (potassium channels), human defensins were first identified by our group as endogenous potassium channel inhibitors. Encouraged by the discovery of human defensins as potassium channel inhibitors, defensins from invertebrates and fungi were successively found by our group to be potassium channel inhibitors. In addition, a plant defensin was reported to be a potassium channel inhibitor. Since defensins are widely produced by vertebrate, invertebrate, plant and fungi species, the recent work established a new research field on defensin-potassium channel interactions. Here, we review the current work on defensins from vertebrate, invertebrate, plant and fungi species as inhibitors of potassium channels and discuss future work in this research field.

Efficient enzymatic cyclization of an inhibitory cystine knot-containing peptide

Disulfide-rich peptides isolated from cone snails are of great interest as drug leads due to their high specificity and potency toward therapeutically relevant ion channels and receptors. They commonly contain the inhibitor cystine knot (ICK) motif comprising three disulfide bonds forming a knotted core. Here we report the successful enzymatic backbone cyclization of an ICK-containing peptide κ-PVIIA, a 27-amino acid conopeptide from Conus purpurascens, using a mutated version of the bacterial transpeptidase, sortase A. Although a slight loss of activity was observed compared to native κ-PVIIA, cyclic κ-PVIIA is a functional peptide that inhibits the Shaker voltage-gated potassium (Kv) channel. Molecular modeling suggests that the decrease in potency may be related to the loss of crucial, but previously unidentified electrostatic interactions between the N-terminus of the peptide and the Shaker channel. This hypothesis was confirmed by testing an N-terminally acetylated κ-PVIIA, which shows a similar decrease in activity. We also investigated the conformational dynamics and hydrogen bond network of cyc-PVIIA, both of which are important factors to be considered for successful cyclization of peptides. We found that cyc-PVIIA has the same conformational dynamics, but different hydrogen bond network compared to those of κ-PVIIA. The ability to efficiently cyclize ICK peptides using sortase A will enable future protein engineering for this class of peptides and may help in the development of novel therapeutic molecules. Biotechnol. Bioeng. 2016;113: 2202-2212. © 2016 Wiley Periodicals, Inc.

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.

Engineered nanoparticles mimicking cell membranes for toxin neutralization

Protein toxins secreted from pathogenic bacteria and venomous animals rely on multiple mechanisms to overcome the cell membrane barrier to inflict their virulence effect. A promising therapeutic concept toward developing a broadly applicable anti-toxin platform is to administer cell membrane mimics as decoys to sequester these virulence factors. As such, lipid membrane-based nanoparticulates are an ideal candidate given their structural similarity to cellular membranes. This article reviews the virulence mechanisms employed by toxins at the cell membrane interface and highlights the application of cell-membrane mimicking nanoparticles as toxin decoys for systemic detoxification. In addition, the implication of particle/toxin nanocomplexes in the development of toxoid vaccines is discussed.

Toxin acidic residue evolutionary function-guided design of de novo peptide drugs for the immunotherapeutic target, the Kv1.3 channel

During the long-term evolution of animal toxins acting on potassium channels, the acidic residues can orientate the toxin binding interfaces by adjusting the molecular polarity. Based on the evolutionary function of toxin acidic residues, de novo peptide drugs with distinct binding interfaces were designed for the immunotherapeutic target, the Kv1.3 channel. Using a natural basic toxin, BmKTX, as a template, which contains 2 acidic residues (Asp19 and Asp33), we engineered two new peptides BmKTX-19 with 1 acidic residue (Asp33), and BmKTX-196 with 2 acidic residues (Asp6 and Asp33) through only adjusting acidic residue distribution for reorientation of BmKTX binding interface. Pharmacological experiments indicated that BmKTX-19 and BmKTX-196 peptides were specific inhibitors of the Kv1.3 channel and effectively suppressed cytokine secretion. In addition to the structural similarity between the designed and native peptides, both experimental alanine-scanning mutagenesis and computational simulation further indicated that the binding interface of wild-type BmKTX was successfully reoriented in BmKTX-19 and BmKTX-196, which adopted distinct toxin surfaces as binding interfaces. Together, these findings indicate not only the promising prospect of BmKTX-19 and BmKTX-196 as drug candidates but also the desirable feasibility of the evolution-guided peptide drug design for discovering numerous peptide drugs for the Kv1.3 channel.

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.

Computational studies of marine toxins targeting ion channels

Toxins from marine animals offer novel drug leads for treatment of diseases involving ion channels. Computational methods could be very helpful in this endeavour in several ways, e.g., (i) constructing accurate models of the channel-toxin complexes using docking and molecular dynamics (MD) simulations; (ii) determining the binding free energies of toxins from umbrella sampling MD simulations; (iii) predicting the effect of mutations from free energy MD simulations. Using these methods, one can design new analogs of toxins with improved affinity and selectivity properties. Here we present a review of the computational methods and discuss their applications to marine toxins targeting potassium and sodium channels. Detailed examples from the potassium channel toxins—ShK from sea anemone and κ-conotoxin PVIIA—are provided to demonstrate capabilities of the computational methods to give accurate descriptions of the channel-toxin complexes and the energetics of their binding. An example is also given from sodium channel toxins (µ-conotoxin GIIIA) to illustrate the differences between the toxin binding modes in potassium and sodium channels.

Why the Drosophila Shaker K+ channel is not a good model for ligand binding to voltage-gated Kv1 channels

The Drosophila Shaker K(+) channel is the first cloned voltage-gated potassium channel and has, therefore, played an important role in structural and functional studies of those channels. While such a role is well justified for ion permeation, it is not clear whether this also extends to ligand binding. Despite the high degree of homology among Shaker and Kv1 channels, κ-conotoxin PVIIA (κ-PVIIA) binds to Shaker with high affinity but not to Kv1 channels. Here we address this issue by studying binding of κ-PVIIA to Shaker and Kv1 channels using molecular dynamics (MD) simulations. The structures of the channel-toxin complexes are constructed via docking and refinement with MD. The binding mode of each complex is characterized and compared to available mutagenesis data to validate the complex models. The potential of mean force for dissociation of the Shaker-κ-PVIIA complex is calculated from umbrella sampling MD simulations, and the corresponding binding free energy is determined, which provides further validation of the complex structure. Comparison of the Shaker and Kv1 complex models shows that a few mutations in the turret and extended regions are sufficient to abolish the observed sensitivity of Shaker to κ-PVIIA. This study demonstrates that Shaker is not always a good model for Kv1 channels for ligand binding. It also provides insights into the binding of the toxin to potassium channels that will be useful for improving affinity and selectivity properties of Kv1 channels.

Use-dependent block of the voltage-gated Na(+) channel by tetrodotoxin and saxitoxin: effect of pore mutations that change ionic selectivity

Voltage-gated Na(+) channels (NaV channels) are specifically blocked by guanidinium toxins such as tetrodotoxin (TTX) and saxitoxin (STX) with nanomolar to micromolar affinity depending on key amino acid substitutions in the outer vestibule of the channel that vary with NaV gene isoforms. All NaV channels that have been studied exhibit a use-dependent enhancement of TTX/STX affinity when the channel is stimulated with brief repetitive voltage depolarizations from a hyperpolarized starting voltage. Two models have been proposed to explain the mechanism of TTX/STX use dependence: a conformational mechanism and a trapped ion mechanism. In this study, we used selectivity filter mutations (K1237R, K1237A, and K1237H) of the rat muscle NaV1.4 channel that are known to alter ionic selectivity and Ca(2+) permeability to test the trapped ion mechanism, which attributes use-dependent enhancement of toxin affinity to electrostatic repulsion between the bound toxin and Ca(2+) or Na(+) ions trapped inside the channel vestibule in the closed state. Our results indicate that TTX/STX use dependence is not relieved by mutations that enhance Ca(2+) permeability, suggesting that ion-toxin repulsion is not the primary factor that determines use dependence. Evidence now favors the idea that TTX/STX use dependence arises from conformational coupling of the voltage sensor domain or domains with residues in the toxin-binding site that are also involved in slow inactivation.

In silico analysis of conformational changes induced by mutation of aromatic binding residues: consequences for drug binding in the hERG K+ channel

Pharmacological inhibition of cardiac hERG K⁺ channels is associated with increased risk of lethal arrhythmias. Many drugs reduce hERG current by directly binding to the channel, thereby blocking ion conduction. Mutation of two aromatic residues (F656 and Y652) substantially decreases the potency of numerous structurally diverse compounds. Nevertheless, some drugs are only weakly affected by mutation Y652A. In this study we utilize molecular dynamics simulations and docking studies to analyze the different effects of mutation Y652A on a selected number of hERG blockers. MD simulations reveal conformational changes in the binding site induced by mutation Y652A. Loss of π-π-stacking between the two aromatic residues induces a conformational change of the F656 side chain from a cavity facing to cavity lining orientation. Docking studies and MD simulations qualitatively reproduce the diverse experimentally observed modulatory effects of mutation Y652A and provide a new structural interpretation for the sensitivity differences.

Molecular mechanism of δ-dendrotoxin-potassium channel recognition explored by docking and molecular dynamic simulations

δ-Dendrotoxin, isolated from mamba snake venom, has 57 residues cross-linked by three disulfide bridges. The protein shares a pharmacological activity with other animal toxins, the potent blockade of potassium channels, but is structurally unrelated to toxins of different species. We employed alanine-scanning mutagenesis to explore the molecular mechanism of δ-dendrotoxin binding to potassium channels, using protein-protein docking and molecular dynamic simulations. In our reasonable model of the δ-dendrotoxin-ShaKv1.1 complex, δ-dendrotoxin interacted mainly with the N-terminal region and the turn of two antiparallel β-sheets of the channel. This binding mode could well explain the functional roles of critical residues in δ-dendrotoxin and the ShaKv1.1 channel. Structural analysis indicated that the critical Lys6 residue of δ-dendrotoxin plugged its side chain into a channel selectivity filter. Another two critical δ-dendrotoxin residues, Lys3 and Arg10, were found to contact channel residues through strong polar and nonpolar interactions, especially salt-bridge interactions. As for the ShaKv1.1 channel, the channel turrets were found in the "half-open state," and two of four Glu423 in the turrets of the channel B and D chains could interact, respectively, with Lys3 and Lys26 of δ-dendrotoxin through electrostatic interactions. The essential Asp431 channel residue was found to associate electrostatically with Arg10 of δ-dendrotoxin, and a critical Tyr449 channel residue was just under the channel-interacting surface of δ-dendrotoxin. Together, these novel data may accelerate the structure-function research of toxins in the dendrotoxin family and be of significant value in revealing the diverse interactions between animal toxins and potassium channels.

Use of venom peptides to probe ion channel structure and function

Venoms of snakes, scorpions, spiders, insects, sea anemones, and cone snails are complex mixtures of mostly peptides and small proteins that have evolved for prey capture and/or defense. These deadly animals have long fascinated scientists and the public. Early studies isolated lethal components in the search for cures and understanding of their mechanisms of action. Ion channels have emerged as targets for many venom peptides, providing researchers highly selective and potent molecular probes that have proved invaluable in unraveling ion channel structure and function. This minireview highlights molecular details of their toxin-receptor interactions and opportunities for development of peptide therapeutics.

The outer vestibule of the Na+ channel-toxin receptor and modulator of permeation as well as gating

The outer vestibule of voltage-gated Na(+) channels is formed by extracellular loops connecting the S5 and S6 segments of all four domains ("P-loops"), which fold back into the membrane. Classically, this structure has been implicated in the control of ion permeation and in toxin blockage. However, conformational changes of the outer vestibule may also result in alterations in gating, as suggested by several P-loop mutations that gave rise to gating changes. Moreover, partial pore block by mutated toxins may reverse gating changes induced by mutations. Therefore, toxins that bind to the outer vestibule can be used to modulate channel gating.

Spontaneous conformational change and toxin binding in alpha7 acetylcholine receptor: insight into channel activation and inhibition

Nicotinic AChRs (nAChRs) represent a paradigm for ligand-gated ion channels. Despite intensive studies over many years, our understanding of the mechanisms of activation and inhibition for nAChRs is still incomplete. Here, we present molecular dynamics (MD) simulations of the alpha7 nAChR ligand-binding domain, both in apo form and in alpha-Cobratoxin-bound form, starting from the respective homology models built on crystal structures of the acetylcholine-binding protein. The toxin-bound form was relatively stable, and its structure was validated by calculating mutational effects on the toxin-binding affinity. However, in the apo form, one subunit spontaneously moved away from the conformation of the other four subunits. This motion resembles what has been proposed for leading to channel opening. At the top, the C loop and the adjacent beta7-beta8 loop swing downward and inward, whereas at the bottom, the F loop and the C terminus of beta10 swing in the opposite direction. These swings appear to tilt the whole subunit clockwise. The resulting changes in solvent accessibility show strong correlation with experimental results by the substituted cysteine accessibility method upon addition of acetylcholine. Our MD simulation results suggest a mechanistic model in which the apo form, although predominantly sampling the "closed" state, can make excursions into the "open" state. The open state has high affinity for agonists, leading to channel activation, whereas the closed state upon distortion has high affinity for antagonists, leading to inhibition.

Molecular basis of inhibitory peptide maurotoxin recognizing Kv1.2 channel explored by ZDOCK and molecular dynamic simulations

Inhibitory peptide-channel interactions have been utilized to characterize both channels and peptides; however, the fundamental basis for these interactions remains elusive. Here, combined computation methods were employed to study the specific binding of maurotoxin (MTX) peptide to Kv1.2 channel. In the first stage, numerous predicted complexes were generated by docking an ensemble of all 35 NMR conformations of MTX to Kv1.2 channel with ZDOCK program. Then the resulted complexes were clustered and classified into four main binding modes, based on experimental information and interaction energy analysis after the energy minimization and molecular dynamics (MD) simulations. By examining the stability of the plausible candidates through unrestrained MD simulations and calculation of the binding free energies, a final reasonable MTX-Kv1.2 complex was identified, with an overall high degree of correlation between the calculation and experiment on mutational effects. In the obtained complex structure model, MTX mainly used its beta-sheet domains to associate the channel mouth instead of the well-recognized functionally important S5P linkers of Kv1.2 channel. Structure analysis characterized that the most essential Tyr(32) residue of MTX was surrounded by a "pocket" formed by many nonpolar and polar residues of Kv1.2 channel, and revealed a pore-blocking Lys(23) and an important Lys(7) stabilized by strong electrostatic interactions with Asp(379) of Kv1.2. Furthermore, a stepwise structural arrangement for both ligand and receptor was found to accompany the tighter interaction of MTX into the target channel. The starting conformation of MTX, the side-chain conformation of the most important residue Tyr(32), and proper introduction of flexibility for candidate complexes were demonstrated to be considerably important factors for obtaining the final reasonable complex structure model. All these findings should not only be helpful for identifying more plausible K(+) channel-inhibitory peptide complex structures, but also provide intrinsically valuable structural biology information to interpret binding affinities, specificities, and diversity of K(+) channel-nature toxin interactions.

On the Dielectric Boundary in Poisson-Boltzmann Calculations

In applying the Poisson-Boltzmann (PB) equation for calculating the electrostatic free energies of solute molecules, an open question is how to specify the boundary between the low-dielectric solute and the high-dielectric solvent. Two common specifications of the dielectric boundary, as the molecular surface (MS) or the van der Waals (vdW) surface of the solute, give very different results for the electrostatic free energy of the solute. With the same atomic radii, the solute is more solvent-exposed in the vdW specification. One way to resolve the difference is to use different sets of atomic radii for the two surfaces. The radii for the vdW surface would be larger in order to compensate for the higher solvent exposure. Here we show that radius re-parameterization required for bringing MS-based and vdW-based PB results to agreement is solute-size dependent. The difference in atomic radii for individual amino acids as solutes is only 2-5% but increases to over 20% for proteins with ~200 residues. Therefore two sets of radii that yield identical MS-based and vdW-based PB results for small solutes will give very different PB results for large solutes. This finding raises issues about two common practices. The first is the use of atomic radii, which are parameterized against either experimental solvation data or data obtained from explicit-solvent simulations on small compounds, for PB calculations on proteins. The second is the parameterization of vdW-based generalized Born models against MS-based PB results.

A holistic approach to protein docking

Docking of unbound protein structures into a complex has gained significant progress in recent years, but nonetheless still poses a great challenge. We have pursued a holistic approach to docking which brings together effective methods at different stages. First, protein-protein interaction sites are predicted or obtained from experimental studies in the literature. Interface prediction/experimental data are then used to guide the generation of docked poses or to rank docked poses generated from an unbiased search. Finally, selected models are refined by lengthy molecular dynamics (MD) simulations in explicit water. For CAPRI target T27, we used information on interaction sites as input to drive docking and as a filter to rank docked poses. Lead candidates were then clustered according to RMSD among them. From the clustering, 10 models were selected and subject to refinement by MD simulations. Our Model 7 is rated number one among all submissions according to L_rmsd. Six of our other submissions are rated acceptable. As scorer, eight of our submissions are rated acceptable.

Prediction of salt and mutational effects on the association rate of U1A protein and U1 small nuclear RNA stem/loop II

We have developed a computational approach for predicting protein-protein association rates (Alsallaq and Zhou, Structure 2007, 15, 215). Here we expand the range of applicability of this approach to protein-RNA binding and report the first results for protein-RNA binding rates predicted from atomistic modeling. The system studied is the U1A protein and stem/loop II of the U1 small nuclear RNA. Experimentally it was observed that the binding rate is significantly reduced by increasing salt concentration while the dissociation changes little with salt concentration, and charges distant from the binding site make marginal contribution to the binding rate. These observations are rationalized. Moreover, predicted effects of salt and charge mutations are found to be in quantitative agreement with experimental results.

Speeding the recovery from ultraslow inactivation of voltage-gated Na+ channels by metal ion binding to the selectivity filter: a foot-on-the-door?

Slow inactivated states in voltage-gated ion channels can be modulated by binding molecules both to the outside and to the inside of the pore. Thus, external K(+) inhibits C-type inactivation in Shaker K(+) channels by a "foot-in-the-door" mechanism. Here, we explore the modulation of a very long-lived inactivated state, ultraslow inactivation (I(US)), by ligand binding to the outer vestibule in voltage-gated Na(+) channels. Blocking the outer vestibule by a mutant mu-conotoxin GIIIA substantially accelerated recovery from I(US). A similar effect was observed if Cd(2+) was bound to a cysteine engineered to the selectivity filter (K1237C). In K1237C channels, exposed to 30 microM Cd(2+), the time constant of recovery from I(US) was decreased from 145.0 +/- 10.2 s to 32.5 +/- 3.3 s (P < 0.001). Recovery from I(US) was only accelerated if Cd(2+) was added to the bath solution during recovery (V = -120 mV) from I(US), but not when the channels were selectively exposed to Cd(2+) during the development of I(US) (-20 mV). These data could be explained by a kinetic model in which Cd(2+) binds with high affinity to a slow inactivated state (I(S)), which is transiently occupied during recovery from I(US). A total of 50 microM Cd(2+) produced an approximately 8 mV hyperpolarizing shift of the steady-state inactivation curve of I(S), supporting this kinetic model. Binding of lidocaine to the internal vestibule significantly reduced the number of channels entering I(US), suggesting that I(US) is associated with a conformational change of the internal vestibule of the channel. We propose a molecular model in which slow inactivation (I(S)) occurs by a closure of the outer vestibule, whereas I(US) arises from a constriction of the internal vestibule produced by a widening of the selectivity filter region. Binding of Cd(2+) to C1237 promotes the closure of the selectivity filter region, thereby hastening recovery from I(US). Thus, Cd(2+) ions may act like a foot-on-the-door, kicking the I(S) gate to close.

Quantitative computer simulations of biomolecules: A snapshot

A recent workshop titled "Quantitative Computational Biophysics" at Florida State University provided an overview of the state of the art in quantitative modeling of biomolecular systems. The presentations covered a wide range of interrelated topics, including the development and validation of force fields, the modeling of protein-protein interactions, the sampling of conformational space, and the assessment of equilibration and statistical errors. Substantial progress in all these areas was reported.

Prediction of protein-protein association rates from a transition-state theory

We recently developed a theory for the rates of protein-protein association. The theory is based on the concept of a transition state, which separates the bound state, with numerous short-range interactions but restricted translational and rotational freedom, and the unbound state, with, at most, a small number of interactions but expanded configurational freedom. When not accompanied by large-scale conformational changes, protein-protein association becomes diffusion limited. The association rate is then predicted as k(a)=k(a)(0)exp(-DeltaG(el)(double dagger)/k(B)T), where DeltaG(el)(double dagger) is the electrostatic interaction free energy in the transition state, k(a)(0) is the rate in the absence of electrostatic interactions, and k(B)T is thermal energy. Here, this transition-state theory is used to predict the association rates of four protein complexes. The predictions for the wild-type complexes and 23 mutants are found to agree closely with experimental data over wide ranges of ionic strength.

Energy landscape and transition state of protein-protein association

Formation of a stereospecific protein complex is favored by specific interactions between two proteins but disfavored by the loss of translational and rotational freedom. Echoing the protein folding process, we have previously proposed a transition state for protein-protein association. Here we clarify the specification of the transition state by working with two types of toy models for protein association. A "hemisphere" model consists of two matching hemispheres as associating proteins, and a "crater" model consists of a spherical protein with a crater to which another spherical protein fits snugly. Short-range pairwise interactions between sites across the interface hold together the bound complex. Small relative translation and rotation between the subunits quickly destroy the interactions, leading to a sharp transition between the bound state with numerous short-range interactions but restricted translation and rotational freedom and the unbound state with, at most, a small number of interactions but expanded configurational freedom. This transition sets the outer boundary of the bound state as well as the transition state for association. The energy landscape is funnel-like, with the deep well of the bound state surrounded by a broad shallow basin. Calculations with the toy models suggest that mutational change in the interaction energy in the x-ray structure of a protein-protein complex, commonly used to approximate the effect on the association constant, overestimates the effect of mutation by 10-20%. With an eye toward specifying the transition states of actual protein complexes, we find that the total number of contacts between the subunits serves as a good surrogate of the interaction energy. This formalism of protein-protein association is applied to the barnase-barstar complex, reproducing the experimental results for the association rate over a wide range of ionic strength.

Mechanism of Ion Permeation in a Model Channel: Free Energy Surface and Dynamics of K+Ion Transport in an Anion-Doped Carbon Nanotube

The mechanism of the ion permeation is investigated for an anion-doped carbon nanotube, as a model of the K+ channel, by analyzing the free energy surface and the dynamics of the ion permeation through the model channel. It is found that the main rate-determining step is how an ion enters the channel. The entrance of the ion is mostly blocked by a water molecule located at this entrance. Only about 10% of K+ ions which reach the mouth of the channel can really enter the channel. The rejection rate sensitively depends on the location of this water molecule, which is easily controlled by the charge of the carbon nanotube; for example, the maximum permeation is obtained when the anion charge is at a certain value, -5.4e in the present model. At this charge, the facile translocation of the ion inside the channel is also induced due to the number of fluctuations of the ions inside the channel. Therefore, the so-called "Newton's balls", a toy model, combined with a simple ion diffusion model for explaining the fast ion permeation should be modified. The present analysis thus suggests that there exists an optimum combination of the length and the charge of the carbon nanotube for the most efficient ion permeation.

BcIV, a new paralyzing peptide obtained from the venom of the sea anemone Bunodosoma caissarum. A comparison with the Na+ channel toxin BcIII

Sea anemones produce a wide variety of biologically active compounds, such as the proteinaceous neurotoxins and cytolysins. Herein we report a new peptide, purified to homogeneity from the neurotoxic fraction of B. caissarum venom, by using gel filtration followed by rp-HPLC, naming it as BcIV. BcIV is a 41 amino acid peptide (molecular mass of 4669 amu) possessing 6 cysteines covalently linked by three disulfide bonds. This toxin has 45 and 48% of identity when compared to APETx1 and APETx2 from Anthopleura elegantissima, respectively, and 42% of identity with Am-II and BDS-I and-II obtained from Antheopsis maculata and Anemonia sulcata, respectively. This neurotoxin presents only a weak-paralyzing action (minimal Lethal Dose close to 2000 microg/kg) in swimming crabs Callinectes danae. This appears to be a different effect to that caused by the type 1 sea anemone toxin BcIII that is lethal to the same animals at lower doses (LD50=219 microg/kg). Circular dichroism spectra of BcIII and BcIV show a high content of beta-strand secondary structure in both peptides, very similar to type 1 sodium channel toxins from various sea anemones, and to APETx1 and APETx2 from A. elegantissima, a HERG channel modulator and an ASIC3 inhibitor, respectively. Interestingly, BcIII and BcIV have similar effects on the action potential of the crab leg nerves, suggesting the same target in this tissue. As BcIII was previously reported as a Na+ channel effector and BcIV is inactive over Na+ currents of mammalian GH3 cells, we propose a species-specific action for this new molecule. A molecular model of BcIV was constructed using the structure of the APETx1 as template and putative key residues are discussed.

Solution structure of PcFK1, a spider peptide active against Plasmodium falciparum

Psalmopeotoxin I (PcFK1) is a 33-amino-acid residue peptide isolated from the venom of the tarantula Psalmopoeus cambridgei. It has been recently shown to possess strong antiplasmodial activity against the intra-erythrocyte stage of Plasmodium falciparum in vitro. Although the molecular target for PcFK1 is not yet determined, this peptide does not lyse erythrocytes, is not cytotoxic to nucleated mammalian cells, and does not inhibit neuromuscular function. We investigated the structural properties of PcFK1 to help understand the unique mechanism of action of this peptide and to enhance its utility as a lead compound for rational development of new antimalarial drugs. In this paper, we have determined the three-dimensional solution structure by (1)H two-dimensional NMR means of recombinant PcFK1, which is shown to belong to the ICK structural superfamily with structural determinants common to several neurotoxins acting as ion channels effectors.

I-conotoxin superfamily revisited

The I-conotoxin superfamily (I-Ctx) is known to have four disulfide bonds with the cysteine arrangement C-C-CC-CC-C-C, and the members inhibit or modify ion channels of nerve cells. Recently, Olivera and co-workers (FEBS J. 2005; 272: 4178-4188) have suggested that the previously described I-Ctx should now be divided into two different gene superfamilies, namely, I1 and I2, in view of their having two different types of signal peptides and exhibiting distinct functions. We have revisited the 28 entries presently grouped as I-Ctx in UniProt Swiss-Prot knowledgebase, and on the basis of in silico analysis have divided them into I1 and I2 superfamilies. The sequence analysis has provided a framework for in silico annotation enabling us to carry out computer-based functional characterization of the UniProtKB/TrEMBL entry Q59AA4 from Conus miles and to predict it as a member of the I2 superfamily. Furthermore, we have predicted the mature toxin of this entry and have proposed that it may be an inhibitor of voltage-gated potassium channels.