Structure Of Messenger RNA Research Articles

Novel eIF4A1 inhibitors with anti‐tumor activity in lymphoma

BackgroundDeregulated translation initiation is implicated extensively in cancer initiation and progression. It is actively pursued as a viable target that circumvents the dependency on oncogenic signaling, a significant factor in current strategies. Eukaryotic translation initiation factor (eIF) 4A plays an essential role in translation initiation by unwinding the secondary structure of messenger RNA (mRNA) upstream of the start codon, enabling active ribosomal recruitment on the downstream genes. Several natural product molecules with similar scaffolds, such as Rocaglamide A (RocA), targeting eIF4A have been reported in the last decade. However, their clinical utilization is still elusive due to several pharmacological limitations. In this study we identified new eIF4A1 inhibitors and their possible mechanisms.MethodsIn this report, we conducted a pharmacophore-based virtual screen of RocA complexed with eIF4A and a polypurine RNA strand for novel eIF4A inhibitors from commercially available compounds in the MolPort Database. We performed target-based screening and optimization of active pharmacophores. We assessed the effects of novel compounds on biochemical and cell-based assays for efficacy and mechanistic evaluation.ResultsWe validated three new potent eIF4A inhibitors, RBF197, RBF 203, and RBF 208, which decreased diffuse large B-cell lymphoma (DLBCL) cell viability. Biochemical and cellular studies, molecular docking, and functional assays revealed that thosenovel compounds clamp eIF4A into mRNA in an ATP-independent manner. Moreover, we found that RBF197 and RBF208 significantly depressed eIF4A-dependent oncogene expression as well as the colony formation capacity of DLBCL. Interestingly, exposure of these compounds to non-malignant cells had only minimal impact on their growth and viability.ConclusionsIdentified compounds suggest a new strategy for designing novel eIF4A inhibitors.

Cold Shock Response in Bacteria

Bacteria often encounter temperature fluctuations in their natural habitats and must adapt to survive. The molecular response of bacteria to sudden temperature upshift or downshift is termed the heat shock response (HSR) or the cold shock response (CSR), respectively. Unlike the HSR, which activates a dedicated transcription factor that predominantly copes with heat-induced protein folding stress, the CSR is mediated by a diverse set of inputs. This review provides a picture of our current understanding of the CSR across bacteria. The fundamental aspects of CSR involved in sensing and adapting to temperature drop, including regulation of membrane fluidity, protein folding, DNA topology, RNA metabolism, and protein translation, are discussed. Special emphasis is placed on recent findings of a CSR circuitry in Escherichia coli mediated by cold shock family proteins and RNase R that monitors and modulates messenger RNA structure to facilitate global translation recovery during acclimation.

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population.

BackgroundThe role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of “silent” genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking.ResultsWe developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population.ConclusionsThese findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a “Structural Predictivity Index” (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

Encapsulation state of messenger RNA inside lipid nanoparticles

Understanding the structure of messenger RNA (mRNA) lipid nanoparticles, and specifically the microenvironment of the mRNA molecules within these entities, is fundamental to advancing their biomedical potential. Here, we show that a permeating cationic dye, thionine, can serve as a cryogenic electron microscopy contrasting agent by binding selectively to encapsulated mRNA without disturbing lipid nanoparticle morphology. Cryo-electron microscopy images identify the mRNA location, revealing that mRNA may exist within solvent-filled cavities or may be substantially lipid associated.

Advanced sequence optimization for the high efficient yield of human group A rotavirus VP6 recombinant protein in Escherichia coli and its use as immunogen.

Rotavirus is the important etiological agents of infectious diarrhea among children under 5 years old. Rotaviruses are divided into 10 serogroups (A-J) and each group is based on genetic properties of major structural protein VP6. We designed a novel VP6 sequence optimization to increase the expression level of this protein. Numerous factors such as codon adaptation index, codon pair bias, and guanine-cytosine content were adapted based on Escherichiacoli codon usage. In addition, the ribosome binding site (RBS) of pET-15b was redesigned by the RBS calculator and the secondary structure of VP6 messenger RNA was optimized in the whole length of the coding sequence. Various factors including isopropyl beta- d-thiogalactoside (IPTG) concentration, temperature, and induction time were analyzed for the optimization of the best expression in E. coli by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and western blotting. The recombinant VP6 (rVP6) protein was purified by the Ni-sepharose and then the hyperimmune sera were generated against rVP6 in rabbits. Among three different temperatures, IPTG concentrations, and postinductions, the level of rVP6 was higher at 37°C, 1 mM of IPTG, and 8 h, respectively. Also, the high expression level of rVP6 was obtained in the insoluble aggregate form (43.8 g/L). After purification, the yield of rVP6 was 10.83 g/L. The rVP6 specific antiserum was confirmed by both immunofluorescent and western blotting. The versatile sequence optimization was the reason to produce a high level of rVP6 compared to other reports and can potentially apply to produce cheaper commercial kits to diagnose serological tests and new rotavirus vaccines.

Investigating Viruses During the Transformation of Molecular Biology: Part II.

My scientific career started at an extraordinary time, shortly after the discoveries of the helical structure of DNA, the central dogma of DNA to RNA to protein, and the genetic code. Part I of this series emphasizes my education and early studies highlighted by the isolation and characterization of numerous vaccinia virus enzymes, determination of the cap structure of messenger RNA, and development of poxviruses as gene expression vectors for use as recombinant vaccines. Here I describe a shift in my research focus to combine molecular biology and genetics for a comprehensive understanding of poxvirus biology. The dominant paradigm during the early years was to select a function, isolate the responsible proteins, and locate the corresponding gene, whereas later the common paradigm was to select a gene, make a mutation, and determine the altered function. Motivations, behind-the-scenes insights, importance of new technologies, and the vital roles of trainees and coworkers are emphasized.

Multiple upstream start codons (AUG) in 5' untranslated region enhance translation efficiency of cry2Ac11 without helper protein.

Cry2Ac11, a 65 kDa insecticidal protein produced by Bacillus thuringiensis, shows toxicity against dipteran and lepidopteran larvae. It is encoded by cry2Ac11 gene ( orf3), which is part of an operon comprising orf1, orf2, and orf3. Orf2, a helper protein, helps in proper folding and prevents aberrant aggregation of newly produced molecules. In this study, we have elucidated the effect of different mutations in translation initiation region (TIR), particularly the ribosomal binding site and the start codon (RBS-ATG) on cry2Ac11 gene expression without helper protein. All recombinant constructs were expressed in acrystalliferous B. thuringiensissubsp israelensis 4Q7 under the control of strong chimeric promoter cyt1AP/STAB. Of all the mutants, mut/RBS2, with two consecutive AUGs after the spacer region in TIR, exhibited 89- and 2246-fold higher transcript levels compared with 4Q7-operSalI/RBS ( cry2Ac11 operon) and 4Q7-w-RBS ( cry2Ac11 gene), respectively. The analysis of mut/RBS2 messenger RNA (mRNA) structure in the RBS-AUG region showed the presence of RBS in the single-stranded part of the moderately stable hairpin loop. The high expression efficiency of Cry2Ac11 mutant without helper protein is a cumulative and cooperative result of chimeric promoter cyt1AP/STAB-SD with the optimal context of RBS-AUG region provided by multiple AUGs and stabilizer sequence at 3' ends.

RNA Therapeutics in Cardiovascular Precision Medicine.

Since our knowledge on structure and function of messenger RNA (mRNA) has expanded from merely being an intermediate molecule between DNA and proteins to the notion that RNA is a dynamic gene regulator that can be modified and edited, RNA has become a focus of interest into developing novel therapeutic schemes. Therapeutic modulation of RNA molecules by DNA- and RNA-based therapies has broadened the scope of therapeutic targets in infectious diseases, cancer, neurodegenerative diseases and most recently in cardiovascular diseases as well. Currently, antisense oligonucleotides (ASO), small interfering RNAs (siRNAs), and microRNAs are the most widely applied therapeutic strategies to target RNA molecules and regulate gene expression and protein production. However, a number of barriers have to be overcome including instability, inadequate binding affinity and delivery to the tissues, immunogenicity, and off-target toxicity in order for these agents to evolve into efficient drugs. As cardiovascular diseases remain the leading cause of mortality worldwide, a large number of clinical trials are under development investigating the safety and efficacy of RNA therapeutics in clinical conditions such as familial hypercholesterolemia, diabetes mellitus, hypertriglyceridemia, cardiac amyloidosis, and atrial fibrillation. In this review, we summarize the clinical trials of RNA-targeting therapies in cardiovascular disease and critically discuss the advances, the outcomes, the limitations and the future directions of RNA therapeutics in precision transcriptomic medicine.

Messenger RNA Structure Regulates Translation Initiation: A Mechanism Exploited from Bacteria to Humans.

Messenger RNA Structure Regulates Translation Initiation: A Mechanism Exploited from Bacteria to Humans.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure.

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1. Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Tuning the transcription and translation of L-amino acid deaminase in Escherichia coli improves α-ketoisocaproate production from L-leucine.

α-Ketoisocaproate (KIC) is used widely in the pharmaceutical and nutraceutical industries. In previous studies, we achieved a one-step biosynthesis of KIC from l-leucine, using an Escherichia coli whole-cell biocatalyst expressing an l-amino acid deaminase (l-AAD) from Proteus vulgaris. Herein, we report the fine-tuning of l-AAD gene expression in E. coli BL21 (DE3) at the transcriptional and translational levels to improve the KIC titer. By optimizing the plasmid origin with different copy numbers, modulating messenger RNA structure downstream of the initiation codon, and designing the sequences at the ribosome binding site, we increased biocatalyst activity to 31.77%, 24.89%, and 30.20%, respectively, above that achieved with BL21/pet28a-lad. The highest KIC titers reached 76.47 g·L-1, 80.29 g·L-1, and 81.41 g·L-1, respectively. Additionally, the integration of these three engineering strategies achieved an even higher KIC production of 86.55 g·L-1 and a higher l-leucine conversion rate of 94.25%. The enzyme-engineering strategies proposed herein may be generally applicable to the construction of other biocatalysts.

EF-G catalyzed translocation dynamics in the presence of ribosomal frameshifting stimulatory signals.

Programmed -1 ribosomal frameshifting (-1PRF) is tightly regulated by messenger RNA (mRNA) sequences and structures in expressing two or more proteins with precise ratios from a single mRNA. Using single-molecule fluorescence resonance energy transfer (smFRET) between (Cy5)EF-G and (Cy3)tRNALys, we studied the translational elongation dynamics of -1PRF in the Escherichia coli dnaX gene, which contains three frameshifting signals: a slippery sequence (A AAA AAG), a Shine-Dalgarno (SD) sequence and a downstream hairpin. The frameshift promoting signals mostly impair the EF-G-catalyzed translocation step of the two tRNALys and the slippery codons from the A- and P- sites. The hairpin acts as a road block slowing the translocation rate. The upstream SD sequence together with the hairpin promotes dissociation of futile EF-G and thus causes multiple EF-G driven translocation attempts. A slippery sequence also helps dissociation of the EF-G by providing alternative base-pairing options. These results indicate that frameshifting takes place during the repetitive ribosomal conformational changes associated with EF-G dissociation upon unsuccessful translocation attempts of the second slippage codon from the A- to the P- sites.

MRNA secondary structure optimization using a correlated stem-loop prediction.

Secondary structure of messenger RNA plays an important role in the bio-synthesis of proteins. Its negative impact on translation can reduce the yield of protein by slowing or blocking the initiation and movement of ribosomes along the mRNA, becoming a major factor in the regulation of gene expression. Several algorithms can predict the formation of secondary structures by calculating the minimum free energy of RNA sequences, or perform the inverse process of obtaining an RNA sequence for a given structure. However, there is still no approach to redesign an mRNA to achieve minimal secondary structure without affecting the amino acid sequence. Here we present the first strategy to optimize mRNA secondary structures, to increase (or decrease) the minimum free energy of a nucleotide sequence, without changing its resulting polypeptide, in a time-efficient manner, through a simplistic approximation to hairpin formation. Our data show that this approach can efficiently increase the minimum free energy by >40%, strongly reducing the strength of secondary structures. Applications of this technique range from multi-objective optimization of genes by controlling minimum free energy together with CAI and other gene expression variables, to optimization of secondary structures at the genomic level.

Mutation in the 3’untranslated region of APP as a genetic determinant of cerebral amyloid angiopathy

Aβ-related cerebral amyloid angiopathy (CAA) is a major cause of primary non-traumatic brain hemorrhage. In families with an early onset of the disease, CAA can be due to amyloid precursor protein (APP) pathogenic variants or duplications. APP duplications lead to a ~1.5-fold increased APP expression, resulting in Aβ overproduction and deposition in the walls of leptomeningeal vessels. We hypothesized that rare variants in the 3'untranslated region (UTR) of APP might lead to APP overexpression in patients with CAA and no APP pathogenic variant or duplication. We performed direct sequencing of the whole APP 3'UTR in 90 patients with CAA and explored the functional consequences of one previously unreported variant. We identified three sequence variants in four patients, of which a two-base pair deletion (c.*331_*332del) was previously unannotated and absent from 175 controls of same ethnicity. This latter variant was associated with increased APP expression in vivo and in vitro. Bioinformatics and functional assays showed that the APP c.*331_*332del variant affected APP messenger RNA (mRNA) structure and binding of two microRNAs (miR-582-3p and miR-892b), providing a mechanism for the observed effects on APP expression. These results identify APP 3'UTR sequence variants as genetic determinants of Aβ-CAA.

Effects of messenger RNA structure and other translational control mechanisms on major histocompatibility complex-I mediated antigen presentation.

Effective T‐cell surveillance of antigen‐presenting cells is dependent on the expression of an array of antigenic peptides bound to major histocompatibility complex (MHC) class I (MHC‐I) or class II (MHC‐II) molecules. Pathogens co‐evolving with their hosts exploit crucial translational regulatory mechanisms in order to evade host immune recognition and thereby sustain their infection. Evasion strategies that downregulate viral protein synthesis and thereby restrict antigen presentation to cytotoxic T‐cells through the endogenous MHC‐I pathway have been implicated in the pathogenesis of viral‐associated malignancies. An understanding of the mechanisms by which messenger RNA (mRNA) structure modulates both viral mRNA translation and the antigen processing machinery to escape immune surveillance, will stimulate the development of alternative therapeutic strategies focused on RNA‐directed drugs designed to enhance immune responses against infected cells. In this review, we discuss regulatory aspects of the MHC‐I pathway and summarize current knowledge of the role attributed by mRNA structure and other translational regulatory mechanisms in immune evasion. In particular we highlight the impact of recently identified G‐quadruplex structures within virally encoded transcripts as unique regulatory signals for translational control and antigen presentation. WIREs RNA 2015, 6:157–171. doi: 10.1002/wrna.1262This article is categorized under: 1RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems2Translation > Translation Regulation

Single-Molecule Profiling of Ribosome Translational Phenomena

Translation elongation is a heterogeneous process, involving multiple compositional factors stochastically binding to the ribosome to direct protein synthesis, which in turn regulates the conformation of the ribosome. The rate of translation is often regulated by the underlying messenger RNA (mRNA) sequence and structure. Here, we use single-molecule fluorescence resonant energy transfer (FRET) and colocalization with zero-mode waveguides (ZMWs) to correlate directly ribosome conformations and compositions of thousands of ribosomes simultaneously during multiple rounds of elongation. This allows us to profile global translational rates while delineating mechanistic details of the dynamics with codon resolution. We first studied translation of a canonical mRNA with uniform translation rates to establish the tight interplay between compositional factors and conformational dynamics of the ribosome during elongation. We then determined how mRNA sequences and structures, such as hairpins and possible ribosome-mRNA base pairing, as expressed in the the dnaX −1 frameshifting sequence, perturb the basal elongation process. Our results show how mRNAs can modulate and uncouple ribosomal conformational and compositional dynamics.

Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo

RNA plays a dual role as an informational molecule and a direct effector of biological tasks. The latter function is enabled by RNA’s ability to adopt complex secondary and tertiary folds and thus has motivated extensive computational1–2 and experimental3–8 efforts for determining RNA structures. Existing approaches for evaluating RNA structure have been largely limited to in vitro systems, yet the thermodynamic forces which drive RNA folding in vitro may not be sufficient to predict stable RNA structures in vivo5. Indeed, the presence of RNA binding proteins and ATP-dependent helicases can influence which structures are present inside cells. Here we present an approach for globally monitoring RNA structure in native conditions in vivo with single nucleotide precision. This method is based on in vivo modification with dimethyl sulfate (DMS), which reacts with unpaired adenine and cytosine residues9, followed by deep sequencing to monitor modifications. Our data from yeast and mammalian cells are in excellent agreement with known mRNA structures and with the high-resolution crystal structure of the Saccharomyces cerevisiae ribosome10. Comparison between in vivo and in vitro data reveals that in rapidly dividing cells there are vastly fewer structured mRNA regions in vivo than in vitro. Even thermostable RNA structures are often denatured in cells, highlighting the importance of cellular processes in regulating RNA structure. Indeed, analysis of mRNA structure under ATP-depleted conditions in yeast reveals that energy-dependent processes strongly contribute to the predominantly unfolded state of mRNAs inside cells. Our studies broadly enable the functional analysis of physiological RNA structures and reveal that, in contrast to the Anfinsen view of protein folding, thermodynamics play an incomplete role in determining mRNA structure in vivo.

Most RNAs regulating ribosomal protein biosynthesis in Escherichia coli are narrowly distributed to Gammaproteobacteria.

In Escherichia coli, 12 distinct RNA structures within the transcripts encoding ribosomal proteins interact with specific ribosomal proteins to allow autogenous regulation of expression from large multi-gene operons, thus coordinating ribosomal protein biosynthesis across multiple operons. However, these RNA structures are typically not represented in the RNA Families Database or annotated in genomic sequences databases, and their phylogenetic distribution is largely unknown. To investigate the extent to which these RNA structures are conserved across eubacterial phyla, we created multiple sequence alignments representing 10 of these messenger RNA (mRNA) structures in E. coli. We find that while three RNA structures are widely distributed across many phyla of bacteria, seven of the RNAs are narrowly distributed to a few orders of Gammaproteobacteria. To experimentally validate our computational predictions, we biochemically confirmed dual L1-binding sites identified in many Firmicute species. This work reveals that RNA-based regulation of ribosomal protein biosynthesis is used in nearly all eubacterial phyla, but the specific RNA structures that regulate ribosomal protein biosynthesis in E. coli are narrowly distributed. These results highlight the limits of our knowledge regarding ribosomal protein biosynthesis regulation outside of E. coli, and the potential for alternative RNA structures responsible for regulating ribosomal proteins in other eubacteria.

MRNA secondary structure optimization using a correlated stem–loop prediction

Secondary structure of messenger RNA plays an important role in the bio-synthesis of proteins. Its negative impact on translation can reduce the yield of protein by slowing or blocking the initiation and movement of ribosomes along the mRNA, becoming a major factor in the regulation of gene expression. Several algorithms can predict the formation of secondary structures by calculating the minimum free energy of RNA sequences, or perform the inverse process of obtaining an RNA sequence for a given structure. However, there is still no approach to redesign an mRNA to achieve minimal secondary structure without affecting the amino acid sequence. Here we present the first strategy to optimize mRNA secondary structures, to increase (or decrease) the minimum free energy of a nucleotide sequence, without changing its resulting polypeptide, in a time-efficient manner, through a simplistic approximation to hairpin formation. Our data show that this approach can efficiently increase the minimum free energy by >40%, strongly reducing the strength of secondary structures. Applications of this technique range from multi-objective optimization of genes by controlling minimum free energy together with CAI and other gene expression variables, to optimization of secondary structures at the genomic level.

Mutational analysis of AXIN2, MSX1, and PAX9 in two Mexican oligodontia families

The genes for axin inhibition protein 2 (AXIN2), msh homeobox 1 (MSX1), and paired box gene 9 (PAX9) are involved in tooth root formation and tooth development. Mutations of the AXIN2, MSX1, and PAX9 genes are associated with non-syndromic oligodontia. In this study, we investigated phenotype and AXIN2, MSX1, and PAX9 gene variations in two Mexican families with non-syndromic oligodontia. Individuals from two families underwent clinical examinations, including an intra-oral examination and panoramic radiograph. Retrospective data were reviewed, and peripheral blood samples were collected. The exons and exon-intronic boundaries of the AXIN2, MSX1, and PAX9 genes were sequenced and analyzed. Protein and messenger RNA structures were predicted using bioinformative software programs. Clinical and oral examinations revealed isolated non-syndromic oligodontia in the two Mexican families. The average number of missing teeth was 12. The sequence analysis of exons and exon-intronic regions of AXIN2, MSX1, and PAX9 revealed 11 single-nucleotide polymorphisms (SNPs), including seven in AXIN2, two in MSX1, and three in PAX9. One novel SNP of MSX1, c.476T>G (Leu159Arg), was found in all of the studied patients in the families. MSX1 Leu159Arg and PAX9 Ala240Pro change protein and messenger RNA structures. Our findings suggested that a combined reduction of MSX1 and PAX9 gene dosages increased the risk for oligodontia in the Mexican families, as in vivo investigation has indicated that interaction between Msx1 and Pax9 is required for tooth development.