Sequence-structure relationships in yeast mRNAs

Probing RNA structural landscapes across Candida yeast genomes

Understanding the intricate roles of RNA molecules in virulence and host-pathogen interactions can provide valuable insights into combatting infections and improving human health. Although much progress has been achieved in understanding transcriptional regulation during host-pathogen interactions in diverse species, more is needed to know about the structure of pathogen RNAs. This is particularly true for fungal pathogens, including pathogenic yeasts of the Candida genus, which are the leading cause of hospital-acquired fungal infections. Our work addresses the gap between RNA structure and their biology by employing genome-wide structure probing to comprehensively explore the structural landscape of mRNAs and long non-coding RNAs (lncRNAs) in the four major Candida pathogens. Specifically focusing on mRNA, we observe a robust correlation between sequence conservation and structural characteristics in orthologous transcripts, significantly when sequence identity exceeds 50%, highlighting structural feature conservation among closely related species. We investigate the impact of single nucleotide polymorphisms (SNPs) on mRNA secondary structure. SNPs within 5' untranslated regions (UTRs) tend to occur in less structured positions, suggesting structural constraints influencing transcript regulation. Furthermore, we compare the structural properties of coding regions and UTRs, noting that coding regions are generally more structured than UTRs, consistent with similar trends in other species. Additionally, we provide the first experimental characterization of lncRNA structures in Candida species. Most lncRNAs form independent subdomains, similar to human lncRNAs. Notably, we identify hairpin-like structures in lncRNAs, a feature known to be functionally significant. Comparing hairpin prevalence between lncRNAs and protein-coding genes, we find enrichment in lncRNAs across Candida species, humans, and Arabidopsis thaliana, suggesting a conserved role for these structures. In summary, our study offers valuable insights into the interplay between RNA sequence, structure, and function in Candida pathogens, with implications for gene expression regulation and potential therapeutic strategies against Candida infections.

Modeling and computational characterization of a Xanthomonas sp. Hypothetical protein identifies a remote ortholog of Burkholderia lethal factor 1

Burkholderia Lethal Factor 1 (BLF1) is a deamidase first characterized in Burkholderia pseudomallei. This enzyme inhibits cellular protein synthesis by deamidating a glutamine residue to a glutamic acid in its target protein, the eukaryotic translation initiation factor 4 A (eIF4A). In this work, we present the characterization of a hypothetical protein from Xanthomonas sp. Leaf131 as the first report of a BLF1 family ortholog outside of the Burkholderia genus. Although standard sequence similarity searches such as BLAST were not able to detect the homology between the Xanthomonas sp. Leaf131 hypothetical protein sequence and BLF1, our computed structure model for the Xanthomonas sp. hypothetical protein revealed structural similarities with an RMSD of 2.7 Å/164 Cα atoms and a TM-score of 0.72 when superposed. Structural comparisons of the Xanthomonas model structure against BLF1 and Escherichia coli cytotoxic necrotizing factor 1 (CNF1) revealed that the conserved signature LXGC motif and putative catalytic residues are structurally aligned thus signifying a level of functional or mechanistic similarity. Protein-protein docking analysis and molecular dynamics simulations also demonstrated that eIF4A could still be a possible target substrate for deamidation by XLF1 as it is for BLF1. We therefore propose that this Xanthomonas hypothetical protein be renamed as Xanthomonas Lethal Factor 1 (XLF1). Our work also provides further evidence of the utility of programs such as AlphaFold in bridging the computational function annotation transfer gap despite very low sequence identities of under 20%.Communicated by Ramaswamy H. Sarma.

Conserved RNA structures in the intergenic regions of ambisense viruses

Ambisense viruses are negative-sense single-stranded RNA viruses that use a unique expression strategy. Their genome contains at least one ambisense RNA segment that carries two oppositely oriented reading frames separated by an intergenic region. It is believed that a structural RNA element within the intergenic region is involved in transcription termination. However, a general overview over the structural repertoire of ambisense intergenic regions is currently lacking. In this study we investigated the structural potential of the intergenic regions of all known ambisense viruses and compared their structural repertoire by structure-guided clustering. Intergenic regions of most ambisense viruses possess a high potential to build stable secondary structures and many viruses share common structural motifs in the intergenic regions of their ambisense segments. We demonstrate that (i) within the phylogenetic virus groups sets of conserved functional structures are present, but that (ii) between the groups conservation is low to non-existent. These results reflect a high degree of freedom to regulate ambisense transcription termination and also imply that the genetic strategy of having an ambisense RNA genome has evolved several times independently.

Melting temperature highlights functionally important RNA structure and sequence elements in yeast mRNA coding regions

Secondary structure elements in the coding regions of mRNAs play an important role in gene expression and regulation, but distinguishing functional from non-functional structures remains challenging. Here we investigate the dependence of sequence–structure relationships in the coding regions on temperature based on the recent PARTE data by Wan et al. Our main finding is that the regions with high and low thermostability (high T_m and low T_m regions) are under evolutionary pressure to preserve RNA secondary structure and primary sequence, respectively. Sequences of low T_m regions display a higher degree of evolutionary conservation compared to high T_m regions. Low T_m regions are under strong synonymous constraint, while high T_m regions are not. These findings imply that high T_m regions contain thermo-stable functionally important RNA structures, which impose relaxed evolutionary constraint on sequence as long as the base-pairing patterns remain intact. By contrast, low thermostability regions contain single-stranded functionally important conserved RNA sequence elements accessible for binding by other molecules. We also find that theoretically predicted structures of paralogous mRNA pairs become more similar with growing temperature, while experimentally measured structures tend to diverge, which implies that the melting pathways of RNA structures cannot be fully captured by current computational approaches.

Adaptation of mRNA structure to control protein folding

Comparison of mRNA and protein structures shows that highly structured mRNAs typically encode compact protein domains suggesting that mRNA structure controls protein folding. This function is apparently performed by distinct structural elements in the mRNA, which implies 'fine tuning' of mRNA structure under selection for optimal protein folding. We find that, during evolution, changes in the mRNA folding energy follow amino acid replacements, reinforcing the notion of an intimate connection between the structures of a mRNA and the protein it encodes, and the double encoding of protein sequence and folding in the mRNA.

Protein functional features are reflected in the patterns of mRNA translation speed

BACKGROUND

The degeneracy of the genetic code makes it possible for the same amino acid string to be coded by different messenger RNA (mRNA) sequences. These "synonymous mRNAs" may differ largely in a number of aspects related to their overall translational efficiency, such as secondary structure content and availability of the encoded transfer RNAs (tRNAs). Consequently, they may render different yields of the translated polypeptides. These mRNA features related to translation efficiency are also playing a role locally, resulting in a non-uniform translation speed along the mRNA, which has been previously related to some protein structural features and also used to explain some dramatic effects of "silent" single-nucleotide-polymorphisms (SNPs). In this work we perform the first large scale analysis of the relationship between three experimental proxies of mRNA local translation efficiency and the local features of the corresponding encoded proteins.

RESULTS

We found that a number of protein functional and structural features are reflected in the patterns of ribosome occupancy, secondary structure and tRNA availability along the mRNA. One or more of these proxies of translation speed have distinctive patterns around the mRNA regions coding for certain protein local features. In some cases the three patterns follow a similar trend. We also show specific examples where these patterns of translation speed point to the protein's important structural and functional features.

CONCLUSIONS

This support the idea that the genome not only codes the protein functional features as sequences of amino acids, but also as subtle patterns of mRNA properties which, probably through local effects on the translation speed, have some consequence on the final polypeptide. These results open the possibility of predicting a protein's functional regions based on a single genomic sequence, and have implications for heterologous protein expression and fine-tuning protein function.

Evolutionary Origin and Conserved Structural Building Blocks of Riboswitches and Ribosomal RNAs: Riboswitches as Probable Target Sites for Aminoglycosides Interaction

PURPOSE

Riboswitches, as noncoding RNA sequences, control gene expression through direct ligand binding. Sporadic reports on the structural relation of riboswitches with ribosomal RNAs (rRNA), raises an interest in possible similarity between riboswitches and rRNAs evolutionary origins. Since aminoglycoside antibiotics affect microbial cells through binding to functional sites of the bacterial rRNA, finding any conformational and functional relation between riboswitches/rRNAs is utmost important in both of medicinal and basic research.

METHODS

Analysis of the riboswitches structures were carried out using bioinformatics and computational tools. The possible functional similarity of riboswitches with rRNAs was evaluated based on the affinity of paromomycin antibiotic (targeting "A site" of 16S rRNA) to riboswitches via docking method.

RESULTS

There was high structural similarity between riboswitches and rRNAs, but not any particular sequence based similarity between them was found. The building blocks including "hairpin loop containing UUU", "peptidyl transferase center conserved hairpin A loop"," helix 45" and "S2 (G8) hairpin" as high identical rRNA motifs were detected in all kinds of riboswitches. Surprisingly, binding energies of paromomycin with different riboswitches are considerably better than the binding energy of paromomycin with "16S rRNA A site". Therefore the high affinity of paromomycin to bind riboswitches in comparison with rRNA "A site" suggests a new insight about riboswitches as possible targets for aminoglycoside antibiotics.

CONCLUSION

These findings are considered as a possible supporting evidence for evolutionary origin of riboswitches/rRNAs and also their role in the exertion of antibiotics effects to design new drugs based on the concomitant effects via rRNA/riboswitches.

Solving nucleic acid structures by molecular replacement: examples from group II intron studies

Structured RNA molecules are key players in ensuring cellular viability. It is now emerging that, like proteins, the functions of many nucleic acids are dictated by their tertiary folds. At the same time, the number of known crystal structures of nucleic acids is also increasing rapidly. In this context, molecular replacement will become an increasingly useful technique for phasing nucleic acid crystallographic data in the near future. Here, strategies to select, create and refine molecular-replacement search models for nucleic acids are discussed. Using examples taken primarily from research on group II introns, it is shown that nucleic acids are amenable to different and potentially more flexible and sophisticated molecular-replacement searches than proteins. These observations specifically aim to encourage future crystallographic studies on the newly discovered repertoire of noncoding transcripts.

Specific temperature-induced perturbations of secondary mRNA structures are associated with the cold-adapted temperature-sensitive phenotype of influenza A virus

For decades, cold-adapted, temperature-sensitive (ca/ts) strains of influenza A virus have been used as live attenuated vaccines. Due to their great public health importance it is crucial to understand the molecular mechanism(s) of cold adaptation and temperature sensitivity that are currently unknown. For instance, secondary RNA structures play important roles in influenza biology. Thus, we hypothesized that a relatively minor change in temperature (32-39°C) can lead to perturbations in influenza RNA structures and, that these structural perturbations may be different for mRNAs of the wild type (wt) and ca/ts strains. To test this hypothesis, we developed a novel in silico method that enables assessing whether two related RNA molecules would undergo (dis)similar structural perturbations upon temperature change. The proposed method allows identifying those areas within an RNA chain where dissimilarities of RNA secondary structures at two different temperatures are particularly pronounced, without knowing particular RNA shapes at either temperature. We identified such areas in the NS2, PA, PB2 and NP mRNAs. However, these areas are not identical for the wt and ca/ts mutants. Differences in temperature-induced structural changes of wt and ca/ts mRNA structures may constitute a yet unappreciated molecular mechanism of the cold adaptation/temperature sensitivity phenomena.

Flexible programming of cell-free protein synthesis using magnetic bead-immobilized plasmids

The use of magnetic bead-immobilized DNA as movable template for cell-free protein synthesis has been investigated. Magnetic microbeads containing chemically conjugated plasmids were used to direct cell-free protein synthesis, so that protein generation could be readily programmed, reset and reprogrammed. Protein synthesis by using this approach could be ON/OFF-controlled through repeated addition and removal of the microbead-conjugated DNA and employed in sequential expression of different genes in a same reaction mixture. Since the incubation periods of individual template plasmids are freely controllable, relative expression levels of multiple proteins can be tuned to desired levels. We expect that the presented results will find wide application to the flexible design and execution of synthetic pathways in cell-free chassis.