Biochemical characterization of mRNA capping enzyme from Faustovirus
- S. Hong Chan,
- Christa N. Molé,
- Dillon Nye,
- Lili Mitchell,
- Nan Dai,
- Jackson Buss,
- Daniel W. Kneller,
- Joseph M. Whipple and
- G. Brett Robb
- Corresponding author: chan{at}neb.com
Abstract
The mammalian mRNA 5′ cap structures play important roles in cellular processes such as nuclear export, efficient translation, and evading cellular innate immune surveillance and regulating 5′-mediated mRNA turnover. Hence, installation of the proper 5′ cap is crucial in therapeutic applications of synthetic mRNA. The core 5′ cap structure, Cap-0, is generated by three sequential enzymatic activities: RNA 5′ triphosphatase, RNA guanylyltransferase, and cap N7-guanine methyltransferase. Vaccinia virus RNA capping enzyme (VCE) is a heterodimeric enzyme that has been widely used in synthetic mRNA research and manufacturing. The large subunit of VCE D1R exhibits a modular structure where each of the three structural domains possesses one of the three enzyme activities, whereas the small subunit D12L is required to activate the N7-guanine methyltransferase activity. Here, we report the characterization of a single-subunit RNA capping enzyme from an amoeba giant virus. Faustovirus RNA capping enzyme (FCE) exhibits a modular array of catalytic domains in common with VCE and is highly efficient in generating the Cap-0 structure without an activation subunit. Phylogenetic analysis suggests that FCE and VCE are descended from a common ancestral capping enzyme. We found that compared to VCE, FCE exhibits higher specific activity, higher activity toward RNA containing secondary structures and a free 5′ end, and a broader temperature range, properties favorable for synthetic mRNA manufacturing workflows.
Keywords
INTRODUCTION
The rapid and global deployment of synthetic mRNA vaccines against SARS-CoV-2 was the result of decades of research on several fronts (Dolgin 2021). The core technologies of synthetic mRNA are in vitro transcription (IVT) using T7 RNA polymerase (Davanloo et al. 1984) and 5′ cap incorporation using vaccinia virus RNA capping enzyme (VCE) (Martin and Moss 1976; Shuman and Moss 1990; Fuchs et al. 2016) or cap analogs (Ishikawa et al. 2009; Henderson et al. 2021). The 5′ cap structures perform important biological functions such as efficient translation of the mRNA, evasion of the innate immune surveillance systems, and 5′-mediated mRNA turnover (Ramanathan et al. 2016). Therefore, 5′ cap incorporation is considered one of the critical quality attributes (CQAs) in mRNA manufacturing (Daniel et al. 2022). Methods to assess the extent of 5′ cap incorporation have been published (Cairns et al. 2003; Beverly et al. 2016; Chan et al. 2022; Vlatkovic et al. 2022).
The mRNA 5′ cap structures were first identified by multiple research groups in the 1970s (Desrosiers et al. 1975; Furuichi et al. 1975; Moyer et al. 1975; Wei et al. 1975; Furuichi 2015). In lower eukaryotes such as fungi, mRNA contains a 5′ cap structure (Cap-0) where a N7-methyl-guanosine is linked to the 5′ end of the RNA via a 5′-to-5′ triphosphate linkage (Fig. 1A). In metazoans, the Cap-0 structure is further modified to form the Cap-1 structure where a methyl group is added to the 2′-O group of the first nucleotide of the RNA (Barbosa and Moss 1978). The Cap-1 structure has been shown to be an important factor in evading the innate immune surveillance systems in humans (Goubau et al. 2014; Hyde and Diamond 2015). The Cap-1 structure can be further modified to the Cap-2 structure, where the 2′-O groups of the first and second nucleotides of the RNA are methylated (Fig. 1A; Werner et al. 2011; Smietanski et al. 2014; Drazkowska et al. 2022). The Cap-2 structure has been implicated in further evasion of innate immune response over the lifetime of mRNA molecules (Despic and Jaffrey 2023).
(A) Enzymatic reactions involved in RNA 5′ capping. First, the RNA triphosphatase (TPase) hydrolyzes the 5′ triphosphate group (pppN−) to form a 5′ diphosphate (ppN−) (Reaction 1). Then, an RNA 5′ guanylyltransferase adds a GMP group to the 5′ diphosphate to form an unmethylated G-capped structure (GpppN−) through a reversible mechanism (Reaction 2). This is followed by an RNA cap guanosine N7-methyltransferase activity that adds a methyl group to the N7 position of the G cap and forms the Cap-0 structure (m7GpppN−) (Reaction 3). In metazoans, a fourth enzyme RNA cap 2′-O-methyltransferase (2′-OMTase) further adds a methyl group to the 2′-O position of the initiating nucleotide of the transcript, forming the Cap-1 structure (m7GpppNm−) (Reaction 4). (B) Catalytic domains of FCE. (Top panels) Individual domains of the AlphaFold2 prediction of FCE were aligned to published atomic structures. Insets show the alignment of highly conserved residues labeled with FCE residue numbering. (Bottom panels) Amino acid residues of the motifs were aligned with indicated enzymes. Cet1, Ceg1, and Abd1 are the TPase, GTase, and N7MTase of Saccharomyces cerevisiae, respectively. hRNMT is the human RNA cap N7MTase. MCE refers to Acanthamoeba polyphaga mimivirus capping enzyme. Conserved residues are highlighted in gray. FCE residues that are mutated in this study and the corresponding residues in orthologous enzymes are highlighted in black. (i) The predicted structure of FCE's TPase domain is aligned to the VCE TPase domain-ppRNA complex extracted from the CryoEM cotranscriptional capping complex (PDB ID 6RIE). The inset shows the alignment of the conserved ELE–ELE motif. (ii) The predicted structure of FCE's GTase domain is aligned to the X-ray crystal structure of Paramecium bursaria chlorella virus (PBCV-1) GTase with covalent GMP intermediate (PDB ID 1CKN). Inset shows the phosphamide covalent linkage between the catalytic Lys of PBCV-1 GTase and GMP. (iii) The predicted FCE N7MTase domain structure is aligned to the X-ray crystal structure of the N7MTase domain of VCE D1 subunit in complex with co-product SAH (PDB ID 4CKB). Inset shows the conserved SAM-binding site.
The first mRNA capping enzyme was purified from vaccinia virus particles (Martin and Moss 1975; Martin et al. 1975; Shuman et al. 1980) around the same time mRNA cap structures were discovered. As the first mRNA capping enzyme available in recombinant form (Shuman 1990; Shuman and Moss 1990), VCE became the reference for subsequent research on cellular and viral RNA capping enzymes (Ghosh and Lima 2010; Shuman 2015; De Vlugt et al. 2018; Park et al. 2022). For example, RNA capping enzymes in general refer to proteins that generate the Cap-0 structure through three enzyme activities that VCE possesses: the RNA triphosphatase activity (TPase), the RNA guanylyltransferase activity (GTase), and the RNA cap N7-guanosine methyltransferase activity (N7MTase) (Fig. 1A, Reactions 1–3). The TPase domain consists of the triphosphate tunnel metalloenzyme (TTM) fold, where the 5′ triphosphate RNA is inserted into a long β-barrel structure for hydrolysis (Kyrieleis et al. 2014). The GTase domain shares the adenosyl transferase active site of DNA and RNA ligases that involves a covalent Lys-GMP intermediate and has the propensity to catalyze the forward and reverse reactions (Soulière et al. 2008). The RNA cap N7MTase domain exhibits the characteristic core fold of the class I MTase family (De la Peña et al. 2007; Kyrieleis et al. 2014).
The large subunit of VCE (D1R) adopts a modular structure where each of the three structural domains carries one of the component enzyme activities (Shuman and Morham 1990; Kyrieleis et al. 2014). These functional domains are largely conserved in eukaryotes, but their organization can be variable. For example, three individual proteins, Cet1, Ceg1, and Abd1, carry each of the three activities in S. cerevisiae (Shuman et al. 1994; Mao et al. 1995; Lima et al. 1999). In mammals, RNMT possesses the N7MTase activity, and RNGTT possesses the TPase and GTase activities where the TTM domain is replaced by a cysteinyl-phosphatase domain (Wen et al. 1998). In metazoans, the Cap-0 structure is further modified to the Cap-1 structure by a 2′OMTase (Fig. 1A, Reaction 4). The enzymes CMTR1 and VP39 possess the 2′OMTase activity in humans and vaccinia virus, respectively. RNA viruses have evolved a diverse set of mechanisms for mRNA capping: formation of an enzyme–pRNA covalent intermediate in vesicular stomatitis virus (Ogino and Banerjee 2007) and coronavirus (Park et al. 2022) through the kinase-like Nidovirus RdRp-associated nucleotidyltransferase domain (NiRAN); direct transfer of m7GMP to pp-RNA in Alphavirus (Ahola et al. 1995; Tomar et al. 2011), and cap-snatching from host mRNA in influenza virus (Shaw and Lamb 1984; Reich et al. 2014) and the yeast Totivirus (Fujimura and Esteban 2012). Viral RNA capping mechanisms have been exploited for antiviral drug development. For example, cap-snatching inhibitors are an effective class of anti-influenza drugs (Hsu et al. 2012; Jones et al. 2016; García-Sastre 2019; Todd et al. 2021).
Vaccinia virus belongs to a monophyletic class of viruses called the nucleocytoplasmic large DNA viruses (NCLDVs), most of which encode their own capping enzymes. In this communication, we characterize a single-subunit mRNA capping enzyme from another NCLDV, Faustovirus, that infects amoeba. Unlike VCE, a heterodimeric enzyme that requires a small subunit to generate Cap-0 efficiently (Mao and Shuman 1994; Kyrieleis et al. 2014), the single-chain Faustovirus mRNA capping enzyme (FCE) generates Cap-0 from 5′ triphosphate RNA efficiently. Similar to mimivirus, an NCLDV that also infects amoeba, Faustovirus does not seem to encode a small RNA capping enzyme subunit (Benarroch et al. 2008). FCE and the VCE large subunit (D1R) share a common organization of the catalytic domains. The two capping enzymes exhibit the same enzyme activities under similar optimal reaction conditions. Interestingly, our initial biochemical characterization reveals that FCE exhibits a higher specific activity in Cap-0 synthesis, operates under a broader active temperature range, and is more active on structured RNA.
RESULTS
Identification and predicted features
Faustovirus mRNA capping enzyme (FCE; GenBank Accession AMN83561.1; 873 aa) is a single-subunit RNA enzyme from Faustovirus, a giant virus that infects vermamoeba. It shares a very low sequence identity with the large subunit of VCE (D1R; GenBank Accession AAO89385.1; 844 aa) (14.8%) and has no BLAST hits within poxviruses (data not shown), a double-stranded DNA virus family that vaccinia virus belongs to. FCE also exhibits very low sequence identity (13.8%) to another single-subunit RNA capping enzyme from A. polyphaga mimivirus (GenBank Accession AAV50651.1; 1170 aa). However, manual inspection of the amino acids sequence showed that FCE possesses the conserved catalytic amino acids residues in capping enzymes from organisms such as vaccinia virus, S. cerevisiae, PBCV-1, mimivirus, and human (Fig. 1B). AlphaFold2 prediction suggests that FCE is organized into three domains which are highly similar to the experimentally determined structures of corresponding domains of VCE D1R, the PBCV-1 GTase, and the mimivirus capping enzyme TPase domain (Fig. 1B; Supplemental Fig. 1; Håkansson et al. 1997; Benarroch et al. 2008; Kyrieleis et al. 2014; Hillen et al. 2019).
Enzymatic activities
To verify that FCE generates the Cap-0 structure, we used purified FCE under standard VCE capping reaction conditions to cap a 25-nt RNA oligonucleotide containing a triphosphate group at the 5′ end. The reactions were then analyzed by nucleoside analysis using LC–MS. The results showed that FCE generates the same Cap-0 structure on RNA oligonucleotides with an adenosine or a guanosine as the initiating nucleotide (Supplemental Fig. 2A).
We next verified the presence of the three individual enzyme activities. Using capillary electrophoresis as readout (Chan et al. 2022), we subjected FCE to conditions that facilitated the TPase activity, TPase and GTase activities, or all three activities. Using substrate RNA oligonucleotides containing a 5′ triphosphate group (ppp−), FCE generated a 5′ diphosphate end (pp-) in the absence of GTP and SAM, unmethylated G-capped 5′ end (Gppp−) in the presence of GTP only, and Cap-0 N7-methyl G-capped 5′ end (m7Gppp−) in the presence of both GTP and SAM (Supplemental Fig. 2B). A time course of the FCE capping reaction in the presence of GTP and SAM showed that the conversion of 5′ ppp− to pp−, Gppp−, and then m7Gppp− mirrored that of VCE reactions (Fig. 2).
FCE and VCE generate the m7Gppp− (Cap-0) structure through the same pathway. Under identical conditions (1× RNA capping buffer, 0.5 mM GTP, 0.1 mM SAM, 1 mM DTT, and 5 nM enzyme at 25°C), FCE and VCE underwent the same reaction pathway to generate the m7Gppp− cap from 5′ triphosphate. As the concentration of 5′ ppp-RNA (blue circles) decreased over time, the concentration of the TPase product pp− (orange circles) increased initially and decreased over time, followed by a similar trend with a smaller amplitude for the GTase product Gppp− (gray circles) and an increase in N7MTase product m7Gppp− (yellow circles).
To verify the active sites, we mutated FCE residues that correspond to important active site residues of VCE and compared the mutants to wild-type (WT) FCE activity (Fig. 1B, bottom panels). Assays were designed to interrogate each of the three enzyme activities quantitatively using substrate RNA containing a 5′ triphosphate, diphosphate, or unmethylated G-cap group. Results for WT FCE are shown in Figure 3A. The FCE TPase domain hydrolyzed over 90% of 5′ ppp-RNA to pp-RNA in Reaction 1 (in the absence of GTP and SAM). The GTase activity guanylylated ∼90% of 5′ pp-RNA to the Gppp− capped form in Reaction 2 (with GTP but no SAM). The FCE N7MTase activity methylated ∼80% of Gppp-RNA to m7Gppp− in Reaction 3 (in the presence of SAM without GTP). The no enzyme controls showed that the ppp-RNA preparation contained a small amount of pp-RNA, the pp-RNA substrate preparation contained a small amount of ppp-RNA, and the Gppp-RNA preparation contained a small amount of pp-RNA. These impurities were side products in the manufacturing of the oligonucleotide RNA molecules.
Catalytic amino acid residues of FCE. (A) Under standard conditions (see Materials and Methods), 2.5 nM of FCE WT enzyme converted ≥90% of ppp-RNA into pp-RNA in the absence of GTP and SAM (Reaction 1), ∼90% pp− into Gppp− in the presence of GTP (Reaction 2) and ∼80% of Gppp− into m7Gppp− in the presence of SAM (Reaction 3). (B) Mutations of the Glu residues of the ELE motifs knocked out TPase activity but did not affect GTase or N7MTase activities. Mutants E26A/E28A, E184A/E186A, and E26A/E28A/E184A/E186A exhibited the same trend. (C) Mutation of Lys282 of the conserved KTDG motif to Met eliminated the formation of the GMP-GTase covalent intermediate (left panel). The mutant K282M exhibited undetectable GTase activity (Reaction 2) but exhibited high level of TPase (Reaction 1) and N7MTase activity (Reaction 3). (D) Mutation of Tyr712 in the conserved SAM-binding motif to Ala significantly decreased the N7MTase activity (Reaction 3) compared to the WT FCE enzyme but not TPase (Reaction 1) or GTase activity (Reaction 2).
TPase
Alanine mutations were introduced to the two conserved metal-coordinating glutamic acids in the TPase domain ELE motifs (Fig. 1B, left bottom panel). Double mutants E26A/E28A and E184A/E186A and quadruple mutant E26A/E28A/E184A/E186A were generated and tested for the presence of the three enzyme activities. All three mutants exhibited a similar activity profile: the mutants were inactive in Reaction 1 but displayed activity similar to the WT enzyme in Reactions 2 and 3 (Fig. 3B), showing that FCE contains the conserved ELE motifs required for its TPase reaction.
GTase
FCE contains the conserved 282KTDG285 motif of adenosyl/guanylyl transferases where Lys282 is expected to form a covalent intermediate with GMP (Fig. 1B, center bottom panel; Soulière et al. 2008). We introduced a conservative K282M mutation and tested the enzyme's ability to self guanylylate and guanylylate pp-RNA. Incubation with GTP in the presence of magnesium showed that mutant K282M is incapable of forming the covalent GMP intermediate compared to the WT enzyme (Fig. 3C, left panel). Capping enzyme activity assays showed that the mutant was unable to convert pp-RNA into Gppp− in Reaction 2 but exhibited TPase and N7MTase activities in Reactions 1 and 3, respectively (Fig. 3C, right panel), showing that Lys282 is the catalytic lysine residue required for the guanylyltransferase activity in FCE.
N7MTase
The FCE N7MTase domain contains a conserved 707Q/NFAIHY712 SAM-binding motif, where the tyrosine residue in VCE D1R stacks with the adenosine ring of SAM and the rest of the motif spans the length of SAM-binding pocket (Fig. 1B, right bottom panel). Mutation of Tyr712 to Ala resulted in a significant decrease of N7MTase activity (Fig. 3D, Reaction 3) but did not affect the TPase (Reaction 1) or GTase (Reaction 2) activities, showing the importance of Tyr712 in N7MTase activity in FCE.
Biochemical characterization
FCE exhibits higher specific activity, defined as enzyme concentration to achieve 50% Cap-0 formation (Cap50), than VCE (Fig. 4A). FCE capped 50% of the 5′ ppp-RNA substrate at ∼1 nM compared to VCE that capped 50% of the substrate at ∼2.5 nM at 37°C (Fig. 4A). Similar to VCE, FCE does not exhibit a preference toward RNA initiated with an adenosine or a guanosine (Fig. 4B), the two most abundant 5′ nt of eukaryotic and viral mRNA (Wang et al. 2019). Protein purity and concentrations were verified using SDS-PAGE and densitometry (Supplemental Fig. 3).
(A) FCE exhibits higher specific activity than VCE. RNA capping reactions were performed using decreasing concentrations of two separate preparations of both FCE and VCE under standard conditions at 37°C for 30 min. FCE achieved 50% capping at ∼1 nM enzyme concentration, while VCE achieved the same level of capping at ∼2.5 nM enzyme concentration. (B) FCE and VCE do not discriminate between RNA initiating with guanosine or adenosine. RNA capping reactions were performed on 5′ triphosphate RNA oligonucleotides initiating with a guanosine (+1G) or an adenosine (+1A) under standard conditions at 37°C for 30 min. FCE and VCE exhibited similar level of capping activity on both +1G and +1A RNA. Data points are the average value of triplicated experiments. Error bars represent standard deviations.
pH
FCE exhibits the highest Cap-0 generating activity at pH 7.0–8.0 (Fig. 5A). Interestingly, our data showed that different enzyme activities were affected at either end of the pH range. We found that low pH (pH 6.0) impedes the TPase activity as the level of ppp-RNA remains high. High pH (pH 9.0) impedes the GTase and/or N7MTase activity as the level of pp-RNA remains high (Fig. 5B). VCE exhibits a similar differential response.
Effect of pH on capping activities. (A) FCE and VCE exhibit the highest capping activity at pH 7.0 and pH 8.0. (B) Effect of pH on component enzyme activities. Both enzymes exhibit similar responses outside of optimal pH values. At pH 6.0, the accumulation of ppp− substrate indicates that the TPase activity of both enzymes is negatively affected. At pH 9.0, the accumulation of pp− product indicates that the forward GTase activity of both enzymes is negatively affected. Data points are the average value of triplicated experiments. Error bars represent standard deviations.
NaCl
Similar to VCE, FCE exhibits full Cap-0 capping activity up to 50 mM NaCl at low enzyme concentrations (Fig. 6A). With 100 mM NaCl, both enzymes exhibited lower m7Gppp-generating activity but achieved close to 100% m7Gppp-RNA capping at 5 nM enzyme concentration. Interestingly, the TPase activity of FCE is more sensitive to increasing NaCl concentration than that of VCE, as indicated by the higher fraction of ppp-RNA at 200 mM NaCl and above in the FCE reactions (Fig. 6B).
Effect of NaCl concentration on capping activities. (A) FCE and VCE exhibit the highest capping activity in zero or 50 mM NaCl. At enzyme concentrations of 5 nM or higher, both enzymes exhibit high capping activity in 100 mM NaCl. (B) Differential effect of NaCl on FCE and VCE. In 200 mM NaCl, the TPase activity of FCE decreased significantly whereas the GTase, but not TPase, activity of VCE was affected. The TPase activity of VCE started to be negatively impacted by 300 mM of NaCl. Data points are the average value of triplicated experiments. Error bars represent standard deviations.
MgCl2 concentration
Magnesium ions are required for the TPase and GTase activities. FCE and VCE exhibit the highest Cap-0 generating activity at 0.5–2 mM MgCl2 (Fig. 7). The GTase activity of FCE was significantly impeded at 5 mM MgCl2, as evident by the accumulation of pp-RNA as the reaction product. At the lower end of MgCl2 concentration, the TPase activity of FCE was greatly hampered at 0.1 mM MgCl2, as shown by the accumulation of ppp-RNA substrate, whereas the TPase activity of VCE was less affected under the same conditions. Our results show that FCE and VCE share similar optimal MgCl2 concentrations, but their component enzyme activities respond differently outside of this optimal range.
Differential effect of MgCl2 concentration on capping activities. FCE and VCE exhibit the highest capping activity in 0.5–2 mM MgCl2, but individual activities of the two enzymes are differentially affected outside of the optimal MgCl2 concentration range. At 0.1 mM MgCl2, the TPase activity of FCE is impeded whereas in VCE the GTase activity was affected. For FCE, the GTase activity was negatively impacted at 5 mM MgCl2 with TPase starting to be affected at 10 mM MgCl2. For VCE, the GTase activity started to be affected at 5 mM MgCl2 as indicated by pp-RNA accumulation. However, the TPase activity of VCE was not affected up to 20 mM MgCl2 due to the lack of ppp− accumulation. In each experiment, 5 nM of enzymes and 0.5 mM of RNA substrate were used. Data points are the average value of triplicated experiments. Error bars represent standard deviations.
Reaction temperature
FCE remains highly active for Cap-0 incorporation between 25°C and 55°C, compared to VCE which is most active between 35°C and 45°C (Fig. 8A, blue and orange, respectively). Inspection of intermediate products shows that 5′ pp-RNA began to accumulate at 45°C and became the dominant product at 50°C and above for VCE, indicating a defect in GTase activity (Fig. 8B, bottom panel) and consistent with a previous report that capping activity of VCE can be inactivated by preheating at 55°C (Shatzer et al. 2008). For FCE, accumulation of 5′ pp-RNA product only began at 60°C, where a large fraction of reaction products was Gppp-RNA, suggesting that the N7MTase activity of FCE was compromised (Fig. 8B, top panel). On the other hand, the absence of Gppp-RNA product in the VCE high-temperature reactions suggests that both the N7MTase and the GTase activities were compromised, or the N7MTase activity was compromised and the Gppp− products were rapidly converted back to pp− by the reversible GTase activity (Soulière et al. 2008). It should be noted that for both enzymes, very little-to-no ppp− substrate was detected at all the temperatures tested, showing that the TPase activities are active from 15°C to 60°C.
RNA capping activity of FCE and VCE across a temperature range. RNA capping assays were performed on 0.5 µM of a 20 nt poly(A) synthetic oligonucleotide with a 3′ FAM group using 5 nM of FCE or VCE under standard conditions (see Materials and Methods). (A) FCE generated more m7Gppp-cap at all temperatures tested. Data points are the average of triplicated reactions. Error bars represent standard deviations. (B) A breakdown of the intermediate products across the temperature range after 30 min.
RNA capping activity on structured 5′ ends
A series of RNA oligonucleotides (23–25 nt) were designed to investigate the effect of secondary structure and the availability of a free 5′ end on FCE and VCE capping activity (Fig. 9). A low enzyme concentration (2 nM) was used to assess differences in capping activity of the two enzymes. Under these conditions, FCE generated more Cap-0 product than VCE on an RNA oligo with two or more unpaired nucleotides at the 5′ end within hairpin structures with varying degrees of thermodynamic stability (Oligos 1, 2, and 3). Intriguingly, FCE generated less Cap-0 product on hairpin structures containing a 1-nt 5′ overhang or no overhang compared to VCE (Oligos 4 and 5). Furthermore, FCE generated very little-to-no Cap-0 product on Oligos 5–7 where the 5′ nt was base-paired or recessive within a duplex structure. A good explanation of why each capping enzyme prefers specific structured models of RNA remains elusive without rigorous methods to evaluate the folding dynamics of the RNA molecules and the enzyme–RNA binding kinetics. However, we found that FCE reactions on Oligos 5, 6, and 7 possessed a larger fraction of unreacted ppp-RNA (20%–50%) than on Oligo 4 (>5%) (Supplemental Fig. 4, top panel), suggesting that the TPase activity was less active on the 5′ end with a blunt or recessive structure. Except for Oligo 6, we did not observe a major accumulation of unreacted ppp-RNA in any VCE reactions (Supplemental Fig. 4, bottom panel). This suggests that the TPase activity of VCE is more resilient to 5′ structures than that of FCE. It should be noted that in cases where VCE generated more Cap-0 product than FCE, only 10%–60% of the model hairpin RNA was capped, whereas 80% or more of the RNA was capped when FCE outperformed VCE. In general, our data suggests that FCE exhibits higher capping activity than VCE on stem–looped RNA structures with free 5′ ends (Fuchs et al. 2016), and that the design of the 5′ sequence of synthetic mRNA can have a major impact on the effectiveness of enzymatic RNA capping.
RNA capping activity of FCE and VCE on structured model RNA. Synthetic oligonucleotides are designed to adopt hairpin structures with different folding energies and availability of the 5′ triphosphate. All RNA oligos contain a FAM-dT for capillary electrophoresis readout and their theoretical minimum free energies (MFEs) of unfolding are listed. Capping reactions were carried out on 0.5 µM of the RNA oligos using 2 nM of FCE or VCE at 37°C under standard conditions (see Materials and Methods). Data points are the average of triplicated reactions. Error bars represent standard deviations.
RNA capping activity on long in vitro transcripts with different 5′ sequences
To assess the RNA capping activity of FCE and VCE on in vitro transcripts, 20 nM FCE or VCE was used to cap 3.5 µM (2 µg/µL) of FLuc transcripts (1.8 kb) containing the 5′ sequence of the first 50 nt from human β globin and synthetic transcripts derived from pRNA21, FLuc, and Comirnaty in 100 µL reactions for 1 h at 37°C. The extent of cap incorporation was measured using targeted RNase H cleavage and LC–MS (Chan et al. 2022). Overall, transcripts containing different 5′ sequences were capped to a different extent by both enzymes, and FCE generated a higher percentage of m7Gppp− cap than VCE in the in vitro transcripts tested (Fig. 10). Finally, we showed that FCE and vaccinia virus cap 2′-OMTase can be used in a one-pot format to generate the Cap-1 structure on in vitro transcripts (Supplemental Fig. 5).
RNA capping activity of FCE and VCE on FLuc transcript (1.8 kb) containing a 5′sequence derived from the first 50 nt of human β globin variant 1, synthetic 5′ UTR of pRNA21, synthetic 5′ UTR of pFLuc, and Comirnaty. FCE or VCE (20 nM) was used to cap 3.5 µM of the in vitro transcripts for 1 h at 37°C. Following capping reactions, mRNA capping was measured using targeted RNase H cleavage and LC–MS. Data points are the average of triplicated experiments. Error bars represent standard deviations.
Phylogenetic analysis
Vaccinia virus and Faustovirus belong to a monophyletic group of eukaryotic double-stranded DNA viruses called the NCLDVs. As a group of dsDNA viruses that are diverse in genome size (ranges from 100 kb to 2.5 Mb), virion morphology, and hosts which include amoeba, unicellular algae, insects, and humans, NCLDVs only share three conserved genes—genes that code for a family B DNA polymerase, a primase–helicase, and a predicted transcription factor required for late transcription in poxvirus VLTF3 (Koonin and Yutin 2019). However, NCLDVs also share a core set of genes necessary for cytoplasmic viral transcription and replication and capsid formation (Koonin and Yutin 2019). Cytoplasmic viruses do not have access to the host transcription and RNA capping machineries in the nucleus and face selective pressure to maintain an effective transcription and RNA capping mechanism. Hence, the mRNA capping machinery is broadly conserved among the NCLDV, with the exception of those that have adapted nuclear transcription and replication (Allen et al. 2006; Chinchar et al. 2017).
To examine the conservation of FCE and survey the diversity of RNA capping enzymes within the virosphere, viral genomes were searched for the presence of the cap N7MTase domain (Pfam PF03291). Eighty-three genomes, all belonging to the NCLDV, were identified and clustered based on proteome similarity using the program vContact2 (Fig. 11A; Bin Jang et al. 2019). The amino acid sequences of representative members of the identified capping enzyme gene clusters were aligned, and their sequence similarity is represented in the form of a phylogenetic tree (Fig. 11B). These capping enzyme genes were then mapped to the proteome similarity networks. We found that phylogenetic analysis of the capping enzymes is largely consistent with proteome similarity networks consistent with the notion that broad conservation and correlation with proteome similarity highlights the importance of capping enzyme to the lifestyles of these viruses (Fig. 11A; Koonin and Yutin 2019). The capping enzyme from Faustoviruses is most similar to those from African swine fever viruses and Pacmanviruses (Fig. 11B, blue dots) and these three families of viruses cluster in the same proteome similarity network (Fig. 11A, blue dots). Most of the clades, except Panodoraviruses (orange dots) and Phycodnaviruses (red dots), possess putative RNA capping enzymes as multifunctional proteins containing all three enzyme activities. Poxvirus capping enzymes are the only group of enzymes that possess an obvious activation subunit (similar to the D12L of vaccinia virus) and form a distinct branch on the phylogenetic tree (Fig. 11B, green dots). Similar to a previous study, the poxviruses form a distinct group in our proteome similarity networks analysis (Fig. 11A, green dots; Guglielmini et al. 2019). Interestingly, our analysis suggests that the phylogenetic relationship of RNA capping enzymes does not always correlate with proteome similarity. For example, Pandoravirus and Phycodnavirus capping enzymes have separated TPase, MTase, and GTase genes and are clustered under a common ancestor distinct from the multifunctional capping enzymes (Fig. 11B, orange and red dots), even though the two groups of NCLDVs are clustered to different proteome similarity networks. On the other hand, although the mimivirus capping enzyme belongs to the clade (Fig. 11B, yellow dots) that is phylogenetically distinct from the clade containing the Phycodnavirus capping enzyme (Fig. 11B, red dots), the two groups are clustered to connected networks in proteome similarity network analysis (Fig. 11A).
(A) The CheckV database of complete viral genomes was searched for viruses which contain a cap M7MTase domain using hmmsearch and the PFAM profile PF03291. The 83 identified viruses were clustered by proteome similarity using the vContact2 program. The program was run using DIAMOND (Buchfink et al. 2021) for the all-vs-all protein alignment, and the figure was prepared by excluding the reference virus database. (B) A subset of the 83 capping enzymes was selected to represent the diversity of the data set, aligned and used to construct a phylogenetic tree. The tree is built from an MAFFT alignment of the N7MTase domains; for viruses which contain separate genes from TPase, GTase and N7MTase domains, only the N7MTase protein was used. The viruses were mapped to the clustering in (A) (color dots), showing that the clustering based on the whole proteome is consistent with the phylogeny of the cap N7MTase genes. On the right, putative TPase (TP), GTase (GT), and N7MTase (MT) domains, as well as homologs of the regulatory subunit (D12L) of VCE are shown as white hexagons. Domains attached by solid lines are found within one gene.
DISCUSSION
In this communication, we present biochemical and enzymatic characterization of a single-subunit RNA capping enzyme from Faustovirus. Faustovirus has been shown to infect Vermamoeba vermiformis, an amoeba that has been found in human-associated environments such as hospital water networks, drinking water, and human stool samples (Bradbury 2014; Hajialilo et al. 2015; Pagnier et al. 2015; Delafont et al. 2018; Rolland et al. 2019). As a giant virus (NCLDV that has a genome size of 400 kb or larger) with a 460 kb dsDNA genome, our phylogenetic analysis, based on two universally conserved genes (Family B DNA polymerase and major capsid protein), as well as a previously published analysis (Koonin and Yutin 2019), suggests that Faustovirus is more closely related to Asfavirus (genome size of 170–190 kb) than giant viruses of similar genome sizes (such as Cedratvirus of 590 kb or Marseillevirus of 360 kb) (Koonin and Yutin 2019). Amino acid sequence analysis also showed that FCE is most closely related to the Asfavirus RNA capping enzyme (Fig. 11B, blue dots). Our phylogenetic analysis suggests that NCLDVs derive their RNA capping enzymes from a common ancestral enzyme which is a single-subunit protein with multiple activities. It is also likely that the functional dependence of an activation subunit of poxvirus RNA capping enzymes (e.g., vaccinia virus D12L) occurred after the speciation event that produced the ancestors of the Poxviridae and Asfarviridae such that poxvirus RNA capping enzymes remains the only NCLDV capping enzyme that requires an activation subunit.
Similar to NCLDV RNA capping enzymes, RNA capping enzymes in eukaryotes adopt a diverse configuration. Fungal capping enzymes are composed of separate TPase, GTase, and N7MTase proteins (e.g., Cet1, Ceg1, and Abd1 of S. cerevisiae), whereas higher eukaryotes have a bifunction TPase and GTase protein (e.g., RNGTT of humans and Drosophila), which contains a cysteinyl-phosphatase domain for TPase activity, and a separate N7MTase protein (e.g., RNMT of humans) (Yue et al. 1997; Changela et al. 2001). It is not clear whether the separation of these enzyme activities on individual proteins contributes to any biological functions or regulation. But it is fascinating that the modular nature of RNA capping enzymes allows the evolution of so many variations throughout eukaryotes and their viruses (see Ramanathan et al. 2016, for an overview of the RNA capping machineries of RNA viruses), while each being functionally adapted to the organisms as an essential component of the transcription–translation machinery.
With higher specific activity and broader temperature range than VCE, FCE may be a preferred choice for mRNA manufacturing applications and may enable major improvements in mRNA manufacturing workflows. For example, low-temperature capping may help prevent the autohydrolysis of long synthetic RNA molecules. Capping at higher temperatures may increase capping efficiency on structured RNA. On the other hand, VCE accepts some modified GTP as cap donor (Ohno et al. 2023), such that VCE is a better choice for applications such as selective enrichment of transcription start sites using Cappable-Seq and ReCappable-Seq (Ettwiller et al. 2016; Ugolini et al. 2022; Yan et al. 2022), labeling of in vitro transcripts using a fluorescent derivative of GTP (Gunawardana et al. 2015), or applications that require derivatization of the 5′ end of in vitro transcripts (Ohno et al. 2023).
As with most enzymes, in vitro activities of FCE can vary based on reaction conditions such as enzyme, substrate and reagent concentrations, supplementation of additives, temperature, and accessibility of the substrate. For RNA capping enzymes, the true substrate is the 5′ triphosphate group of newly synthesized RNA. Our data and previously published results on model RNA oligonucleotides showed that the secondary structure of the 5′ region plays a significant role in the effectiveness of enzymatic capping (Fig. 9; Fuchs et al. 2016). We further demonstrated that variation of the 5′ sequence of long in vitro transcripts can have a significant impact on enzymatic capping (Fig. 10). Therefore, exploring the sequence space for highly “cappable” 5′ ends, in addition to their in vivo expression levels (Sample et al. 2019; Jia et al. 2020; Leppek et al. 2022), would be beneficial for manufacturing high-quality mRNA that offers high translation efficiency and low immunogenicity for therapeutic applications.
MATERIALS AND METHODS
Enzymes and RNA substrates
Faustovirus RNA capping enzyme (FCE) and VCE were obtained from New England Biolabs, Inc. FCE mutant constructs were generated by GenScript by site-directed mutagenesis. His-tagged FCE mutant proteins were expressed from Escherichia coli BL21(DE3)star and purified by Ni-column chromatography and size-exclusion chromatography by GenScript to over 90% homogeneity (determined by Coomassie blue-stained SDS-PAGE). RNA oligonucleotides containing 5′ triphosphate, 5′ diphosphate, 5′ Gppp− cap, and 5′ m7Gppp− cap were synthesized by Bio-Synthesis, Inc.
Enzyme assays
Unless specified otherwise, RNA capping enzyme assays were performed in 10 µL of 1× RNA capping buffer (New England Biolabs, Inc.) containing 0.5 µM substrate RNA, 0.1 mM SAM, 0.5 mM GTP, 1 mM DTT, and a specified concentration of capping enzymes. The reactions were carried out in a thermocycler for 30 min at a specific temperature and quenched by adding an equal volume of a 2× stop solution containing 20 mM EDTA and 2% SDS. The quenched reactions were diluted for capillary electrophoresis on an Applied Biosystems 3730xl Genetic Analyzer (96 capillary array) using POP-7 polymer and GeneScan 120 LIZ dye Size Standard (Greenough et al. 2016; Wulf et al. 2019; Chan et al. 2022). Peak detection and quantification were performed using the Peak Scanner software v.1.0 (Thermo Fisher Scientific) and an in-house data analysis suite. Validation of the capillary electrophoresis platform for RNA cap analysis has been published (Chan et al. 2022).
Verification of the identity of 5′ cap
RNA oligonucleotides starting with a +1G (pppG-25mer) and +1A (pppA-25mer) were capped using FCE under standard conditions: 0.5 µM substrate RNA, 0.1 mM SAM, 0.5 mM GTP, 1 mM DTT, 100 nM FCE in 10 µL of 1× RNA capping buffer for 30 min at 37°C. The reactions were treated by Nucleoside Digestion Mix (New England Biolabs, Inc.) to release the nucleosides and the dinucleotide cap structure (m7GpppG or m7GpppA). Specifically, 4 µL of 10× Nucleoside Digestion Mix reaction buffer, 24 µL of nuclease-free water, and 2 µL of the Nuclease Digestion Mix were added to 10 µL of the capping reactions. After the reactions were incubated at 37°C for 1 h, the samples were subjected to LC–MS nucleoside analysis using an Agilent 1200 Series LC–MS System equipped with a G1315D diode array detector and a 6120 Single Quadrupole Mass Detector in both positive (+ESI) and negative (−ESI) electrospray ionization modes. LC was performed on a Waters Atlantis T3 column (4.6 × 150 mm, 3 µm) with a gradient mobile phase consisting of aqueous ammonium acetate (10 mM, pH 4.5) and methanol. The identity of each peak was confirmed by MS. The relative abundance of nucleosides was determined by comparing the respective peak area at 260 nm or their respective UV absorption maxima to that of cytidine in each sample.
Self-guanylylation assay
In a 20 µL reaction, 100 nM of WT FCE or mutant K282M was incubated with 1 µM of 32P-α-GTP (PerkinElmer) in a 20 µL reaction containing 1× RNA capping buffer at room temperature for 1 min. The reactions were quenched by adding 10 µL of 3× SDS loading dye (New England Biolabs, Inc.). The quenched reactions were analyzed by SDS-PAGE on a 10%–20% Tris-Glycine SDS polyacrylamide mini-gel (Thermo Fisher). The presence of the GMP-GTase covalent intermediate was visualized by autoradiography using a Phosphor Storage Screen (Cytiva) and the Amersham Typhoon RGB scanner (Cytiva).
Capping activity assay on 5′ secondary structures
A series of RNA oligonucleotides was designed to possess stem–loop structures of different stability and different accessibility to the 5′ triphosphate group (Fig. 9). Oligo 1 is expected to form a short 3-bp stem with an unstructured 5′ end. Oligo 2 is expected to form a longer stem structure with mismatches such that the theoretical MFE of unfolding is a moderate −3.5 kcal/mol with the first two nucleotides free from base-pairing. Oligo 3 is also expected to have two free nucleotides at the 5′ end but adopt a more stable stem–loop structure with an MFE of −7.2 kcal/mol. Oligos 4 and 5 have the same secondary structure but are expected to have one free nucleotide at the 5′ end and a fully base-paired 5′ end, respectively. Oligos 6 and 7 are expected to have a slightly more stable stem–loop structure with MFE of −8.9 kcal/mol and one- or two-base recessive 5′ ends, respectively. The RNA oligonucleotides contained a FAM-conjugated dT in the loop region to facilitate readout using capillary electrophoresis. To facilitate differentiating capping activity of FCE and VCE on the structured RNA oligonucleotides, a lower concentration of enzymes (2 nM) was used in the assay with 0.5 µM of the structured RNA oligos at 37°C for 30 min under standard conditions. The reactions were quenched and analyzed by capillary electrophoresis as described above.
RNA capping activity of FCE and VCE on long in vitro transcripts with different 5′ sequences
To assess the RNA capping activity of FCE and VCE on long in vitro transcripts, a plasmid containing a firefly luciferase gene (FLuc) under the control of the standard T7 phi10 promoter was constructed. To investigate the effect of the 5′ sequence on RNA capping enzymes, the first 50 nt of the FLuc transcribing region was replaced with the first 50 nt of human globin subunit β (HBB) mRNA (GenBank accession NM_000518.5), synthetic transcript from pRNA21 or Comirnaty. T7 promoter phi2.5 was used in the constructs carrying HBB and Comirnaty to facilitate efficient transcription of the transcripts initiating with an adenosine. IVT was done using the HiScribe T7 High Yield RNA Synthesis Kit (NEB) following the manufacturer's instructions. The in vitro transcripts (1.8 kb) were purified using the Monarch RNA Clean-up Kit (50 µg) (NEB). RNA capping reactions were carried out using 20 nM of FCE or VCE in 20 µL reactions containing 1× FCE reaction buffer, 0.1 mM SAM, 0.5 mM GTP, 1 mM DTT, and 3.5 µM (2 µg/µL) in vitro transcripts at 37°C for 1 h. The extent of cap incorporation was assessed by an RNase H-based cap analysis method published previously (Chan et al. 2022).
One-pot Cap-1 incorporation using FCE and vaccinia RNA cap 2′-O-methyltransferase
To determine whether FCE is compatible with vaccinia RNA cap 2′-OMTase in a one-pot format, 50 or 100 nM of FCE was incubated with 0.5 µM of the synthetic FLuc in vitro transcript (1.8 kb) (above) in the presence or absence of 5 U/µL vaccinia cap 2′-OMTase (NEB) in 20 µL reactions containing 1× RNA capping buffer (New England Biolabs, Inc.), 0.2 mM SAM, 0.5 mM GTP, and 1 mM DTT at 37°C for 1 h. The extent of Cap-1 incorporation was analyzed by an RNase H-based cap analysis method previously described (Chan et al. 2022).
Phylogenetic analysis
The CheckV database of viral genomes (Nayfach et al. 2021) was searched for viruses encoding capping enzymes using the HMMER program (Wheeler and Eddy 2013) and the PFam profile for the Poxvirus N7MTase domain (Pfam: PF03291). The resulting viral genomes were all members of the NCLDV, and these were clustered on the basis of proteome similarity using vContact2 (Bin Jang et al. 2019). Associated capping enzymes were aligned with MAFFT (Madeira et al. 2022) using the G-INS-I algorithm and used as an input for FastTree (Price et al. 2010) with default parameters to produce a phylogenetic tree which was rooted according to Guglielmini et al. (2019). The initial alignment of capping enzymes was used together with the VCE crystal structure (Kyrieleis et al. 2014) and predicted structures (Mirdita et al. 2022) to define boundaries between TPase, GTase, and MTase domains. Profiles for these domains were constructed and used to identify distant homologs of the poxvirus capping enzyme domains in viral genomes for which they occur in separate genes. To identify vaccinia D12L homologs, NCLDV-translated proteomes were searched with the existing profile hidden Markov model (HMM) (PFAM: PF03341).
SUPPLEMENTAL MATERIAL
Supplemental material is available for this article.
COMPETING INTEREST STATEMENT
S.H.C., C.N.M., D.N., L.M., N.D., J.B., D.W.K., J.M.W., and G.B.R. are current employees of New England Biolabs that markets molecular biology reagents. The affiliation does not affect the authors’ impartiality, adherence to journal standards and policies, or availability of data.
ACKNOWLEDGMENTS
We would like to thank the NEB sequencing core for the capillary electrophoresis analysis as well as all the facilities personnel for daily support of the operation. This work is supported by New England Biolabs, Inc.
Footnotes
-
Article is online at http://www.rnajournal.org/cgi/doi/10.1261/rna.079738.123.
-
Freely available online through the RNA Open Access option.
- Received June 5, 2023.
- Accepted August 7, 2023.
This article, published in RNA, is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), as described at http://creativecommons.org/licenses/by-nc/4.0/.





















