CodonBert large language model for mRNA vaccines.

mRNA-based vaccines and therapeutics are gaining popularity and usage across a wide range of conditions. One of the critical issues when designing such mRNAs is sequence optimization. Even small proteins or peptides can be encoded by an enormously large number of mRNAs. The actual mRNA sequence can have a large impact on several properties including expression, stability, immunogenicity, and more. To enable the selection of an optimal sequence, we developed CodonBERT, a large language model (LLM) for mRNAs. Unlike prior models, CodonBERT uses codons as inputs which enables it to learn better representations. CodonBERT was trained using more than 10 million mRNA sequences from a diverse set of organisms. The resulting model captures important biological concepts. CodonBERT can also be extended to perform prediction tasks for various mRNA properties. CodonBERT outperforms previous mRNA prediction methods including on a new flu vaccine dataset.

Understanding structure activity relationships of Good HEPES lipids for lipid nanoparticle mRNA vaccine applications.

Lipid nanoparticles (LNPs) have shown great promise as delivery vehicles to transport messenger ribonucleic acid (mRNA) into cells and act as vaccines for infectious diseases including COVID-19 and influenza. The ionizable lipid incorporated within the LNP is known to be one of the main driving factors for potency and tolerability. Herein, we describe a novel family of ionizable lipids synthesized with a piperazine core derived from the HEPES Good buffer. These ionizable lipids have unique asymmetric tails and two dissimilar degradable moieties incorporated within the structure. Lipids tails of varying lengths, degrees of unsaturation, branching, and the inclusion of additional ester moieties were evaluated for protein expression. We observed several key lipid structure activity relationships that correlated with improved protein production in vivo, including lipid tails of 12 carbons on the ester side and the effect of carbon spacing on the disulfide arm of the lipids. Differences in LNP physical characteristics were observed for lipids containing an extra ester moiety. The LNP structure and lipid bilayer packing, visualized through Cryo-TEM, affected the amount of protein produced in vivo. In non-human primates, the Good HEPES LNPs formulated with an mRNA encoding an influenza hemagglutinin (HA) antigen successfully generated functional HA inhibition (HAI) antibody titers comparable to the industry standards MC3 and SM-102 LNPs, demonstrating their promise as a potential vaccine.

Pan-SARS neutralizing responses after third boost vaccination in non-human primate immunogenicity model.

The emergence of SARS-CoV-2 variants, especially Beta and Delta, has raised concerns about the reduced protection from previous infection or vaccination based on the original Wuhan-Hu-1 (D614) virus. To identify promising regimens for inducing neutralizing titers towards new variants, we evaluated monovalent and bivalent mRNA vaccines either as primary vaccination or as a booster in nonhuman primates (NHPs). Two mRNA vaccines, D614-based MRT5500 and Beta-based MRT5500ß, tested in sequential regimens or as a bivalent combination in naïve NHPs produced modest neutralizing titers to heterologous variants. However, when mRNA vaccines were administered as a booster to pre-immune NHPs, we observed a robust increase in neutralizing titers with expanded breadth towards all tested variants, and notably SARS-CoV-1. The breadth of the neutralizing response was independent of vaccine sequence or modality, as we further showed either MRT5500 or recombinant subunit Spike protein (with adjuvant) can serve as boosters to induce broadly neutralizing antibodies in the NHPs primed with MRT5500. The data support the notion that a third vaccination is key to boosting existing titers and improving the breadth of antibodies to address variants of concern, including those with an E484K mutation in the Receptor Binding Domain (RBD) (Beta, Gamma).

Development of multivalent mRNA vaccine candidates for seasonal or pandemic influenza.

Recent approval of mRNA vaccines for emergency use against COVID-19 is likely to promote rapid development of mRNA-based vaccines targeting a wide range of infectious diseases. Compared to conventional approaches, this vaccine modality promises comparable potency while substantially accelerating the pace of development and deployment of vaccine doses. Already demonstrated successfully for single antigen vaccines such as for COVID-19, this technology could be optimized for complex multi-antigen vaccines. Herein, utilizing multiple influenza antigens, we demonstrated the suitability of the mRNA therapeutic (MRT) platform for such applications. Seasonal influenza vaccines have three or four hemagglutinin (HA) antigens of different viral subtypes. In addition, influenza neuraminidase (NA), a tetrameric membrane protein, is identified as an antigen that has been linked to protective immunity against severe viral disease. We detail the efforts in optimizing formulations of influenza candidates that use unmodified mRNA encoding full-length HA or full-length NA encapsulated in lipid nanoparticles (LNPs). HA and NA mRNA-LNP formulations, either as monovalent or as multivalent vaccines, induced strong functional antibody and cellular responses in non-human primates and such antigen-specific antibody responses were associated with protective efficacy against viral challenge in mice.

Immunogenicity and efficacy of mRNA COVID-19 vaccine MRT5500 in preclinical animal models.

Emergency use authorization of COVID vaccines has brought hope to mitigate pandemic of coronavirus disease 2019 (COVID-19). However, there remains a need for additional effective vaccines to meet the global demand and address the potential new viral variants. mRNA technologies offer an expeditious path alternative to traditional vaccine approaches. Here we describe the efforts to utilize an mRNA platform for rational design and evaluations of mRNA vaccine candidates based on the spike (S) glycoprotein of SARS-CoV-2. Several mRNA constructs of S-protein, including wild type, a pre-fusion stabilized mutant (2P), a furin cleavage-site mutant (GSAS) and a double mutant form (2P/GSAS), as well as others, were tested in animal models for their capacity to elicit neutralizing antibodies (nAbs). The lead 2P/GSAS candidate was further assessed in dose-ranging studies in mice and Cynomolgus macaques, and for efficacy in a Syrian golden hamster model. The selected 2P/GSAS vaccine formulation, designated MRT5500, elicited potent nAbs as measured in neutralization assays in all three preclinical models and more importantly, protected against SARS-CoV-2-induced weight loss and lung pathology in hamsters. In addition, MRT5500 elicited TH1-biased responses in both mouse and non-human primate (NHP), thus alleviating a hypothetical concern of potential vaccine-associated enhanced respiratory diseases known associated with TH2-biased responses. These data position MRT5500 as a viable vaccine candidate for entering clinical development.

Premature Termination Codons Are Recognized in the Nucleus in A Reading-Frame Dependent Manner.

mRNAs containing premature termination codons (PTCs) are known to be degraded via nonsense-mediated mRNA decay (NMD). Unexpectedly, we found that mRNAs containing any type of PTC (UAA, UAG, UGA) are detained in the nucleus whereas their wild-type counterparts are rapidly exported. This retention is strictly reading-frame dependent. Strikingly, our data indicate that translating ribosomes in the nucleus proofread the frame and detect the PTCs in the nucleus. Moreover, the shuttling NMD protein Upf1 specifically associates with PTC+ mRNA in the nucleus and is required for nuclear retention of PTC+ mRNA. Together, our data lead to a working model that PTCs are recognized in the nucleus by translating ribosomes, resulting in recruitment of Upf1, which in turn functions in nuclear retention of PTC+ mRNA. Nuclear PTC recognition adds a new layer of proofreading for mRNA and may be vital for ensuring the extraordinary fidelity required for protein production.

Export and stability of naturally intronless mRNAs require specific coding region sequences and the TREX mRNA export complex.

A great deal is known about the export of spliced mRNAs, but little is known about the export of mRNAs encoded by human cellular genes that naturally lack introns. Here, we investigated the requirements for export of three naturally intronless mRNAs (HSPB3, IFN-α1, and IFN-ß1). Significantly, we found that all three mRNAs are stable and accumulate in the cytoplasm, whereas size-matched random RNAs are unstable and detected only in the nucleus. A portion of the coding region confers this stability and cytoplasmic localization on the naturally intronless mRNAs and a cDNA transcript, which is normally retained in the nucleus and degraded. A polyadenylation signal, TREX mRNA export components, and the mRNA export receptor TAP are required for accumulation of the naturally intronless mRNAs in the cytoplasm. We conclude that naturally intronless mRNAs contain specific sequences that result in efficient packaging into the TREX mRNA export complex, thereby supplanting the splicing requirement for efficient mRNA export.

A role for TREX components in the release of spliced mRNA from nuclear speckle domains.

The TREX complex, which functions in mRNA export, is recruited to mRNA during splicing. Both the splicing machinery and the TREX complex are concentrated in 20-50 discrete foci known as nuclear speckle domains. In this study, we use a model system where DNA constructs are microinjected into HeLa cell nuclei, to follow the fates of the transcripts. We show that transcripts lacking functional splice sites, which are inefficiently exported, do not associate with nuclear speckle domains but are instead distributed throughout the nucleoplasm. In contrast, pre-mRNAs containing functional splice sites accumulate in nuclear speckles, and our data suggest that splicing occurs in these domains. When the TREX components UAP56 or Aly are knocked down, spliced mRNA, as well as total polyA+ RNA, accumulates in nuclear speckle domains. Together, our data raise the possibility that pre-mRNA undergoes splicing in nuclear speckle domains, before their release by TREX components for efficient export to the cytoplasm.

Human DDX3 functions in translation and interacts with the translation initiation factor eIF3.

The conserved RNA helicase DDX3 is of major medical importance due to its involvement in numerous cancers, human hepatitis C virus (HCV) and HIV. Although DDX3 has been reported to have a wide variety of cellular functions, its precise role remains obscure. Here, we raised a new antibody to DDX3 and used it to show that DDX3 is evenly distributed throughout the cytoplasm at steady state. Consistent with this observation, HA-tagged DDX3 also localizes to the cytoplasm. RNAi of DDX3 in both human and Drosophila cells shows that DDX3 is required for cell viability. Moreover, using RNAi, we show that DDX3 is required for expression of protein from reporter constructs. In contrast, we did not detect a role for DDX3 in nuclear steps in gene expression. Further insight into the function of DDX3 came from the observation that its major interaction partner is the multi-component translation initiation factor eIF3. We conclude that a primary function for DDX3 is in protein translation, via an interaction with eIF3.

Splicing promotes rapid and efficient mRNA export in mammalian cells.

The numerous steps in protein gene expression are extensively coupled to one another through complex networks of physical and functional interactions. Indeed, >25 coupled reactions, often reciprocal, have been documented among such steps as transcription, capping, splicing, and polyadenylation. Coupling is usually not essential for gene expression, but instead enhances the rate and/or efficiency of reactions and, physiologically, may serve to increase the fidelity of gene expression. Despite numerous examples of coupling in gene expression, whether splicing enhances mRNA export still remains controversial. Although splicing was originally reported to promote export in both mammalian cells and Xenopus oocytes, it was subsequently concluded that this was not the case. These newer conclusions were surprising in light of the observations that the mRNA export machinery colocalizes with splicing factors in the nucleus and that splicing promotes recruitment of the export machinery to mRNA. We therefore reexamined the relationship between splicing and mRNA export in mammalian cells by using FISH, in combination with either transfection or nuclear microinjection of plasmid DNA. Together, these analyses indicate that both the kinetics and efficiency of mRNA export are enhanced 6- to 10-fold (depending on the construct) for spliced mRNAs relative to their cDNA counterparts. We conclude that splicing promotes mRNA export in mammalian cells and that the functional coupling between splicing and mRNA export is a conserved and general feature of gene expression in higher eukaryotes.

The signal sequence coding region promotes nuclear export of mRNA.

In eukaryotic cells, most mRNAs are exported from the nucleus by the transcription export (TREX) complex, which is loaded onto mRNAs after their splicing and capping. We have studied in mammalian cells the nuclear export of mRNAs that code for secretory proteins, which are targeted to the endoplasmic reticulum membrane by hydrophobic signal sequences. The mRNAs were injected into the nucleus or synthesized from injected or transfected DNA, and their export was followed by fluorescent in situ hybridization. We made the surprising observation that the signal sequence coding region (SSCR) can serve as a nuclear export signal of an mRNA that lacks an intron or functional cap. Even the export of an intron-containing natural mRNA was enhanced by its SSCR. Like conventional export, the SSCR-dependent pathway required the factor TAP, but depletion of the TREX components had only moderate effects. The SSCR export signal appears to be characterized in vertebrates by a low content of adenines, as demonstrated by genome-wide sequence analysis and by the inhibitory effect of silent adenine mutations in SSCRs. The discovery of an SSCR-mediated pathway explains the previously noted amino acid bias in signal sequences and suggests a link between nuclear export and membrane targeting of mRNAs.

Expression and molecular characterization of ZmMYB-IF35 and related R2R3-MYB transcription factors.

R2R3-MYB transcription factors play many important roles in higher plants including the regulation of secondary metabolism, the control of cell shape, and in the response to various stress conditions. In spite of their large number and significance, very few of these genes have been functionally characterized in monocots. Here, we describe the characterization of ZmMYB-IF35 from maize. Using GAL4 fusion constructs, we show that ZmMYB-IF35 possesses the ability to bind DNA in a sequence specific manner and activate transcription in yeast. We also show that ZmMYB-IF35 is capable of binding to the a1 promoter in planta, but it is not sufficient for activation of a1 transcription. Interestingly, a chimeric protein consisting of the MYB domain from ZmMYB-IF35 and the non-MYB C-terminal region of P1, a closely related R2R3-MYB protein, activated transcription from the a1 promoter in planta, suggesting that regions outside the conserved R2R3-MYB domain contribute to regulatory specificity. In situ hybridization experiments demonstrate that ZmMYB-IF35 expresses primarily in epidermal and vascular cells, while its rice ortholog, OsMYB-IF, displays a broad expression pattern in aerial parts of the plant. Together, our results provide novel insights on the participation of ZmMYB-IF35 and related genes in the regulation of secondary metabolic pathways in the grasses.

The human Shwachman-Diamond syndrome protein, SBDS, associates with ribosomal RNA.

Shwachman-Diamond syndrome (SDS) is an autosomal recessive disorder characterized by bone marrow failure, exocrine pancreatic dysfunction, and leukemia predisposition. Mutations in the SBDS gene are identified in most patients with SDS. SBDS encodes a highly conserved protein of unknown function. Data from SBDS orthologs suggest that SBDS may play a role in ribosome biogenesis or RNA processing. Human SBDS is enriched in the nucleolus, the major cellular site of ribosome biogenesis. Here we report that SBDS nucleolar localization is dependent on active rRNA transcription. Cells from patients with SDS or Diamond-Blackfan anemia are hypersensitive to low doses of actinomycin D, an inhibitor of rRNA transcription. The addition of wild-type SBDS complements the actinomycin D hypersensitivity of SDS patient cells. SBDS migrates together with the 60S large ribosomal subunit in sucrose gradients and coprecipitates with 28S ribosomal RNA (rRNA). Loss of SBDS is not associated with a discrete block in rRNA maturation or with decreased levels of the 60S ribosomal subunit. SBDS forms a protein complex with nucleophosmin, a multifunctional protein implicated in ribosome biogenesis and leukemogenesis. Our studies support the addition of SDS to the growing list of human bone marrow failure syndromes involving the ribosome.

Human mRNA export machinery recruited to the 5' end of mRNA.

Pre-mRNAs undergo splicing to remove introns, and the spliced mRNA is exported to the cytoplasm for translation. Here we investigated the mechanism for recruitment of the conserved mRNA export machinery (TREX complex) to mRNA. We show that the human TREX complex is recruited to a region near the 5' end of mRNA, with the TREX component Aly bound closest to the 5' cap. Both TREX recruitment and mRNA export require the cap, and these roles for the cap are splicing dependent. CBP80, which is bound to the cap, associates efficiently with TREX, and Aly mediates this interaction. Together, these data indicate that the CBP80-Aly interaction results in recruitment of TREX to the 5' end of mRNA, where it functions in mRNA export. As a consequence, the mRNA would be exported in a 5' to 3' direction through the nuclear pore, as observed in early electron micrographs of giant Balbiani ring mRNPs.

Metabolite profiling as a functional genomics tool.

Plants accumulate a very large number of small molecules (phytochemicals) with important functions in the ecology of plants and in the protection against biotic and abiotic stress conditions. Little is known on how phytochemical biosynthetic pathways are regulated, which is a key step to successfully engineering plant metabolism. Plant natural products are usually not essential, and genetic analyses often fail to identify phenotypes associated with the absence of these compounds. We have investigated the use of metabolite profiling of plant cells in culture to establish the function of transcription factors suspected to control plant metabolic pathways.

Recently duplicated maize R2R3 Myb genes provide evidence for distinct mechanisms of evolutionary divergence after duplication.

R2R3 Myb genes are widely distributed in the higher plants and comprise one of the largest known families of regulatory proteins. Here, we provide an evolutionary framework that helps explain the origin of the plant-specific R2R3 Myb genes from widely distributed R1R2R3 Myb genes, through a series of well-established steps. To understand the routes of sequence divergence that followed Myb gene duplication, we supplemented the information available on recently duplicated maize (Zea mays) R2R3 Myb genes (C1/Pl1 and P1/P2) by cloning and characterizing ZmMyb-IF35 and ZmMyb-IF25. These two genes correspond to the recently expanded P-to-A group of maize R2R3 Myb genes. Although the origins of C1/Pl1 and ZmMyb-IF35/ZmMyb-IF25 are associated with the segmental allotetraploid origin of the maize genome, other gene duplication events also shaped the P-to-A clade. Our analyses indicate that some recently duplicated Myb gene pairs display substantial differences in the numbers of synonymous substitutions that have accumulated in the conserved MYB domain and the divergent C-terminal regions. Thus, differences in the accumulation of substitutions during evolution can explain in part the rapid divergence of C-terminal regions for these proteins in some cases. Contrary to previous studies, we show that the divergent C termini of these R2R3 MYB proteins are subject to purifying selection. Our results provide an in-depth analysis of the sequence divergence for some recently duplicated R2R3 Myb genes, yielding important information on general patterns of evolution for this large family of plant regulatory genes.