Evolution of protein inhibitors of serine proteinases: positive Darwinian selection or compositional effects?

Evolutionary Aspects of the Structural Convergence and Functional Diversification of Kunitz-Domain Inhibitors

Kunitz-type domains are ubiquitously found in natural systems as serine protease inhibitors or animal toxins in venomous animals. Kunitz motif is a cysteine-rich peptide chain of ~ 60 amino acid residues with alpha and beta fold, stabilized by three conserved disulfide bridges. An extensive dataset of amino acid variations is found on sequence analysis of various Kunitz peptides. Kunitz peptides show diverse biological activities like inhibition of proteases of other classes and/or adopting a new function of blocking or modulating the ion channels. Based on the amino acid residues at the functional site of various Kunitz-type inhibitors, it is inferred that this 'flexibility within the structural rigidity' is responsible for multiple biological activities. Accelerated evolution of functional sites in response to the co-evolving molecular targets of the hosts of venomous animals or parasites, gene sharing, and gene duplication have been discussed as the most likely mechanisms responsible for the functional heterogeneity of Kunitz-domain inhibitors.

DOES EVOLUTIONARY PLASTICITY EVOLVE?

During the development of a multicellular organism from a zygote, a large number of epigenetic interactions take place on every level of suborganismal organization. This raises the possibility that the system of epigenetic interactions may compensate or "buffer" some of the changes that occur as mutations on its lowest levels, and thus stabilize the phenotype with respect to mutations. This hypothetical phenomenon will be called "epigenetic stability." Its potential importance stems from the fact that phenotypic variation with a genetic basis is an essential prerequisite for evolution. Thus, variation in epigenetic stability might profoundly affect attainable rates of evolution. While representing a systemic property of a developmental system, epigenetic stability might itself be genetically determined and thus be subject to evolutionary change. Whether or not this is the case should ideally be answered directly, that is, by experimentation. The time scale involved and our insufficient quantitative understanding of developmental pathways will probably preclude such an approach in the foreseeable future. Preliminary answers are sought here by using a biochemically motivated model of a small but central part of a developmental pathway. Modeled are sets of transcriptional regulators that mutually regulate each other's expression and thereby form stable gene expression patterns. Such gene-expression patterns, crucially involved in determining developmental pattern formation events, are most likely subject to strong stabilizing natural selection. After long periods of stabilizing selection, the fraction of mutations causing changes in gene-expression patterns is substantially reduced in the model. Epigenetic stability has increased. This phenomenon is found for widely varying regulatory scenarios among transcription factor genes. It is discussed that only epistatic (nonlinear) gene interactions can cause such change in epigenetic stability. Evidence from paleontology, molecular evolution, development, and genetics, consistent with the existence of variation in epigenetic stability, is discussed. The relation of epigenetic stability to developmental canalization is outlined. Experimental scenarios are suggested that may provide further evidence.

Differential transcriptome analysis supports Rhodnius montenegrensis and Rhodnius robustus (Hemiptera, Reduviidae, Triatominae) as distinct species

Chagas disease is one of the main parasitic diseases found in Latin America and it is estimated that between six and seven million people are infected worldwide. Its etiologic agent, the protozoan Trypanosoma cruzi, is transmitted by triatomines, some of which from the genus Rhodnius. Twenty species are currently recognized in this genus, including some closely related species with low levels of morphological differentiation, such as Rhodnius montenegrensis and Rhodnius robustus. In order to investigate genetic differences between these two species, we generated large-scale RNA-sequencing data (consisting of four RNA-seq libraries) from the heads and salivary glands of males of R. montenegrensis and R. robustus. Transcriptome assemblies produced for each species resulted in 64,952 contigs for R. montenegrensis and 70,894 contigs for R. robustus, with N50 of approximately 2,100 for both species. SNP calling based on the more complete R. robustus assembly revealed 3,055 fixed interspecific differences and 216 transcripts with high levels of divergence which contained only fixed differences between the two species. A gene ontology enrichment analysis revealed that these highly differentiated transcripts were enriched for eight GO terms related to AP-2 adaptor complex, as well as other interesting genes that could be involved in their differentiation. The results show that R. montenegrensis and R. robustus have a substantial quantity of fixed interspecific polymorphisms, which suggests a high degree of genetic divergence between the two species and likely corroborates the species status of R. montenegrensis.

Serine protease inhibitors of parasitic helminths

Serine protease inhibitors (serpins) are a superfamily of structurally conserved proteins that inhibit serine proteases and play key physiological roles in numerous biological systems such as blood coagulation, complement activation and inflammation. A number of serpins have now been identified in parasitic helminths with putative involvement in immune regulation and in parasite survival through interference with the host immune response. This review describes the serpins and smapins (small serine protease inhibitors) that have been identified in Ascaris spp., Brugia malayi, Ancylostoma caninum Onchocerca volvulus, Haemonchus contortus, Trichinella spiralis, Trichostrongylus vitrinus, Anisakis simplex, Trichuris suis, Schistosoma spp., Clonorchis sinensis, Paragonimus westermani and Echinococcus spp. and discusses their possible biological functions, including roles in host-parasite interplay and their evolutionary relationships.

Purification and biochemical characterization of ovine alpha-1-proteinase inhibitor: mechanistic adaptations and role of Phe350 and Met356

The glycoprotein alpha-1-proteinase inhibitor (alpha-1-PI) is a member of the serpin super family that causes rapid and irreversible inhibition of redundant serine protease activity. A homogenous preparation of ovine alpha-1-PI, a 60 kDa protein was obtained by serially subjecting ovine serum to 40-70% (NH(4))(2)SO(4) precipitation, Blue Sepharose, size-exclusion, and concanavalin-A chromatography. Extensive insights into the trypsin, chymotrypsin, and elastase interaction with ovine alpha-1-PI, point towards the involvement of Phe(350) besides the largely conserved Met(356) in serine protease recognition and consequent inhibition. The N-terminal of C-terminal peptides cleaved on interaction with elastase, trypsin, and chymotrypsin prove the presence of diffused sub-sites in the vicinity of Met(356) and the strategically positioned Pro anchored peptide stretch. Further, human alpha-1-PI is more thermolabile compared to ovine alpha-1-PI, higher thermolability is mainly attributed to poorer glycosylation. The enzymatic deglycosylation of human and ovine alpha-1-PI results in diminished thermostability of the inhibitors, with sharp decrease in thermal transition temperatures but retaining their inhibitory potency. Homology modeling of the deduced amino acid sequence of ovine alpha-1-PI using the human alpha-1-PI template has been used to explain the observed inhibitor-protease interactions.

Serpins: structure, function and molecular evolution

The superfamily of serine proteinase inhibitors (serpins) are involved in a number of fundamental biological processes such as blood coagulation, complement activation, fibrinolysis, angiogenesis, inflammation and tumor suppression and are expressed in a cell-specific manner. The average protein size of a serpin family member is 350-400 amino acids, but gene structure varies in terms of number and size of exons and introns. Previous studies of all known serpins identified 16 clades and 10 orphan sequences. Vertebrate serpins can be conveniently classified into six sub-groups. We provide additional data that updates the phylogenetic analysis in the context of structural and functional properties of the proteins. From these, we can conclude that the functional classification of serpins relies on their protein structure and not on sequence similarity.

Additivity-based prediction of equilibrium constants for some protein-protein associations

For many protein families, such as serine proteinases or serine proteinase inhibitors, the family assignment predicts reactivity only in general terms. Both detailed specificity and quantitative reactivity are lacking. We believe that, for many such protein families, algorithms can be devised by defining the subset of n functionally important sequence positions, making the 19n possible single mutants and measuring their reactivity. Given the assumption that the contributions of the n positions are additive, the reactivities of the 20(n) variants can be predicted. This is illustrated by an almost complete algorithm for the Kazal family of protein inhibitors of serine proteinases.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

Molecular evolution of interleukin-3

Chimpanzee, tamarin, and marmoset interleukin-3 (IL-3) genes were cloned, sequenced, and expressed. Western blot analysis demonstrated that functional genes were isolated. IL-3 sequences were compared with those of mouse, rat, rhesus monkey, gibbon, and man. Multiple alignment of the IL-3 coding regions showed that only a few regions had been conserved during mammalian evolution, which are likely associated with functional domains of the IL-3 protein. Substitution rates for the various lineages were calculated and the numbers of synonymous and nonsynonymous substitutions were estimated separately. Distance matrices of the IL-3 coding regions were used to construct phylogenetic trees which revealed large differences in IL-3 evolution rate as well as a more rapid substitution rate for rodents and a rate slowdown during hominoid evolution. Extremes were rhesus monkey IL-3, which accumulated few synonymous substitutions, and gibbon IL-3, which had almost exclusively synonymous substitutions. In rhesus monkey IL-3, nonsynonymous substitutions outnumbered synonymous substitutions, which could not be readily explained by a random process of substitutions. We assume that during evolution of IL-3, the majority of the amino acid replacements and the impaired interspecies functional cross-reactivity originate from selection mechanisms with the most likely selective force being the structure of the heterodimeric IL.3 cell-surface receptor. Insight into IL-3 architecture and structural analysis of the IL-3 receptor are needed to analyze the unusually fast evolution of IL-3 in more detail.

Nucleotide composition as a driving force in the evolution of retroviruses

All complete retrovirus sequences in the GenEMBL database were examined with the goal of assessing possible relationships between the nucleotide composition of retroviral genomes, the amino acid composition of retroviral proteins, and evolutionary strategies used by retroviruses. The results demonstrated that the genome of each viral lineage has a characteristic base composition and that the variations between groups are related to retroviral phylogeny. By analogy to microbial species, we suggest that the variations arise from group-specific patterns of directional mutations where the bias can be exerted on any of the four nucleotides. It is most likely that the mutational patterns are introduced during reverse transcription, and a direct participation of reverse transcriptase in the process is suspected. A straightforward strategy was used to analyze the compositional relationship between nucleotides and encoded amino acids. The procedure entailed calculations of amino acid frequencies from nucleotide content and the comparison of the calculated values to the observed amino acid frequencies in retroviruses. The results revealed an excellent correspondence between variation in genomic base composition and variation in amino acid composition of proteins with the compositional differences extending into all major coding regions of the viruses. Because of the magnitude and dispersion of these effects, and because of the nonconservative nature of many of the substitutions between groups with different genomic biases, we suggest that the variations in protein composition driven by biased nucleotide frequencies are an important factor in shaping the characteristic phenotypes of the different viral lineages. A clue to the nature of the evolutionary forces that are responsible for the generation of nucleotide biases was provided by the observation that viruses with radically different base frequencies most often inhabit the same cell type. This observation, along with analysis of amino acid and nucleotide replacement patterns between and within reverse transcriptase sequences from the various groups, permitted us to advance a model for the evolution of retroviruses. According to the model, speciation could initiate when daughter virions from a single progenitor vary in the direction of their mutational bias. These variations would exert a pleiotropic effect on the frequencies of nucleotides in all viral genes and consequently on the frequencies of amino acids in the encoded proteins. The variants with the most extreme compositional differences would have a selective advantage because their different precursor requirements would enable them to occupy different ecological niches within a single cell.(ABSTRACT TRUNCATED AT 400 WORDS)

Evolutionary origin of a Kunitz-type trypsin inhibitor domain inserted in the amyloid beta precursor protein of Alzheimer's disease

The Kunitz-type protease inhibitor is one of the serine protease inhibitors. It is found in blood, saliva, and all tissues in mammals. Recently, a Kunitz-type sequence was found in the protein sequence of the amyloid beta precursor protein (beta APP). It is known that beta APP accumulates in the neuritic plaques and cerebrovascular deposits of patients with Alzheimer's disease. Collagen type VI in chicken also has an insertion of a Kunitz-type sequence. To elucidate the evolutionary origin of these insertion sequences, we constructed a phylogenetic tree by use of all the available sequences of Kunitz-type inhibitors. The tree shows that the ancestral gene of the Kunitz-type inhibitor appeared about 500 million years ago. Thereafter, this gene duplicated itself many times, and some of the duplicates were inserted into other protein-coding genes. During this process, the Kunitz-type sequence in the present beta APP gene diverged from its ancestral gene about 270 million years ago and was inserted into the gene soon after duplication. Although the function of the insertion sequences is unknown, our molecular evolutionary analysis shows that these insertion sequences in beta APP have an evolutionarily close relationship with the inter-alpha-trypsin inhibitor or trypstatin, which inhibits the activity of tryptase, a novel membrane-bound serine protease in human T4+ lymphocytes.

Extreme divergence of a novel wheat thionin generated by a mutational burst specifically affecting the mature protein domain of the precursor

A new type of neutral thionin (type V), specifically expressed in developing wheat endosperm, has been found to be encoded by a set of single-copy genes located in the long arms of chromosomes 1A, 1B and 1D, within less than 10,000 base-pairs of those corresponding to the highly basic type-I thionins. Divergence between types I and V has occurred through a process of accelerated evolution that has affected the amino acid sequence of the mature thionin but not the precursor domains corresponding to the N-terminal signal peptide and the long C-terminal acidic peptide. This process involved a deletion and a non-synonymous nucleotide substitution rate equal to the synonymous rate in the thionin sequence.

Extreme variations in the ratios of non-synonymous to synonymous nucleotide substitution rates in signal peptide evolution

Nucleotide sequences encoding signal peptides from the precursors of alpha-amylase/trypsin inhibitors from cereals are homologous to those corresponding to the precursors of thaumatin II and of plastocyanins. Non-synonymous (KA) and synonymous (KS) rates of nucleotide substitutions have been calculated for all possible binary combinations. Extreme variation in KA/KS ratios has been observed; from the 0.167 average found within the plastocyanin family to an average of 1.90 calculated for the inhibitors/thaumatin II transition. A similar calculation has been carried out for the signal peptide sequences of thionins, which are unrelated to those of the alpha-amylase/trypsin inhibitor family, and an average KA/KS of 0.12 has been obtained. This variation can be largely explained in terms of an empirical index of stability related to amino acid composition and seems to be independent of functional constraints.

Mammalian alpha 1-antitrypsins: comparative biochemistry and genetics of the major plasma serpin

1. Human alpha 1-antitrypsin (alpha 1 AT) has been extensively characterized and reviewed. It is the archetypal member of the superfamily of serine proteinase inhibitors, the serpins. As human alpha 1-antitrypsin exhibits a relatively high concentration in plasma and is usually the highest concentration serpin, it can be referred to as the major plasma serpin. 2. alpha 1-Antitrypsin from species other than man has been characterized for two major reasons: (1) for use in a model animal system to assist with the study of the human alpha 1 AT deficiency disease; and (2) to find polymorphism for use in gene mapping and linkage studies or for parentage analysis. 3. The diverse range of reasons for studying alpha 1AT has yielded a vast array of literature that is often not well cross-referenced. 4. The characteristic features of alpha 1AT in all species examined to date will be presented with a view to examining which features are important structurally and functionally from an evolutionary perspective. 5. In mouse, horse, rabbit and guinea pig, multigene families which appear to have arisen from alpha 1AT have been found. The functional and evolutionary implications of these paralogous genes will also be discussed.

Evolution of sequences encoding the principal neutralization epitope of human immunodeficiency virus 1 is host dependent, rapid, and continuous

The principal neutralization epitope of human immunodeficiency virus 1 is localized in the third variable (V3) domain of the external envelope and has been shown to bind isolate-specific antibodies. Therefore, the extent of variation within the nucleic acid sequence encoding this epitope was studied in DNA directly obtained from peripheral blood mononuclear cells of six children and their plasma donor. This revealed that the quasi-species distribution of sequences obtained after cloning varied from recipient to recipient and that the distance from the donor sequences increased over time. V3 nucleotide evolution rates averaged 9.5 x 10(-3) per site per year for silent sites and 11.4 x 10(-3) per site per year for nonsilent sites (vs. 9.7 and 9.8 x 10(-3) per site per year for a control region 5' adjacent to the V3 region) and, although individual differences were observed, did not correlate with the serum antigen levels or disease progression. Sequences of both the epitope coding region itself (V3) and the control region upstream diverted more from the donor sequence among children not progressing to AIDS than among children progressing to AIDS. The evolution of V3 sequences is apparently host dependent, rapid, and independent of the level of antigen expression.

Amino acid sequences of ovomucoid third domain from 25 additional species of birds

Ovomucoids were isolated from 25 avian species other than the 101 studied in Laskowski et al. (1987, Biochemistry 26, 202-221). These were subjected to limited proteolysis with an appropriate enzyme, and connecting peptide extended ovomucoid third domains were isolated and sequenced to the end in a protein sequencer. Of the 25 new sequences, 13 duplicate ones were already known, and 12 are unique. Probably the most striking findings are a Pro14----Ser14 replacement in weka, an Ala14----Thr15 replacement in Bulwer's pheasant, the discovery of two additional amino acid residues Ile18 and Gly18 at the P1 reactive site position in Kalij pheasant and tawny frogmouth, respectively, and the first finding of a negative (Glu34) rather than positive (Lys34 or Arg34) amino acid residue at the NH2 terminus of the alpha helix in caracara ovomucoid third domain. These results complete the determination of all the sequences of ovomucoid third domains in the four species genus Gallus, in the five species genus Syrmaticus, and in the two species genera Aix and Pavo.

A maximum likelihood approach to the detection of selection from a phylogeny

A large amount of information is contained within the phylogenetic relationships between species. In addition to their branching patterns it is also possible to examine other aspects of the biology of the species. The influence that deleterious selection might have is determined here. The likelihood of different phylogenies in the presence of selection is explored to determine the properties of such a likelihood surface. The calculation of likelihoods for a phylogeny in the presence and absence of selection, permits the application of a likelihood ratio test to search for selection. It is shown that even a single selected site can have a strong effect on the likelihood. The method is illustrated with an example from Drosophila melanogaster and suggests that deleterious selection may be acting on transposable elements.