Evolutionary Conservation of Protein Backbone Flexibility

Internal protein dynamics is essential for biological function. During evolution, protein divergence is functionally constrained: properties more relevant for function vary more slowly than less important properties. Thus, if protein dynamics is relevant for function, it should be evolutionary conserved. In contrast with the well-studied evolution of protein structure, the evolutionary divergence of protein dynamics has not been addressed systematically before, apart from a few case studies. X-Ray diffraction analysis gives information not only on protein structure but also on B-factors, which characterize the flexibility that results from protein dynamics. Here we study the evolutionary divergence of protein backbone dynamics by comparing the C_α flexibility (B-factor) profiles for a large dataset of homologous proteins classified into families and superfamilies. We show that C_α flexibility profiles diverge slowly, so that they are conserved at family and superfamily levels, even for pairs of proteins with nonsignificant sequence similarity. We also analyze and discuss the correlations among the divergences of flexibility, sequence, and structure. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. David Pollock]     

On the evolutionary conservation of hydrogen bonds made by buried polar amino acids: the hidden joists,braces and trusses of protein architecture

Background  

The hydrogen bond patterns between mainchain atoms in protein structures not only give rise to regular secondary structures but also satisfy mainchain hydrogen bond potential. However, not all mainchain atoms can be satisfied through hydrogen bond interactions that arise in regular secondary structures; in some locations sidechain-to-mainchain hydrogen bonds are required to provide polar group satisfaction. Buried polar residues that are hydrogen-bonded to mainchain amide atoms tend to be highly conserved within protein families, confirming that mainchain architecture is a critical restraint on the evolution of proteins. We have investigated the stabilizing roles of buried polar sidechains on the backbones of protein structures by performing an analysis of solvent inaccessible residues that are entirely conserved within protein families and superfamilies and hydrogen bonded to an equivalent mainchain atom in each family member.     

Molecular Evolution of the Plant R Regulatory Gene Family

Anthocyanin pigmentation patterns in different plant species are controlled in part by members of the myc-like R regulatory gene family. We have examined the molecular evolution of this gene family in seven plant species. Three regions of the R protein show sequence conservation between monocot and dicot R genes. These regions encode the basic helix-loop-helix domain, as well as conserved N-terminal and C-terminal domains; mean replacement rates for these conserved regions are 1.02 X 10(-9) nonsynonymous nucleotide substitutions per site per year. More than one-half of the protein, however, is diverging rapidly, with nonsynonymous substitution rates of 4.08 X 10(-9) substitutions per site per year. Detailed analysis of R homologs within the grasses (Poaceae) confirm that these variable regions are indeed evolving faster than the flanking conserved domains. Both nucleotide substitutions and small insertion/deletions contribute to the diversification of the variable regions within these regulatory genes. These results demonstrate that large tracts of sequence in these regulatory loci are evolving at a fairly rapid rate.     

Probing into the role of conserved N-glycosylation sites in the Tyrosinase glycoprotein family

N-linked glycosylation has a profound effect on the proper folding, oligomerization and stability of glycoproteins. These glycans impart many properties to proteins that may be important for their proper functioning, besides having a tendency to exert a chaperone-like effect on them. Certain glycosylation sites in a protein however, are more important than other sites for their function and stability. It has been observed that some N-glycosylation sites are conserved over families of glycoproteins over evolution, one such being the tyrosinase related protein family. The role of these conserved N-glycosylation sites in their trafficking, sorting, stability and activity has been examined here. By scrutinizing the different glycosylation sites on this family of glycoproteins it was inferred that different sites in the same family of polypeptides can perform distinct functions and conserved sites across the paralogues may perform diverse functions.     

Exon-intron structure of outlier tick lipocalins indicate a monophyletic origin within the larger lipocalin family

All tick proteins assigned to the lipocalin family lack the structural conserved regions (SCRs) that are characteristic of the kernel lipocalins and can thus be classified as outliers. These tick proteins have been assigned to the tick lipocalin family based on database searches that indicated homology between tick sequences and the fact that the histamine binding protein (HBP2) from the hard tick Rhipicephalus appendiculatus (Ixodidae) shows structural similarity to the lipocalin fold. Sequence identity between kernel and outlier lipocalins falls below 20% and the question raised is whether the outlier and kernel lipocalins are truly homologous. More specifically in the case of the tick lipocalins, whether their structural fold is derived from the lipocalin fold or whether convergent evolution resulted in the generation of the basic lipocalin-like fold which consists of an eight stranded continuous anti-parallel beta-barrel terminated by a C-terminal alpha-helix that lies parallel to the barrel. The current study determined the gene structure for HBP2 and TSGP1, TSGP2 and TSGP4, lipocalins identified from the soft tick Ornithodoros savignyi (Argasidae). All tick lipocalins have four introns (A-D) with conserved positions and phases within the tick lipocalin sequence alignment. The positions and phase information are also conserved with regard to the rest of the lipocalin family. Phylogenetic analysis using this information shows conclusively that tick lipocalins are evolutionary related to the rest of the lipocalin family. Tick lipocalins are grouped within a monophyletic clade that indicates a monophyletic origin within the tick lineage and also group with the other arthropod lipocalins in a larger clade. Phylogenetic analysis of sequence alignments based on conserved secondary structure of the lipocalin fold support the conclusions from the gene structure trees. These results indicate that exon-intron arrangement can be useful for the inclusion of outlier lipocalins within the larger lipocalin family.     

Evolution of structurally disordered proteins promotes neostructuralization

Protein structure is generally more conserved than sequence, but for regions that can adopt different structures in different environments, does this hold true? Understanding how structurally disordered regions evolve altered secondary structure element propensities as well as conformational flexibility among paralogs are fundamental questions for our understanding of protein structural evolution. We have investigated the evolutionary dynamics of structural disorder in protein families containing both orthologs and paralogs using phylogenetic tree reconstruction, protein structure disorder prediction, and secondary structure prediction in order to shed light upon these questions. Our results indicate that the extent and location of structurally disordered regions are not universally conserved. As structurally disordered regions often have high conformational flexibility, this is likely to have an effect on how protein structure evolves as spatially altered conformational flexibility can also change the secondary structure propensities for homologous regions in a protein family.     

What the evolution of the amyloid protein precursor supergene family tells us about its function

The Alzheimer's disease amyloid protein precursor (APP) gene is part of a multi-gene super-family from which sixteen homologous amyloid precursor-like proteins (APLP) and APP species homologues have been isolated and characterised. Comparison of exon structure (including the uncharacterised APL-1 gene), construction of phylogenetic trees, and analysis of the protein sequence alignment of known homologues of the APP super-family were performed to reconstruct the evolution of the family and to assess the functional significance of conserved protein sequences between homologues. This analysis supports an adhesion function for all members of the APP super family, with specificity determined by those sequences which are not conserved between APLP lineages, and provides evidence for an increasingly complex APP superfamily during evolution. The analysis also suggests that Drosophila APPL and Caenorhabditis elegans APL-1 may be a fourth APLP lineage indicating that these proteins, while not functional homologues of human APP, are similarly likely to regulate cell adhesion. Furthermore, the betaA4 sequence is highly conserved only in APP orthologues, strongly suggesting this sequence is of significant functional importance in this lineage.     

Domain evolution in the GH13 pullulanase subfamily with focus on the carbohydrate-binding module family 48

Glycoside hydrolase (GH) family 13 comprises about 30 different specificities. Four of them have been proposed to form the GH13 pullulanase subfamily: pullulanase, isoamylase, maltooligosyl trehalohydrolase and branching enzyme forming the seven CAZy GH13 subfamilies: GH13 8-GH13 14. Recently, a new family of carbohydrate-binding modules (CBMs), the family CBM48 has been established containing the putative starch-binding domains from the pullulanase subfamily, the β-subunit of AMP-activated protein kinase and some other GH13 enzymes with pullulanase and/or α-amylase-pullulanase specificity. Since all of these enzymes are multidomain proteins and the structure for at least one representative of each enzyme specificity has already been determined, the main goal of the present study was to elucidate domain evolution within this GH13 pullulanase subfamily (84 real enzymes) focusing on the CBM48 module. With regard to CBM48 positioning in the amino acid sequence, the N-terminal end of a protein appears to be a predominant position. This is especially true for isoamylases and maltooligosyl trehalohydrolases. Secondary structure-based alignment of CBM modules from CBM48, CBM20 and CBM21 revealed that several residues known as consensus for CBM20 and CBM21 could also be identified in CBM48, but only branching enzymes possess the aromatic residues that correspond with the two tryptophans forming the evolutionary conserved starch-binding site 1 in CBM20. The evolutionary trees constructed for the individual domains, complete alignment, and the conserved sequence regions of the α-amylase family were found to be comparable to each other (except for the C-domain tree) with two basic parts: (i) branching enzymes and maltooligosyl trehalohydrolases; and (ii) pullulanases and isoamylases. Taxonomy was respected only within clusters with pure specificity, i.e. the evolution of CBM48 reflects the evolution of specificities rather than evolution of species. This is a feature different from the one observed for the starch-binding domain of the family CBM20 where the starch-binding domain evolution reflects the evolution of species.     

Relation between sequence similarity and structural similarity in proteins. Role of important properties of amino acids

In a previous paper we obtained ten (orthogonal) factors, linear combinations of which can express the properties of the 20 naturally occurring amino acids. In this paper, we assume that the most important properties (linear combinations of these ten factors) that determine the three-dimensional structure of a protein are conserved properties, i.e., are those that have been conserved during evolution. Two definitions of a conserved property are presented: (1) a conserved property for an average protein is defined as that linear combination of the ten factors that optimally expresses the similarity of one amino acid to another (hence, little change during evolution), as given by the relatedness odds matrix of Dayhoff et al.; (2) a conserved property for each position in the amino acid sequence (locus) of a specific family of homologous proteins (the cytochromec family or the globin family) is defined as that linear combination of the ten factors that is common among a set of amino acids at a given locus when the sequences are properly aligned. When the specificity at each locus is averaged over all loci, the same features are observed for three expressions of these two definitions, namely the conserved property for an average protein, the average conserved property for the cytochromec family, and the average conserved property for the globin family; we find that bulk and hydrophobicity (information about packing and long-range interactions) are more important than other properties, such as the preference for adopting a specific backbone structure (information about short-range interactions). We also demonstrate that the sequence profile of a conserved property, defined for each locus of a protein family (definition 2), corresponds uniquely to the three-dimensional structure, while the conserved property for an average protein (definition 1) is not useful for the prediction of protein structure. The amino acid sequences of numerous proteins are searched to find those that are similar, in terms of the conserved properties (definition 2), to sequences of the same size from one of the homologous families (cytochromec and globin, respectively) for whose loci the conserved properties were defined. Many similar sequences are found, the number of similarities decreasing with increasing size of the segment. However, the segments must be rather long (15 residues) before the comparisons become meaningful. As an example, one sufficiently large sequence (20 residues) from a protein of known structure (apo-liver alcohol dehydrogenase that is not a member of either family) is found to be similar in the conserved properties to a particular sequence of a member of the family of human hemoglobin chains, and the two sequences have similar structures. This means that, since conserved properties are expected to be structure determinants, we can use the conserved properties to predict an initial protein structure for subsequent energy minimization for a protein for which the conserved properties are similar to those of a family of proteins with a sufficiently large number of homologous amino acid sequences; such a large number of homologous sequences is required to define a conserved property for each locus of the homologous protein family.     

Exploring some of the physico-chemical properties of the LAMMER protein kinase DOA of Drosophila

Members of the highly conserved LAMMER family of protein kinases have been described in all eukaryotes. LAMMER kinases possess markedly similar peptide motifs in their kinase catalytic subdomains that are responsible for phosphotransfer and substrate interaction, suggesting that family members serve similar functions in widely diverged species. This hypothesis is supported by their phosphorylation of SR and SR-related proteins in diverged species. Here we describe a 3-dimensional homology model of the catalytic domain of DOA, a representative LAMMER kinase, encoded by the Drosophila locus Darkener of apricot (Doa). Homology modeling of DOA based on a Sky1p template revealed a highly conserved structural framework within conserved core regions. These adopt typical kinase folding like that of other protein kinases. However, in contrast to Sky1p, some structural features, such as those in helix ?C suggest that the DOA kinase is not a constitutively active enzyme but requires activation. This may occur by phosphorylation within an activation loop that forms a broad turn and in which interactions between the side chains occur across the loop. The fold of the activation loop is stabilized through interactions with residues in the C-terminal tail, which is not part of the conserved kinase core and is variable among protein kinases. Immediately following the activation loop in the segment between the ?9 sheet and helix ?F is a P+1 loop. The electrostatic surface potential of the DOA substrate binding groove is largely negative, as it is in other known SR protein kinases, suggesting that DOA substrates must be basic. All differences between D. melanogaster and other Drosophila species are single amino acid changes situated in regions outside of any ?-helices or ?-sheets, and after modeling these had absolutely no visible effect on protein structure. The absence of evolved amino acid changes among 12 Drosophila species that would cause at least predictable changes in DOA structure indicate that evolution has already selected evolved mutations for having minimal effect on kinase structure.     

Comparison of scales of amino acid side chain properties by conservation during evolution of four proteins

As the amino acid sequence of a given protein changes along the phylogenetic tree, enough of the overall folding pattern must be conserved to ensure that the protein still fulfils its biological function. Eighteen published scales which tabulate various side chain properties are compared here by computing the variance of each scale when applied to each of several protein families. The conservation of each scale of side chain properties is examined for the 20,627 residues in 60 mammalian myoglobins, 31 mammalian ribonucleases, insulin A and B chains (29 sequences each), 29 vertebrate and 28 plant cytochrome c's. Those scales which are the most highly conserved through the evolution of each protein family may well be the best predictors of protein folding patterns. The mean-area-buried scale and the optimized matching hydrophobicities scale are more conserved than other scales. An additional result is the relatively poor conservation across evolution of the Chou-Fasman secondary structure predictors.     

The characterisation of six ADAMTS proteases in the basal chordate Ciona intestinalis provides new insights into the vertebrate ADAMTS family

ADAMTS, constituting a recently discovered family of secreted zinc-dependent metalloproteases, have been shown to have critical physiological roles through identification of a number of natural animal and human gene mutations. The identification of six ADAMTS genes in the basal chordate Ciona intestinalis provides new insight into how, when and in what order the vertebrate orthologues have evolved. The phylogenetic assignments, based on sequences conserved across all genes, are supported by conserved domain structures within defined sub-families. The phylogeny and the frequent localisation of ADAMTS genes in paralogous regions of the genome are consistent with the vertebrate lineages having arisen by large scale or genome duplication. The high level of conservation in the protease active site of vertebrate orthologues within some sub-families suggests subfunctionalisation, whereas the greater divergence in others would favour the evolution of novel substrate specificities and these observations are borne-out where substrate-specificity is known. The expansion and sub-specialization of the ADAMTS family is a component of the increased complexity of extracellular matrix that is associated with the evolution of vertebrates.     

Exploring the factors determining the dynamics of different protein folds

Normal mode analyses of homologous proteins at the family and superfamily level show that slow dynamics are similar and are preserved through evolution. This study investigates how the slow dynamics of proteins is affected by variation in the protein architecture and fold. For this purpose, we have used computer-generated protein models based on idealized protein structures with varying folds. These are shown to be protein-like in their behavior, and they are used to investigate the influence of architecture and fold on the slow dynamics. We compared the dynamics of models having different folds but similar architecture and found the architecture to be the dominant factor for the slow dynamics.     

The PUF family of RNA-binding proteins: does evolutionarily conserved structure equal conserved function?

Drosophila Pumilio (Pum) protein is a founder member of a novel family of RNA-binding proteins, known as the PUF family. The PUF proteins constitute an evolutionarily highly conserved family of proteins present from yeast to humans and plants, and are characterized by a highly conserved C-terminal RNA-binding domain, composed of eight tandem repeats. The conserved biochemical features and genetic function of PUF family members have emerged from studies of model organisms. PUF proteins bind to related sequence motifs in the 3' untranslated region (3'UTR) of specific target mRNAs and repress their translation. Frequently, PUF proteins function asymmetrically to create protein gradients, thus causing asymmetric cell division and regulating cell fate specification. Thus, it was recently proposed that the primordial role of PUF proteins is to sustain mitotic proliferation of stem cells. Here we review the evolution, conserved genetic and biochemical properties of PUF family of proteins, and discuss protein interactions, upstream regulators and downstream targets of PUF proteins. We also suggest that a conserved mechanism of PUF function extends to the newly described mammalian members of the PUF family (human PUM1 and PUM2, and mouse Pum1 and Pum2), that show extensive homology to Drosophila Pum, and could have an important role in cell development, fate specification and differentiation.     

Structural Determinants of Oligomerization of the Aquaporin-4 Channel

The aquaporin (AQP) family of integral membrane protein channels mediate cellular water and solute flow. Although qualitative and quantitative differences in channel permeability, selectivity, subcellular localization, and trafficking responses have been observed for different members of the AQP family, the signature homotetrameric quaternary structure is conserved. Using a variety of biophysical techniques, we show that mutations to an intracellular loop (loop D) of human AQP4 reduce oligomerization. Non-tetrameric AQP4 mutants are unable to relocalize to the plasma membrane in response to changes in extracellular tonicity, despite equivalent constitutive surface expression levels and water permeability to wild-type AQP4. A network of AQP4 loop D hydrogen bonding interactions, identified using molecular dynamics simulations and based on a comparative mutagenic analysis of AQPs 1, 3, and 4, suggest that loop D interactions may provide a general structural framework for tetrameric assembly within the AQP family.     

Evolutionary analysis of the mammalian M1 aminopeptidases reveals conserved exon structure and gene death

The members of the M1 aminopeptidase family share conserved domains, yet show functional divergence within the family as a whole. In order to better understand this family, this study analyzed the mammalian members in depth at exon, gene, and protein levels. The twelve human members, eleven rat members, and eleven mouse members were first analyzed in multiple alignments to visualize both reported and unreported conserved domains. Phylogenetic trees were then generated for humans, rats, mice, and all mammals to determine how closely related the homologs were and to gain insight to the divergence in the family members. This produced three groups with similarity within the family. Next, a synteny study was completed to determine the present locations of the genes and changes that had occurred. It became apparent that gene death likely resulted in the lack of one member in mouse and rat. Finally, an in-depth analysis of the exon structure revealed that nine members of the human family and eight in mouse, are highly conserved within the exon structure. Taken together, these results indicate that the M1 aminopeptidase family is a divergent family with three subgroups and that genetic evidence mirrors categorization of the family by enzymatic function.