Phylogenetic analysis of 70 kD heat shock protein sequences suggests a chimeric origin for the eukaryotic cell nucleus

Background: The evolutionary relationships between archaebacteria, eubacteria and eukaryotic cells are of central importance in biology. The current view is that each of these three groups of organisms constitutes a monophyletic domain, and that eukaryotic cells have evolved from an archaebacterial ancestor. Recent studies on a number of highly conserved protein sequences do not, however, support this view and raise important questions concerning the evolutionary relationships between all extant organisms, particularly regarding the origin of eukaryotic cells.Results We have used sequences of 70 kD heat shock protein (hsp70) — the most conserved protein found to date in all species — to examine the evolutionary relationship between various species. We have obtained two new archaebacterial hsp70 sequences from the species, Thermoplasma acidophilum and Halobacterium cutirubrum. A global comparison of hsp70 sequences, including our two new sequences, shows that all known archaebacterial homologs share a number of sequence signatures with the Gram-positive group of bacteria that are not found in any other prokaryotic or eukaryotic species. In contrast, the eukaryotic homologs are shown to share a number of unique sequence features with the Gram-negative bacteria that are not present in any archaebacteria. Detailed phylogenetic analyses of hsp70 sequences strongly support a specific evolutionary relationship between archaebacteria and Gram-positive bacteria on the one hand, and Gram-negative bacteria and eukaryotes on the other. The phylogenetic analyses also indicate a polyphyletic branching of archaebacteria within the Gram-positive species. The possibility that the observed relationships are due to horizontal gene transfers can be excluded on the basis of sequence characteristics of different groups of homologs.Conclusion Our results do not support the view that archaebacteria constitute a monophyletic domain, but instead suggest a close evolutionary linkage between archaebacteria and Gram-positive bacteria. Furthermore, in contrast to the presently accepted view, eukaryotic hsp70s show a close and specific relationship to those from Gram-negative species. To explain the phylogenies based on different gene sequences, a chimeric model for the origin of the eukaryotic cell nucleus involving fusion between an archaebacterium and a Gram-negative eubacterium is proposed. Several predictions from the chimeric model are discussed.     

Exploration of the Arrest Peptide Sequence Space Reveals Arrest-enhanced Variants

Translational arrest peptides (APs) are short stretches of polypeptides that induce translational stalling when synthesized on a ribosome. Mechanical pulling forces acting on the nascent chain can weaken or even abolish stalling. APs can therefore be used as in vivo force sensors, making it possible to measure the forces that act on a nascent chain during translation with single-residue resolution. It is also possible to score the relative strengths of APs by subjecting them to a given pulling force and ranking them according to stalling efficiency. Using the latter approach, we now report an extensive mutagenesis scan of a strong mutant variant of the Mannheimia succiniciproducens SecM AP and identify mutations that further increase the stalling efficiency. Combining three such mutations, we designed an AP that withstands the strongest pulling force we are able to generate at present. We further show that diproline stretches in a nascent protein act as very strong APs when translation is carried out in the absence of elongation factor P. Our findings highlight critical residues in APs, show that certain amino acid sequences induce very strong translational arrest and provide a toolbox of APs of varying strengths that can be used for in vivo force measurements.     

Powerful Sequence Similarity Search Methods and In-Depth Manual Analyses Can Identify Remote Homologs in Many Apparently “Orphan” Viral Proteins

The genome sequences of new viruses often contain many “orphan” or “taxon-specific” proteins apparently lacking homologs. However, because viral proteins evolve very fast, commonly used sequence similarity detection methods such as BLAST may overlook homologs. We analyzed a data set of proteins from RNA viruses characterized as “genus specific” by BLAST. More powerful methods developed recently, such as HHblits or HHpred (available through web-based, user-friendly interfaces), could detect distant homologs of a quarter of these proteins, suggesting that these methods should be used to annotate viral genomes. In-depth manual analyses of a subset of the remaining sequences, guided by contextual information such as taxonomy, gene order, or domain cooccurrence, identified distant homologs of another third. Thus, a combination of powerful automated methods and manual analyses can uncover distant homologs of many proteins thought to be orphans. We expect these methodological results to be also applicable to cellular organisms, since they generally evolve much more slowly than RNA viruses. As an application, we reanalyzed the genome of a bee pathogen, Chronic bee paralysis virus (CBPV). We could identify homologs of most of its proteins thought to be orphans; in each case, identifying homologs provided functional clues. We discovered that CBPV encodes a domain homologous to the Alphavirus methyltransferase-guanylyltransferase; a putative membrane protein, SP24, with homologs in unrelated insect viruses and insect-transmitted plant viruses having different morphologies (cileviruses, higreviruses, blunerviruses, negeviruses); and a putative virion glycoprotein, ORF2, also found in negeviruses. SP24 and ORF2 are probably major structural components of the virions.     

Identification of a SecM segment required for export-coupled release from elongation arrest

SecM in Escherichia coli has two functionally crucial regions. The arrest motif near the C-terminus interacts with the ribosomal exit tunnel to arrest its own translational elongation. The signal sequence at the N-terminus directs the SecM nascent polypeptide to the Sec-mediated export pathway to release the arrested state of translation. Here, we addressed the importance of the central region of SecM. Characterization of internal substitution and deletion mutants revealed that a segment from residue 100 to residue 109 is required for the export-coupled release of the SecM nascent chain from the elongation-arrested state. Thus, the central region of SecM is not just a geometric linker but it participates actively in the regulation of translation arrest.     

Nascent SecM Chain Outside the Ribosome Reinforces Translation Arrest

SecM, a bacterial secretion monitor protein, contains a specific amino acid sequence at its C-terminus, called arrest sequence, which interacts with the ribosomal tunnel and arrests its own translation. The arrest sequence is sufficient and necessary for stable translation arrest. However, some previous studies have suggested that the nascent chain outside the ribosome affects the stability of translation arrest. To clarify this issue, we performed in vitro translation assays with HaloTag proteins fused to the C-terminal fragment of E. coli SecM containing the arrest sequence or the full-length SecM. We showed that the translation of HaloTag proteins, which are fused to the fragment, is not effectively arrested, whereas the translation of HaloTag protein fused to full-length SecM is arrested efficiently. In addition, we observed that the nascent SecM chain outside the ribosome markedly stabilizes the translation arrest. These results indicate that changes in the nascent polypeptide chain outside the ribosome can affect the stability of translation arrest; the nascent SecM chain outside the ribosome stabilizes the translation arrest.     

Emergence of protein fold families through rational design

Diverse proteins with similar structures are grouped into families of homologs and analogs, if their sequence similarity is higher or lower, respectively, than 20%–30%. It was suggested that protein homologs and analogs originate from a common ancestor and diverge in their distinct evolutionary time scales, emerging as a consequence of the physical properties of the protein sequence space. Although a number of studies have determined key signatures of protein family organization, the sequence-structure factors that differentiate the two evolution-related protein families remain unknown. Here, we stipulate that subtle structural changes, which appear due to accumulating mutations in the homologous families, lead to distinct packing of the protein core and, thus, novel compositions of core residues. The latter process leads to the formation of distinct families of homologs. We propose that such differentiation results in the formation of analogous families. To test our postulate, we developed a molecular modeling and design toolkit, Medusa, to computationally design protein sequences that correspond to the same fold family. We find that analogous proteins emerge when a backbone structure deviates only 1–2 Å root-mean-square deviation from the original structure. For close homologs, core residues are highly conserved. However, when the overall sequence similarity drops to ~25%–30%, the composition of core residues starts to diverge, thereby forming novel families of protein homologs. This direct observation of the formation of protein homologs within a specific fold family supports our hypothesis. The conservation of amino acids in designed sequences recapitulates that of the naturally occurring sequences, thereby validating our computational design methodology.     

The ribosomal exit tunnel functions as a discriminating gate

Translation of SecM stalls unless its N-terminal part is "pulled" by the protein export machinery. Here we show that the sequence motif FXXXXWIXXXXGIRAGP that includes a specific arrest point (Pro) causes elongation arrest within the ribosome. Mutations that bypass the elongation arrest were isolated in 23S rRNA and L22 r protein. Such suppressor mutations occurred at a few specific residues of these components, which all face the narrowest constriction of the ribosomal exit tunnel. Thus, we suggest that this region of the exit tunnel interacts with nascent translation products and functions as a discriminating gate.     

Electrostatics in the ribosomal tunnel modulate chain elongation rates

Electrostatic potentials along the ribosomal exit tunnel are nonuniform and negative. The significance of electrostatics in the tunnel remains relatively uninvestigated, yet they are likely to play a role in translation and secondary folding of nascent peptides. To probe the role of nascent peptide charges in ribosome function, we used a molecular tape measure that was engineered to contain different numbers of charged amino acids localized to known regions of the tunnel and measured chain elongation rates. Positively charged arginine or lysine sequences produce transient arrest (pausing) before the nascent peptide is fully elongated. The rate of conversion from transiently arrested to full-length nascent peptide is faster for peptides containing neutral or negatively charged residues than for those containing positively charged residues. We provide experimental evidence that extraribosomal mechanisms do not account for this charge-specific pausing. We conclude that pausing is due to charge-specific interactions between the tunnel and the nascent peptide.     

Staphylococcus aureus Sortase A-Mediated Incorporation of Peptides: Effect of Peptide Modification on Incorporation

The endogenous Staphylococcus aureus sortase A (SrtA) transpeptidase covalently anchors cell wall-anchored (CWA) proteins equipped with a specific recognition motif (LPXTG) into the peptidoglycan layer of the staphylococcal cell wall. Previous in situ experiments have shown that SrtA is also able to incorporate exogenous, fluorescently labelled, synthetic substrates equipped with the LPXTG motif (K(FITC)LPETG-amide) into the bacterial cell wall, albeit at high concentrations of 500 μM to 1 mM. In the present study, we have evaluated the effect of substrate modification on the incorporation efficiency. This revealed that (i) by elongation of LPETG-amide with a sequence of positively charged amino acids, derived from the C-terminal domain of physiological SrtA substrates, the incorporation efficiency was increased by 20-fold at 10 μM, 100 μM and 250 μM; (ii) Substituting aspartic acid (E) for methionine increased the incorporation of the resulting K(FITC)LPMTG-amide approximately three times at all concentrations tested; (iii) conjugation of the lipid II binding antibiotic vancomycin to K(FITC)LPMTG-amide resulted in the same incorporation levels as K(FITC)LPETG-amide, but much more efficient at an impressive 500-fold lower substrate concentration. These newly developed synthetic substrates can potentially find broad applications in for example the in situ imaging of bacteria; the incorporation of antibody recruiting moieties; the targeted delivery and covalent incorporation of antimicrobial compounds into the bacterial cell wall.     

Endosymbiotic Bacteroidales Bacteria of the Flagellated Protist Pseudotrichonympha grassii in the Gut of the Termite Coptotermes formosanus

A unique lineage of bacteria belonging to the order Bacteroidales was identified as an intracellular endosymbiont of the protist Pseudotrichonympha grassii (Parabasalia, Hypermastigea) in the gut of the termite Coptotermes formosanus. We identified the 16S rRNA, gyrB, elongation factor Tu, and groEL gene sequences in the endosymbiont and detected a very low level of sequence divergence (<0.9% of the nucleotides) in the endosymbiont population within and among protist cells. The Bacteroidales endosymbiont sequence was affiliated with a cluster comprising only sequences from termite gut bacteria and was not closely related to sequences identified for members of the Bacteroidales attached to the cell surfaces of other gut protists. Transmission electron microscopy showed that there were numerous rod-shaped bacteria in the cytoplasm of the host protist, and we detected the endosymbiont by fluorescence in situ hybridization (FISH) with an oligonucleotide probe specific for the 16S rRNA gene identified. Quantification of the abundance of the Bacteroidales endosymbiont by sequence-specific cleavage of rRNA with RNase H and FISH cell counting revealed, surprisingly, that the endosymbiont accounted for 82% of the total bacterial rRNA and 71% of the total bacterial cells in the gut community. The genetically nearly homogeneous endosymbionts of Pseudotrichonympha were very abundant in the gut symbiotic community of the termite.     

Proteogenomic and functional analysis of chromate reduction in Acidiphilium cryptum JF-5, an Fe(III)-respiring acidophile

Acidiphilium cryptum JF-5, an acidophilic iron-respiring Alphaproteobacterium, has the ability to reduce chromate under aerobic and anaerobic conditions, making it an intriguing and useful model organism for the study of extremophilic bacteria in bioremediation applications. Genome sequence annotation suggested two potential mechanisms of Cr(VI) reduction, namely, a number of c-type cytochromes, and a predicted NADPH-dependent Cr(VI) reductase. In laboratory studies using pure cultures of JF-5, an NADPH-dependent chromate reductase activity was detected primarily in soluble protein fractions, and a periplasmic c-type cytochrome (ApcA) was also present, representing two potential means of Cr(VI) reduction. Upon further examination, it was determined that the NADPH-dependent activity was not specific for Cr(VI), and the predicted proteins were not detected in Cr(VI)-grown cultures. Proteomic data did show measureable amounts of ApcA in cells grown with Cr(VI). Purified ApcA is reducible by menadiol, and in turn can reduce Cr(VI), suggesting a means to obtain electrons from the respiratory chain and divert them to Cr(VI). Electrochemical measurements confirm that Cr reduction by ApcA is pH dependent, with low pH being favored. Homology modeling of ApcA and comparison to a known Cr(VI)-reducing c-type cytochrome structure revealed basic amino acids which could interact with chromate ion. From these studies, it can be concluded that A. cryptum has the physiologic and genomic capability to reduce Cr(VI) to the less toxic Cr(III). However, the expected chromate reductase mechanism may not be the primary means of Cr(VI) reduction in this organism.     

Co-translational protein targeting catalyzed by the Escherichia coli signal recognition particle and its receptor.

The Ffh-4.5S ribonucleoprotein particle (RNP) and FtsY from Escherichia coli are homologous to essential components of the mammalian signal recognition particle (SRP) and SRP receptor, respectively. The ability of these E. coli components to function in a bona fide co-translational targeting pathway remains unclear. Here we demonstrate that the Ffh-4.5S RNP and FtsY can efficiently replace their mammalian counterparts in targeting nascent secretory proteins to microsomal membranes in vitro. Targeting in the heterologous system requires a hydrophobic signal sequence, utilizes GTP and, moreover, occurs co-translationally. Unlike mammalian SRP, however, the Ffh-4.5S RNP is unable to arrest translational elongation, which results in a narrow time window for the ribosome nascent chain to interact productively with the membrane-bound translocation machinery. The highly negatively charged N-terminal domain of FtsY, which is a conserved feature among prokaryotic SRP receptor homologs, is important for translocation and acts to localize the protein to the membrane. Our data illustrate the extreme functional conservation between prokaryotic and eukaryotic SRP and SRP receptors and suggest that the basic mechanism of co-translational protein targeting is conserved between bacteria and mammals.     

Divergent stalling sequences sense and control cellular physiology

Recent studies have identified several amino acid sequences that interact with the ribosomal interior components and arrest their own elongation. Whereas stalling of the inducible class depends on specific low-molecular weight compounds, that of the intrinsic class is released when the nascent chain is transported across or inserted into the membrane. The stalled ribosome alters messenger RNA secondary structure and thereby contributes to regulation of the cis-located target gene expression at different levels. The stalling sequences are divergent but likely to utilize non-uniform nature of the peptide bond formation reactions and are recruited relatively recently to different biological systems, possibly including those to be identified in forthcoming studies.     

Comparative Genomics of DNA Fragments from Six Antarctic Marine Planktonic Bacteria

Six environmental fosmid clones from Antarctic coastal water bacterioplankton were completely sequenced. The genome fragments harbored small-subunit rRNA genes that were between 85 and 91% similar to those of their nearest cultivated relatives. The six fragments span four phyla, including the Gemmatimonadetes, Proteobacteria (α and γ), Bacteroidetes, and high-G+C gram-positive bacteria. Gene-finding and annotation analyses identified 244 total open reading frames. Amino acid comparisons of 123 and 113 Antarctic bacterial amino acid sequences to mesophilic homologs from G+C-specific and SwissProt/UniProt databases, respectively, revealed widespread adaptation to the cold. The most significant changes in these Antarctic bacterial protein sequences included a reduction in salt-bridge-forming residues such as arginine, glutamic acid, and aspartic acid, reduced proline contents, and a reduction in stabilizing hydrophobic clusters. Stretches of disordered amino acids were significantly longer in the Antarctic sequences than in the mesophilic sequences. These characteristics were not specific to any one phylum, COG role category, or G+C content and imply that underlying genotypic and biochemical adaptations to the cold are inherent to life in the permanently subzero Antarctic waters.     

Direct probing of the interaction between the signal sequence of nascent preprolactin and the signal recognition particle by specific cross-linking

We have studied the interaction between the signal sequence of nascent preprolactin and the signal recognition particle (SRP) during the initial events in protein translocation across the endoplasmic reticulum membrane. A new method of affinity labeling was used, whereby lysine residues, carrying the photoreactive group 4-(3-trifluoromethyldiazirino) benzoic acid in their side chains, are incorporated into a protein by means of modified lysyl-tRNA, and cross-linking to the interacting component is induced by irradiation. SRP interacts through its Mr 54,000 polypeptide component with the signal sequences of nascent preprolactin chains containing about 70 residues, and with decreasing affinity with longer chains as well; it causes inhibition of elongation. Binding of SRP is reversible and requires the nascent chain to be bound to a functional ribosome. SRP cross-linked to the signal sequence still inhibits elongation but does not prevent it completely. We conclude that SRP does not block the exit site of the polypeptide chain on the ribosome. The SRP receptor of the endoplasmic reticulum membrane displaces the signal sequence from SRP and, even if SRP is cross-linked, releases elongation arrest.     

Ribosome stalling during translation elongation induces cleavage of mRNA being translated in Escherichia coli

Recently, it has been found that ribosome pausing at stop codons caused by certain nascent peptides induces cleavage of mRNA in Escherichia coli cells (1, 2). The question we addressed in the present study is whether mRNA cleavage occurs when translation elongation is prevented. We focused on a specific peptide sequence (AS17), derived from SecM, that is known to cause elongation arrest. When the crp-crr fusion gene encoding CRP-AS17-IIA(Glc) was expressed, cAMP receptor protein (CRP) proteins truncated around the arrest sequence were efficiently produced, and they were tagged by the transfer-messenger RNA (tmRNA) system. Northern blot analysis revealed that both truncated upstream crp and downstream crr mRNAs were generated along with reduced amounts of the full-length crp-crr mRNA. The truncated crp mRNA dramatically decreased in the presence of tmRNA due to rapid degradation. The 3' ends of truncated crp mRNA correspond well to the C termini of the truncated CRP proteins. We conclude that ribosome stalling by the arrest sequence induces mRNA cleavage near the arrest point, resulting in nonstop mRNAs that are recognized by tmRNA. We propose that the mRNA cleavage induced by ribosome stalling acts in concert with the tmRNA system as a way to ensure quality control of protein synthesis and possibly to regulate the expression of certain genes.     

Genetic Identification of Nascent Peptides That Induce Ribosome Stalling

Several nascent peptides stall ribosomes during their own translation in both prokaryotes and eukaryotes. Leader peptides that induce stalling can regulate downstream gene expression. Interestingly, stalling peptides show little sequence similarity and interact with the ribosome through distinct mechanisms. To explore the scope of regulation by stalling peptides and to better understand the mechanism of stalling, we identified and characterized new examples from random libraries. We created a genetic selection that ties the life of Escherichia coli cells to stalling at a specific site. This selection relies on the natural bacterial system that rescues arrested ribosomes. We altered transfer-messenger RNA, a key component of this rescue system, to direct the completion of a necessary protein if and only if stalling occurs. We identified three classes of stalling peptides: C-terminal Pro residues, SecM-like peptides, and the novel stalling sequence FXXYXIWPP. Like the leader peptides SecM and TnaC, the FXXYXIWPP peptide induces stalling efficiently by inhibiting peptidyl transfer. The nascent peptide exit tunnel and peptidyltransferase center are implicated in this stalling event, although mutations in the ribosome affect stalling on SecM and FXXYXIWPP differently. We conclude that ribosome stalling can be caused by numerous sequences and is more common than previously believed.     

“Protein aggregates” contain RNA and DNA,entrapped by misfolded proteins but largely rescued by slowing translational elongation

All neurodegenerative diseases feature aggregates, which usually contain disease‐specific diagnostic proteins; non‐protein constituents, however, have rarely been explored. Aggregates from SY5Y‐APP_Sw neuroblastoma, a cell model of familial Alzheimer''s disease, were crosslinked and sequences of linked peptides identified. We constructed a normalized “contactome” comprising 11 subnetworks, centered on 24 high‐connectivity hubs. Remarkably, all 24 are nucleic acid‐binding proteins. This led us to isolate and sequence RNA and DNA from Alzheimer''s and control aggregates. RNA fragments were mapped to the human genome by RNA‐seq and DNA by ChIP‐seq. Nearly all aggregate RNA sequences mapped to specific genes, whereas DNA fragments were predominantly intergenic. These nucleic acid mappings are all significantly nonrandom, making an artifactual origin extremely unlikely. RNA (mostly cytoplasmic) exceeded DNA (chiefly nuclear) by twofold to fivefold. RNA fragments recovered from AD tissue were ~1.5‐to 2.5‐fold more abundant than those recovered from control tissue, similar to the increase in protein. Aggregate abundances of specific RNA sequences were strikingly differential between cultured SY5Y‐APP_Sw glioblastoma cells expressing APOE3 vs. APOE4, consistent with APOE4 competition for E‐box/CLEAR motifs. We identified many G‐quadruplex and viral sequences within RNA and DNA of aggregates, suggesting that sequestration of viral genomes may have driven the evolution of disordered nucleic acid‐binding proteins. After RNA‐interference knockdown of the translational‐procession factor EEF2 to suppress translation in SY5Y‐APP_Sw cells, the RNA content of aggregates declined by >90%, while reducing protein content by only 30% and altering DNA content by ≤10%. This implies that cotranslational misfolding of nascent proteins may ensnare polysomes into aggregates, accounting for most of their RNA content.     

Stability of RNA quadruplex in open reading frame determines proteolysis of human estrogen receptor α

mRNAs encodes not only information that determines amino acid sequences but also additional layers of information that regulate the translational processes. Notably, translational halt at specific position caused by rare codons or stable RNA structures is one of the potential factors regulating the protein expressions and structures. In this study, a quadruplex-forming potential (QFP) sequence derived from an open reading frame of human estrogen receptor α (hERα) mRNA was revealed to form parallel G-quadruplex and halt the translation elongation in vitro. Moreover, when the full-length hERα and variants containing synonymous mutations in the QFP sequence were expressed in cells, translation products cleaved at specific site were observed in quantities dependent on the thermodynamic stability of the G-quadruplexes. These results suggest that the G-quadruplex formation in the coding region of the hERα mRNA impacts folding and proteolysis of hERα protein by slowing down or temporarily stalling the translation elongation.