The Role of RNA Secondary Structure in Regulation of Gene Expression in Bacteria

Genome-Wide Computational Prediction and Analysis of Noncoding RNAs in Oleidesulfovibrio alaskensis G20

Noncoding RNAs (ncRNAs) play key roles in the regulation of important pathways, including cellular growth, stress management, signaling, and biofilm formation. Sulfate-reducing bacteria (SRB) contribute to huge economic losses causing microbial-induced corrosion through biofilms on metal surfaces. To effectively combat the challenges posed by SRB, it is essential to understand their molecular mechanisms of biofilm formation. This study aimed to identify ncRNAs in the genome of a model SRB, Oleidesulfovibrio alaskensis G20 (OA G20). Three in silico approaches revealed genome-wide distribution of 37 ncRNAs excluding tRNAs in the OA G20. These ncRNAs belonged to 18 different Rfam families. This study identified riboswitches, sRNAs, RNP, and SRP. The analysis revealed that these ncRNAs could play key roles in the regulation of several pathways of biosynthesis and transport involved in biofilm formation by OA G20. Three sRNAs, Pseudomonas P10, Hammerhead type II, and sX4, which were found in OA G20, are rare and their roles have not been determined in SRB. These results suggest that applying various computational methods could enrich the results and lead to the discovery of additional novel ncRNAs, which could lead to understanding the "rules of life of OA G20" during biofilm formation.

Effects of mRNA conformational switching on translational noise in gene circuits

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

Identification of an N-terminal tag (580N) that improves the biosynthesis of fluorescent proteins in Francisella tularensis and other Gram-negative bacteria

Utilization of fluorescent proteins is widespread for the study of microbial pathogenesis and host-pathogen interactions. Here, we discovered that linkage of the 36 N-terminal amino acids of FTL_0580 (a hypothetical protein of Francisella tularensis) to fluorescent proteins increases the fluorescence emission of bacteria that express these recombinant fusions. This N-terminal peptide will be referred to as 580N. Western blotting revealed that the linkage of 580N to Emerald Green Fluorescent Protein (EmGFP) in F. tularensis markedly improved detection of this protein. We therefore hypothesized that transcripts containing 580N may be translated more efficiently than those lacking the coding sequence for this leader peptide. In support, expression of emGFPFt that had been codon-optimized for F. tularensis, yielded significantly enhanced fluorescence than its non-optimized counterpart. Furthermore, fusing emGFP with coding sequence for a small N-terminal peptide (Serine-Lysine-Isoleucine-Lysine), which had previously been shown to inhibit ribosomal stalling, produced robust fluorescence when expressed in F. tularensis. These findings support the interpretation that 580N enhances the translation efficiency of fluorescent proteins in F. tularensis. Interestingly, expression of non-optimized 580N-emGFP produced greater fluorescence intensity than any other construct. Structural predictions suggested that RNA secondary structure also may be influencing translation efficiency. When expressed in Escherichia coli and Klebsiella pneumoniae bacteria, 580N-emGFP produced increased green fluorescence compared to untagged emGFP (neither allele was codon optimized for these bacteria). In conclusion, fusing the coding sequence for the 580N leader peptide to recombinant genes might serve as an economical alternative to codon optimization for enhancing protein expression in bacteria.

Regulation of bacterial gene expression by non-coding RNA: It is all about time!

Commensal and pathogenic bacteria continuously evolve to survive in diverse ecological niches by efficiently coordinating gene expression levels in their ever-changing environments. Regulation through the RNA transcript itself offers a faster and more cost-effective way to adapt than protein-based mechanisms and can be leveraged for diagnostic or antimicrobial purposes. However, RNA can fold into numerous intricate, not always functional structures that both expand and obscure the plethora of roles that regulatory RNAs serve within the cell. Here, we review the current knowledge of bacterial non-coding RNAs in relation to their folding pathways and interactions. We posit that co-transcriptional folding of these transcripts ultimately dictates their downstream functions. Elucidating the spatiotemporal folding of non-coding RNAs during transcription therefore provides invaluable insights into bacterial pathogeneses and predictive disease diagnostics. Finally, we discuss the implications of co-transcriptional folding andapplications of RNAs for therapeutics and drug targets.

Comparative Study of Single-stranded Oligonucleotides Secondary Structure Prediction Tools

BACKGROUND

Single-stranded nucleic acids (ssNAs) have important biological roles and a high biotechnological potential linked to their ability to bind to numerous molecular targets. This depends on the different spatial conformations they can assume. The first level of ssNAs spatial organisation corresponds to their base pairs pattern, i.e. their secondary structure. Many computational tools have been developed to predict the ssNAs secondary structures, making the choice of the appropriate tool difficult, and an up-to-date guide on the limits and applicability of current secondary structure prediction tools is missing. Therefore, we performed a comparative study of the performances of 9 freely available tools (mfold, RNAfold, CentroidFold, CONTRAfold, MC-Fold, LinearFold, UFold, SPOT-RNA, and MXfold2) on a dataset of 538 ssNAs with known experimental secondary structure.

RESULTS

The minimum free energy-based tools, namely mfold and RNAfold, and some tools based on artificial intelligence, namely CONTRAfold and MXfold2, provided the best results, with [Formula: see text] of exact predictions, whilst MC-fold seemed to be the worst performing tool, with only [Formula: see text] of exact predictions. In addition, UFold and SPOT-RNA are the only options for pseudoknots prediction. Including in the analysis of mfold and RNAfold results 5-10 suboptimal solutions further improved the performances of these tools. Nevertheless, we could observe issues in predicting particular motifs, such as multiple-ways junctions and mini-dumbbells, or the ssNAs whose structure has been determined in complex with a protein. In addition, our benchmark shows that some effort has to be paid for ssDNA secondary structure predictions.

CONCLUSIONS

In general, Mfold, RNAfold, and MXfold2 seem to currently be the best choice for the ssNAs secondary structure prediction, although they still show some limits linked to specific structural motifs. Nevertheless, actual trends suggest that artificial intelligence has a high potential to overcome these remaining issues, for example the recently developed UFold and SPOT-RNA have a high success rate in predicting pseudoknots.

c-di-GMP regulates activity of the PlzA RNA chaperone from the Lyme disease spirochete

PlzA is a c-di-GMP-binding protein crucial for adaptation of the Lyme disease spirochete Borrelia (Borreliella) burgdorferi during its enzootic life cycle. Unliganded apo-PlzA is important for vertebrate infection, while liganded holo-PlzA is important for survival in the tick; however, the biological function of PlzA has remained enigmatic. Here, we report that PlzA has RNA chaperone activity that is inhibited by c-di-GMP binding. Holo- and apo-PlzA bind RNA and accelerate RNA annealing, while only apo-PlzA can strand displace and unwind double-stranded RNA. Guided by the crystal structure of PlzA, we identified several key aromatic amino acids protruding from the N- and C-terminal domains that are required for RNA-binding and unwinding activity. Our findings illuminate c-di-GMP as a switch controlling the RNA chaperone activity of PlzA, and we propose that complex RNA-mediated modulatory mechanisms allow PlzA to regulate gene expression during both the vector and host phases of the B. burgdorferi life cycle.

Riboswitch-mediated regulation of riboflavin biosynthesis genes in prokaryotes

Prokaryotic organisms frequently use riboswitches to quantify intracellular metabolite concentration via high-affinity metabolite receptors. Riboswitches possess a metabolite-sensing system that controls gene regulation in a cis-acting fashion at the initiation of transcriptional/translational level by binding with a specific metabolite and controlling various biochemical pathways. Riboswitch binds with flavin mononucleotide (FMN), a phosphorylated form of riboflavin and controls gene expression involved in riboflavin biosynthesis and transport pathway. The first step of the riboflavin biosynthesis pathway is initiated by the conversion of guanine nucleotide triphosphate (GTP), which is an intermediate of the purine biosynthesis pathway. An alternative pentose phosphate pathway of riboflavin biosynthesis includes the enzymatic conversion of ribulose-5-phosphate into 3, 4 dihydroxy-2-butanone-4-phosphates by DHBP synthase. The product of ribAB interferes with both GTP cyclohydrolase II as well as DHBP synthase activities, which catalyze the cleavage of GTP and converts DHBP Ribu5P in the initial steps of both riboflavin biosynthesis branches. Riboswitches are located in the 5' untranslated region (5' UTR) of messenger RNAs and contain an aptamer domain (highly conserved in sequence) where metabolite binding leads to a conformational change in an aptamer domain, which modulate the regulation of gene expression located on bacterial mRNA. In this review, we focus on how riboswitch regulates the riboflavin biosynthesis pathway in Bacillus subtilis and Lactobacillus plantarum.

Insight rifampicin-resistant (rpoB) mutation in Pseudomonas stutzeri leads to enhance the biosynthesis of secondary metabolites to survive against harsh environments

In this study, a wild-type and five distinct rifampicin-resistant (Rif^r) rpoB mutants of Pseudomonas stutzeri (i.e., Q518R, D521Y, D521V, H531R and I614T) ability were investigated against harsh environments (particularly nutritional complexity). Among these, the robust Rif^r phenotype of P. Stutzeri was associated only with base replacements of the amino deposits. The use of carboxylic and amino acids significantly increased in various Rif^r mutants than that of wild type of P. stutzeri. The assimilation of carbon and nitrogen (N) sources of Rif^r mutants' confirmed that the organism maintains the adaptation in nutritionally complex environments. Acetylene reduction assay at different times also found the variability for N-fixation in all strains. Among them, the highest nitrogenase activity was determined in mutant 'D521V'. The assimilation of carbon and nitrogen sources of P. stutzeri and its Rif^r mutants ensures that the organism maintains the adaptability in nutritionally complex environments through fixing more nitrogen.

Biophysical characterisation of human LincRNA-p21 sense and antisense Alu inverted repeats

Human Long Intergenic Noncoding RNA-p21 (LincRNA-p21) is a regulatory noncoding RNA that plays an important role in promoting apoptosis. LincRNA-p21 is also critical in down-regulating many p53 target genes through its interaction with a p53 repressive complex. The interaction between LincRNA-p21 and the repressive complex is likely dependent on the RNA tertiary structure. Previous studies have determined the two-dimensional secondary structures of the sense and antisense human LincRNA-p21 AluSx1 IRs using SHAPE. However, there were no insights into its three-dimensional structure. Therefore, we in vitro transcribed the sense and antisense regions of LincRNA-p21 AluSx1 Inverted Repeats (IRs) and performed analytical ultracentrifugation, size exclusion chromatography, light scattering, and small angle X-ray scattering (SAXS) studies. Based on these studies, we determined low-resolution, three-dimensional structures of sense and antisense LincRNA-p21. By adapting previously known two-dimensional information, we calculated their sense and antisense high-resolution models and determined that they agree with the low-resolution structures determined using SAXS. Thus, our integrated approach provides insights into the structure of LincRNA-p21 Alu IRs. Our study also offers a viable pipeline for combining the secondary structure information with biophysical and computational studies to obtain high-resolution atomistic models for long noncoding RNAs.

How does RNA fold dynamically?

Recent advances in interrogating RNA folding dynamics have shown the classical model of RNA folding to be incomplete. Here, we pose three prominent questions for the field that are at the forefront of our understanding of the importance of RNA folding dynamics for RNA function. The first centers on the most appropriate biophysical framework to describe changes to the RNA folding energy landscape that a growing RNA chain encounters during transcriptional elongation. The second focuses on the potential ubiquity of strand displacement - a process by which RNA can rapidly change conformations - and how this process may be generally present in broad classes of seemingly different RNAs. The third raises questions about the potential importance and roles of cellular protein factors in RNA conformational switching. Answers to these questions will greatly improve our fundamental knowledge of RNA folding and function, drive biotechnological advances that utilize engineered RNAs, and potentially point to new areas of biology yet to be discovered.

Probing RNA structures and functions by solvent accessibility: an overview from experimental and computational perspectives

Characterizing RNA structures and functions have mostly been focused on 2D, secondary and 3D, tertiary structures. Recent advances in experimental and computational techniques for probing or predicting RNA solvent accessibility make this 1D representation of tertiary structures an increasingly attractive feature to explore. Here, we provide a survey of these recent developments, which indicate the emergence of solvent accessibility as a simple 1D property, adding to secondary and tertiary structures for investigating complex structure–function relations of RNAs.

Influence of Flagellin Polymorphisms, Gene Regulation, and Responsive Memory on the Motility of Xanthomonas Species That Cause Bacterial Spot Disease of Solanaceous Plants

Increasingly, new evidence has demonstrated variability in the epitope regions of bacterial flagellin, including in regions harboring the microbe-associated molecular patterns flg22 and flgII-28 that are recognized by the pattern recognition receptors FLS2 and FLS3, respectively. Additionally, because bacterial motility is known to contribute to pathogen virulence and chemotaxis, reductions in or loss of motility can significantly reduce bacterial fitness. In this study, we determined that variations in flg22 and flgII-28 epitopes allow some but not all Xanthomonas spp. to evade both FLS2- and FLS3-mediated oxidative burst responses. We observed variation in the motility for many isolates, regardless of their flagellin sequence. Instead, we determined that past growth conditions may have a significant impact on the motility status of isolates, because we could minimize this variability by inducing motility using chemoattractant assays. Additionally, motility could be significantly suppressed under nutrient-limited conditions, and bacteria could "remember" its prior motility status after storage at ultracold temperatures. Finally, we observed larger bacterial populations of strains with flagellin variants predicted not to be recognized by either FLS2 or FLS3, suggesting that these bacteria can evade flagellin recognition in tomato plants. Although some flagellin variants may impart altered motility and differential recognition by the host immune system, external growth parameters and gene expression regulation appear to have more significant impacts on the motility phenotypes for these Xanthomonas spp.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.