Length Of RNA Research Articles

Engineering a DNA polymerase for modifying large RNA at specific positions.

The synthesis of large RNA with precise modifications at specific positions is in high demand for both basic research and therapeutic applications, but efficient methods are limited. Engineered DNA polymerases have recently emerged as attractive tools for RNA labelling, offering distinct advantages over conventional RNA polymerases. Here, through semi-rational designs, we engineered a DNA polymerase variant and used it to precisely incorporate a diverse range of modifications, including base modifications, 2'-ribose modifications and backbone modifications, into desired positions within RNA. We achieved efficiencies exceeding 85% in the majority of modification cases, demonstrating success in introducing 2'-O-methyl, phosphorothioate, N4-acetylcytidine and a fluorophore to specific sites in eGFP and Firefly luciferase messenger RNA. Our mRNA products with N4-acetylcytidine, 2'-O-methyl and/or phosphorothioate have demonstrated the ability to enhance stability and affect protein production. This method presents a promising tool for the comprehensive functionalization of RNA, enabling the introduction of plentiful modifications irrespective of RNA lengths and sequences.

CParty: hierarchically constrained partition function of RNA pseudoknots.

Biologically relevant RNA secondary structures are routinely predicted by efficient dynamic programming algorithms that minimize their free energy. Starting from such algorithms, one can devise partition function algorithms, which enable stochastic perspectives on RNA structure ensembles. As the most prominent example, McCaskill's partition function algorithm is derived from pseudoknot-free energy minimization. While this algorithm became hugely successful for the analysis of pseudoknot-free RNA structure ensembles, as of yet there exists only one pseudoknotted partition function implementation, which covers only simple pseudoknots and comes with a borderline-prohibitive complexity of O(n5) in the RNA length n. Here, we develop a partition function algorithm corresponding to the hierarchical pseudoknot prediction of HFold, which performs exact optimization in a realistic pseudoknot energy model. In consequence, our algorithm CParty carries over HFold's advantages over classical pseudoknot prediction in characterizing the Boltzmann ensemble at equilibrium. Given an RNA sequence S and a pseudoknot-free structure G, CParty computes the partition function over all possibly pseudoknotted density-2 structures G∪G' of S that extend the fixed G by a disjoint pseudoknot-free structure G'. Thus, CParty follows the common hypothesis of hierarchical pseudoknot formation, where pseudoknots form as tertiary contacts only after a first pseudoknot-free "core" G and we call the computed partition function hierarchically constrained (by G). Like HFold, the dynamic programming algorithm CParty is very efficient, achieving the low complexity of the pseudoknot-free algorithm, i.e. cubic time and quadratic space. Finally, by computing pseudoknotted ensemble energies, we unveil kinetics features of a therapeutic target in SARS-CoV-2. CParty is available at https://github.com/HosnaJabbari/CParty.

Full-length direct RNA sequencing uncovers stress granule-dependent RNA decay upon cellular stress.

Cells react to stress by triggering response pathways, leading to extensive alterations in the transcriptome to restore cellular homeostasis. The role of RNA metabolism in shaping the cellular response to stress is vital, yet the global changes in RNA stability under these conditions remain unclear. In this work, we employ direct RNA sequencing with nanopores, enhanced by 5' end adapter ligation, to comprehensively interrogate the human transcriptome at single-molecule and -nucleotide resolution. By developing a statistical framework to identify robust RNA length variations in nanopore data, we find that cellular stress induces prevalent 5' end RNA decay that is coupled to translation and ribosome occupancy. Unlike typical RNA decay models in normal conditions, we show that stress-induced RNA decay is dependent on XRN1 but does not depend on deadenylation or decapping. We observed that RNAs undergoing decay are predominantly enriched in the stress granule transcriptome while inhibition of stress granule formation via genetic ablation of G3BP1 and G3BP2 rescues RNA length. Our findings reveal RNA decay as a key component of RNA metabolism upon cellular stress that is dependent on stress granule formation.

Silencing ACE1 Gene with dsRNA of Different Lengths Impairs Larval Development in Leptinotarsa decemlineata.

In the search for effective strategies to control the Colorado Potato Beetle, RNA interference technology has emerged as a promising method due to its capacity to suppress genes selectively. Factors such as the target gene and double-stranded RNA (dsRNA) length are critical for optimizing gene silencing efficiency. In this study, we designed and synthesized in vitro dsRNAs of varying lengths targeting the ACE1 gene, which encodes the AChE1 isoform of acetylcholinesterase in the beetle. All tested dsRNA lengths (222 bp, 543 bp, 670 bp, and 870 bp) promoted transcript reduction. The 670 bp dsRNA was the most effective, reducing transcript levels by approximately 40% by day seven, followed by the 543 bp dsRNA. No significant differences were observed between the 222 bp and 870 bp dsRNAs. Furthermore, all of the dsRNA lengths resulted in reduced weight gain and increased mortality in larvae, with the 670 bp dsRNA showing the highest mortality rate, leaving only 63% larval survival, a trend that persisted through day nine. These findings emphasize that dsRNA length is a key factor in the silencing response, underscoring the importance of selecting the optimal length while considering the gene's target, stability, and delivery methods. This study contributes to establishing design criteria for dsRNA, aiding in the development of more effective and sustainable pest management strategies.

Full-length direct RNA sequencing uncovers stress-granule dependent RNA decay upon cellular stress.

Cells react to stress by triggering response pathways, leading to extensive alterations in the transcriptome to restore cellular homeostasis. The role of RNA metabolism in shaping the cellular response to stress is vital, yet the global changes in RNA stability under these conditions remain unclear. In this work, we employ direct RNA sequencing with nanopores, enhanced by 5' end adaptor ligation, to comprehensively interrogate the human transcriptome at single-molecule and nucleotide resolution. By developing a statistical framework to identify robust RNA length variations in nanopore data, we find that cellular stress induces prevalent 5' end RNA decay that is coupled to translation and ribosome occupancy. Unlike typical RNA decay models in normal conditions, we show that stress-induced RNA decay is dependent on XRN1 but does not depend on deadenylation or decapping. We observed that RNAs undergoing decay are predominantly enriched in the stress granule transcriptome while inhibition of stress granule formation via genetic ablation of G3BP1 and G3BP2 rescues RNA length. Our findings reveal RNA decay as a key determinant of RNA metabolism upon cellular stress and dependent on stress-granule formation.

Did the exposure of coacervate droplets to rain make them the first stable protocells?

Membraneless coacervate microdroplets have long been proposed as model protocells as they can grow, divide, and concentrate RNA by natural partitioning. However, the rapid exchange of RNA between these compartments, along with their rapid fusion, both within minutes, means that individual droplets would be unable to maintain their separate genetic identities. Hence, Darwinian evolution would not be possible, and the population would be vulnerable to collapse due to the rapid spread of parasitic RNAs. In this study, we show that distilled water, mimicking rain/freshwater, leads to the formation of electrostatic crosslinks on the interface of coacervate droplets that not only suppress droplet fusion indefinitely but also allow the spatiotemporal compartmentalization of RNA on a timescale of days depending on the length and structure of RNA. We suggest that these nonfusing membraneless droplets could potentially act as protocells with the capacity to evolve compartmentalized ribozymes in prebiotic environments.

In vitro higher-order oligomeric assembly of the respiratory syncytial virus M2-1 protein with longer RNAs.

The respiratory syncytial virus (RSV) M2-1 protein is a transcriptional antitermination factor crucial for efficiently synthesizing multiple full-length viral mRNAs. During RSV infection, M2-1 exists in a complex with mRNA within cytoplasmic compartments called inclusion body-associated granules (IBAGs). Prior studies showed that M2-1 can bind along the entire length of viral mRNAs instead of just gene-end (GE) sequences, suggesting that M2-1 has more sophisticated RNA recognition and binding characteristics. Here, we analyzed the higher oligomeric complexes formed by M2-1 and RNAs in vitro using size exclusion chromatography (SEC), electrophoretic mobility shift assays (EMSA), negative stain electron microscopy (EM), and mutagenesis. We observed that the minimal RNA length for such higher oligomeric assembly is about 14 nucleotides for polyadenine sequences, and longer RNAs exhibit distinct RNA-induced binding modality to M2-1, leading to enhanced particle formation frequency and particle homogeneity as the local RNA concentration increases. We showed that particular cysteine residues of the M2-1 cysteine-cysteine-cystine-histidine (CCCH) zinc-binding motif are essential for higher oligomeric assembly. Furthermore, complexes assembled with long polyadenine sequences remain unaffected when co-incubated with ribonucleases or a zinc chelation agent. Our study provided new insights into the higher oligomeric assembly of M2-1 with longer RNA.IMPORTANCERespiratory syncytial virus (RSV) causes significant respiratory infections in infants, the elderly, and immunocompromised individuals. The virus forms specialized compartments to produce genetic material, with the M2-1 protein playing a pivotal role. M2-1 acts as an anti-terminator in viral transcription, ensuring the creation of complete viral mRNA and associating with both viral and cellular mRNA. Our research focuses on understanding M2-1's function in viral mRNA synthesis by modeling interactions in a controlled environment. This approach is crucial due to the challenges of studying these compartments in vivo. Reconstructing the system in vitro uncovers structural and biochemical aspects and reveals the potential functions of M2-1 and its homologs in related viruses. Our work may contribute to identifying targets for antiviral inhibitors and advancing RSV infection treatment.

RiboDiffusion: tertiary structure-based RNA inverse folding with generative diffusion models.

RNA design shows growing applications in synthetic biology and therapeutics, driven by the crucial role of RNA in various biological processes. A fundamental challenge is to find functional RNA sequences that satisfy given structural constraints, known as the inverse folding problem. Computational approaches have emerged to address this problem based on secondary structures. However, designing RNA sequences directly from 3D structures is still challenging, due to the scarcity of data, the nonunique structure-sequence mapping, and the flexibility of RNA conformation. In this study, we propose RiboDiffusion, a generative diffusion model for RNA inverse folding that can learn the conditional distribution of RNA sequences given 3D backbone structures. Our model consists of a graph neural network-based structure module and a Transformer-based sequence module, which iteratively transforms random sequences into desired sequences. By tuning the sampling weight, our model allows for a trade-off between sequence recovery and diversity to explore more candidates. We split test sets based on RNA clustering with different cut-offs for sequence or structure similarity. Our model outperforms baselines in sequence recovery, with an average relative improvement of 11% for sequence similarity splits and 16% for structure similarity splits. Moreover, RiboDiffusion performs consistently well across various RNA length categories and RNA types. We also apply in silico folding to validate whether the generated sequences can fold into the given 3D RNA backbones. Our method could be a powerful tool for RNA design that explores the vast sequence space and finds novel solutions to 3D structural constraints. The source code is available at https://github.com/ml4bio/RiboDiffusion.

In Vivo Delivery of Spherical and Cylindrical In Vitro Reconstituted Virus-like Particles Containing the Same Self-Amplifying mRNA.

The dramatic effectiveness of recent mRNA (mRNA)-based COVID vaccines delivered in lipid nanoparticles has highlighted the promise of mRNA therapeutics in general. In this report, we extend our earlier work on self-amplifying mRNAs delivered in spherical in vitro reconstituted virus-like particles (VLPs), and on drug delivery using cylindrical virus particles. In particular, we carry out separate in vitro assemblies of a self-amplifying mRNA gene in two different virus-like particles: one spherical, formed with the capsid protein of cowpea chlorotic mottle virus (CCMV), and the other cylindrical, formed from the capsid protein of tobacco mosaic virus (TMV). The mRNA gene is rendered self-amplifying by genetically fusing it to the RNA-dependent RNA polymerase (RdRp) of Nodamura virus, and the relative efficacies of cell uptake and downstream protein expression resulting from their CCMV- and TMV-packaged forms are compared directly. This comparison is carried out by their transfections into cells in culture: expressions of two self-amplifying genes, enhanced yellow fluorescent protein (EYFP) and Renilla luciferase (Luc), packaged alternately in CCMV and TMV VLPs, are quantified by fluorescence and chemiluminescence levels, respectively, and relative numbers of the delivered mRNAs are measured by quantitative real-time PCR. The cellular uptake of both forms of these VLPs is further confirmed by confocal microscopy of transfected cells. Finally, VLP-mediated delivery of the self-amplifying-mRNA in mice following footpad injection is shown by in vivo fluorescence imaging to result in robust expression of EYFP in the draining lymph nodes, suggesting the potential of these plant virus-like particles as a promising mRNA gene and vaccine delivery modality. These results establish that both CCMV and TMV VLPs can deliver their in vitro packaged mRNA genes to immune cells and that their self-amplifying forms significantly enhance in situ expression. Choice of one VLP (CCMV or TMV) over the other will depend on which geometry of nucleocapsid is self-assembled more efficiently for a given length and sequence of RNA, and suggests that these plant VLP gene delivery systems will prove useful in a wide variety of medical applications, both preventive and therapeutic.

Ribosomal RNA expansion segments and their role in ribosome biology.

Ribosomes are universally conserved cellular machines that catalyze protein biosynthesis. The active sites underly immense evolutionary conservation resulting in virtually identical core structures of ribosomes in all domains of life including organellar ribosomes. However, more peripheral structures of cytosolic ribosomes changed during evolution accommodating new functions and regulatory options. The expansion occurred at the riboprotein level, including more and larger ribosomal proteins and at the RNA level increasing the length of ribosomal RNA. Expansions within the ribosomal RNA occur as clusters at conserved sites that face toward the periphery of the cytosolic ribosome. Recent biochemical and structural work has shed light on how rRNA-specific expansion segments (ESs) recruit factors during translation and how they modulate translation dynamics in the cytosol. Here we focus on recent work on yeast, human and trypanosomal cytosolic ribosomes that explores the role of two specific rRNA ESs within the small and large subunit respectively. While no single regulatory strategy exists, the absence of ESs has consequences for proteomic stability and cellular fitness, rendering them fascinating evolutionary tools for tailored protein biosynthesis.

Repurposing CRISPR-Cas13 systems for robust mRNA trans-splicing

Type VI CRISPR enzymes have been developed as programmable RNA-guided Cas proteins for eukaryotic RNA editing. Notably, Cas13 has been utilized for site-targeted single base edits, demethylation, RNA cleavage or knockdown and alternative splicing. However, the ability to edit large stretches of mRNA transcripts remains a significant challenge. Here, we demonstrate that CRISPR-Cas13 systems can be repurposed to assist trans-splicing of exogenous RNA fragments into an endogenous pre-mRNA transcript, a method termed CRISPR Assisted mRNA Fragment Trans-splicing (CRAFT). Using split reporter-based assays, we evaluate orthogonal Cas13 systems, optimize guide RNA length and screen for optimal trans-splicing site(s) across a range of intronic targets. We achieve markedly improved editing of large 5’ and 3’ segments in different endogenous mRNAs across various mammalian cell types compared to other spliceosome-mediated trans-splicing methods. CRAFT can serve as a versatile platform for attachment of protein tags, studying the impact of multiple mutations/single nucleotide polymorphisms, modification of untranslated regions (UTRs) or replacing large segments of mRNA transcripts.

Multiple Oligo assisted RNA Pulldown via Hybridization followed by Mass Spectrometry (MORPH-MS) for exploring the RNA-Protein interactions

ABSTRACT Understanding RNA-protein interactions is crucial for deciphering the cellular functions and molecular mechanisms of regulatory RNAs. Consequently, there is a constant need to develop innovative and cost-effective methods to uncover such interactions. We developed a simple and cost-effective technique called Multiple Oligo assisted RNA Pulldown via Hybridization (MORPH) to identify proteins interacting with a specific RNA. MORPH employs a tiling array of antisense oligos (ASOs) to efficiently capture the RNA of interest along with proteins associated with it. Unlike existing techniques that rely on multiple individually biotinylated oligos spanning the entire RNA length, MORPH stands out by utilizing a single biotinylated oligo to capture all the ASOs. To evaluate MORPH’s efficacy, we applied this technique combined with mass spectrometry to identify proteins interacting with lncRNA NEAT1, which has previously been studied using various methods. Our results demonstrate that despite being a simple and inexpensive procedure, MORPH performs on par with existing methods. Abbreviations: ASO, Antisense oligo; lncRNA, long non-coding RNA; MORPH, Multiple Oligo assisted RNA Pulldown via Hybridization

Factors Affecting Stability of RNA – Temperature, Length, Concentration, pH, and Buffering Species

RNA is prone to both chemical degradation and/or physical instability. Some of the factors affecting stability of RNA in solution are its length, 3′ poly A tail and 5′ cap integrity, excipients, buffering species, pH of the solution, nucleases, and divalent cations. In this work, we showed the effect of temperature, messenger RNA (mRNA) length, buffering species, pH of the solution, and the concentration of mRNA on its chemical and physical stability. Our thermodynamic analysis of a 4000 nucleotide-long mRNA measured an activation energy of 31.5 kcal/mol normalized per phosphodiester backbone. We found mRNA length to be negatively correlated to its stability. Buffering species and pH of the solution affected mRNA integrity along with affecting the onset temperature of melting obtained by Differential Scanning Calorimetry (DSC) thermograms. It was also found that increasing the concentration of mRNA in solution increased its stability.

Profiling stress-triggered RNA condensation with photocatalytic proximity labeling

Stress granules (SGs) are highly dynamic cytoplasmic membrane-less organelles that assemble when cells are challenged by stress. RNA molecules are sorted into SGs where they play important roles in maintaining the structural stability of SGs and regulating gene expression. Herein, we apply a proximity-dependent RNA labeling method, CAP-seq, to comprehensively investigate the content of SG-proximal transcriptome in live mammalian cells. CAP-seq captures 457 and 822 RNAs in arsenite- and sorbitol-induced SGs in HEK293T cells, respectively, revealing that SG enrichment is positively correlated with RNA length and AU content, but negatively correlated with translation efficiency. The high spatial specificity of CAP-seq dataset is validated by single-molecule FISH imaging. We further apply CAP-seq to map dynamic changes in SG-proximal transcriptome along the time course of granule assembly and disassembly processes. Our data portray a model of AU-rich and translationally repressed SG nanostructure that are memorized long after the removal of stress.

Comparison of Proline-glutamate-proline-glutamate-polymorphic GC-rich Sequences Family Protein Wag22 (Rv1759c), PE_PGRS31 (Rv1768), PE_PGRS32 (Rv1803), and PE_PGRS33 gene (Rv1818c) in Exponential State and Under In vitro Model of Latency in Same Clinical Isolates of Mycobacterium tuberculosis: Frameshift Mutation in Extensively Drug-resistant and Totally Drug-resistant tuberculosis Bacilli

AbstractBackground:Proline-glutamate (PE)/proline-PE (PPE) proteins play an important role in the development of mycobacterial pathogenicity by modulating the host immune system. In the present investigation, the structural changes in PE-polymorphic GC-rich sequences (PGRS) family protein Wag22 (Rv1759c), PE_PGRS31 (Rv1768), PE_PGRS32 (Rv1803), and PE_PGRS33 gene (Rv1818c) were compared and analyzed in exponential state and underin vitromodel of latency in same clinical isolates ofMycobacterium tuberculosis(MTB).Methods:MTB strains were isolated from clinically and laboratory-confirmed cases of tuberculosis (TB). The TB isolates were subjected to the Xpert MTB/rifampin test and then, further susceptibility testing using proportional methods was performed on them. The isolates were characterized using both 16S–23S RNA andhsp65 genes spacer polymerase chain reaction-restriction fragment length polymorphism. Selected isolates studied at two experimental set–up at exponential phase OD 600 = 0.05 (5 cfu/mL × 106 cfu/mL) and under zero oxygen and nutrition for 26 months to selected isolates studied at two experimental setup in exponential phase OD600 = 0.05 (5 cfu/mL × 106 cfu/mL) and under zero oxygen and nutrition after 26 months. Whole-genome sequencing was performed on studied isolates and the protein structures were analyzed using a bioinformatics web server.Results:No deletion, insertion, or substation occurred in susceptible, mono-drug and multidrug resistant-TB isolates were observed at PE-PGRS family protein Wag22 (Rv1759c) and PE_PGRS31 (Rv1768) at exponential phase. Although, a large deletion (at Rv1759c; Rv1768) was observed in extensively drug-resistant (XDR) and totally drug-resistant (TDR) TB isolates at the exponential phase. All studied TDR-TB isolates had a common deletion position from amino acid 1 (methionine) to amino acid 83 (glycine) and from amino acid 725 (proline) to amino acid 914 (threonine) at PE-PGRS family protein Wag22 (Rv1759c). At PE_PGRS32 (Rv1803), deletion occurred from amino acid 1 (methionine) to amino acid 212 (glycine) in latent TDR-TB bacilli. No changes in Rv1803 were observed in other studied isolates. In contrast, 66.6% of studied isolates had either insertion, deletion, substitution, or combination of changes at PE_PGRS33 (Rv1818c). However, the majority of changes at Rv1818c occurred in drug-resistant isolates. We also documented the region of deletion and insertion at PE_PGRS33 (Rv1818c) is different in active and latent TDR-TB isolates.Conclusions:Changes in these PE-PGRS family protein was associated with drug susceptibility patterns of individual isolates. Our result showed a total frameshift mutation of protein that had a different length in comparison to the original protein. These changes might disturb the interactions between XDR and TDR-TB isolates and immune responses, which needs further investigation.

ATP-induced cross-linking of a biomolecular condensate

DEAD-box helicases are important regulators of biomolecular condensates. However, the mechanisms through which these enzymes affect the dynamics of biomolecular condensates have not been systematically explored. Here, we demonstrate the mechanism by which the mutation of a DEAD-box helicase’s catalytic core alters ribonucleoprotein condensate dynamics in the presence of ATP. Through altering RNA length within the system, we are able to attribute the altered biomolecular dynamics and material properties to physical cross-linking of RNA facilitated by the mutant helicase. These results suggest that mutant condensates approach a gel transition when RNA length is increased to lengths comparable to eukaryotic mRNA. Lastly, we show that this cross-linking effect is tunable with ATP concentration, uncovering a system whose RNA mobility and material properties vary with enzyme activity. More generally, these findings point to a fundamental mechanism for modulating condensate dynamics and emergent material properties through nonequilibrium, molecular-scale interactions.

The effect of polymer length in liquid-liquid phase separation

Understanding the thermodynamics that drive liquid-liquid phase separation (LLPS) is quite important given the number of diverse biomolecular systems undergoing this phenomenon. Many studies have focused on condensates of long polymers, but very few systems of short-polymer condensates have been observed and studied. Here, we study a short-polymer system of various lengths of poly-adenine RNA and peptides formed by the RGRGG sequence repeats to understand the underlying thermodynamics of LLPS. Using the recently developed COCOMO coarse-grained (CG) model, we predicted condensates for lengths as short as 5-10 residues, which was then confirmed by experiment, making this one of the smallest LLPS systems yet observed. A free-energy model reveals that the length dependence of condensation is driven primarily by entropy of confinement. The simplicity of this system will provide the basis for understanding more biologically realistic systems.

Surfactants or scaffolds? RNAs of varying lengths control the thermodynamic stability of condensates differently

Biomolecular condensates, thought to form via liquid-liquid phase separation of intracellular mixtures, are multicomponent systems that can include diverse types of proteins and RNAs. RNA is a critical modulator of RNA–protein condensate stability, as it induces an RNA concentration-dependent reentrant phase transition—increasing stability at low RNA concentrations and decreasing it at high concentrations. Beyond concentration, RNAs inside condensates can be heterogeneous in length, sequence, and structure. Here, we use multiscale simulations to understand how different RNA parameters interact with one another to modulate the properties of RNA–protein condensates. To do so, we perform residue/nucleotide resolution coarse-grained molecular dynamics simulations of multicomponent RNA–protein condensates containing RNAs of different lengths and concentrations, and either FUS or PR25 proteins. Our simulations reveal that RNA length regulates the reentrant phase behavior of RNA–protein condensates: increasing RNA length sensitively rises the maximum value that the critical temperature of the mixture reaches, and the maximum concentration of RNA that the condensate can incorporate before beginning to become unstable. Strikingly, RNAs of different lengths are organized heterogeneously inside condensates, which allows them to enhance condensate stability via two distinct mechanisms: shorter RNA chains accumulate at the condensate’s surface acting as natural biomolecular surfactants, while longer RNA chains concentrate inside the core to saturate their bonds and enhance the density of molecular connections in the condensate. Using a patchy particle model, we additionally demonstrate that the combined impact of RNA length and concentration on condensate properties is dictated by the valency, binding affinity, and polymer length of the various biomolecules involved. Our results postulate that diversity on RNA parameters within condensates allows RNAs to increase condensate stability by fulfilling two different criteria: maximizing enthalpic gain and minimizing interfacial free energy; hence, RNA diversity should be considered when assessing the impact of RNA on biomolecular condensates regulation.

Random and Natural Non-Coding RNA Have Similar Structural Motif Patterns but Differ in Bulge, Loop, and Bond Counts.

An important question in evolutionary biology is whether (and in what ways) genotype-phenotype (GP) map biases can influence evolutionary trajectories. Untangling the relative roles of natural selection and biases (and other factors) in shaping phenotypes can be difficult. Because the RNA secondary structure (SS) can be analyzed in detail mathematically and computationally, is biologically relevant, and a wealth of bioinformatic data are available, it offers a good model system for studying the role of bias. For quite short RNA (length L≤126), it has recently been shown that natural and random RNA types are structurally very similar, suggesting that bias strongly constrains evolutionary dynamics. Here, we extend these results with emphasis on much larger RNA with lengths up to 3000 nucleotides. By examining both abstract shapes and structural motif frequencies (i.e., the number of helices, bonds, bulges, junctions, and loops), we find that large natural and random structures are also very similar, especially when contrasted to typical structures sampled from the spaces of all possible RNA structures. Our motif frequency study yields another result, where the frequencies of different motifs can be used in machine learning algorithms to classify random and natural RNA with high accuracy, especially for longer RNA (e.g., ROC AUC 0.86 for L = 1000). The most important motifs for classification are the number of bulges, loops, and bonds. This finding may be useful in using SS to detect candidates for functional RNA within 'junk' DNA regions.

Optimal Conditions for In Vitro Assembly of Respiratory Syncytial Virus Nucleocapsid-like Particles.

The nucleocapsids (NCs) of the respiratory syncytial virus (RSV) can display multiple morphologies in vivo, including spherical, asymmetric, and filamentous conformations. Obtaining homogeneous ring-like oligomers in vitro is significant since they structurally represent one turn of the characteristic RSV NC helical filament. Here, we analyzed and optimized conditions for forming homogenous, recombinant nucleocapsid-like particles (NCLPs) of RSV in vitro. We examined the effects of modifying the integrated RNA length and sequence, altering incubation time, and varying buffer parameters, including salt concentration and pH, on ring-like NCLPs assembly using negative stain electron microscopy (EM) imaging. We showed that high-quality, homogeneous particles are assembled when incubating short, adenine-rich RNA sequences with RNA-free N associated with P (N0P). Further, we reported that a co-incubation duration greater than 3 days, a NaCl concentration between 100 mM and 200 mM, and a pH between 7 and 8 are optimal for N-RNA ring assembly with polyadenine RNA sequences. We believe assembling high-quality, homogeneous NCLPs in vitro will allow for further analysis of RSV RNA synthesis. This work may also lend insights into obtaining high-resolution nucleocapsid homogeneous structures for in vitro analysis of antiviral drug candidates against RSV and related viruses.