Cited by 409RNA and protein synthesis mechanismsRNA Interference and Gene DeliveryRNA Research and Splicing

The DOI matches its Crossref record; 115 of the 117 references it cites were found in a registry or an index; 22 of the 29 sources its claims cite were checked for support.

Combinatorial optimization of mRNA structure, stability, and translation for RNA-based therapeutics

Open access

Published version of a preprint

Abstract

Therapeutic mRNAs and vaccines are being developed for a broad range of human diseases, including COVID-19. However, their optimization is hindered by mRNA instability and inefficient protein expression. Here, we describe design principles that overcome these barriers. We develop a new RNA sequencing-based platform called PERSIST-seq to systematically delineate in-cell mRNA stability, ribosome load, as well as in-solution stability of a library of diverse mRNAs. We find that, surprisingly, in-cell stability is a greater driver of protein output than high ribosome load. We further introduce a method called In-line-seq, applied to thousands of diverse RNAs, that reveals sequence and structure-based rules for mitigating hydrolytic degradation. Our findings show that "superfolder" mRNAs can be designed to improve both stability and expression that are further enhanced through pseudouridine nucleoside modification. Together, our study demonstrates simultaneous improvement of mRNA stability and protein expression and provides a computational-experimental platform for the enhancement of mRNA medicines.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper develops PERSIST-seq and In-line-seq to systematically optimize therapeutic mRNA stability, structure, and translation. It finds that in-cell mRNA stability is a greater driver of total protein output than high ribosome load, and that structure-optimized CDSs plus pseudouridine modification can improve both in-solution stability and protein expression.

Key findings

  • In a 233-mRNA library, PERSIST-seq showed in-cell half-lives ranging from <5 h to >15 h; a model integrating ribosome load and in-cell stability predicted measured luciferase activity with r = 0.71, supporting in-cell stability as the dominant predictor of total protein output.
  • In-solution stability was driven mainly by mRNA length and structure: the shortest MEV-encoding mRNAs had half-lives of 3.4 ± 0.6 h, while the longest eGFP-encoding mRNAs had half-lives of 1.1 ± 0.08 h; the ratio 3.0 ± 0.7 matched the inverse length ratio 958 nt/250 nt = 3.8.
  • In-line-seq on 3,030 RNA fragments (2,165 after filtering) yielded DegScore, which predicted Nluc CDS in-solution half-lives with Pearson R = −0.63, p < 0.0001, outperforming dGMFE (R = −0.53) and summed unpaired probability (R = −0.60).
  • Pseudouridine (ψ) or m1ψ substitution suppressed in-line hydrolysis and gave 1.2–2.7-fold stabilization for longer Nluc mRNAs; a RiboTree/LinearDesign/DegScore-optimized mRNA had t1/2 = 2.4 relative to Nluc start and 1.7-fold higher protein expression, and after 2 h accelerated degradation most ψ-modified optimized mRNA remained functional while the starting mRNA gave negligible activity.

Methods

  • PERSIST-seq: pooled full-length IVT mRNA library of 233 constructs with varied 5′ UTRs, 3′ UTRs, and CDSs, barcoded and transfected into HEK293T cells; polysome sucrose-gradient fractionation and time-course sampling measured ribosome load, in-cell stability, and in-solution stability by barcode sequencing.
  • In-line-seq: Eterna “Roll-your-own-structure” challenge generated 68-nt RNAs; 3,030 fragments were profiled by single-nucleotide-resolution in-line hydrolysis and SHAPE, then DegScore was trained as a windowed ridge regression model and validated against independent PERSIST-seq in-solution half-lives.
  • Follow-up one-by-one assays used luciferase activity after transfection, capillary electrophoresis for in-solution half-lives, SHAPE/DMS chemical mapping, and ψ or m1ψ substitutions; design algorithms included LinearDesign, GC-rich codon sampling, RiboTree, and Eterna participant designs.

Limitations

  • Stated: experiments were limited to HEK293T cells, so optimal UTRs are likely cell-type dependent; the study focused on naked, non-formulated mRNAs, and formulated mRNA libraries may reveal additional effects; DegScore was trained on unmodified nucleotides and does not account for ψ stabilization; blinding and randomization were not used.
  • Evident from the text: the main pooled library contained 233 constructs, with some design rules derived from Eterna-designed subsets and accelerated degradation conditions (pH 10 with Mg2+), and the protein-output model assumes first-order translation and decay kinetics with negligible Nluc protein degradation.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41467-022-28776-w

Cites

115 of 115 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Figure 1. PERSIST-seq overview and illustrative ribosome load insights. (A) Overview of the mRNA optimization workflow. Literature mined and rationally designed 5’ and 3’ UTRs were combined with Eterna and algorithmically designed coding sequences. All sequences were then experimentally tested in parallel for in-solution and in-cell stability as well as ribosome load. The mRNA design included unique, 6–9 nt barcodes in the 3’ UTR for tag counting by short-read sequencing. (B) Experimental design for testing in-solution and in-cell stability and ribosome load in parallel. mRNAs were in vitro transcribed, 5’ capped, and polyadenylated in a pooled format before transfection into HEK293T cells or being subjected to in-solution degradation. Transfected cells were then harvested for sucrose gradient fractionation or in-cell degradation analysis. (C) Polysome trace from transfected HEK293T cells with 233-mRNA pool. (D) 5’ UTR variants display a higher variance in mean ribosome load per construct as determined from polysome sequencing. The formula for ribosome load is given. (E) Heatmaps from polysome profiles of mRNA designs selected from the top, middle, and bottom five mRNAs (by ribosome load) from each design category. (F) Secondary structure model of the SARS-CoV-2 5’ UTR. Introduced mutations and substitutions are highlighted. (G) Heatmaps of SARS-CoV-2 5’ UTR variants’ polysome profiles sorted by ribosome load.
Figure 1. PERSIST-seq overview and illustrative ribosome load insights. (A) Overview of the mRNA optimization workflow. Literature mined and rationally designed 5’ and 3’ UTRs were combined with Eterna and algorithmically designed coding sequences. All sequences were then experimentally tested in parallel for in-solution and in-cell stability as well as ribosome load. The mRNA design included unique, 6–9 nt barcodes in the 3’ UTR for tag counting by short-read sequencing. (B) Experimental design for testing in-solution and in-cell stability and ribosome load in parallel. mRNAs were in vitro transcribed, 5’ capped, and polyadenylated in a pooled format before transfection into HEK293T cells or being subjected to in-solution degradation. Transfected cells were then harvested for sucrose gradient fractionation or in-cell degradation analysis. (C) Polysome trace from transfected HEK293T cells with 233-mRNA pool. (D) 5’ UTR variants display a higher variance in mean ribosome load per construct as determined from polysome sequencing. The formula for ribosome load is given. (E) Heatmaps from polysome profiles of mRNA designs selected from the top, middle, and bottom five mRNAs (by ribosome load) from each design category. (F) Secondary structure model of the SARS-CoV-2 5’ UTR. Introduced mutations and substitutions are highlighted. (G) Heatmaps of SARS-CoV-2 5’ UTR variants’ polysome profiles sorted by ribosome load.CC BY-NC · PubMed Central
Figure 2. In-cell RNA stability drives downstream protein expression levels. (A) In-cell half-life of each mRNA design in HEK293T cells. (B) Higher polysome load correlates with decreased in-cell half-life. Correlation between in-cell half-life and mean ribosome load across the entire profile (left), monosome-to-free subunit ratio (center), or polysome-to-monosome ratio (right). (C) In-cell half-life and mean ribosome load for individual mRNA designs with varying UTRs. (D) Kinetic model for predicting protein expression from mRNA half-life and ribosome load. P(t) is protein quantity at time t; m0 is the mass of mRNA present at t=0; l is mRNA length; kt is translation rate; and km and kp are rates of mRNA and protein decay, respectively. (E) Protein expression predicted using the kinetic model in (D) on the basis of mRNA half-life and ribosome load. Top: predicted protein expression of each UTR variant; note closer similarity to in-cell half-life data than to ribosome load in (C). Bottom: predicted protein expression normalized by mRNA length (corresponding to transfecting equal masses of each mRNA). (F) Correlation of predicted protein expression and Nluc/Fluc activity in HEK293T cells. (G) In-solution half-life of various mRNA design variants. mRNA lifetimes are strongly dependent on mRNA length and designed structures, revealed by time courses of mRNA degradation under accelerated aging conditions (10 mM MgCl2, 50 mM Na-CHES, pH 10.0). (H) Nucleotide-resolution in vitro DMS mapping confirms large differences in structural accessibility between a highly structured JEV-HA-Nluc mRNA construct, “LinearDesign-1” and a highly unstructured construct “Yellowstone”. The 5’ and 3’ UTRs (hHBB) were kept constant between designs. (I) Nucleotide DMS accessibility mapped on to structures from DMS-directed structure prediction.
Figure 2. In-cell RNA stability drives downstream protein expression levels. (A) In-cell half-life of each mRNA design in HEK293T cells. (B) Higher polysome load correlates with decreased in-cell half-life. Correlation between in-cell half-life and mean ribosome load across the entire profile (left), monosome-to-free subunit ratio (center), or polysome-to-monosome ratio (right). (C) In-cell half-life and mean ribosome load for individual mRNA designs with varying UTRs. (D) Kinetic model for predicting protein expression from mRNA half-life and ribosome load. P(t) is protein quantity at time t; m0 is the mass of mRNA present at t=0; l is mRNA length; kt is translation rate; and km and kp are rates of mRNA and protein decay, respectively. (E) Protein expression predicted using the kinetic model in (D) on the basis of mRNA half-life and ribosome load. Top: predicted protein expression of each UTR variant; note closer similarity to in-cell half-life data than to ribosome load in (C). Bottom: predicted protein expression normalized by mRNA length (corresponding to transfecting equal masses of each mRNA). (F) Correlation of predicted protein expression and Nluc/Fluc activity in HEK293T cells. (G) In-solution half-life of various mRNA design variants. mRNA lifetimes are strongly dependent on mRNA length and designed structures, revealed by time courses of mRNA degradation under accelerated aging conditions (10 mM MgCl2, 50 mM Na-CHES, pH 10.0). (H) Nucleotide-resolution in vitro DMS mapping confirms large differences in structural accessibility between a highly structured JEV-HA-Nluc mRNA construct, “LinearDesign-1” and a highly unstructured construct “Yellowstone”. The 5’ and 3’ UTRs (hHBB) were kept constant between designs. (I) Nucleotide DMS accessibility mapped on to structures from DMS-directed structure prediction.CC BY-NC · PubMed Central
Figure 3. High-throughput in-line hydrolysis uncovers principles of in-solution RNA degradation. (A) Eterna participants were asked to design 68-nucleotide RNA fragments maximizing sequence and structure diversity. 3030 constructs were characterized and probed using high-throughput in-line degradation (In-line-seq). (B) Nucleotide-resolution degradation of 2165 68-nt RNA sequences (filtered for signal quality), probed by In-line-seq, sorted by hierarchical clustering on degradation profiles. (C) Sequences span a diverse set of secondary structure motifs, revealing patterns in degradation based on both sequence (i.e., linkages ending at 3’ uridine are particularly reactive) and structure (symmetric internal loops, circled, have suppressed hydrolytic degradation compared to asymmetric internal loops). (D) The ridge regression model “DegScore” was trained to predict per-nucleotide degradation from sequence and loop assignment information. Coefficients with the largest magnitude corresponded to sequence identity immediately after the link, with U being most disfavored. (E) DegScore showed improved predictive power on mRNAs over two other metrics previously posited to predict RNA stability. (F) Introduction of pseudouridine (ψ) modifications stabilizes selected short RNAs at U nucleotides in both loop motifs and in fully unstructured RNAs. (G) Capillary electrophoresis characterization of fragmentation time courses of Nluc mRNA molecules designed with extensive structure (LinearDesign-1) and relatively less structure (Yellowstone), synthesized with standard nucleotides and with ψ modifications. The full-length mRNA band is indicated with a red asterisk. The Tetrahymena ribozyme P4-P6 domain RNA was included after degradation as a control. (H) Exponential fits of capillary electrophoresis measurements of intact RNA over ten time points confirm >3x difference between in-solution lifetimes of LinearDesign-1 and Yellowstone Nluc mRNAs. Inset: Calculated half-lives. Error bars represent standard deviations from bootstrapped exponential fits (n = 1000).
Figure 3. High-throughput in-line hydrolysis uncovers principles of in-solution RNA degradation. (A) Eterna participants were asked to design 68-nucleotide RNA fragments maximizing sequence and structure diversity. 3030 constructs were characterized and probed using high-throughput in-line degradation (In-line-seq). (B) Nucleotide-resolution degradation of 2165 68-nt RNA sequences (filtered for signal quality), probed by In-line-seq, sorted by hierarchical clustering on degradation profiles. (C) Sequences span a diverse set of secondary structure motifs, revealing patterns in degradation based on both sequence (i.e., linkages ending at 3’ uridine are particularly reactive) and structure (symmetric internal loops, circled, have suppressed hydrolytic degradation compared to asymmetric internal loops). (D) The ridge regression model “DegScore” was trained to predict per-nucleotide degradation from sequence and loop assignment information. Coefficients with the largest magnitude corresponded to sequence identity immediately after the link, with U being most disfavored. (E) DegScore showed improved predictive power on mRNAs over two other metrics previously posited to predict RNA stability. (F) Introduction of pseudouridine (ψ) modifications stabilizes selected short RNAs at U nucleotides in both loop motifs and in fully unstructured RNAs. (G) Capillary electrophoresis characterization of fragmentation time courses of Nluc mRNA molecules designed with extensive structure (LinearDesign-1) and relatively less structure (Yellowstone), synthesized with standard nucleotides and with ψ modifications. The full-length mRNA band is indicated with a red asterisk. The Tetrahymena ribozyme P4-P6 domain RNA was included after degradation as a control. (H) Exponential fits of capillary electrophoresis measurements of intact RNA over ten time points confirm >3x difference between in-solution lifetimes of LinearDesign-1 and Yellowstone Nluc mRNAs. Inset: Calculated half-lives. Error bars represent standard deviations from bootstrapped exponential fits (n = 1000).CC BY-NC · PubMed Central
Figure 4. Integration of 5’/3’ UTRs, structure-optimized CDSs, and pseudouridine (ψ) together enhance mRNA stability and translational output. (A) CDS and 5’/3’UTR combinations differentially impact protein synthesis. Six mRNA constructs were in vitro synthesized and luciferase activity was measured 6 or 24 hrs post-transfection. Inclusion of ψ was tested on two selected constructs. (B) Workflow for different approaches to design the CDS variants tested in (C). (C) Variations in CDS design facilitates high in-solution stability and differential protein expression. In vitro transcribed mRNAs (24 in total) were subjected to in-solution degradation or transfected into HEK293T cells for 6 and 24 hrs. In-solution half-lives and luciferase activity are normalized to the Nluc START reference construct. Predicted secondary structures are shown for select constructs with colors indicating DegScore at each nucleotide. Designs derived from LinearDesign solutions are marked with a purple triangle. (D) Predicted secondary structure overview of Ribotree_LinearDesign_degscoreall_1. Zoomed boxes indicate sequence optimizations and subsequent structural changes made by DegScore to the reference LinearDesign construct. (E) Increased in-solution half-life and enhanced luciferase expression at 24 hrs, but not 6 hrs, correlate with DegScore. (F) Schematic for testing the synergy between RNA modifications and mRNA design rules on downstream stability and protein output. mRNAs were in vitro synthesized with or without ψ and subjected to degradation conditions. Samples were collected overtime and the RNA was purified before being transfected into HEK293T cells. Luciferase activity was measured 24 hrs after transfection. (G) Luciferase activity of the reference Nluc sequence and DegScore-optimized CDS with or without ψ after being subjected to in-solution degradation.
Figure 4. Integration of 5’/3’ UTRs, structure-optimized CDSs, and pseudouridine (ψ) together enhance mRNA stability and translational output. (A) CDS and 5’/3’UTR combinations differentially impact protein synthesis. Six mRNA constructs were in vitro synthesized and luciferase activity was measured 6 or 24 hrs post-transfection. Inclusion of ψ was tested on two selected constructs. (B) Workflow for different approaches to design the CDS variants tested in (C). (C) Variations in CDS design facilitates high in-solution stability and differential protein expression. In vitro transcribed mRNAs (24 in total) were subjected to in-solution degradation or transfected into HEK293T cells for 6 and 24 hrs. In-solution half-lives and luciferase activity are normalized to the Nluc START reference construct. Predicted secondary structures are shown for select constructs with colors indicating DegScore at each nucleotide. Designs derived from LinearDesign solutions are marked with a purple triangle. (D) Predicted secondary structure overview of Ribotree_LinearDesign_degscoreall_1. Zoomed boxes indicate sequence optimizations and subsequent structural changes made by DegScore to the reference LinearDesign construct. (E) Increased in-solution half-life and enhanced luciferase expression at 24 hrs, but not 6 hrs, correlate with DegScore. (F) Schematic for testing the synergy between RNA modifications and mRNA design rules on downstream stability and protein output. mRNAs were in vitro synthesized with or without ψ and subjected to degradation conditions. Samples were collected overtime and the RNA was purified before being transfected into HEK293T cells. Luciferase activity was measured 24 hrs after transfection. (G) Luciferase activity of the reference Nluc sequence and DegScore-optimized CDS with or without ψ after being subjected to in-solution degradation.CC BY-NC · PubMed Central

Verification

Claims the paper makes

93 Own result · 38 Own interpretation · 15 Own method · 8 No basis given · 5 Couldn’t confirm · 3 Partly supported · 2 Supported so far · 1 Supported · 1 Not supported · 1 No full text of the sources · accuracy not yet measured

  • It may state two separate claims as one. 22
  • Probably not a substantial claim. 15
  • The quote only partly states this claim. 1
36 of 36
  • The authors state this as an assumption.

    Nucleotides that are not base paired in structure are uniformly prone to hydrolytic degradation.

    “the assumption that nucleotides that are not base paired in structure would be uniformly prone to hydrolytic degradation7.”
  • The authors state this as an assumption.

    Translation follows first-order kinetics and mRNA and protein decay are exponential.

    “assuming first-order rate of translation and exponential decay for mRNAs and for proteins.”
  • The authors state this as an assumption.

    The high pH buffer containing Mg2+ mimics the high effective pH and positively charged environment that can arise in lipid nanoparticles, protamine, and other formulations for mRNA therapeutics.

    “To mimic the high effective pH and positively charged environment that can arise in lipid nanoparticles, protamine, and other formulations for mRNA therapeutics6, we used a high pH buffer containing Mg2+ to accelerate degradation (10 mM MgCl2, 50 mM Na-CHES, pH 10.0, 24 °C); conditions without Mg2+ or at lower pH lead to similar conclusions on relative stabilities across RNA variants (see below).”
  • The paper relies on this without support.

    The effects of ψ modifications on mRNA stability in solution have not been reported.

    “While ψ and derivatives have been reported to stabilize mRNAs against decay in cells67, the effects of these modifications on mRNA stability in solution have not been reported.”
  • The paper relies on this without support.

    High mRNA structure adversely impacts the cellular translational apparatus.

    “the expectation that high mRNA structure should adversely impact the cellular translational apparatus.”
  • The paper relies on this without support.

    Once delivered into patient cells, the candidate mRNA must outcompete other cytoplasmic mRNAs for the translation machinery and avoid the cellular mRNA degradation machinery to express maximum amounts of the desired encoded protein.

    “Moreover, once delivered into patient cells, the candidate mRNA must outcompete other cytoplasmic mRNAs for the translation machinery and avoid the cellular mRNA degradation machinery to express maximum amounts of the desired encoded protein.”
  • The authors state this as an assumption.

    mRNA translation and mRNA/protein decay follow first-order kinetics.

    “To calculated the expected protein levels assuming first order kinetics of mRNA translation and mRNA/protein decay, the following differential equations are used:”
  • The paper relies on this without support.

    The luciferase assay relies on accumulation of luciferase in the cytosol.

    “Because this assay relies on accumulation of luciferase in the cytosol, any signal peptides sequences (Table S1) were removed from the CDS for templates and mRNA for these transfection and luciferase activity experiments.”
  • The authors state this as an assumption.

    The majority of the other RNA in the lane were degradation products from the mRNA of interest.

    “assuming that the majority of the other RNA in the lane were degradation products from the mRNA of interest”
9 of 9

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

3 of 3

Citations

Cites

Show all 117 references

Cited by

Sorted by citations

See all 409

Metadata

Authors
Leppek, K., Byeon, G. W., Kladwang, W., Wayment-Steele, H. K., Kerr, C. H., Xu, A. F., Kim, D. S., Topkar, V. V., Choe, C., Rothschild, D., Tiu, G. C., Wellington-Oguri, R., Fujii, K., Sharma, E., Watkins, A. M., Nicol, J. J., Romano, J., Tunguz, B., Diaz, F., Cai, H., Guo, P., Wu, J., Meng, F., Shi, S., Participants, E., Dormitzer, P. R., Solórzano, A., Barna, M., & Das, R.
Affiliations
  • Stanford University (US)
  • Eterna Massive Open Laboratory (US)
  • Nvidia (United States) (US)
  • Pfizer (United States) (US)
Corresponding author
  • Barna, Maria
  • Das, Rhiju
Venue
Nature Communications
Published
2022-03-22
Volume
13
Issue
1
Pages
1536-1536
ISSN
2041-1723 (ISSN-L)
OpenAlex
W4220991830
Cited by
399
References
117
Citations by year
2026: 88 · 2025: 127 · 2024: 99 · 2023: 53 · 2022: 23
FWCI
24.66
Type
journal article
Language
en
Open access
Locations
  • Nature Communications (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • PubMed Central (submittedVersion, cc-by-nc)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
Funders
  • National Science Foundation
  • Damon Runyon Cancer Research Foundation (DRG-2237-15)
  • Pfizer
  • Oregon State University
  • Baidu
  • National Institutes of Health
  • University of California, San Francisco
  • Stanford Bio-X
  • Chemistry, Engineering and Medicine for Human Health, Stanford University
  • Canadian Institutes of Health Research
  • National Cancer Institute (R21CA219847)
  • National Institute of General Medical Sciences (R35GM122579)
  • Eunice Kennedy Shriver National Institute of Child Health and Human Development (F30HD100123, R01HD086634)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA and protein synthesis mechanisms
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Interference and Gene Delivery
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Research and Splicing
Keywords
mRNA stability, mRNA structure, mRNA vaccines, RNA therapeutics, protein expression, hydrolytic degradation, pseudouridine modification, mRNA design, COVID-19 vaccines

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-11-05, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-11-05 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review