Studies on sea snake venom.

Erabutoxins a and b are neurotoxins isolated from venom of a sea snake Laticauda semifasciata (erabu-umihebi). Amino acid sequences of the toxins indicated that the toxins are members of a superfamily consisting of short and long neurotoxins and cytotoxins found in sea snakes and terrestrial snakes. The short neurotoxins to which erabutoxins belong act by blocking the nicotinic acetylcholine receptor on the post synaptic membrane in a manner similar to that of curare. X-ray crystallography and NMR analyses showed that the toxins have a three-finger structure, in which three fingers made of three loops emerging from a dense core make a gently concave surface of the protein. The sequence comparison and the location of essential residues on the protein suggested the mechanism of binding of the toxin to the acetylcholine receptor. Classification of snakes by means of sequence comparison and that based on different morphological features were inconsistent, which led the authors to propose a hypothesis "Evolution without divergence."

Tertiary structure of erabutoxin b in aqueous solution as elucidated by two-dimensional nuclear magnetic resonance.

The three-dimensional structure of erabutoxin b, a short-chain neurotoxic peptide purified from the venom of the sea snake Laticauda semifasciata, was determined in aqueous solution by two-dimensional proton nuclear magnetic resonance and simulated annealing-based calculations. On the basis of 883 assigned nuclear Overhauser effect (NOE) connectivities, 676 final distance constraints were derived and used together with 38 torsion angle (phi, chi 1) constraints, four distance constraints derived from disulfide bridges and 30 distance constraints derived from hydrogen bonds. A total of 14 converged structures were obtained from 50 runs of calculations. The atomic root-mean-square difference about the mean coordinate positions (excluding the residues 18 to 22) is 0.60 A for backbone atoms (N, C alpha and C'). The protein consists of a core region from which three finger-like loops emerge outwards. It includes a short, two-stranded antiparallel beta-sheet of residues 2 to 5 and 13 to 16, a three-stranded antiparallel beta-sheet involving residues 23 to 30, 35 to 41 and 50 to 56, and four disulfide bridges in the core region. Comparison with two crystal structures of erabutoxin b at 1.4 A and 1.7 A resolution indicated that the solution and the crystal structures were very similar, but less defined regions were observed at the localized region of the tip of the central loop and the outside of the third loop in solution. Other short-chain alpha-neurotoxins showed structural characteristics similar to those of erabutoxin b.

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

The three-dimensional structure in solution of the alpha-neurotoxin from the black mamba (Dendroaspis polylepis polylepis) has been determined by nuclear magnetic resonance spectroscopy. A high quality structure for this 60-residue protein was obtained from 656 NOE distance constraints and 143 dihedral angle constraints, using the distance geometry program DIANA for the structure calculation and AMBER for restrained energy minimization. For a group of 20 conformers used to represent the solution structure, the average root-mean-square deviation value calculated for the polypeptide backbone heavy atoms relative to the mean structure was 0.45 A. The protein consists of a core region from which three finger-like loops extend outwards. It includes a short, two-stranded antiparallel beta-sheet of residues 1-5 and 13-17, a three-stranded antiparallel beta-sheet involving residues 23-31, 34-42 and 51-55, and four disulfide bridges in the core region. There is also extensive non-regular hydrogen bonding between the carboxy-terminal tail of the polypeptide chain and the rest of the core region. Comparison with the crystal structure of erabutoxin-b indicates that the structure of alpha-neurotoxin is quite similar to other neurotoxin structures, but that local structural differences are seen in regions thought to be important for binding of neurotoxins to the acetylcholine receptor. For two regions of the alpha-neurotoxin structure there is evidence for an equilibrium between multiple conformations, which might be related to conformational rearrangements upon binding to the receptor. Overall, the alpha-neurotoxin presents itself as a protein with a stable core and flexible surface areas that interact with the acetylcholine receptor in such a way that high affinity binding is achieved by conformational rearrangements of the deformable regions of the neurotoxin structure.

Theoretical analysis of the structure of the peptide fasciculin and its docking to acetylcholinesterase.

The fasciculins are a family of closely related peptides that are isolated from the venom of mambas and exert their toxic action by inhibiting acetylcholinesterase (AChE). Fasciculins belong to the structural family of three-fingered toxins from Elapidae snake venoms, which include the alpha-neurotoxins that block the nicotinic acetylcholine receptor and the cardiotoxins that interact with cell membranes. The features unique to the known primary and tertiary structures of the fasciculin molecule were analyzed. Loop I contains an arginine at position 11, which is found only in the fasciculins and could form a pivotal anchoring point to AChE. Loop II contains five cationic residues near its tip, which are partly charge-compensated by anionic side chains in loop III. By contrast, the other three-fingered toxins show full charge compensation within loop II. The interaction of fasciculin with the recognition site on acetylcholinesterase was investigated by estimating a precollision orientation followed by determination of the buried surface area of the most probable complexes formed, the electrostatic field contours, and the detailed topography of the interaction surface. This approach has led to testable models for the orientation and site of bound fasciculin.

Three-dimensional crystal structure of recombinant erabutoxin a at 2.0 A resolution.

Recombinant erabutoxin a (Ea(r)) has been crystallized by vapour diffusion in hanging drops. The crystals belong to space group P2(1)2(1)2(1) with cell dimensions a = 55.8 A, b = 53.4 A, c = 40.8 A. Diffraction data have been recorded on a FAST detector up to 2.0 A. The atomic crystal structure of Ea(r) has been determined by initial refinement of the structure of the isotoxin erabutoxin b (Eb) the crystals of which were grown under identical conditions. The R-factor was 23% at 2.0 A resolution. The secondary and tertiary structures of Ea(r) are shown to be identical with that of wild-type Eb, within the experimental error.

Functional determinants by which snake and cone snail toxins block the alpha 7 neuronal nicotinic acetylcholine receptors.

Snakes and cone snails produce toxins which block muscular and/or neuronal nicotinic acetylcholine receptors (AChRs). This paper mostly focuses on the determinants by which a snake long chain curaremimetic toxin and the cone snail toxin ImI bind specifically to the alpha 7 neuronal receptor. In both cases, the site involves a small turn-like structure constrained by two half-cystines.

Mechanism of protein folding. III. Disulfide bonding.

It was shown in lysozyme and phospholipase, and generally in proteins with disulfide bonds, that after the formation of secondary structures the hydrophobic interactions between the key pairs responsible for folding tertiary structures bring several cysteine residues close together. Among the possible combination of cysteine residues some definite pairs are realized in the tertiary structure. In the Appendix to this paper an algebraic relation is given which must be satisfied for two cysteine residues to make a disulfide bond. This relation is too strict to be applied to real problems, where the two cysteines come close together, but the distance is still too great to yield a disulfide bond. In this case the two residues can attract each other by disulfide formation potential. A geometrical graphic representation is given for the locus of the H atom of the SH group in the cysteine residue. This looks like a lampshade and provides us with a guide to select the correct choice among cysteine pairs. This method is applied to lysozyme and phospholipase to supplement the discussion of the preceding paper (T. Yoshimura, H. Noguchi, T. Inoue and N. Saitô, Biophys. Chem. 40 (1991) 277).

Comparison of refolding patterns of erabutoxin b and cardiotoxin 3.10.2 from snake venom.

The refolding patterns of erabutoxin b (a neurotoxin) and cardiotoxin 3.10.2 (from Naja naja siamensis venom) have been studied by reducing both the proteins by treatment with reduced dithiothreitol followed by renaturation by treatment with oxidised dithiothreitol. Isoelectric focusing of the samples trapped at varying time intervals during renaturation of the proteins reveals formation of intermediates in the folding pathway with cardiotoxin 3.10.2. having fewer intermediates than erabutoxin b and faster rate of refolding (1 hr and 3 hr respectively).

Structure of dimeric and monomeric erabutoxin a refined at 1.5 A resolution.

Erabutoxin a has been crystallized in its monomeric and dimeric forms. The structures were refined at 1.50 and 1.49 A resolution, respectively, using synchrotron radiation data. The crystals belong to space group P212121, with cell dimensions a = 49.84, b = 46.62, c = 21.22 A for the monomer and a = 55.32, b = 53.54, c = 40.76 A for the dimer. Using starting models from earlier structure determinations, the monomeric structure refined to an R value of 16.7% (8004 unique reflections, 17.0-1.50 A resolution range), while the dimeric structure has been solved by the molecular-replacement method with a final R value of 16.9% (19 444 unique reflections, 17.4-1.49 A resolution range). The high-resolution electron-density maps clearly revealed significant discrete disorder in the proteins and allowed an accurate determination of the solvent structure. For the monomer, the side chains of six residues were modelled with alternate conformers and 106 sites for water molecules and one site for a sulfate ion were included in the final model, whereas for the dimer, 206 sites for water molecules were included and both C-terminal residues together with the side chains of 11 residues adopted alternative conformations. A comparison was made with earlier structure determinations. The features of the solvent structure of the erabutoxin molecules are discussed in detail.

Multistep modeling of protein structure: application to bungarotoxin.

Modelling of bungarotoxin in atomic details is presented in this article. The model-building procedure utilizes the low-resolution crystal coordinates of the c-alpha atoms of bungarotoxin, sequence homology within the neurotoxin family, as well as high-resolution x-ray diffraction data of cobratoxin and erabutoxin. Our model-building procedure involves: (a) principles of comparative modelling, (b) embedding procedures of distance geometry, and (c) use of molecular mechanics for optimizing packing. The model is not only consistent with the c-alpha coordinates of crystal structure, but also agrees with solution conformational features of the triple-stranded beta sheet as observed by NOE measurements.

Three-dimensional structures of proteins determined by two-dimensional NMR and distance geometry calculations.

Due to the recent development of NMR spectroscopy, the three-dimensional structures of proteins up to 10 Kd in aqueous solution can be determined in atomic levels. In the present article, we show the intuitive description on the theoretical background of structural determination by NMR, where we take erabutoxin b as an example and compare the structures in aqueous solution and in the crystalline state. Then, we show the three-dimensional structure of mouse epidermal growth factor and human transforming growth factor alpha determined by NMR. On the basis of the results, we discuss the molecular mechanism of recognition between the epidermal growth factor and its receptor.

Solution structure of neuronal bungarotoxin determined by two-dimensional NMR spectroscopy: calculation of tertiary structure using systematic homologous model building, dynamical simulated annealing, and restrained molecular dynamics.

Neuronal bungarotoxin has previously been shown, using two-dimensional 1H NMR spectroscopy, to have a triple-stranded antiparallel beta-sheet structure which dimerizes in solution [Oswald, R.E., Sutcliffe, M.J., Bamberger, M., Loring, R.H., Braswell, E., & Dobson, C.M. (1991) Biochemistry 30, 4901-4909]. In this paper, structural calculations are described which use the 582 experimentally measured NOE restraints in conjunction with 27 phi-angle restraints from J-value measurements. The positions of the N-terminal region and C-terminal region were poorly defined in the calculated structures with respect to the remainder of the structure. The region of the structure containing the triple-stranded beta-sheet was, however, well defined and similar to that found in the structure of homologous alpha-bungarotoxin (45% amino acid identity). The experimental restraints did not result in a well-defined dimer interface region because of the small number of NOEs which could be identified in this region. An approach was therefore adopted which produced model structures based to varying degrees on the alpha-bungarotoxin structure. Fourteen different structures were generated in this manner and subsequently used as starting points for refinement using dynamical simulated annealing followed by restrained molecular dynamics. This approach, which combines NMR data and homologous model building, has enabled a family of structures to be proposed for the dimeric molecule. In particular, Phe49 has been identified as possibly playing an important role in dimer formation, this residue in one chain interacting with the corresponding residue in the adjacent chain.

Solution conformation of cobrotoxin: a nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing study.

The solution conformation of cobrotoxin has been determined by using proton nuclear magnetic resonance spectroscopy. With the combination of various two-dimensional NMR techniques, the 1H-NMR spectrum of cobrotoxin was completely assigned (Yu et al., 1990). A set of 435 approximate interproton distance restraints was derived from nuclear Overhauser enhancement (NOE) measurements. These NOE constraints, in addition to the 29 dihedral angle constraints (from coupling constant measurements) and 26 hydrogen bonding restraints (from the pattern of short-range NOEs), form the basis of 3-D structure determination by the hybrid distance geometry-dynamical simulated annealing method. The 23 structures that were obtained satisfy the experimental restraints, display small deviation from idealized covalent geometry, and possess good nonbonded contacts. Analysis of converged structures indicated that there are two antiparallel beta sheets (double and triple stranded), duly confirming our earlier observations. These are well defined in terms of both atomic root mean square (RMS) differences and backbone torsional angles. The average backbone RMS deviation between the calculated structures and the mean structure, for the beta-sheet regions, is 0.92 A. The mean solution structure was compared with the X-ray crystal structure of erabutoxin b, the homologous protein. This yielded information that both structures resemble each other except at the exposed loop/surface regions, where the solution structure seems to possess more flexibility.

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.

Structure determination of a dimeric form of erabutoxin-b, crystallized from a thiocyanate solution.

Erabutoxin-b, M(r) = 6861.1, a single 62 amino-acid chain folded by four disulfide bridges, was crystallized in a new orthorhombic form by using thiocyanate as crystallizing agent. The space group is P2(1)2(1)2(1) with a = 53.36 (4), b = 40.89 (4), c = 55.71 (5) A, V = 121533.1 A and Z = 8. X-ray diffraction data were recorded at the LURE synchrotron facility (lambda = 1.405 A). The structure was solved by molecular replacement and shows a dimeric association through an anti-parallel beta-sheet around the twofold non-crystallographic axis. The two independent molecules, one SCN- ion and 97 associated water molecules were refined by molecular dynamics and annealing techniques to R = 19.6% (10,913 Fobs, resolution 5-1.7 A). The thiocyanate ion is located at the interface of the dimer and close to the non-crystallographic twofold axis.

The length of a single turn controls the overall folding rate of "three-fingered" snake toxins.

Snake curaremimetic toxins are short all-beta proteins, containing several disulfide bonds which largely contribute to their stability. The four disulfides present in snake toxins make a "disulfide beta-cross"-fold that was suggested to be a good protein folding template. Previous studies on the refolding of snake toxins (Ménez, A. et al. (1980) Biochemistry 19, 4166-4172) showed that this set of natural homologous proteins displays different rates of refolding. These studies suggested that the observed different rates could be correlated to the length of turn 2, one out of five turns present in the toxins structure and close to the "disulfide beta-cross". To demonstrate this hypothesis, we studied the refolding pathways and kinetics of two natural isotoxins, toxin alpha (Naja nigricollis) and erabutoxin b (Laticauda semifasciata), and two synthetic homologues, the alpha mutants, alpha60 and alpha62. These mutants were designed to probe the peculiar role of the turn 2 on the refolding process by deletion or insertion of one residue in the turn length that reproduced the natural heterogeneity at that locus. The refolding was studied by electrospray mass spectrometry (ESMS) time-course analysis. This analysis permitted both the identification and quantitation of the population of intermediates present during the process. All toxins were shown to share the same sequential scheme for disulfide bond formation despite large differences in their refolding rates. The results presented here demonstrate definitely that no residues except those forming turn 2 accounted for the observed differences in the refolding rate of toxins.

Genetic engineering of snake toxins. The functional site of Erabutoxin a, as delineated by site-directed mutagenesis, includes variant residues.

Using site-directed mutagenesis, we previously identified some residues that probably belong to the site by which Erabutoxin a (Ea), a sea snake toxin, recognizes the nicotinic acetylcholine receptor (AcChoR) (Pillet, L., Trémeau, O., Ducancel, F. Drevet, P., Zinn-Justin, S., Pinkasfeld, S., Boulain, J.-C., and Ménez, A. (1993) J. Biol. Chem. 268, 909-916). We have now studied the effect of mutating 26 new positions on the affinity of Ea for AcChoR. The mutations are F4A, N5V, H6A, Q7L, S9G, Q10A, P11N, Q12A, T13V, T14A, K15A, T16A, delta S18, E21A, Y25F, Q28A, S30A, T35A, I36R, P44V, T45A, V46A, K47A, P48Q, I50Q, and S53A. Binding affinity decreases upon mutation at Gln-7, Gln-10 and to a lesser extent at His-6, Ser-9 and Tyr-25 whereas it increases upon mutation at Ile-36. Other mutations have no effect on Ea affinity. In addition, new mutations of the previously explored Ser-8, Asp-31, Arg-33, and Glu-38 better explain the functional role of these residues in Ea. The previous and present mutational analysis suggest that the "functional" site of Ea covers a homogeneous surface of at least 680 A2, encompassing the three toxin loops, and includes both conserved and variant residues. The variable residues might contribute to the selectivity of Ea for some AcChoRs, including those from fish, the prey of sea snakes.

Three-dimensional solution structure of a curaremimetic toxin from Naja nigricollis venom: a proton NMR and molecular modeling study.

The solution conformation of toxin alpha from Naja nigricollis (61 amino acids and four disulfides), a snake toxin which specifically blocks the activity of the nicotinic acetylcholine receptor (AcChoR), has been determined using nuclear magnetic resonance spectroscopy and molecular modeling. The solution structures were calculated using 409 distance and 73 dihedral angle restraints. The average atomic rms deviation between the eight refined structures and the mean structure is approximately 0.5 A for the backbone atoms. The overall folding of toxin alpha consists of three major loops which are stabilized by three disulfide bridges and one short C terminal loop stabilized by a fourth disulfide bridge. All the disulfides are grouped in the same region of the molecule, forming a highly constrained structure from which the loops protrude. As predicted, this structure appears to be very similar to the 1.4-A resolution crystal structure of another snake neurotoxin, namely, erabutoxin b from Laticauda semifasciata. The atomic rms deviation for the backbone atoms between the solution and crystal structures is approximately 1.7 A. The minor differences which are observed between the two structures are partly related to the deletion of one residue from the chain of toxin alpha. It is notable that, although the two toxins differ from each other by 16 amino acid substitutions, their side chains have an essentially similar spatial organization. However, most of the side chains which constitute the presumed AcChoR binding site for the curaremimetic toxins are poorly resolved in toxin alpha.(ABSTRACT TRUNCATED AT 250 WORDS)

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.

Genomic structures of cardiotoxin 4 and cobrotoxin from Naja naja atra (Taiwan cobra).

Two genomic DNAs with the size of 2.3 kb and 2.4 kb, which were isolated from the liver of Naja naja atra (Taiwan cobra), encoded the precursors of cardiotoxin 4 and cobrotoxin, respectively. Both genes shared virtually identical overall organization with three exons separated by two introns, which were inserted in the similar positions of the gene's coding regions. Moreover, their nucleotide sequences shared approximately 84.2% identity. This result reveals the evolutionary relationship between cardiotoxin and cobrotoxin. The exon/intron structures of cardiotoxin 4 and cobrotoxin genes were similar to that reported for erabutoxin c gene, a neurotoxin genomic DNA from a sea snake (Laticauda semifasciata). However, in contrast to the finding that the intron 2 of these genes had a similar size, a notable variation with the size of intron 1 was observed (1233 bp, 1269 bp and 197 bp for cardiotoxin 4, cobrotoxin and erabutoxin c genes, respectively). The different size with intron 1 is due to the middle region at the first intron of cardiotoxin 4 and cobrotoxin genes, which encoded small nucleolar RNA (snoRNA), being absent in that of erabutoxin c gene. These results, together with the finding of the potential mobility of snoRNA genes during evolution, suggest that intron insertions or deletions of snoRNA genes occur with the evolutionary divergence of snake neurotoxins and cardiotoxins.