Optimization of 5' Untranslated Region of Modified mRNA for Use in Cardiac or Hepatic Ischemic Injury.

Modified mRNA (modRNA) is a gene-delivery platform for transiently introducing a single gene or several genes of interest to different cell types and tissues. modRNA is considered to be a safe vector for gene transfer, as it negligibly activates the innate immune system and does not compromise the genome integrity. The use of modRNA in basic and translational science is rising, due to the clinical potential of modRNA. We are currently using modRNA to induce cardiac regeneration post-ischemic injury. Major obstacles in using modRNA for cardiac ischemic disease include the need for the direct and single administration of modRNA to the heart and the inefficient translation of modRNA due to its short half-life. Modulation of the 5' untranslated region (5' UTR) to enhance translation efficiency in ischemic cardiac disease has great value, as it can reduce the amount of modRNA needed per delivery and will achieve higher and longer protein production post-single delivery. Here, we identified that 5' UTR, from the fatty acid metabolism gene carboxylesterase 1D (Ces1d), enhanced the translation of firefly luciferase (Luc) modRNA by 2-fold in the heart post-myocardial infarction (MI). Moreover, we identified, in the Ces1d, a specific RNA element (element D) that is responsible for the improvement of modRNA translation and leads to a 2.5-fold translation increment over Luc modRNA carrying artificial 5' UTR, post-MI. Importantly, we were able to show that 5' UTR Ces1d also enhances modRNA translation in the liver, but not in the kidney, post-ischemic injury, indicating that Ces1d 5' UTR and element D may play a wider role in translation of protein under an ischemic condition.

Introduction

Ischemic heart disease is the leading cause of death for both men and women in the United States, killing about 610,000 people per year., Many scientists around the world are studying different approaches to induce cardiac regeneration in hope to unravel new treatments that can improve cardiac function post-ischemic injury. Genetic medicine is one of the avenues of research to treat a failing heart post-myocardial infarction (MI). In this research, we aimed to adjust gene expression in the heart using viral vectors, small molecules, or an RNA-based approach to promote cardiac protection and cardiovascular or cardiac regeneration in ischemic cardiac disease. Modified messenger RNA (modRNA) is a novel gene therapy platform that can be used to alter the levels of proteins in mammalian cells and tissues, and to treat heart diseases.6, 7, 8 The concept of therapeutically altering mRNA expression has great potential in the treatment of human diseases. To date, several therapeutic approaches using small interfering RNA (siRNA) and antisense oligonucleotides have been successfully tested to reduce mRNA levels in the cells., Yet, the upregulation of proteins in tissues is challenging, mostly to the high amount of mRNA needed to treat human tissue. The supply of a large amount of mRNA in vivo can subsequently elicit a higher immune response. Preclinical studies with modRNA showed that due to the transient expression of modRNA (gene expression return to baseline is a few days), modRNA will be needed, both in direct or intravenous delivery methods, to be redelivered to achieve higher gene expression.,12, 13, 14, 15

The expression of genes is controlled intricately at the post-transcriptional level. The level of an individual mRNA type inside of a cell does not ensure the synthesis of comparable amounts of respective proteins, both positive and negative modulators that influence translation and maintain certain levels of protein. Eukaryotic gene translation is regulated in the translation level by several components, including 5′ untranslated region (UTR),18, 19, 20 3′ UTR,21, 22, 23 poly(A) tail,24, 25, 26, 27 cap structure,28, 29, 30, 31 etc. There are multiple regulatory elements within the untranslated regions (UTRs) of the mRNA, which are critical for the stability and translation of mRNA into protein., Moreover, the length and secondary structure of 5′ UTR and mutations within the 5′ UTR have been reported to be associated with certain human diseases. The 5′ UTR plays a significant role in the regulation of translational efficiency by helping the ribosome to bind the messenger RNA (mRNA) in the proximity of the start codon and thus, is a main contributor to the cellular proteome. Additionally, certain RNA elements within the 5′ UTR may change its secondary structure (e.g., internal ribosome entry sites [IRESs], upstream start codon (AUG)s, or upstream open reading frames [uORFs]) and can be an important contributor of the entire translation rate., Besides, 5′ UTRs can contain sequence elements that can function as binding sites for regulatory proteins.

Until this date, modRNA, which had been used preclinically in cardiac research using artificial 5′ UTR (36 nt), was first described by Warren et al.

Asrani et al. has used in vitro screening for optimization mRNA 5′ UTR to improve expression of arginase 1 (ARG1). They showed that plasmid-based screening methods do not correlate with protein expression driven by exogenously expressed mRNA and that improved 5′ UTR but not 3′ UTR appears to be the key driver in protein expression for exogenously delivered mRNA.

We recently optimized the synthesis and the delivery of modRNA for cardiac delivery use. In the present study, we compared transcriptomic and proteomic analysis of similar tissues to identify, in vivo, a potential 5′ UTR and the elements within 5′ UTR that can enhance the translation post-cardiac and hepatic ischemic injury condition.

Results

Characterizing the Ischemic Heart Transcriptome and Proteome

To characterize the dynamics of heart left-ventricle (LV) transcriptome and proteome post-MI, we analyzed changes in genes expression and protein level in the LV of mice, 4 and 24 h post-MI and compared it to LV from sham-operated mice (Figure 1A). In total, we detected 14,552 genes and 2,397 proteins in our samples. With the comparison of the two datasets, we found 2,272 genes with corresponding proteins. Out of all genes with corresponding proteins, 239 genes and 120 proteins were differentially expressed (q value < 0.05), 4 h post-MI. 24 h post-MI, 1,702 genes and 272 proteins were differentially expressed. Hierarchical clustering dendrogram of gene expression (Figures 1B and S1A) and protein levels (Figure 1C) shows that in both cases, the experimental groups are clustered together, demonstrating significant differences in the transcriptome and proteome post-MI. Pearson correlation analysis revealed a positive correlation between changes in gene expression and protein expression, both 4 h post-MI (r squared = 0.02; Figure 1D) and 24 h post-MI (r squared = 0.13; Figure 1E). Whereas we found a significant correlation between the changes in gene expression and protein levels in the ischemic heart post-MI and in our search for a 5′ UTR element that may elevate translation of modRNA in the heart post-MI, we were interested in identifying genes with a significant noncorrelation relationship between mRNA and protein expression. Therefore, we screened for genes that encode for proteins with elevated levels in 4 or 24 h post-MI (fold change > 2) and mRNA downregulated at 4 or 24 h post-MI (fold change < 0.64) and 5′ UTR shorter than 100 bp. We were able to identify 3 genes in 4 h (Figure 1F) and 18 genes in 24 h (Figure 1G) that displayed high protein expression accompanied with lower or unchanged mRNA level, 24 h post-MI, compared to sham hearts. We identify 5 genes (gelsolin [Gsn], pregnancy-zone protein [Pzp], Serpina, fructosamine 3 kinase [Fn3k], and carboxylesterase 1D [Ces1d], marked in yellow, in Figures 1F and 1G) that had the shortest 5′ UTR among those 19 genes (as FERM domain containing 5 [Frmd5] and Serpina 1b are present in both our 4- and 24-h screen results), with an upregulated protein expression not related to their mRNA expression post-MI. In addition, we validate Ces1d expression results using qPCR and western blot, showing similar mRNA and protein expression as evaluated by RNA sequencing (RNA-seq) and proteomic analysis (Figures S1B–S1D).

Figure 1

Characterizing the Ischemic Heart Transcriptome and Proteome

(A) Sham-operated heart or heart, 4 or 24 h post-MI was collected, and the ischemic area tissue (or equivalent area in sham hearts) was divided into two equal pieces. One-half of the ischemic heart was sequenced for transcriptomic analysis (n = 10 total, Sham n = 3, 4 h post-MI n = 3, or 24 h post-MI n = 4), whereas the other one-half of the ischemic heart was evaluated for protein level using mass spectrometry (n = 12 total, Sham n = 4, 4 h post-MI n = 4, or 24 h post-MI n = 4). (B and C) Hierarchical clustering dendrogram of 2,272 genes with corresponding mRNA level (B) or 2,397 protein intensity (C) in Sham, 4 h post-MI or 24 h post-MI hearts. (D and E) Correlation analysis between changes in levels of proteins and mRNA in the LV, 4 (D) or 24 (E) h post-MI. The bottom-right shaded rectangles include genes that show static or reduced mRNA levels post-MI, whereas their encoded protein levels increased with comparison to sham. (F and G) A list of genes that encode for proteins with elevated protein levels (fold change > 2), while showing mRNA downregulated (fold change < 0.64), 4 h (F) or 24 h (G) post-MI. Genes with yellow backgrounds have a 5′ UTR that is shorter than 100 base pairs (see Table S1 for the full sequences of the 5′ UTR of the 5 genes).

To evaluate the translational efficiency of their 5′ UTR, we designed a Luc modRNA construct carrying different 5′ UTRs taken from the selected genes (Figure 2A). We then compared the expression of 5 newly designed, different Luc modRNAs with different 5′ UTRs with our Luc modRNAs that carry artificial 5′ UTRs (36 nt, Luc-Control) and are commonly used for in vitro modRNA production (Figure 2B). Our in vitro screening showed that all Luc modRNAs carrying different 5′ UTRs allowed protein translation of modRNA in rat neonatal cardiomyocytes (RNCMs). However, in comparison to the Luc-Control, Gsn, Pzp, Serpina, and Fn3k had significantly lower translation, except Luc modRNA carrying 5′ UTR of Ces1d (Luc-Ces1d), which showed a significant increase of 23% in modRNA translation in comparison to Luc-Control (Figure 2C). To evaluate that this result is not due to different transfection efficiency, a parallel in vitro experiment was done with cotransfection of firefly Luc modRNA carrying different 5′ UTR modRNAs and renilla Luc modRNA as an internal control that carries control 5′ UTR in neonatal rat CMs. 24 h post-transfection, we used IVIS system to measure simultaneously the translation of the two different luciferase (firefly and renilla) modRNA in vitro (Figure S2). Our results show, similar to Figure 2, that Ces1d modRNA gave significantly higher firefly Luc modRNA translation with comparison to firefly Luc-Control modRNA, without significant change in renilla Luc modRNA translation. These results indicate that the high translation is due to the use of Ces1d 5′ UTR and not due to different transfection efficiency.

Figure 2

Evaluating the Translation Efficiency of modRNA Carrying Different Potential 5′ UTRs in Neonatal Rat CMs or in a Mouse Heart Ischemic Model

(A) Schematic representation of modRNA structure and the replacement of different potential 5′ UTRs in the Luc modRNA. (B) Experimental plan to evaluate the translation efficiency of Luc or GFP modRNA carrying different potential 5′ UTRs in neonatal rat CMs using IVIS or western blot analysis, respectively. (C and D) Quantification of the IVIS (C; n = 4) and western blot (D) experiments that were described in (B). (E) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying different potential 5′ UTRs in mouse hearts, 24 h post-MI using IVIS analysis. (F) Quantification of the IVIS experiment that was described in (E) (n = 4). One-way ANOVA and Tukey’s multiple comparison test were used for (C) and (F). ∗∗∗p < 0.001, ∗p < 0.05; N.S, not significant.

Additionally, to confirm our finding, we made EGFP modRNA, carrying 5′ UTR of Ces1d (EGFP-Ces1d), or EGFP modRNA, carrying artificial 5′ UTR (EGFP-Control) and compared their translation level of EGFP using western blot in neonatal rat CMs. Similar to our previous results, EGFP-Ces1d showed a 22% increase in modRNA translation in comparison to EGFP-Control (Figure 2D).

We decided to evaluate further our Luc modRNA carrying different 5′ UTRs in an MI mouse model. We selected the Luc-5′ UTR modRNAs (Pzp, Serpina 1b, and Ces1d) that showed the highest translation in our in vitro assay (Figures 2B–2D), and we measured their translation at 1, 2, and 3 days post-MI (Figure 2E). Our analysis revealed that Luc-Ces1d modRNA had a significantly (2-fold) higher expression in comparison to Luc-Control or Luc modRNA carrying 5′ UTR of Pzp or Serpina (Figure 2F).

To identify if Ces1d 5′ UTR regulates protein translation only in the heart ischemic mouse model, we evaluated both Luc-Ces1d and Luc-Control modRNA translation in a nonischemic heart (Figures 3A and 3B) or ischemic heart mouse models (Figures 3C, 3D, S3A, and S3B). We show that whereas there were no significant changes in the translation level of Luc-Ces1d in comparison with Luc-Control in the heart without MI, 1, 2, and 3 days post-injection, there is a significantly higher translation of Luc-Ces1d, 2 days post-MI in comparison to the Luc-Control. This may indicate that Ces1d enhances modRNA translation in the heart only under an ischemic condition, such as MI.

Figure 3

5′ UTR of Ces1d Enhances Significantly mRNA Translation in Ischemic Heart but Not in a Nonischemic Mouse Model

(A) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying 5′ UTR of Ces1d or artificial (control) 5′ UTR in a nonischemic heart model. (B) Quantification of the experiment that was described in (A) (n = 15). (C) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying 5′ UTR of Ces1d or artificial (control) 5′ UTR in an ischemic heart model. (D) Quantification of the experiment that was described in (C) (n = 15). Two-way ANOVA and Tukey’s multiple comparison test were used for (B) and (D). ∗p < 0.05; N.S, not significant.

We then wanted to evaluate the role of Ces1d in modRNA translation under an ischemic condition in other organs besides the heart. As acute hepatic or renal ischemia may lead to hepatic or renal failure that may be fatal, we decided to choose both liver and kidney as two representative organs for ischemic disease. Similar to the heart, ischemic injury in the liver significantly increases the expression of Luc-Ces1d (4 h post-delivery) in comparison to Luc-Control (Figures 4A–4D, S3C, and S3D). Again, similar to the heart, no significant differences were seen between groups in the liver without ischemic injury (Figures 4A and 4B). On the contrary, Luc-Ces1d had no significant translation differences in comparison to Luc-Control in the ischemic and nonischemic condition in the kidney (Figures 4E–4H, S3E, and S3F).

Figure 4

5′ UTR of Ces1d Enhances Significantly Luc modRNA Translation in Liver but Not in a Nonischemic or Kidney Ischemic Mouse Model

(A) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying 5′ UTR of Ces1d or artificial (control) 5′ UTR in a nonischemic liver model. (B) Quantification of the experiment that was described in (A) (n = 6). (C) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying 5′ UTR of Ces1d or artificial (control) 5′ UTR in an ischemic liver model. (D) Quantification of the experiment that was described in (C) (n = 6). (E) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying 5′ UTR of Ces1d or artificial (control) 5′ UTR in a nonischemic kidney model. (F) Quantification of the experiment that was described in (E) (n = 6). (G) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying 5′ UTR of Ces1d or artificial (control) 5′ UTR in an ischemic kidney model. (H) Quantification of the experiment that was described in (G) (n = 6). Two-way ANOVA and Tukey’s multiple comparison test were used for (B), (D), (F), and (H). ∗p < 0.05; N.S, not significant.

To identify the RNA element in Ces1d 5′ UTR that is responsible for the significant enhancement of translation of modRNA carrying Ces1d 5′ UTR, we examined the Ces1d 5′ UTR consensus elements that had been conserved among different species (e.g., mouse, rat, pig, gerbil, and human). Interestingly, 4 out of the 5 elements (elements B, C, D, and E) were conserved among species (Figure 5A). Based on this information, we designed a Luc modRNA construct carrying different 5′ UTR elements of Ces1d (elements A–E) and compared its translation ability in neonatal rat CM to Luc-Ces1d (Figures 5B and 5C). Our results indicated a significant reduction in translation ability in elements A, B, C, and E but not in D (Figure 5C). We hypothesized that element D is the RNA element that is responsible for the higher translation ability of Luc-Ces1d. To test our hypothesis, we compared in an ischemic heart model the expression of Luc-Ces1d or Luc-element D with Luc-Control over 3 days (Figure 5D). Our results show that Luc modRNA, carrying 5′ UTR of Ces1d or element D, has higher translation in comparison to Luc-Control (Figure 5E). Intriguingly, element D alone had a significantly higher translation in the heart at day 3 post-MI (Figure 5D). Overall, with the combination of the 3-day readout, Luc-element D has a 2.5-fold higher translation in the heart post-MI in comparison to the widely used artificial 5′ UTR (Luc-Control). However, Luc-element D failed to increase translation over Luc-Ces1d in liver ischemic injury (Figures S4A and S4B) and over Luc-Ces1d or Luc-Control in heart nonischemic injury (Figures S5A and S5B).

Figure 5

A Specific RNA Element in 5′ UTR of Ces1d Enhances Significantly Luc modRNA Translation in a Heart Ischemic Mouse Model

(A) List of the different RNA elements in 5′ UTR of Ces1d that has been conserved across different species. (B) Experimental plan to evaluate the translation efficiency of different RNA elements in 5′ UTR of Ces1d in neonatal rat CMs using Luc modRNA and IVIS analysis. (C) Quantification of the experiments that were described in (B) (n = 4). (D) Experimental plan to evaluate the translation efficiency of Luc modRNA carrying full-length 5′ UTR of Ces1d, only element D from Ces1d 5′ UTR, or artificial (control) 5′ UTR in an ischemic heart model. (E) Quantification of the experiments that were described in (D) (n = 5). One-way ANOVA and Tukey’s multiple comparison test were used for (C). Two-way ANOVA and Tukey’s multiple comparison test were used for (E). ∗∗∗p < 0.001, ∗∗p < 0.01, ∗p < 0.05; N.S, not significant.

Discussion

The use of modRNA as a gene-delivery tool is growing in the field of therapeutic medicine. Our laboratory aims to use modRNA to induce cardiac protection and cardiovascular or cardiac regeneration post-myocardial infarction. Zangi et al. have shown that immediate delivery of vascular endothelial growth factor (VEGF)-A modRNA in a MI mouse model led to induction of cardiovascular regeneration. A follow-up study in large animals shows a significant improvement of cardiac function when VEGF-A modRNA was delivered 1 week post-MI. Recently, it was shown that intradermal delivery of VEGF-A modRNA in patients suffering from type 2 diabetes is safe and may promote angiogenesis. VEGF-A modRNA has now been evaluated in phase 2a human clinical trials to improve cardiac function in heart-failure patients. In parallel, other groups are using the modRNA platform in preclinical studies of different liver diseases to deliver different target genes. The different liver disease models that have been used are for factor IX deficiency hemophilia B, acute intermittent porphyria, glycogen storage disease type 1A, thrombotic thrombocytopenic purpura, alpha-1 antitrypsin deficiency, Crigler-Najjar syndrome type 1, and urea cycle disorder with target genes of FIX, PBGD, glucose-6-phosphatase, ADAMTS13, SERPINA1, bilirubin-UGT, and ornithine transcarbamoylase, respectively. As was described above, both the heart and liver disease have been heavily studied in search for new treatments. One of the obstacles in moving to large animals and clinical trials is the need for a large amount of modRNA to transfect human or large-animal heart and liver. In addition, due to the short expression of modRNA, there might be a need to deliver the modRNA several times. In order to reduce the need for a large amount of modRNA for delivery, we aim to improve modRNA translation so the amount of protein that has been received from the same amount of modRNA will be higher.

Several elements within RNA are responsible for the regulation of translation, including type of nucleotides, poly(A) tail length, or the structure of 5′ UTR, 3′ UTR, and cap analog. We showed that replacement of pseudouridine with N1-Methyl-Pseudouridine-5'-Triphosphate (1-M-pseudouridine) results in a higher modRNA translation in the heart of the latter. In addition we found that a longer poly(A) tail can increase the translation of the modRNA in vivo (unpublished data). Here, we describe a unique first-of-this-kind screen, comparing proteomic and transcriptomic analysis to identify the 5′ UTR sequence that can increase the translation of modRNA in an ischemic condition with comparison to the widely used artificial 5′ UTR that has been used in previous modRNA publications.4, 5, 6,,,

Our data show that there is a positive correlation at both 4 and 24 h post-MI in mRNA levels and protein intensity in the ischemic heart. However, we were able to identify 19 different genes with negative correlation, in which their mRNA level was reduced or unchanged, whereas their protein levels were upregulated, 4 and 24 h post-MI. The average length of 5′ UTRs is ∼100 to ∼220 nt across species. In vertebrates, longer 5′ UTRs tend to be associated with poor translation. Therefore, we decided to select 5 genes with the shortest 5′ UTR (<100 base pairs) (Figure 1).

Negative correlation between mRNA levels and protein expression has been reported, especially upon internal or external stimuli that trigger alteration in translation of specific genes. VEGF has been shown to be a stress-induced protein in many conditions, such as hypoxia and hypoglycemia., Other genes in which translation changes in response to stimulation of internal or external signals are platelet-derived growth factor 2 (PDGF2) and transforming growth factor β (TGF-β). During the embryonic stage, when most of the organ development and cell differentiation takes place, translation regulation plays a significant role by altering the levels of expression of specific mRNA subsets during a certain time frame, whereas the bulk of transcripts remain unaffected.54, 55, 56

In our study, we were able to identify several 5′ UTRs that may allow translation of modRNA in the heart or liver post-ischemic injury. Both Pzp and Serpina 1b 5′ UTR show nonsignificant and similar translation ability as the control 5′ UTR sequence in vivo (Figures 2E and 2F). As our control 5′ UTR (artificial 5′ UTR) is the top and most known and commonly used 5′ UTR in the mRNA field,,,57, 58, 59, 60, 61 we believe that these results indicate how well, relative to this premium control, our selected 5′ UTR sequences performed.

In this particular work, we focused on the 5′ UTR of Ces1d (Figures 2, 3, and 4) as an enhancer of modRNA translation in cardiac and hepatic ischemic conditions. Ces1d belongs to a family of carboxylesterases that have important roles in lipid metabolism, and they hydrolyze endogenous esters and thioesters. Carboxylesterases are known for involvement in the detoxification of environmental toxicants, as well as the prodrug metabolism. Ces1d is the functional mouse ortholog of human CES1 and has similar protein-expression profiles at different cells or tissues as CES1. Several roles of Ces1d have been reported to be directly associated with lipid metabolism.62, 63, 64 As lipid metabolism is an important process for normal heart function, fatty acids are the preferred substrates under aerobic conditions. As Ces1d association with lipid metabolism and as MI lead to alteration of the lipid metabolism, we hypothesize that Ces1d mRNA has been triggered by the ischemic condition in the heart post-MI, which leads it to translate better. We showed that element D is the RNA element in Ces1d that is responsible for the elevation in mRNA translation post-MI (Figure 4). It will be interesting to evaluate the element D and Ces1d in different ischemic conditions and different organs bedsides the heart, liver, and kidney. The fact that Ces1d gave higher translation in the ischemic heart and liver but not the kidney is interesting, as all three organs primarily use fatty acid oxidation for energy. This may indicate that each organ and physiological condition will need a separate evaluation using a similar approach as ours.

In conclusion, we have identified the 5′ UTR of Ces1d and RNA element D in Ces1d 5′ UTR as RNA elements for improving modRNA translation in the heart and liver post-ischemic injury. This may have clinical applications, as both organs have been heavily targeted with modRNA in different cardiac and hepatic diseases. Our results may contribute to designing superior modRNA for preclinical studies in ischemic cardiac and liver diseases, carrying 5′ UTR of Ces1d or RNA element D.

Methods

Mice

All animal procedures were performed according to protocols approved by the Icahn School of Medicine at Mount Sinai Institutional Care and Use Committee. Swiss Webster (CFW) mice were used for the study. Before surgery, mice were anesthetized with a cocktail of 100 mg/kg ketamine and 10 mg/kg xylazine. For protein expression, mice were injected with 25 μg of modRNA in citrate buffer directly into the myocardium during open-chest surgery. When required, 25 μg of modRNA was injected into the border zone with three injections immediately after left anterior descending artery (LAD) ligation.

modRNA Synthesis

Clean PCR products generated with plasmid templates (GeneArt; Thermo Fisher Scientific) were used as the template for mRNA. modRNAs were generated by transcription in vitro with a customized ribonucleoside blend of anti-reverse cap analog, 3′-O-Me-m7G(5′)ppp(5′)G (6 mM; TriLink Biotechnologies), guanosine triphosphate (1.5 mM; Life Technologies), adenosine triphosphate (7.5 mM; Life Technologies), cytidine triphosphate (7.5 mM; Life Technologies), and N1-Methyl-Pseudouridine-5'-Triphosphate (7.5 mM; TriLink Biotechnologies). The mRNA was purified with the MEGAclear kit (Life Technologies) and treated with Antarctic Phosphatase (New England Biolabs). It was then repurified with the MEGAclear kit. The mRNA was quantified on a NanoDrop spectrometer (Thermo Scientific), precipitated with ethanol and ammonium acetate, and resuspended in 10 mM Tris-HCl and 1 mM EDTA.

Rat Neonatal CM Isolation

Ventricular RNCMs were isolated from 3- to 4-day-old Sprague-Dawley rat pups by a Pierce primary cardiomyocytes isolation kit (Thermo Fisher Scientific; catalog number 88282). After isolation, cells were incubated in 10% horse serum DMEM and change media, and a cardiomyocyte growth supplement was added after 16 h and transfected with modRNA.

In Vitro modRNA Transfection

2.5 μg/well of a 24-well plate of Luc modRNA or 10 μg/well of a 6-well plate of nuclear GFP (nGFP) modRNA was transfected into neonatal rat CMs using the transfection reagent jetPEI (Polyplus).The transfection mixture was prepared according to the manufacturer’s protocol, and then it was added to cells cultured in DMEM containing 10% fetal bovine serum (FBS) and antibiotic-antimycotic (anti-anti). Then, 24 h post-transfection, the expression level of cells was imaged and measured in IVIS, or cell lysates were collected and analyzed by western blot.

Heart Ischemic Injury

MI was induced by permanent ligation of the LAD. The left thoracic region was shaved and sterilized. After intubation, the heart was exposed by left thoracotomy. The LAD was ligated with a suture. Mouse hearts without (sham) or with ischemic injury (MI) were collected at 4 and 24 h post-MI (Figure 1A). The ischemic area tissue (or equivalent area in sham hearts) was collected and divided into two pieces and snap frozen quickly. One-half of the ischemic heart tissues was sent for RNA-seq, whereas the other one-half was sent for proteomics analysis to evaluate gene and protein fold change, respectively, between MI and sham hearts. When required, we injected 25 μg modRNA into the infarct border zone immediately after LAD ligation. The thoracotomy and skin were sutured closed in layers. Excess air was removed from the thoracic cavity, and the mouse was removed from ventilation when normal breathing was established.

Liver Ischemic Injury

Liver ischemia was induced by closing the left lateral lobe and median for 1 h, and 25 μg modRNA was injected into the left lateral lobe in three different injections immediately after clip was removed.

Kidney Ischemic Injury

Kidney ischemia was induced by applying a micro clip onto the renal artery and renal vein. Successful ischemia can be visually confirmed by a gradual uniform darkening of the kidney. The clip was removed after 30 min, and 25 μg of modRNA was injected in the kidney with three injections.

In Vivo modRNA Delivery

Luc modRNA (25 μg) in a total volume of 60 μL in Tris-borate (TB) buffer was delivered via direct injection to the myocardium. The sucrose-citrate buffer contains 20 μL sucrose in nuclease-free water (0.3 g/mL) and 20 μL citrate (0.1 M, pH 7; Sigma) mixed with 20 μL modRNA. The transfection mixture was directly injected (three individual injections, 20 μL each) into the heart, kidney, or liver.

Detection of Luciferase Expression Using the IVIS System

Bioluminescence imaging of the transfected cells (24–72 h) or injected mice was taken at different time points (4–144 h) in the IVIS system. To visualize cells expressing firefly luciferase in vitro, D-luciferin was added to the cell-culture plate, and an image was taken in the IVIS system (IVIS Spectrum National Center for Research Resources [NCRR] S10-RR026561-01 at the Preclinical Small Imaging Core at Mount Sinai Medical Center). To visualize cells expressing renilla luciferase in vitro, cells were washed twice with media, renilla luciferase substrate was added to the cell-culture plate, and an image was taken in emission filter 500. To visualize tissues expressing Luc in vivo, mice were anesthetized with isoflurane (Abbott Laboratories), and luciferin (150 μg/g body weight; Sigma) was injected intraperitoneally. Mice were imaged using an IVIS imaging system (IVIS Spectrum NCRR S10-RR026561-01 at the Preclinical Small Imaging Core at Mount Sinai Medical Center) every 2 min until the Luc signal reached a plateau. Imaging data were analyzed and quantified with Living Image software.

Western Blotting

Cell lysates were collected and subjected to SDS-PAGE in 12% precast NuPAGE Bis/Tris gels (Invitrogen) under reducing conditions in MES running buffer (Invitrogen). The resulting bands were transferred onto a nitrocellulose membrane (Bio-Rad) by blotting in a semidry transfer apparatus with NuPAGE-MOPS (3-(N-morpholino)propanesulfonic acid) transfer buffer (Invitrogen). The membrane was blocked by incubation with Tris-buffered saline (TBS)/Tween containing 5% dry milk powder and incubated with specific primary antibodies overnight at 4°C. It was then washed in TBS/Tween and incubated with rabbit or goat secondary antibodies conjugated to horseradish peroxidase for 1 h at room temperature. Antibody binding was detected with an enhanced chemiluminescence (ECL) detection system (Pierce). We used prestained protein standards (Amersham) to determine molecular weight.

RNA Isolation

Total RNA was isolated with an RNeasy Mini Kit (QIAGEN), and DNA was degraded by treatment with TURBO DNase (Invitrogen). RNA quality was checked by a bioanalyzer.

RNA Sequencing

Poly(A)-tailed RNA was prepared by the Epigenomics Core at Cornell Medical College, with the mRNA Seq Sample Prep Kit (Illumina) and used to create libraries for HiSeq 2000 sequencing (Illumina). We used single, 50-bp reads for sequencing. We obtained a mean of 30 million reads per sample, with a mean quality score of 35.2. We used Partek flow software for data analysis. RNA-seq reads were aligned to mm10 with STAR version 2.53a. Read counts were generated by the application of the Partek expectation-maximization (E/M) algorithm to University of California, Santa Cruz (UCSC), RefSeq 2017-08-02. Counts were normalized with trimmed mean of maximum (TMM)-values algorithm, and the Partek flow gravitational search algorithm (GSA) was used for statistical analysis. The RNA-seq data used in this study is available using GenBank: GSE138201.

Protein Mass Spectrometry

All solvents were high-performance liquid chromatography (HPLC) grade from Sigma-Aldrich, and all chemicals, where not stated otherwise, were obtained from Sigma-Aldrich. For ample preparation, samples were lysed in Biognosys’ lysis buffer with a TissueLyser II bead mill (QIAGEN) using stainless-steel grinding beads for 3 min at 30 Hz. Samples were treated with benzonase after lysis to reduce DNA contamination. Protein concentrations of the lysates were estimated using a UV-VIS (visible) spectrometer (SPECTROstar Nano, BMG Labtech). Approximately 100 μg of protein from each sample was denatured using Biognosys’ denature buffer, reduced using Biognosys’ reduction solution for 60 min at 37°C and alkylated using Biognosys’ alkylation solution for 30 min at room temperature in the dark. Subsequently, digestion to peptides was carried out using 3 μg of trypsin (Promega) overnight at 37°C.

Peptides were desalted using C18 MacroSpin columns (The Nest Group), according to the manufacturer’s instructions, and dried down using a SpeedVac system. Peptides were resuspended in 50 μL liquid chromatography (LC) solvent A (1% acetonitrile, 0.1% formic acid [FA]) and spiked with Biognosys’ iRT (indexed retention time) kit calibration peptides. Peptide concentrations were determined using a UV-VIS spectrometer (SPECTROstar Nano, BMG Labtech).

For high pH reverse-phase (HPRP) fractionation, equal volumes of samples from each condition were pooled. Ammonium hydroxide was added to all pools to a pH value > 10. The fractionation was performed using a Dionex UltiMate 3000 RS Pump (Thermo Scientific) on an Acquity UPLC CSH (ultra-performance liquid chromatography charged surface hybrid) C18 1.7 μm, 2.1 × 150 mm column (Waters). The gradient was 2% to 35% solvent B in 10 min; solvents were the following: A, 20 mM ammonium formate in H2O, and B, acetonitrile. Fractions were taken every 15 s and sequentially pooled to 4 fraction pools. These were dried down and resolved in 20 μL solvent A. Prior to mass spectrometric analyses, they were spiked with Biognosys’ iRT kit calibration peptides. Peptide concentrations were determined using a UV-VIS spectrometer (SPECTROstar Nano, BMG Labtech).

For the liquid chromatography-tandem mass spectrometry (LC-MS/MS; shotgun) measurements, 2 μg of peptides was injected to an in-house packed C18 column (ReproSil-Pur, Dr. Maisch GmbH; 1.9 μm particle size, 120 Å pore size, 75 μm inner diameter, 50 cm length; New Objective) on a Thermo Scientific EASY-nLC 1200 nano-liquid chromatography system connected to a Thermo Scientific Q Exactive HF mass spectrometer equipped with a standard nano-electrospray source. LC solvents were the following: A, 1% acetonitrile in water with 0.1% FA, and B, 15% water in acetonitrile with 0.1% FA. The nonlinear LC gradient was 1%–55% solvent B in 60 min, followed by 55%–90% B in 10 s, 90% B for 10 min, 90%−100% B in 10 s, and 1% B for 5 min. A modified top 15 method was used.

The mass spectrometric data were analyzed using Biognosys’ search engine Pulsar (version 1.0.16105), and the false-discovery rate on peptide and protein level was set to 1%. A mouse UniProt FASTA database (Mus musculus, 2017-07-01) was used for the search engine, allowing for 2 missed cleavages and variable modifications (N-terminal acetylation, methionine oxidation).

For the LC-MS/MS heart-rate monitor (HRM) measurements, 2 μg of peptides per sample was injected to the C18 column (ReproSil-Pur, Dr. Maisch GmbH; 1.9 μm particle size, 120 Å pore size, 75 μm inner diameter, 50 cm length; New Objective) on a Thermo Scientific EASY-nLC 1200 nano-liquid chromatography system connected to a Thermo Scientific Q Exactive HF mass spectrometer equipped with a standard nano-electrospray source. LC solvents were the following: A, 1% acetonitrile in water with 0.1% formic acid, and B, 15% water in acetonitrile with 0.1% formic acid. The nonlinear LC gradient was 1%–52% solvent B in 60 min, followed by 52%−90% B in 10 s, and 90% B for 10 min. A data-independent acquisition (DIA) method with one full-range survey scan and 14 DIA windows was used.

HRM mass spectrometric data were analyzed using Spectronaut Pulsar software (Biognosys). The false-discovery rate on peptide and protein level was set to 1%, and data were filtered using row-based extraction. The assay library (protein inventory) generated in this project was used for the analysis. The HRM measurements analyzed with Spectronaut were normalized using local regression normalization.

Merging Proteomics and Transcriptomics Data

For proteomics/transcriptomics merging, gene IDs contained in the transcriptomics datasets were matched to gene names in the mouse UniProt Swiss-Prot proteome.

Statistical Analysis for Luciferase Activity Assay and RNA Protein Correlation

All statistical analysis was performed with GraphPad Prism software. Values are reported as mean ± SD. One-way ANOVA with Bonferroni correction (∗p < 0.05 considered significant) was used for comparisons among groups. Pearson’s R correlation was calculated for correlation between changes in mRNA expression and protein levels.

Author Contributions

N.S. designed the study, carried out most of the experiments, analyzed most of the data, and wrote the manuscript. Y.H. designed, carried out, and analyzed the RNA-seq and proteomics experiments and revised the manuscript. M.T.K.S. prepared modRNAs and performed experiments. K.K. performed experiments and revised the manuscript. A.M. revised the manuscript. N.H., B.A., and S.A. performed experiments and helped analyze the data. E.C. carried out all surgery and modRNA injections in mouse models. L.Z. designed the experiments, analyzed data, and wrote the manuscript.

Conflicts of Interest

L.Z., N.S., and Y.H. are inventors of a provisional patent, which covers the results in this manuscript.