Minimum-free-energy distribution of RNA secondary structures: Entropic and thermodynamic properties of rare events

Probing the large deviations for the beta random walk in random medium

We consider a discrete-time random walk on a one-dimensional lattice with space- and time-dependent random jump probabilities, known as the beta random walk. We are interested in the probability that, for a given realization of the jump probabilities (a sample), a walker starting at the origin at time t=0 is at position beyond ξsqrt[T/2] at time T. This probability fluctuates from sample to sample and we study the large-deviation rate function, which characterizes the tails of its distribution at large time T≫1. It is argued that, up to a simple rescaling, this rate function is identical to the one recently obtained exactly by two of the authors for the continuum version of the model. That continuum model also appears in the macroscopic fluctuation theory of a class of lattice gases, e.g., in the so-called KMP model of heat transfer. An extensive numerical simulation of the beta random walk, based on an importance sampling algorithm, is found in good agreement with the detailed analytical predictions. A first-order transition in the tilted measure, predicted to occur in the continuum model, is also observed in the numerics.

Influence of the mRNA initial region on protein production: a case study using recombinant detoxified pneumolysin as a model

Recombinant proteins are of great importance in modern society, mostly as biopharmaceutical products. However, challenging and complex processes with low production yield are major drawbacks. Normally, the optimization to overcome these obstacles is focused on bioreactor and purification processes, and the biomolecular aspects are neglected, seen as less important. In this work, we present how the 5' mRNA secondary structure region can be relevant for translation and, therefore, protein production. For this, Escherichia coli BL21(DE3) clones, producing recombinant detoxified pneumolysin (PdT) with and without the N-terminal His-tag, were cultivated in 10-L bioreactors. Another version of the pdt gene (version 2) with synonymous changes in the 5'-end nucleotide sequence was also obtained. Protein production, plasmid stability, carbon sources, and acetic acid were quantified during the cultures. Furthermore, in silico mRNA analyses were performed using TIsigner and RNAfold. The results showed that the His-tag presence at the N-terminus generated a minimum 1.5-fold increase in target protein synthesis, which was explained by the in silico mRNA analyses that returned an mRNA secondary structure easier to translate and, therefore, higher protein production than without the His-tag. The pdt gene version 2 showed lower 5' mRNA opening energy than version 1, allowing higher PdT production even without a tag. This work reveals that simple mRNA analyses during heterologous gene design and production steps can help reach high-recombinant protein titers in a shorter time than using only traditional bioprocess optimization strategies.

Signaling landscape of mitochondrial non-coding RNAs

Human mitochondria are the vital cell organelle acting as a storehouse of energy generation and diverse regulatory functions. Mitochondrial DNA comprises 93% coding region and 7% non-coding regions, in which the non-coding region hypothesized as responsible for signaling is our specific interest. Here, we explored the unknown functions of mitochondrial non-coding RNAs by studying their respective signaling pathways. We retrieved conserved motifs of interactions from known experimental protein-RNA complexes to model unknown mitochondrial ncRNA sequences. Our results provide the ncRNAs list and show their involvement in four crucial pathways, such as (i) Processing of Capped Intron-Containing Pre-mRNA, (ii) Spliceosome, (iii) Spliceosomal assembly, and (iv) RNA Polymerase II Transcription, respectively. The interactome analysis revealed that the SRSF2 and U2AF2 proteins interact with ncRNAs. Further, we have simulated the selected ncRNA-protein complexes in cell-like environmental conditions and found them stable in terms of energetics. Through our study, we have identified an apparent interaction of mitochondrial ncRNAs with proteins and their role in critical signaling pathways, providing new insights into mitochondrial ncRNA-dependent gene regulation.Communicated by Ramaswamy H. Sarma.

Investigation into the prevalent CRISPR-Cas systems among the Aeromonas genus

Clustered regularly interspaced short palindromic repeats (CRISPR) is a prokaryotic adaptive immune system that checks invasion by mobile genetic elements through nuclease targeting. In this study, we investigated the occurrence, diversity, and features of the CRISPR system in the genus Aeromonas using bioinformatics tools. Only 13 out of 122 complete genomes (10.66%) of the genus Aeromonas from the NCBI GenBank database harbored the CRISPR system. The Type I-F system was the most prevalent CRISPR system among the Aeromonads, followed by the Type I-E system. Only one strain harbored a Type I-C CRISPR system. Among the Aeromonads, Aeromonas caviae (22.7%) and Aeromonas veronii (20%) had a higher prevalence rate of the complete CRISPR system. The analysis of direct repeat (DR) sequences showed that all could form stable RNA secondary structures. A phylogenetic tree generated for the Cas1 protein classified CRISPR subtypes into three distinct clusters. Among the 748 spacers investigated, 41.98% and 17.25% showed perfect homology to phage and plasmid sequences, respectively. Some arrays had duplicated spacers. The CRISPR loci are closely linked to antibiotic resistance genes in most strains. Collectively, our results would contribute to research on antibiotic resistance in the Aeromonas group, and provide new insights into the diversity and evolution of the CRISPR-Cas system.

Probing large deviations of the Kardar-Parisi-Zhang equation at short times with an importance sampling of directed polymers in random media

The one-point distribution of the height for the continuum Kardar-Parisi-Zhang equation is determined numerically using the mapping to the directed polymer in a random potential at high temperature. Using an importance sampling approach, the distribution is obtained over a large range of values, down to a probability density as small as 10^{-1000} in the tails. The short-time behavior is investigated and compared with recent analytical predictions for the large-deviation forms of the probability of rare fluctuations, showing a spectacular agreement with the analytical expressions. The flat and stationary initial conditions are studied in the full space, together with the droplet initial condition in the half-space.

Analysis of in vitro evolution reveals the underlying distribution of catalytic activity among random sequences

The emergence of catalytic RNA is believed to have been a key event during the origin of life. Understanding how catalytic activity is distributed across random sequences is fundamental to estimating the probability that catalytic sequences would emerge. Here, we analyze the in vitro evolution of triphosphorylating ribozymes and translate their fitnesses into absolute estimates of catalytic activity for hundreds of ribozyme families. The analysis efficiently identified highly active ribozymes and estimated catalytic activity with good accuracy. The evolutionary dynamics follow Fisher's Fundamental Theorem of Natural Selection and a corollary, permitting retrospective inference of the distribution of fitness and activity in the random sequence pool for the first time. The frequency distribution of rate constants appears to be log-normal, with a surprisingly steep dropoff at higher activity, consistent with a mechanism for the emergence of activity as the product of many independent contributions.

Chromatin accessibility and guide sequence secondary structure affect CRISPR-Cas9 gene editing efficiency

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated protein 9 (CRISPR-Cas9) systems have emerged as the method of choice for genome editing, but large variations in on-target efficiencies continue to limit their applicability. Here, we investigate the effect of chromatin accessibility on Cas9-mediated gene editing efficiency for 20 gRNAs targeting 10 genomic loci in HEK293T cells using both SpCas9 and the eSpCas9(1.1) variant. Our study indicates that gene editing is more efficient in euchromatin than in heterochromatin, and we validate this finding in HeLa cells and in human fibroblasts. Furthermore, we investigate the gRNA sequence determinants of CRISPR-Cas9 activity using a surrogate reporter system and find that the efficiency of Cas9-mediated gene editing is dependent on guide sequence secondary structure formation. This knowledge can aid in the further improvement of tools for gRNA design.

Novel Approach to Analyzing MFE of Noncoding RNA Sequences

Genomic studies have become noncoding RNA (ncRNA) centric after the study of different genomes provided enormous information on ncRNA over the past decades. The function of ncRNA is decided by its secondary structure, and across organisms, the secondary structure is more conserved than the sequence itself. In this study, the optimal secondary structure or the minimum free energy (MFE) structure of ncRNA was found based on the thermodynamic nearest neighbor model. MFE of over 2600 ncRNA sequences was analyzed in view of its signal properties. Mathematical models linking MFE to the signal properties were found for each of the four classes of ncRNA analyzed. MFE values computed with the proposed models were in concordance with those obtained with the standard web servers. A total of 95% of the sequences analyzed had deviation of MFE values within ±15% relative to those obtained from standard web servers.

Studying the features of 57 confirmed CRISPR loci in 29 strains of Escherichia coli

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) system is a novel type of innate defense system in prokaryotes for destruction of exogenous elements. To gain further insight into behavior and organization of the system, the extensive analysis of the available sequenced genomes is necessary. The dynamic nature of CRISPR loci is possibly valuable for typing and relative analyses of strains and microbial population. There are a few orderly bioinformatics investigations about the structure of CRISPR sequences in the Escherichia coli strains. In this study, 57 CRISPR loci were selected from 32 Escherichia coli strains to investigate their structural characteristics and potential functions using bioinformatics tools. Our results showed that most strains contained several loci that mainly included conserved direct repeats, while the spacers were highly variable. Moreover, RNA analysis of the sequences indicated that all loci could form stable RNA secondary structures and showed homology mostly with phages compared to plasmids. Only three strains included cas genes around their loci.