Proline codon pair selection determines ribosome pausing strength and translation efficiency in bacteria

Synonymous and non-synonymous codon substitutions can alleviate dependence on GroEL for folding

The Escherichia coli GroEL/ES chaperonin system facilitates protein folding in an ATP-driven manner. There are <100 obligate clients of this system in E. coli although GroEL can interact and assist the folding of a multitude of proteins in vitro. It has remained unclear, however, which features distinguish obligate clients from all the other proteins in an E. coli cell. To address this question, we established a system for selecting mutations in mouse dihydrofolate reductase (mDHFR), a GroEL interactor, that diminish its dependence on GroEL for folding. Strikingly, both synonymous and non-synonymous codon substitutions were found to reduce mDHFR's dependence on GroEL. The non-synonymous substitutions increase the rate of spontaneous folding whereas computational analysis indicates that the synonymous substitutions appear to affect translation rates at specific sites.

Revealing Molecular Patterns of Alzheimer's Disease Risk Gene Expression Signatures in COVID-19 Brains

Background

Various virus infections are known to predispose to Alzheimer's disease (AD), and a linkage between COVID-19 and AD has been established. COVID-19 infection modulates the gene expression of the genes implicated in progression of AD.

Objective

Determination of molecular patterns and codon usage and context analysis for the genes that are modulated during COVID-19 infection and are implicated in AD was the target of the study.

Methods

Our study employed a comprehensive array of research methods, including relative synonymous codon usage, Codon adaptation index analysis, Neutrality and parity analysis, Rare codon analyses, and codon context analysis. This meticulous approach was crucial in determining the molecular patterns present in genes up or downregulated during COVID-19 infection.

Results

G/C ending codons were preferred in upregulated genes while not in downregulated genes, and in both gene sets, longer genes have high expressivity. Similarly, T over A nucleotide was preferred, and selection was the major evolutionary force in shaping codon usage in both gene sets. Apart from stops codons, codons CGU - Arg, AUA - Ile, UUA - Leu, UCG - Ser, GUA - Val, and CGA - Arg in upregulated genes, while CUA - Leu, UCG - Ser, and UUA - Leu in downregulated genes were present below the 0.5%. Glutamine-initiated codon pairs have high residual values in upregulated genes. Identical codon pairs GAG-GAG and GUG-GUG were preferred in both gene sets.

Conclusions

The shared and unique molecular features in the up- and downregulated gene sets provide insights into the complex interplay between COVID-19 infection and AD. Further studies are required to elucidate the relationship of these molecular patterns with AD pathology.

Evolving precision: rRNA expansion segment 7S modulates translation velocity and accuracy in eukaryal ribosomes

Ribosome-enhanced translational miscoding of the genetic code causes protein dysfunction and loss of cellular fitness. During evolution, open reading frame length increased, necessitating mechanisms for enhanced translation fidelity. Indeed, eukaryal ribosomes are more accurate than bacterial counterparts, despite their virtually identical, conserved active centers. During the evolution of eukaryotic organisms ribosome expansions at the rRNA and protein level occurred, which potentially increases the options for translation regulation and cotranslational events. Here we tested the hypothesis that ribosomal RNA expansions can modulate the core function of the ribosome, faithful protein synthesis. We demonstrate that a short expansion segment present in all eukaryotes' small subunit, ES7S, is crucial for accurate protein synthesis as its presence adjusts codon-specific velocities and guarantees high levels of cognate tRNA selection. Deletion of ES7S in yeast enhances mistranslation and causes protein destabilization and aggregation, dramatically reducing cellular fitness. Removal of ES7S did not alter ribosome architecture but altered the structural dynamics of inter-subunit bridges thus affecting A-tRNA selection. Exchanging the yeast ES7S sequence with the human ES7S increases accuracy whereas shortening causes the opposite effect. Our study demonstrates that ES7S provided eukaryal ribosomes with higher accuracy without perturbing the structurally conserved decoding center.

Assessment of the broad-spectrum host targeting antiviral efficacy of halofuginone hydrobromide in human airway, intestinal and brain organotypic models

Halofuginone hydrobromide has shown potent antiviral efficacy against a variety of viruses such as SARS-CoV-2, dengue, or chikungunya virus, and has, therefore, been hypothesized to have broad-spectrum antiviral activity. In this paper, we tested this broad-spectrum antiviral activity of Halofuginone hydrobomide against viruses from different families (Picornaviridae, Herpesviridae, Orthomyxoviridae, Coronaviridae, and Flaviviridae). To this end, we used relevant human models of the airway and intestinal epithelium and regionalized neural organoids. Halofuginone hydrobomide showed antiviral activity against SARS-CoV-2 in the airway epithelium with no toxicity at equivalent concentrations used in human clinical trials but not against any of the other tested viruses.

VeloPro: A pipeline integrating Ribo-seq and AlphaFold deciphers association patterns between translation velocity and protein structure features

VeloPro integrates Ribo-seq data and AlphaFold2-predicted 3D protein structure information for characterization of the association patterns between translation velocity and many protein structure features in prokaryotic and eukaryotic organisms across different taxonomical clades such as bacteria, fungi, protozoa, nematode, plants, insect, and mammals. We illustrated that association patterns between translation velocity and protein structure features differ across organisms, partially reflecting their taxonomical relationship.

Domains and residues of the Saccharomyces cerevisiae hnRNP protein Hrp1 important for transcriptional autoregulation and noncoding RNA termination

Proteins that bind the nascent transcript exiting RNA polymerase II can regulate transcription elongation. The essential Saccharomyces cerevisiae hnRNP protein Hrp1 is one such protein and participates in both cleavage and polyadenylation-coupled and Nrd1-Nab3-Sen1-dependent RNA polymerase II termination. Prior evidence that Hrp1 is a positive RNA polymerase II elongation factor suggests that its release from the elongation complex promotes termination. Here we report the effects of deletions and substitutions in Hrp1 on its autoregulation via an Nrd1-Nab3-Sen1-dependent transcription attenuator in the 5'-UTR of its mRNA and on the function of an Hrp1-dependent Nrd1-Nab3-Sen1 terminator in the SNR82 snoRNA gene. Deletion of either of two central RNA recognition motifs or either of the flanking low-sequence complexity domains is lethal. Smaller, viable deletions in the amino-terminal low-sequence complexity domain cause readthrough of both the HRP1 attenuator and SNR82 terminator. Substitutions that cause readthrough localized mostly to the RNA recognition motifs, although not always to the RNA-binding face. We found that autoregulation of Hrp1 mRNA synthesis is surprisingly robust, overcoming the expected lethal effects of the start codon and frameshift mutations via overexpression of the mRNA up to 40-fold. Our results suggest a model in which binding of attenuator or terminator elements in the nascent transcript by RNA recognition motifs 1 and 2 disrupts interactions between RNA recognition motif 2 and the RNA polymerase II elongation complex, increasing its susceptibility to termination.

Analysis of 3.5 million SARS-CoV-2 sequences reveals unique mutational trends with consistent nucleotide and codon frequencies

BACKGROUND

Since the onset of the SARS-CoV-2 pandemic, bioinformatic analyses have been performed to understand the nucleotide and synonymous codon usage features and mutational patterns of the virus. However, comparatively few have attempted to perform such analyses on a considerably large cohort of viral genomes while organizing the plethora of available sequence data for a month-by-month analysis to observe changes over time. Here, we aimed to perform sequence composition and mutation analysis of SARS-CoV-2, separating sequences by gene, clade, and timepoints, and contrast the mutational profile of SARS-CoV-2 to other comparable RNA viruses.

METHODS

Using a cleaned, filtered, and pre-aligned dataset of over 3.5 million sequences downloaded from the GISAID database, we computed nucleotide and codon usage statistics, including calculation of relative synonymous codon usage values. We then calculated codon adaptation index (CAI) changes and a nonsynonymous/synonymous mutation ratio (dN/dS) over time for our dataset. Finally, we compiled information on the types of mutations occurring for SARS-CoV-2 and other comparable RNA viruses, and generated heatmaps showing codon and nucleotide composition at high entropy positions along the Spike sequence.

RESULTS

We show that nucleotide and codon usage metrics remain relatively consistent over the 32-month span, though there are significant differences between clades within each gene at various timepoints. CAI and dN/dS values vary substantially between different timepoints and different genes, with Spike gene on average showing both the highest CAI and dN/dS values. Mutational analysis showed that SARS-CoV-2 Spike has a higher proportion of nonsynonymous mutations than analogous genes in other RNA viruses, with nonsynonymous mutations outnumbering synonymous ones by up to 20:1. However, at several specific positions, synonymous mutations were overwhelmingly predominant.

CONCLUSIONS

Our multifaceted analysis covering both the composition and mutation signature of SARS-CoV-2 gives valuable insight into the nucleotide frequency and codon usage heterogeneity of SARS-CoV-2 over time, and its unique mutational profile compared to other RNA viruses.

Omics profiling identifies the regulatory functions of the MAPK/ERK pathway in nephron progenitor metabolism

Nephron endowment is defined by fetal kidney growth and crucially dictates renal health in adults. Defects in the molecular regulation of nephron progenitors contribute to only a fraction of reduced nephron mass cases, suggesting alternative causative mechanisms. The importance of MAPK/ERK activation in nephron progenitor maintenance has been previously demonstrated, and here, we characterized the metabolic consequences of MAPK/ERK deficiency. Liquid chromatography/mass spectrometry-based metabolomics profiling identified 42 reduced metabolites, of which 26 were supported by in vivo transcriptional changes in MAPK/ERK-deficient nephron progenitors. Among these, mitochondria, ribosome and amino acid metabolism, together with diminished pyruvate and proline metabolism, were the most affected pathways. In vitro cultures of mouse kidneys demonstrated a dosage-specific function for pyruvate in controlling the shape of the ureteric bud tip, a regulatory niche for nephron progenitors. In vivo disruption of proline metabolism caused premature nephron progenitor exhaustion through their accelerated differentiation in pyrroline-5-carboxylate reductases 1 (Pycr1) and 2 (Pycr2) double-knockout kidneys. Pycr1/Pycr2-deficient progenitors showed normal cell survival, indicating no changes in cellular stress. Our results suggest that MAPK/ERK-dependent metabolism functionally participates in nephron progenitor maintenance by monitoring pyruvate and proline biogenesis in developing kidneys.

The Bacterial ClpXP-ClpB Family Is Enriched with RNA-Binding Protein Complexes

In the matrix of bacteria/mitochondria/chloroplasts, Lon acts as the degradation machine for soluble proteins. In stress periods, however, proteostasis and survival depend on the strongly conserved Clp/Hsp100 family. Currently, the targets of ATP-powered unfoldases/disaggregases ClpB and ClpX and of peptidase ClpP heptameric rings are still unclear. Trapping experiments and proteome profiling in multiple organisms triggered confusion, so we analyzed the consistency of ClpP-trap targets in bacteria. We also provide meta-analyses of protein interactions in humans, to elucidate where Clp family members are enriched. Furthermore, meta-analyses of mouse complexomics are provided. Genotype–phenotype correlations confirmed our concept. Trapping, proteome, and complexome data retrieved consistent coaccumulation of CLPXP with GFM1 and TUFM orthologs. CLPX shows broad interaction selectivity encompassing mitochondrial translation elongation, RNA granules, and nucleoids. CLPB preferentially attaches to mitochondrial RNA granules and translation initiation components; CLPP is enriched with them all and associates with release/recycling factors. Mutations in CLPP cause Perrault syndrome, with phenotypes similar to defects in mtDNA/mtRNA. Thus, we propose that CLPB and CLPXP are crucial to counteract misfolded insoluble protein assemblies that contain nucleotides. This insight is relevant to improve ClpP-modulating drugs that block bacterial growth and for the treatment of human infertility, deafness, and neurodegeneration.

Translation stalling proline motifs are enriched in slow-growing, thermophilic, and multicellular bacteria

Rapid bacterial growth depends on the speed at which ribosomes can translate mRNA into proteins. mRNAs that encode successive stretches of proline can cause ribosomes to stall, substantially reducing translation speed. Such stalling is especially detrimental for species that must grow and divide rapidly. Here, we focus on di-prolyl motifs (XXPPX) and ask whether their prevalence varies with growth rate. To find out we conducted a broad survey of such motifs in >3000 bacterial genomes across 35 phyla. Indeed, fast-growing species encode fewer motifs than slow-growing species, especially in highly expressed proteins. We also found many di-prolyl motifs within thermophiles, where prolines can help maintain proteome stability. Moreover, bacteria with complex, multicellular lifecycles also encode many di-prolyl motifs. This is especially evident in the slow-growing phylum Myxococcota. Bacteria in this phylum encode many serine-threonine kinases, and many di-prolyl motifs at potential phosphorylation sites within these kinases. Serine-threonine kinases are involved in cell signaling and help regulate developmental processes linked to multicellularity in the Myxococcota. Altogether, our observations suggest that weakened selection on translational rate, whether due to slow or thermophilic growth, may allow di-prolyl motifs to take on new roles in biological processes that are unrelated to translational rate.

Exceptionally versatile take II: post-translational modifications of lysine and their impact on bacterial physiology

Among the 22 proteinogenic amino acids, lysine sticks out due to its unparalleled chemical diversity of post-translational modifications. This results in a wide range of possibilities to influence protein function and hence modulate cellular physiology. Concomitantly, lysine derivatives form a metabolic reservoir that can confer selective advantages to those organisms that can utilize it. In this review, we provide examples of selected lysine modifications and describe their role in bacterial physiology.