Unstructured regions in elongation factors EF1A from three overkingdom of the living world

A new characteristic for classification of the Living World which based on the ability of amino acid sequences to form unstructured regions that appear as loops in their 3D structure is described. Our approach is in principle different from RNA and protein phylogenies that are based on the alignment of amino acid sequences from different organisms. Introduction of new structural-functional characteristic in itself is of undoubted interest because megataxonomy and macrophylogeny lack features that may be resolve evolutionary relation between different groups of organisms though apparent abundance of such characteristics are present. We used the program FoldUnfold to search for unstructured regions in the elongation factors EF1A. The reliability of loop prediction was checked against five factors whose structure is known from X-ray analysis. In addition to two cross-bridges between three structural domains in the elongation factors, the program predicts extra loops. Not counting the effector loop that is inherent to all factors, there are six. Three (A, B and C) of the six different loops are revealed in the first domain, one loop (D) in the second, and two loops (E and F) in the third domain of the factor, all six of which are never found in the same factor. Signatures of elongation factors for each Superkingdom of the Living World have been found for several dozen typical representatives from each Superkingdom. These signatures lead to the variation of the number of loops and their localization within the factor domains. The obtained data leads us to believe that the approach based on the prediction of unstructured protein loops--up to six--must have higher resolution than the method based on the indels (insertion + deletion), the number of which equals one for the same elongation factors. In our analysis, the specificity of sequences it is important, in addition to the existence of loops. Since the total number of loops predicted in the factors increases with the complexity of an organism, we propose the following about the role of the loops in evolution: holding to the principle of "thrifty inventiveness", Nature operates with different universal inserts (loops) adapting their number and location among the factor domains as well as their amino acid composition so that the protein will perform special functions: one in protozoa and several in higher organisms.     

Evolution of the macroglobulin protein family: from bacteria to primates

Proteins of the macroglobulin family are thioether-containing glycoproteins that act as inhibitors of a wide range of hydrolases, transporters and regulators of cytokine, hormone, lipid and oligonucleotide synthesis. As ancient components of innate immunity, these proteins are involved in folding of endogenous proteins as well as recognition and presentation of exogenous antigens. Interaction of macroglobulins with transmembrane receptors triggers cascades of reactions that regulate energy metabolism, cell division and apoptosis, participate in reproduction and cancerogenesis. A broad spectrum of conformational and functional states of molecules, depending on the type of ligands, and an appropriate set of implemented functions allow us to consider these proteins as key regulators of proteostasis. This review addresses the structure and function of macroglobulin proteins during evolution of organisms staying at different phylogenetic levels.     

Exploring metazoan evolution through dynamic and holistic changes in protein families and domains

ABSTRACT: BACKGROUND: Proteins convey the majority of biochemical and cellular activities in organisms. Over the course of evolution, proteins undergo normal sequence mutations as well as large scale mutations involving domain duplication and/or domain shuffling. These events result in the generation of new proteins and protein families. Processes that affect proteome evolution drive species diversity and adaptation. Herein, change over the course of metazoan evolution, as defined by birth/death and duplication/deletion events within protein families and domains, was examined using the proteomes of 9 metazoan and two outgroup species. RESULTS: In studying members of the three major metazoan groups, the vertebrates, arthropods, and nematodes, we found that the number of protein families increased at the majority of lineages over the course of metazoan evolution where the magnitude of these increases was greatest at the lineages leading to mammals. In contrast, the number of protein domains decreased at most lineages and at all terminal lineages. This resulted in a weak correlation between protein family birth and domain birth; however, the correlation between domain birth and domain member duplication was quite strong. These data suggest that domain birth and protein family birth occur via different mechanisms, and that domain shuffling plays a role in the formation of protein families. The ratio of protein family birth to protein domain birth (domain shuffling index) suggests that shuffling had a more demonstrable effect on protein families in nematodes and arthropods than in vertebrates. Through the contrast of high and low domain shuffling indices at the lineages of Trichinella spiralis and Gallus gallus, we propose a link between protein redundancy and evolutionary changes controlled by domain shuffling; however, the speed of adaptation among the different lineages was relatively invariant. Evaluating the functions of protein families that appeared or disappeared at the last common ancestors (LCAs) of the three metazoan clades supports a correlation with organism adaptation. Furthermore, bursts of new protein families and domains in the LCAs of metazoans and vertebrates are consistent with whole genome duplications. CONCLUSION: Metazoan speciation and adaptation were explored by birth/death and duplication/deletion events among protein families and domains. Our results provide insights into protein evolution and its bearing on metazoan evolution.     

Structural study of the catalytic domain of PKCzeta using infrared spectroscopy and two-dimensional infrared correlation spectroscopy

The secondary structure of the catalytic domain from protein kinase C zeta was studied using IR spectroscopy. In the presence of the substrate MgATP, there was a significant change in the secondary structure. After heating to 80 degrees C, a 14% decrease in the alpha-helix component was observed, accompanied by a 6% decrease in the beta-pleated sheet; no change was observed in the large loops or in 3(10)-helix plus associated loops. The maximum increase with heating was observed in the aggregated beta-sheet component, with an increase of 14%. In the presence of MgATP, and compared with the sample heated in its absence, there was a substantial decrease in the 3(10)-helix plus associated loops and an increase in alpha-helix. Synchronous 2D-IR correlation showed that the main changes occurred at 1617 cm(-1), which was assigned to changes in the intermolecular aggregated beta-sheet of the denaturated protein. This increase was mainly correlated with the change in alpha-helix. In the presence of MgATP, the main correlation was between aggregated beta-sheet and the large loops component. The asynchronous 2D-correlation spectrum indicated that a number of components are transformed in intermolecularly aggregated beta-sheet, especially the alpha-helix and beta-sheet components. It is interesting that changes in 3(10)-helix plus associated loops and in alpha-helix preceded changes in large loops, which suggests that the open loops structure exists as an intermediate state during denaturation. In summary, IR spectroscopy revealed an important effect of MgATP on the secondary structure and on the thermal unfolding process when this was induced, whereas 2D-IR correlation spectroscopy allowed us to show the establishment of the denaturation pathway of this protein.     

Small but mighty: Functional landscape of the versatile geminivirus-encoded C4 protein

The fast-paced evolution of viruses enables them to quickly adapt to the organisms they infect by constantly exploring the potential functional landscape of the proteins encoded in their genomes. Geminiviruses, DNA viruses infecting plants and causing devastating crop diseases worldwide, produce a limited number of multifunctional proteins that mediate the manipulation of the cellular environment to the virus’ advantage. Among the proteins produced by the members of this family, C4, the smallest one described to date, is emerging as a powerful viral effector with unexpected versatility. C4 is the only geminiviral protein consistently subjected to positive selection and displays a number of dynamic subcellular localizations, interacting partners, and functions, which can vary between viral species. In this review, we aim to summarize our current knowledge on this remarkable viral protein, encompassing the different aspects of its multilayered diversity, and discuss what it can teach us about geminivirus evolution, invasion requirements, and virulence strategies.     

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an "activated" form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.     

Zinc proteomes, phylogenetics and evolution

Evolution has not been studied in detail with reference to the changing environment. This requires a study of the inorganic chemistry of organisms, especially metalloproteins. The evolution of organisms has been analysed many times previously using comparative studies, fossils, and molecular sequences of proteins, DNA and 16s rRNA (Zhang and Gladyshev, Chem. Rev., 2009, 109, 4828). These methods have led to the confirmation of Darwin's original proposal that evolution followed from natural selection in a changing environment often pictured as a tree. In all cases, the main tree in its upper later reaches has been well studied but its lower earlier parts are not so well defined. To approach this topic we have treated evolution as due to the intimate combination of the effect of chemical changes in the environment and in the organisms (Williams and da Silva, The Chemistry of Evolution, 2006, Elsevier). The best chemicals to examine are inorganic ions as they are common to both. As a more detailed example of the chemical study of organisms we report in this paper a bioinformatic approach to the characterization of the zinc proteomes. We deduce them from the 821 totally sequenced DNA of organisms available on NCBI, exploiting a published method developed by one of us (Andreini, Bertini and Rosato, Acc. Chem. Res., 2009, 42, 1471). Comparing the derived zinc-finger-containing proteins and zinc hydrolytic enzymes in organisms of different complexity there is a correlation in their changes during evolution related to environmental change.     

Chromatin loops,illegitimate recombination,and genome evolution

Chromosomal rearrangements frequently occur at specific places (“hot spots”) in the genome. These recombination hot spots are usually separated by 50–100 kb regions of DNA that are rarely involved in rearrangements. It is quite likely that there is a correlation between the above‐mentioned distances and the average size of DNA loops fixed at the nuclear matrix. Recent studies have demonstrated that DNA loop anchorage regions can be fairly long and can harbor DNA recombination hot spots. We previously proposed that chromosomal DNA loops may constitute the basic units of genome organization in higher eukaryotes. In this review, we consider recombination between DNA loop anchorage regions as a possible source of genome evolution.     

Protein Interaction Networks

The study of protein interactions is playing an ever increasing role in our attempts to understand cells and diseases on a system-wide level. This article reviews several experimental approaches that are currently being used to measure protein–protein, protein–DNA and gene–gene interactions. These techniques have now been scaled up to produce extensive genome-wide data sets that are providing us with a first glimpse of global interaction networks. Complementing these experimental approaches, several computational methodologies to predict protein interactions are also reviewed. Existing databases that serve as repositories for protein interaction information and how such databases are used to analyze high-throughput data from a pathway perspective is also addressed. Finally, current efforts to combine multiple data types to obtain more accurate and comprehensive models of protein interactions are discussed. It is clear that the evolution of these experimental and computational approaches is rapidly changing our view of biology, and promises to provide us with an unprecedented ability to model cells and organisms at a system-wide level.     

Sequence variation at two eosinophil-associated ribonuclease loci in humans

Host defense against invading pathogens is of great importance to the survival of higher organisms. We have been studying the evolution of mammalian eosinophil-associated ribonucleases (EARs), which are members of the ribonuclease A superfamily with known antipathogen activities. Earlier studies showed that positive selection promoted rapid diversification of paralogous EAR genes in both primates and rodents. Intraspecifically, however, it is unknown whether these genes also have divergent alleles. The recent discovery that the gene repertoire of the EAR family is much larger in rodents than in primates has led us to consider the possibility that primates maintain a large number of polymorphic alleles to compensate for a smaller gene repertoire. Here we present sequences of 2417 nucleotides at the two EAR loci, the eosinophil-derived neurotoxin (EDN, RNase 2) and eosinophil cationic protein (ECP, RNase 3), from >50 human individuals. Our data demonstrate that the nucleotide diversities (0.06-0.11%) at these loci are typical for human nuclear genes, thus permitting us to reject this polymorphism hypothesis. No significant departure from neutrality is noted and no signs of overdominant selection are observed. Similar patterns were observed in a preliminary study of chimpanzees. In summary, our results suggest that the antipathogen functions of the primate EARs are conserved after they are established and that these proteins are not currently undergoing rapid diversification in response to challenge from invading microorganisms.     

Disordered regions in elongation factors EF1A in the three superkingdoms of life

A new criterion proposed for classification of the living world is based on the ability of the protein amino acid sequence to form disordered regions, appearing as loops in the 3D structure. The approach used fundamentally differs from the approaches based on comparisons of certain RNA or protein sequences of different organisms. Introduction of any new structural-functional criterion that could resolve the evolutionary relationships between the main groups of origin organisms is of interest in itself, as megasystematics and macrophylogeny lack informative criteria despite the apparent abundance of molecular characteristics. The specialized program FoldUnfold was used to search for disordered regions in the elongation factors EF1A (EFs). The reliability of loop prediction was verified against five EFs with the structures known from X-ray analysis. It was demonstrated with the example of several dozens of typical representatives of the living world that the program predicts extra loops in addition to two linkers between three structural domains in EFs. Besides the effector loop, contained in all EFs, six loops were detected at maximum. Of them, three loops (A, B, and C) are in domain I, one (D) is in domain II, and two (E and F) are in domain III. Moreover, all six loops are never present in the same EF. The EF signatures were determined for each of the superkingdoms of life. Each superkingdom displayed variations in the number of loops and their location within the EF domains. Not only the presence of a particular loop was important in the analysis, but also the specificity of its amino acid sequence. As the total number of predicted loops in EFs increases with the increasing complexity of organisms, the following evolutionary role was postulated for the loops. Following the principle of thrifty inventiveness, nature operates with different universal inserts (loops), adapting their number, location within the EF domains, and amino acid composition so that the protein performs specialized functions—single in protozoa and several in higher organisms.     

Dynameomics: large-scale assessment of native protein flexibility

Structure is only the first step in understanding the interactions and functions of proteins. In this paper, we explore the flexibility of proteins across a broad database of over 250 solvated protein molecular dynamics simulations in water for an aggregate simulation time of approximately 6 micros. These simulations are from our Dynameomics project, and these proteins represent approximately 75% of all known protein structures. We employ principal component analysis of the atomic coordinates over time to determine the primary axis and magnitude of the flexibility of each atom in a simulation. This technique gives us both a database of flexibility for many protein fold families and a compact visual representation of a particular protein's native-state conformational space, neither of which are available using experimental methods alone. These tools allow us to better understand the nature of protein motion and to describe its relationship to other structural and dynamical characteristics. In addition to reporting general properties of protein flexibility and detailing many dynamic motifs, we characterize the relationship between protein native-state flexibility and early events in thermal unfolding and show that flexibility predicts how a protein will begin to unfold. We provide evidence that fold families have conserved flexibility patterns, and family members who deviate from the conserved patterns have very low sequence identity. Finally, we examine novel aspects of highly inflexible loops that are as important to structural integrity as conventional secondary structure. These loops, which are difficult if not impossible to locate without dynamic data, may constitute new structural motifs.     

Comparative study of ammonium transporters in different organisms by study of a large number of structural protein features via data mining algorithms

Ammonium is an excellent nitrogen source, and ammonium transfer is a fundamental process in most organisms. Membrane transport of ammonium is the key component of nitrogen metabolism mediated by Ammonium Transporter/Methylamine Permease/Rhesus (AMT/MEP/Rh) protein family. Ammonium transporters play different physiological roles in various organisms. Here, we looked at the protein characteristics of ammonium transporters in different organisms to create a link between protein characteristics and the organism. In order to increase the accuracy and precision of the employed models, for the first time, an attempt was made to cover all structural aspects of ammonium transporters in animals, bacteria, fungi, plants, and human by extracting and calculating 874 protein attributes of primary, secondary, and tertiary structures for each ammonium transporter. Then, various weighting and modeling algorithms were applied to determine how structural protein features change between organisms. Considering a large number of protein attributes made it possible to detect key protein characteristics in the structure of ammonium transporters. The results, for the first time, indicated that His-based features including count/frequency of His and frequency/count of Ile-His were the most significant features generating different types of ammonium transporters within organisms. Within different tested models, the C5.0 model was the most efficient and precise model for discrimination of organism type, based on ammonium transporter sequence, with the precision of 94.85%. The determination of protein characteristics of ammonium transporters in different organisms provides a new vista for understanding the evolution of transporters based on the modulation of protein characteristics and facilitates engineering of new transporters. In our point of view, dissecting a large number of structural protein characteristics through data mining algorithms provides a novel functional strategy for studying evolution and phylogeny. This research will serve as a basis for future studies on engineering novel ammonium transporters.     

Protein interaction networks

The study of protein interactions is playing an ever increasing role in our attempts to understand cells and diseases on a system-wide level. This article reviews several experimental approaches that are currently being used to measure protein-protein, protein-DNA and gene-gene interactions. These techniques have now been scaled up to produce extensive genome-wide data sets that are providing us with a first glimpse of global interaction networks. Complementing these experimental approaches, several computational methodologies to predict protein interactions are also reviewed. Existing databases that serve as repositories for protein interaction information and how such databases are used to analyze high-throughput data from a pathway perspective is also addressed. Finally, current efforts to combine multiple data types to obtain more accurate and comprehensive models of protein interactions are discussed. It is clear that the evolution of these experimental and computational approaches is rapidly changing our view of biology, and promises to provide us with an unprecedented ability to model cells and organisms at a system-wide level.     

The evolution of RNA editing and pentatricopeptide repeat genes

The pentatricopeptide repeat (PPR) is a degenerate 35-amino-acid structural motif identified from analysis of the sequenced genome of the model plant Arabidopsis thaliana. From the wealth of sequence information now available from plant genomes, the PPR protein family is now known to be one of the largest families in angiosperm species, as most genomes encode 400-600 members. As the number of PPR genes is generally only c. 10-20 in other eukaryotic organisms, including green algae, the family has obviously greatly expanded during land plant evolution. This provides a rare opportunity to study selection pressures driving a 50-fold expansion of a single gene family. PPR proteins are sequence-specific RNA-binding proteins involved in many aspects of RNA processing in organelles. In this review, we will summarize our current knowledge about the evolution of PPR genes, and will discuss the relevance of the dramatic expansion in the family to the functional diversification of plant organelles, focusing primarily on RNA editing.     

Stealth proteins: in silico identification of a novel protein family rendering bacterial pathogens invisible to host immune defense

There are a variety of bacterial defense strategies to survive in a hostile environment. Generation of extracellular polysaccharides has proved to be a simple but effective strategy against the host's innate immune system. A comparative genomics approach led us to identify a new protein family termed Stealth, most likely involved in the synthesis of extracellular polysaccharides. This protein family is characterized by a series of domains conserved across phylogeny from bacteria to eukaryotes. In bacteria, Stealth (previously characterized as SacB, XcbA, or WefC) is encoded by subsets of strains mainly colonizing multicellular organisms, with evidence for a protective effect against the host innate immune defense. More specifically, integrating all the available information about Stealth proteins in bacteria, we propose that Stealth is a D-hexose-1-phosphoryl transferase involved in the synthesis of polysaccharides. In the animal kingdom, Stealth is strongly conserved across evolution from social amoebas to simple and complex multicellular organisms, such as Dictyostelium discoideum, hydra, and human. Based on the occurrence of Stealth in most Eukaryotes and a subset of Prokaryotes together with its potential role in extracellular polysaccharide synthesis, we propose that metazoan Stealth functions to regulate the innate immune system. Moreover, there is good reason to speculate that the acquisition and spread of Stealth could be responsible for future epidemic outbreaks of infectious diseases caused by a large variety of eubacterial pathogens. Our in silico identification of a homologous protein in the human host will help to elucidate the causes of Stealth-dependent virulence. At a more basic level, the characterization of the molecular and cellular function of Stealth proteins may shed light on fundamental mechanisms of innate immune defense against microbial invasion.     

Chaperonin 60: a paradoxical,evolutionarily conserved protein family with multiple moonlighting functions

Chaperonin 60 is the prototypic molecular chaperone, an essential protein in eukaryotes and prokaryotes, whose sequence conservation provides an excellent basis for phylogenetic analysis. Escherichia coli chaperonin 60 (GroEL), the prototype of this family of proteins, has an established oligomeric‐structure‐based folding mechanism and a defined population of folding partners. However, there is a growing number of examples of chaperonin 60 proteins whose crystal structures and oligomeric composition are at variance with GroEL, suggesting that additional complexities in the protein‐folding function of this protein should be expected. In addition, many organisms have multiple chaperonin 60 proteins, some of which have lost their protein‐folding ability. It is emerging that this highly conserved protein has evolved a bewildering variety of additional biological functions – known as moonlighting functions – both within the cell and in the extracellular milieu. Indeed, in some organisms, it is these moonlighting functions that have been left after the loss of the protein‐folding activity. This highlights the major paradox in the biology of chaperonin 60. This article reviews the relationship between the folding and non‐folding (moonlighting) activities of the chaperonin 60 family and discusses current knowledge on their molecular evolution focusing on protein domains involved in the non‐folding chaperonin functions in an attempt to understand the emerging biology of this evolutionarily ancient protein family.     

Insights into the molecular evolution of the PDZ/LIM family and identification of a novel conserved protein motif

The PDZ and LIM domain-containing protein family is encoded by a diverse group of genes whose phylogeny has currently not been analyzed. In mammals, ten genes are found that encode both a PDZ- and one or several LIM-domains. These genes are: ALP, RIL, Elfin (CLP36), Mystique, Enigma (LMP-1), Enigma homologue (ENH), ZASP (Cypher, Oracle), LMO7 and the two LIM domain kinases (LIMK1 and LIMK2). As conventional alignment and phylogenetic procedures of full-length sequences fell short of elucidating the evolutionary history of these genes, we started to analyze the PDZ and LIM domain sequences themselves. Using information from most sequenced eukaryotic lineages, our phylogenetic analysis is based on full-length cDNA-, EST-derived- and genomic- PDZ and LIM domain sequences of over 25 species, ranging from yeast to humans. Plant and protozoan homologs were not found. Our phylogenetic analysis identifies a number of domain duplication and rearrangement events, and shows a single convergent event during evolution of the PDZ/LIM family. Further, we describe the separation of the ALP and Enigma subfamilies in lower vertebrates and identify a novel consensus motif, which we call 'ALP-like motif' (AM). This motif is highly-conserved between ALP subfamily proteins of diverse organisms. We used here a combinatorial approach to define the relation of the PDZ and LIM domain encoding genes and to reconstruct their phylogeny. This analysis allowed us to classify the PDZ/LIM family and to suggest a meaningful model for the molecular evolution of the diverse gene architectures found in this multi-domain family.     

Optimal growth temperature of prokaryotes correlates with class II amino acid composition

Partitioning of aminoacyl-tRNA synthetases and their associated amino acids into two classes allows us to distinguish between thermophilic and mesophilic species based only on amino acids composition. The CLASSDB program has been developed for amino acid content analysis in organisms treated individually or pooled together to form a pattern of characteristic properties. A strong correlation has been observed between optimal growth temperature (OGT) of organisms and class II amino acids content. Amino acid composition in organisms closely related phylogenetically but dissimilar in their OGT testifies that thermo-adaptation happens rather rapidly on the time scale of evolution.     

The conserved regulation of mitochondrial uncoupling proteins: From unicellular eukaryotes to mammals

Uncoupling proteins (UCPs) belong to the mitochondrial anion carrier protein family and mediate regulated proton leak across the inner mitochondrial membrane. Free fatty acids, aldehydes such as hydroxynonenal, and retinoids activate UCPs. However, there are some controversies about the effective action of retinoids and aldehydes alone; thus, only free fatty acids are commonly accepted positive effectors of UCPs. Purine nucleotides such as GTP inhibit UCP-mediated mitochondrial proton leak. In turn, membranous coenzyme Q may play a role as a redox state-dependent metabolic sensor that modulates the complete activation/inhibition of UCPs. Such regulation has been observed for UCPs in microorganisms, plant and animal UCP1 homologues, and UCP1 in mammalian brown adipose tissue. The origin of UCPs is still under debate, but UCP homologues have been identified in all systematic groups of eukaryotes. Despite the differing levels of amino acid/DNA sequence similarities, functional studies in unicellular and multicellular organisms, from amoebae to mammals, suggest that the mechanistic regulation of UCP activity is evolutionarily well conserved. This review focuses on the regulatory feedback loops of UCPs involving free fatty acids, aldehydes, retinoids, purine nucleotides, and coenzyme Q (particularly its reduction level), which may derive from the early stages of evolution as UCP first emerged.