Nuclear magnetic resonance solution structure of the alpha-neurotoxin from the black mamba (Dendroaspis polylepis polylepis)

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.

Solution structure of gamma-bungarotoxin: the functional significance of amino acid residues flanking the RGD motif in integrin binding

Gamma-bungarotoxin, a snake venom protein isolated from Bungarus multicinctus, contains 68 amino acids, including 10 cysteine residues and a TAVRGDGP sequence at positions 30-37. The solution structure of gamma-bungarotoxin has been determined by nuclear magnetic resonance (NMR) spectroscopy. The structure is similar to that of the short-chain neurotoxins that contain three loops extending from a disulfide-bridged core. The tripeptide Arg-Gly-Asp (RGD) sequence is located at the apex of the flexible loop and is similar to that of other RGD-containing proteins. However, gamma-bungarotoxin only inhibits platelet aggregations with an IC50 of 34 microM. To understand its weak activity in inhibiting platelet aggregation, we mutated the RGD loop sequences of rhodostomin, a potent platelet aggregation inhibitor, from RIPRGDMP to TAVRGDGP, resulting in a 196-fold decrease in activity. In addition, the average Calpha-to-Calpha distance between R33 and G36 of gamma-bungarotoxin is 6.02 A, i.e., shorter than that of other RGD-containing proteins that range from 6.55 to 7.46 A. These results suggested that the amino acid residues flanking the RGD motif might control the width of the RGD loop. This structural difference may be responsible for its decrease in platelet aggregation inhibition compared with other RGD-containing proteins.

Structurally conserved α-neurotoxin genes encode functionally diverse proteins in the venom of Naja sputatrix 1

The structure and organization of the genes encoding the long-chain neurotoxins and four other isoforms of weak neurotoxins in the venom of Naja sputatrix are reported. The genes contained three exons interrupted by two introns, a structure similar to other members of the three-finger toxin family. The proteins encoded by these genes, however, show varied affinity towards nicotinic acetylcholine receptors. Phylogenetic analysis of these genes showed that the weak neurotoxin gene is confined to a distinct group. We also observe that specific mutations of the gene provide the diversity in function in these toxins while maintaining a common structural scaffold. This forms the first report where the molecular basis of evolution of postsynaptic neurotoxins from an ancestral gene can be demonstrated using the same species of snake.

Cytotoxic potency of cardiotoxin from Naja sputatrix: development of a new cytolytic assay

The possible involvement of specific regions/loops of cardiotoxin from Naja sputatrix venom in mediating its cytolytic activity is evaluated using a new cytolytic assay. In this assay, the amount of chloramphenicol acetyltransferase (CAT) that is released upon lysis of the cellular membranes by the cytotoxin has been measured as an index of cytolysis. This newly developed CAT system is more sensitive than the traditional haemolysis method utilizing red blood cells or the lactate dehydrogenase assay for cytolysis. Series of chimaeric toxin molecules have been constructed by swapping the loops between highly hydrophilic neurotoxin and highly hydrophobic cardiotoxin molecules from Naja sputatrix, which are known to exhibit structural similarity (three-finger conformation) but to have different functional properties. Comparison of the cytolytic activities of the recombinant chimaeric toxins demonstrated the possible involvement of all three loops of cardiotoxin in its cytolytic potency. However, the first two loops of the protein appear to make the major contribution to its lytic activity. cDNAs encoding cardiotoxin and the chimaeric toxins, when expressed in transfected cultured Chinese hamster ovary cells, resulted in cell lysis, indicating that these cDNAs can be developed as useful cytolytic agents.

Classification of protein disulphide-bridge topologies

The preferential occurrence of certain disulphide-bridge topologies in proteins has prompted us to design a method and a program, KNOT-MATCH, for their classification. The program has been applied to a database of proteins with less than 65% homology and more than two disulphide bridges. We have investigated whether there are topological preferences that can be used to group proteins and if these can be applied to gain insight into the structural or functional relationships among them. The classification has been performed by Density Search and Hierarchical Clustering Techniques, yielding thirteen main protein classes from the superimposition and clustering process. It is noteworthy that besides the disulphide bridges, regular secondary structures and loops frequently become correctly aligned. Although the lack of significant sequence similarity among some clustered proteins precludes the easy establishment of evolutionary relationships, the program permits us to find out important structural or functional residues upon the superimposition of two protein structures apparently unrelated. The derived classification can be very useful for finding relationships among proteins which would escape detection by current sequence or topology-based analytical algorithms.

Modeling the third loop of short-chain snake venom neurotoxins: roles of the short-range and long-range interactions

The influence of long-range interactions on local structures is an important issue in understanding protein folding process and protein structure stability. Using short-chain snake venom neurotoxin as a model system, we have studied the conformational properties of eight different loop III sequences either in the environment of one of the short-chain neurotoxin, erabutoxin b (PDB ID 1nxb), or in free state by Monte Carlo simulated annealing method. The surrounding protein structure was found to be crucial in stabilizing the loop conformation. Although all the eight peptides prefer type V beta turn in solution, three of them (KPGI, KPGV, KSGI) turn to type II beta turn and the other five (KKGI, KKGV, KNGI, KQGI, and KRGV) are confined to more rigid type V beta turn conformation in the protein structure. Using flexible tetra-glycine-peptide to screen the backbone conformational space in the protein environment also validates the results. This study shows that long-range interactions do contribute to the stability and the types of conformation for a surface loop in protein, while short-range interactions may only provide candidate conformations, which then have to be filtered by the long-range interactions further.

Do structural deviations between toxins adopting the same fold reflect functional differences?

Three-finger proteins form a structurally related family of compounds that exhibit a great variety of biological properties. To address the question of the prediction of functional areas on their surfaces, we tentatively conferred the acetylcholinesterase inhibitory activity of fasciculins on a short-chain curaremimetic toxin. For this purpose, we assimilated the three-dimensional structure of fasciculin 2 with the one of toxin alpha. This comparison revealed that the tips of the first and second loops, together with the C terminus residue, deviated most. A first recombinant fasciculin/toxin alpha chimera was designed by transferring loop 1 in its entirety together with the tip of loop 2 of fasciculin 2 into the toxin alpha scaffold. A second chimera (rChII) was obtained by adding the point Asn-61 --> Tyr substitution. Comparison of functional and structural properties of both chimeras show that rChII can accommodate the imposed modifications and displays nearly all the acetylcholinesterase-blocking activities of fasciculins. The three-dimensional structure of rChII demonstrates that rChII adopts a typical three-fingered fold with structural features of both parent toxins. Taken together, these results emphasize the great structural flexibility and functional adaptability of that fold and confirm that structural deviations between fasciculins and short-chain neurotoxins do indeed reflect functional diversity.

High resolution x-ray analysis of two mutants of a curaremimetic snake toxin

A previous mutational analysis of erabutoxin a (Ea), a curaremimetic toxin from sea snake venom, showed that the substitutions S8G and S8T caused, respectively, 176-fold and 780-fold affinity decreases for the nicotinic acetylcholine receptor (AchR). In view of the fact that the side-chain of Ser8 is buried in the wild-type toxin, we wondered whether these affinity changes reflect a direct binding contribution of S8 to the receptor and/or conformational changes that could have occurred in Ea as a result of the introduced mutations. To approach this question, we solved X-ray structures of the two mutants S8G and S8T at high resolution (0.18 nm and 0.17 nm, with R factors of 18.0% and 17.9%, respectively). The data show that none of the mutations significantly modified the toxin structure. Even within the site where the toxin binds to the receptor the backbone conformation remained unchanged. Therefore, the low affinities of the mutants S8T and S8G cannot be explained by a large conformational change of the toxin structure. Although we cannot exclude the possibility that undetectable structural changes have occurred in the toxin mutants, our data support the view that, although buried between loop I and II, S8 is part of the functional epitope of the toxin.

Comparison of three classes of snake neurotoxins by homology modeling and computer simulation graphics

We present a systematic structure comparison of three major classes of postsynaptic snake toxins, which include short and long chain alpha-type neurotoxins plus one angusticeps-type toxin of black mamba snake family. Two novel alpha-type neurotoxins isolated from Taiwan cobra (Naja naja atra) possessing distinct primary sequences and different postsynaptic neurotoxicities were taken as exemplars for short and long chain neurotoxins and compared with the major lethal short-chain neurotoxin in the same venom, i.e., cobrotoxin, based on the derived three-dimensional structure of this toxin in solution by NMR spectroscopy. A structure comparison among these two alpha-neurotoxins and angusticeps-type toxin (denoted as FS2) was carried out by the secondary-structure prediction together with computer homology-modeling based on multiple sequence alignment of their primary sequences and established NMR structures of cobrotoxin and FS2. It is of interest to find that upon pairwise superpositions of these modeled three-dimensional polypeptide chains, distinct differences in the overall peptide flexibility and interior microenvironment between these toxins can be detected along the three constituting polypeptide loops, which may reflect some intrinsic differences in the surface hydrophobicity of several hydrophobic peptide segments present on the surface loops of these toxin molecules as revealed by hydropathy profiles. Construction of a phylogenetic tree for these structurally related and functionally distinct toxins corroborates that all long and short toxins present in diverse snake families are evolutionarily related to each other, supposedly derived from an ancestral polypeptide by gene duplication and subsequent mutational substitutions leading to divergence of multiple three-loop toxin peptides.

Protein similarities beyond disulphide bridge topology

Structural superimposition is an important procedure to analyse the relationships between proteins. A new approach and program, KNOT-MATCH, has been developed for automated structural superimposition of proteins by means of their disulphide bridge topology. As a result of the superimposition, regular secondary structures, loops and clusters of residues become correctly aligned. This fact allows us to find out important structural overlaps of residues, sometimes with functional significance, not only among proteins belonging to the same family but also between apparently non-related proteins. Different disulphide-rich protein families, such as EGF-like, defensin-like and plant protease inhibitors, have been self or cross analysed with this approach. Some amino acids that have been experimentally determined to be structural and/or functional key residues for these proteins are conserved in the three-dimensional space after superimposition by KNOT-MATCH. The program can be very useful for finding relationships among proteins that would be hidden to the current alignment methods based on sequence and on main-chain topology.

Proline brackets and identification of potential functional sites in proteins: toxins to therapeutics

Protein toxins induce their specific pharmacological effects through protein protein interaction. Identification of these protein-protein interaction sites could lead to prototypes of highly specific therapeutic agents. However, deciphering the structure function relationships of protein toxins and locating the functional sites is a difficult, tedious and cumbersome task. We recently developed a novel predictive method to identify potential protein protein interaction sites directly from the amino acid sequence of a protein (R. M. Kini and H. J. Evans (1996) FEBS Lett. 385, 81-86) based on the presence of proline residues, a common residue found predominantly in the flanking segments of protein-protein interaction sites (R. M. Kini and H. J. Evans (1996) Biochem. Biophys. Res. Commun. 212, 1115-1124). It is a simple and straight-forward method. This review describes the new method and its application to solve structure function relationships of protein toxins. The method is useful in identifying functional sites in toxins, despite the subtle and complex nature of their structure function relationships and saves significant amounts of time and money.

Flanking proline residues identify the L-type Ca2+ channel binding site of calciseptine and FS2

Calciseptine and FS2 are 60-amino acid polypeptides, isolated from venom of the black mamba (Dendroaspis polylepis polylepis), that block voltage-dependent L-type Ca2+ channels. We predicted that these polypeptides contain an identical functional site between residues 43 and 46 by searching for proline residues that mark the flanks of protein-protein interaction sites [Kini, R. M., and Evans, H. J. (1966) FEBS Lett. 385, 81-86]. The predicted Ca2+ channel binding site also occurs in closely related toxins, C10S2C2 and S4C8. Therefore, it is likely that these toxins also will block L-type Ca2+ channels. To test the proposed binding site on calciseptine and FS2, an eight-residue peptide, named L-calchin (L-type calcium channel inhibitor), was synthesized and examined for biological activity. As expected for an L-type Ca2+ channel blocker, L-calchin reduced peak systolic and developed pressure in isolated rat heart Langendorff preparations without affecting diastolic pressure or heart rate. Furthermore, L-calchin caused a voltage-independent block of L-type Ca2+ channel currents in whole-cell patch-clamped rabbit ventricular myocytes. Thus the synthetic peptide exhibits the L-type Ca2+ channel blocking properties of the parent molecules, calciseptine and FS2, but with a lower potency. These results strongly support the identification of a site in calciseptine and FS2 that is important for binding to L-type Ca2+ channels and reinforce the importance of proline brackets flanking protein-protein interaction sites.

Snake venom cardiotoxins-structure, dynamics, function and folding

Snake cardiotoxins are highly basic (pI > 10) small molecular weight (approximately 6.5 kDa), all beta-sheet proteins. They exhibit a broad spectrum of interesting biological activities. The secondary structural elements in these toxins include antiparallel double and triple stranded beta-sheets. The three dimensional structures of these toxins reveal an unique asymmetric distribution of the hydrophobic and hydrophilic amino acids. The 3D structures of closely related snake venom toxins such as neurotoxins and cardiotoxin-like basic proteins (CLBP) fail to show similar pattern(s) in the distribution of polar and nonpolar residues. Recently, many novel biological activities have been reported for cardiotoxins. However, to-date, there is no clear structure-function correlation(s) available for snake venom cardiotoxins. The aim of this comprehensive review is to summarize and critically evaluate the progress in research on the structure, dynamics, function and folding aspects of snake venom cardiotoxins.

Solution structure of toxin b, a long neurotoxin from the venom of the king cobra (Ophiophagus hannah)

The solution structure of toxin b, a long neurotoxin (73 amino acids and 5 disulfides) from the venom of Ophiophagus hannah (king cobra), has been determined using 1H NMR and dynamical simulated annealing techniques. The structures were calculated using 485 distance constraints and 52 dihedral angle restraints. The 21 structures that were obtained satisfy the experimental restraints and possess good nonbonded contacts. Analysis of the converged structures revealed that the protein consists of a core region from which three finger-like loops extend outwards. The regular secondary structure in toxin b includes a double and a triple stranded antiparallel beta sheet. Comparison with the solution structures of other long neurotoxins reveals that although the structure of toxin b is similar to those of previously reported long neurotoxins, clear local structural differences are observed in regions proposed to be involved in binding to the acetylcholine receptor. A positively charged cluster is found in the C-terminal tail, in Loop III, and in the tip of Loop II. This cationic cluster could be crucial for the binding of the long neurotoxins to the acetylcholine receptor.

The role of acidic amino acid residues in the structural stability of snake cardiotoxins

We have recently shown that membrane-related activities of cardiotoxin V from Naja naja atra (CTX A5) are diminished at acidic pH although the overall beta-sheet structure of the molecule is maintained. In order to understand more about the mechanism of inactivation of CTX at acidic pH, we studied the effect of pH and denaturing reagents on the structural stability of CTX. We found, first, pH-induced structural transitions occurred in CTX A5 at two pH values as judged by the CD ellipticity around 195 nm: an increase in the beta-sheet content occurred around pH 4 and followed by a decrease, therein, around pH 2. The pKa of three acidic amino acid residues in CTX A5, i.e., Glu-17, Asp-42, and Asp-59, were determined to be 4.0, 3.2, and below 2.3, respectively, by NMR spectroscopy. The low pKa value of Asp-59 implies salt bridge formation between Lys-2 and Asp-59. Thus, electrostatic interaction may stabilize the three loop structure in addition to the hydrogen bonds between N- and C-termini of CTX molecule. Second, 2,2,2-trifluoroethanol (TFE) and guanidinium chloride (GdmHCI) were found to induce alpha-helical and random coil formation, respectively, in CTX A5 and eight other beta-sheet CTXs. Comparison of the relative potencies of TFE and GdmHCI to induce structural changes suggests that the amino acid residue located at position 17 plays a role in the structural stability. Specifically, CTXs containing negatively charged Glu-17 are least stable. It is suggested that Glu-17 may perturb the interaction between Lys-2 and Asp-59, and thus the overall stability of beta-sheet, in the presence of denaturing reagent. In conclusion, the perturbed structural stability of CTXs may partially explain the lower activity CTX exhibits at acidic pH. A structural model to account for the unfolding and refolding of CTX molecules without the breaking of disulfide bonds is also proposed.

Biology and treatment of the mamba snakebite

Mambas are venomous African snakes that are capable of inflicting fatal envenomation. The mambas (genus Dendrosaspis) are members of the family Elapidae [1]. Four species of mamba inhabit equatorial and southern Africa. Without medical treatment, mamba bites are frequently fatal. First aid treatment includes lymphatic retardation with immobilization and pressure wrap. Medical management requires the intravenous administration of mamba-specific antivenin.

Amino acid residue: is it structural or functional?

A new approach is suggested for delineating the structural and functional amino acid residues in proteins with known three-dimensional structure, basing on the involvement of residues in intramolecular hydrophobic and hydrophilic interactions and additional information about the conservativity of the residues. The approach is applied to the families of homologous neurotoxins and cardiotoxins. The results obtained concerning the role of amino acid residues in both families of toxins accord well with the similarity of their fold, but different mechanisms of action. Current approach can be used for detailed characterization of protein spatial structures, as well as for rational protein engineering.

Acetylcholinesterase inhibition by fasciculin: crystal structure of the complex

The crystal structure of the snake toxin fasciculin, bound to mouse acetylcholinesterase (mAChE), at 3.2 A resolution reveals a synergistic three-point anchorage consistent with the picomolar dissociation constant of the complex. Loop II of fasciculin contains a cluster of hydrophobic residues that interact with the peripheral anionic site of the enzyme and sterically occlude substrate access to the catalytic site. Loop I fits in a crevice near the lip of the gorge to maximize the surface area of contact of loop II at the gorge entry. The fasciculin core surrounds a protruding loop on the enzyme surface and stabilizes the whole assembly. Upon binding of fasciculin, subtle structural rearrangements of AChE occur that could explain the observed residual catalytic activity of the fasciculin-enzyme complex.

Detailed assessment of spatial hydrophobic and electrostatic properties of 2D NMR-derived models of neurotoxin II

2D NMR-derived spatial structures of neurotoxin II (NtII) and several homologous toxins in solution were assessed by comparison with their own amino acid sequences using a three-dimensional (3D) profile method. 3D profiles of all the toxin models match the sequences well and, therefore, the method of 3D profile was demonstrated to work correctly for these well-resolved NMR structures in aqueous solution. At the same time, the profile window plots reveal low scores in the bottom tip of loop II (residues 22-34 in NtII) and in beta-strand of loop III (residues 49-52). Some residues in the first poor-scoring region are of functional importance being involved in binding with nicotinic acetylcholine receptor (AChR). Furthermore, the second segment participates in intermolecular hydrogen bonding upon dimerization of postsynaptic neurotoxins in solution resulting in increasing of the 3D-1D score for residues at the interface between monomers. Therefore, the 3D profile method can be useful for detection functionally-important regions in well-resolved protein structures.

1H-NMR assignments and secondary structure of dendroaspin, an RGD-containing glycoprotein IIb-IIIa (alpha IIb-beta 3) antagonist with a neurotoxin fold

Dendroaspin, also referred to as mambin, was originally isolated from the venom of the Elapidae snake Dendroaspis jamesoni kaimose. It shares a high level of sequence similarity with the short-chain neurotoxins found in other Elapidae but displays approximately 1000-fold lower neurotoxin activity than the closely related protein erabutoxin b. However, unlike neurotoxins, it contains an RGD (Arg-Gly-Asp) motif and functions as an antagonist of platelet aggregation and cell-cell adhesion of comparable potency to the disintegrins from the venoms of Viperidae. We have determined the secondary structure of dendroaspin using 1H-NMR spectroscopy. Its structure resembles that of the short-chain neurotoxins, with three loops extending from a disulphide-bridged core; however, the strands of the triple-stranded beta-sheet are shorter and the loop containing the RGD sequence is moved away from this sheet. The structure bears little resemblance to that of the disintegrins, except in the RGD-containing loop, suggesting that this loop may be of prime importance in its inhibitory function. Comparison of this preliminary structure with that of the neurotoxins and disintegrins furthers our understanding of the mechanism of integrin antagonists and shows how the neurotoxin fold can be manipulated to give a variety of inhibitors.

Solution structure of omega-conotoxin GVIA using 2-D NMR spectroscopy and relaxation matrix analysis

We report here the solution structure of omega-conotoxin GVIA, a peptide antagonist of the N-type neuronal voltage-sensitive calcium channel. The structure was determined using two-dimensional NMR in combination with distance geometry and restrained molecular dynamics. The full relaxation matrix analysis program MARDIGRAS was used to generate maximum and minimum distance restraints from the crosspeak intensities in NOESY spectra. The 187 restraints obtained were used in conjunction with 23 angle restraints from vicinal coupling constants as input for the structure calculations. The backbones of the best 21 structures match with an average pairwise RMSD of 0.58 A. The structures contain a short segment of triple-stranded beta-sheet involving residues 6-8, 18-21, and 24-27, making this the smallest published peptide structure to contain a triple-stranded beta-sheet. Conotoxins have been shown to be effective neuroprotective agents in animal models of brain ischemia. Our results should aid in the design of novel nonpeptide compounds with potential therapeutic utility.

Determination of the nuclear-magnetic-resonance solution structure of cardiotoxin CTX IIb from Naja mossambica mossambica

The NMR structure of cardiotoxin CTX IIb from Naja mossambica mossambica in aqueous solution was determined from a total of 593 nuclear Overhauser enhancement distance constraints and 135 dihedral angle constraints, which were collected using two-dimensional homonuclear 1H-NMR experiments. Structure calculations were performed with the program DIANA, using the redundant dihedral angle constraints strategy for improved convergence, followed by restrained energy minimization with the program FANTOM and a modified version of the program AMBER. The CTX IIb structure is represented by a group of 20 conformers with an average root-mean-square deviation relative to the mean solution structure of 0.072 nm for the backbone atoms, and 0.116 nm for all heavy atoms. The molecular structure of CTX IIb is characterized by a three-stranded beta-sheet made up of residues 20-26, 32-39 and 48-54, and a two-stranded beta-sheet composed of residues 1-5 and 10-14. A cluster of four disulfide bonds, 3-21, 14-38, 42-53 and 54-59, form the core of the molecule and crosslink the individual polypeptide strands. The NMR structure is similar to the previously reported X-ray crystal structure of the cardiotoxin CTX VII4 from the same species. Differences between the two structures were noted in the tips of the two loops formed by residues 6-9 and 27-31, which connect the beta-strand 1-5 with 10-14, and 20-26 with 32-39, respectively. For these loops the NMR data also indicate significantly increased dynamic disorder in the solution structure. These observations are discussed with respect to earlier suggestions by others that these two loops are essential structural elements for function and specificity of a wide variety of homologous toxins.