Evolution of Snake Venom Disintegrins by Positive Darwinian Selection

Abstract

PII-disintegrins, cysteine-rich polypeptides broadly distributed in the venoms of geographically diverse species of vipers and rattlesnakes, antagonize the adhesive functions of beta(1) and beta(3) integrin receptors. PII-disintegrins evolved in Viperidae by neofunctionalization of disintegrin-like domains of duplicated PIII-snake venom hemorrhagic metalloproteinase (SVMP) genes recruited into the venom proteome before the radiation of the advanced snakes. Minimization of the gene (loss of introns and coding regions) and the protein structures (successive loss of disulfide bonds) underpins the postduplication divergence of disintegrins. However, little is known about the underlying genetic mechanisms that have generated the structural and functional diversity among disintegrins. Phylogenetic inference and maximum likelihood-based codon substitution approaches were used to analyze the evolution of the disintegrin family. The topology of the phylogenetic tree does not parallel that of the species tree. This incongruence is consistent with that expected for a multigene family undergoing a birth-and-death process in which the appearance and disappearance of loci are being driven by selection. Cysteine and buried residues appear to be under strong purifying selection due to their role in maintaining the active conformation of disintegrins. Divergence of disintegrins is strongly influenced by positive Darwinian selection causing accelerated rate of substitution in a substantial proportion of surface-exposed disintegrin residues. Global and lineage-specific sites evolving under diversifying selection were identified. Several sites are located within the integrin-binding loop and the C-terminal tail, two regions that form a conformational functional epitope. Arginine-glycine-aspartic acid (RGD) was inferred to represent the ancestral integrin-recognition motif, which emerged from the subgroup of PIII-SVMPs bearing the RDECD sequence. The most parsimonious nucleotide substitution model required for the emergence of all known disintegrin's integrin inhibitory motifs from an ancestral RGD sequence involves a minimum of three mutations. The adaptive advantage of the emergence of motifs targeting beta(1) integrins and the role of positively selected sites located within nonfunctional disintegrin regions appear to be difficult to rationalize in the context of a predator-prey arms race. Perhaps, this represents a consequence of the neofunctionalization potential of the disintegrin domain, a feature that may underlie its recruitment into the venom proteome followed by its successful transformation into a toxin.