Robustness of predictions of extremely thermally stable proteins in ancient organisms

A number of studies have addressed the environmental temperatures experienced by ancient life. Computational studies using a nonhomogeneous evolution model have estimated ancestral G + C contents of ribosomal RNAs and the amino acid compositions of ancestral proteins, generating hypotheses regarding the mesophilic last universal common ancestor. In contrast, our previous study computationally reconstructed ancestral amino acid sequences of nucleoside diphosphate kinases using a homogeneous model and then empirically resurrected the ancestral proteins. The thermal stabilities of these ancestral proteins were equivalent to or greater than those of extant homologous thermophilic proteins, supporting the thermophilic universal ancestor theory. In this study, we reinferred ancestral sequences using a dataset from which hyperthermophilic sequences were excluded. We also reinferred ancestral sequences using a nonhomogeneous evolution model. The newly reconstructed ancestral proteins are still thermally stable, further supporting the hypothesis that the ancient organisms contained thermally stable proteins and therefore that they were thermophilic.     

Amino Acid Compositional Shifts During Streptophyte Transitions to Terrestrial Habitats

Across the streptophyte lineage, which includes charophycean algae and embryophytic plants, there have been at least four independent transitions to the terrestrial habitat. One of these involved the evolution of embryophytes (bryophytes and tracheophytes) from a charophycean ancestor, while others involved the earliest branching lineages, containing the monotypic genera Mesostigma and Chlorokybus, and within the Klebsormidiales and Zygnematales lineages. To overcome heat, water stress, and increased exposure to ultraviolet radiation, which must have accompanied these transitions, adaptive mechanisms would have been required. During periods of dehydration and/or desiccation, proteomes struggle to maintain adequate cytoplasmic solute concentrations. The increased usage of charged amino acids (DEHKR) may be one way of maintaining protein hydration, while increased use of aromatic residues (FHWY) protects proteins and nucleic acids by absorbing damaging UV, with both groups of residues thought to be important for the stabilization of protein structures. To test these hypotheses we examined amino acid sequences of orthologous proteins representing both mitochondrion- and plastid-encoded proteomes across streptophytic lineages. We compared relative differences within categories of amino acid residues and found consistent patterns of amino acid compositional fluxuation in extra-membranous regions that correspond with episodes of terrestrialization: positive change in usage frequency for residues with charged side-chains, and aromatic residues of the light-capturing chloroplast proteomes. We also found a general decrease in the usage frequency of hydrophobic, aliphatic, and small residues. These results suggest that amino acid compositional shifts in extra-membrane regions of plastid and mitochondrial proteins may represent biochemical adaptations that allowed green plants to colonize the land.     

Closed Loops of TIM Barrel Protein Fold

Abstract

The closed loops within the proteins of the TIM-barrel fold family are analyzed and compared sequence- and structure-wise. The size distribution of the closed loops of the TIM-barrels confirms universal preference to the standard size of 25–30 residues. 3D structural RMSD comparisons of the closed loops and presentation of their sequences in binary form suggest that the TIM-barrel proteins are built from descendants of several types of basic closed loop prototypes. Comparison of these prototypes points to a likely common ancestor—the alpha helix containing closed loops of 28 amino acids. The presumed ancestor is characterized by specific binary consensus sequence.     

Closed loops of TIM barrel protein fold

The closed loops within the proteins of the TIM-barrel fold family are analyzed and compared sequence- and structure-wise. The size distribution of the closed loops of the TIM-barrels confirms universal preference to the standard size of 25-30 residues. 3D structural RMSD comparisons of the closed loops and presentation of their sequences in binary form suggest that the TIM-barrel proteins are built from descendants of several types of basic closed loop prototypes. Comparison of these prototypes points to a likely common ancestor--the alpha helix containing closed loops of 28 amino acids. The presumed ancestor is characterized by specific binary consensus sequence.     

In Vitro Evolution Reveals Noncationic Protein–RNA Interaction Mediated by Metal Ions

RNA–peptide/protein interactions have been of utmost importance to life since its earliest forms, reaching even before the last universal common ancestor (LUCA). However, the ancient molecular mechanisms behind this key biological interaction remain enigmatic because extant RNA–protein interactions rely heavily on positively charged and aromatic amino acids that were absent (or heavily under-represented) in the early pre-LUCA evolutionary period. Here, an RNA-binding variant of the ribosomal uL11 C-terminal domain was selected from an approximately 10¹⁰ library of partially randomized sequences, all composed of ten prebiotically plausible canonical amino acids. The selected variant binds to the cognate RNA with a similar overall affinity although it is less structured in the unbound form than the wild-type protein domain. The variant complex association and dissociation are both slower than for the wild-type, implying different mechanistic processes involved. The profile of the wild-type and mutant complex stabilities along with molecular dynamics simulations uncovers qualitative differences in the interaction modes. In the absence of positively charged and aromatic residues, the mutant uL11 domain uses ion bridging (K⁺/Mg²⁺) interactions between the RNA sugar-phosphate backbone and glutamic acid residues as an alternative source of stabilization. This study presents experimental support to provide a new perspective on how early protein–RNA interactions evolved, where the lack of aromatic/basic residues may have been compensated by acidic residues plus metal ions.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Resonant recognition model and protein topography. Model studies with myoglobin, hemoglobin and lysozyme

This study describes the further extension of the resonant recognition model for the analysis and prediction of protein--protein and protein--DNA structure/function dependencies. The model is based on the significant correlation between spectra of numerical presentations of the amino acid or nucleotide sequences of proteins and their coded biological activity. According to this physico-mathematical method, it is possible to define amino acids in the sequence which are predicted to be the most critical for protein function. Using sperm whale myoglobin, human hemoglobin and hen egg white lysozyme as model protein examples, sets of predicted amino acids, or so-called 'hot spots', have been identified within the tertiary structure. It was found for each protein that the predicted 'hot spots', which are distributed along the primary sequence, are spatially grouped in a dome-like arrangement over the active site. The identified amino acids did not correspond to the amino acid residues which are involved in the chemical reaction site of these proteins. It is thus proposed that the resonant recognition model helps to identify amino acid residues which are important for the creation of the molecular structure around the catalytic active site and also the associated physical field conditions required for biorecognition, docking of the specific substrate and full biological activity.     

Selective constraints, amino acid composition, and the rate of protein evolution

What are the major forces governing protein evolution? A common view is that proteins with strong structural and functional requirements evolve more slowly than proteins with weak constraints, because a stringent negative selection pressure limits the number of substitutions. In contrast, Graur claimed that the substitution rate of a protein is mainly determined by its amino acid composition and the changeabilities of amino acids. In this paper, however, we found that the relative changeabilities of amino acids in mammalian proteins are different for transmembranal and nontransmembranal segments, which have very distinct structural requirements. This indicates that the changeability of a given residue is influenced by the structural and functional context. We also reexamined the relationship between substitution rate and amino acid composition. Indeed, the two kinds of segments exhibit contrasting amino acid compositions: transmembranal regions are made up mainly of hydrophobic residues (a total frequency of approximately 60%) and are very poor in polar amino acids (<5%), whereas nontransmembranal segments have frequencies of 30% and 22%, respectively. Interestingly, we found that within a given integral membrane protein, nontransmembranal segments accumulate, on average, twice as many substitutions as transmembranal regions. However, regression analyses showed that the variability in amino acid frequencies among proteins cannot explain more than 30% of the variability in substitution rate for the transmembranal and nontransmembranal data sets. Furthermore, transmembranal and nontransmembranal segments evolving at the same rate in different proteins have different compositions, and the compositions of slowly evolving and rapidly evolving segments of the same type are similar. From these observations, we conclude that the rate of protein evolution is only weakly affected by amino acid composition but is mostly determined by the strength of functional requirements or selective constraints.     

Evolution of amino acid frequencies in proteins over deep time: inferred order of introduction of amino acids into the genetic code

To understand more fully how amino acid composition of proteins has changed over the course of evolution, a method has been developed for estimating the composition of proteins in an ancestral genome. Estimates are based upon the composition of conserved residues in descendant sequences and empirical knowledge of the relative probability of conservation of various amino acids. Simulations are used to model and correct for errors in the estimates. The method was used to infer the amino acid composition of a large protein set in the Last Universal Ancestor (LUA) of all extant species. Relative to the modern protein set, LUA proteins were found to be generally richer in those amino acids that are believed to have been most abundant in the prebiotic environment and poorer in those amino acids that are believed to have been unavailable or scarce. It is proposed that the inferred amino acid composition of proteins in the LUA probably reflects historical events in the establishment of the genetic code.     

Role of mutational bias and natural selection on genome-wide nucleotide bias in prokaryotic organisms

Correlations between genomic GC contents and amino acid frequencies were studied in the homologous sequences of 12 eubacterial genomes. Results show that amino acids encoded by GC-rich codons increases significantly with genomic GC contents, whereas opposite trend was observed in case of amino acids encoded by GC-poor codons. Further studies show all the amino acids do not change in the predicted direction according to their genomic GC pressure, suggesting that protein evolution is not entirely dictated by their nucleotide frequencies. Amino acid substitution matrix calculated among hydrophobic, amphipathic and hydrophilic amino acid groups' shows that amphipathic and hydrophilic amino acids are more frequently substituted by hydrophobic amino acids than from hydrophobic to hydrophilic or amphipathic amino acids. This indicates that nucleotide bias induces a directional changes in proteome composition in such a way that underwent strong changes in hydropathy values. In fact, significant increases in hydrophobicity values have also been observed with the increase of genomic GC contents. Correlations between GC contents and amino acid compositions in three different predicted protein secondary structures show that hydropathy values increases significantly with GC contents in aperiodic and helix structures whereas strand structure remains insensitive with the genomic GC levels. The relative importance of mutation and selection on the evolution of proteins have been discussed on the basis of these results.     

Primordia Vita. Deconvolution from Modern Sequences.

Evolution of the triplet code is reconstructed on the basis of consensus temporal order of appearance of amino acids. Several important predictions are confirmed by computational sequence analyses. The earliest amino acids, alanine and glycine, have been encoded by GCC and GGC codons, as today. They were succeeded, respectively, by A- and G-series of amino acids, encoded by pyrimidine-central and purine-central codons. The length of the earliest proteins is estimated to be 6–7 residues. The earliest mRNAs were short G+C-rich molecules. These short sequences could have formed hairpins. This is confirmed by analysis of modern prokaryotic mRNA sequences. Predominant size of detected ancient hairpins also corresponds to 6–7 amino acids, as above. Vestiges of last common ancestor can be found in extant proteins in form of entirely conserved short sequences of size six to nine residues present in all or almost all sequenced prokaryotic proteomes (omnipresent motifs). The functions of the topmost conserved octamers are not involved in the basic elementary syntheses. This suggests an initial abiotic supply of amino acids, bases and sugars. Presented at: National Workshop on Astrobiology: Search for Life in the Solar System, Capri, Italy, 26 to 28 October, 2005.     

Early fixation of an optimal genetic code

The evolutionary forces that produced the canonical genetic code before the last universal ancestor remain obscure. One hypothesis is that the arrangement of amino acid/codon assignments results from selection to minimize the effects of errors (e.g., mistranslation and mutation) on resulting proteins. If amino acid similarity is measured as polarity, the canonical code does indeed outperform most theoretical alternatives. However, this finding does not hold for other amino acid properties, ignores plausible restrictions on possible code structure, and does not address the naturally occurring nonstandard genetic codes. Finally, other analyses have shown that significantly better code structures are possible. Here, we show that if theoretically possible code structures are limited to reflect plausible biological constraints, and amino acid similarity is quantified using empirical data of substitution frequencies, the canonical code is at or very close to a global optimum for error minimization across plausible parameter space. This result is robust to variation in the methods and assumptions of the analysis. Although significantly better codes do exist under some assumptions, they are extremely rare and thus consistent with reports of an adaptive code: previous analyses which suggest otherwise derive from a misleading metric. However, all extant, naturally occurring, secondarily derived, nonstandard genetic codes do appear less adaptive. The arrangement of amino acid assignments to the codons of the standard genetic code appears to be a direct product of natural selection for a system that minimizes the phenotypic impact of genetic error. Potential criticisms of previous analyses appear to be without substance. That known variants of the standard genetic code appear less adaptive suggests that different evolutionary factors predominated before and after fixation of the canonical code. While the evidence for an adaptive code is clear, the process by which the code achieved this optimization requires further attention.     

The very early evolution of protein translocation across membranes

In this study, we used a computational approach to investigate the early evolutionary history of a system of proteins that, together, embed and translocate other proteins across cell membranes. Cell membranes comprise the basis for cellularity, which is an ancient, fundamental organizing principle shared by all organisms and a key innovation in the evolution of life on Earth. Two related requirements for cellularity are that organisms are able to both embed proteins into membranes and translocate proteins across membranes. One system that accomplishes these tasks is the signal recognition particle (SRP) system, in which the core protein components are the paralogs, FtsY and Ffh. Complementary to the SRP system is the Sec translocation channel, in which the primary channel-forming protein is SecY. We performed phylogenetic analyses that strongly supported prior inferences that FtsY, Ffh, and SecY were all present by the time of the last universal common ancestor of life, the LUCA, and that the ancestor of FtsY and Ffh existed before the LUCA. Further, we combined ancestral sequence reconstruction and protein structure and function prediction to show that the LUCA had an SRP system and Sec translocation channel that were similar to those of extant organisms. We also show that the ancestor of Ffh and FtsY that predated the LUCA was more similar to FtsY than Ffh but could still have comprised a rudimentary protein translocation system on its own. Duplication of the ancestor of FtsY and Ffh facilitated the specialization of FtsY as a membrane bound receptor and Ffh as a cytoplasmic protein that could bind nascent proteins with specific membrane-targeting signal sequences. Finally, we analyzed amino acid frequencies in our ancestral sequence reconstructions to infer that the ancestral Ffh/FtsY protein likely arose prior to or just after the completion of the canonical genetic code. Taken together, our results offer a window into the very early evolutionary history of cellularity.     

生物细胞中颗粒状多聚磷酸盐细胞器的结构与功能

多聚磷酸盐(polyphosphat,poly P)是一种由数十个或上百个磷酸根聚合而成的生物大分子,以颗粒状、胶体状和溶解状等多种状态存在于各类生物细胞中。生物体中的poly P能够通过分解提供能量;鳌合金属离子来调节细胞内渗透压,维持质膜稳定;与蛋白质或DNA结合稳定其结构,减轻细胞应激损伤。颗粒状多聚磷酸盐细胞器主要指细胞中用于贮存颗粒状poly P、金属阳离子以及蛋白质、氨基酸和少量水等物质的细胞器。在寄生虫细胞中颗粒状聚磷细胞器常称为酸性钙体,而细菌或者其他微生物细胞中则称为异染颗粒,但是随着研究的不断深入,发现酸性钙体和异染颗粒都具有相似的结构特征,遂将其统一定义为颗粒状多聚磷酸盐细胞器。颗粒状聚磷细胞器的发现拓展了生物共同祖先(last universal common ancestor,LUCA)的学说,丰富了原核生物细胞器认知,我们相信该细胞器在生命起源、抗环境胁迫、生物互作和代谢调控等方面具有重要功能,在疾病治疗以及磷生物地球化学循环过程中发挥重要作用。     

Nucleotide sequence of the gene for the ferrienterochelin receptor FepA in Escherichia coli. Homology among outer membrane receptors that interact with TonB

We have determined the nucleotide sequence of the Escherichia coli fepA gene, which codes for the outer membrane receptor for ferrienterochelin and colicins B and D. The predicted FepA polypeptide has a molecular weight of 79,908 and consists of 723 amino acids. A 22-amino acid leader or signal peptide preceded the mature protein. With respect to overall composition, hydropathy, net charge and distribution of nonpolar segments, the FepA polypeptide was typical of other E. coli outer membrane proteins, except that FepA contained 2 cysteine residues. Comparison of the deduced amino acid sequence of FepA with that of three other TonB-dependent receptors (BtuB, FhuA, and IutA) revealed only a few regions of sequence homology; one of these included the amino-termini. An amino acid substitution within the conserved amino-terminal region of BtuB resulted in production of a receptor that had normal binding functions but was incapable of energy-dependent transport of vitamin B12. This result suggests that the amino-terminal end of these four polypeptides is involved in interaction with the TonB protein or another step of energy transduction. Three other regions of homology were shared among the four proteins: one near residues 50 to 70, another at about residue 100 to 140, and the last between 20 and 40 amino acid residues from the carboxyl terminus. The function of these three regions remains speculative.     

Identification of amino acid residues required for Ras p21 target activation.

The Krev-1 gene has been shown to suppress ras-mediated transformation in vitro. Both ras and Krev-1 proteins have identical effector domains (ras residues 32 to 40), which are required for biological activity and for the interaction of Ras p21 with Ras GTPase-activating protein (GAP). In this study, five amino acid residues flanking the ras effector domain, which are not conserved with the Krev-1 protein, were shown to be required for normal protein-protein interactions and biological activity. The substitution of Krev-1 p21 residues 26, 27, 30, 31, and 45 with the corresponding amino acid residues from Ras p21 resulted in a Krev-1 protein which had ras function in both mammalian and yeast biological assays. Replacement of these residues in Ras p21 with the corresponding Krev-1 p21 amino acids resulted in ras proteins which were impaired biologically or reduced in their affinity for in vitro GAP binding. Evaluation of these mutant ras proteins have implications for Ras p21-GAP interactions in vivo.     

Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex

It is important to identify hotspot residues that determine protein–protein interactions in interfaces of macromolecular complexes. We have applied a combination of ancestral sequence reconstruction and protein design to identify hotspots within imidazole glycerol phosphate synthase (ImGPS). ImGPS is a key metabolic enzyme complex, which links histidine and de novo purine biosynthesis and consists of the cyclase subunit HisF and the glutaminase subunit HisH. Initial fluorescence titration experiments showed that HisH from Zymomonas mobilis (zmHisH) binds with high affinity to the reconstructed HisF from the last universal common ancestor (LUCA‐HisF) but not to HisF from Pyrobaculum arsenaticum (paHisF), which differ by 103 residues. Subsequent titration experiments with a reconstructed evolutionary intermediate linking LUCA‐HisF and paHisF and inspection of the subunit interface of a contemporary ImGPS allowed us to narrow down the differences crucial for zmHisH binding to nine amino acids of HisF. Homology modeling and in silico mutagenesis studies suggested that at most two of these nine HisF residues are crucial for zmHisH binding. These computational results were verified by experimental site‐directed mutagenesis, which finally enabled us to pinpoint a single amino acid residue in HisF that is decisive for high‐affinity binding of zmHisH. Our work shows that the identification of protein interface hotspots can be very efficient when reconstructed proteins with different binding properties are included in the analysis. Proteins 2017; 85:312–321. © 2016 Wiley Periodicals, Inc.     

Physicochemical characterization and functional analysis of some snake venom toxin proteins and related non-toxin proteins of other chordates

Snake venom contains a diverse array of proteins and polypeptides. Cytotoxins and short neurotoxins are non-enzymatic polypeptide components of snake venom. The three-dimensional structure of cytotoxin and short neurotoxin resembles a three finger appearance of three-finger protein super family. Different family members of three-finger protein super family are employed in diverse biological functions. In this work we analyzed the cytotoxin, short neurotoxin and related non-toxin proteins of other chordates in terms of functional analysis, amino acid compositional (%) profile, number of amino acids, molecular weight, theoretical isoelectric point (pI), number of positively charged and negatively charged amino acid residues, instability index and grand average of hydropathy with the help of different bioinformatical tools. Among all interesting results, profile of amino acid composition (%) depicts that all sequences contain a conserved cysteine amount but differential amount of different amino acid residues which have a family specific pattern. Involvement in different biological functions is one of the driving forces which contribute the vivid amino acid composition profile of these proteins. Different biological system dependent adaptation gives the birth of enriched bio-molecules. Understanding of physicochemical properties of these proteins will help to generate medicinally important therapeutic molecules for betterment of human lives.     

The nature of the last universal ancestor and the root of the tree of life, still open questions.

The nature of the last universal ancestor to all extent cellular organisms and the rooting of the universal tree of life are fundamental questions which can now be addressed by molecular evolutionists. Several scenarios have been proposed during the last years, based on the phylogenies of ribosomal RNA and of duplicated proteins, which suggest that the last universal ancestor was either an RNA progenote or an hyperthermophilic prokaryote. We discuss these hypotheses in the light of new data on the evolution of DNA metabolizing enzymes and of contradictions between different protein phylogenies. We conclude that the last universal ancestor was a member of the DNA world already containing several DNA polymerases and DNA topoisomerases. Furthermore, we criticize current data which suggest that the rooting of the universal tree of life is located in the eubacterial branch and we conclude that both rooting the universal tree and the nature of the last universal ancestor are still open questions.     

Identification of a common hyaluronan binding motif in the hyaluronan binding proteins RHAMM, CD44 and link protein.

We have previously identified two hyaluronan (HA) binding domains in the HA receptor, RHAMM, that occur near the carboxyl-terminus of this protein. We show here that these two HA binding domains are the only HA binding regions in RHAMM, and that they contribute approximately equally to the HA binding ability of this receptor. Mutation of domain II using recombinant polypeptides of RHAMM demonstrates that K423 and R431, spaced seven amino acids apart, are critical for HA binding activity. Domain I contains two sets of two basic amino acids, each spaced seven residues apart, and mutation of these basic amino acids reduced their binding to HA--Sepharose. These results predict that two basic amino acids flanking a seven amino acid stretch [hereafter called B(X7)B] are minimally required for HA binding activity. To assess whether this motif predicts HA binding in the intact RHAMM protein, we mutated all basic amino acids in domains I and II that form part of these motifs using site-directed mutagenesis and prepared fusion protein from the mutated cDNA. The altered RHAMM protein did not bind HA, confirming that the basic amino acids and their spacing are critical for binding. A specific requirement for arginine or lysine residues was identified since mutation of K430, R431 and K432 to histidine residues abolished binding. Clustering of basic amino acids either within or at either end of the motif enhanced HA binding activity while the occurrence of acidic residues between the basic amino acids reduced binding. The B(X7)B motif, in which B is either R or K and X7 contains no acidic residues and at least one basic amino acid, was found in all HA binding proteins molecularly characterized to date. Recombinant techniques were used to generate chimeric proteins containing either the B(X7)B motifs present in CD44 or link protein, with the amino-terminus of RHAMM (amino acids 1-238) that does not bind HA. All chimeric proteins containing the motif bound HA in transblot analyses. Site-directed mutations of these motifs in CD44 sequences abolished HA binding. Collectively, these results predict that the motif of B(X7)B as a minimal binding requirement for HA in RHAMM, CD44 and link protein, and occurs in all HA binding proteins described to date.