Estimating uncertainty in predicted folding free energy changes of RNA secondary structures

Structural Impact of 3-methylcytosine Modification on the Anticodon Stem-loop of a Neuronally-enriched Arginine tRNA

All tRNAs undergo a series of chemical modifications to fold and function correctly. In mammals, the C32 nucleotide in the anticodon loop of tRNA-Arg-CCU and UCU is methylated to form 3-methylcytosine (m3C). Deficiency of m3C in arginine tRNAs has been linked to human neurodevelopmental disorders, indicating a critical biological role for m3C modification. However, the structural repercussions of m3C modification are not well understood. Here, we examine the structural effects of m3C32 modification on the anticodon stem loop (ASL) of human tRNA-Arg-UCU-4-1, a unique tRNA with enriched expression in the central nervous system. Optical melting experiments demonstrate that m3C modification can locally disrupt nearby base pairing within the ASL while simultaneously stabilizing the ASL electrostatically, resulting in little net change thermodynamically. The isoenergetic nature of the C32-A38 pair versus the m3C32-A38 pair may help discriminate against structures not adopting canonical C32-A38 pairings, as most other m3C pairings are unfavorable. Furthermore, multidimensional NMR reveals that after m3C modification there are changes in hairpin loop structure and dynamics, the structure of A37, and the neighboring A31-U39 base pair. However, these structural changes after modification are made while maintaining the shape of the C32-A38 pairing, which is essential for efficient tRNA function in translation. These findings suggest that m3C32 modification could alter interactions of tRNA-Arg isodecoders with one or more binding partners while simultaneously maintaining the tRNA's ability to function in translation.

Conformational dynamics of the hepatitis C virus 3'X RNA

The 3' end of the hepatitis C virus genome is terminated by a highly conserved, 98 nt sequence called 3'X. This untranslated structural element is thought to regulate several essential RNA-dependent processes associated with infection. 3'X has two proposed conformations comprised of either three or two stem-loop structures that result from the different base-pairing interactions within the first 55 nt. Here, we used single-molecule Förster resonance energy transfer spectroscopy to monitor the conformational status of fluorescently labeled constructs that isolate this region of the RNA (3'X55). We observed that 3'X55 can adopt both proposed conformations and the relative abundance of them can be modulated by either solution conditions or nucleotide deletions. Furthermore, interconversion between the two conformations takes place over the course of several hours. The simultaneous existence of two slowly interconverting conformations may help prime individual copies of the viral genome for either viral protein or RNA synthesis, thereby minimizing conflicts between these two competing processes.

Estimating RNA Secondary Structure Folding Free Energy Changes with efn2

A number of analyses require estimates of the folding free energy changes of specific RNA secondary structures. These predictions are often based on a set of nearest neighbor parameters that models the folding stability of a RNA secondary structure as the sum of folding stabilities of the structural elements that comprise the secondary structure. In the software suite RNAstructure, the free energy change calculation is implemented in the program efn2. The efn2 program estimates the folding free energy change and the experimental uncertainty in the folding free energy change. It can be run through the graphical user interface for RNAstructure, from the command line, or a web server. This chapter provides detailed protocols for using efn2.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Active secretion of a thermostable transglutaminase variant in Escherichia coli

Background

Streptomyces mobaraenesis transglutaminase (smTG) is widely used to generate protein crosslinking or attachment of small molecules. However, the low thermostability is a main obstacle for smTG application. In addition, it is still hard to achieve the secretory expression of active smTG in E. coli, which benefits the enzyme evolution. In this study, a combined strategy was conducted to improve the thermostability and secretory expression of active smTG in E. coli.

Results

First, the thermostable S. mobaraenesis transglutaminase variant S2P-S23V-Y24N-S199A-K294L (TGm1) was intracellularly expressed in pro-enzyme form in E. coli. Fusing the pro-region of Streptomyces hygroscopicus transglutaminase (proH) and TrxA achieved a 9.78 U/mL of intracellular smTG activity, 1.37-fold higher than the TGm1 fused with its native pro-region. After in vitro activation by dispase, the TGm1 with proH yielded FRAPD-TGm1, exhibiting 0.95 ℃ and 94.25% increases in melting temperature and half-life at 60 ℃ compared to FRAP-TGm1 derived from the expression using its native pro-region, respectively. Second, the TGm1 with proH was co-expressed with transglutaminase activating protease and chaperones (DnaK, DnaJ, and GrpE) in E. coli, achieving 9.51 U/mL of intracellular FRAPD-TGm1 without in vitro activation. Third, the pelB signal peptide was used to mediate the secretory expression of active TGm in E. coli, yielding 0.54 U/mL of the extracellular FRAPD-TGm1. A script was developed to shuffle the codon of pelB and calculate the corresponding mRNA folding energy. A 1.8-fold increase in the extracellular expression of FRAPD-TGm1 was achieved by the Top-9 pelB sequence derived from the coding sequences with the lowest mRNA folding energy. Last, deleting the gene of Braun’s lipoprotein further increased the extracellular yield of FRAPD-TGm1 by 31.2%, reached 1.99 U/mL.

Conclusions

The stabilized FRAPD-smTG here could benefit the enzyme application in food and non-food sectors, while the E. coli system that enables secretory expression of active smTG will facilitate the directed evolution for further improved catalytic properties. The combined strategy (N-terminal modification, co-expression with chaperones, mRNA folding energy optimization of signal peptide, and lipoprotein deletion) may also improve the secretory expression of other functional proteins in E. coli.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01801-9.

Specific length and structure rather than high thermodynamic stability enable regulatory mRNA stem-loops to pause translation

Translating ribosomes unwind mRNA secondary structures by three basepairs each elongation cycle. Despite the ribosome helicase, certain mRNA stem-loops stimulate programmed ribosomal frameshift by inhibiting translation elongation. Here, using mutagenesis, biochemical and single-molecule experiments, we examine whether high stability of three basepairs, which are unwound by the translating ribosome, is critical for inducing ribosome pauses. We find that encountering frameshift-inducing mRNA stem-loops from the E. coli dnaX mRNA and the gag-pol transcript of Human Immunodeficiency Virus (HIV) hinders A-site tRNA binding and slows down ribosome translocation by 15-20 folds. By contrast, unwinding of first three basepairs adjacent to the mRNA entry channel slows down the translating ribosome by only 2-3 folds. Rather than high thermodynamic stability, specific length and structure enable regulatory mRNA stem-loops to stall translation by forming inhibitory interactions with the ribosome. Our data provide the basis for rationalizing transcriptome-wide studies of translation and searching for novel regulatory mRNA stem-loops.

An in-silico approach to design potential siRNAs against the ORF57 of Kaposi's sarcoma-associated herpesvirus

Kaposi's sarcoma-associated herpesvirus (KSHV) is one of the few human oncogenic viruses, which causes a variety of malignancies, including Kaposi's sarcoma, multicentric Castleman disease, and primary effusion lymphoma, particularly in human immunodeficiency virus patients. The currently available treatment options cannot always prevent the invasion and dissemination of this virus. In recent times, siRNA-based therapeutics are gaining prominence over conventional medications as siRNA can be designed to target almost any gene of interest. The ORF57 is a crucial regulatory protein for lytic gene expression of KSHV. Disruption of this gene translation will inevitably inhibit the replication of the virus in the host cell. Therefore, the ORF57 of KSHV could be a potential target for designing siRNA-based therapeutics. Considering both sequence preferences and target site accessibility, several online tools (i-SCORE Designer, Sfold web server) had been utilized to predict the siRNA guide strand against the ORF57. Subsequently, off-target filtration (BLAST), conservancy test (fuzznuc), and thermodynamics analysis (RNAcofold, RNAalifold, and RNA Structure web server) were also performed to select the most suitable siRNA sequences. Finally, two siRNAs were identified that passed all of the filtration phases and fulfilled the thermodynamic criteria. We hope that the siRNAs predicted in this study would be helpful for the development of new effective therapeutics against KSHV.