Cited by 203RNA and protein synthesis mechanismsSARS-CoV-2 and COVID-19 ResearchRNA Interference and Gene Delivery

The DOI matches its Crossref record; 93 of the 93 references it cites were found in a registry or an index; 70 of the 136 sources its claims cite were checked for support.

Detailed Dissection and Critical Evaluation of the Pfizer/BioNTech and Moderna mRNA Vaccines

Open access
Xia, Xuhua·2021·Vaccines

Abstract

The design of Pfizer/BioNTech and Moderna mRNA vaccines involves many different types of optimizations. Proper optimization of vaccine mRNA can reduce dosage required for each injection leading to more efficient immunization programs. The mRNA components of the vaccine need to have a 5′-UTR to load ribosomes efficiently onto the mRNA for translation initiation, optimized codon usage for efficient translation elongation, and optimal stop codon for efficient translation termination. Both 5′-UTR and the downstream 3′-UTR should be optimized for mRNA stability. The replacement of uridine by N1-methylpseudourinine (Ψ) complicates some of these optimization processes because Ψ is more versatile in wobbling than U. Different optimizations can conflict with each other, and compromises would need to be made. I highlight the similarities and differences between Pfizer/BioNTech and Moderna mRNA vaccines and discuss the advantage and disadvantage of each to facilitate future vaccine improvement. In particular, I point out a few optimizations in the design of the two mRNA vaccines that have not been performed properly.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

Xia dissects the design of Pfizer/BioNTech BNT-162b2 and Moderna mRNA-1273, comparing 5′-UTRs, codon optimization, translation initiation, termination signals, and 3′-UTRs. He argues several optimizations are suboptimal—especially Moderna’s CGG codon choice and both vaccines’ UGA-containing stop signals—while noting some interpretations remain uncertain.

Key findings

  • Both vaccines increased CGN Arg codons and reduced AGR codons: BNT-162b2 reduced AGR by 8 and mRNA-1273 by 28; highly expressed ribosomal protein genes encode 64.2% of Arg by CGN. However, the paper argues CGG is not the optimal CGN codon: highly expressed skeletal muscle genes prefer CGC over CGG (paired t-test, t = 3.075, DF = 25, p = 0.0034), and mRNA-1273 overused CGG.
  • BNT162b2 uses the human α-globin 5′-UTR with Kozak GCCACCAUG and little secondary structure flanking the start codon; mRNA-1273 has a GC-rich element and secondary structure with MFE −12.3 for the 40-nt window around the start codon, versus −11.1 for α-globin. The paper says BNT162b2’s 5′-UTR appears superior but cautions that 8.6% of ribosomal protein genes have MFE equal to or smaller than −12.3.
  • BNT162b2 ends with UGAUGA and mRNA-1273 with UGAUAAUAG; the paper argues the optimal stop signal should be UAAA, because UGA is leakiest and highly expressed RP genes favor +4R (likelihood ratio chi-square = 931.0514, DF = 17, p < 0.000001). Ψ substitution may further increase stop-codon misreading and +1 frameshifting.
  • Given 100 μg mRNA-1273 versus 30 μg BNT-162b2 and equivalent efficacy, the paper infers Pfizer/BioNTech mRNA may produce about 3.3 times as many proteins per mRNA. CAI values are 0.94925 for BNT-162b2 and 0.97939 for mRNA-1273, but higher CAI is not necessarily better because Moderna used a “fundamentalist” and Pfizer/BioNTech a “liberal” codon strategy.

Methods

  • Sequence data included 33 RPS and 50 RPL ubiquitously expressed ribosomal protein genes, 11,327 annotated CDSs from human chromosome 1, intron nucleotide frequencies from chromosomes 18–22, SARS-CoV/SARS-CoV-2 and related coronavirus genomes, and vaccine mRNA sequences (BNT-162b2 public; mRNA-1273 putative from Jeong et al., with UTRs matching Moderna patents).
  • Tissue expression used Human Protein Atlas data from 62 tissues; top 50 skeletal-muscle highly expressed genes were selected, with 26 having CDSs ≥300 codons for Table 2; ZAP expression was extracted from three HPA files.
  • Analyses used DAMBE for CDS extraction, codon frequencies, CAI, ITE, PWM, and MFE (Vienna RNA fold); SARS-CoV-2 transcriptomic SRA SRR11886744 was analyzed for subgenomic spike mRNA; statistical tests included paired t-test and likelihood ratio test.

Limitations

  • Stated: the mRNA-1273 sequence is putative and “should be considered putative,” though its UTRs match Moderna patents; there is no systematic evaluation of the fundamentalist versus liberal codon optimization strategies in translation efficiency; the 5′-UTR secondary-structure interpretation is cautioned because some highly expressed genes have comparable MFE; and Ψ replacement complicates stop-codon reasoning.
  • Evident: the study is computational/in silico and reports no new wet-lab translation or immunogenicity experiments; inferences about dosage, protein output, and translation efficiency are indirect.
  • Evident: the paper itself notes that tRNA gene copy number is problematic as a proxy for tRNA abundance because nearly half of human tRNA genes are not expressed.

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

Integrity signals

Registry

Crossref · 10.3390/vaccines9070734

Cites

93 of 93 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Figure 1 Tissue-specific gene expression of ZAP (ZC3HAV1 and its long form ZC3HAV1L) extracted from three gene expression files from Human Protein Atlas (www.proteinatlas.org, accessed on 20 June 2021) representing three independent transcriptomic experiments: (A) proteinatlas.tsv, (B) rna_tissue_gtex.tsv and (C) rna_tissue.fantom.txv. ZAP expression is low in skeletal muscle (pointed to by the red arrow) in all three data sets. The horizontal axis is alphabetically sorted. The three data sets do not include the same types of tissues.
Figure 1 Tissue-specific gene expression of ZAP (ZC3HAV1 and its long form ZC3HAV1L) extracted from three gene expression files from Human Protein Atlas (www.proteinatlas.org, accessed on 20 June 2021) representing three independent transcriptomic experiments: (A) proteinatlas.tsv, (B) rna_tissue_gtex.tsv and (C) rna_tissue.fantom.txv. ZAP expression is low in skeletal muscle (pointed to by the red arrow) in all three data sets. The horizontal axis is alphabetically sorted. The three data sets do not include the same types of tissues.CC BY · PubMed Central
Figure 2 5′-UTR and secondary structure of the two vaccine mRNAs. (A) 5′-UTR of BNT162b2 was partially taken from human α-globin gene where the 37-nt 5′-UTR is shared by HBA1 and HBA2 mRNAs. The upstream TATA box is colored in red, and the Kozak consensus highlighted in bold. (B) Secondary structure stability measured by MFE (minimum folding energy) over a sliding windows of 40 nt, with start codon AUG at Mid-window sites 55–57. A strong secondary structure (small MFE) is visible in sequences flanking the start codon in mRNA-1273. (C) Visualization of the secondary structure embedding the start codon in mRNA-1273.
Figure 2 5′-UTR and secondary structure of the two vaccine mRNAs. (A) 5′-UTR of BNT162b2 was partially taken from human α-globin gene where the 37-nt 5′-UTR is shared by HBA1 and HBA2 mRNAs. The upstream TATA box is colored in red, and the Kozak consensus highlighted in bold. (B) Secondary structure stability measured by MFE (minimum folding energy) over a sliding windows of 40 nt, with start codon AUG at Mid-window sites 55–57. A strong secondary structure (small MFE) is visible in sequences flanking the start codon in mRNA-1273. (C) Visualization of the secondary structure embedding the start codon in mRNA-1273.CC BY · PubMed Central
Figure 3 MFE (minimum fold energy) distribution of sequences flanking the start codon AUG in ribosomal protein genes. MFE is calculated from 40-nt with 20 nt upstream of AUG and 20 nt downstream including AUG), based on 186 mRNA variants of 34 RPS and 54 RPL genes. A total of 8.6% of MFE values are equal to, or smaller (more stable secondary structure) than, MFE for mRNA-1273.
Figure 3 MFE (minimum fold energy) distribution of sequences flanking the start codon AUG in ribosomal protein genes. MFE is calculated from 40-nt with 20 nt upstream of AUG and 20 nt downstream including AUG), based on 186 mRNA variants of 34 RPS and 54 RPL genes. A total of 8.6% of MFE values are equal to, or smaller (more stable secondary structure) than, MFE for mRNA-1273.CC BY · PubMed Central
Figure 4 Discontinuous transcription of subgenomic SARS-2-S mRNA. (A) SARS-CoV-2 genome (NC_045512) with the leader sequence at sites 1–70, and two shaded transcription regulatory sequence (TRS), with TRS-L at the end of the leader sequence and TRS-B upstream of the coding sequence for the spike protein colored in red. The transcription of the negative strand progresses up to TRS-B and pauses. (B) The negative strand shifts location from base-pairing with TRS-B to base-pairing with TRS-L. Transcription of the negative strand continues. (C) SARS-2-S mRNA transcribed from the negative subgenomic strand. (D) Three frequently used TRS in coronavirus, with the first two likely being degenerate forms of the third.
Figure 4 Discontinuous transcription of subgenomic SARS-2-S mRNA. (A) SARS-CoV-2 genome (NC_045512) with the leader sequence at sites 1–70, and two shaded transcription regulatory sequence (TRS), with TRS-L at the end of the leader sequence and TRS-B upstream of the coding sequence for the spike protein colored in red. The transcription of the negative strand progresses up to TRS-B and pauses. (B) The negative strand shifts location from base-pairing with TRS-B to base-pairing with TRS-L. Transcription of the negative strand continues. (C) SARS-2-S mRNA transcribed from the negative subgenomic strand. (D) Three frequently used TRS in coronavirus, with the first two likely being degenerate forms of the third.CC BY · PubMed Central

Verification

Claims the paper makes

47 Own interpretation · 39 No basis given · 28 Partly supported · 20 Own result · 13 No full text of the sources · 8 Supported · 7 Not checked yet · 5 Supported so far · 2 Mixed support · 2 Couldn’t confirm · 1 Not supported · 1 Own method · accuracy not yet measured

Verify all 12 references — free account

Show all 118 claims

  • It may state two separate claims as one. 41
  • Probably not a substantial claim. 16
  • The quote only partly states this claim. 1
55 of 55
  • The paper relies on this without support.

    Proper optimization of vaccine mRNA can reduce dosage required for each injection leading to more efficient immunization programs.

    “Proper optimization of vaccine mRNA can reduce dosage required for each injection leading to more efficient immunization programs.”
  • The authors state this as an assumption.

    The putative sequence of mRNA-1273 is the real mRNA-1273.

    “In this context, it may be safe to assume that the putative sequence is the real mRNA-1273.”
  • The paper relies on this without support.

    The two mRNA vaccines are administered through muscle injection.

    “As the two mRNA vaccines are administered through muscle injection, it is relevant to characterize features of highly expressed muscle genes.”
  • The authors state this as an assumption.

    tRNA abundance is well correlated with tRNA gene copy number.

    “Assuming that tRNA abundance is well correlated with tRNA gene copy number, which is true for Saccharomyces cerevisiae [34] but not known for other eukaryotes”
  • The authors state this as an assumption.

    An injection of Pfizer/BioNTech or Moderna vaccine produces the same number of the encoded spike proteins.

    “one may assume that an injection of Pfizer/BioNTech or Moderna vaccine produces the same number of the encoded spike proteins.”
  • The paper relies on this without support.

    All U nucleotides in the two mRNA vaccines were replaced by N1-methylpseudouridines (Ψ).

    “All U nucleotides in the two mRNA vaccines were replaced by N1-methylpseudouridines (Ψ)”
  • The paper relies on this without support.

    The Pfizer/BioNTech and Moderna vaccines are equivalent in efficacy despite the difference in dosage.

    “Given the difference in dosage (100 μg with mRNA-1273 and 30 μg with BNT-162b2) and the equivalence in efficacy”
  • The paper relies on this without support.

    There is no systematic evaluation of the two codon optimization strategies in translation efficiency.

    “There is no systematic evaluation of these two codon optimization strategies in translation efficiency.”
  • The paper relies on this without support.

    In E. coli, AAC is a major codon for Asn and AAU is a rarely used minor codon.

    “AAC is a major codon (heavily used by highly expressed genes and decoded by the most abundant isoacceptor tRNA {XE “tRNA: isoacceptor”}{XE “isoacceptor tRNA”}) whereas AAU is a rarely used minor codon.”
  • The paper relies on this without support.

    Highly expressed E. coli genes use AAC almost exclusively to encode Asn.

    “Highly expressed E. coli genes use AAC almost exclusively to encode Asn”
  • The authors state this as an assumption.

    Conserved amino acid sites are functionally important and variable amino acid sites are of limited importance.

    “He classified amino acid sites into conserved sites (assumed to be functionally important) and variable sites (assumed to be of limited importance).”
  • The authors state this as an assumption.

    The optimization done by natural selection can be extrapolated to translation in muscle cells.

    “These first two strategies assume that the optimization done by natural selection can be extrapolated to translation in muscle cells.”
  • The paper relies on this without support.

    The 5′-UTR serves two key functions: to stabilize mRNA and to facilitate scanning by the small ribosome subunit to localize the start codon.

    “5′-UTR serves two key functions: to stabilize mRNA and to facilitate scanning by small ribosome subunit to localize the start codon.”
  • The paper relies on this without support.

    G/U base pairs are usually too weak to close a stem.

    “G/U base pairs are usually too weak to close a stem.”
  • The paper relies on this without support.

    Functionally important and highly expressed genes such as ribosomal protein genes are expected to evolve strong termination signals.

    “Functionally important and highly expressed genes such as RP genes are expected to evolve strong termination signals.”
  • The authors state this as an assumption.

    The observation that glycine interacts more favorably with purine than pyrimidine in protein-DNA interactions can be extrapolated to protein-RNA interactions.

    “If one may extrapolate this observation from DNA to RNA, then one would predict that G626 would favor purine at the +4 site, i.e., a purine at +4 site is more likely to enhance the stop signal by interacting with G626 in eRF1 than a pyrimidine.”
16 of 16

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?

Citations

Cites

Show all 92 references

Research works printed but not linked: 3

  1. Xia, X. (2018). Bioinformatics and Translation Elongation. Bioinformatics and the Cell: Modern Computational Approaches in Genomics, Proteomics and Transcriptomics, Springer. · DOI · Paper page
  2. Xia, X. (2018). Bioinformatics and Translation Initiation. Bioinformatics and the Cell: Modern Computational Approaches in Genomics, Proteomics and Transcriptomics, Springer. · DOI · Paper page
  3. Xia, X. (2018). Bioinformatics and Translation Termination in Bacteria. Bioinformatics and the Cell: Modern Computational Approaches in Genomics, Proteomics and Transcriptomics, Springer. · DOI · Paper page

Cited by

Sorted by citations

See all 203

Metadata

Authors
Xia, X.
Affiliations
University of Ottawa (CA)
Corresponding author
Xia, Xuhua
Venue
Vaccines
Publisher
MDPI AG
Published
2021-07-03
Volume
9
Issue
7
Pages
734
ISSN
2076-393X (ISSN-L)
OpenAlex
W3177545102
Cited by
201
References
92
Citations by year
2026: 28 · 2025: 40 · 2024: 42 · 2023: 49 · 2022: 34 · 2021: 8
FWCI
10.44
Type
journal article
Language
en
Open access
Locations
  • Vaccines (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • uO Research (University of Ottawa) (submittedVersion, cc-by-nd)
  • University of Ottawa - Library
Funders
Natural Sciences and Engineering Research Council of Canada (RGPIN/2018-03878)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA and protein synthesis mechanisms
  • Health Sciences › Medicine › Infectious Diseases › SARS-CoV-2 and COVID-19 Research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Interference and Gene Delivery
Keywords
mRNA vaccines, BNT162b2, Moderna vaccine, SARS-CoV-2, 5' untranslated region, codon usage, stop codon, 3' untranslated region, mRNA stability, N1-methylpseudouridine, wobble pairing, translation initiation, vaccine design

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