Widespread Influence of 3'-End Structures on Mammalian mRNA Processing and Stability.

The physiological relevance of structures within mammalian mRNAs has been elusive, as these mRNAs are less folded in cells than in vitro and have predicted secondary structures no more stable than those of random sequences. Here, we investigate the possibility that mRNA structures facilitate the 3'-end processing of thousands of human mRNAs by juxtaposing poly(A) signals (PASs) and cleavage sites that are otherwise too far apart. We find that RNA structures are predicted to be more prevalent within these extended 3'-end regions than within PAS-upstream regions and indeed are substantially more folded within cells, as determined by intracellular probing. Analyses of thousands of ectopically expressed variants demonstrate that this folding both enhances processing and increases mRNA metabolic stability. Even folds with predicted stabilities resembling those of random sequences can enhance processing. Structure-controlled processing can also regulate neighboring gene expression. Thus, RNA structure has widespread roles in mammalian mRNA biogenesis and metabolism.

Magnesium ions mitigate metastable states in the regulatory landscape of mRNA elements.

Residing in the 5' untranslated region of the mRNA, the 2'-deoxyguanosine (2'-dG) riboswitch mRNA element adopts an alternative structure upon binding of the 2'-dG molecule, which terminates transcription. RNA conformations are generally strongly affected by positively charged metal ions (especially Mg2+). We have quantitatively explored the combined effect of ligand (2'-dG) and Mg2+ binding on the energy landscape of the aptamer domain of the 2'-dG riboswitch with both explicit solvent all-atom molecular dynamics simulations (99 µsec aggregate sampling for the study) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments. We show that both ligand and Mg2+ are required for the stabilization of the aptamer domain; however, the two factors act with different modalities. The addition of Mg2+ remodels the energy landscape and reduces its frustration by the formation of additional contacts. In contrast, the binding of 2'-dG eliminates the metastable states by nucleating a compact core for the aptamer domain. Mg2+ ions and ligand binding are required to stabilize the least stable helix, P1 (which needs to unfold to activate the transcription platform), and the riboswitch core formed by the backbone of the P2 and P3 helices. Mg2+ and ligand also facilitate a more compact structure in the three-way junction region.

mRNA Fragmentation Pattern Detected by SHAPE.

The success of messenger RNA (mRNA) vaccines in controlling COVID-19 has warranted further developments in new technology. Currently, their quality control process largely relies on low-resolution electrophoresis for detecting chain breaks. Here, we present an approach using multi-primer reverse transcription sequencing (MPRT-seq) to identify degradation fragments in mRNA products. Using this in-house-made mRNA containing two antigens and untranslated regions (UTRs), we analyzed the mRNA completeness and degradation pattern at a nucleotide resolution. We then analyzed the sensitive base sequence and its correlation with the secondary structure. Our MPRT-seq mapping shows that certain sequences on the 5' of bulge-stem-loop structures can result in preferential chain breaks. Our results agree with commonly used capillary electrophoresis (CE) integrity analysis but at a much higher resolution, and can improve mRNA stability by providing information to remove sensitive structures or sequences in the mRNA sequence design.

Unraveling the Role of BRCA1 variants in Dysregulation of Transcriptional and Post-Transcriptional Mechanisms in Breast Cancer.

Objective: To screen BRCA1 gene variants and predict potential role of the identified variants in breast cancer. Method: This case-control study included two hundred and fifty breast cancer patients and equal healthy individuals from the Federal Breast Cancer Screening Centre, Pakistan Institute of Medical Sciences, Islamabad from March 2021- January 2023. Demographic data was collected through questionnaires and clinical data was assessed using mammograms, ultrasound, histopathology and immunohistochemistry reports. Polymerase chain reaction and Sanger sequencing approach were used to detect variants in BRCA1 gene. In-silico analyses were carried out to predict mutation effect, miRNA binding site alterations and change in mRNA structure and stability. Results: Invasive ductal carcinoma was the most prevalent type of breast cancer. Old age [OR: 2.8149 (1.5995 to 4.9538) p value = 0.0003] and family history [OR: 4.3186 (1.7336 to 10.7581) p value = 0.001] were significant breast cancer risk. Six variants were identified. Two novel missense variants, Chr17:43082553A>T and Chr17:43093710A>T were predicted deleterious as these disrupted interaction with PALB2 and importin alpha's NLS2 site, respectively. In silico analysis predicted the loss of hsa-miR-1179 binding site due to variant Chr17:43093220T>C. Moreover, four variants were predicted to affect the mRNA structure and stability. Conclusion: Two novel variants were predicted to be pathogenic. In-silico analysis predicted the loss of miRNA binding site along with change in mRNA secondary structure plus stability, possible mechanisms for oncogenesis. Further, expressional studies are required to confirm BRCA1 gene dysregulation in breast cancer due to these variants.

Initiator AUGs Are Discriminated from Elongator AUGs Predominantly through mRNA Accessibility in C. crescentus.

The initiation of translation in bacteria is thought to occur upon base pairing between the Shine-Dalgarno (SD) site in the mRNA and the anti-SD site in the rRNA. However, in many bacterial species, such as Caulobacter crescentus, a minority of mRNAs have SD sites. To examine the functional importance of SD sites in C. crescentus, we analyzed the transcriptome and found that more SD sites exist in the coding sequence than in the preceding start codons. To examine the function of SD sites in initiation, we designed a series of mutants with altered ribosome accessibility and SD content in translation initiation regions (TIRs) and in elongator AUG regions (EARs). A lack of mRNA structure content is required for initiation in TIRs, and, when introduced into EARs, can stimulate initiation, thereby suggesting that low mRNA structure content is a major feature that is required for initiation. SD sites appear to stimulate initiation in TIRs, which generally lack structure content, but SD sites only stimulate initiation in EARs if RNA secondary structures are destabilized. Taken together, these results suggest that the difference in secondary structure between TIRs and EARs directs ribosomes to start codons where SD base pairing can tune the efficiency of initiation, but SDs in EARs do not stimulate initiation, as they are blocked by stable secondary structures. This highlights the importance of studying translation initiation mechanisms in diverse bacterial species. IMPORTANCE Start codon selection is an essential process that is thought to occur via the base pairing of the rRNA to the SD site in the mRNA. This model is based on studies in E. coli, yet whole-genome sequencing revealed that SD sites are absent at start codons in many species. By examining the transcriptome of C. crescentus, we found more SD-AUG pairs in the CDS of mRNAs than preceding start codons, yet these internal sites do not initiate. Instead, start codon regions have lower mRNA secondary structure content than do internal SD-AUG regions. Therefore, we find that start codon selection is not controlled by the presence of SD sites, which tune initiation efficiency, but by lower mRNA structure content surrounding the start codon.

Fully sequencing the cassava full-length cDNA library reveals unannotated transcript structures and alternative splicing events in regions with a high density of single nucleotide variations, insertions-deletions, and heterozygous sequences.

The primary transcript structure provides critical insights into protein diversity, transcriptional modification, and functions. Cassava transcript structures are highly diverse because of alternative splicing (AS) events and high heterozygosity. To precisely determine and characterize transcript structures, fully sequencing cloned transcripts is the most reliable method. However, cassava annotations were mainly determined according to fragmentation-based sequencing analyses (e.g., EST and short-read RNA-seq). In this study, we sequenced the cassava full-length cDNA library, which included rare transcripts. We obtained 8,628 non-redundant fully sequenced transcripts and detected 615 unannotated AS events and 421 unannotated loci. The different protein sequences resulting from the unannotated AS events tended to have diverse functional domains, implying that unannotated AS contributes to the truncation of functional domains. The unannotated loci tended to be derived from orphan genes, implying that the loci may be associated with cassava-specific traits. Unexpectedly, individual cassava transcripts were more likely to have multiple AS events than Arabidopsis transcripts, suggestive of the regulated interactions between cassava splicing-related complexes. We also observed that the unannotated loci and/or AS events were commonly in regions with abundant single nucleotide variations, insertions-deletions, and heterozygous sequences. These findings reflect the utility of completely sequenced FLcDNA clones for overcoming cassava-specific annotation-related problems to elucidate transcript structures. Our work provides researchers with transcript structural details that are useful for annotating highly diverse and unique transcripts and alternative splicing events.

A Comprehensive Review of mRNA Vaccines.

mRNA vaccines have been demonstrated as a powerful alternative to traditional conventional vaccines because of their high potency, safety and efficacy, capacity for rapid clinical development, and potential for rapid, low-cost manufacturing. These vaccines have progressed from being a mere curiosity to emerging as COVID-19 pandemic vaccine front-runners. The advancements in the field of nanotechnology for developing delivery vehicles for mRNA vaccines are highly significant. In this review we have summarized each and every aspect of the mRNA vaccine. The article describes the mRNA structure, its pharmacological function of immunity induction, lipid nanoparticles (LNPs), and the upstream, downstream, and formulation process of mRNA vaccine manufacturing. Additionally, mRNA vaccines in clinical trials are also described. A deep dive into the future perspectives of mRNA vaccines, such as its freeze-drying, delivery systems, and LNPs targeting antigen-presenting cells and dendritic cells, are also summarized.

mRNAs sequestered in stress granules recover nearly completely for translation.

Stress granules (SGs) are membrane-less condensates composed of RNA and protein that assemble in response to stress stimuli and disassemble when stress is lifted. Both assembly and disassembly are tightly controlled processes, yet, it remains elusive whether mRNAs in SGs completely recover for translation following stress relief. Using RNA-seq of translating fractions in human cell line, we found that higher fraction of the m6A-modified mRNAs recovered for translation compared to unmodified mRNAs, i.e. 95% vs 84%, respectively. Considering structural mRNA analysis, we found that the m6A modification enhances structuring at nucleotides in its close vicinity. Our results suggest that SG-sequestered mRNAs disassemble nearly completely from SGs and the m6A modification may display some advantage to the mRNAs in their recovery for translation likely by m6A-driven structural stabilization.

Decoding on the ribosome depends on the structure of the mRNA phosphodiester backbone.

During translation, the ribosome plays an active role in ensuring that mRNA is decoded accurately and rapidly. Recently, biochemical studies have also implicated certain accessory factors in maintaining decoding accuracy. However, it is currently unclear whether the mRNA itself plays an active role in the process beyond its ability to base pair with the tRNA. Structural studies revealed that the mRNA kinks at the interface of the P and A sites. A magnesium ion appears to stabilize this structure through electrostatic interactions with the phosphodiester backbone of the mRNA. Here we examined the role of the kink structure on decoding using a well-defined in vitro translation system. Disruption of the kink structure through site-specific phosphorothioate modification resulted in an acute hyperaccurate phenotype. We measured rates of peptidyl transfer for near-cognate tRNAs that were severely diminished and in some instances were almost 100-fold slower than unmodified mRNAs. In contrast to peptidyl transfer, the modifications had little effect on GTP hydrolysis by elongation factor thermal unstable (EF-Tu), suggesting that only the proofreading phase of tRNA selection depends critically on the kink structure. Although the modifications appear to have no effect on typical cognate interactions, peptidyl transfer for a tRNA that uses atypical base pairing is compromised. These observations suggest that the kink structure is important for decoding in the absence of Watson-Crick or G-U wobble base pairing at the third position. Our findings provide evidence for a previously unappreciated role for the mRNA backbone in ensuring uniform decoding of the genetic code.

Deciphering the rules of mRNA structure differentiation in Saccharomyces cerevisiae in vivo and in vitro with deep neural networks.

The structure of mRNA in vivo is unwound to some extent in response to multiple factors involved in the translation process, resulting in significant differences from the structure of the same mRNA in vitro. In this study, we have proposed a novel application of deep neural networks, named DeepDRU, to predict the degree of mRNA structure unwinding in vivo by fitting five quantifiable features that may affect mRNA folding: ribosome density (RD), minimum folding free energy (MFE), GC content, translation initiation ribosome density (INI) and mRNA structure position (POS). mRNA structures with adjustment of the simulated structural features were designed and then fed into the trained DeepDRU model. We found unique effect regions of these five features on mRNA structure in vivo. Strikingly, INI is the most critical factor affecting the structure of mRNA in vivo, and structural sequence features, including MFE and GC content, have relatively smaller effects. DeepDRU provides a new paradigm for predicting the unwinding capability of mRNA structure in vivo. This improved knowledge about the mechanisms of factors influencing the structural capability of mRNA to unwind will facilitate the design and functional analysis of mRNA structure in vivo.

Exposing synonymous mutations.

Synonymous codon changes, which do not alter protein sequence, were previously thought to have no functional consequence. Although this concept has been overturned in recent years, there is no unique mechanism by which these changes exert biological effects. A large repertoire of both experimental and bioinformatic methods has been developed to understand the effects of synonymous variants. Results from this body of work have provided global insights into how biological systems exploit the degeneracy of the genetic code to control gene expression, protein folding efficiency, and the coordinated expression of functionally related gene families. Although it is now clear that synonymous variants are important in a variety of contexts, from human disease to the safety and efficacy of therapeutic proteins, there is no clear consensus on the approaches to identify and validate these changes. Here, we review the diverse methods to understand the effects of synonymous mutations.

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.

Computational Methods for Modeling Aptamers and Designing Riboswitches.

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Modulation of protein synthesis by polyamines.

Polyamines are ubiquitous small basic molecules that play important roles in cell growth and viability. Since polyamines mainly exist as a polyamine-RNA complex, we looked for proteins whose synthesis is preferentially stimulated by polyamines at the level of translation, and thus far identified 17 proteins in Escherichia coli and 6 proteins in eukaryotes. The mechanisms of polyamine stimulation of synthesis of these proteins were investigated. In addition, the role of eIF5A, containing hypusine formed from spermidine, on protein synthesis is described. These results clearly indicate that polyamines and eIF5A contribute to cell growth and viability through modulation of protein synthesis.

Causal signals between codon bias, mRNA structure, and the efficiency of translation and elongation.

Ribosome profiling data report on the distribution of translating ribosomes, at steady-state, with codon-level resolution. We present a robust method to extract codon translation rates and protein synthesis rates from these data, and identify causal features associated with elongation and translation efficiency in physiological conditions in yeast. We show that neither elongation rate nor translational efficiency is improved by experimental manipulation of the abundance or body sequence of the rare AGG tRNA. Deletion of three of the four copies of the heavily used ACA tRNA shows a modest efficiency decrease that could be explained by other rate-reducing signals at gene start. This suggests that correlation between codon bias and efficiency arises as selection for codons to utilize translation machinery efficiently in highly translated genes. We also show a correlation between efficiency and RNA structure calculated both computationally and from recent structure probing data, as well as the Kozak initiation motif, which may comprise a mechanism to regulate initiation.

Translational halt during elongation caused by G-quadruplex formed by mRNA.

mRNA forms various secondary and tertiary structures that affect gene expression. Although structures formed in the untranslated regions (UTRs) of mRNAs that inhibit translation have been characterized, stable mRNA structures in open reading frames (ORFs) may also cause translational halt or slow translation elongation. We previously established a method, termed a synchronized translation assay, that enables time course analysis of single turnover translation elongation. In this method, translation initiation, which is a rate determining step of the translation procedure, can be ignored because all ribosomes are synchronized on a specific position of mRNA before translation elongation is restarted from this position. In this paper, we used the synchronized translation assay to evaluate the effects of a G-quadruplex structure located at various positions within the mRNA ORF on translational halt.

Computational design of mRNA vaccines.

mRNA technology has emerged as a successful vaccine platform that offered a swift response to the COVID-19 pandemic. Accumulating evidence shows that vaccine efficacy, thermostability, and other important properties, are largely impacted by intrinsic properties of the mRNA molecule, such as RNA sequence and structure, both of which can be optimized. Designing mRNA sequence for vaccines presents a combinatorial problem due to an extremely large selection space. For instance, due to the degeneracy of the genetic code, there are over 10632 possible mRNA sequences that could encode the spike protein, the COVID-19 vaccines' target. Moreover, designing different elements of the mRNA sequence simultaneously against multiple objectives such as translational efficiency, reduced reactogenicity, and improved stability requires an efficient and sophisticated optimization strategy. Recently, there has been a growing interest in utilizing computational tools to redesign mRNA sequences to improve vaccine characteristics and expedite discovery timelines. In this review, we explore important biophysical features of mRNA to be considered for vaccine design and discuss how computational approaches can be applied to rapidly design mRNA sequences with desirable characteristics.

Strategies to overcome the challenges of low or no expression of heterologous proteins in Escherichia coli.

Protein expression is a critical process in diverse biological systems. For Escherichia coli, a widely employed microbial host in industrial catalysis and healthcare, researchers often face significant challenges in constructing recombinant expression systems. To maximize the potential of E. coli expression systems, it is essential to address problems regarding the low or absent production of certain target proteins. This article presents viable solutions to the main factors posing challenges to heterologous protein expression in E. coli, which includes protein toxicity, the intrinsic influence of gene sequences, and mRNA structure. These strategies include specialized approaches for managing toxic protein expression, addressing issues related to mRNA structure and codon bias, advanced codon optimization methodologies that consider multiple factors, and emerging optimization techniques facilitated by big data and machine learning.

Advancements of in vitro transcribed mRNA (IVT mRNA) to enable translation into the clinics.

The accelerated progress and approval of two mRNA-based vaccines to address the SARS-CoV-2 virus were unprecedented. This record-setting feat was made possible through the solid foundation of research on in vitro transcribed mRNA (IVT mRNA) which could be utilized as a therapeutic modality. Through decades of thorough research to overcome barriers to implementation, mRNA-based vaccines or therapeutics offer many advantages to rapidly address a broad range of applications including infectious diseases, cancers, and gene editing. Here, we describe the advances that have supported the adoption of IVT mRNA in the clinics, including optimization of the IVT mRNA structural components, synthesis, and lastly concluding with different classes of IVT RNA. Continuing interest in driving IVT mRNA technology will enable a safer and more efficacious therapeutic modality to address emerging and existing diseases.

Molecular cloning, expression, mRNA secondary structure and immunological characterization of mussel foot proteins (Mfps) (Mollusca: Bivalvia).

The macroscale production of mussel foot proteins (Mfps) in the expression system has not succeeded to date. The principal reasons for this are low levels of expression and yield of Mfps, lack of post-translational modifications (PTMs), and immunological toxic effects on the host system. Identification of post-translational modification sites, suitable expression hosts, and immunological responses through an experimental approach is very costly and time-consuming. However, in the present study, in silico post-translation modification, antigenicity, allergenicity, and the immunological reaction of all available Mfps were characterized. Furthermore, all Mfps were codon optimized in three different expression systems to determine the best expression host. Finally, we performed the in-silico cloning of all codon-optimized Mfps in a suitable host (E. coli K12, pET28a(+) vector) and analyzed the secondary structure of mRNA and its structural stability. Among the 78 Mfps, six fps are considered potential allergenic proteins, six fps are considered non-allergenic proteins, and all other fps are probably allergenic. High antigenicity was observed in bacterial cells as compared to yeast and tumor cells. Nevertheless, the predicted expression of Mfps in a bacterial host is higher than in other expression hosts. Important to note that all Mfps showed significant immunological activity in the human system, and we concluded that these antigenic, allergenic, and immunological properties are directly correlated with their amino acid composition. The study's major goal is to provide a comprehensive understanding of Mfps and aid in the future genetic engineering and expression of Mfps and its diverse applications in different fields.Communicated by Ramaswamy H. Sarma.