BioSim Lab @ UB

Most amino acids are encoded by multiple codons, making the genetic code degenerate. Synonymous mutations affect protein translation and folding, but their impact on RNA itself is often neglected. We developed a genetic algorithm that introduces synonymous mutations to control the diversity of structures sampled by an mRNA. The behavior of the designed mRNAs reveals a physical code layered in the genetic code. We find that mRNA conformational heterogeneity directs physical properties and functional outputs of RNA-protein complexes and biomolecular condensates. The role of structure and disorder of proteins in biomolecular condensates is well appreciated, but we find that RNA conformational heterogeneity is equally important. This feature of RNA enables both evolution and engineers to build cellular structures with specific material and responsive properties.

Drosophila germ granules enrich mRNAs critical for fly development. Within germ granules, mRNAs form multi-transcript clusters marked by increased mRNA concentration, creating an elevated potential for intermolecular base pairing. However, the type and abundance of intermolecular base pairing in mRNA clusters is poorly characterized. Using single-molecule super-resolution microscopy, chemical probing for base accessibility, phase separation assays, and simulations, we demonstrated that mRNAs remain well-folded upon localization to germ granules. While most base pairing is intramolecular, mRNAs still display the ability for intermolecular base pairing, facilitating clustering without high sequence complementarity or significant melting of secondary structure. This base pairing among mRNAs is driven by scattered and discontinuous stretches of bases appearing on the surface of folded RNAs, providing multivalency to clustering but exhibits low probability for sustained interactions. Notably, engineered germ granule mRNAs with exposed GC-rich complementary sequences (CSs) presented within stable stem loops induce sustained base pairing in vitro and enhanced intermolecular interactions in vivo. However, the presence of these stem loops alone disrupts fly development, and the addition of GC-rich CSs exacerbates this phenotype. Although germ granule mRNAs contain numerous GC-rich CSs capable of stable intermolecular base pairing, they are primarily embedded by RNA folding. This study emphasizes the role of RNA folding in controlling the type and abundance of intermolecular base pairing, thereby preserving the functional integrity of mRNAs within the germ granules.

Repeat RNA sequences self-associate to form condensates. Simulations of a coarse-grained single-interaction site model for (CAG)n (n = 30 and 31) show that the salt-dependent free energy gap, ΔGS, between the ground (perfect hairpin) and the excited state (slipped hairpin (SH) with one CAG overhang) of the monomer for (n even) is the primary factor that determines the rates and yield of self-assembly. For odd n, the free energy (GS) of the ground state, which is an SH, is used to predict the self-association kinetics. As the monovalent salt concentration, CS, increases, ΔGS and GS increase, which decreases the rates of dimer formation. In contrast, ΔGS for shuffled sequences, with the same length and sequence composition as (CAG)31, is larger, which suppresses their propensities to aggregate. Although demonstrated explicitly for (CAG) polymers, the finding of inverse correlation between the free energy gap and RNA aggregation is general.

Low-complexity nucleotide repeat sequences, which are implicated in several neurological disorders, undergo liquid–liquid phase separation (LLPS) provided the number of repeat units, n, exceeds a critical value. Here, we establish a link between the folding landscapes of the monomers of trinucleotide repeats and their propensity to self-associate. Simulations using a coarse-grained Self-Organized Polymer (SOP) model for (CAG)n repeats in monovalent salt solutions reproduce experimentally measured melting temperatures, which are available only for small n. By extending the simulations to large n, we show that the free-energy gap, ΔGS, between the ground state (GS) and slipped hairpin (SH) states is a predictor of aggregation propensity. The GS for even n is a perfect hairpin (PH), whereas it is a SH when n is odd. The value of ΔGS (zero for odd n) is larger for even n than for odd n. As a result, the rate of dimer formation is slower in (CAG)30 relative to (CAG)31, thus linking ΔGS to RNA–RNA association. The yield of the dimer decreases dramatically, compared to the wild type, in mutant sequences in which the population of the SH decreases substantially. Association between RNA chains is preceded by a transition to the SH even if the GS is a PH. The finding that the excitation spectrum—which depends on the exact sequence, n, and ionic conditions—is a predictor of self-association should also hold for other RNAs (mRNA for example) that undergo LLPS.

Although it is known that RNA undergoes liquid–liquid phase separation, the interplay between the molecular driving forces and the emergent features of the condensates, such as their morphologies and dynamic properties, is not well understood. We introduce a coarse-grained model to simulate phase separation of trinucleotide repeat RNAs, which are implicated in neurological disorders. After establishing that the simulations reproduce key experimental findings, we show that once recruited inside the liquid droplets, the monomers transition from hairpin-like structures to extended states. Interactions between the monomers in the condensates result in the formation of an intricate and dense intermolecular network, which severely restrains the fluctuations and mobilities of the RNAs inside large droplets. In the largest densely packed high-viscosity droplets, the mobility of RNA chains is best characterized by reptation, reminiscent of the dynamics in polymer melts. Our work provides a microscopic framework for understanding liquid–liquid phase separation in RNA, which is not easily discernible in current experiments.