Step-wise thermal denaturation of cobrotoxin, a snake venom neurotoxin fromNaja naja atra: A proton nuclear magnetic resonance study

Temperature dependence of proton nuclear magnetic resonance spectra has been followed for cobrotoxin, a postsynaptic neurotoxin fromNaja naja atra venom. Several aromatic amino-acid residues, including the functionally essential Trp-29 located at the tip of the central loop of the molecule, have been found to undergo a thermal structural transition above the global thermal denaturation temperature. It is suggested that a local structure around these residues behaves somehow independently of the rest of the molecule, and that such structural organization may be favorable for a conformational change of a neurotoxin molecule on binding to acetylcholine receptor.     

Effect of photooxidation on biological properties of botulinum neurotoxin A

Photooxidation of botulinic neurotoxin A in the presence of methylene blue is associated with a decrease in toxicity down to complete detoxication. During neurotoxin photooxidation, when the toxicity makes up to 1 to 3% of the original one, the conformation of the neurotoxin molecule and its antigenic properties remain unchanged. Under these conditions, using diethylpyrocarbonate, a specific reagent for histidine, the photooxidized neurotoxin was found to contain 5-6 oxidized histidine residues per molecule of neurotoxin; this was accompanied by changes in the UV absorbance spectrum around 280 nm. It was assumed that the main decrease in neurotoxin toxicity during photooxidation is probably due to oxidation of tryptophane, since the differential UV spectra suggest that the higher the extremum around 280 nm, the greater the decrease of toxicity; chemical modification of histidine residues alone causes no noticeable detoxication.     

Role of arginine and histidine residues in the biological activity of botulinic neurotoxin A

Treatment of botulinic neurotoxin A with cyclohexanedione demonstrated that modification of 5 to 10 arginine residues does not change the neurotoxin toxicity, while after modification of 15-20 arginine residues the toxicity is decreased by 40-50% of the original value. Butanedione exerts a stronger detoxicating effect on neurotoxin than cyclohexanedione. The molecular conformation of the modified toxin derivatives and their precipitability upon interaction with antisera against toxin and toxin fragments does not change thereby. The non-toxic derivatives of toxin containing 40 modified arginine residues possess a partial serological affinity for the original toxin in a reaction with antiserum against toxin but do not interact with the antifragment sera. The molecular conformation of these preparations is changed considerably. It is assumed that one or two arginine residues are located near the toxic site of the neurotoxin molecule and are also components of its antigenic determinants. Modification of histidine residues in the neurotoxin molecule by diethylpyrocarbonate is accompanied by a decrease of its toxicity. An additional 10% toxicity is revealed upon modification of 11-13 histidine residues. The molecular conformation of the modified derivatives of neurotoxin and their precipitability do not change thereby. It is probable that 1 or 2 histidine residues are located at or near the toxic site. The data obtained suggest that histidine residues are not localized in antigenic determinants of the neurotoxin molecule.     

Structurally homologous toxins isolated from the Taiwan cobra (Naja naja atra) differ significantly in their structural stability

Cardiotoxin and neurotoxin analogues isolated from snake venom sources are highly homologous proteins (>50% homology) with similar three-dimensional structures but exhibit drastically different biological properties. In the present study, we compare the conformational stability of cardiotoxin analogue III (CTX III) and cobrotoxin (CBTX), a neurotoxin analogue, from the Taiwan cobra (Naja naja atra), using circular dichroism spectroscopy and hydrogen-deuterium (H/D) exchange techniques in conjunction with two-dimensional NMR methods. Contrary to expectations, it is found that CTX III and CBTX differ significantly in their structural stabilities. The three-dimensional structure of CBTX is less stable than that of CTX III. The amide protons of residues at the N- and C-terminal ends of the CTX III molecule are strongly protected against H/D exchange, implying that the terminal ends of the molecule are bridged together by significant numbers of hydrogen bonds. However, in CBTX, amide protons at the terminal ends of the molecule do not exhibit an significant protection against H/D exchange. Comparison of the protection factors of the various amide protons in CTX III and CBTX reveals that the extraordinary stability of CTX III stems from the strong network of interactions among the residues at the N- and C-terminal ends and also due to the tight and ordered packing of the nonpolar residues involved in the triple-stranded, anti-parallel, beta-sheet segment of the molecule.     

The role of glycine residues at the C-terminal peptide segment in antinociceptive activity: a molecular dynamics simulation

Elucidating structural determinants in the functional regions of toxins can provide useful knowledge for designing novel analgesic peptides. Glycine residues at the C-terminal region of the neurotoxin BmK AGP-SYPU2 from the scorpion Buthus martensii Karsch (BmK) have been shown to be crucial to its analgesic activity. However, there has been no research on the structure–function relationship between the C-terminal segment of this toxin and its analgesic activity. To address this issue, we performed three MD simulations: one on the native structure and the other two on mutants of that structure. Results of these calculations suggest that the existence of glycine residues at the C-terminal segment stabilizes the protruding topology of the NC domain, which is considered an important determinant of the analgesic activity of BmK AGP-SYPU2.     

Structural and biochemical characterization of the therapeutic Anabaena variabilis phenylalanine ammonia lyase

We have recently observed promising success in a mouse model for treating the metabolic disorder phenylketonuria with phenylalanine ammonia lyase (PAL) from Rhodosporidium toruloides and Anabaena variabilis. Both molecules, however, required further optimization in order to overcome problems with protease susceptibility, thermal stability, and aggregation. Previously, we optimized PAL from R. toruloides, and in this case we reduced aggregation of the A. variabilis PAL by mutating two surface cysteine residues (C503 and C565) to serines. Additionally, we report the structural and biochemical characterization of the A. variabilis PAL C503S/C565S double mutant and carefully compare this molecule with the R. toruloides engineered PAL molecule. Unlike previously published PAL structures, significant electron density is observed for the two active-site loops in the A. variabilis C503S/C565S double mutant, yielding a complete view of the active site. Docking studies and N-hydroxysuccinimide-biotin binding studies support a proposed mechanism in which the amino group of the phenylalanine substrate is attacked directly by the 4-methylidene-imidazole-5-one prosthetic group. We propose a helix-to-loop conformational switch in the helices flanking the inner active-site loop that regulates accessibility of the active site. Differences in loop stability among PAL homologs may explain the observed variation in enzyme efficiency, despite the highly conserved structure of the active site. A. variabilis C503S/C565S PAL is shown to be both more thermally stable and more resistant to proteolytic cleavage than R. toruloides PAL. Additional increases in thermal stability and protease resistance upon ligand binding may be due to enhanced interactions among the residues of the active site, possibly locking the active-site structure in place and stabilizing the tetramer. Examination of the A. variabilis C503S/C565S PAL structure, combined with analysis of its physical properties, provides a structural basis for further engineering of residues that could result in a better therapeutic molecule.     

Sequential1H-NMR assignments of neurotoxin III from the sea anemoneHeteractis macrodactylus and structural comparison with related toxins

The complete sequence-specific assignment of resonances in the¹H-NMR spectrum of the polypeptide neurotoxin III (Hm III) from the sea anemoneHeteractis macrodactylus is described. Comparison of the chemical shifts and pattern of NOEs for Hm III with those for the related toxin Hp III fromHeteractis paumotensis, which differs only in the substitution of Asn for Tyr at position 11, shows that the overall secondary and tertiary structures are conserved. The largest differences in chemical shift caused by the substitution at position 11 are observed for the NH resonances of Arg-13, Thr-14, Ala-15, Leu-17, and Cys-26. The C^αH resonances influenced most are those of ASP-6, Gly-9, Leu-17, and Glu-42, while the most affected C^βH resonances are from Leu-17, Glu-28, and Lys-32. The absence of long-range NOEs to the aromatic ring of Tyr-11 as well as the lack of significant chemical shift effects on residues outside the loop comprising residues 7–16 confirm that this part of the loop makes no long-lived contacts with the rest of the molecule. The deviations from random coil shifts of Hm III are compared with those of the related anemone toxins Hp II, Hp III, and toxin I fromStichodactyla helianthus (Sh I). The similarity in deviations in chemical shift as a function of sequence position for these four toxins emphasizes the overall structural homology among these polypeptides.     

Characterization of the subunits of botulinum neurotoxin type A

A comparative amino acid analysis of botulinum neurotoxin type A and its subunits has been carried out. The heavy and light chains of neurotoxin have the same ratios of polar and non-polar amino acids (1.3:1), the amount of tryptophan residues in the heavy chain is 4 times as much as that in the light chain, and the number of SH-groups exceeds that in the light chains 2-fold. In neurotoxin, two N-terminal amino acid residues--alanine and leucine--were identified. Alanine was found to be the N-terminus of the heavy chain. The fluorescence spectra of neurotoxin subunits indicate differences in the conformational state of the polypeptide chains. The antigenic non-identity of botulinum neurotoxin A subunits suggests the presence in the neurotoxin molecule of at least two antigenic determinants, corresponding to the heavy and light chains.     

Structure and refinement of penicillopepsin at 1.8 A resolution

Penicillopepsin, the aspartyl protease from the mould Penicillium janthinellum, has had its molecular structure refined by a restrained-parameter least-squares procedure at 1.8 Å resolution to a conventional R-factor of 0.136. The estimated co-ordinate accuracy for the majority of the 2363 atoms of the enzyme is better than 0.12 Å. The average atomic thermal vibration parameter, B, for the atoms of the enzyme is 14.5 Ų. One determining factor of this low average B value is the large central hydrophobic core, in which there are two prominent clusters of aromatic residues, one of nine, the other of seven residues. The N and C-terminal domains of penicillopepsin display an approximate 2-fold symmetry: 70 residue pairs are topologically equivalent, related by a rotation of 177 ° and a translation of 1.2 Å. The analysis of the secondary structural features of the molecule reveals non-linear hydrogen bonding. In penicillopepsin, there is no difference in the mean hydrogen-bond parameters for the elements of α-helix, parallel or antiparallel β-pleated sheet. The mean values for these structural elements are: NO, 2.90 Å; NHO, 1.95 Å; N?O, 160 °. The average hydrogen-bond parameters of the reverse β-turns and the 3₁₀ helices are distinctly different from the above values. The analysis of sidechain conformational angles χ¹ and χ² penicillopepsin and other enzyme structures refined in this laboratory shows much narrower distributions as compared with those compiled from unrefined protein structures. The close proximity of the carboxyl groups of Asp33 and Asp213 suggests that they share a proton in a tight hydrogen-bonded environment (Asp33OD2 to Asp213OD1 is 2.87 Å). There are several solvent molecules in the active site region and, in particular, O39 forms hydrogen-bonded interactions with both aspartate residues. The disposition of the two carboxyl groups suggests that neither is likely to be involved in a direct nucleophilic attack on the scissile bond of a substrate. The average atomic B-factors of the residues in this region of the molecule are between 5 and 8 Ų, confirming the proposal that conformational mobility of the active site residues has no role in the enzymatic mechanism. However, conformational mobility of neighbouring regions of the molecule e.g. the “flap” containing Tyr75, is verified by the high B-factors for those residues. The positions of 319 solvent sites per asymmetric unit have been selected from difference electron density maps and refined. Thirteen have been classified as internal, and several of these may have key roles during catalysis. The positively charged N^ζ atom of Lys304 forms hydrogen bonds to the carboxylate of Asp14 (internal ion pair) and to two internal water molecules O5 and O25. The protonated side-chain of Asp300 forms a hydrogen bond to Thr214O, 2.78 Å, and is the recipient of a hydrogen bond from a surface pocket water molecule O46. There is no possibility for direct interaction between Asp300 and Lys304 without large conformational changes of their environment. The intermolecular packing involves many protein-protein contacts (66 residues) with a large number of solvent molecules involved in bridging between polar residues at the contact surface. The penicillopepsin molecules resemble an approximate hexagonal close-packing of spheres with each molecule having 12 “nearest” neighbours.     

Spatial organization of catalase proteins

The three-dimensional structure of the heme-containing fungal catalase fromPenicillium vitale (m.m. 2,80,000) has been studied by X-ray analysis at 2.0 A resolution. The molecule is tetramer, each subunit contains 670 aminoacid residues identified to construct “X-ray” primary structure. The subunit is built of three compact domains and their connections. The first domain of about 350 residues contains aβ-barrel flanked by helices, the second domain of 70 residues is formed by four helices and the third one is composed of 150 residues and is topologically similar to flavodoxin. The active site including heme is deeply buried near theβ-barrel. A comparison of the structure of catalase fromPenicillium vitale with that of beef liver catalase revealed very close structural homology of the first and the second domain, but the third domain is entirely absent in beef liver catalase. A catalase from thermophillic bacteriaThermus thermophilus (m.m. 2,10,000) has been first isolated, crystallized and studied by X-ray analysis. Crystals are cubic, space group is P2₁3, a = 133.4 Å. The molecule is a hexamer with trigonal symmetry 32. The electron density map at 3 Å resolution made it possible to trace the polypeptide chain. The main structural motif is formed by four near parallel helices. There is no heme inThermus thermophilus catalase, the active site is between the four helices and contains two manganese ions.     

Structure and biological activity of hagfish insulin

An isomorphously phased electron density map of hagfish (Myxine glutinosa) insulin has been calculated at a resolution of 3·1 Å spacing. The molecule crystallises with one molecule per asymmetric unit but is organised as a symmetric dimer lying on a 2-fold crystal axis. The structure of the hagfish insulin monomer is much more similar to that of pig insulin molecule 2 than molecule 1 of the dimer that constitutes one third of the 2 Zn insulin hexamer. There are different conformations however at the N and C termini of the B-chain. At the C terminus the two final residues on hagfish insulin partially obscure the A1 glycine residue, which in pig insulin is exposed. This structural difference has been shown, however, not to be responsible for the reduced activity of the hagfish insulin.     

Structure and biological activity of botulinum neurotoxin

Botulinum neurotoxin appears to undergo structural alterations after synthesis and also before it inhibits neurotransmitter release. Discussions and conjectures are presented in this context: 1. At what sites on the 150 kDa neurotoxin does posttranslational proteolytic processing occur? 2. Does neurotransmitter inhibition depend on separation of a segment of the neurotoxin from the rest of the molecule? 3. At what step in the intoxication pathway does activation of neurotoxin (enhanced lethality following limited proteolysis) manifest? 4. Can the receptor binding parameters (based on bovine brain synaptosome and lipid membrane), channel forming property (lipid bilayer membrane) and intracellular inhibitory activity (based on permeabilized chromaffin and PC 12 cells) provide clues to differences in the lethal potency between the neurotoxin serotypes? In addition, the following issues are considered: 5. The spontaneous fragmentation of isolated 50 kDa light chain, after its separation from 100 kDa heavy chain, 6. Effect of specific chemical modification of Arg, His, Lys, Trp, Tyr and Asp/Glu residues of types A, B and E neurotoxins on lethality and antigenicity, and 7. Development of second generation toxoids.     

Conformational changes associated with the nicking and activation of botulinum neurotoxin type E

Secondary and tertiary structural parameters of type E botulinum neurotoxin in the unactivated single-chain and activated two-chain (i.e., after proteolytic cleavage) forms were analyzed using circular dichroism, derivative absorption and fluorescence spectroscopy. The estimated secondary structures (22 and 20% alpha-helix, 44 and 44% beta-pleated sheets, and 34 and 36% random coils for the single- and two-chain neurotoxins, respectively) indicated that virtually no change occurred upon nicking of the single-chain neurotoxin. About 57% of the 70 Tyr residues were exposed in the single-chain form, which increased to 62% in the two-chain form. Fluorescence quenching experiments with neutral, anionic and cationic quenchers indicated that about 40% of the maximum accessible fluorescent Trp residues were exposed on the surface of the single-chain neurotoxin as compared to only 20% in the case of the two-chain neurotoxin. Acrylamide was the most effective quencher with a fraction accessibility of 0.56 and 0.48 of maximum accessible Trp fluorescence residues in the single and two-chain forms of the neurotoxin, respectively. Native polyacrylamide gel electrophoresis of the two forms of the neurotoxin revealed greater mobility for the two chain form. This indicates that the surface charges in the single-chain neurotoxin were altered upon nicking. These observations suggest that nicking of the single-chain type E neurotoxin results in refolding and redistribution of the surface charges of the neurotoxin.     

Aptamer Selection Express: A Novel Method for Rapid Single-Step Selection and Sensing of Aptamers

Here we describe a new DNA capture element (DCE) sensing system, based on the quenching and dequenching of a double-stranded aptamer. This system shows very good sensitivity and thermal stability. While quenching, dequenching, and separating the DCE systems made from different aptamers (all selected by SELEX), an alternative method to rapidly select aptamers was developed—the Aptamer Selection Express (ASExp). This process has been used to select aptamers against different types of targets (Bacillus anthracis spores, Bacillus thuringiensis spores, MS-2 bacteriophage, ovalbumin, and botulinum neurotoxin). The DCE systems made from botulinum neurotoxin aptamers selected by ASExp have been investigated. The results of this investigation indicate that ASExp can be used to rapidly select aptamers for the DCE sensing system.     

Statistics and Physical Origins of pK and Ionization State Changes upon Protein-Ligand Binding

This work investigates statistical prevalence and overall physical origins of changes in charge states of receptor proteins upon ligand binding. These changes are explored as a function of the ligand type (small molecule, protein, and nucleic acid), and distance from the binding region. Standard continuum solvent methodology is used to compute, on an equal footing, pK changes upon ligand binding for a total of 5899 ionizable residues in 20 protein-protein, 20 protein-small molecule, and 20 protein-nucleic acid high-resolution complexes. The size of the data set combined with an extensive error and sensitivity analysis allows us to make statistically justified and conservative conclusions: in 60% of all protein-small molecule, 90% of all protein-protein, and 85% of all protein-nucleic acid complexes there exists at least one ionizable residue that changes its charge state upon ligand binding at physiological conditions (pH = 6.5). Considering the most biologically relevant pH range of 4-8, the number of ionizable residues that experience substantial pK changes (ΔpK > 1.0) due to ligand binding is appreciable: on average, 6% of all ionizable residues in protein-small molecule complexes, 9% in protein-protein, and 12% in protein-nucleic acid complexes experience a substantial pK change upon ligand binding. These changes are safely above the statistical false-positive noise level. Most of the changes occur in the immediate binding interface region, where approximately one out of five ionizable residues experiences substantial pK change regardless of the ligand type. However, the physical origins of the change differ between the types: in protein-nucleic acid complexes, the pK values of interface residues are predominantly affected by electrostatic effects, whereas in protein-protein and protein-small molecule complexes, structural changes due to the induced-fit effect play an equally important role. In protein-protein and protein-nucleic acid complexes, there is a statistically significant number of substantial pK perturbations, mostly due to the induced-fit structural changes, in regions far from the binding interface.     

Spatial structure of apamin in solution]

The calculation of the complete spatial structure of the bee venom peptide neurotoxin apamin has been carried out by means of a method elaborated earlier. It is based on the joint utilization of the molecular mechanics algorithms and NMR spectroscopy data. It was established that the molecule backbone conformation in solution may be represented as the combination of the beta-turn III (residues 2-5) and alpha-helical segment (9-18) both separated by the non-standard bend IV (5-8). The most probable system of the intramolecular hydrogen bonds in the apamin polypeptide backbone was proposed. Certain amino acid residues have been shown to be characterized by the lack of strict determination of the conformations of their side chains which may be realized in a few states providing approximately equal stabilization of the same form of the main chain. The conformational parameters of the proposed apamin structural model are appropriate to the NMR spectroscopy data derived from the literature and used in the calculations and are not contradictory to other experimental information.     

The Unfolding MD Simulations of Cyclophilin: Analyzed by Surface Contact Networks and Their Associated Metrics

Currently, considerable interest exists with regard to the dissociation of close packed aminoacids within proteins, in the course of unfolding, which could result in either wet or dry moltenglobules. The progressive disjuncture of residues constituting the hydrophobic core ofcyclophilin from L. donovani (LdCyp) has been studied during the thermal unfolding of the molecule, by molecular dynamics simulations. LdCyp has been represented as a surface contactnetwork (SCN) based on the surface complementarity (S_m) of interacting residues within themolecular interior. The application of S_m to side chain packing within proteins make it a very sensitive indicator of subtle perturbations in packing, in the thermal unfolding of the protein. Network based metrics have been defined to track the sequential changes in the disintegration ofthe SCN spanning the hydrophobic core of LdCyp and these metrics prove to be highly sensitive compared to traditional metrics in indicating the increased conformational (and dynamical) flexibility in the network. These metrics have been applied to suggest criteria distinguishing DMG, WMG and transition state ensembles and to identify key residues involved in crucial conformational/topological events during the unfolding process.     

Purification and characterization of a neurotoxin from the venom of Ophiophagus hannah (king cobra)

A neurotoxin, Oh9-1, from the venom of Ophiophagus hannah was isolated by a combination of ion-exchange chromatography and reverse phase HPLC. Amino acid sequence analysis revealed that Oh9-1 consists of 57 amino acids and eight cysteine residues. This protein was mainly constituted with beta-sheet as evidenced by CD spectrum. Oh9-1 inhibited carbachol-induced muscle contraction in an irreversible manner and the dose for achieving 50% inhibition was approximately fourfold that of alpha-bungarotoxin. Since the residues in alpha-neurotoxins closely involve in the binding with acetylcholine receptors are not highly conserved in this toxin molecule, Oh9-1 represents a novel type of neurotoxin structurally distinct from alpha-neurotoxins.     

An approach for the assessment of the order of disruption of the elements of protein structure upon protein unfolding: A study of carbonic anhydrase B

An experimental approach named μ-analysis has been developed in order to elucidate the sequence of the loss of ordered structure by elements of a protein during the denaturation of the molecule. This approach is applicable for the analysis of proteins that fold (unfold) in a multistep process that involve the formation (destruction) of a range of intermediate states. The concept of the approach consists in systematic analysis of mutagenized forms of the protein with point substitutions of hydrophobic amino-acid residues and additional cysteine bridges. Importantly, the substitutions of the amino-acid residues must be localized to the same structural elements of the protein. Point substitutions of hydrophobic amino-acid residues mainly provide information on the structural elements of the protein that are disrupted at the final stages of protein denaturation. The addition of cysteine bridges to the surface of the protein molecule allows investigation of structural elements of the protein that are the first to unfold upon protein denaturation. Calorimetric studies of non-equilibrium melting of bovine carbonic anhydrase B yielded information on the rate constants of the unfolding of ten mutant forms of the protein. The analysis of the effects of mutations on the rates of different stages of protein unfolding allowed for elucidation of the order of disruption of structural elements of carbonic anhydrase B upon thermal denaturation.