Inventing an arsenal: adaptive evolution and neofunctionalization of snake venom phospholipase A2 genes
© Lynch; licensee BioMed Central Ltd. 2007
Received: 02 October 2006
Accepted: 18 January 2007
Published: 18 January 2007
Gene duplication followed by functional divergence has long been hypothesized to be the main source of molecular novelty. Convincing examples of neofunctionalization, however, remain rare. Snake venom phospholipase A2 genes are members of large multigene families with many diverse functions, thus they are excellent models to study the emergence of novel functions after gene duplications.
Here, I show that positive Darwinian selection and neofunctionalization is common in snake venom phospholipase A2 genes. The pattern of gene duplication and positive selection indicates that adaptive molecular evolution occurs immediately after duplication events as novel functions emerge and continues as gene families diversify and are refined. Surprisingly, adaptive evolution of group-I phospholipases in elapids is also associated with speciation events, suggesting adaptation of the phospholipase arsenal to novel prey species after niche shifts. Mapping the location of sites under positive selection onto the crystal structure of phospholipase A2 identified regions evolving under diversifying selection are located on the molecular surface and are likely protein-protein interactions sites essential for toxin functions.
These data show that increases in genomic complexity (through gene duplications) can lead to phenotypic complexity (venom composition) and that positive Darwinian selection is a common evolutionary force in snake venoms. Finally, regions identified under selection on the surface of phospholipase A2 enzymes are potential candidate sites for structure based antivenin design.
Phospholipase A2s (PLA2; EC 184.108.40.206) are esterolytic enzymes that hydrolyze glycerophospholipids at the sn-2 fatty acyl bond, releasing lysophospholipids and fatty acids. PLA2s play key roles in various biological processes in mammals including signal transduction, lipid digestion, host defense and production of eicosanoids and other lysophospholipid derivates with potent biological activities . PLA2 enzymes are also the major components of snake venoms where they function to immobilize and rapidly kill prey . PLA2s from elapid venoms (group-I) are structurally similar to pancreatic secretions while PLA2s from viper venoms (group-II) are structurally similar to inflammatory secretions . A third group of PLA2s (group-III) have been identified from the venom of bees, jellyfish, scorpions and lizards  indicating that PLA2s have been recruited into a toxic function multiple times in diverse lineages.
Snake venom PLA2s are members of large multigene families with diverse pharmacological activities including neurotoxic, myotoxic, cardiotoxic, anticoagulant and hemolytic effects . These diverse activities evolved from an ancestral nontoxic PLA2 by a process of repeated gene duplication followed by functional divergence. PLA2 toxicity is independent of enzymatic activity [3, 4] and is mediated through "pharmacological sites" on the protein surface that directly interact with ligands on the cell membrane . Thus, the surface of PLA2s forms a scaffold for adaptive modification that has been used to generate a diverse array of pharmacological effects through a process of neofunctionalization (the generation of new protein functions that were not the primary function of the ancestral protein).
Previous studies of PLA2 genes identified accelerated evolution of group-I genes from Naja naja  and group-II genes from Trimeresurus  and Vipera  consistent with positive Darwinian selection, however these studies focused on one or two species, included relatively few genes and used methods that lack power to detect episodic adaptive evolution. A larger study of group-I and -II genes found that amino acid substitutions were correlated with surface accessibility , suggesting that modifications of surface residues and positive selection play important roles in generating toxin diversity. To further explore this possibility I compiled an extensive dataset of full length snake venom PLA2 genes from public databases, inferred the gene trees for these toxins and tested for episodes of positive Darwinian selection coincident with the origin of novel pharmacological effects and recurrent diversifying selection on specific sites. In addition, I used a larger amino acids dataset to test if conclusions drawn from the smaller nucleotide dataset were robust to phylogenetic inference.
The results indicate that adaptive evolution is common in snake venom PLA2 genes and is associated with the evolution of new toxin functions and speciation events, demonstrating that molecular adaptation has played a pervasive role in the evolution of snakes and their venom arsenal. This analysis has identified the mutational pathway leading from non-toxic to highly toxic PLA2 enzymes, reconstructing the processes of mutation and adaptation. Finally, increases in genomic complexity gained through gene duplications has promoted the evolution an increasingly complex phenotype (venom composition), providing a link between molecular, phenotypic and organismal evolution.
Results and discussion
Gene duplication and speciation history
Even though branch support for the group-I and group-II nucleotide trees was high in this analysis, nucleotide data are only about a third of the gene sequences that are available, the majority are amino acid data (and thus not suitable for codon-based selection analysis discussed later). To test if the topology of nucleotide-based gene trees was sensitive to taxon sampling I inferring group-I and group-II phylogenies using larger amino acids datasets (245 group-I and 271 group-II genes, respectively). Although the trees inferred from amino acid data (Figs. 3 and 4) had lower support for recent lineages than the nucleotide data, perhaps because protein sequences have not accumulated enough phylogenetically informative substitutions to accurately reconstruct recent branching orders, the deeper nodes were well supported and the overall topology was congruent between amino acid and nucleotide data indicating that inferences based on the nucleotide datasets are reliable.
Origin and elaboration of toxic genes
Snake venom PLA2s have diverse pharmacological activities including neurotoxic, myotoxic, cardiotoxic, anticoagulant and hemolytic effects , which must have originated after they diverged from nontoxic ancestors. Although uncertainty in the tree topology and the richness of toxin functions makes assessment of specific ancestral and derived functions difficult for all lineages, it is clear from these phylogenies that many novel functions have originated in PLA2 genes after gene duplications. For example, three nontoxic group-I PLA2 genes isolated from Laticuadata semifasciata pancreas  branch near the base of the Elapinae group in the nucleotide tree, but are the most basal snake group-I genes in the amino acids tree (Figs. 1 and 3). These pancreatic genes have been suggested to be intermediates between nontoxic and toxic enzymes , suggesting that duplication of an ancestral nontoxic gene originally expressed in the pancreas was followed by recruitment into the venom gland and the emergence of toxic functions. Also within group-I, a clade containing neurotoxins from Laticuadata has emerged from antiplatelet enzymes; the nested position of this clade indicates that the neurotoxic effect of these enzymes is derived from more ancestral antiplatelet enzymes. The origin of group-II toxins that target muscle is also associated with a gene duplication event. Group-II myotoxins share a unique amino acid substitution at residue 49 (aspartate to lysine) that abolishes enzymatic activity [13, 14]. Thus, Lys49-myotoxins evolved a novel nonhydrolytic mechanism to induce membrane damage [15, 16] after duplication of an ancestral Asp49-PLA2 that was not myotoxic.
I used maximum likelihood models of coding-sequence evolution [17, 18] to test the hypothesis that functional diversification of snake venom PLA2 genes is driven by positive Darwinian selection. This method determines the strength and direction of selection by estimating the nonsynonymous-to-synonymous substitution rate (d N /d S = ω), with ω = 1, <1, and >1 indicating neutral evolution, purifying selection and directional selection, respectively. The branch-specific one-ratio model is the simplest, estimating the same ω for all branches in the phylogeny. The estimate of ω for group-I genes under this model, 1.28, is an average over all codons and lineages, highlighting the dominant role of positive selection on elapid venom phospholipases. The estimate of ω for group-II genes under the one-ratio model, 0.686, indicates that group-II genes are generally under purifying selection, however, this estimate is higher than reported from most genes.
The one-ratio model can only detect adaptive evolution when the majority of amino acids and lineages under study have been under positive selection (such as in group-I genes). If adaptive evolution is primarily episodic, then short episodes of positive selection, which are followed by long periods of purifying selection, will not be detected. To test for episodes of positive selection in group-I and group-II gene lineages, I used a free-ratios model that estimates separate d N /d S ratios for all lineages in the tree. These models fit the data significantly better than either the one-ratio model or a constrained one-ratio model with ω forced to be 1 (group-I genes) or a free-ratio model with lineages previously identified with ω >1 constrained to be 1 (group-II genes), indicating that episodes of directional selection are common in snake venom PLA2 evolution with nearly 32% and 21% of group-I and -II gene lineages, respectively, having been under directional selection (Figs. 1 and 2). Moreover, there are several branches with extremely high ω values, including two group-I and three group-II branches with ω >3, one group-I branch with ω >5 and one group-I branch with ω = 9.06 (Figs. 1 and 2).
Ohno's model  of post-duplication divergence predicts an increase in the nonsynonymous substitution rate following duplication as positive Darwinian selection drives the fixation of mutations that confer new or modified functions on gene duplicates. To test for accelerated evolution after duplication I used smaller datasets for which speciation and duplication events could be unambiguously assigned to each branch and a two-ratios model that estimated different ω parameters for post-duplication (PD) and post-speciation (PS) branches. Surprisingly, in group-I genes PD and PS branches have nearly identical ω values (ωPD = 1.12, ωPS = 1.22), indicating that positive selection is associated with both gene duplication and speciation. In contrast to group-I genes, group-II gene PD branches evolve much faster than PS branches (ωPD = 1.4, ωPS = 0.63), consistent with the classical model of neofunctionalization.
Here, neofunctionalization is defined as the emergence of a new toxic effect from an ancestral enzyme that did not posses that effect as its main toxin function (for example, neurotoxic Laticuadata genes and Lys49-myotoxins discussed above). Strikingly, positive selection occurred in the stem-lineage of 67% (4/6) of group-I functional groups and 88% (7/8) of group-II functional groups (Figs. 1 and 2) indicating that positive selection played a pervasive role in the origin of novel toxin functions during the diversification of vipers and elapids and their venoms. It also suggest lineages which can be targeted for ancestral sequence reconstructions for characterization of ancestral toxin functions to compare extant functions to.
The importance of gene duplication to the evolution of species-specific traits is relatively unknown, but duplications resulting in species-specific adaptations have been demonstrated for some genes [20, 21]. The unexpectedly high group-I ωPS may be the result of enzyme adaptation to new prey preference after speciation. Indeed, the three semi-aquatic kraits (Laticaudata sp.) prey primarily on moray and conger eels and assorted fishes [22, 23] while Australian copperheads (Austrelaps) prey on frogs and lizards . In the Elapinae group, the king cobra (Ophiophagus hannah) and kraits (Bungarus sp.) feed almost exclusively on snakes and other reptiles , while the true cobras (Naja sp.) and the Eastern brown snake (Pseudonaja textilis) feed on small mammals, amphibians and birds . This pattern suggest a scenario where dietary shifts after speciation runs the PLA2 gene repertoire through a "selective sieve"; those genes which are no longer effective in subduing new prey species are lost, while genes that are still effective adapt to the new prey type and subsequently diversify.
Maximum Likelihood Parameter Estimates for Group-I PLA2 Genes.
M0: one ratio...
= ω 0
P < 0.01
M0: one ratio-C
ω 0 constrained to 1
M0: one ratio-2
= ω 0
ω PD = 1.22, ω PS = 1.12
see Figure 1
P << 0.001
p0 = 0.148, ω0 = 0
p1 = 0.582, ω2 = 1
p0 = 0.145, ω0 = 0
P << 0.001
36 (PP ≥ 0.99)
p1 = 0.5, ω2 = 1
4 (0.95 ≤ PP < 0.99)
p2 = 0.355, ω 2 = 4.53
10 (PP < 0.95)
p0 = 0.343, ω0 = 0.115
P << 0.001
74 (PP ≥ 0.99)
p1 = 0.435, ω 2 = 1.43
7 (0.95 ≤ PP < 0.99)
p2 = 0.222, ω 2 = 4.63
4 (PP < 0.95)
p = 0.266, q= 0.182
p0 = 0.70, p = 0.276, q = 0.228
P << 0.001
28 (PP ≥ 0.99)
p1 = 0.297, ω = 3.41
2 (0.95 ≤ PP < 0.99)
9 (PP < 0.95)
Maximum Likelihood Parameter Estimates for Group-II PLA2 Genes.
M0: one ratio...
M0: one ratio-2
ω PD = 1.41, ωPS = 0.63
P = 0.014
see Figure 1
P << 0.001
ω 0 constrained to 1
P << 0.01
p0 = 0.124, ω0 = 0
p1 = 0.876, ω2 = 1
p0 = 0.124, ω0 = 0
P << 0.001
30 (PP ≥ 0.99)
p1 = 0.605, ω2 = 1
2 (0.95 ≤ PP < 0.99)
p2 = 0.271, ω 2 = 3.81
9 (PP < 0.95)
p0 = 0.346, ω0 = 0.078
P << 0.001
21 (PP ≥ 0.99)
p1 = 0.42, ω2 = 0.83
5 (0.95 ≤ PP < 0.99)
p2 = 0.234, ω 2 = 2.55
8 (PP < 0.95)
p = 0.332, q = 0.293
p0 = 0.81, p = 0.352, q = 0.348
P << 0.001
15 (PP ≥ 0.99)
p1 = 0.19, ω = 2.58
7 (0.95 ≤ PP < 0.99)
6 (PP < 0.95)
The molecular evolution of group-I and group-II PLA2 genes, such as the birth and death like and "selective sieve" processes of gene duplication, divergence and loss are similar to evolution of other snake venom proteins, particularly the elapid three-finger toxins . Indeed, there is even evidence for species-specific toxin adaptation to prey type within the three-finger toxins and maintenance of a well-ordered tertiary structure  similar to that seen in PLA2 genes, suggesting that this mode of molecular evolution may be common in venom genes.
Kini and Evans  have proposed that 'target sites' on the surface of prey cells are recognized by 'pharmacological sites' on PLA2 enzymes. These protein-protein interactions determine PLA2 specificity by having complementary charges, hydrophobicities, and Van der Waals contact surfaces. This model, combined with the analyses above, suggest that entirely new functions originate after duplication through substitutions in pharmacological sites that alter binding specificities. Although most substitutions will likely disrupt binding specificity for the current target site, a few may create new interaction sites leading to the emergence of novel functions.
The extraordinary level of positive selection acting on snake venom phospholipase A2 genes indicates that adaptive molecular evolution plays an important role in the emergence of these novel functions, continues as functions are diversified and refined, and may contribute to niche differentiation after speciation. Interestingly, mapping sites under positive selection onto the structure of PLA2 enzymes has identified regions that are attractive candidates for structure-based drug design. These data also demonstrate that increases in genomic complexity gained through gene duplications has lead to an increase in phenotypic complexity (venom composition) and likely the ability of venomous snakes to adapt to new prey types.
Sequence alignment and phylogenetic reconstruction
PLA2 genes were obtained from public database (GenBank GI's for each gene are shown in Figs. 1 and 2 for nucleotide data and in Figs. 3 and 4 for amino acid data). Partial sequences, sequences with insertions-deletions that caused reading frame shifts and sequences with premature stop codons were excluded from analysis as likely pseudogenes in the nucleotide analysis. Group-I and -II protein sequences were aligned in Clustal W  and adjusted by eye using Se-Al. Bayesian phylogenetic analyses were performed using MrBayes v3.0 . Tree searches were run using four Markov chains for 3,000,000 generations saving every 100th tree and a codon-based GTR+Γ+I model of sequence evolution or a JTT+Γ+I model for amino acid data. Models are nucleotide and amino acid data were selected using ModelGenerator. After stationarity, the final 15000 trees were used to build a consensus tree. Each analysis was performed three times to ensure convergence of tree topologies.
Tests for selection and ancestral sequence reconstructions
I used codon-based maximum likelihood models implemented in CODEML in the PAML package of programs (version 3.14)  to test for lineages under positive selection using the one-ratio and free-ratios models; this package of programs was also used for ancestral sequence reconstructions. To test for differences in post-duplication (PD) and post-speciation (PS) branches I used a smaller dataset of group-I and -II genes that included at least five representatives of each species/pharmacological group and for which gene duplication and speciation events could accurately be assigned for each branch. Twice the log likelihood difference between models, 2Δℓ = 2(ℓ1-ℓ0), is compared to a χ2-distribution with the degrees of freedom equal to the number of parameter differences between the models to test whether the alternative model (free-ratio or PD-PS) fits the data significantly better than the null model (one-ratio). If a lineage has a d N /d S > 1 and the likelihood ratio test is significant, than the neutral model of evolution is violated and positive selection is suggested. To explicitly test for positive selection I used constrained models that fixed ω at 1.
I used three pairs of site-specific models  to identify specific amino acids under diversifying selection: M0 and M3 (discrete), M1 (neutral) and M2 (selection), and M7 (beta) and M8 (beta & ω). Model M0 estimates a single ω parameter for all sites and branches, while model M3 (discrete) estimates three independent ω parameters and the proportion of sites belonging to each ω-class directly from the data. Model M1 (neutral) assumes two classes of sites in the protein with the proportion of conserved sites (ω = 0) and neutral sites (ω = 1) estimated from the data, while model M2 (selection) adds a third site class with an additional ω estimated as a free parameter allowing for sites with ω > 1. Model M7 (beta) uses a beta distribution B(p, q) with ω restricted to the interval (0,1) while Model M8 (beta & ω) adds a site class with the ω and the proportion of sites with that ω estimated from the data, allowing for sites with ω > 1. Twice the log likelihood difference between the models is compared to the χ2 distribution and tests for variation in ω among sites. After maximum likelihood parameter estimates are calculated, the Bayes theorem is used to calculate the posterior probability of belonging to a site class, when a site is identified with ω > 1 than positive selection is indicated. Sites with posterior probabilities of > 0.5 are reported here. The approximate posterior mean ω for each site from model M3 with two site classes were mapped onto the crystal structures of group-I (PDB ID: 1A3D) and -II (PDB ID: 1OZ6) PLA2 proteins using ICM-Browser (available from http://www.molsoft.com) or Chimera. The three-dimensonal space filling structures were generated with Deep View – spdbv v3.7 .
I thank G. P. Wagner for helpful comments and discussion of the manuscript and Z. Yang for answering questions regarding the use of PAML. I also thank M. J. McCarron for careful reading of the manuscript and the anonymous reviewers for their comments. Support for this research was provided by departmental graduate student research funds provided to VJL.
- Dennis EA: Phospholipases. The Enzymes. Edited by: Boyer PD. 1983, New York, Academic Press, 16: 3rdGoogle Scholar
- Kini RM: Venom Phospholipase A2 enzymes: Structure, Function and Mechanism. 1997, Chichester, Wiley, 511-Google Scholar
- Kini RM, Evans HJ: The role of enzymatic activity in inhibition of the extrinsic tenase complex by phospholipase A2 isoenzymes from Naja nigricollis venom. Toxicon. 1995, 33: 1585-1590. 10.1016/0041-0101(95)00103-4.View ArticlePubMedGoogle Scholar
- Rufini S: The proliferative effects of ammodytin L from the venom of Vipera ammodytes on 208F rat fibroblasts in culture. Biochemical Journal. 1996, 320: 318-326.View ArticleGoogle Scholar
- Kini RM, Evans HJ: A model to explain the pharmacological effects of snake venom phospholipases A2. Toxicon. 1989, 27: 613-635. 10.1016/0041-0101(89)90013-5.View ArticlePubMedGoogle Scholar
- Chuman Y, Nobuhisa I, Ogawa T, Deshimaru M, Chijiwa T, Tan NH, Fukumaki Y, Shimohigashi Y, Ducancel F, Boulain JC: Regional and accelerated molecular evolution in group I snake venom gland phospholipase A2 isozymes. Toxicon. 2000, 38: 449-462. 10.1016/S0041-0101(99)00165-8.View ArticlePubMedGoogle Scholar
- Nakashima K, Nobuhisa I, Deshimaru M, Nakai M, Ogawa T, Shimohigashi Y, Fukumaki Y, Hattori M, Sakaki Y, Hattori S, Ohno M: Accelerated Evolution in the Protein-Coding Regions is Universal in Crotalinae Snake Venom Gland Phospholipase A2 Isozyme Genes. PNAS. 1995, 92: 5605-5609. 10.1073/pnas.92.12.5605.PubMed CentralView ArticlePubMedGoogle Scholar
- Kordis D, Bdolah A, Gubensek F: Positive Darwinian Selection inVipera palaestinaePhospholipase A2Genes Is Unexpectedly Limited to the Third Exon. Biochemical and Biophysical Research Communications. 1998, 251: 613-619. 10.1006/bbrc.1998.9528.View ArticlePubMedGoogle Scholar
- Kini RM: Accelerated evolution and molecular surface of venom phospholipase A2 enzymes. Journal of Molecular Evolution. 1999, 48: 125-132. 10.1007/PL00006450.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Slowinski JB, Knight A, Rooney AP: Inferring Species Trees from Gene Trees: A Phylogenetic Analysis of the Elapidae (Serpentes) Based on the Amino Acid Sequences of Venom Proteins. Molecular Phylogenetics and Evolution. 1997, 8: 349-362. 10.1006/mpev.1997.0434.View ArticlePubMedGoogle Scholar
- Fujimi TJ, Kariya Y, Tsuchiya T, Tamiya T: Nucleotide sequence of phospholipase A2 gene expressed in snake pancreas reveals the molecular evolution of toxic phospholipase A2 genes. Gene. 2002, 292: 225-231. 10.1016/S0378-1119(02)00682-0.View ArticlePubMedGoogle Scholar
- Maraganore JM, Heinrikson RL: The lysine-49 phospholipase A2 from the venom of Agkistrodon piscivorus piscivorus. Relation of structure and function to other phospholipases A2 [published erratum appears in J Biol Chem 1993 Mar 15;268(8):6064]. J Biol Chem. 1986, 261: 4797-4804.PubMedGoogle Scholar
- van den Bergh CJ: The role of aspartic acid-49 in the active site of phospholipase A2. A site-specific mutagenesis study of porcine pancreatic phospholipase A2 and the rationale of the enzymatic activity of [lysine49]phospholipase A2 from Agkistrodon piscivorus piscivorus' venom. European Journal of Biochemistry. 1998, 176: 353-357. 10.1111/j.1432-1033.1988.tb14288.x.View ArticleGoogle Scholar
- Diaz C: The effect of myotoxins isolated from Bothrops snake venoms on multilamellar liposomes: relationship to phospholipase A2, anticoagulant and myotoxic activities. Biochim Biophys Acta. 1991, 1070: 455-460. 10.1016/0005-2736(91)90086-N.View ArticlePubMedGoogle Scholar
- Rufini S: Calcium ion independent membrane leakage induced by phospholipase-like myotoxins. Biochemistry. 1992, 31: 12424-12430. 10.1021/bi00164a018.View ArticlePubMedGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11: 725-736.PubMedGoogle Scholar
- Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15: 568-573.View ArticlePubMedGoogle Scholar
- Ohno S: Evolution by Gene Duplication. 1970, Berlin, SpringerView ArticleGoogle Scholar
- Zhang J: Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nature Genetics. 2002, 30: 411-415. 10.1038/ng852.View ArticlePubMedGoogle Scholar
- Riehle MM, Bennett AF, Long AD: Genetic architecture of thermal adaptation in Escherichia coli. PNAS. 2001, 98: 525-530. 10.1073/pnas.021448998.PubMed CentralView ArticlePubMedGoogle Scholar
- Shine R: Relationships between sexual dimophism and niche partitioning within a clade of sea snakes (Laticaudinae). Oecologia. 2002, 133: 45-53. 10.1007/s00442-002-1012-7.View ArticlePubMedGoogle Scholar
- Su Y: Food Habits of the Sea Snake, Laticauda semifasciata. Zoological Studies. 2005, 44: 403-408.Google Scholar
- Shine R: Habitats, diets, and sympatry in snakes: a study from Australia. Canadian Journal of Zoology. 1977, 55: 1118-1128.View ArticleGoogle Scholar
- Oriov N: Venomous snakes of southern China. Reptilia. 1992, 31: 22-29.Google Scholar
- Yang Z: Statistical Methods for detecting molecular adaptation. Trends in Ecology and Evoltuion. 2000, 15: 1994-1997.Google Scholar
- Kini RM: Excitement ahead: structure, function and mechanism of snake venom phospholipase A2 enzymes. Toxicon. 2003, 42: 827-840. 10.1016/j.toxicon.2003.11.002.View ArticlePubMedGoogle Scholar
- Chioato L, Ward RJ: Mapping structural determinants of biological activities in snake venom phospholipases A2 by sequence analysis and site directed mutagenesis. Toxicon. 2003, 42: 869-883. 10.1016/j.toxicon.2003.11.027.View ArticlePubMedGoogle Scholar
- Soares AM, Giglio JR: Chemical modifications of phospholipases A2 from snake venoms: effects on catalytic and pharmacological properties. Toxicon. 2003, 42: 855-868. 10.1016/j.toxicon.2003.11.004.View ArticlePubMedGoogle Scholar
- Fry BG WW: Molecular evolution and phylogeny of elapid snake venom three-finger toxins. Journal of Molecular Evolution. 2003, 57: 110-129. 10.1007/s00239-003-2461-2.View ArticlePubMedGoogle Scholar
- Higgins D: CLUSTAL W: improving the sensitivity of progressivemultiple sequence alignment through sequence weighting,position-specific gap penalties and weight matrix choice. Nucleic Acids Research. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Z: PAML. a program package for phylogenetic analysis by maximum likelihood. CABIOS. 1997, 13: 555-556.PubMedGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AMK: Codon-Substitution Models for Heterogeneous Selection Pressure at Amino Acid Sites. Genetics. 2000, 155: 431-449.PubMed CentralPubMedGoogle Scholar
- Guex N: SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modeling. Electrophoresis. 1997, 18: 2714-2723. 10.1002/elps.1150181505.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.