5'-UTR G-quadruplex structures acting as translational repressors

Optical Probes for Cellular Imaging of G-quadruplexes: Beyond Fluorescence Intensity Probes

The study of G-quadruplex (G4) structures that form in DNA and RNA is a rapidly growing field, which has evolved from in vitro studies of isolated G4 sequences to genome-wide detection of G4s in a cellular context. This work has revealed the tangible and significant effects that G4s may have on biological regulation. This minireview describes recent progress in the design of photoluminescent intensity-independent optical probes for G4s. We discuss the design and use of probes based on fluorescence or phosphorescence lifetime, rather than intensity-based detection; spectral ratiometric probes; and fluorescent probes for single-molecule G4-detection. We argue that each of these modalities improve unbiased G4 detection in cellular experiments, overcoming problems associated with unknown cellular uptake of probes or their organelle concentration. We discuss the improvements offered by these types of probes, as well as limitations and future research directions needed to facilitate more robust research into G4 biology.

Multi-objective computational optimization of human 5' UTR sequences

The computational design of messenger RNA (mRNA) sequences is a critical technology for both scientific research and industrial applications. Recent advances in prediction and optimization models have enabled the automatic scoring and optimization of $5^\prime $ UTR sequences, key upstream elements of mRNA. However, fully automated design of $5^\prime $ UTR sequences with more than two objective scores has not yet been explored. In this study, we present a computational pipeline that optimizes human $5^\prime $ UTR sequences in a multi-objective framework, addressing up to four distinct and conflicting objectives. Our work represents an important advancement in the multi-objective computational design of mRNA sequences, paving the way for more sophisticated mRNA engineering.

Emerging drugs for Epstein-Barr virus associated-diseases

Epstein-Barr virus (EBV) is the first identified oncogenic virus. It causes three types of diseases: lymphomas, carcinomas, and autoimmune diseases. It is estimated that two hundred thousand deaths are due to EBV each year. After a primary infection, EBV can remain latent lifelong. Reactivation to lytic phase can be induced by various drugs including small organic molecules, biologics, or a combination of both. In this review, we identified the most relevant results obtained with small organic compounds against Epstein-Barr virus-associated diseases. Specific treatments targeting Epstein-Barr Nuclear Antigen 1 are emerging concerning small organic molecules and showed promising results against several EBV-related malignancies.

Identification of G-quadruplex nucleic acid structures by high-throughput sequencing: A review

G-quadruplexes (G4s) are non-canonical nucleic acid secondary structures formed by guanine-rich DNA or RNA sequences. These structures play pivotal roles in cellular processes, including DNA replication, transcription, RNA splicing, and protein translation. High-throughput sequencing has significantly advanced the study of G4s by enabling genome-wide mapping and detailed characterization. This review provides a comprehensive overview of current methods for G4 identification using high-throughput sequencing, focusing on key techniques such as G4-seq, G4-ChIP-seq, G4-CUT&Tag, LiveG4ID-seq, G4assess, HepG4-seq, rG4-seq, RT-stop profiling with DMS-m7G footprinting, G4RP-seq, Keth-seq, and SHALIPE-seq. We discuss the principles, advantages, limitations, and applications of these methods, highlighting their contribution to our understanding of G4 biology. The review also emphasizes the need for improved tools to explore the dynamic behavior of G4s, particularly in living organisms.

5'-UTR G-Quadruplex-Mediated Translation Regulation in Eukaryotes: Current Understanding and Methodological Challenges

RNA G-quadruplexes (rG4s) in 5'-UTRs represent complex regulatory elements capable of both inhibiting and activating mRNA translation through diverse mechanisms in eukaryotes. This review analyzes the evolution of our understanding of 5'-UTR rG4-mediated translation regulation, from early discoveries of simple translation inhibitors to the current recognition of their multifaceted regulatory roles. We discuss canonical and non-canonical rG4 structures, their interactions with regulatory proteins, including helicases and FMRP, and their function in both cap-dependent and IRES-mediated translation. Special attention is given to the synergistic effects between rG4s and upstream open reading frames (uORFs), stress-responsive translation regulation, and their role in repeat-associated non-AUG (RAN) translation linked to neurodegenerative diseases. We critically evaluate methodological challenges in the field, including limitations of current detection methods, reporter system artifacts, and the necessity to verify rG4 presence in endogenous transcripts. Recent technological advances, including genome editing and high-throughput sequencing approaches, have revealed that rG4 effects are more complex and context-dependent than initially thought. This review highlights the importance of developing more robust methodologies for studying rG4s at endogenous levels and carefully reevaluating previously identified targets, while emphasizing their potential as therapeutic targets in various diseases.

Protein binding in an mRNA 5'-UTR sterically hinders translation

Structures in the 5' untranslated regions (UTRs) of mRNAs can physically modulate translation efficiency by impeding the scanning ribosome or by sequestering the translational start site. We assessed the impact of stable protein binding in 5'- and 3'-UTRs on translation efficiency by targeting the MS2 coat protein to a reporter RNA via its hairpin recognition site. Translation was assessed from the reporter RNA when coexpressed with MS2 coat proteins of varying affinities for the RNA, and at different expression levels. Binding of high-affinity proteins in the 5'-UTR hindered translation, whereas no effect was observed when the coat protein was targeted to the 3'-UTR. Inhibition of translation increased with coat protein concentration and affinity, reaching a maximum of 50%-70%. MS2 proteins engineered to bind two reporter mRNA sites had a stronger effect than those binding a single site. Our findings demonstrate that protein binding in an mRNA 5'-UTR physically impedes translation, with the effect governed by affinity, concentration, and sterics.

An interpretable RNA foundation model for exploring functional RNA motifs in plants

The complex 'language' of plant RNA encodes a vast array of biological regulatory elements that orchestrate crucial aspects of plant growth, development and adaptation to environmental stresses. Recent advancements in foundation models (FMs) have demonstrated their unprecedented potential to decipher complex 'language' in biology. In this study, we introduced PlantRNA-FM, a high-performance and interpretable RNA FM specifically designed for plants. PlantRNA-FM was pretrained on an extensive dataset, integrating RNA sequences and RNA structure information from 1,124 distinct plant species. PlantRNA-FM exhibits superior performance in plant-specific downstream tasks. PlantRNA-FM achieves an F1 score of 0.974 for genic region annotation, whereas the current best-performing model achieves 0.639. Our PlantRNA-FM is empowered by our interpretable framework that facilitates the identification of biologically functional RNA sequence and structure motifs, including both RNA secondary and tertiary structure motifs across transcriptomes. Through experimental validations, we revealed translation-associated RNA motifs in plants. Our PlantRNA-FM also highlighted the importance of the position information of these functional RNA motifs in genic regions. Taken together, our PlantRNA-FM facilitates the exploration of functional RNA motifs across the complexity of transcriptomes, empowering plant scientists with capabilities for programming RNA codes in plants.

RNA G-quadruplex formation in biologically important transcribed regions: can two-tetrad intramolecular RNA quadruplexes be formed?

G-quadruplexes (G4s) formed within RNA are emerging as promising targets for therapeutic intervention in cancer, neurodegenerative disorders and infectious diseases. Sequences containing a succession of short GG blocks, or uneven G-tract lengths unable to form three-tetrad G4s (GG motifs), are overwhelmingly more frequent than canonical motifs involving multiple GGG blocks. We recently showed that DNA is not able to form stable two-tetrad intramolecular parallel G4s. Whether RNA GG motifs can form intramolecular G4s under physiological conditions and play regulatory roles remains a burning question. In this study, we performed a systematic analysis and experimental evaluation of a number of biologically important RNA regions involving RNA GG motifs. We show that most of these motifs do not form stable intramolecular G4s but need to dimerize to form stable G4 structures. The strong tendency of RNA GG motif G4s to associate may participate in RNA-based aggregation under conditions of cellular stress.

The role of structure in regulatory RNA elements

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e., RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Neuronal activity rapidly reprograms dendritic translation via eIF4G2:uORF binding

Learning and memory require activity-induced changes in dendritic translation, but which mRNAs are involved and how they are regulated are unclear. In this study, to monitor how depolarization impacts local dendritic biology, we employed a dendritically targeted proximity labeling approach followed by crosslinking immunoprecipitation, ribosome profiling and mass spectrometry. Depolarization of primary cortical neurons with KCl or the glutamate agonist DHPG caused rapid reprogramming of dendritic protein expression, where changes in dendritic mRNAs and proteins are weakly correlated. For a subset of pre-localized messages, depolarization increased the translation of upstream open reading frames (uORFs) and their downstream coding sequences, enabling localized production of proteins involved in long-term potentiation, cell signaling and energy metabolism. This activity-dependent translation was accompanied by the phosphorylation and recruitment of the non-canonical translation initiation factor eIF4G2, and the translated uORFs were sufficient to confer depolarization-induced, eIF4G2-dependent translational control. These studies uncovered an unanticipated mechanism by which activity-dependent uORF translational control by eIF4G2 couples activity to local dendritic remodeling.

Examination of a novel expression-based gene-SNP annotation strategy to identify tissue-specific contributions to heritability in multiple traits

Complex traits show clear patterns of tissue-specific expression influenced by single nucleotide polymorphisms (SNPs), yet current strategies aggregate SNP effects to genes by employing simple physical proximity-based windows. Here, we examined whether incorporating SNPs with effects on tissue-specific cis-expression would improve our ability to detect trait-relevant tissues across 31 complex traits using stratified linkage disequilibrium score regression (S-LDSC). We found that a physical proximity annotation produced more significant tissue enrichments and larger S-LDSC regression coefficients, as compared to an expression-based annotation. Furthermore, we showed that our expression-based annotation did not outperform an annotation strategy in which an equal number of randomly chosen SNPs were annotated to genes within the same genomic window, suggesting extensive redundancy among SNP effect estimates due to linkage disequilibrium. That said, current sample sizes limit estimation of cis-genetic SNP effects; therefore, we recommend reexamination of the expression-based annotation when larger tissue-specific expression datasets become available. To examine the influence of sample size, we used a large whole blood eQTL reference panel (N = 31,684) applying a similar expression-based annotation strategy. We found that significant cis-expression QTLs in whole blood did not outperform the physical proximity annotation when estimating tissue-specific SNP heritability enrichment for either high- or low-density lipoprotein phenotypes but performed similarly for inflammatory bowel disease. Finally, we report new and updated tissue enrichment estimates across 31 complex traits, such as significant heritability enrichment of the frontal cortex for cognitive performance, educational attainment, and intelligence, providing further evidence of this structure's importance in higher cognitive function.

Stable RNA G-Quadruplex in the 5'-UTR of Human cIAP1 mRNA Promotes Translation in an IRES-Independent Manner

RNA G-quadruplex (rG4) structures can influence the fate and functions of mRNAs, especially the translation process. The presence of rG4 structures in 5'-untranslated regions (5'-UTRs) of mRNAs generally represses translation. However, rG4 structures can also promote internal ribosome entry site (IRES)-mediated translation as one of its determinants. Here, we report the identification of an evolutionary conserved rG4-forming sequence motif at the extreme 5'-end of the unusually long 5'-UTR (1.7 kb) in the transcript of human cIAP1 gene encoding the cellular inhibitor of apoptosis protein-1 that promotes cell survival by suppressing apoptosis and is overexpressed in various cancer cells. Expectedly, NMR study, CD spectroscopy, and UV melting assay confirm the formation of a potassium ion-dependent intramolecular and parallel rG4 structure at the sequence stretch. Moreover, the G4-RNA-specific precipitation using biotin-linked biomimetic BioCyTASQ validates the formation of the rG4 structure in the cIAP1 5'-UTR in cells. Interestingly, disruption of the rG4 structure in the cIAP1 5'-UTR results in a dramatic reduction in translation of the downstream luciferase reporter in cells, suggesting a translation-promoting effect of the rG4 structure, contrary to many earlier reports. Furthermore, enhancement of translation by the cIAP1 rG4 structure occurs in an IRES-independent manner.

Discovery of Novel Coumarin-quinolinium Derivatives as Pan-KRAS Translation Inhibitors by Targeting 5'-UTR RNA G-Quadruplexes

Hyperactivated KRAS mutations fuel tumorigenesis and represent attractive targets for cancer treatment. While covalent inhibitors have shown clinical benefits against the KRASG12C mutant, advancements for non-G12C mutants remain limited, highlighting the urgent demand for pan-KRAS inhibitors. RNA G-quadruplexes (rG4s) in the 5'-untranslated region of KRAS mRNA can regulate KRAS translation, making them promising targets for pan-KRAS inhibitor development. Herein, we designed and synthesized 50 novel coumarin-quinolinium derivatives, leveraging our previously developed rG4-specific ligand, QUMA-1. Notably, several compounds exhibited potent antiproliferative activity against cancer cells as pan-KRAS translation inhibitors. Among them, 15a displayed exceptional capability in stabilizing KRAS rG4s, suppressing KRAS translation, and consequently modulating MAPK and PI3K-AKT pathways. 15a induced cell cycle arrest, prompted apoptosis in KRAS-driven cancer cells, and effectively inhibited tumor growth in a KRAS mutant xenograft model. These findings underscore the potential of 15a as a pan-KRAS translation inhibitor, offering a novel and promising approach to target various KRAS-driven cancers.

TMPRSS2 isoform 1 downregulation by G-quadruplex stabilization induces SARS-CoV-2 replication arrest

BACKGROUND

SARS-CoV-2 infection depends on the host cell factors angiotensin-converting enzyme 2, ACE2, and the transmembrane serinprotease 2, TMPRSS2. Potential inhibitors of these proteins would be ideal targets against severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) infection. Our data opens the possibility that changes within TMPRSS2 can modulate the outcome during a SARS-CoV-2 infection.

RESULTS

We reveal that TMPRSS2 acts not only during viral entry but has also an important role during viral replication. In addition to previous functions for TMPRSS2 during viral entry, we determined by specific downregulation of distinct isoforms that only isoform 1 controls and supports viral replication. G-quadruplex (G4) stabilization by chemical compounds impacts TMPRSS2 gene expression. Here we extend and in-depth characterize these observations and identify that a specific G4 in the first exon of the TMPRSS2 isoform 1 is particular targeted by the G4 ligand and affects viral replication. Analysis of potential single nucleotide polymorphisms (SNPs) reveals that a reported SNP at this G4 in isoform 1 destroys the G4 motif and makes TMPRSS2 ineffective towards G4 treatment.

CONCLUSION

These findings uncover a novel mechanism in which G4 stabilization impacts SARS-CoV-2 replication by changing TMPRSS2 isoform 1 gene expression.

Genetic variations in G-quadruplex forming sequences affect the transcription of human disease-related genes

Guanine-rich DNA strands can fold into non-canonical four-stranded secondary structures named G-quadruplexes (G4s). G4s folded in proximal promoter regions (PPR) are associated either with positive or negative transcriptional regulation. Given that single nucleotide variants (SNVs) affecting G4 folding (G4-Vars) may alter gene transcription, and that SNVs are associated with the human diseases' onset, we undertook a novel comprehensive study of the G4-Vars genome-wide (G4-variome) to find disease-associated G4-Vars located into PPRs. We developed a bioinformatics strategy to find disease-related SNVs located into PPRs simultaneously overlapping with putative G4-forming sequences (PQSs). We studied five G4-Vars disturbing in vitro the folding and stability of the G4s located into PPRs, which had been formerly associated with sporadic Alzheimer's disease (GRIN2B), a severe familiar coagulopathy (F7), atopic dermatitis (CSF2), myocardial infarction (SIRT1) and deafness (LHFPL5). Results obtained in cultured cells for these five G4-Vars suggest that the changes in the G4s affect the transcription, potentially contributing to the development of the mentioned diseases. Collectively, data reinforce the general idea that G4-Vars may impact on the different susceptibilities to human genetic diseases' onset, and could be novel targets for diagnosis and drug design in precision medicine.

Analysis of Nucleotide Variations in Human G-Quadruplex Forming Regions Associated with Disease States

While the role of G quadruplex (G4) structures has been identified in cancers and metabolic disorders, single nucleotide variations (SNVs) and their effect on G4s in disease contexts have not been extensively studied. The COSMIC and CLINVAR databases were used to detect SNVs present in G4s to identify sequence level changes and their effect on the alteration of the G4 secondary structure. A total of 37,515 G4 SNVs in the COSMIC database and 2378 in CLINVAR were identified. Of those, 7236 COSMIC (19.3%) and 457 (19%) of the CLINVAR variants result in G4 loss, while 2728 (COSMIC) and 129 (CLINVAR) SNVs gain a G4 structure. The remaining variants potentially affect the folding energy without affecting the presence of a G4. Analysis of mutational patterns in the G4 structure shows a higher selective pressure (3-fold) in the coding region on the template strand compared to the reverse strand. At the same time, an equal proportion of SNVs were observed among intronic, promoter, and enhancer regions across strands.

Effects of Combinations of Untranslated-Region Sequences on Translation of mRNA

mRNA-based therapeutics have been found to be a promising treatment strategy in immunotherapy, gene therapy, and cancer treatments. Effectiveness of mRNA therapeutics depends on the level and duration of a desired protein's expression, which is determined by various cis- and trans-regulatory elements of the mRNA. Sequences of 5' and 3' untranslated regions (UTRs) are responsible for translational efficiency and stability of mRNA. An optimal combination of the regulatory sequences allows researchers to significantly increase the target protein's expression. Using both literature data and previously obtained experimental data, we chose six sequences of 5'UTRs (adenoviral tripartite leader [TPL], HBB, rabbit β-globin [Rabb], H4C2, Moderna, and Neo2) and five sequences of 3'UTRs (mtRNR-EMCV, mtRNR-AES, mtRNR-mtRNR, BioNTech, and Moderna). By combining them, we constructed 30 in vitro transcribed RNAs encoding firefly luciferase with various combinations of 5'- and 3'UTRs, and the resultant bioluminescence was assessed in the DC2.4 cell line at 4, 8, 24, and 72 h after transfection. The cellular data enabled us to identify the best seven combinations of 5'- and 3'UTRs, whose translational efficiency was then assessed in BALB/c mice. Two combinations of 5'- and 3'UTRs (5'Rabb-3'mtRNR-EMCV and 5'TPL-3'Biontech) led to the most pronounced increase in the luciferase amount in the in vivo experiment in mice. Subsequent analysis of the stability of the mRNA indicated that the increase in luciferase expression is explained primarily by the efficiency of translation, not by the number of RNA molecules. Altogether, these findings suggest that 5'UTR-and-3'UTR combinations 5'Rabb-3'mtRNR- EMCV and 5'TPL-3'Biontech lead to high expression of target proteins and may be considered for use in preventive and therapeutic modalities based on mRNA.

Alternative splicing of BCL-x is controlled by RBM25 binding to a G-quadruplex in BCL-x pre-mRNA

BCL-x is a master regulator of apoptosis whose pre-mRNA is alternatively spliced into either a long (canonical) anti-apoptotic Bcl-xL isoform, or a short (alternative) pro-apoptotic Bcl-xS isoform. The balance between these two antagonistic isoforms is tightly regulated and overexpression of Bcl-xL has been linked to resistance to chemotherapy in several cancers, whereas overexpression of Bcl-xS is associated to some forms of diabetes and cardiac disorders. The splicing factor RBM25 controls alternative splicing of BCL-x: its overexpression favours the production of Bcl-xS, whereas its downregulation has the opposite effect. Here we show that RBM25 directly and specifically binds to GQ-2, an RNA G-quadruplex (rG4) of BCL-x pre-mRNA that forms at the vicinity of the alternative 5' splice site leading to the alternative Bcl-xS isoform. This RBM25/rG4 interaction is crucial for the production of Bcl-xS and depends on the RE (arginine-glutamate-rich) motif of RBM25, thus defining a new type of rG4-interacting domain. PhenDC3, a benchmark G4 ligand, enhances the binding of RBM25 to the GQ-2 rG4 of BCL-x pre-mRNA, thereby promoting the alternative pro-apoptotic Bcl-xS isoform and triggering apoptosis. Furthermore, the screening of a combinatorial library of 90 putative G4 ligands led to the identification of two original compounds, PhenDH8 and PhenDH9, superior to PhenDC3 in promoting the Bcl-xS isoform and apoptosis. Thus, favouring the interaction between RBM25 and the GQ-2 rG4 of BCL-x pre-mRNA represents a relevant intervention point to re-sensitize cancer cells to chemotherapy.

Development of a highly optimized procedure for the discovery of RNA G-quadruplexes by combining several strategies

RNA G-quadruplexes (rG4s) are non-canonical secondary structures that are formed by the self-association of guanine quartets and that are stabilized by monovalent cations (e.g. potassium). rG4s are key elements in several post-transcriptional regulation mechanisms, including both messenger RNA (mRNA) and microRNA processing, mRNA transport and translation, to name but a few examples. Over the past few years, multiple high-throughput approaches have been developed in order to identify rG4s, including bioinformatic prediction, in vitro assays and affinity capture experiments coupled to RNA sequencing. Each individual approach had its limits, and thus yielded only a fraction of the potential rG4 that are further confirmed (i.e., there is a significant level of false positive). This report aims to benefit from the strengths of several existing approaches to identify rG4s with a high potential of being folded in cells. Briefly, rG4s were pulled-down from cell lysates using the biotinylated biomimetic G4 ligand BioTASQ and the sequences thus isolated were then identified by RNA sequencing. Then, a novel bioinformatic pipeline that included DESeq2 to identify rG4 enriched transcripts, MACS2 to identify rG4 peaks, rG4-seq to increase rG4 formation probability and G4RNA Screener to detect putative rG4s was performed. This workflow uncovers new rG4 candidates whose rG4-folding was then confirmed in vitro using an array of established biophysical methods. Clearly, this workflow led to the identification of novel rG4s in a highly specific and reliable manner.

The hide-and-seek game of the oncogenic Epstein-Barr virus-encoded EBNA1 protein with the immune system: An RNA G-quadruplex tale

The Epstein-Barr virus (EBV) is the first oncogenic virus described in human. EBV infects more than 90% of the human population worldwide, but most EBV infections are asymptomatic. After the primary infection, the virus persists lifelong in the memory B cells of the infected individuals. Under certain conditions the virus can cause several human cancers, that include lymphoproliferative disorders such as Burkitt and Hodgkin lymphomas and non-lymphoid malignancies such as 100% of nasopharyngeal carcinoma and 10% of gastric cancers. Each year, about 200,000 EBV-related cancers emerge, hence accounting for at least 1% of worldwide cancers. Like all gammaherpesviruses, EBV has evolved a strategy to escape the host immune system. This strategy is mainly based on the tight control of the expression of its Epstein-Barr nuclear antigen-1 (EBNA1) protein, the EBV-encoded genome maintenance protein. Indeed, EBNA1 is essential for viral genome replication and maintenance but, at the same time, is also highly antigenic and T cells raised against EBNA1 exist in infected individuals. For this reason, EBNA1 is considered as the Achilles heel of EBV and the virus has seemingly evolved a strategy that employs the binding of nucleolin, a host cell factor, to RNA G-quadruplex (rG4) within EBNA1 mRNA to limit its expression to the minimal level required for function while minimizing immune recognition. This review recapitulates in a historical way the knowledge accumulated on EBNA1 immune evasion and discusses how this rG4-dependent mechanism can be exploited as an intervention point to unveil EBV-related cancers to the immune system.

G-quadruplex in hepatitis B virus pregenomic RNA promotes its translation

Hepatitis B virus (HBV) is a hepatotropic DNA virus that has a very compact genome. Due to this genomic density, several distinct mechanisms are used to facilitate the viral life cycle. Recently, accumulating evidence show that G-quadruplex (G4) in different viruses play essential regulatory roles in key steps of the viral life cycle. Although G4 structures in the HBV genome have been reported, their function in HBV replication remains elusive. In this study, we treated an HBV replication-competent cell line and HBV-infected cells with the G4 structure stabilizer pyridostatin (PDS) and evaluated different HBV replication markers to better understand the role played by the G4. In both models, we found PDS had no effect on viral precore RNA (pcRNA) or pre-genomic RNA (pgRNA), but treatment did increase HBeAg/HBc ELISA reads and intracellular levels of viral core/capsid protein (HBc) in a dose-dependent manner, suggesting post-transcriptional regulation. To further dissect the mechanism of G4 involvement, we used in vitro-synthesized HBV pcRNA and pgRNA. Interestingly, we found PDS treatment only enhanced HBc expression from pgRNA but not HBeAg expression from pcRNA. Our bioinformatic analysis and CD spectroscopy revealed that pgRNA harbors a conserved G4 structure. Finally, we introduced point mutations in pgRNA to disrupt its G4 structure and observed the resulting mutant failed to respond to PDS treatment and decreased HBc level in in vitro translation assay. Taken together, our data demonstrate that HBV pgRNA contains a G4 structure that plays a vital role in the regulation of viral mRNA translation.

The use of RNA-based treatments in the field of cancer immunotherapy

Over the past several decades, mRNA vaccines have evolved from a theoretical concept to a clinical reality. These vaccines offer several advantages over traditional vaccine techniques, including their high potency, rapid development, low-cost manufacturing, and safe administration. However, until recently, concerns over the instability and inefficient distribution of mRNA in vivo have limited their utility. Fortunately, recent technological advancements have mostly resolved these concerns, resulting in the development of numerous mRNA vaccination platforms for infectious diseases and various types of cancer. These platforms have shown promising outcomes in both animal models and humans. This study highlights the potential of mRNA vaccines as a promising alternative approach to conventional vaccine techniques and cancer treatment. This review article aims to provide a thorough and detailed examination of mRNA vaccines, including their mechanisms of action and potential applications in cancer immunotherapy. Additionally, the article will analyze the current state of mRNA vaccine technology and highlight future directions for the development and implementation of this promising vaccine platform as a mainstream therapeutic option. The review will also discuss potential challenges and limitations of mRNA vaccines, such as their stability and in vivo distribution, and suggest ways to overcome these issues. By providing a comprehensive overview and critical analysis of mRNA vaccines, this review aims to contribute to the advancement of this innovative approach to cancer treatment.

Investigating the NRAS 5' UTR as a target for small molecules

Neuroblastoma RAS (NRAS) is an oncogene that is deregulated and highly mutated in cancers including melanomas and acute myeloid leukemias. The 5' untranslated region (UTR) (5' UTR) of the NRAS mRNA contains a G-quadruplex (G4) that regulates translation. Here we report a novel class of small molecule that binds to the G4 structure located in the 5' UTR of the NRAS mRNA. We used a small molecule microarray screen to identify molecules that selectively bind to the NRAS-G4 with submicromolar affinity. One compound inhibits the translation of NRAS in vitro but showed only moderate effects on the NRAS levels in cellulo. Rapid Amplification of cDNA Ends and RT-PCR analysis revealed that the predominant NRAS transcript does not possess the G4 structure. Thus, although NRAS transcripts lack a G4 in many cell lines the concept of targeting folded regions within 5' UTRs to control translation remains a highly attractive strategy.

Structural and Functional Classification of G-Quadruplex Families within the Human Genome

G-quadruplexes (G4s) are short secondary DNA structures located throughout genomic DNA and transcribed RNA. Although G4 structures have been shown to form in vivo, no current search tools that examine these structures based on previously identified G-quadruplexes and filter them based on similar sequence, structure, and thermodynamic properties are known to exist. We present a framework for clustering G-quadruplex sequences into families using the CD-HIT, MeShClust, and DNACLUST methods along with a combination of Starcode and BLAST. Utilizing this framework to filter and annotate clusters, 95 families of G-quadruplex sequences were identified within the human genome. Profiles for each family were created using hidden Markov models to allow for the identification of additional family members and generate homology probability scores. The thermodynamic folding energy properties, functional annotation of genes associated with the sequences, scores from different prediction algorithms, and transcription factor binding motifs within a family were used to annotate and compare the diversity within and across clusters. The resulting set of G-quadruplex families can be used to further understand how different regions of the genome are regulated by factors targeting specific structures common to members of a specific cluster.

Demir G, Demir D. Quantum tunneling time delay investigation of K K + ion in human telomeric G-quadruplex systems

Guanine-rich quadruplex DNA (G-quadruplex) is of interest both in cell biology and nanotechnology. Its biological functions necessitate a G-quadruplex to be stabilized against escape of the monovalent metal cations. The potassium ion (K K + is particularly important as it experiences a potential energy barrier while it enters and exits the G-quadruplex systems which are normally found in human telomere. In the present work, we analyzed the time it takes for theK K + cations to get in and out of the G-quadruplex. Our time estimate is based on entropic tunneling time-a time formula which gave biologically relevant results for DNA point mutation by proton tunneling. The potential energy barrier experienced byK K + ions is determined from a quantum mechanical simulation study, Schrodinger equation is solved using MATLAB, and the computed eigenfunctions and eigenenergies are used in the entropic tunneling time formula to compute the time delay and charge accumulation rate during the tunneling ofK K + in G-quadruplex. The computations have shown that ion tunneling takes picosecond times. In addition, averageK K + accumulation rate is found to be in the picoampere range. Our results show that time delay during theK K + ion tunneling is in the ballpark of the conformational transition times in biological systems, and it could be an important parameter for understanding its biological role in human DNA as well as for the possible applications in biotechnology. To our knowledge, for the first time in the literature, time delay during the ion tunneling from and into G-quadruplexes is computed.

Optimization of G-Quadruplex Ligands through a SAR Study Combining Parallel Synthesis and Screening of Cationic Bis(acylhydrazones)

G-quadruplexes (G4s), secondary structures adopted by guanine-rich DNA and RNA sequences, are implicated in numerous biological processes and have been suggested as potential drug targets. Accordingly, there is an increasing interest in developing high-throughput methods that allow the generation of congeneric series of G4-targeting molecules ("ligands") and investigating their interactions with the targets. We have developed an operationally simple method of parallel synthesis to generate "ready-to-screen" libraries of cationic acylhydrazones, a motif that we have previously identified as a promising scaffold for potent, biologically active G4 ligands. Combined with well-established screening techniques, such as fluorescence melting, this method enables the rapid synthesis and screening of combinatorial libraries of potential G4 ligands. Following this protocol, we synthesized a combinatorial library of 90 bis(acylhydrazones) and screened it against five different nucleic acid structures. This way, we were able to analyze the structure-activity relationships within this series of G4 ligands, and identified three novel promising ligands whose interactions with G4-DNAs of different topologies were studied in detail by a combination of several biophysical techniques, including native mass spectrometry, and molecular modeling.

G4Atlas: a comprehensive transcriptome-wide G-quadruplex database

RNA G-quadruplex (rG4) is a vital RNA tertiary structure motif that involves the base pairs on both Hoogsteen and Watson-Crick faces of guanines. rG4 is of great importance in the post-transcriptional regulation of gene expression. Experimental technologies have advanced to identify in vitro and in vivo rG4s across diverse transcriptomes. Building on these recent advances, here we present G4Atlas, the first transcriptome-wide G-quadruplex database, in which we have collated, classified, and visualized transcriptome rG4 experimental data, generated from rG4-seq, chemical profiling and ligand-binding methods. Our comprehensive database includes transcriptome-wide rG4s generated from 82 experimental treatments and 238 samples across ten species. In addition, we have also included RNA secondary structure prediction information across both experimentally identified and unidentified rG4s to enable users to display any potential competitive folding between rG4 and RNA secondary structures. As such, G4Atlas will enable users to explore the general functions of rG4s in diverse biological processes. In addition, G4Atlas lays the foundation for further data-driven deep learning algorithms to examine rG4 structural features.

Ribo-uORF: a comprehensive data resource of upstream open reading frames (uORFs) based on ribosome profiling

Upstream open reading frames (uORFs) are typically defined as translation sites located within the 5' untranslated region upstream of the main protein coding sequence (CDS) of messenger RNAs (mRNAs). Although uORFs are prevalent in eukaryotic mRNAs and modulate the translation of downstream CDSs, a comprehensive resource for uORFs is currently lacking. We developed Ribo-uORF (http://rnainformatics.org.cn/RiboUORF) to serve as a comprehensive functional resource for uORF analysis based on ribosome profiling (Ribo-seq) data. Ribo-uORF currently supports six species: human, mouse, rat, zebrafish, fruit fly, and worm. Ribo-uORF includes 501 554 actively translated uORFs and 107 914 upstream translation initiation sites (uTIS), which were identified from 1495 Ribo-seq and 77 quantitative translation initiation sequencing (QTI-seq) datasets, respectively. We also developed mRNAbrowse to visualize items such as uORFs, cis-regulatory elements, genetic variations, eQTLs, GWAS-based associations, RNA modifications, and RNA editing. Ribo-uORF provides a very intuitive web interface for conveniently browsing, searching, and visualizing uORF data. Finally, uORFscan and UTR5var were developed in Ribo-uORF to precisely identify uORFs and analyze the influence of genetic mutations on uORFs using user-uploaded datasets. Ribo-uORF should greatly facilitate studies of uORFs and their roles in mRNA translation and posttranscriptional control of gene expression.

Aberrant splicing in human cancer: An RNA structural code point of view

Alternative splicing represents an essential process that occurs widely in eukaryotes. In humans, most genes undergo alternative splicing to ensure transcriptome and proteome diversity reflecting their functional complexity. Over the last decade, aberrantly spliced transcripts due to mutations in cis- or trans-acting splicing regulators have been tightly associated with cancer development, largely drawing scientific attention. Although a plethora of single proteins, ribonucleoproteins, complexed RNAs, and short RNA sequences have emerged as nodal contributors to the splicing cascade, the role of RNA secondary structures in warranting splicing fidelity has been underestimated. Recent studies have leveraged the establishment of novel high-throughput methodologies and bioinformatic tools to shed light on an additional layer of splicing regulation in the context of RNA structural elements. This short review focuses on the most recent available data on splicing mechanism regulation on the basis of RNA secondary structure, emphasizing the importance of the complex RNA G-quadruplex structures (rG4s), and other specific RNA motifs identified as splicing silencers or enhancers. Moreover, it intends to provide knowledge on newly established techniques that allow the identification of RNA structural elements and highlight the potential to develop new RNA-oriented therapeutic strategies against cancer.

Type I arginine methyltransferases are intervention points to unveil the oncogenic Epstein-Barr virus to the immune system

The oncogenic Epstein-Barr virus (EBV) evades the immune system but has an Achilles heel: its genome maintenance protein EBNA1. Indeed, EBNA1 is essential for viral genome maintenance but is also highly antigenic. Hence, EBV seemingly evolved a system in which the glycine-alanine repeat (GAr) of EBNA1 limits the translation of its own mRNA to the minimal level to ensure its essential function, thereby, at the same time, minimizing immune recognition. Therefore, defining intervention points at which to interfere with GAr-based inhibition of translation is an important step to trigger an immune response against EBV-carrying cancers. The host protein nucleolin (NCL) plays a critical role in this process via a direct interaction with G-quadruplexes (G4) formed in the GAr-encoding sequence of the viral EBNA1 mRNA. Here we show that the C-terminal arginine-glycine-rich (RGG) motif of NCL is crucial for its role in GAr-based inhibition of translation by mediating interaction of NCL with G4 of EBNA1 mRNA. We also show that this interaction depends on the type I arginine methyltransferase family, notably PRMT1 and PRMT3: drugs or small interfering RNA that target these enzymes prevent efficient binding of NCL on G4 of EBNA1 mRNA and relieve GAr-based inhibition of translation and of antigen presentation. Hence, this work defines type I arginine methyltransferases as therapeutic targets to interfere with EBNA1 and EBV immune evasion.

GAIA: G-quadruplexes in alive creature database

G-quadruplexes (G4) are 3D structures that are found in both DNA and RNA. Interest in this structure has grown over the past few years due to both its implication in diverse biological mechanisms and its potential use as a therapeutic target, to name two examples. G4s in humans have been widely studied; however, the level of their study in other species remains relatively minimal. That said, progress in this field has resulted in the prediction of G4s structures in various species, ranging from bacteria to eukaryotes. These predictions were analysed in a previous study which revealed that G4s are present in all living kingdoms. To date, eleven different databases have grouped the various G4s depending on either their structures, on the proteins that might bind them, or on their location in the various genomes. However, none of these databases contains information on their location in the transcriptome of many of the implicated species. The GAIA database was designed so as to make this data available online in a user-friendly manner. Through its web interface, users can query GAIA to filter G4s, which, we hope, will help the research in this field. GAIA is available at: https://gaia.cobius.usherbrooke.ca.

Small molecule-based detection of non-canonical RNA G-quadruplex structures that modulate protein translation

Tandem repeats of guanine-rich sequences in RNA often form thermodynamically stable four-stranded RNA structures. Such RNA G-quadruplexes have long been considered to be linked to essential biological processes, yet their physiological significance in cells remains unclear. Here, we report a approach that permits the detection of RNA G-quadruplex structures that modulate protein translation in mammalian cells. The approach combines antibody arrays and RGB-1, a small molecule that selectively stabilizes RNA G-quadruplex structures. Analysis of the protein and mRNA products of 84 cancer-related human genes identified Nectin-4 and CapG as G-quadruplex-controlled genes whose mRNAs harbor non-canonical G-quadruplex structures on their 5′UTR region. Further investigations revealed that the RNA G-quadruplex of CapG exhibits a structural polymorphism, suggesting a possible mechanism that ensures the translation repression in a KCl concentration range of 25–100 mM. The approach described in the present study sets the stage for further discoveries of RNA G-quadruplexes.

Role of single nucleotide polymorphisms of the HSD3B1 gene (rs6203 and rs33937873) in the prediction of prostate cancer risk

3-β-hydroxysteroid dehydrogenase 1 (HSD3B1) is shown to affect dihydrotestosterone level in prostatic tissue which is a risk factor for prostate cancer (PC). The present study aimed to determine whether rs33937873 (G313A) and rs6203 (C338T) single nucleotide polymorphisms (SNP) in HSD3B1 gene was a potential risk factor for PC susceptibility and can predict the recurrence of PC in Egyptian patients. A total of 186 Egyptian patients were selected with incident primary PC and compared with 180 age healthy controls. The frequencies and the main effect of rs33937873 and rs6203 in HSD3B1 were compared and investigated between the patients and control using genotyping technique and statistical analysis. The mutant GA genotype of G313A in rs33937873 SNP was considered as an independent risk for PC in the multivariate regression analysis [odds ratio (OR)=2.7, 95% confidence intervals (CI): 1.2-5.5, P=0.01] together with positive history of hypertension (HTN) (OR=6.2, 95% CI: 3.2-12.1, P=0.0001) and begin prostatic hyperplasia (BPH; OR=8.9, 95% CI: 4.5-17.5, P=0.0001). Conversely, in rs6203 (C338T), C allele is considered as major risk allele in the development of PC (OR=1.8, 95% CI: 1.3-2.4, P=0.0003). The univariate logistic regression analyses indicated that CC genotype of rs6203 was a PC risk factor (OR=1.9, 95% CI: 1.3-2.9, P=0.002). In addition, the frequency of the A-C haplotype established by rs33937873-rs6203 was also significantly higher for PC (P=0.013). The predication of PC recurrence was associated only with positive family history (OR=7.7, 95% CI: 2.3-25.9, P=0.001) and not for The G313A and C338T SNPs. These results suggested that the two HSD3B1 polymorphisms rs33937873 and rs6203 may modify the risk of PC, particularly among patients with HTN and history of BPH, suggesting them as prominent future markers for prediction of PC risk.

To unwind the biological knots: The DNA/RNA G-quadruplex resolvase RHAU (DHX36) in development and disease

The G-quadruplex (G4) sequences are short fragments of 4-interval triple guanine (G) with frequent and ubiquitous distribution in the genome and RNA transcripts. The G4 sequences are usually folded into secondary "knot" structure via Hoogsteen hydrogen bond to exert negative regulation on a variety of biological processes, including DNA replication and transcription, mRNA translation, and telomere maintenance. Recent structural biological and mouse genetics studies have demonstrated that RHAU (DHX36) can bind and unwind the G4 "knots" to modulate embryonic development and postnatal organ function. Deficiency of RHAU gives rise to embryonic lethality, impaired organogenesis, and organ dysfunction. These studies uncovered the pivotal G4 resolvase function of RHAU to release the G4 barrier, which plays fundamental roles in development and physiological homeostasis. This review discusses the latest advancements and findings in deciphering RHAU functions using animal models.

G-quadruplex occurrence and conservation: more than just a question of guanine–cytosine content

G-quadruplexes are motifs found in DNA and RNA that can fold into tertiary structures. Until now, they have been studied experimentally mainly in humans and a few other species. Recently, predictions have been made with bacterial and archaeal genomes. Nevertheless, a global comparison of predicted G4s (pG4s) across and within the three living kingdoms has not been addressed. In this study, we aimed to predict G4s in genes and transcripts of all kingdoms of living organisms and investigated the differences in their distributions. The relation of the predictions with GC content was studied. It appears that GC content is not the only parameter impacting G4 predictions and abundance. The distribution of pG4 densities varies depending on the class of transcripts and the group of species. Indeed, we have observed that, in coding transcripts, there are more predicted G4s than expected for eukaryotes but not for archaea and bacteria, while in noncoding transcripts, there are as many or fewer predicted G4s in all species groups. We even noticed that some species with the same GC content presented different pG4 profiles. For instance, Leishmania major and Chlamydomonas reinhardtii both have 60% of GC content, but the former has a pG4 density of 0.07 and the latter 1.16.

Reversal of G-Quadruplexes' Role in Translation Control When Present in the Context of an IRES

G-quadruplexes (GQs) are secondary nucleic acid structures that play regulatory roles in various cellular processes. G-quadruplex-forming sequences present within the 5' UTR of mRNAs can function not only as repressors of translation but also as elements required for optimum function. Based upon previous reports, the majority of the 5' UTR GQ structures inhibit translation, presumably by blocking the ribosome scanning process that is essential for detection of the initiation codon. However, there are certain mRNAs containing GQs that have been identified as positive regulators of translation, as they are needed for translation initiation. While most cellular mRNAs utilize the 5' cap structure to undergo cap-dependent translation initiation, many rely on cap-independent translation under certain conditions in which the cap-dependent initiation mechanism is not viable or slowed down, for example, during development, under stress and in many diseases. Cap-independent translation mainly occurs via Internal Ribosomal Entry Sites (IRESs) that are located in the 5' UTR of mRNAs and are equipped with structural features that can recruit the ribosome or other factors to initiate translation without the need for a 5' cap. In this review, we will focus only on the role of RNA GQs present in the 5' UTR of mRNAs, where they play a critical role in translation initiation, and discuss the potential mechanism of this phenomenon, which is yet to be fully delineated.

Identification and targeting of G-quadruplex structures in MALAT1 long non-coding RNA

RNA G-quadruplexes (rG4s) have functional roles in many cellular processes in diverse organisms. While a number of rG4 examples have been reported in coding messenger RNAs (mRNA), so far only limited works have studied rG4s in non-coding RNAs (ncRNAs), especially in long non-coding RNAs (lncRNAs) that are of emerging interest and significance in biology. Herein, we report that MALAT1 lncRNA contains conserved rG4 motifs, forming thermostable rG4 structures with parallel topology. We also show that rG4s in MALAT1 lncRNA can interact with NONO protein with high specificity and affinity in vitro and in nuclear cell lysate, and we provide cellular data to support that NONO protein recognizes MALAT1 lncRNA via rG4 motifs. Notably, we demonstrate that rG4s in MALAT1 lncRNA can be targeted by the rG4-specific small molecule, peptide, and L-aptamer, leading to the dissociation of MALAT1 rG4-NONO protein interaction. Altogether, this study uncovers new and important rG4s in MALAT1 lncRNAs, reveals their specific interactions with NONO protein, offers multiple strategies for targeting MALAT1 and its RNA-protein complex via its rG4 structure and illustrates the prevalence and significance of rG4s in ncRNAs.

The different activities of RNA G-quadruplex structures are controlled by flanking sequences

The role of G-quadruplex (G4) RNA structures is multifaceted and controversial. Here, we have used as a model the EBV-encoded EBNA1 and the Kaposi's sarcoma-associated herpesvirus (KSHV)-encoded LANA1 mRNAs. We have compared the G4s in these two messages in terms of nucleolin binding, nuclear mRNA retention, and mRNA translation inhibition and their effects on immune evasion. The G4s in the EBNA1 message are clustered in one repeat sequence and the G4 ligand PhenDH2 prevents all G4-associated activities. The RNA G4s in the LANA1 message take part in similar multiple mRNA functions but are spread throughout the message. The different G4 activities depend on flanking coding and non-coding sequences and, interestingly, can be separated individually. Together, the results illustrate the multifunctional, dynamic and context-dependent nature of G4 RNAs and highlight the possibility to develop ligands targeting specific RNA G4 functions. The data also suggest a common multifunctional repertoire of viral G4 RNA activities for immune evasion.

Sequence determinants as key regulators in gene expression of T cells

T cell homeostasis, T cell differentiation, and T cell effector function rely on the constant fine-tuning of gene expression. To alter the T cell state, substantial remodeling of the proteome is required. This remodeling depends on the intricate interplay of regulatory mechanisms, including post-transcriptional gene regulation. In this review, we discuss how the sequence of a transcript influences these post-transcriptional events. In particular, we review how sequence determinants such as sequence conservation, GC content, and chemical modifications define the levels of the mRNA and the protein in a T cell. We describe the effect of different forms of alternative splicing on mRNA expression and protein production, and their effect on subcellular localization. In addition, we discuss the role of sequences and structures as binding hubs for miRNAs and RNA-binding proteins in T cells. The review thus highlights how the intimate interplay of post-transcriptional mechanisms dictate cellular fate decisions in T cells.

Interaction of 9-Methoxyluminarine with Different G-Quadruplex Topologies: Fluorescence and Circular Dichroism Studies

The interactions of G-quadruplexes of different topologies with highly fluorescent 9-methoxyluminarine ligand 9-MeLM were investigated by fluorescence and circular dichroism spectroscopy. The results showed that 9-methoxyluminarine was able to interact and did not destabilize any investigated molecular targets. The studied compound was selectively quenched by parallel c-MYC G-quadruplex DNA, whereas hybrid and antiparallel G4 topology caused only a negligible decrease in the fluorescence of the ligand. A high decrease of fluorescence of the ligand after binding with c-MYC G-quadruplex suggests that this molecule can be used as a selective probe for parallel G-quadruplexes.

Unwinding the roles of RNA helicase MOV10

MOV10 is an RNA helicase that associates with the RNA‐induced silencing complex component Argonaute (AGO), likely resolving RNA secondary structures. MOV10 also binds the Fragile X mental retardation protein to block AGO2 binding at some sites and associates with UPF1, a principal component of the nonsense‐mediated RNA decay pathway. MOV10 is widely expressed and has a key role in the cellular response to viral infection and in suppressing retrotransposition. Posttranslational modifications of MOV10 include ubiquitination, which leads to stimulation‐dependent degradation, and phosphorylation, which has an unknown function. MOV10 localizes to the nucleus and/or cytoplasm in a cell type‐specific and developmental stage‐specific manner. Knockout of Mov10 leads to embryonic lethality, underscoring an important role in development where it is required for the completion of gastrulation. MOV10 is expressed throughout the organism; however, most studies have focused on germline cells and neurons. In the testes, the knockdown of Mov10 disrupts proliferation of spermatogonial progenitor cells. In brain, MOV10 is significantly elevated postnatally and binds mRNAs encoding cytoskeleton and neuron projection proteins, suggesting an important role in neuronal architecture. Heterozygous Mov10 mutant mice are hyperactive and anxious and their cultured hippocampal neurons have reduced dendritic arborization. Zygotic knockdown of Mov10 in Xenopus laevis causes abnormal head and eye development and mislocalization of neuronal precursors in the brain. Thus, MOV10 plays a vital role during development, defense against viral infection and in neuronal development and function: its many roles and regulation are only beginning to be unraveled.

This article is categorized under:

RNA Interactions with Proteins and Other Molecules > RNA‐Protein Complexes

RNA Interactions with Proteins and Other Molecules > Protein‐RNA Interactions: Functional Implications

Detecting and Profiling Endogenous RNA G-Quadruplexes in the Human Transcriptome

G-quadruplexes are the non-canonical nucleic acid structures that are preferentially formed in G-rich regions. This structure has been shown to be associated with many biological functions. Regardless of the broad efforts on DNA G-quadruplexes, we still have limited knowledge on RNA G-quadruplexes, especially in a transcriptome-wide manner. Herein, by integrating the DMS-seq and the bioinformatics pipeline, we profiled and depicted the RNA G-quadruplexes in the human transcriptome. The genes that contain RNA G-quadruplexes in their specific regions are significantly related to immune pathways and the COVID-19-related gene sets. Bioinformatics analysis reveals the potential regulatory functions of G-quadruplexes on miRNA targeting at the scale of the whole transcriptome. In addition, the G-quadruplexes are depleted in the putative, not the real, PAS-strong poly(A) sites, which may weaken the possibility of such sites being the real cleaved sites. In brief, our study provides insight into the potential function of RNA G-quadruplexes in post-transcription.

The Evolution of G-quadruplex Structure in mRNA Untranslated Region

The RNA G-quadruplex (rG4) is a kind of non-canonical high-order secondary structure with important biological functions and is enriched in untranslated regions (UTRs) of protein-coding genes. However, how rG4 structures evolve is largely unknown. Here, we systematically investigated the evolution of RNA sequences around UTR rG4 structures in 5 eukaryotic organisms. We found universal selection on UTR sequences, which facilitated rG4 formation in all the organisms that we analyzed. While G-rich sequences were preferred in the rG4 structural region, C-rich sequences were selectively not preferred. The selective pressure acting on rG4 structures in the UTRs of genes with higher G content was significantly smaller. Furthermore, we found that rG4 structures experienced smaller evolutionary selection near the translation initiation region in the 5' UTR, near the polyadenylation signals in the 3' UTR, and in regions flanking the miRNA targets in the 3' UTR. These results suggest universal selection for rG4 formation in the UTRs of eukaryotic genomes and the selection may be related to the biological functions of rG4s.

RNA G-quadruplexes (rG4s): genomics and biological functions

G-quadruplexes (G4s) are non-classical DNA or RNA secondary structures that have been first observed decades ago. Over the years, these four-stranded structural motifs have been demonstrated to have significant regulatory roles in diverse biological processes, but challenges remain in detecting them globally and reliably. Compared to DNA G4s (dG4s), the study of RNA G4s (rG4s) has received less attention until recently. In this review, we will summarize the innovative high-throughput methods recently developed to detect rG4s on a transcriptome-wide scale, highlight the many novel and important functions of rG4 being discovered in vivo across the tree of life, and discuss the key biological questions to be addressed in the near future.

G-Quadruplex Structures in Bacteria: Biological Relevance and Potential as an Antimicrobial Target

DNA strands consisting of multiple runs of guanines can adopt a noncanonical, four-stranded DNA secondary structure known as G-quadruplex or G4 DNA. G4 DNA is thought to play an important role in transcriptional and translational regulation of genes, DNA replication, genome stability, and oncogene expression in eukaryotic genomes. In other organisms, including several bacterial pathogens and some plant species, the biological roles of G4 DNA and G4 RNA are starting to be explored. Recent investigations showed that G4 DNA and G4 RNA are generally conserved across plant species. In silico analyses of several bacterial genomes identified putative guanine-rich, G4 DNA-forming sequences in promoter regions. The sequences were particularly abundant in certain gene classes, suggesting that these highly diverse structures can be employed to regulate the expression of genes involved in secondary metabolite synthesis and signal transduction. Furthermore, in the pathogen Mycobacterium tuberculosis, the distribution of G4 motifs and their potential role in the regulation of gene transcription advocate for the use of G4 ligands to develop novel antitubercular therapies. In this review, we discuss the various roles of G4 structures in bacterial DNA and the application of G4 DNA as inhibitors or therapeutic agents to address bacterial pathogens.

HIV-1 genomic RNA U3 region forms a stable quadruplex-hairpin structure

The U3 promoter region of the HIV-1 long terminal repeat (LTR) has previously been shown to fold into a series of dynamic G-quadruplex structures. Among the G-quadruplexes identified in the LTR sequence, LTR-III was shown to be the most stable in vitro. NMR studies of this 28-nucleotide (nt) DNA revealed a unique quadruplex-hairpin structure. Whether the hairpin forms in RNA element is unknown and the role of the hairpin in the structure and stability of quadruplexes has not been characterized. Here, we used optical and thermodynamic studies to address these questions. The wild-type LTR-III RNA formed a monomolecular quadruplex with a parallel topology using only propeller loops, including the hairpin loop element. By comparison to the WT and variant RNAs, LTR-III DNA structures were more heterogeneous and less stable. Increased stability of the RNA suggests that the RNA quadruplex-hairpin structure may be a more attractive therapeutic target than the analogous DNA element.

CNBP Binds and Unfolds In Vitro G-Quadruplexes Formed in the SARS-CoV-2 Positive and Negative Genome Strands

The Coronavirus Disease 2019 (COVID-19) pandemic has become a global health emergency with no effective medical treatment and with incipient vaccines. It is caused by a new positive-sense RNA virus called severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2). G-quadruplexes (G4s) are nucleic acid secondary structures involved in the control of a variety of biological processes including viral replication. Using several G4 prediction tools, we identified highly putative G4 sequences (PQSs) within the positive-sense (+gRNA) and negative-sense (−gRNA) RNA strands of SARS-CoV-2 conserved in related betacoronaviruses. By using multiple biophysical techniques, we confirmed the formation of two G4s in the +gRNA and provide the first evidence of G4 formation by two PQSs in the −gRNA of SARS-CoV-2. Finally, biophysical and molecular approaches were used to demonstrate for the first time that CNBP, the main human cellular protein bound to SARS-CoV-2 RNA genome, binds and promotes the unfolding of G4s formed by both strands of SARS-CoV-2 RNA genome. Our results suggest that G4s found in SARS-CoV-2 RNA genome and its negative-sense replicative intermediates, as well as the cellular proteins that interact with them, are relevant factors for viral genes expression and replication cycle, and may constitute interesting targets for antiviral drugs development.

Learning to Model G-Quadruplexes: Current Methods and Perspectives

G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.

MD-TSPC4: Computational Method for Predicting the Thermal Stability of I-Motif

I-Motif is a tetrameric cytosine-rich DNA structure with hemi-protonated cytosine: cytosine base pairs. Recent evidence showed that i-motif structures in human cells play regulatory roles in the genome. Therefore, characterization of novel i-motifs and investigation of their functional implication are urgently needed for comprehensive understanding of their roles in gene regulation. However, considering the complications of experimental investigation of i-motifs and the large number of putative i-motifs in the genome, development of an in silico tool for the characterization of i-motifs in the high throughput scale is necessary. We developed a novel computation method, MD-TSPC4, to predict the thermal stability of i-motifs based on molecular modeling and molecular dynamic simulation. By assuming that the flexibility of loops in i-motifs correlated with thermal stability within certain temperature ranges, we evaluated the correlation between the root mean square deviations (RMSDs) of model structures and the thermal stability as the experimentally obtained melting temperature (Tm). Based on this correlation, we propose an equation for Tm prediction from RMSD. We expect this method can be useful for estimating the overall structure and stability of putative i-motifs in the genome, which can be a starting point of further structural and functional studies of i-motifs.