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.     

Computational analysis of evolutionary divergence of scorpion toxins

Conserved and consensus sequences of several members of scorpion toxin families have been analyzed. Multiple sequence alignment define the highly conserved residues that play important structural role in formation of toxin subgroups. Scorpion toxins have characteristic feature in signature pattern {(GA)}- K-C {(LIVM)}-X (2)- K- C- C-X-C) }. where lysine and cysteine residues are well conserved in the polypeptide. The toxins show two major groups, one with conserved active region GKCMNKGKC and second with CTPK. There is variability in position of lysine and arginine in different toxins. Most of the scorpion toxins contain six conserved cysteines involved in disulphide bonds. From the evolutionary analysis it is concluded that Leiurus is a primitive ancestral species from which probably various toxins have evolved during the evolution.     

Adaptive evolution of animal toxin multigene families

Animal toxins comprise a diverse array of proteins that have a variety of biochemical and pharmacological functions. A large number of animal toxins are encoded by multigene families. From studies of several toxin multigene families at the gene level the picture is emerging that most have been functionally diversified by gene duplication and adaptive evolution. The number of pharmacological activities in most toxin multigene families results from their adaptive evolution. The molecular evolution of animal toxins has been analysed in some multigene families, at both the intraspecies and interspecies levels. In most toxin multigene families, the rate of non-synonymous to synonymous substitutions (dN/dS) is higher than one. Thus natural selection has acted to diversify coding sequences and consequently the toxin functions. The selection pressure for the rapid adaptive evolution of animal toxins is the need for quick immobilization of the prey in classical predator and prey interactions. Currently available evidence for adaptive evolution in animal toxin multigene families will be considered in this review.     

Concerted evolution of sea anemone neurotoxin genes is revealed through analysis of the Nematostella vectensis genome

Gene families, which encode toxins, are found in many poisonous animals, yet there is limited understanding of their evolution at the nucleotide level. The release of the genome draft sequence for the sea anemone Nematostella vectensis enabled a comprehensive study of a gene family whose neurotoxin products affect voltage-gated sodium channels. All gene family members are clustered in a highly repetitive approximately 30-kb genomic region and encode a single toxin, Nv1. These genes exhibit extreme conservation at the nucleotide level which cannot be explained by purifying selection. This conservation greatly differs from the toxin gene families of other animals (e.g., snakes, scorpions, and cone snails), whose evolution was driven by diversifying selection, thereby generating a high degree of genetic diversity. The low nucleotide diversity at the Nv1 genes is reminiscent of that reported for DNA encoding ribosomal RNA (rDNA) and 2 hsp70 genes from Drosophila, which have evolved via concerted evolution. This evolutionary pattern was experimentally demonstrated in yeast rDNA and was shown to involve unequal crossing-over. Through sequence analysis of toxin genes from multiple N. vectensis populations and 2 other anemone species, Anemonia viridis and Actinia equina, we observed that the toxin genes for each sea anemone species are more similar to one another than to those of other species, suggesting they evolved by manner of concerted evolution. Furthermore, in 2 of the species (A. viridis and A. equina) we found genes that evolved under diversifying selection, suggesting that concerted evolution and accelerated evolution may occur simultaneously.     

Diversity of bacterial type II toxin-antitoxin systems: a comprehensive search and functional analysis of novel families

Type II toxin-antitoxin (TA) systems are generally composed of two genes organized in an operon, encoding a labile antitoxin and a stable toxin. They were first discovered on plasmids where they contribute to plasmid stability by a phenomenon denoted as 'addiction', and subsequently in bacterial chromosomes. To discover novel families of antitoxins and toxins, we developed a bioinformatics approach based on the 'guilt by association' principle. Extensive experimental validation in Escherichia coli of predicted antitoxins and toxins increased significantly the number of validated systems and defined novel toxin and antitoxin families. Our data suggest that toxin families as well as antitoxin families originate from distinct ancestors that were assembled multiple times during evolution. Toxin and antitoxin families found on plasmids tend to be promiscuous and widespread, indicating that TA systems move through horizontal gene transfer. We propose that due to their addictive properties, TA systems are likely to be maintained in chromosomes even though they do not necessarily confer an advantage to their bacterial hosts. Therefore, addiction might play a major role in the evolutionary success of TA systems both on mobile genetic elements and in bacterial chromosomes.     

Diverse evolutionary mechanisms shape the type III effector virulence factor repertoire in the plant pathogen Pseudomonas syringae

Many gram-negative pathogenic bacteria directly translocate effector proteins into eukaryotic host cells via type III delivery systems. Type III effector proteins are determinants of virulence on susceptible plant hosts; they are also the proteins that trigger specific disease resistance in resistant plant hosts. Evolution of type III effectors is dominated by competing forces: the likely requirement for conservation of virulence function, the avoidance of host defenses, and possible adaptation to new hosts. To understand the evolutionary history of type III effectors in Pseudomonas syringae, we searched for homologs to 44 known or candidate P. syringae type III effectors and two effector chaperones. We examined 24 gene families for distribution among bacterial species, amino acid sequence diversity, and features indicative of horizontal transfer. We assessed the role of diversifying and purifying selection in the evolution of these gene families. While some P. syringae type III effectors were acquired recently, others have evolved predominantly by descent. The majority of codons in most of these genes were subjected to purifying selection, suggesting selective pressure to maintain presumed virulence function. However, members of 7 families had domains subject to diversifying selection.     

Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.     

Novel genes encoding six kinds of three-finger toxins in Ophiophagus hannah (king cobra) and function characterization of two recombinant long-chain neurotoxins

Three-finger toxins are a family of low-molecular-mass toxins (<10 kDa) having very similar three-dimensional structures. In the present study, 19 novel cDNAs coding three-finger toxins were cloned from the venom gland of Ophiophagus hannah (king cobra). Alignment analysis showed that the putative peptides could be divided into six kinds of three-finger toxins: LNTXs (long-chain neurotoxins), short-chain neurotoxins, cardiotoxins (CTXs), weak neurotoxins, muscarinic toxins and a toxin with a free SH group. Furthermore, a phylogenetic tree was established on the basis of the toxin cDNAs and the previously reported similar nucleotide sequences from the same source venom. It indicated that three-finger-toxin genes in O. hannah diverged early in the course of evolution by long- and short-type pathways. Two LNTXs, namely rLNTX1 (recombinant LNTX1) and rLNTX3, were expressed and showed cytolytic activity in addition to their neurotoxic function. By comparing the functional residues, we offer some possible explanations for the differences in their neurotoxic function. Moreover, a plausible elucidation of the additonal cytolytic activity was achieved by hydropathy-profile analysis. This, to our knowledge, is the first observation that recombinant long chain alpha-neurotoxins have a CTX-like cytolytic activity.     

Analysis of active site residues of botulinum neurotoxin E by mutational, functional, and structural studies: Glu335Gln is an apoenzyme

Clostridial neurotoxins comprising the seven serotypes of botulinum neurotoxins and tetanus neurotoxin are the most potent toxins known to humans. Their potency coupled with their specificity and selectivity underscores the importance in understanding their mechanism of action in order to develop a strategy for designing counter measures against them. To develop an effective vaccine against the toxin, it is imperative to achieve an inactive form of the protein which preserves the overall conformation and immunogenicity. Inactive mutants can be achieved either by targeting active site residues or by modifying the surface charges farther away from the active site. The latter affects the long-range forces such as electrostatic potentials in a subtle way without disturbing the structural integrity of the toxin causing some drastic changes in the activity/environment. Here we report structural and biochemical analysis on several mutations on Clostridium botulinum neurotoxin type E light chain with at least two producing dramatic effects: Glu335Gln causes the toxin to transform into a persistent apoenzyme devoid of zinc, and Tyr350Ala has no hydrolytic activity. The structural analysis of several mutants has led to a better understanding of the catalytic mechanism of this family of proteins. The residues forming the S1' subsite have been identified by comparing this structure with a thermolysin-inhibitor complex structure.     

Domain loss facilitates accelerated evolution and neofunctionalization of duplicate snake venom metalloproteinase toxin genes

Gene duplication is a key mechanism for the adaptive evolution and neofunctionalization of gene families. Large multigene families often exhibit complex evolutionary histories as a result of frequent gene duplication acting in concordance with positive selection pressures. Alterations in the domain structure of genes, causing changes in the molecular scaffold of proteins, can also result in a complex evolutionary history and has been observed in functionally diverse multigene toxin families. Here, we investigate the role alterations in domain structure have on the tempo of evolution and neofunctionalization of multigene families using the snake venom metalloproteinases (SVMPs) as a model system. Our results reveal that the evolutionary history of viperid (Serpentes: Viperidae) SVMPs is repeatedly punctuated by domain loss, with the single loss of the cysteine-rich domain, facilitating the formation of P-II class SVMPs, occurring prior to the convergent loss of the disintegrin domain to form multiple P-I SVMP structures. Notably, the majority of phylogenetic branches where domain loss was inferred to have occurred exhibited highly significant evidence of positive selection in surface-exposed amino acid residues, resulting in the neofunctionalization of P-II and P-I SVMP classes. These results provide a valuable insight into the mechanisms by which complex gene families evolve and detail how the loss of domain structures can catalyze the accelerated evolution of novel gene paralogues. The ensuing generation of differing molecular scaffolds encoded by the same multigene family facilitates gene neofunctionalization while presenting an evolutionary advantage through the retention of multiple genes capable of encoding functionally distinct proteins.     

Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution

Efforts to describe toxins from the two major families of venomous snakes (Viperidae and Elapidae) usually reveal proteins belonging to few structural types, particular of each family. Here we carried on an effort to determine uncommon cDNAs that represent possible new toxins from Lachesis muta (Viperidae). In addition to nine classes of typical toxins, atypical molecules never observed in the hundreds of Viperidae snakes studied so far are highly expressed: a diverging C-type lectin that is related to Viperidae toxins but appears to be independently originated; an ohanin-like toxin, which would be the third member of the most recently described class of Elapidae toxins, related to human butyrophilin and B30.2 proteins; and a 3FTx-like toxin, a new member of the widely studied three-finger family of proteins, which includes major Elapidae neurotoxins and CD59 antigen. The presence of these common and uncommon molecules suggests that the repertoire of toxins could be more conserved between families than has been considered, and their features indicate a dynamic process of venom evolution through molecular mechanisms, such as multiple recruitments of important scaffolds and domain exchange between paralogs, always keeping a minimalist nature in most toxin structures in opposition to their nontoxin counterparts.     

一个新的东亚钳蝎毒素(BmKT_1)全长cDNA的克隆和分析

首先构建了东亚钳蝎毒腺组织 c DNA文库 ;根据已知的东亚钳蝎哺乳动物毒素氨基酸序列保守区设计引物 ,并用 PCR从 c DNA文库中扩增出一个 c DNA片段作为筛选 c DNA文库的探针 ;从 c DNA文库中筛选到二个编码同一个新的蝎毒素多肽的 c DNA,它们除 3′- UTR外 ,其余序列完全一致 .它们均含有 2 55bp长的开放阅读框 ,编码 85肽的前体毒素 ,包括 1 9个氨基酸残基的信号肽 ,66个残基的成熟毒素 (命名为 Bm KT1) ;Bm KT1氨基酸序列与已知的蝎毒素具有较大的同源性 ,与 Bm KM1,Lqq ,Lqhα IT和 Bm K M10 的同源性分别为 77%、67%、67%和 65% .Bm KT1的 C端不存在末端修饰步骤且具有一个与这些毒素不相同的特征结构 ,即在末端延伸了两个氨基酸残基 - P- S,推测 Bm KT1具有新的活性功能特征 .     

Domain organization and evolution of multifunctional autoprocessing repeats-in-toxin (MARTX) toxin in Vibrio vulnificus

The objective of this study was to analyze multifunctional autoprocessing repeats-in-toxin (MARTX) toxin domain organization within the aquatic species Vibrio vulnificus as well as to study the evolution of the rtxA1 gene. The species is subdivided into three biotypes that differ in host range and geographical distribution. We have found three different types (I, II, and III) of V. vulnificus MARTX (MARTX(Vv)) toxins with common domains (an autocatalytic cysteine protease domain [CPD], an α/β-hydrolase domain, and a domain resembling that of the LifA protein of Escherichia coli O127:H6 E2348/69 [Efa/LifA]) and specific domains (a Rho-GTPase inactivation domain [RID], a domain of unknown function [DUF], a domain resembling that of the rtxA protein of Photorhabdus asymbiotica [rtxA(PA)], and an actin cross-linking domain [ACD]). Biotype 1 isolates harbor MARTX(Vv) toxin types I and II, biotype 2 isolates carry MARTX(Vv) toxin type III, and biotype 3 isolates have MARTX(Vv) toxin type II. The analyzed biotype 2 isolates harbor two identical copies of rtxA1, one chromosomal and the other plasmidic. The evolutionary history of the gene demonstrates that MARTX(Vv) toxins are mosaics, comprising pieces with different evolutionary histories, some of which have been acquired by intra- or interspecific horizontal gene transfer. Finally, we have found evidence that the evolutionary history of the rtxA1 gene for biotype 2 differs totally from the gene history of biotypes 1 and 3.     

Evolutionary Analyses Suggest a Function of MxB Immunity Proteins Beyond Lentivirus Restriction

Viruses impose diverse and dynamic challenges on host defenses. Diversifying selection of codons and gene copy number variation are two hallmarks of genetic innovation in antiviral genes engaged in host-virus genetic conflicts. The myxovirus resistance (Mx) genes encode interferon-inducible GTPases that constitute a major arm of the cell-autonomous defense against viral infection. Unlike the broad antiviral activity of MxA, primate MxB was recently shown to specifically inhibit lentiviruses including HIV-1. We carried out detailed evolutionary analyses to investigate whether genetic conflict with lentiviruses has shaped MxB evolution in primates. We found strong evidence for diversifying selection in the MxB N-terminal tail, which contains molecular determinants of MxB anti-lentivirus specificity. However, we found no overlap between previously-mapped residues that dictate lentiviral restriction and those that have evolved under diversifying selection. Instead, our findings are consistent with MxB having a long-standing and important role in the interferon response to viral infection against a broader range of pathogens than is currently appreciated. Despite its critical role in host innate immunity, we also uncovered multiple functional losses of MxB during mammalian evolution, either by pseudogenization or by gene conversion from MxA genes. Thus, although the majority of mammalian genomes encode two Mx genes, this apparent stasis masks the dramatic effects that recombination and diversifying selection have played in shaping the evolutionary history of Mx genes. Discrepancies between our study and previous publications highlight the need to account for recombination in analyses of positive selection, as well as the importance of using sequence datasets with appropriate depth of divergence. Our study also illustrates that evolutionary analyses of antiviral gene families are critical towards understanding molecular principles that govern host-virus interactions and species-specific susceptibility to viral infection.     

Orthologs, paralogs and genome comparisons

During the past decade, ancient gene duplications were recognized as one of the main forces in the generation of diverse gene families and the creation of new functional capabilities. New tools developed to search data banks for homologous sequences, and an increased availability of reliable three-dimensional structural information led to the recognition that proteins with diverse functions can belong to the same superfamily. Analyses of the evolution of these superfamilies promises to provide insights into early evolution but are complicated by several important evolutionary processes. Horizontal transfer of genes can lead to a vertical spread of innovations among organisms, therefore finding a certain property in some descendants of an ancestor does not guarantee that it was present in that ancestor. Complete or partial gene conversion between duplicated genes can yield phylogenetic trees with several, apparently independent gene duplications, suggesting an often surprising parallelism in the evolution of independent lineages. Additionally, the breakup of domains within a protein and the fusion of domains into multifunctional proteins makes the delineation of superfamilies a task that remains difficult to automate.     

Revealing the paradox of drug reward in human evolution

Neurobiological models of drug abuse propose that drug use is initiated and maintained by rewarding feedback mechanisms. However, the most commonly used drugs are plant neurotoxins that evolved to punish, not reward, consumption by animal herbivores. Reward models therefore implicitly assume an evolutionary mismatch between recent drug-profligate environments and a relatively drug-free past in which a reward centre, incidentally vulnerable to neurotoxins, could evolve. By contrast, emerging insights from plant evolutionary ecology and the genetics of hepatic enzymes, particularly cytochrome P450, indicate that animal and hominid taxa have been exposed to plant toxins throughout their evolution. Specifically, evidence of conserved function, stabilizing selection, and population-specific selection of human cytochrome P450 genes indicate recent evolutionary exposure to plant toxins, including those that affect animal nervous systems. Thus, the human propensity to seek out and consume plant neurotoxins is a paradox with far-reaching implications for current drug-reward theory. We sketch some potential resolutions of the paradox, including the possibility that humans may have evolved to counter-exploit plant neurotoxins. Resolving the paradox of drug reward will require a synthesis of ecological and neurobiological perspectives of drug seeking and use.     

Insights into the evolutionary origins of clostridial neurotoxins from analysis of the <Emphasis Type="Italic">Clostridium botulinum</Emphasis> strain A neurotoxin gene cluster

Background  

Clostridial neurotoxins (CNTs) are the most deadly toxins known and causal agents of botulism and tetanus neuroparalytic diseases. Despite considerable progress in understanding CNT structure and function, the evolutionary origins of CNTs remain a mystery as they are unique to Clostridium and possess a sequence and structural architecture distinct from other protein families. Uncovering the origins of CNTs would be a significant contribution to our understanding of how pathogens evolve and generate novel toxin families.