Evolutionary relationship between the bacterial HPr kinase and the ubiquitous PEP-carboxykinase: expanding the P-loop nucleotidyl transferase superfamily

Similarities between protein three-dimensional structures can reveal evolutionary and functional relationships not apparent from sequence comparison alone. Here we report such a similarity between the metabolic enzymes histidine phosphocarrier protein kinase (HPrK) and phosphoenolpyruvate carboxykinase (PCK), suggesting that they are evolutionarily related. Current structure classifications place PCK and other P-loop containing nucleotidyl-transferases into different folds. Our comparison of both HPrK and PCK to other P-loop containing proteins reveals that all share a common structural motif consisting of an alphabeta segment containing the P-loop flanked by an additional beta-strand that is adjacent in space, but far apart along the sequence. Analysis also shows that HPrK/PCK differ from other P-loop containing structures no more than they differ from each other. We thus suggest that HPrK and PCK should be classified with other P-loop containing proteins, and that all probably share a common ancestor that probably contained a simple P-loop motif with different protein segments being added or lost over the course of evolution. We used the structure-based sequence alignment containing residues specific to HPrK/PCK to identify additional members of this P-loop containing family.     

Selective deletion of the NH2-terminal variable region of cardiac troponin T in ischemia reperfusion by myofibril-associated mu-calpain cleavage

The structure of the NH2-terminal region of troponin T (TnT) is hypervariable among the muscle type-specific isoforms and is also regulated by alternative RNA splicing. This region does not contain binding sites for other thin filament proteins, but alteration of its structure affects the Ca2+ regulation of muscle contraction. Here we report a truncated cardiac TnT produced during myocardial ischemia reperfusion. Amino acid sequencing and protein fragment reconstruction determined that it is generated by a posttranslational modification selectively removing the NH2-terminal variable region and preserving the conserved core structure of TnT. Triton X-100 extraction of cardiac muscle fibers promoted production of the NH2-terminal truncated cardiac TnT (cTnT-ND), indicating a myofibril-associated proteolytic activity. Mu-calpain is a myofibril-associated protease and is known to degrade TnT. Supporting a role of mu-calpain in producing cTnT-ND in myocardial ischemia reperfusion, calpain inhibitors decreased the level of cTnT-ND in Triton-extracted myofibrils. Mu-calpain treatment of the cardiac myofibril and troponin complex specifically reproduced cTnT-ND. In contrast, mu-calpain treatment of isolated cardiac TnT resulted in nonspecific degradation, suggesting that this structural modification is relevant to physiological structures of the myofilament. Triton X-100 treatment of transgenic mouse cardiac myofibrils overexpressing fast skeletal muscle TnT produced similar NH2-terminal truncations of the endogenous and exogenous TnT, despite different amino acid sequences at the cleavage site. With the functional consequences of removing the NH2-terminal variable region of TnT, the mu-calpain-mediated proteolytic modification of TnT may act as an acute mechanism to adjust muscle contractility under stress conditions.     

Protein evolution: structure-function relationships of the oncogene beta-catenin in the evolution of multicellular animals

Beta-catenin functions as a cytoskeletal linker protein in cadherin-mediated adhesion and as a signal mediator in wnt-signal transduction pathways. We use a novel integrative approach, combining evolutionary, genomic, and three-dimensional structural data to analyze and trace the structural and functional evolution of beta-catenin genes. This approach also enabled us to examine the effects of gene duplication on the structure and function of beta-catenin genes in Drosophila, C. elegans, and vertebrates. By sampling a large number of different taxa, we identified both ancestral and derived motifs and residues within the different regions of the beta-catenin proteins. Projecting amino acid substitutions onto the three- dimensional structure established for mouse beta-catenin, we identified specific domains that exhibit loss and gain of selective constraints during beta catenin evolution. Structural changes, changes in the amino acid substitution rate, and the appearance of novel functional domains in beta-catenin can be mapped to specific branches on the metazoan tree. Together, our analyses suggest that a single, beta-catenin gene fulfilled both adhesion and signaling functions in the last common ancestor of metazoans some 700 million years ago. In addition, gene duplications facilitated the evolution of beta-catenins with novel functions and allowed the evolution of multiple, single-function proteins (cell adhesion or wnt-signaling) from the ancestral, dual-function protein. Integrative methods such as those we have applied here, utilizing the 'natural experiments' present in animal diversity, can be employed to identify novel and shared functional motifs and residues in virtually any protein among the proteomes of model systems and humans.     

Fold change in evolution of protein structures

Typically, protein spatial structures are more conserved in evolution than amino acid sequences. However, the recent explosion of sequence and structure information accompanied by the development of powerful computational methods led to the accumulation of examples of homologous proteins with globally distinct structures. Significant sequence conservation, local structural resemblance, and functional similarity strongly indicate evolutionary relationships between these proteins despite pronounced structural differences at the fold level. Several mechanisms such as insertions/deletions/substitutions, circular permutations, and rearrangements in beta-sheet topologies account for the majority of detected structural irregularities. The existence of evolutionarily related proteins that possess different folds brings new challenges to the homology modeling techniques and the structure classification strategies and offers new opportunities for protein design in experimental studies.     

DOM-fold: a structure with crossing loops found in DmpA, ornithine acetyltransferase, and molybdenum cofactor-binding domain

Understanding relationships between sequence, structure, and evolution is important for functional characterization of proteins. Here, we define a novel DOM-fold as a consensus structure of the domains in DmpA (L-aminopeptidase D-Ala-esterase/amidase), OAT (ornithine acetyltransferase), and MocoBD (molybdenum cofactor-binding domain), and discuss possible evolutionary scenarios of its origin. As shown by a comprehensive structure similarity search, DOM-fold distinguished by a two-layered beta/alpha architecture of a particular topology with unusual crossing loops is unique to those three protein families. DmpA and OAT are evolutionarily related as indicated by their sequence, structural, and functional similarities. Structural similarity between the DmpA/OAT superfamily and the MocoBD domains has not been reported before. Contrary to previous reports, we conclude that functional similarities between DmpA/OAT proteins and N-terminal nucleophile (Ntn) hydrolases are convergent and are unlikely to be inherited from a common ancestor.     

Modeling a minimal ribosome based on comparative sequence analysis

We have determined the three-dimensional organization of ribosomal RNAs and proteins essential for minimal ribosome function. Comparative sequence analysis identifies regions of the ribosome that have been evolutionarily conserved, and the spatial organization of conserved domains is determined by mapping these onto structures of the 30S and 50S subunits determined by X-ray crystallography. Several functional domains of the ribosome are conserved in their three-dimensional organization in the Archaea, Bacteria, Eucaryotic nuclear, mitochondria and chloroplast ribosomes. In contrast, other regions from both subunits have shifted their position in three-dimensional space during evolution, including the L11 binding domain and the alpha-sarcin-ricin loop (SRL). We examined conserved bridge interactions between the two ribosomal subunits, giving an indication of which contacts are more significant. The tRNA contacts that are conserved were also determined, highlighting functional interactions as the tRNA moves through the ribosome during protein synthesis. To augment these studies of a large collection of comparative structural models sampled from all major branches on the phylogenetic tree, Caenorhabditis elegans mitochondrial rRNA is considered individually because it is among the smallest rRNA sequences known. The C.elegans model supports the large collection of comparative structure models while providing insight into the evolution of mitochondrial ribosomes.     

Ion binding to calmodulin

Summary Over the past few years calcium has emerged as an important bioregulator. Upon external stimulation, the cell generates a transient Ca2+ increase, which is transformed into a cellular event through a molecular cascade. The first step in this cascade is the binding of calcium to proteins present in the cytosol. These proteins capable of binding Ca²⁺ under physiological conditions all belong to the same evolutionary family that evolved from a common ancestor. However, they strongly differ in the properties of their calcium binding sites. Calmodulin, the ubiquitous calcium binding protein present in all eukaryotic cells, is very close to the ancestor protein, presents four calcium binding sites which bind calcium, magnesium and monovalent ions competitively and is involved in the triggering of cellular processes. Parvalbumin, another member of the family, is more specialized and found mostly in fast-twitch skeletal muscle. It binds calcium and magnesium with high affinity and seems to be involved in muscle relaxation. On the other hand, troponin C which confers Ca²⁺ sensitivity to acto-myosin interaction exhibits both triggering and relaxing sites. The study of intracellular Ca^2– binding proteins has shown that calcium binding proteins have evolved from a simple common structure to fulfill different functions.Abbreviations CaBP calcium-binding protein - ICaBP the vitamin D-dependent intestinal Cat+binding protein - S-100 the glial S-100 protein - RLC the phosphorylatable myosin regulatory light chain - CaM calmodulin - Pa parvalbumin - TnC troponin C - TnI troponin I - Hepes N-2-hydroxyethylpipezarine, N-2-ethane-sulfonic acid - W7 N-(6-Aminohexyl)-5-chloro-l-Naphtalene sulfonamide - SDS sodium dodecyl sulfate - NMR nuclear magnetic resonance     

Modularity in the evolution of yeast protein interaction network

Protein interaction networks are known to exhibit remarkable structures: scale-free and small-world and modular structures. To explain the evolutionary processes of protein interaction networks possessing scale-free and small-world structures, preferential attachment and duplication-divergence models have been proposed as mathematical models. Protein interaction networks are also known to exhibit another remarkable structural characteristic, modular structure. How the protein interaction networks became to exhibit modularity in their evolution? Here, we propose a hypothesis of modularity in the evolution of yeast protein interaction network based on molecular evolutionary evidence. We assigned yeast proteins into six evolutionary ages by constructing a phylogenetic profile. We found that all the almost half of hub proteins are evolutionarily new. Examining the evolutionary processes of protein complexes, functional modules and topological modules, we also found that member proteins of these modules tend to appear in one or two evolutionary ages. Moreover, proteins in protein complexes and topological modules show significantly low evolutionary rates than those not in these modules. Our results suggest a hypothesis of modularity in the evolution of yeast protein interaction network as systems evolution.     

Structural comparison of cupredoxin domains: domain recycling to construct proteins with novel functions.

The three-dimensional structures of the copper-containing enzymes ascorbate oxidase, ceruloplasmin, and nitrite reductase, comprised of multiple domains with a cupredoxin fold, are consistent with having evolved from a common ancestor. The presence or absence of copper sites has complicated ascertaining the structural and evolutionary relationship among these and related proteins. Simultaneous structural superposition of the enzyme domains and their known cupredoxin relatives shows clearly that there are at least six cupredoxin classes, and that the evolution of the conserved core of these domains is independent of the presence or absence of copper sites. Relationships among the variable loops in these structures show that the two-domain ancestor of the blue oxidases contained a trinuclear-copper interface but could not have functioned in a monomeric state. Comparison of the sequence of the copper-containing, iron-regulating protein. Ferrous transport (Fet3) from yeast to the structurally defined core and loop residues of the cupredoxins suggests specific residues that could be involved in the ferroxidase activity of Fet3.     

Protein structure search and local structure characterization

Background  

Structural similarities among proteins can provide valuable insight into their functional mechanisms and relationships. As the number of available three-dimensional (3D) protein structures increases, a greater variety of studies can be conducted with increasing efficiency, among which is the design of protein structural alphabets. Structural alphabets allow us to characterize local structures of proteins and describe the global folding structure of a protein using a one-dimensional (1D) sequence. Thus, 1D sequences can be used to identify structural similarities among proteins using standard sequence alignment tools such as BLAST or FASTA.     

Systematic analysis of short internal indels and their impact on protein folding

Background  

Protein sequence insertions/deletions (indels) can be introduced during evolution or through alternative splicing (AS). Alternative splicing is an important biological phenomenon and is considered as the major means of expanding structural and functional diversity in eukaryotes. Knowledge of the structural changes due to indels is critical to our understanding of the evolution of protein structure and function. In addition, it can help us probe the evolution of alternative splicing and the diversity of functional isoforms. However, little is known about the effects of indels, in particular the ones involving core secondary structures, on the folding of protein structures. The long term goal of our study is to accurately predict the protein AS isoform structures. As a first step towards this goal, we performed a systematic analysis on the structural changes caused by short internal indels through mining highly homologous proteins in Protein Data Bank (PDB).     

Orthologs, paralogs and genome comparisons

During the past decade, ancient gene duplications were recognized as one of the main forces in the generation of diverse gene families and the creation of new functional capabilities. New tools developed to search data banks for homologous sequences, and an increased availability of reliable three-dimensional structural information led to the recognition that proteins with diverse functions can belong to the same superfamily. Analyses of the evolution of these superfamilies promises to provide insights into early evolution but are complicated by several important evolutionary processes. Horizontal transfer of genes can lead to a vertical spread of innovations among organisms, therefore finding a certain property in some descendants of an ancestor does not guarantee that it was present in that ancestor. Complete or partial gene conversion between duplicated genes can yield phylogenetic trees with several, apparently independent gene duplications, suggesting an often surprising parallelism in the evolution of independent lineages. Additionally, the breakup of domains within a protein and the fusion of domains into multifunctional proteins makes the delineation of superfamilies a task that remains difficult to automate.     

Solution Structure and Calcium-Binding Properties of M-Crystallin, A Primordial βγ-Crystallin from Archaea

The lens βγ-crystallin superfamily has many diverse but topologically related members belonging to various taxa. Based on structural topology, these proteins are considered to be evolutionarily related to lens crystallins, suggesting their origin from a common ancestor. Proteins with βγ-crystallin domains, although found in some eukaryotes and eubacteria, have not yet been reported in archaea. Sequence searches in the genome of the archaebacterium Methanosarcina acetivorans revealed the presence of a protein annotated as a βγ-crystallin family protein, named M-crystallin. Solution structure of this protein indicates a typical βγ-crystallin fold with a paired Greek-key motif. Among the known structures of βγ-crystallin members, M-crystallin was found to be structurally similar to the vertebrate lens βγ-crystallins. The Ca^2 +-binding properties of this primordial protein are somewhat more similar to those of vertebrate βγ-crystallins than to those of bacterial homologues. These observations, taken together, suggest that amphibian and vertebrate βγ-crystallin domains are evolutionarily more related to archaeal homologues than to bacterial homologues. Additionally, identification of a βγ-crystallin homologue in archaea allows us to demonstrate the presence of this domain in all the three domains of life.     

Natural history of ABC systems: not only transporters

In recent years, our understanding of the functioning of ABC (ATP-binding cassette) systems has been boosted by the combination of biochemical and structural approaches. However, the origin and the distribution of ABC proteins among living organisms are difficult to understand in a phylogenetic perspective, because it is hard to discriminate orthology and paralogy, due to the existence of horizontal gene transfer. In this chapter, I present an update of the classification of ABC systems and discuss a hypothetical scenario of their evolution. The hypothetical presence of ABC ATPases in the last common ancestor of modern organisms is discussed, as well as the additional possibility that ABC systems might have been transmitted to eukaryotes, after the two endosymbiosis events that led to the constitution of eukaryotic organelles. I update the functional information of selected ABC systems and introduce new families of ABC proteins that have been included recently into this vast superfamily, thanks to the availability of high-resolution three-dimensional structures.     

Phylogenetic divergence of fish and mammalian metallothionein: relationships with structural diversification and organismal temperature

Metallothioneins (MTs) are nonenzymatic low molecular weight proteins, that play an important role in the homeostasis and detoxification of heavy metals in a large variety of organisms. These proteins are endowed with striking features, including an unusual amino acid composition characterized by the presence of 20 cysteines out of a total of 60 residues and absence of secondary structure elements. It is generally accepted that MTs underwent few modifications during evolution because of these structural and functional constraints. Such a conclusion is founded on the studies carried out mostly on MTs of mammalian origin. For such a reason, we have decided to compare the MTs of homeothermic and poikilothermic organisms, such as mammals and fish, with the specific aim to put in relation phylogenetic divergence and structural/functional adaptation to temperature. We have included in our analysis also Antarctic Notothenioids, a fish group characterized by genetic isolation and cold-adaptation to a particular harsh environment. We have determined the average hydropathic index of ancestral MT sequences and used them to infer the temperatures of the environment housing the hypothetical ancestor organisms. Finally, we have derived phylogenetic relationships of MT molecules from the pairwise comparison of their three-dimensional structures.     

Comparison of sequence-based and structure-based phylogenetic trees of homologous proteins: Inferences on protein evolution

Several studies based on the known three-dimensional (3-D) structures of proteins show that two homologous proteins with insignificant sequence similarity could adopt a common fold and may perform same or similar biochemical functions. Hence, it is appropriate to use similarities in 3-D structure of proteins rather than the amino acid sequence similarities in modelling evolution of distantly related proteins. Here we present an assessment of using 3-D structures in modelling evolution of homologous proteins. Using a dataset of 108 protein domain families of known structures with at least 10 members per family we present a comparison of extent of structural and sequence dissimilarities among pairs of proteins which are inputs into the construction of phylogenetic trees. We find that correlation between the structure-based dissimilarity measures and the sequence-based dissimilarity measures is usually good if the sequence similarity among the homologues is about 30% or more. For protein families with low sequence similarity among the members, the correlation coefficient between the sequence-based and the structure-based dissimilarities are poor. In these cases the structure-based dendrogram clusters proteins with most similar biochemical functional properties better than the sequence-similarity based dendrogram. In multi-domain protein families and disulphide-rich protein families the correlation coefficient for the match of sequence-based and structure-based dissimilarity (SDM) measures can be poor though the sequence identity could be higher than 30%. Hence it is suggested that protein evolution is best modelled using 3-D structures if the sequence similarities (SSM) of the homologues are very low.     

Functional and evolutionary relationships of troponin C.

Striated muscle contraction is initiated when, following membrane depolarization, Ca(2+) binds to the low-affinity Ca(2+) binding sites of troponin C (TnC). The Ca(2+) activation of this protein results in a rearrangement of the components (troponin I, troponin T, and tropomyosin) of the thin filament, resulting in increased interaction between actin and myosin and the formation of cross bridges. The functional properties of this protein are therefore critical in determining the active properties of striated muscle. To date there are 61 known TnCs that have been cloned from 41 vertebrate and invertebrate species. In vertebrate species there are also distinct fast skeletal muscle and cardiac TnC proteins. While there is relatively high conservation of the amino acid sequence of TnC homologs between species and tissue types, there is wide variation in the functional properties of these proteins. To date there has been extensive study of the structure and function of this protein and how differences in these translate into the functional properties of muscles. The purpose of this work is to integrate these studies of TnC with phylogenetic analysis to investigate how changes in the sequence and function of this protein, integrate with the evolution of striated muscle.     

Conservation of residue interactions in a family of Ca-binding proteins

In the TNC family of Ca-binding proteins (calmodulin, parvalbumin, intestinal calcium binding protein and troponin C) approximately 70 well-conserved amino acid sequences and six crystal structures are known. We find a clear correlation between residue contacts in the structures and residue conservation in the sequences: residues with strong sidechain-sidechain contacts in the three-dimenesional structure tend to be the more conserved in the sequence. This is one way to quantify the intuitive notion of the importance of sidechain interactions for maintaining protein three-dimensional structure in evolution and may usefully be taken into account in planning point mutations in protein engineering.     

Functional architecture of the CFTR chloride channel

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ATP-binding cassette (ABC) family of membrane transport proteins. CFTR is unique among ABC proteins in that it functions not as an active transporter but as an ATP-gated Cl^? channel. As an ion channel, the function of the CFTR transmembrane channel pore that mediates Cl^? movement has been studied in great detail. On the other hand, only low resolution structural data is available on the transmembrane parts of the protein. The structure of the channel pore has, however, been modeled on the known structure of active transporter ABC proteins. Currently, significant barriers exist to building a unified view of CFTR pore structure and function. Reconciling functional data on the channel with indirect structural data based on other proteins with very different transport functions and substrates has proven problematic. This review summarizes current structural and functional models of the CFTR Cl^? channel pore, including a comprehensive review of previous electrophysiological investigations of channel structure and function. In addition, functional data on the three-dimensional arrangement of pore-lining helices, as well as contemporary hypotheses concerning conformational changes in the pore that occur during channel opening and closing, are discussed. Important similarities and differences between different models of the pore highlight current gaps in our knowledge of CFTR structure and function. In order to fill these gaps, structural and functional models of the membrane-spanning pore need to become better integrated.     

Collective dynamics differentiates functional divergence in protein evolution

Protein evolution is most commonly studied by analyzing related protein sequences and generating ancestral sequences through Bayesian and Maximum Likelihood methods, and/or by resurrecting ancestral proteins in the lab and performing ligand binding studies to determine function. Structural and dynamic evolution have largely been left out of molecular evolution studies. Here we incorporate both structure and dynamics to elucidate the molecular principles behind the divergence in the evolutionary path of the steroid receptor proteins. We determine the likely structure of three evolutionarily diverged ancestral steroid receptor proteins using the Zipping and Assembly Method with FRODA (ZAMF). Our predictions are within ∼2.7 Å all-atom RMSD of the respective crystal structures of the ancestral steroid receptors. Beyond static structure prediction, a particular feature of ZAMF is that it generates protein dynamics information. We investigate the differences in conformational dynamics of diverged proteins by obtaining the most collective motion through essential dynamics. Strikingly, our analysis shows that evolutionarily diverged proteins of the same family do not share the same dynamic subspace, while those sharing the same function are simultaneously clustered together and distant from those, that have functionally diverged. Dynamic analysis also enables those mutations that most affect dynamics to be identified. It correctly predicts all mutations (functional and permissive) necessary to evolve new function and ∼60% of permissive mutations necessary to recover ancestral function.