Identification of Cis-Acting Elements in the 3′-Untranslated Region of the Dengue Virus Type 2 RNA That Modulate Translation and Replication*

Using the massively parallel genetic algorithm for RNA folding, we show that the core region of the 3′-untranslated region of the dengue virus (DENV) RNA can form two dumbbell structures (5′- and 3′-DBs) of unequal frequencies of occurrence. These structures have the propensity to form two potential pseudoknots between identical five-nucleotide terminal loops 1 and 2 (TL1 and TL2) and their complementary pseudoknot motifs, PK2 and PK1. Mutagenesis using a DENV2 replicon RNA encoding the Renilla luciferase reporter indicated that all four motifs and the conserved sequence 2 (CS2) element within the 3′-DB are important for replication. However, for translation, mutation of TL1 alone does not have any effect; TL2 mutation has only a modest effect in translation, but translation is reduced by ∼60% in the TL1/TL2 double mutant, indicating that TL1 exhibits a cooperative synergy with TL2 in translation. Despite the variable contributions of individual TL and PK motifs in translation, WT levels are achieved when the complementarity between TL1/PK2 and TL2/PK1 is maintained even under conditions of inhibition of the translation initiation factor 4E function mediated by LY294002 via a noncanonical pathway. Taken together, our results indicate that the cis-acting RNA elements in the core region of DENV2 RNA that include two DB structures are required not only for RNA replication but also for optimal translation.

severely affected viral replication (16,28). Downstream of the 5Ј-UTR in the capsid-coding region is the capsid-coding region hairpin RNA (cHP) (Fig. 1A), which was shown to play a role for the efficient selection of the start codon for initiation of translation of the single polyprotein (29,30). Moreover, it is also required for DENV and WNV replication (29,30). Recently, another base-paired interaction between the motifs known as the "5Ј-3Ј downstream of AUG region," or DAR, was also identified to be required for viral RNA replication (31,32).
The 3Ј-UTR of flaviviruses exhibits a great deal of sequence divergence and size heterogeneity but also contains conserved sequences that are required for replication (9 -14, 16, 17, 33). It consists of the variable region, core region (CR), and terminal stem-loop (3Ј-SL) region (Fig. 1A).
The 3Ј-SL region is formed within ϳ100 nt from the 3Ј-end (34,35), which has been shown to play a significant role in replication (33, 36 -39). The CR includes the 3Ј-cyclization sequence (3Ј-CS1) containing a short motif that is complementary to the 5Ј cyclization sequence (5Ј-CS) located in the capsid coding region (Fig. 1A). A long range RNA-RNA interaction between these complementary motifs resulting in circularization of the genome was postulated (40) and demonstrated by atomic force microscopy (13). The physical and functional RNA-RNA interaction between 5Ј-and 3Ј-CS1 was shown to be essential for RNA synthesis in vitro (10,41) and in replication of subgenomic replicons or infectious clones in cultured mammalian cells (11,12,14,42,43). In addition to the 5Ј-and 3Ј-CS1, two regions of complementarity in the 5Ј-and 3Ј-terminal sequences, termed the upstream of AUG region (UAR), were also identified as a requirement for cyclization and RNA replication (13,17,44). The 5Ј-UAR is located upstream of the 5Ј-CS, and the 3Ј-UAR is downstream of 3Ј-CS1 (Fig. 1A).
The CR exhibits relatively high sequence conservation among members of MBFV and is believed to fold into well defined secondary structures independent of other regions of the 3Ј-UTR (45,46). Within this region, upstream from the 3Ј-CS1, the CS2 is present in all three subgroups in the MBFV, whereas the repeat conserved sequence 2 (RCS2) is present only in JEV and DENV subgroups (40) (Fig. 1B).
The regions containing CS2 and RCS2 are part of two almost identical dumbbell-like (DB) secondary structures that have been postulated to form in DENV (Fig. 1B). It is thought that this duplication of dumbbell-like structures arose because of the repetition of CS2 and RCS2 sequences in the 3Ј-and 5Ј-DB, respectively. The leftmost loops of the 5Ј-and 3Ј-DBs contain a set of identical sequences (underlined in 5Ј-GAAGCUGUA-3Ј) (45). Sequence analysis of different serotypes of DENV showed that the internal five nucleotides in these loops (referred to as TL1 and TL2) and the two complementary pentanucleotide sequences downstream of each DB, PK2 (5Ј-GCAGC-3Ј) and PK1 (5Ј-ACAGC-3Ј), are highly conserved, suggesting possible base pair interactions between TL1/PK2 and TL2/PK1 to form two pseudoknots (5Јand 3Ј-) (45).
The role of the CR region in flavivirus replication or translation is not understood at present. In this study, we sought to examine the functional roles of the TLs and PKs by using a Renilla luciferase (Rluc) reporter replicon in BHK21 cells. We show here that all four motifs played a crucial role in replication in a sequence-dependent manner. In contrast, there seems to be a differential role for TL1/PK2 and TL2/PK1 in translation. The lack of functional similarity between TL1 and TL2 despite their identical sequences within two DBs could be explained by results of our analysis of RNA secondary structures, their stabilities, and frequency of occurrence computed by the massively parallel genetic algorithm (MPGAfold) (47)(48)(49)(50)(51). The results of our study, taken together, provide a new insight into the role of the cis-acting elements in the CR region in translation and replication.

Cells and Primers
BHK-21 cells (ATCC, Manassas, VA) were maintained at 37°C with 5% CO 2 in complete Earle's minimum essential medium (EMEM) (CellGro, Manassas, VA), containing high glucose (1 g/liter) and L-glutamine, and supplemented with 10% fetal bovine serum and antibiotics (Invitrogen). All primers were obtained from Integrated DNA Technologies (Coralville, IA) and resuspended in H 2 O. Primer sequences are listed in supplemental Table 1.

Replicon Construction
DENV2 Rluc Replicon-The Rluc gene was amplified from phRL-SV40 (Promega, Madison, WI) placing SalI sites on both the 5Ј-and 3Ј-ends using RlucF and RlucR primers (supplemental Table 1). The resulting fragment was cloned into pGEMTEasy (Promega). The resulting insert vector was used as a template using RlucF and RlucIresR primers (supplemental Table 1) to create a fragment containing Rluc with a SalI site at the 5Ј-end and the 5Ј-terminus of EMCV IRES at the 3Ј-end. Next, using the GFP-expressing replicon, pRS424GFPIresDEN2 (52), as a template and primers IresF and NS1BstEIIR, a fragment containing the entire sequence of EMCV IRES followed by the 5Ј-end of NS1 up to the BstEII site was created. Overlap PCR was performed using the products from the first two PCRs and primers RlucF and NS1BstEIIR. The entire fragment was cloned into pGEMTEasy. The insert was subcloned into the replicon plasmid pRS424GFPIresDEN2 using SalI and BstEII (New England Biolabs, Ipswich, MA) to create the DENV2 Rluc replicon.
GND Replicon-A replication-defective replicon clone was created by engineering a point mutation in the conserved GDD motif of the RdRP domain in NS5 (nucleotides 9553-9561). First, an AatII/ApaI restriction fragment from the C-terminal region of NS5 to the ApaI site at the 3Ј-UTR was isolated from pRS424GFPIresDEN2 and subcloned into pBR322 vector. The GDD 3 GND mutation was introduced by site-directed mutagenesis using primers GNDF and GNDR (supplemental Table 1). The fragment was released with PmlI and ApaI and subcloned into XhoI-SacI sites of pBR322, which contained the replicon sequence from the XhoI site in NS3 to the SacI site at the end of the 3Ј-UTR. The sequence between XhoI and SacI was cloned into the WT DENV2 Rluc replicon to yield the GND mutant.

Replicon Mutants
The 5Ј-or 3Ј-UTR fragments containing the desired mutations were first amplified by overlap extension site-directed mutagenesis PCR using primers containing desired mutations (supplemental Table 1) and were purified from agarose gel using the Zymoclean gel DNA recovery kit (Zymo Research, Irvine, CA). Three purified PCRs were combined in a single tube and vacuum-dried. Twenty micrograms of WT replicon were digested overnight with ϳ40 units of SalI for 5Ј-UTR mutations or BbvCI for 3Ј-UTR mutations. These were purified by phenol/chloroform/isoamyl alcohol (25:24:1), ethanol-precipitated, and resuspended in 5 l of water. The linearized WT or the mutant replicon plasmid was mixed with the vacuumdried mutant PCR fragments, and the mixture was introduced into Saccharomyces cerevisiae YPH857 (50 l) made chemically competent using the S.c. EasyComp Transformation Kit (Invitrogen) (53). The 5Ј-or 3Ј-UTRs of mutant replicons were amplified from the colonies, and their sequences were verified (MC Lab, South San Francisco, CA). Plasmid DNAs harboring the mutations were recovered using the Zymoprep II yeast plasmid minipreparation kit (Zymo Research), and 2 l were used to transform 20 l of MAX Efficiency Stbl2-competent cells (Invitrogen). Replicon plasmids were propagated in Escherichia coli from single colonies, and their sequences within 5Ј-and 3Ј-UTRs were again verified. Details of sequences for mutant replicons used in this study are shown in supplemental Table 1 and Figs. 3A, 4A, and 5A.

Templates for in Vitro RdRP Assay
DNA templates for in vitro transcription harboring the desired mutations in a 719-nt template containing the first 230 nt from the 5Ј-end and the last 489 nt of the genome (54) were prepared by an overlap PCR amplification strategy as follows. In the first PCR, the region spanning the 5Ј-UTR and the NS5-3Ј-UTR junction was amplified using the pSY2 plasmid as the template and the primers pSY2 EcoRI and NS5-3Ј-UTRr. In the second PCR, the full-length 3Ј-UTR was amplified from the replicon mutants using the start of 3Ј-UTR and 3Ј-end DEN2 primers (see supplemental Table 1 for primer sequences). The two fragments were joined together by an overlap PCR and purified from the agarose gel.

In Vitro Transcription of RNA Templates
Approximately 50 g of WT and mutant replicons were linearized with EcoICRI and purified twice with phenol/chloroform/isoamyl alcohol and once with chloroform to remove excess phenol. DNA was ethanol-precipitated and resuspended in 20 l of RNase-free water. Three microliters were used to make RNA with the Ampliscribe SP6 high yield transcription kit (Epicenter Biotechnologies, Madison, WI). The reaction mixtures (20 l) contained 3.13 mM GTP, 8 mM m 7 G(5Ј)ppp(5Ј)G cap analog (New England Biolabs) and 80 units of RNase inhibitor (New England Biolabs) and incubated for 3-4 h at 37°C. Following transcription, DNAse I (Epicenter) was directly added to the mixture and incubated for another 1 h at 37°C. 0.5 l of each reaction was visualized on a 0.5% agarose gel to verify the integrity of the RNA prior to storage or use. The plasmid encoding the rabbit ␣-globin 5Јand 3Ј-UTRs flanking the firefly luciferase gene (pGLGpA) (15) (kindly provided by Dr. Theo Dreher of Oregon State University) was linearized by Acc65I digestion and used for the in vitro transcription. The GLGpA RNA was used as an internal control for transfections of the WT and mutant replicons. GLGpA and (ϩ) templates for the in vitro RdRP assays were made using the Ampliscribe T7 high yield transcription kit.

RNA Transfection
For each reaction, ϳ10 6 BHK-21 cells were resuspended in 100 l of Ingenio solution (Mirus Bio, Madison, WI), and ϳ3 g of replicon RNA and 0.1 g of GLGpA RNA were added. Transfections were performed using a Nucleofector II electroporator (Amaxa Biosystems, Cologne, Germany) using program A031. Following pulsing, cells were carefully transferred in a tube containing 1 ml of prewarmed complete EMEM. Cells were allowed to recover at 37°C for 5 min and then transferred to another tube with 3 ml of complete EMEM. Cells were plated in 8 wells of a 48-well plate, each containing 400 l. The rest of the cells were plated in a 6-well plate containing 3 ml of medium. Each transfection experiment was done in quadruplicate and was repeated at least three times using at least two different batches of in vitro transcribed RNA.

Cap-independent Translation Assay
In vitro transcribed mutant RNAs were transfected into BHK21 cells as described above. After transfection, cells were allowed to recover in 1 ml of complete EMEM at 37°C for 5 min. 100 l was seeded in each of 10 wells of a 48-well plate containing 100 l of medium. Half of the wells were used as no-treatment controls, whereas the other half contained LY294002 with a final concentration of 40 M (Cayman Chemical, Ann Arbor, MI). Cells were lysed at 2 h post-transfection (hpt) and assayed as described below.

Luciferase Assays
At the desired time points, cells were lysed with 60 l of 1ϫ Renilla luciferase lysis buffer to determine the Rluc signal using a kit (Promega). Three additional wells were also lysed at 2 hpt with 1ϫ Passive Lysis Buffer (Promega) to measure the GLGpA Fluc activity. To ensure complete lysis, plates were put on an orbital shaker for 30 min. Rluc activity was measured using a Centro LB 960 luminometer (Berthold Technologies) by injecting 100 l of 1ϫ Renilla luciferase substrate and reading for 10 s after a 2-s delay. Fluc readings were similarly measured using 1ϫ luciferase assay reagent (Promega).

Immunofluorescence
Following transfection, cells were plated into wells of a 2-well chamber slides (BD Falcon). At 96 hpt, the culture medium was removed, and the cells were washed with PBS. The cells were fixed in 100% cold methanol and placed at Ϫ20°C for 30 min. After fixing, the cells were blocked in 1ϫ PBS with 1% nonfat dry milk for 1 h with rocking. Cells were washed three times with PBS and then incubated with rabbit anti-NS5 IgG (1:250) in PBS for 2 h with rocking. After washing, cells were incubated with goat anti-rabbit IgG conjugated to FITC (ICN, Solon, OH) (1:250) for 2 h with rocking. The secondary antibody was removed, and cells were washed six times in PBS, mounted using ProLong reagent (Molecular Probes), and visualized using an Olympus Fluoview FV300 laser confocal microscope.

Quantitative PCR
Early Time Points-To determine the RNA stabilities of replicon RNAs, we transfected 6 g of each RNA by electroporation into 2 ϫ 10 6 BHK-21 cells as described above, and the transfected cells were plated into 24-well plates (2 ϫ 10 5 cells/ well). At the desired time points (1, 2, 4, and 6 hpt), cells were harvested. Total RNA was extracted using 250 l of TRIzol reagent (Invitrogen) following the manufacturer's protocol, resuspended in 40 l of diethylpyrocarbonate-treated water, and stored immediately at Ϫ80°C until used for qPCR analysis. Each RNA sample was analyzed in duplicate wells of a 96-well plate, and the assay was performed on the Applied Biosystems 7900HT fast real-time PCR system. Primers and probe were targeted to amplify nucleotides 9725-9820 in the NS5 gene. The qPCR reaction mixture (50 l) contained 100 ng of extracted RNA, a 0.2 M concentration of primers qPCR NS5F and qPCR NS5R (supplemental Table 1), 0.1 M qPCR NS5 probe (supplemental Table 1) (from Dr. Robin Levis, Food and Drug Administration), 1ϫ rTth EZ buffer, 300 M each dNTP, 3 mM Mn(OAc) 2 , and 5 units of rTth DNA polymerase (Applied Biosystems). The reactions were subjected to 60°C for 30 min for reverse transcription, followed by 94°C for 1 min, and PCRs were amplified through 40 cycles of 94°C for 15 s and 60°C for 1 min. RNA was quantified by reference to RNA extracted from a virus stock with a known titer.
96 hpt Time Point-Transfected cells grown in 6-well plates for 96 hpt were washed with PBS and lysed with 1 ml of TRIzol reagent, and total RNA was purified following the manufacturer's protocol. cDNA was synthesized from 1.5 g of total RNA using the iScript cDNA synthesis kit (Bio-Rad). qPCR analysis was done on a iQ5 multicolor real-time PCR detection system (Bio-Rad) on 2 l of the cDNA reaction together with a 0.25 M concentration of qPCR NS1F and qPCR NS1R primers (supplemental Table 1; primer sequences from R. Takhampunya) in a 20-l reaction using iQ SYBR Green Supermix (Bio-Rad). Each data point was done in triplicate. Thermocycling conditions were as follows: 95°C for 3 min, 35 cycles of 95°C for 20 s, and 30 s each at 55°C and 72°C. Measurements of SYBR Green signal were done at the annealing step. RNA copy numbers were determined using purified PCR products of the NS1 region of known concentrations. Values were further normalized from GAPDH qPCR results generated using the same protocol and cDNA batch but with primers qPCR GAPDHF and qPCR GAPDHR (supplemental Table 1).
After induction, cells were harvested and lysed with 30 ml of lysis buffer (50 mM HEPES, 0.3 M NaCl, 10% glycerol, 1% Non-idet P-40, pH 7) supplemented with 500 l of Protease Inhibitor Mixture VII (Calbiochem) and lysozyme (Pierce) on ice for 30 min. The lysates were sonicated on ice for 20 min total processing time (15 s on, 45 s off) and were centrifuged at 18,000 ϫ g for 45 min. Supernatant was allowed to bind to TALON resin (Clontech) in a 50-ml tube for 1 h in a cold room on a shaker.
The resin was washed five times with 10 ml of Buffer A (50 mM NaH 2 PO 4 , 0.3 M NaCl, 10% glycerol, pH 7). At the last wash, the slurry was transferred and packed in a 0.8 ϫ 4-cm Poly-Prep chromatography column (Bio-Rad). The protein was eluted with 10 ml of buffer A containing 150 mM imidazole. Fractions of 500 l were collected, and the presence of protein was monitored by a Bradford assay on a 96-well plate. The eluted fractions containing the most protein were combined in a single tube, and 20 units of SUMO protease I (Life Sensors, Malvern, PA) was added to cleave the N-terminal histidine tag to produce a full-length NS5 with an authentic N terminus. Digestion was allowed to proceed overnight on ice in the cold room. The purified protein was dialyzed using Spectra/Por 6 dialysis tubing (50,000 molecular weight cut-off; Spectrum Laboratories, Rancho Dominguez, CA) for 6 h in 1 liter of enzyme buffer (50 mM Tris-Cl, 50 mM NaCl, 5 mM MgCl 2 , 40% glycerol, 1 mM fresh dithiothreitol, pH 7.5).

In Vitro RdRP Assay
Each 25-l RdRP reaction contained 150 ng of RNA template, 4 Ci of [␣-32 P]GTP, 0.5 mM ATP, CTP, UTP, 12.5 M GTP, and 20 nM NS5 in 50 mM HEPES, 10 mM KCl, 5 mM MgCl 2 , 2 mM MnCl 2 , pH 8. The reaction was incubated at 37°C for 30 min. The RdRP products were purified as described previously (10) and resuspended in 20 l of diethylpyrocarbonate-treated water. A 5-l aliquot was mixed with an equal amount of Gel Loading Buffer II (Ambion) and was heat-denatured at 70°C for 10 min and then flash-cooled on ice for 5 min. Samples were run on a 8 M urea, 4% polyacrylamide gel in 1ϫ TBE buffer at 125 V and visualized using a Storm 840 PhosphorImager (Amersham Biosciences). Band intensities were quantified using ImageJ version 1.45 software (55).

MPGAfold Analysis
For secondary structure prediction analysis, we used the MPGAfold (47-51, 56, 57) to fold a 719-nt subgenomic RNA (minigenome) from the DENV2 (New Guinea C strain) containing the 5Ј-terminal 226 nt, 42 nt from the C-terminal coding region of NS5, including the UAG termination codon, and the 451-nt 3Ј-UTR derived from pSY2 (10). Previous studies have shown that this RNA molecule contains the essential cis-acting elements required for efficient translation (15,29,30,58) as well as negative strand RNA synthesis in vitro (10,41,54,59). MPGAfold utilizes the genetic operators of mutation, recombination, and selection in parallel on a high performance computer. The program uses previously reported energy rules (60), including efn2 coaxial energy calculations at run time, to drive a population (e.g. 16,384) of RNA structures toward a set of structures that have high fitness (low free energy). MPGAfold is capable of predicting RNA secondary structure folding dynamics, focusing on the formation of significant final structures and folding intermediates (61-65). Thus, MPGAfold is especially suited for identifying and visualizing, in addition to highly stable RNA conformations with low free energy, metastable states that potentially occur in a folding pathway. Because the algorithm is nondeterministic, 100 independent runs at a 16,384 population level were performed for each WT and mutant RNA to determine the structure or structural components that form a consensus. StructureLab (56,66) and, more specifically, the StemTrace (51) component of StructureLab were used to analyze the results and to obtain the frequencies of structural motifs (supplemental Table 2).

Statistical Methods
Rluc readings were normalized against the initial amount of RNA transfected and the corresponding Fluc signal. These were expressed as a percentage of WT. Because our objective is to determine the effect of each mutation on translation and replication, the values for each group were compared with WT using unpaired Student's t test (␣ ϭ 0.01), using GraphPad Prism version 5.03 (GraphPad Software, San Diego, CA).

RESULTS
The 3Ј-DB Forms a Stable Secondary Structure, Whereas the 5Ј-DB Occurs Mostly in RNA Folding Intermediates as Revealed by MPGAfold-A number of studies have indicated that secondary structures in both 5Ј-and 3Ј-terminal regions are required for efficient viral RNA translation (15,29,30,58) and viral RNA synthesis (10 -12, 16, 41). In this study, we focused on the role of 5Ј-and 3Ј-DBs and potential pseudoknot base pair interactions within the CR of DENV2 RNA in translation and replication, which is currently not understood. This knowledge is a prerequisite for a detailed analysis of the trans-acting factors that interact with these elements in mediating their effects in these processes.
Earlier RNA secondary structure prediction studies have revealed two almost identical DB-shaped structures (5Ј-and 3Ј-DBs) in the CR of DENV2 3Ј-UTR (45,46,67) (Fig. 1, A and  B). These structural predictions and phylogenetic analyses were based on the 3Ј-UTR sequences alone from RNA folding algorithms. However, in previous studies using a subgenomic RNA of 719 nt and the viral replicase complex from DENV2-infected mammalian cells or purified NS5, we showed that long range interactions between 5Ј-and 3Ј-terminal regions of DENV2 RNA are important for RNA synthesis in vitro. For example, the 3Ј-UTR itself is not an active template for RNA synthesis in vitro by the viral polymerase unless the 5Ј-UTR is also added in trans or is present in the same RNA (10,16,41,59). Therefore, in this study, we used the "minigenome," a subgenomic RNA sequence containing both 5Ј-and 3Ј-terminal sequences, for  (53), we have replaced the structural proteins with Rluc and EMCV IRES, retaining the first 75 nt of C containing the cHP and the 5Ј-CS, and the last 73 amino acids of E for proper endoplasmic reticulum translocation of NS1. Mutant replicons were made using this backbone for yeast recombination. VR, variable region. B, the secondary structure analysis of the core region of the 3Ј-UTR showed two DB structures (5Ј-and 3Ј-DB) (45). In the left arm of both DBs are identical pentanucleotide sequences, TL1 (nt 10,474 -10,478) and TL2 (nt 10,562-10,566), which are predicted to base-pair with downstream sequences PK2 (nt 10530 -10534) and PK1 (nt 10617-10621), respectively, to form two pseudoknot structures, 5Ј-and 3Ј- (45). On the right arm of the DBs are duplicated conserved sequences RCS2 and CS2. C, a minigenome of DENV2 RNA sequence containing both 5Ј-and 3Ј-UTR sequences and other essential elements for RNA synthesis in vitro (10,54) was used for analysis of secondary structures by MPGAfold as described under "Experimental Procedures." Based on this analysis, a stable best fit structure (E ϭ Ϫ226.2 kcal/mol) has only the 3Ј-DB, whereas the appearance of the 5Ј-DB is visible in metastable structures such as in D (Ϫ225.6 kcal/mol). The frequency of occurrence of individual stems among the final predicted structures is color-coded, which increases from left to right (shown in the bottom left). secondary structure prediction analysis using the MPGAfold algorithm. The MPGAfold algorithm (47,50,56,57), starting from a pool of randomly generated stem-loop structures from a given RNA sequence, recombines these parent structures over multiple generations until the population of RNA molecules evolves and reaches a consensus RNA structure. Because of the stochastic nature of MPGAfold, multiple independent runs must be done to determine the overall consensus structure for the given RNA. Fig. 1C shows the best fit structure (Ϫ226.2 kcal/mol) from 100 different runs of the 719-nt DENV2 WT RNA sequence. The frequency of occurrence of each local structure is shown by a color scheme. The appearance of the already identified and functionally characterized DENV2 RNA secondary structures like the SLA, SLB (16,27,28), cHP (29,30), 3Ј-SL (33,35,37,68), and 5Ј-3Ј-CS1 base pairing (10 -12, 14, 40, 41, 59) confirms the reliability of the MPGAfold-predicted RNA structures.
MPGAfold shows the formation of the 3Ј-DB with 97% frequency of occurrence ( Fig. 1C and supplemental Table 2). However, the sequences involved in the proposed 5Ј-DB do not readily form the expected DB structure; instead, its upstream half, including TL1, is involved in a long stem-loop structure, whereas the lower half folds into another stem-loop with a frequency similar to that of the 3Ј-DB (Ͼ90%) (Fig. 1C).
In some runs, we noticed that the same RNA did form both 5Ј-and 3Ј-DBs as an energetically close suboptimal structure (Fig. 1D, Ϫ225.6 kcal/mol). The frequency of occurrence of the 5Ј-DB among the final structures reached only 6%. However, in an additional 62% of runs, it was found in intermediate structures.
This difference in the folding behavior of the two DBs led us to hypothesize that it is the 3Ј-DB that plays a significant role in the viral life cycle. MPGAfold analysis showed that TL2 is part of an unpaired loop in the 3Ј-DB in the best fit structure, whereas TL1 is alternatively base-paired with upstream sequences (Fig. 1B). Thus, despite their identical sequences of TL1 and TL2, it is TL2 that is available to bind PK1 in the 3Ј-DB to form the putative 3Ј-. This notion is consistent with sequence analysis of different flavivirus 3Ј-UTRs showing that the exact sequence of TL2 is conserved in DENV, WNV, Murray Valley encephalitis virus, and JEV, whereas TL1 is only conserved within the DENVs (45).
Construction and Characterization of the DENV2 Rluc Replicon-To study the effect of different RNA motifs on viral translation and replication, we constructed a DENV2 New Guinea C strain replicon expressing Rluc (Fig. 1A). This was cloned into the yeast shuttle vector pRS424 directly following the SP6 promoter, which was used for in vitro transcription. The viral structural proteins were replaced by the Rluc gene fused with an EMCV IRES element, which directs cap-independent translation of the nonstructural proteins, NS1-NS5. Specific sequences from the N terminus of the capsid-coding region and the last 73 amino acids of the envelope protein were retained because these are required in cis for RNA replication. For specific details of the construction, see "Experimental Procedures." In addition, we constructed the replication-deficient GND replicon by substitution of the highly conserved GDD 3 GND in the polymerase domain (69,70).
Kinetics of the Rluc activity following electroporation of the WT and GND replicon RNAs in BHK21 cells show a peak of Rluc activity at 2 hpt, indicating translation of the input WT and GND replicon RNAs ( Fig. 2A). This activity declined over time, which was then followed by a second increase of luciferase activity in WT replicon up to 96 hpt due to translation of the newly replicated RNA; this second burst of Rluc activity was not observed in the GND mutant. This was confirmed by directly quantifying replicon RNA at 96 hpt through qPCR (Fig. 2B) and by immunofluorescent staining for detection of NS5 (Fig. 2C). These results indicated that the WT and GND replicons could be used in our analysis of the effects of cis-acting elements in translation and replication by quantifying the luciferase activities at 2 and 96 hpt, respectively.
The Roles of Conserved TL and PK Motifs in RNA Translation-To verify whether the RNA folding inequality and the frequency of occurrence between the 5Ј-and the 3Ј-DBs based on the MPGAfold analysis reflect differences in TL1 and TL2 functionalities (Fig. 1, C and D), we analyzed the effects of mutations engineered to disrupt the potential base pair interactions among TL1/PK2 in the 5Ј-DB and TL2/PK1 in the 3Ј-DB (Fig. 1B). We first deleted the conserved sequence of 5 nt in one or both TLs, creating mutants ⌬ 5 TL1, ⌬ 5 TL2, and ⌬ 5 TL1⌬ 5 TL2 (Fig. 3A). The results at 2 hpt indicated that ⌬ 5 TL2 modestly affected translation (72% of WT, p Ͻ 0.0001), whereas deleting TL1 did not have any effect (Fig. 3B). However, when both were deleted in tandem, translation further decreased to 40% of WT (p Ͻ 0.0001), suggesting that TL1 exhibits a cooperative synergy with TL2 in enhancement of translation of DENV2 RNA.
To further investigate the role of TL1 and TL2, we employed a second approach to disrupt the TL1/PK2 and TL2/PK1 base pair interactions. We introduced point mutations by "flipping" the TL sequences (5Ј-GCUGU-3Ј 3 5Ј-UGUCG-3Ј), generating TL1 Flip , TL2 Flip , and TL1 Flip TL2 Flip (Fig. 3A). At 2 hpt, TL1 Flip was efficiently translated at levels similar to WT, confirming the results obtained with ⌬ 5 TL1 (Fig. 3B, leftmost  panel). In the TL2 Flip mutant, translation was reduced to 80% of WT (p ϭ 0.0006) and 33% in TL1 Flip TL2 Flip (p Ͻ 0.0001). The results obtained with Flip mutants thus confirm our results obtained with the deletion mutants described above. Next, we examined whether the reduction in translation of some of the mutant RNAs is due to differences in RNA stabilities. RNA stabilities of the WT and mutants were determined by qPCR. The results showed that there was no appreciable change in RNA stabilities of these mutants from that of the WT (Fig. 3C).
We sought to examine whether the putative partners of TL1 and TL2 would also have a similar effect on translation. Conserved sequences PK2 and PK1 were flipped to make PK2 Flip and PK1 Flip . These mutations unexpectedly did not affect translation (Fig. 3B, second panel). Another set of PK2 and PK1 mutants, PK2 mut and PK1 mut (Fig. 3A), were constructed to verify this result. Both translated efficiently like WT, PK2 Flip , and PK1 Flip (Fig. 3B, rightmost panel). To eliminate the possibility that a mutation in one PK might be rescued by the other non-mutated PK, we mutated both PK2 and PK1 in the same replicon RNA (PK2 mut PK1 mut ). Fig. 3B shows that PK2 mut PK1 mut still translated as efficiently as WT.
Interestingly, the restoration of the base pairing between the TLs and PKs in TL2 Flip PK1 Flip and TL1 Flip PK2 Flip TL2 Flip PK1 Flip having the 20 and 67% defect in translation, respectively, with nonviral sequences (Fig. 3A) restored translation back to WT levels (Fig. 3B, second panel). Similar restoration of translation to WT levels was observed in another set of mutants that swapped the positions of TL2 and PK1 sequences (TL2 3 PK1/PK1 3 TL2 in Fig. 3B, third panel).
Taken together, our evidence suggests that TL1 and TL2 play an important role in translation with unequal contributions and that their specific sequences are not required. Substitutions with non-viral sequences can still maintain efficient translation at WT levels provided that they retain base pairing with PK2 and PK1.
The TLs and PKs Are Critical for Viral Replication-Replication of WT and mutant replicons was assayed by measuring the luciferase signal at 96 hpt. The results indicated that replication of all of the mutants was attenuated (Fig. 3, D and E). TL1 and TL2 have different contributions to replication. A deletion or substitution mutation of the TL1 motif resulted in a 70 -80% decrease in Rluc signal (Fig. 3D) and a 50% drop in RNA levels (Fig. 3E). Mutation of TL2, on the other hand, had a much greater effect, with replication at ϳ10% of WT. Mutating both sequences in one RNA generated a more severe phenotype. In contrast to their phenotypes in translation, PK2 and PK1 mutants are defective in replication (Fig. 3, D and E).
We also sought to examine whether restoring the base pairing in TL and PK mutants with non-viral sequences could rescue the FIGURE 2. Characterization of DENV2 replicon encoding Rluc reporter. A, the kinetics of luciferase expression of WT and a replication-deficient mutant (GDD 3 GND) were followed at different time points post-transfection (hpt). 2 and 96 hpt were determined to be the measure of translation and replication, respectively, of input RNA electroporated into BHK21 cells. RLU, relative light units. B, direct quantification of replicated RNA by quantitative RT-PCR. Quantitative RT-PCR was carried out by RNA extraction from DENV2 replicon-transfected BHK21 cells at 96 hpt as described under "Experimental Procedures." C, detection of DENV2 NS5 by immunofluorescence staining for NS5. Expression of viral proteins was confirmed by immunofluorescence for detection of DENV2 NS5 in BHK21 cells transfected with DENV2 replicon as described under "Experimental Procedures." Cells transfected with infectious RNA of DENV2 were used for a positive control (DENV i.c.). Error bars, S.E. defect in replication. Unlike in translation, the base pairing between mutant TL and PK motifs did not restore RNA synthesis (TL1 Flip PK2 Flip , TL2 Flip PK1 Flip , TL1 Flip PK2 Flip TL2 Flip PK1 Flip , and TL2 3 PK1/PK1 3 TL2 in Fig. 3, D and E). Because PK1 overlaps with the 3Ј-CS1 (Fig. 1B), we sought to determine whether the loss of replication was due to the decrease in base pairing between the cyclization sequences. We introduced compensatory mutations in the 5Ј-CS to restore the loss of base pairing between the 5Ј-3Ј-CS1 in PK1 Flip and TL2 Flip PK1 Flip (indicated by "-3nt"). MPGAfold analysis showed that the frequency of occurrence of the 5Ј-3Ј-CS1 base pairing is reduced as expected in the PK1 Flip mutant because of a reduction in the number of base pairs between the cyclization sequences. This frequency was restored in the PK1 Flip -3nt and the TL2 Flip PK1 Flip -3nt mutants (supplemental Table 2). These compensatory mutations for restoration of the 5Ј-3Ј-CS1 base pairing in these replicons did not improve their replication (Fig. 3, D and E), suggesting that specific sequences within the 5Ј-3Ј-CS1 are important for replication.
The Conserved Nucleotides Surrounding TL2 Motif Also Play a Role in Translation and Replication-Our data in Fig. 3B indicate the possibility that the TL2/PK1 interaction to form the putative 3Јmight be playing a role in stabilizing a structure suitable for translation. We hypothesized that if we increase the strength of the 3Ј-by increasing the number of base pairs, translation efficiency could be enhanced. We chose to mutate the nucleotides adjacent to TL2 that are also conserved and yet are not predicted to participate in the 3Јbase pairing (45,46) and have not been functionally studied previously.
We generated various mutants that increased the number of base pairs by 1-4 bp in the region surrounding TL2 or PK1 (Fig.  4A). Mutating the A immediately downstream or upstream of TL2 (ϩ1 down or ϩ1 up) while increasing the base pairing to 6 nucleotides has a dramatic effect on translation (45 or 60% of WT, respectively; Fig. 4B). Mutations of other residues resulting in 7-9 bp overall had negative effects on translation except for one mutant containing 7 continuous bp, which did not show an appreciable effect (ϩ2 up in Fig. 4, A and B) (p ϭ 0.1118). The data collectively suggest that the conserved nucleotides adjacent to TL2 are also important for efficient translation.
In contrast to their phenotypes in translation, these mutants exhibit a more pronounced defect in replication (Fig. 4C). It is interesting to note that even one base change (ϩ1 down) can result in a 90% reduction in RNA synthesis. The critical roles of these nucleotides in translation and replication thus explain why these sequences are conserved (45,46). The CS2 Sequence, but Not Secondary Structure, Determines Its Role for Its Functionality-As a part of our overall objective to study the role of 3Ј-DB in translation and replication, we examined the contribution of the CS2 region within the 3Ј-DB in these two processes. One previous study done using WNV replicon reported that the entire deletion of CS2 resulted in a significant impairment of replication (12). In this study, we sought to determine the role of CS2 by creating substitution mutations without altering its secondary structure (CS2 mut ) (Fig. 5A) as predicted by MPGAfold (data not shown). To determine which part of the CS2 structure or sequence is important, we made three additional replicons that introduced mutations in the 3Ј-DB junction (Loop A mut ), stem (Stem mut ), and loop (Loop B mut ) (Fig. 5A). All three subdomains affected both translation and replication to a similar extent (Fig. 5, B and C). The combination of all three mutations did not have an additive effect, suggesting that all subdomains most likely functioned as one structural motif that was essential for CS2 function.
Cis-Acting Elements in the Core Region Are Involved in a Noncanonical Cap-independent Translation-It was reported that DENV2 RNA required both 5Ј-and 3Ј-UTR to initiate translation under conditions of inhibition of translation initiation factor 4E by the addition of LY294002 (58). We sought to determine whether the RNA elements within the 3Ј-DB are required for the viral RNA to be translated by a noncanonical mecha-nism. This compound inhibits the phosphotidylinositol-3 kinase, which leads to the hypophosphorylation of 4E-BP1 (eIF4E-binding protein 1) (71). This form of 4E-BP1 is able to sequester eIF4E, thus suppressing cap-dependent translation of the majority of capped mRNAs in the compound-treated host cells. We selected mutant replicons that would give a good representation of all of the mutants we have done so far: ⌬ 5 TL1⌬ 5 TL2, TL1 Flip TL2 Flip , TL1 Flip PK2 Flip TL2 Flip PK1 Flip , PK2 mut PK1 mut , ϩ4 mut, and CS2 mut . The mRNA containing the 5Ј-and 3Ј-UTR of the ␣-globin gene flanking a firefly luciferase coding sequence with an engineered poly(A) tail (15) was used as an internal control for the transfections of WT and mutant Rluc replicons. The time course of this control RNA was followed in the presence and absence of LY294002. The addition of LY294002 to the transfected cells inhibited the cap-dependent translation of this control mRNA (Fig. 6A). However, WT replicon RNA was translated efficiently in the presence of LY294002 (Fig. 6B, WT) (58). We consistently observed an increase in the Rluc signal obtained from the translation of WT replicon RNA in the presence of the LY294002 compared with the signal obtained in the absence of the drug. Moreover, the 5Јand 3Ј-DB mutants that showed deficient translation in the absence of the drug LY294002 were still defective in the presence of the drug (Fig. 6B). However, the cap-dependent translation of the control RNA GLGpA (15) included in all transfec- tion experiments was inhibited in the presence of LY294002 (Fig. 6C), indicating that the drug was active in these assays. These results indicated that the motifs in the CR are required for the noncanonical cap-independent translation of the viral RNA when the cap-dependent translation is inhibited. Interestingly, restoration of base pairing between TK1/PK2 and TL2/ PK1 in the mutant, TL1 Flip PK2 Flip TL2 Flip PK1 Flip , brought translation back to the WT level (Fig. 6B), suggesting that there is a functional interaction between these motifs.
The CR Elements Do Not Affect (Ϫ)-strand Synthesis-To explore the step at which the conserved sequences in the CR act during replication, we tested the efficiencies of RNA synthesis in vitro of (ϩ)-strand RNA templates containing selected mutations using a full-length NS5 with authentic N and C termini. The RdRP assay showed that all mutants produced (Ϫ)-strand efficiently with the exception of PK1 Flip and TL2 Flip PK1 Flip (Fig.  7). This defect in replication was rescued by restoring the 5Ј-3Ј-CS1 complementarity (PK1 Flip -3nt and TL2 Flip PK1 Flip -3nt). These are in contrast to our replicon data that showed that the 5Ј-3Ј-CS1 and CR sequences were required for replication (Fig.  3, D and E). This suggests that these nucleotides are required at subsequent steps of replication.

DISCUSSION
Previous in silico analyses of various DENV2 3Ј-UTR sequences have predicted that the CR folds into two similar DB-like structures, 5Ј-and 3Ј-DB (45,46,67). Our study is unique in that we take into consideration the influence of 5Ј-terminal sequences and structures on the folding of 3Ј-UTR by using the MPGAfold for predicting RNA secondary structures.
Moreover, MPGAfold is able to capture metastable intermediates in an RNA folding pathway. These metastable intermediates as well as stable and best fit structures predicted by MPGAfold have been shown to have biological relevance in a number of previous studies (61)(62)(63)(64)(65)72). Another advantage of MPGAfold is that besides computing the free energy of the RNA, it also provides us with information on the relative frequencies of occurrence of local secondary structural motifs.
Although MPGAfold reveals the formation of RNA secondary structural elements, such as SLA, SLB, cHP, 5Ј-CS/3Ј-CS1 cyclization, 3Ј-DB, and 3Ј-SL readily (Fig. 1C), the 5Ј-DB and 5Ј-3Ј-UAR interaction is not seen in this stable structure. Surprisingly, the base pairing between 5Ј-and 3Ј-UAR elements occurs only 8% of the time in the final structure with an additional 27% in the intermediate structures (supplemental Fig. 1 and supplemental Table 2). This is consistent with the observation that the 5Ј-3Ј-UAR binding occurs only after the 5Ј-3Ј-CS1 circularization has been established (73), suggesting that the latter is the primary driving force in genomic end-to-end communication. This conclusion is also supported by a stem trace plot of the minigenome RNA maturation in an MPGAfold run (data not shown). However, there is strong experimental evidence that 5Ј-3Ј-UAR base pairing is essential for relieving the repression of RNA synthesis by RdRP in vitro by 3Ј-SL (74). The 3Ј-SL includes a part of 3Ј-UAR element, which upon a conformational change in the presence of the viral polymerase is likely to form the duplex 5Ј-3Ј-UAR (74). Conformational changes occurring within the 3Ј-SL due to 5Ј-3Ј-terminal RNA-RNA interaction as well as the NS5/RNA interactions by structure probing and footprinting methods have also been reported for WNV (75). Thus, these results in light of MPGAfold analysis suggest that the 5Ј-UAR forms part of a long stem loop and that the 3Ј-UAR is a part of the 3Ј-SL in the best fit structure (Fig.  1C); the 5Ј-3Ј-UAR duplex, which occurs only in a metastable structure (supplemental Fig. 1 and supplemental Table 2),

FIGURE 5. Effect of mutations in the conserved CS2 region of 3-DB.
A, subsections of the CS2 region were mutated such that the predicted secondary structure is conserved. Bases in lowercase type denote mutation. Translation (B) and replication (C) of WT and mutant replicons transfected into BHK21 cells were measured at 2 and 96 hpt, respectively. Replication measured by Rluc assay (dark gray) and directly by qPCR (light gray) is shown. Error bars, S.E. *, p Ͻ 0.001 compared with WT.
could possibly be formed upon polymerase binding prior to RNA synthesis.
The previously predicted RNA structures of the CR using 3Ј-UTR sequences showed the formation of both the 5Ј-and 3Ј-DBs (45,46,67). However, the best fit structure produced by MPGAfold after 100 independent runs contained only the 3Ј-DB with a high frequency (Fig. 1C). The simultaneous appearance of both DBs occurs mostly in suboptimal folding intermediates in which the population of RNA structures contains 62% 5Ј-DB (Fig. 1D), although among the final structures, the 5Ј-DB was predicted at a very low frequency (6%). Thus, the long range interaction of the 5Ј-and 3Ј-termini is an important determinant in RNA structure prediction of flaviviruses as revealed in this study. In addition to the folding influence of the 5Ј-terminal sequence that includes 5Ј-CS, the rarity of formation of the 5Ј-DB may also be explained by the nature of the bases in its long arm (Fig. 1B). The base pairing in this arm is interrupted midway by an internal loop, reducing the number of stacked pairs to five. On the other hand, the 3Ј-DB has eight uninterrupted base pairs (Fig. 1B). The position of the break in the 5Ј-DB arm increases the potential for that structure to breathe open and allow its sequences to interact with other nucleotides.
Because most of the final structures show the 5Ј-DB sequences in an alternative stem motif, TL1 remains largely inaccessible to the proposed PK2 interaction to form the 5Ј-.  The stability of the 3Ј-DB structure, however, enables TL2 to be maintained as a single-stranded loop that can readily bind to the PK1 and form the 3Ј-. In addition, the number of base pairs at the stalk of the 3Ј-DB corresponds to roughly a half helix turn, positioning TL2 nearer to PK1, suggesting that the 3Јformation is sterically feasible. The differences in the frequency of formation of the DBs and the accessibility of TL2 for interaction with PK1, respectively, support the notion that the 3Ј-DB elements are more important for the virus. This prediction agrees with the sequence alignment data showing that the 3Ј-DB elements are present in the entire flavivirus genus, whereas the 5Ј-DB is only found in the DENV and JEV serological groups (40,76). Moreover, the TL2 and PK1 are conserved within the flaviviruses, whereas TL1 and PK2 sequences are only conserved within the DENV serotypes (45). Our results show that the contribution of TL1 to translation can only be seen when TL2 is also mutated (Fig. 3B), suggesting that TL1 cooperates with TL2 in translation.
MPGAfold analyses for RNA folding were based on the minigenome, whereas the effects of mutations in translation and replication were analyzed in the context of the luciferase-based replicon RNAs (10,210 nt). Therefore, we sought to assess the potential impact of the luciferase coding region on the structures in the 5Ј-and 3Ј-terminal regions that contain the RNA elements of interest. Due to the computational constraints of the MPGAfold algorithm, we surveyed the full replicon secondary structures with the aid of two fast dynamic programming algorithms: Mfold (60,77) and RNAfold (78). A common selfcontained subdomain (1184 nt in length) bordered by conserved interactions between the 5Ј-end of the Rluc gene and the 3Ј-end of the NS5 gene was identified by both algorithms (data not shown). Independent folds of this subdomain by the two dynamic programming algorithms and MPGAfold indicated a potential for interactions between the 5Ј-end of the Rluc gene and the 5Ј-DB sequences upstream of its head SL (data not shown). The key structural motifs discussed here (SLA, UAR, CS, 5Ј-DB head, 3Ј-DB, and 3Ј-SL) were present in the structures predicted by the dynamic programming algorithms, whereas the full 5Ј-DB was still found in energetically close suboptimal structures, suggesting that the interference of the Rluc gene in the RNA folding pathway is minimal.
Our results of replication assays suggest that mutations of these conserved motifs all affect replication. The replicationdeficient PK1 Flip and TL1 Flip PK2 Flip TL2 Flip PK1 Flip mutations impact mainly the frequency of occurrence of the structural motifs rather than inducing major full secondary structure alterations as seen by MPGAfold and yet function differently (supplemental Fig. 3, A and B). Therefore, there seems to be a strong sequence-specific requirement for their function in replication because these mutations result in attenuation (Figs. 3,  D and E, 4C, and 5C). This sequence specificity for replication is most evident in the TL2 Flip PK1 Flip -3nt mutant because the restoration of the TL2-PK1 and 5Ј-3Ј-CS1 base pairs was not sufficient to rescue the lethal phenotype (Figs. 3, D and E). A similar conclusion was reported in which the 3Ј-CS1 involving the same nucleotides that constitute PK1 was mutated in the WNV replicon (43). Their results showed that specific bases in the 5Јand 3Ј-CS1 are important for replication, regardless of comple-mentarity. This sequence-specific requirement of PK1 in the 3Ј-CS1 and its complementary sequence in the 5Ј-CS in replication is not yet understood. Moreover, the in vitro RdRP assays using (ϩ)-strand templates containing WT and 5Ј-and 3Ј-DB mutant sequences and purified full-length DENV2 NS5 polymerase (54,59) revealed that the minus strand synthesis was not affected by these mutations (Fig. 7). These results suggest that the RNA elements in the 5Ј-and 3Ј-DBs are involved at a subsequent step in viral replication.
Our results, taken together, suggest that TL1 and TL2 have variable modulatory effects in translation. It is noteworthy that the translational defects in TL2 Flip , TL1 Flip TL2 Flip , and TL2 3 PK1 could be rescued by a compensatory mutation in PK2 and PK1, although mutations in PK1 and PK2 motifs per se do not show any effect. One possible explanation is that there are alternative base pair interactions in which TL1 or TL2 could participate when PK1 and PK2 are mutated (supplemental Fig. 2). For example, TL2 could participate in an alternative base pair interaction with a motif within the top loop structure of the 3Ј-SL when PK1 is mutated.
More importantly, the ability of the viral RNA to switch to a noncanonical cap-independent pathway in translation was also restored in the mutant, TL1 Flip PK2 Flip TL2 Flip PK1 Flip (Fig. 6B), supporting the model that there is a functional interaction between these motifs that is important in translation. A previous study examined the folding of the 5Ј-and 3Ј-DB regions of DENV4 3Ј-UTR RNA by RNA probing analysis, which revealed that PK2 and PK1 were susceptible to double strand-specific RNase V1 digestion. However, their binding partners as well as the effect of the 5Ј-terminal region on the folding of the CR were not investigated (79). A more thorough structural analysis is warranted to determine the solution structure of the CR and identify short and long range tertiary interactions between this region and different elements in the genome. Recent advances in structure probing have made it possible to analyze RNA structures in the context of the full viral genome. The high throughput selective 2Ј-hydroxyl acylation analyzed by primer extension has been used to determine the complete structure of the HIV-1 genome (80). A similar method can be applied to DENV to solve the secondary structure of the flaviviral genome and identify the binding partners of the elements in the CR.
Mutations of the surrounding nucleotides of the TL2 loop by increasing the number of base pairs affected translation (Fig. 4,  A and B), suggesting the possibility that increasing the putative base pairing between TL2/PK1 destabilizes the 3ЈDB structure. Another possibility is that these four nucleotides also exhibit a sequence-specific requirement for binding of a protein. CS2 has a similar sequence-specific requirement for translation. It seems that mutation of any subregion results in a similar 50% decrease in translation (Fig. 5, A and B). Because the CS2 sequence is duplicated in the 5Ј-DB as RCS2 (Fig. 1B), it is likely that RCS2 may cooperate with CS2 in translation.
Previous studies revealed that some host proteins bind within the 3Ј-UTR, although the precise roles of these RNAprotein interactions in the virus life cycle are not fully understood at present. For example, EF-1␣ binds to the 3Ј-SL of WNV (81, 82) and DENV (83). The La protein and polypyrimidine tract-binding protein bind to the 3Ј-UTR of DENV4 (83).
The polypyrimidine tract-binding protein also interacts with DENV NS4A protein, and this interaction is required for viral replication (84). The host proteins containing the RNA recognition motif, such as TIAR and TIA-1, were shown to interact with WNV 3Ј-SL of (Ϫ)-strand RNA and increase the efficiency of WNV growth (85). On the other hand, the Y-box-binding protein-1 binds to the 3Ј-SL of DENV RNA and represses translation (86). Only the poly(A)-binding protein, which plays an important role in host mRNA translation, was reported to bind within the DENV CR (87).
There are a number of examples in which the RNA elements within the 3Ј-UTRs of positive strand RNA viruses play an important role in translation. In the Turnip crinkle virus 3Ј-UTR, a tRNA-like structure binds to 80 S ribosome and the 60 S ribosomal subunit (64,65). In addition, in many plant virus RNAs, there are different classes of cap-independent translation elements within the 3Ј-UTRs that take over the function of the 5Ј-cap (for a review, see Ref. 88). Although DENV has a 5Ј-cap for translation of its genome, when the cap-dependent translation is limiting, a similar mechanism may take over under conditions of inhibition of translation initiation factor 4E required for cap-dependent translation, as reported previously (58). The identification of the RNA elements within the CR that seem to be required for both cap-dependent and a noncanonical cap-independent translation warrants further investigation to identify host factors that mediate in this pathway.