WRKY gene family evolution in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

The Arabidopsis thaliana WRKY proteins are characterized by a sequence of 60 amino acids including WRKY domain. It is well established that these proteins are involved in the regulation of various physiological programs unique to plants including pathogen defense, senescence and response to environmental stresses, which attracts attention of the scientific community as to how this family might have evolved. We tried to satisfy this curiosity and analyze reasons for duplications of these gene sequences leading to their diversified gene actions. The WRKY sequences available in Arabidopsis thaliana were used to evaluate selection pressure following duplication events. A phylogenetic tree was constructed and the WRKY family was divided into five sub-families. After that, tests were conducted to decide whether positive or purified selection played key role in these events. Our results suggest that purifying selection played major role during the evolution of this family. Some amino acid changes were also detected in specific branches of phylogeny suggesting that relaxed constraints might also have contributed to functional divergence among sub-families. Sites relaxed from purifying selection were identified and mapped onto the structural and functional regions of the WRKY1 protein. These analyses will enhance our understanding of the precise role played by natural selection to create functional diversity in WRKY family.     

The divergence and positive selection of the plant‐specific BURP‐containing protein family

BURP domain‐containing proteins belong to a plant‐specific protein family and have diverse roles in plant development and stress responses. However, our understanding about the genetic divergence patterns and evolutionary rates of these proteins remain inadequate. In this study, 15 plant genomes were explored to elucidate the genetic origins, divergence, and functions of these proteins. One hundred and twenty‐five BURP protein‐encoding genes were identified from four main plant lineages, including 13 higher plant species. The absence of BURP family genes in unicellular and multicellular algae suggests that this family (1) appeared when plants shifted from relatively stable aquatic environments to land, where conditions are more variable and stressful, and (2) is critical in the adaptation of plants to adverse environments. Promoter analysis revealed that several responsive elements to plant hormones and external environment stresses are concentrated in the promoter region of BURP protein‐encoding genes. This finding confirms that these genes influence plant stress responses. Several segmentally and tandem‐duplicated gene pairs were identified from eight plant species. Thus, in general, BURP domain‐containing genes have been subject to strong positive selection, even though these genes have conformed to different expansion models in different species. Our study also detected certain critical amino acid sites that may have contributed to functional divergence among groups or subgroups. Unexpectedly, all of the critical amino acid residues of functional divergence and positive selection were exclusively located in the C‐terminal region of the BURP domain. In conclusion, our results contribute novel insights into the genetic divergence patterns and evolutionary rates of BURP proteins.     

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.     

Evolution of the class C GPCR Venus flytrap modules involved positive selected functional divergence

Background  

Class C G protein-coupled receptors (GPCRs) represent a distinct group of the GPCR family, which structurally possess a characteristically distinct extracellular domain inclusive of the Venus flytrap module (VFTM). The VFTMs of the class C GPCRs is responsible for ligand recognition and binding, and share sequence similarity with bacterial periplasmic amino acid binding proteins (PBPs). An extensive phylogenetic investigation of the VFTMs was conducted by analyzing for functional divergence and testing for positive selection for five typical groups of the class C GPCRs. The altered selective constraints were determined to identify the sites that had undergone functional divergence via positive selection. In order to structurally demonstrate the pattern changes during the evolutionary process, three-dimensional (3D) structures of the GPCR VFTMs were modelled and reconstructed from ancestral VFTMs.     

Combining specificity determining and conserved residues improves functional site prediction

Background  

Predicting the location of functionally important sites from protein sequence and/or structure is a long-standing problem in computational biology. Most current approaches make use of sequence conservation, assuming that amino acid residues conserved within a protein family are most likely to be functionally important. Most often these approaches do not consider many residues that act to define specific sub-functions within a family, or they make no distinction between residues important for function and those more relevant for maintaining structure (e.g. in the hydrophobic core). Many protein families bind and/or act on a variety of ligands, meaning that conserved residues often only bind a common ligand sub-structure or perform general catalytic activities.     

On the unique function of selenocysteine — Insights from the evolution of selenoproteins

The process of natural selection leaves signatures in our genome that can be used to identify functionally important amino acid changes in proteins. In natural populations, amino acids that are better adapted to local conditions might increase in frequency, whereas moderately to severely deleterious protein mutations tend to be eliminated and do not contribute to protein differences between species. Amino acid mutations with no fitness consequences are, however, lost or fixed without regard to natural selection. Looking for evidence of natural selection is, therefore, an attractive strategy for characterizing the contribution of a residue to protein function. Because the majority of identified selenoproteins have now been found in Cys-form, the extent of exchangeability between Sec and Cys residues can be measured in proteins over long periods of time. The statistical analysis of the pattern of Sec/Cys exchanges in diversity (within species) and divergence (between species) data, provides robust inferences of the strength and mode of natural selection acting on these protein sites. Such inferences inform us not only of the long-term exchangeability between Sec and Cys residues, but also of the nature of the selective factors shaping Sec usage in proteins.     

Predicting functionally important residues from sequence conservation

MOTIVATION: All residues in a protein are not equally important. Some are essential for the proper structure and function of the protein, whereas others can be readily replaced. Conservation analysis is one of the most widely used methods for predicting these functionally important residues in protein sequences. RESULTS: We introduce an information-theoretic approach for estimating sequence conservation based on Jensen-Shannon divergence. We also develop a general heuristic that considers the estimated conservation of sequentially neighboring sites. In large-scale testing, we demonstrate that our combined approach outperforms previous conservation-based measures in identifying functionally important residues; in particular, it is significantly better than the commonly used Shannon entropy measure. We find that considering conservation at sequential neighbors improves the performance of all methods tested. Our analysis also reveals that many existing methods that attempt to incorporate the relationships between amino acids do not lead to better identification of functionally important sites. Finally, we find that while conservation is highly predictive in identifying catalytic sites and residues near bound ligands, it is much less effective in identifying residues in protein-protein interfaces. AVAILABILITY: Data sets and code for all conservation measures evaluated are available at http://compbio.cs.princeton.edu/conservation/     

Functional divergence after gene duplication and sequence-structure relationship: a case study of G-protein alpha subunits

In this article, we use animal G-protein alpha subunit family as an example to illustrate a comprehensive analytical pipeline for detecting different types of functional divergence of protein families, which is phylogeny-dependent, combined with ancestral sequence inference and available protein structure information. In particular, we focus on (i) Type-I functional divergence, or site-specific rate shift, as typically exemplified by amino acid residue highly conserved in a subset of homologous genes but highly variable in a different subset of homologous genes, and (ii) Type-II functional divergence, or the shift of cluster-specific amino acid property, as exemplified by a radical shift of amino acid property between duplicate genes, which is otherwise evolutionally conserved. We utilized the software DIVERGE2 to carry out these analyses. In the case of G-protein alpha subunit gene family, we have predicted amino acid residues that are related to either Type-I or Type-II functional divergence. The inferred ancestral sequences for these sites are helpful to explore the trends of functional divergence. Finally, these predicted residues are mapped to the protein structures to test whether these residues may have 3D structure or solvent accessibility preference.     

Adaptive evolution of crustin antimicrobial peptides in decapods

Antimicrobial peptides (AMPs) are present in a wide range of taxonomic groups and played a crucial role in host adaptation to a diverse array of ever-changing pathogens. Crustin, a cysteine-rich cationic AMP, is known to exhibit antimicrobial activity against Gram-positive and Gram-negative bacteria in decapods. Given their important role in host-immune defense, a large proportion of amino acid substitutions in crustin AMPs are expected to be fixed by natural selection. Utilizing the complete coding nucleotide sequence data of crustin, the present study revealed the pervasive role of positive Darwinian selection in the evolution and divergence of crustin AMPs in decapods. Approximately, 20–35?% of codons in two phylogenetically distinct groups of closely related crustins in Penaeid shrimps are shown to have evolved under positive selection. Interestingly, several of these positively selected sites are located at the carboxyl-terminal region, the region that directly interacts with the invading pathogens. Pathogen-mediated selection pressure could be the likely cause for such an accelerated rate of amino acid substitutions and could have contributed to the structural and functional diversification of crustin AMPs in several taxa.     

Diversity of the C-terminal portion of the biphenyl dioxygenase large subunit

The biphenyl dioxygenase (BPDO) catalyses a stereospecific dioxygenation of biphenyl and analogs of it. Aside from being involved in the destruction and detoxification of toxic pollutants in soil, in the context of the green chemistry concept, this enzyme is a promising biocatalyst to design new more selective and more environmentally friendly approaches to manufacture fine chemicals. At this time, most of our knowledge about the variability of key residues determining the substrate specificity and regiospecificity of the enzyme oxygenase component (BphAE) toward biphenyl analogs and about the effect of altering these residues on catalytic properties is based on investigations made with BphAEs from cultured organisms and engineered enzymes derived from them. The purpose of this work was to examine the diversity of the amino acid sequence patterns of the alpha subunit (BphA) C-terminal domain deduced from PCR products amplified from DNA extracted from cultured bacteria of various phylogenetic lines and from the soil microflora of PCB-contaminated soils. Of special interest were segments of the C-terminal portion called regions I, III and IV. Altogether, the phylogenetic tree obtained from aligning the deduced amino acid sequences of BphAs C-terminal domain from cultured bacteria belonging to various ecological niches and from uncultured soil bacteria reveals that most of the BphAs were linked to the three clusters of BphAs previously reported. However, few belong to new branches that diverge from the previously known branches showing a high diversity of BphAs in natural environment. Furthermore, data show a wide distribution of BphAs with family linkages that not only crosses bacterial taxonomic frontiers but also ecological niches. Nevertheless, in spite of this divergence, the sequence patterns of regions III and IV amino acids that are known to influence substrate specificity and regiospecificity are rather conserved among BphAs and the pattern was independent of the family cluster to which they belong. In most cases, regions III and IV amino acid patterns are closer to those of Pseudomonas pseudoalcaligenes KF707 BphA1 than to the most versatile Burkholderia xenovorans LB400 BphA. This might suggest that the PCB-degrading potency of soil bacteria is closer to the one observed for KF707 BphAE than from LB400 BphAE. However, the fact that among less than 20 PCR products amplified from soil DNA that we have sequenced, one of them was very homologous to that of LB400 BphA and in addition, residues 335 and 336 of LB400 were replaced by residues that previous enzyme engineering had shown to extend the range of PCB substrate used by the enzyme strongly suggest that PCB-degrading bacteria are evolving in soil to optimize their PCB-degrading capacity.     

Kinesin-related genes from diplomonad,sponge, amphioxus,and cyclostomes: divergence pattern of kinesin family and evolution of giardial membrane-bounded organella

To understand the question of whether divergence of eukaryotic genes by gene duplications and domain shufflings proceeded gradually or intermittently during evolution, we have cloned and sequenced Giardia lamblia cDNAs encoding kinesins and kinesin-related proteins and have obtained 13 kinesin-related cDNAs, some of which are likely homologs of vertebrate kinesins involved in vesicle transfer to ER, Golgi, and plasma membrane. A phylogenetic tree of the kinesin family revealed that most gene duplications that gave rise to different kinesin subfamilies with distinct functions have been completed before the earliest divergence of extant eukaryotes. This suggests that the complex endomembrane system has arisen very early in eukaryotic evolution, and the diminutive ER and Golgi apparatus recognized in the giardial cells, together with the absence of mitochondria, might be characters acquired secondarily during the evolution of parasitism. To understand the divergence pattern of the kinesin family in the lineage leading to vertebrates, seven more Unc104-related cDNAs have been cloned from sponge, amphioxus, hagfish, and lamprey. The divergence pattern of the animal Unc104/KIF1 subfamily is characterized by two active periods in gene duplication interrupted by a considerably long period of silence, instead of proceeding gradually: animals underwent extensive gene duplications before the parazoan-eumetazoan split. In the early evolution of vertebrates around the cyclostome-gnathostome split, further gene duplications occurred, by which a variety of genes with similar structures over the entire regions were generated. This pattern of divergence is similar to those of animal genes involved in cell-cell communication and developmental control.     

A new branch in the family: structure of aspartate-beta-semialdehyde dehydrogenase from Methanococcus jannaschii

The structure of aspartate-beta-semialdehyde dehydrogenase (ASADH) from Methanococcus jannaschii has been determined to 2.3 angstroms resolution using multiwavelength anomalous diffraction (MAD) phasing of a selenomethionine-substituted derivative to define a new branch in the family of ASADHs. This new structure has a similar overall fold and domain organization despite less than 10% conserved sequence identity with the bacterial enzymes. However, the entire repertoire of functionally important active site amino acid residues is conserved, suggesting an identical catalytic mechanism but with lower catalytic efficiency. A new coenzyme-binding conformation and dual NAD/NADP coenzyme specificity further distinguish this archaeal branch from the bacterial ASADHs. Several structural differences are proposed to account for the dramatically enhanced thermostability of this archaeal enzyme. Finally, the intersubunit communication channel connecting the active sites in the bacterial enzyme dimer has been disrupted in the archaeal ASADHs by amino acid changes that likely prevent the alternating sites reactivity previously proposed for the bacterial ASADHs.     

The catalytic and lectin domains of UDP-GalNAc:polypeptide alpha-N-Acetylgalactosaminyltransferase function in concert to direct glycosylation site selection

UDP-GalNAc:polypeptide alpha-N-Acetylgalactosaminyltransferases (ppGalNAcTs), a family (EC 2.4.1.41) of enzymes that initiate mucin-type O-glycosylation, are structurally composed of a catalytic domain and a lectin domain. Previous studies have suggested that the lectin domain modulates the glycosylation of glycopeptide substrates and may underlie the strict glycopeptide specificity of some isoforms (ppGalNAcT-7 and -10). Using a set of synthetic peptides and glycopeptides based upon the sequence of the mucin, MUC5AC, we have examined the activity and glycosylation site preference of lectin domain deletion and exchange constructs of the peptide/glycopeptide transferase ppGalNAcT-2 (hT2) and the glycopeptide transferase ppGalNAcT-10 (hT10). We demonstrate that the lectin domain of hT2 directs glycosylation site selection for glycopeptide substrates. Pre-steady-state kinetic measurements show that this effect is attributable to two mechanisms, either lectin domain-aided substrate binding or lectin domain-aided product release following glycosylation. We find that glycosylation of peptide substrates by hT10 requires binding of existing GalNAcs on the substrate to either its catalytic or lectin domain, thereby resulting in its apparent strict glycopeptide specificity. These results highlight the existence of two modes of site selection used by these ppGalNAcTs: local sequence recognition by the catalytic domain and the concerted recognition of distal sites of prior glycosylation together with local sequence binding mediated, respectively, by the lectin and catalytic domains. The latter mode may facilitate the glycosylation of serine or threonine residues, which occur in sequence contexts that would not be efficiently glycosylated by the catalytic domain alone. Local sequence recognition by the catalytic domain differs between hT2 and hT10 in that hT10 requires a pre-existing GalNAc residue while hT2 does not.     

Evolutionary Analysis for Functional Divergence of the Toll-Like Receptor Gene Family and Altered Functional Constraints

The Toll-like receptor (TLR) gene family consists of type 1 transmembrane receptors, which play essential roles in both innate immunity and adaptive immune response by ligand recognition and signal transduction. Using all available vertebrate TLR protein sequences, we inferred the phylogenetic tree and then characterized critical amino acid residues for functional divergence by detecting altered functional constraints after gene duplications. We found that the extracellular domain of TLR genes showed higher functional divergence than that of the cytoplasmic domain, particularly in the region between leucine-rich repeat (LRR) 10 and LRR 15 of TLR 4. Our finding supports the concept that sequence evolution in the extracellular domain may be responsible for the broad diversity of TLR ligand-binding affinity, providing a testable hypothesis for potential targets that could be verified by further experimentation.     

Molecular Evolution of Miraculin-Like Proteins in Soybean Kunitz Super-Family

Miraculin-like proteins (MLPs) belong to soybean Kunitz super-family and have been characterized from many plant families like Rutaceae, Solanaceae, Rubiaceae, etc. Many of them possess trypsin inhibitory activity and are involved in plant defense. MLPs exhibit significant sequence identity (~30-95%) to native miraculin protein, also belonging to Kunitz super-family compared with a typical Kunitz family member (~30%). The sequence and structure-function comparison of MLPs with that of a classical Kunitz inhibitor have demonstrated that MLPs have evolved to form a distinct group within Kunitz super-family. Sequence analysis of new genes along with available MLP sequences in the literature revealed three major groups for these proteins. A significant feature of Rutaceae MLP type 2 sequences is the presence of phosphorylation motif. Subtle changes are seen in putative reactive loop residues among different MLPs suggesting altered specificities to specific proteases. In phylogenetic analysis, Rutaceae MLP type 1 and type 2 proteins clustered together on separate branches, whereas native miraculin along with other MLPs formed distinct clusters. Site-specific positive Darwinian selection was observed at many sites in both the groups of Rutaceae MLP sequences with most of the residues undergoing positive selection located in loop regions. The results demonstrate the sequence and thereby the structure-function divergence of MLPs as a distinct group within soybean Kunitz super-family due to biotic and abiotic stresses of local environment.     

Characterization of a gene encoding a Pichia pastoris protein disulfide isomerase

Protein disulphide isomerases belong to the thioredoxin superfamily of protein-thiol oxidoreductases that have two double-cysteine redox-active sites and take part in protein folding in the endoplasmic reticulum (ER). We report here the cloning of a Pichia pastoris genomic DNA fragment (2919 bp) that encodes the full length of a protein disulphide isomerase (PpPDI). The deduced amino acid sequence of PDI consists of 517 residues and carries the two characteristic PDI-type redox-active domains -CGHC-, separated by 338 residues, and two potential N-glycosylation sites. The N-terminal end forms a putative signal sequence, and an acidic C-terminal region represents a possible calcium-binding domain. Together with the -HDEL ER retrieval sequence at the C-terminus, these features indicate that the gene encodes a redox-active ER-resident protein disulphide isomerase. The nucleotide sequence, which also contains two other open reading frames, has been submitted to the EMBL Nucleotide Sequence Database, Accession No. AJ302014.     

Detecting coevolution in and among protein domains

Correlated changes of nucleic or amino acids have provided strong information about the structures and interactions of molecules. Despite the rich literature in coevolutionary sequence analysis, previous methods often have to trade off between generality, simplicity, phylogenetic information, and specific knowledge about interactions. Furthermore, despite the evidence of coevolution in selected protein families, a comprehensive screening of coevolution among all protein domains is still lacking. We propose an augmented continuous-time Markov process model for sequence coevolution. The model can handle different types of interactions, incorporate phylogenetic information and sequence substitution, has only one extra free parameter, and requires no knowledge about interaction rules. We employ this model to large-scale screenings on the entire protein domain database (Pfam). Strikingly, with 0.1 trillion tests executed, the majority of the inferred coevolving protein domains are functionally related, and the coevolving amino acid residues are spatially coupled. Moreover, many of the coevolving positions are located at functionally important sites of proteins/protein complexes, such as the subunit linkers of superoxide dismutase, the tRNA binding sites of ribosomes, the DNA binding region of RNA polymerase, and the active and ligand binding sites of various enzymes. The results suggest sequence coevolution manifests structural and functional constraints of proteins. The intricate relations between sequence coevolution and various selective constraints are worth pursuing at a deeper level.     

Phylogenetic analysis and evolution of aromatic amino acid hydroxylase

A study was performed to investigate the phylogenetic relationship among AAAH members and to statistically evaluate sequence conservation and functional divergence. In total, 161 genes were identified from 103 species. Phylogenetic analysis showed that well-conserved subfamilies exist. Exon-intron structure analysis showed that the gene structures of AAAH were highly conserved across some different lineage species, while some species-specific introns were also found. The dynamic distribution of ACT domain suggested one gene fusion event has occurred in eukaryota. Significant functional divergence was found between some subgroups. Analysis of the site-specific profiles revealed critical amino acid residues for functional divergence. This study highlights the molecular evolution of this family and may provide a starting point for further experimental verifications.     

A unique spacer domain of synaptotagmin IV is essential for Golgi localization

Synaptotagmin (Syt) family members consist of six separate domains: a short amino terminus, a single transmembrane domain, a spacer domain, a C2A domain, a C2B domain and a short carboxyl (C) terminus. Despite sharing the same domain structures, several synaptotagmin isoforms show distinct subcellular localization. Syt IV is mainly localized at the Golgi, while Syt I, a possible Ca(2+)-sensor for secretory vesicles, is localized at dense-core vesicles and synaptic-like microvesicles in PC12 cells. In this study, we sought to identify the region responsible for the Golgi localization of Syt IV by immunocytochemical and biochemical analyses as a means of defining the distinct subcellular localization of the synaptotagmin family. We found that the unique C-terminus of the spacer domain (amino acid residues 73-144) between the transmembrane domain and the C2A domain is essential for the Golgi localization of Syt IV. In addition, the short C-terminus is probably involved in proper folding of the protein, especially the C2B domain. Without the C-terminus, Syt IVdeltaC proteins are not targeted to the Golgi and seem to colocalize with an endoplasmic reticulum (ER) marker (i.e. induce crystalloid ER-like structures). On the basis of these results, we propose that the divergent spacer domain among synaptotagmin isoforms may contain certain signals that determine the final destination of each isoform.