Structural basis for the activation of plant bunyavirus replication machinery and its dual-targeted inhibition by ribavirin

Structural basis for the activation of plant bunyavirus replication machinery and its dual-targeted inhibition by ribavirin

Main

The discovery of tobacco mosaic virus in the late nineteenth century led to the recognition of viruses as a new class of pathogen and the foundation of virology as a discipline. Subsequent studies on plant viruses have contributed greatly to the development of modern virology and molecular biology fields. For instance, the demonstration of viral RNA (vRNA) itself being infectious, and the discovery as well as elucidation of RNA interference as an antiviral defence strategy, were seminal in popularizing these disciplines1,2. Plant viruses are known to cause severe quality and yield losses in agricultural crops with an estimated economic impact of over US$30 billion annually, posing a major threat to global food security3,4. At present, planting disease-resistant cultivars and applying insecticides to manage insect vectors are the primary methods employed to prevent viral infections and reduce disease losses. However, for plants that have already been infected by viruses, there are currently no effective antiviral drugs for treatment.

In the past decade, numerous efforts have been made towards understanding polymerase structures and the replication mechanism of human and animal viruses, to help in the design and development of compounds targeting this machinery with high efficacy and exhibiting broad-spectrum antiviral activity5,6,7. The RNA-dependent RNA polymerases (RdRps) from the orders Articulavirales and Bunyavirales have received particular interest because some of their members (that is, influenza virus, Crimean–Congo haemorrhagic fever virus, Rift Valley fever virus (RVFV) and Lassa virus (LASV)) have high epidemic potential and we lack antivirals or vaccines for them. In recent years, these studies have provided structures for the RdRp of influenza virus8,9,10,11,12,13, La Crosse virus (LACV)14,15,16, RVFV17, severe fever with thrombocytopenia syndrome virus (SFTSV)18,19,20, Hantaan virus (HTNV)21,22, Sin Nombre virus (SNV)23, LASV and Machupo virus (MACV); for the latter two, studies have also provided their structures in the presence of their cognate regulatory Z proteins24,25,26,27.

The International Committee on Taxonomy of Viruses lists more than 2,100 different species of plant viruses so far4; to our knowledge, no structure has been resolved for any plant virus replication machinery. Gaps in our understanding of the structure of plant virus replication machinery, which plays a crucial role in the life cycles of viruses, are hampering the development of pesticides against plant viruses. More than 80% of plant viruses are RNA viruses and replicate in the cytoplasm via a viral-encoded RdRp4. Tomato spotted wilt orthotospovirus (TSWV) is a global negative-sense plant RNA virus with a segmented genome in the family Tospoviridae of the order Bunyavirales, and its encoded RdRp (also known as the L protein; ~330 kDa) is one of the largest RdRps of all plant viruses28,29,30,31. The TSWV L protein is also larger than most animal-infecting viral polymerases of known structure (typically ranging from ~50 to 250 kDa)5,32,33,34,35. In this study we present cryogenic electron microscopy (cryo-EM) structures of full-length TSWV RNA polymerase and provide insights into its activation mechanism and its dual-targeted inhibition by ribavirin.

Results

Overall structure

We determined the cryo-EM structure of both the apo L (L protein without RNA) and the RNA-bound L (L protein in complex with the viral genomic RNA) following successful expression and purification of the full-length L protein of TSWV in High Five insect cells. In the apo-L grids, the major class of particles represented the apo-core structure, which consisted of the polymerase core exclusively and was resolved at a resolution of 3.9 Å (Fig. 1a,b and Supplementary Fig. 1). Density for the amino-terminal domain and the carboxy-terminal domain (CTD) was missing in the electron maps, probably due to the intrinsic flexibility associated with these domains. However, the addition of the vRNA promoter (5′ vRNA(1–17 nucleotides (nts)) and 3′ vRNA(1–17 nts)) resulted in an ordered endonuclease (endo) region in the major class identified through 3D classification. Furthermore, the resolution of this structure was improved to 3.7 Å (Fig. 1a,c and Supplementary Fig. 2). The polymerase cores of the L–5′3′-vRNA complex with endo and the L–5′3′-vRNA complex without endo showed no significant differences (Supplementary Fig. 3). When viral 5′ vRNA(1–10 nts) was added, an additional density corresponding to the CTD was observed in only 4% of the particles. After refinement and post-processing, the overall resolution of TSWV L containing the CTD was up to 4.0 Å (Fig. 1a,d and Supplementary Fig. 4).

Fig. 1: Overall structure of TSWV L protein.
Structural basis for the activation of plant bunyavirus replication machinery and its dual-targeted inhibition by ribavirin

a, Domain architecture of the TSWV L protein. NTD, N-terminal domain; vRBL, vRNA binding lobe; CBD, cap-binding domain; Mid domain, middle domain. bd, Cryo-EM maps and cartoon representations of the apo L-core (PDB ID: 8KIA) (b), vRNA-promoter-bound L (PDB ID: 8KI7) (c) and apo L containing the CTD (PDB ID: 8KI9) (d). The structures are coloured by domain with the same colour code as in a. 3′ RNA and 5′ RNA are shown in soft red and bright orange, respectively. The 5′ hook is highlighted in lime green.

Full size image

The structure of TSWV L can be divided into three distinct regions, analogous to the influenza virus polymerase subunits: the N-terminal PA-like region (residues 1–1012), the central PB1-like region (RdRp, residues 1013–1760) and the C-terminal PB2-like region (residues 1761–2880) (Fig. 1a). Despite significant differences in primary sequence and molecular weight, the overall architecture of TSWV L exhibits remarkable similarity to other structurally characterized polymerases of segmented negative strand RNA viruses (sNSVs) (Extended Data Fig. 1). According to the alignment with the influenza polymerase subunits, the TSWV RNA polymerase core, which forms a globular structure, covers the central PB1-like region and the C-terminal part of the PA-like region as well as the N-terminal part of the PB2-like region (Fig. 1a,b). The PB1-like region of TSWV L adopts the characteristic right-hand fold with canonical fingers, thumb and palm subdomains commonly observed in polymerase domains (Extended Data Fig. 1). The endo domain of TSWV L is larger than those of other structurally characterized polymerases (~310 residues versus ~200 residues), primarily due to the presence of additional helices at the N-terminal region of the endo domain. However, all these domains, belonging to the PD-(D/E)-K superfamily of nucleases, share a common two-lobed, kidney-shaped architecture (Extended Data Fig. 2a–e). DALI36 blast showed that the N-lobe composed of five helices (alpha 1–4 and 10) shares highly similar structural features with an iron-binding protein (PDB ID: 6N63) (Extended Data Fig. 2f). The CTD domain of TSWV L is notably more distinct in size and folding, comprising approximately 800 residues, in contrast to the 400–500 residues found in other sNSVs. It consists of three distinct domains: a cap-binding domain, a middle domain and a C-terminal extension (Extended Data Figs. 3 and 4). The endo domain and CTD extend outward from the core and face each other across a solvent channel in the superposed structure (Fig. 1c,d and Supplementary Fig. 5).

Conformational changes upon vRNA promoter binding

The highly conserved and complementary 5′ and 3′ termini of sNSVs genomic elements play crucial roles as promoters for both transcription and replication. Here we observed that the 5′ and 3′ ends form a continuous double-stranded region encompassing both a proximal and a distal duplex in the TSWV L structure (Figs. 1c and 2a), in contrast to the presence of only a distal duplex in both influenza virus and the other bunyavirus polymerases (Extended Data Fig. 5). The proximal duplex is formed through complementary base pairing of nucleotides 4–10 of the 5′ end with nucleotides 4–10 of the 3′ end. The distal duplex is formed by complementary base pairing of nucleotides 11–13 of the 5′ end with nucleotides 11–13 of the 3′ end. Nucleotides 1–3 of the 5′ end are not visible, whereas nucleotides 1–3 of the 3′ end are visible and are directed towards the active sites of the polymerase (Fig. 2a). In addition, nucleotides 1–10 of the 5′ end form a stem–loop (‘hook’) in a hook-binding pocket. This structural arrangement primarily relies on a mismatched base pair (A1–A10) and two canonical base pairs (G2–C9 and A3–U8) (Fig. 2a).

Fig. 2: vRNA promoter binding and structural changes.
figure 2

a, Structure and schematic representation of the vRNA promoter, with base pairings indicated by lines. The proximal and distal duplexes are outlined, while RNA nucleotides that are present but not visible are depicted in grey. b, Global rearrangements of L structure upon vRNA promoter binding. The structures of apo L (light grey; PDB ID: 8KIA) and vRNA-bound L (coloured; PDB ID: 8KI7) are superimposed. Movement and ordered domains are indicated. c, View of the right-hand-like RdRp showing the ordered active site in vRNA-bound L (coloured), which is disordered in the apo form (light grey). The active site is surrounded by a dotted box in the left panel and zoomed in the right panel. Conserved RdRp motifs A–H and the priming loop are shown in bright orange, yellow, slightly green, cyan, violet, brown, magenta, dark grey and blue, respectively. d,e, The active sites of LACV (d) and FluA (e) polymerases upon vRNA promoter binding. f, Comparison of conformational differences in the active site between apo L (light grey) and 5′-hook-bound L (PDB ID: 8KI8). g, The hydrogen bond (dashed yellow lines) between C5 of the 5′ hook and Lys1291 of the fingertips stabilizes motif F (fingertips) of the polymerase active site. h,i, Interaction between motif F and 5′ hook RNA in LACV (h) and FluA (i) polymerases. j,k, Conservation analysis of amino acid Lys1291 in the polymerases of tospoviruses (j) and sNSVs (k). Identical residues are indicated with a red background, and similar residues are depicted in blue boxes with red text. A blue triangle highlights the position of Lys1291. l,m, Inoculation of N. benthamiana leaves with wild-type (WT) TSWV infectious clone (l) or its mutants (m), via agro-infiltration. Left panels of l and m, accumulation of eGFP fluorescence expressed from TSWV in N. benthamiana leaves was detected at 60 h post inoculation. Scale bars, 200 μm. Right panels of l and m, immunoblot analysis of eGFP and L protein expression in the leaves shown in the left panels was performed. M, medium segment; opt, codon-optimized; SR, small segment reporter; EV, empty vector. The RuBisCO large subunit was stained with Ponceau S (PS) to indicate sample loading.

Source data

Full size image

The elucidation of the structure of L in its apo state allowed us to compare the conformational changes between vRNA-promoter-bound L and apo L. Superimposition of the two structures reveals a noticeable global rearrangement. The binding of the vRNA promoter (Supplementary Table 1) induces the movement of multiple subdomains, including the lid, bridge and α-ribbon, resulting in a more compact polymerase core and an ordered endo region (Fig. 2b). Furthermore, the vRNA promoter binding directly stabilizes the RNA-binding regions of TSWV L, including the arch and clamp domain (Fig. 2b). Two functional elements critical for enzyme activity, the fingertips (also known as motif F) and the priming loop, which are disordered in apo L, are stabilized in the vRNA-promoter-bound state (Fig. 2c). Similar to the vRNA -bound polymerases of LACV and FluA (Fig. 2d–e), the stabilized motif F, together with other motifs (motifs A–E, G and H), forms a stable catalytic centre.

5′ hook vRNA is sufficient to activate polymerase

Because nucleotides 1–10 of 5′-end vRNA form a special hook structure, we further prepared a complex of TSWV L and the single-stranded 5′ vRNA(1–10 nts) promoter and determined its structure to validate whether this hook structure alone can activate the L protein. The structure of the 5′-vRNA(1–10 nts)–L complex was obtained at 3.2 Å (Extended Data Fig. 6a).

In 5′-vRNA(1–10 nts)-bound L, both motif F and the priming loop are ordered compared with apo L and have a similar conformation to vRNA-promoter-bound L (Fig. 2f and Extended Data Fig. 6b), suggesting that the hook structure is sufficient to activate the L protein and facilitate the formation of the catalytic centre. The hydrogen bond between C5 of the 5′ end and residue Lys1291 of the fingertips plays a crucial role in the interaction between the RNA hook and motif F (Fig. 2g). Analysis of the interactions in other viral polymerase structures revealed that the residue Lys956 on motif F of the LACV L protein and the residue Arg238 on motif F of the FluA polymerase also form hydrogen bonds with the protruding base U6 at the bend of the hook structure (Fig. 2h,i). Moreover, residue Lys1291 is highly conserved among all orthotospovirus polymerases, but not those of animal viruses (Fig. 2j,k).

To further investigate the impact of the elements discerned in the activated L structure on viral replication, we introduced specific mutations in the TSWV infectious clone. More specifically, we mutated amino acid Lys1291 to Ala in the L segment and changed nucleotide C into A at the fifth position of the 5′ end in the S RNA template. Whereas clear fluorescence of the enhanced green fluorescent protein (eGFP) marker gene was observed in the agro-infiltrated leaves with wild-type TSWV at three days post inoculation (dpi), no fluorescence was observed with TSWV carrying the LK1291A mutation (Fig. 2l). However, a significantly lower level of eGFP fluorescence was observed in the rescued TSWV containing the SC5A mutation (Fig. 2m). The eGFP accumulation and production of L protein in each sample was verified by immunoblotting and indicated that the mutation introduced into the L protein had no significant effect on its stability (Fig. 2l,m). We also mutated other amino acid residues (Tyr455, Lys467, Thr468, Thr653, Lys656, Asp775, Arg823, Pro824, Gln825, Arg826, Leu827, Arg831, Thr895, Lys896, Leu897, Val1014, His1015, Thr1018, Lys1191, and Glu1432) in the L protein that interact with hook RNA in the viral infectious clone and found that all of them are essential for the rescue of TSWV (Extended Data Fig. 7).

Formation of a complete catalytic centre

To investigate how the ordered motif F functions, we further analysed the interaction between motif F and other regions of the TSWV L protein. Motif F is primarily involved in direct polar interactions with the fingers and core lobe regions, including hydrogen bond interaction between conserved residue Arg1294 and conserved residue Gln1455 on motif B (Fig. 3a and Supplementary Table 2). Also, residue Arg1289 interacts with residue Lys1008, which we found to be absolutely conserved in the L proteins of all bunyaviruses; we denoted this area as motif I (Fig. 3a,b).

Fig. 3: Newly identified conserved motif I in sNSV polymerases.
figure 3

a, Hydrogen bond (dashed yellow line) between Arg1289 of ordered fingertips and Lys1008 of the core lobe in the TSWV polymerase active site. b, The conservation of lysine (motif I) is shown in an alignment in bunyaviruses (TSWV (Orthotospovirus); BADUV, Badu virus (Phasivirus); EMARaV, European mountain ash ringspot-associated virus (Emaravirus); FERV, Ferak virus (Orthoferavirus); GOLV, Gouleako virus (Goukovirus); HEBV, Herbert virus (Herbevirus); HTNV (Orthohnatavirus); JONV, Jonchet virus (Orthojonvirus); KIGV, Kigluaik phantom virus (Orthophasmavirus); LACV (Orthobunyavirus); LASV (Mammarenavirus); MACV (Mammarenavirus); RSV, rice stripe virus (Tenuivirus); RVFV (Phlebovirus); SFTSV (Phlebovirus); SNV (Orthohnatavirus)) and orthomyxoviruses (FluA). The blue triangle highlights the position of the conserved lysine. cj, The surroundings of motif I are illustrated for HTNV (c), LACV (d), LASV (e), MACV (f), RVFV (g), SFTSV (h), SNV (i) and FluA (j) with the neighbouring motifs A–H. k, N. benthamiana leaves were inoculated with wild-type TSWV infectious clone or its mutant via agro-infiltration. Accumulation of eGFP fluorescence expressed from rescued TSWV was detected in N. benthamiana leaves at 60 h post inoculation. Scale bars, 200 μm (left). Expression levels of eGFP, L and its mutant were analysed via western blot (top right), and the S genomic RNA (red arrowhead) and eGFP mRNA transcripts (green arrowhead) were analysed via northern blot using DIG-labelled antisense eGFP probes (bottom right). Ethidium bromide staining of ribosomal RNA (rRNA) was used as the RNA loading control.

Source data

Full size image

To verify the spatial position of motif I in the polymerase structures of other sNSVs, we analysed the active sites of all structurally resolved polymerases. As shown in Fig. 4c–j, motif I is consistently observed in the polymerase active site across various viruses. In traditional bunyaviruses such as HTNV, LACV, RVFV, SFTSV, SNV and TSWV, the conserved Lys residue resides on a flexible loop, proximal to both motif F and motif G (Fig. 3c,d,g–i). For arenaviruses (newly classified in the order Bunyavirales), including LASV and MACV, the conserved Lys residue is located on an α-helix and is close to motif A, motif H and motif F (Fig. 3e,f). In the case of FluA, the spatial characteristics of the Thr42 amino acid are similar to those of arenaviruses (Fig. 3j), indicating that motif I in influenza viruses and arenaviruses may be evolutionarily closer.

Fig. 4: Mechanism of TSWV inhibition by ribavirin.
figure 4

a, Effects of ribavirin treatment on TSWV infection. N. benthamiana plants inoculated with TSWV were subjected to treatment with either buffer or ribavirin (100 μg ml−1). Left: disease development was photographed at 7, 14 and 21 dpi. Right: the accumulation level of nucleocapsid protein (NP) in the systemic leaves of treated plants was determined via immunoblot analysis using an NP-specific antibody. CK+, TSWV-infected sample. b, Left: cut-away global view of the L-ribavirin complex (PDB ID: 8KI6). The upper right panel provides a close-up and sliced-surface view of ribavirin located in the middle of the 5′-hook-binding pocket and the template-binding channel. The lower right panel offers a magnified view of ribavirin positioned in the active site. Ribavirin is shown as a blue shape. c, vRNA structures within the 5′-hook-binding pocket and template-binding channel of the L–ribavirin complex, achieved by superimposing the vRNA-bound L structure onto the ribavirin-bound L structure described here. The two surface views display the location of ribavirin in vRNA-binding position. d, Illustration of the entrance and exit tunnels and active site of L depicting how ribavirin inhibits activity of the polymerase. The L structure (grey) with vRNA (coloured) is displayed with tunnels (cyan) labelled with arrows as template entry, NTP entry, product and template exit tunnel, as calculated using Caver v.3.0.3 (left). Ribavirin is superimposed on the same tunnels as in the left panel and is depicted in blue (right). e, 5′ hook RNA binding analysis of the L protein in the absence or presence of ribavirin by EMSA. Increasing concentrations of L protein were incubated with FAM-labelled 5′ hook RNA(1–10 nts) in the absence of ribavirin (lanes 1–4), and the L and FAM-labelled 5′ hook RNA mixtures were incubated with decreasing amounts of ribavirin (lanes 5–8). f, Effect of ribavirin on in vitro RNA synthesis directed by TSWV RNP. Digoxin-labelled RNA synthesized in vitro by heat-inactivated virus (lane 2) and untreated virus in the absence (lane 1) and presence (lanes 3–5) of increasing concentrations of ribavirin (25, 50 and 100 μg ml−1, respectively).

Source data

Full size image

To explore the biological role of motif I, we evaluated the effect of the residue Lys1008 on virus replication using TSWV infectious clones. Wild-type and mutant-infected clones were infiltrated into Nicotiana benthamiana leaves. eGFP fluorescence was observed in leaves infiltrated with the wild type at 3 dpi, but no fluorescence was observed in leaves infiltrated with the clone carrying the LK1008A mutation (Fig. 3k). Protein immunoblot analysis indicated that both wild-type and mutant L protein were expressed normally, and the accumulation of eGFP protein was consistent with fluorescence intensity (Fig. 3k). These results confirm the crucial role of Lys1008 in rescuing the virus.

Dual-targeted inhibition of polymerase activation

Currently, there are no chemical antiviral compounds directly targeting plant viruses. The structural similarities between the RdRp structure of TSWV and those of the animal/human-infecting viruses of the Bunyavirales and influenza viruses suggest that human antiviral drugs targeting RdRp may also be effective against TSWV. To test this, we used the TSWV infectious clones to screen the inhibitory effects of three well-known human antiviral drugs—namely, remdesivir, ribavirin and favipiravir—on viral replication. Infectious clones were co-infiltrated with the drugs in N. benthamiana leaves, and eGFP fluorescence was examined at 3 dpi. Leaves treated with remdesivir revealed a significant reduction in eGFP expression, while leaves treated with ribavirin showed no fluorescence, indicating that both drugs, in particular ribavirin, were effective against TSWV replication (Extended Data Fig. 8). Although a small amount of fluorescence was observed in the leaves treated with favipiravir, the leaves exhibited necrosis, indicating potential phytotoxicity.

To further substantiate these findings, we also tested the effect of ribavirin on a normal viral infection. N. benthamiana leaves were mechanically inoculated with TSWV in the absence (buffer control) and presence of a ribavirin solution at a concentration of 100 μg ml−1. In the buffer-control-treated group, typical symptoms of a viral infection, including leaf curling and stunting, were observed in the systemic (newly emerged) leaves at 14 and 21 dpi. The presence of TSWV in the systemic leaves was confirmed by immunoblotting assays from as early as 7 dpi. In contrast, in the ribavirin-treated group, no typical symptoms of viral infection were observed in the systemic leaves at any time point, and the absence of TSWV in ribavirin-treated plants was confirmed by immunoblotting (Fig. 4a). These findings indicate that ribavirin strongly inhibits TSWV infection, which is consistent with previous reports37,38. Given that ribavirin is a nucleoside analogue, we hypothesized that its antiviral activity may be attributed to its interaction with the TSWV polymerase, potentially disrupting viral genomic RNA synthesis. To investigate this and identify the binding site of the nucleoside analogue, we mixed purified TSWV polymerase with ribavirin and determined the complex structure at atomic resolution (Supplementary Fig. 6).

Unexpectedly, compared with the density map of L protein in the apo state, the complex revealed two additional densities for L–ribavirin (Supplementary Fig. 7), indicating that ribavirin targets two sites on TSWV L. Ribavirin binds not only to the active site of the polymerase but also to a key channel located in the RNA-binding site (Fig. 4b). By modelling the viral promoter structure in the complex and calculating the internal tunnel of the L core, we found that ribavirin occupies two key sites that effectively block both the active catalytic formation and the entry of the 5′ end into the hook-binding pocket required for the formation of a stem–loop (hook) structure (Fig. 4c,d). Given the previously mentioned critical role of the hook structure in activating the L protein, ribavirin binding at this channel position prevents the allosteric activation of polymerase.

We also performed an electrophoretic mobility shift assay (EMSA) to investigate the binding affinity change between the L protein and 5′ hook RNA in the presence or absence of ribavirin. As shown in Fig. 4e, a fixed amount of FAM (fluorescein)-labelled 5′ hook RNA (1.5 pmol) was incubated with increasing concentrations of TSWV L protein. The free RNA probe was able to completely shift in the presence of 4 pmol of L protein. Upon the addition of ribavirin, RNA shift was significantly reduced, suggesting that ribavirin is able to prevent the L protein from binding 5′ hook RNA (Fig. 4e). Furthermore, we investigated the inhibition activity of ribavirin by an in vitro RdRp activity experiment as reported in ref. 39. Purified TSWV ribonucleoprotein complex (RNP) was incubated in the presence of NTPs (including digoxin-labelled UTP) to synthesize RNA in vitro. The digoxin-labelled RNA synthesized in vitro by heat-inactivated TSWV RNP and untreated RNP in the absence and presence of ribavirin were resolved by RNA gel electrophoresis and further analysed by northern blot. High amounts of RNA product were observed in the absence of ribavirin, while the total amount of RNA synthesis was decreased considerably in the presence of ribavirin (Fig. 4f). These data suggest that ribavirin can inhibit the RNA synthesis activity of polymerase in vitro.

Because ribavirin is a prodrug, it would be activated by intracellular phosphorylation into mono-, di- and triphosphates. The inhibition mechanism of ribavirin against TSWV in N. benthamiana leaves may depend on the activated form of ribavirin. We further tested the effect of ribavirin triphosphate treatment on TSWV infection in planta. N. benthamiana leaves were mechanically inoculated with TSWV in the absence (buffer control) and presence of a ribavirin triphosphate solution at a concentration of 100 μg ml−1. Typical symptoms of viral infection were observed in the buffer-control-treated plants, whereas no typical symptoms appeared in the ribavirin-triphosphate-treated plants at 14 dpi (Fig. 5a). These data suggest that ribavirin triphosphate also strongly inhibits TSWV infection. In addition, we prepared a complex of TSWV L and ribavirin triphosphate and determined its structure at 3.0 Å (Supplementary Fig. 8). Similar to the dual-target binding pattern of ribavirin on TSWV L protein, ribavirin triphosphate binds to both the active site of the polymerase and the entry channel of 5′ hook RNA. Ribavirin triphosphate located in the entry channel forms hydrogen bonds with residues Lys1191, Met777 and Thr779, while active-site-located ribavirin triphosphate makes hydrogen bonds with Tyr908, Lys1540, Lys1541, Asn1538, Lys1370 and Gln1455 (Fig. 5b).

Fig. 5: The inhibition of TSWV L by ribavirin triphosphate.
figure 5

a, Effects of ribavirin triphosphate treatment on TSWV infection. N. benthamiana plants inoculated with TSWV were subjected to treatment with either buffer or ribavirin triphosphate (100 μg ml−1). Left: disease development was photographed at 14 dpi. Right: the accumulation level of TSWV NP protein in the systemic leaves of treated plants was determined via immunoblot analysis using an NP-protein-specific antibody. Ponceau S staining of the RuBisCO large subunit was used to indicate sample loading. b, Left: cut-away global view of the TSWV L–ribavirin triphosphate complex (PDB ID: 9J8V). The L protein is coloured by domain with the same colour code as in Fig. 1. Ribavirin triphosphate is represented as a blue shape with a dashed outline. The right panel shows detailed interactions between L and ribavirin triphosphate. The hydrogen bonds are shown as yellow dashed lines.

Source data

Full size image

Discussion

The structure of TSWV L reported here reveals the structural features of an RNA polymerase from a plant virus belonging to the sNSVs with commonalities to animal-infecting viruses. Polymerases from picornavirus to influenza virus typically range from ~50 to 250 kDa; therefore, the ~330-kDa TSWV L represents the largest polymerase in terms of molecular weight among all resolved RNA polymerase structures to date. The overall architecture of the TSWV L protein consists of three main components: the N-terminal endo fusing a metal-binding protein-like domain, the central core region containing the RdRp polymerase and the variable CTD comprising the cap-binding domain region (Fig. 1 and Extended Data Fig. 1). This implies a common evolutionary origin between plant and animal RNA polymerases of sNSVs. RNA polymerases larger than TSWV L (for example, the ~450-kDa Nairovirus L protein) may also have a similar overall architecture.

The structure of the viral promoter has been controversial in the literature, with various models such as the ‘panhandle’, ‘fork’, ‘corkscrew’ or ‘hook’ proposed on the basis of structural and functional studies of influenza viruses40. The current paradigm of sNSVs polymerase structures is that the vRNA promoter binds to the polymerase in a hook structure, supported by the presence of a stem–loop (hook) structure formed by the 5′ proximal strand and a distal duplex formed by the 5′ and 3′ distal strands5. However, it is still unclear whether proximal 5′ and 3′ single-stranded RNA ends can form a proximal duplex and assemble into a panhandle-like structure. The two single-stranded terminal ends of viral genomic elements, making up the panhandle structure and comprising the promoter of bunyaviruses, are more complementary than those of influenza virus6. In this study, we observed a proximal double-stranded region formed in the viral promoter on the TSWV L protein (Figs. 1c and 2a) that is not observed in influenza virus, whereas a distal duplex is observed in both influenza virus and bunyavirus polymerases8,11,14,20,21,24 (Extended Data Fig. 5). Although the presence of a duplex consisting of the proximal nucleotides 1–10 of the 5′ end with the 3′ end and simultaneous formation of a hook structure in the 5′ end is not an actual state, it still reveals the possibility of the vRNA promoter binding in a fully double-stranded conformation at the template entrance channel of the polymerase. We propose that the newly synthesized L protein (apo L) in TSWV-infected plants is in an inactive state and initially binds to genomic RNA with a fully double-stranded end (Figs. 2a and 6). The binding of the vRNA promoter allows the 5′ end to enter the hook-binding pocket, where it folds into a stabilized hook structure and simultaneously activates the polymerase’s active site, together with ordered endo and CTD to engage it in replication and transcription processes (Figs. 1 and 6). Upon treatment of TSWV-infected plants with ribavirin, ribavirin or its triphosphate form binds to both the 5′-hook-binding pocket entry tunnel and the enzyme’s active site, blocking both vRNA recognition and the formation of a complete catalytic centre, two key steps in polymerase activation, and ultimately locking the polymerase in an inactivated state (Figs. 4–6).

Fig. 6: Structure-based model of TSWV L activation and its dual-targeted inhibition mechanisms.
figure 6

In TSWV-infected plants, the newly synthesized L protein (apo L) remains in an inactivated state with a disordered active site. When the viral genome with a double-stranded terminus binds to the L protein, the proximal double strand unwinds; the 3′ end becomes single-stranded and enters the polymerase active centre, and nucleotides 1–10 of the 5′ end enter the hook-binding pocket to form a hook-like stem–loop. The protruding base C5 at the bend of the hook structure makes hydrogen bonds with Lys1291 of the fingertips (motif F), resulting in an ordered motif F. New hydrogen bonds (Arg1294–Gln1455 and Arg1289–Lys1008) between motif F and other motifs contribute to the formation of a complete polymerase active centre. These alterations in the active site together with the ordered endo and CTD engage TSWV L in an activated state to initiate replication and transcription processes. In ribavirin-treated plants, ribavirin or its triphosphate disrupts both vRNA recognition and the formation of a complete catalytic centre by occupying both the 5′-hook-binding pocket entry tunnel and the enzyme’s active centre, locking the polymerase in an inactivated state.

Full size image

A set of eight structural motifs (A to H) of RdRp play a key role in catalysing RNA synthesis. Motif F comprises positively charged residues that interact with the phosphate group of incoming NTP. Through the interactions between positively charged residues on motif F and nucleotides on the hook RNA in TSWV polymerase (Fig. 2g), another potential role of motif F emerges—it might function as a ‘sensor’ to detect the unique hook structure formed by viral genomic RNA. This recognition is a critical step for polymerase activation. Furthermore, to our knowledge, no new motifs have been reported since the identification of two previously unknown conserved motifs, G and H, specific to all sNSV polymerases by Gerlach et al. in 201514. In this study, while analysing the intramolecular interaction of the TSWV L protein with motif F, we unexpectedly discovered an absolutely conserved lysine in the bunyavirus L protein (Fig. 3a,b). Further structural alignment revealed that this conserved lysine is present in the active site of all sNSV polymerases (Fig. 4c–j), indicating the identification of a new motif, denoted I, which provides deeper insights into the composition of the active site of the sNSV polymerases. Following the recognition of vRNA by the disordered motif F, the region becomes ordered, taking on the role of an ‘adaptor’, and generates new interactions with other residues (including motif I), thereby placing the polymerase in an activated state. Moreover, the Lys1008 residue makes an ionic contact with the product strand phosphate group at the +3 position from the active site when we model a duplex RNA into the active site. This traps the tip of the motif F finger sequence underneath itself, fully structuring motif F for NTP binding and activity (Extended Data Fig. 9). Therefore, another key mechanistic insight may be that this ordering of motif F is probably dependent on first making a ≥3-base-pair product RNA to establish the Lys1008–phosphate contact, and this may in turn trigger a transition to a stable elongation complex.

In contrast to the rapid development of agricultural fungicides, insecticides and herbicides, there remains a notable gap in antiviral pesticides. To investigate the occurrence of similar interactions between viral polymerases and ribavirin or its triphosphate, we searched the PDB database and found structures of ribavirin complexes with the polymerases of noroviruses and aphthoviruses41,42, both of which are positive-stranded RNA viruses. In these structures, ribavirin or its triphosphate binds to a single target site within the catalytic core of the polymerase (Extended Data Fig. 10). The structures of the TSWV L–ribavirin and TSWV L–ribavirin triphosphate complexes thus reveal unique dual target sites of the antiviral nucleoside analogue. In addition to occupying the active site, they can bind at the 5′-hook-binding pocket entry tunnel of TSWV L, preventing the 5′ end from entering the hook-binding pocket to form a hook structure to activate the polymerase (Fig. 4b–d).

Because the components of human and veterinary drugs are now banned from use as pesticides to avoid the emerging risk of resistance, ribavirin is not currently available for direct application to crops, and therefore plant-virus-specific drugs need to be developed. Considering that interaction between the Lys amino acid at position 1291 of the L protein and the C base at the fifth position of the 5′ end hook is crucial for viral infection and is highly conserved in orthotospoviruses but not conserved in human and animal sNSVs (Fig. 2j,k), the Lys1291–C5 interaction site could serve as a specific target for designing pesticides against plant orthotospoviruses. In addition, since the formation of the 5′ end hook structure is a common feature in the initiation stage of sNSV polymerases, the polymerase hook formation tunnel presents a highly attractive target for the design of broad-spectrum antiviral drugs against sNSVs.

Methods

Protein expression and purification

The codon-optimized sequence for the L protein (GenBank: JF960237.1) of the TSWV-YN isolate was synthesized and cloned into a pFastBac expression plasmid with an N-terminal 10×His–Msyb tag. Recombinant baculoviruses were prepared using SF9 insect cells (11496015, Invitrogen), and protein expression was carried out in High Five cells (B85502, Invitrogen) via the Bac-to-Bac expression system. At 60 h post inoculation, the culture was centrifuged to pellet the cells, which were subsequently resuspended in chilled lysis buffer (50 mM Tris (pH 8.0), 200 mM NaCl, 10 mM imidazole) and then lysed by sonication. After centrifugation at 14,000 g for 1 h at 4 °C, the supernatant was loaded twice onto a Ni-affinity column (GE) pre-equilibrated with lysis buffer. The resin was washed three times with 120 ml of lysis buffer, and the bound proteins were eluted with elution buffer containing imidazole (50 mM Tris (pH 8.0), 200 mM NaCl, 300 mM imidazole) and then captured onto a 5-ml Q column (GE Healthcare). The proteins bound to the column were eluted with a linear concentration gradient of NaCl, and the fractions that peaked at 380 mM NaCl were pooled together. Subsequently, this pooled fraction was loaded onto a Superose 6 10/300 GL column (GE Healthcare) and eluted in 50 mM Tris (pH 8.0), 200 mM NaCl, 1 mM Tris (2-carboxyethyl) phosphine buffer. The sample in the single peak was concentrated, flash-frozen and stored at −80 °C before use.

Cryo-EM sample preparation and data collection

Purified TSWV L protein was concentrated to 1.2 mg ml−1 for sample preparation. L protein was mixed with (1) 5′ vRNA(1–17 nts) 5′-pAGAGCAAUCAGGUACAA-3′ OH) and 3′ vRNA(1–17 nts) (5′-pUGUUACCUGAUUGCUCU-3′ OH), (2) 5′ vRNA(1–10 nts) (5′-pAGAGCAAUCA-3′OH), (3) ribavirin or (4) ribavirin 5′-triphosphate at a 1:1.5 molar ratio and incubated for 5 min on ice. Holy-carbon gold grid (Cflat R1.2/1.3 mesh 300) was freshly glow-discharged with a Solarus 950 plasma cleaner (Gatan) for 30 s. A 3-μl aliquot of the complex was transferred onto the grids, blotted with filter paper at 22 °C and 100% humidity, and plunged into ethane using a Vitrobot Mark IV (FEI). For the apo L and L–5′-vRNA(1-17 nts)–3′-vRNA(1-17 nts) complex, micrographs were collected at 300 kV using a Titan Krios microscope (Thermo Fisher), equipped with a K2 detector (Gatan), using SerialEM automated data collection software43. Videos (32 frames, every 0.2 s, total dose 60 e− Å−2) were recorded at a final pixel size of 1.04 Å with a defocus of between −1.5 and −2.0 μm. For the L–5′-vRNA(1-10 nts) and L–ribavirin complexes, micrographs were collected at 300 kV using a Titan Krios microscope (Thermo Fisher), equipped with a K3 detector (Gatan), using SerialEM automated data collection software43. Videos (32 frames, every 0.2 s, total dose 60 e− Å−2) were recorded at a final pixel size of 1.07 Å with a defocus of between −1.5 and −2.0 μm. For the L–ribavirin triphosphate complex, micrographs were collected at 300 kV using a Titan Krios microscope (Thermo Fisher), equipped with a Falcon4 detector (Gatan), using SerialEM automated data collection software43. Videos (32 frames, every 0.2 s, total dose 60 e− Å−2) were recorded at a final pixel size of 1.07 Å with a defocus of between −1.5 and −2.0 μm.

Image processing

For the apo L, L–5′-vRNA(1–17 nts)–3′-vRNA(1–17 nts) complex, L–5′-vRNA(1–10 nts), L–ribavirin complex and L–ribavirin 5′-triphosphate complex, a total of 4,135, 4,557, 4,109, 3,659 and 3,853 micrographs were recorded, respectively. An almost equivalent image processing strategy was used for all the datasets. First, the raw data were processed using MotionCor2 by aligning them and averaging them into motion-corrected summed images. After that, the defocus value for each micrograph was determined using Gctf44. Next, micrographs of high quality were selected for further processing. Particles were then picked and extracted for two-dimensional alignment45. Well-defined particles were selected from here for initial model reconstruction. The initial model was used as a reference for 3D classification. After refinement and post-processing, the overall resolution of the TSWV L apo state, L–5′3′-vRNA, L–5′-vRNA, L–ribavirin and L–ribavirin triphosphate was up to 3.9 Å, 3.7 Å, 3.2 Å, 3.4 Å and 3.0 Å, respectively, on the basis of the gold-standard Fourier shell correlation (threshold, 0.143)46. The quality of the local resolution was evaluated using ResMap47. In L–5′-vRNA datasets, the CTD of L was observed. Among the five 3D classes, the subset, making up 4% of total particles, displayed clear features of the CTD in L. After refinement and post-processing, the overall resolution of TSWV L containing the CTD improved to 4.0 Å. In addition, local refinement and DeepEMhancer were used for further improving the quality of the CTD map48.

Model building and refinement

The initial structure model was built using the AutoBuild program in Phenix and further corrected manually by real-space refinement in COOT v.1.1.11. An AlphaFold model comprising the CTD region was used to guide the model building in COOT of parts of the CTD (in areas where the map density was present but its quality was too low for ab initio unambiguous building).

Infectious clones and virus inoculation

The Agrobacterium tumefaciens culture carrying the full-length TSWV infectious clone L(+)opt, M(–)opt and SR(+)eGFP was obtained from a previously reported source49. For virus inoculation, N. benthamiana leaf tissues infected with the virus were ground at a ratio of 1:10 (w/v) in a 0.01 M phosphate buffer solution at pH 7.0. The resulting crude leaf extracts were then used for rub-inoculation onto the leaves of N. benthamiana plants. The inoculated tobacco plants were cultivated in a growth chamber under day/night temperatures of 23 °C/25 °C, with a photoperiod of 16 h of light and 8 h of darkness.

Antiviral drug treatment of plant leaves

To prepare 100 μg ml−1 solutions of remdesivir (SF1193-5MG, Beyotime), ribavirin (R9644-10MG, Sigma-Aldrich), ribavirin 5′-triphosphate (sc-358826, Santa Cruz Biotechnology) and favipiravir (PHR9070-50MG, Sigma-Aldrich), 1 mg of powder of each of these drugs was mixed with 200 μl of DMSO solvent separately. After the solids were completely dissolved, 10 ml of ddH2O (containing 1% Tween-20) was added to formulate the solution. One hour post inoculation of the N. benthamiana leaves with the virus, solutions containing the drug were evenly applied to the surface of the leaves by a brush.

Fluorescence microscopy

N. benthamiana leaves were collected between 24 and 72 h post agro-infiltration, and GFP fluorescence was assessed using an inverted fluorescence microscope (IX-71-F22 FL/DIC; Olympus). GFP was excited at 488 nm, and the emitted fluorescence was collected through a GFP barrier filter (Olympus). Images were captured and analysed with Image-Pro software (Olympus).

Immunoblots

Total protein was extracted from 1.0 g of agro-infiltrated N. benthamiana leaf tissue in 2 ml of extraction buffer containing 10% (v/v) glycerol, 25 mM Tris-HCl (pH 7.5), 1 mM EDTA, 150 mM NaCl, 10 mM dithiothreitol, 2% (w/v) polyvinylpolypyrrolidone, 0.5% (v/v) Triton X-100 and 1× protease inhibitor cocktail. Protein samples were separated on SDS–PAGE gels, transferred to PVDF membranes, blocked with a 5% skim milk solution and then incubated with primary antibodies, such as anti-GFP (1:5,000, Sigma-Aldrich), anti-TSWV L (1:500, lab-made) and anti-TSWV NP (1:10,000, lab-made)) at room temperature for 2 h or overnight at 4 °C. After incubation with a secondary antibody, such as HRP-conjugated goat anti-rabbit (1:10,000, Sigma-Aldrich) or goat anti-mouse (1:10,000, Sigma-Aldrich), for 1 h, the blots were detected using the ECL Substrate Kit (Thermo Scientific).

Total RNA extraction and northern blotting analysis

Total RNA was extracted from N. benthamiana plant leaves using a Total RNA Purification Kit from Tiangen. For northern blot analysis, specific DIG-labelled probes for both sense and antisense eGFP were synthesized using a DIG High Prime RNA labelling kit from Roche. Total RNAs were separated on 1% formaldehyde agarose gels and then transferred to Hybond-N+ membranes from GE Healthcare. The membrane blots were subsequently hybridized with the specific DIG-labelled probe and processed following the manufacturer’s protocol using a DIG-High Prime Detection Starter Kit II, also from Roche.

EMSA

In a standard EMSA, varying amounts of L protein were incubated in 10 μl of binding buffer (45 mM Tris, 45 mM boric acid, pH 8.4) along with 1.5 pmol of FAM-labelled RNA probe (GenScript). The reaction mixtures were incubated at room temperature for 10 min, followed by the addition of 5 μl of loading buffer (0.1% bromophenol blue, 0.1% xylene cyanol and 15% Ficoll 400 in 0.5× TBE buffer). The mixtures were then separated by electrophoresis on a 1.5% agarose gel in 0.5× TBE buffer (45 mM Tris, 45 mM boric acid, 1 mM EDTA, pH 8.0). The gel was subsequently electroblotted onto a nylon membrane (Hybond N+; GE Healthcare). The FAM-labelled RNA was detected using ChemiDoc Imaging Systems from BioRad.

In vitro TSWV RdRp activity assay

In vitro RdRp assays were performed as described in ref. 39, with minor modifications. The assays used 10–15 μg of purified TSWV RNP in a final volume of 25 μl, containing 20 mM HEPES (pH 7.4), 0.5 mM magnesium acetate, 5 mM MnCl2, 2.5 mM DTT, 1 mM of each NTP (including digoxin-labelled UTP, 53119620, Roche), 0.1% NP-40 and 0.8 U μl−1 RNasin. The reaction mixtures were incubated at 30 °C for 1.5 h. Digoxin-labelled RNA products were separated by electrophoresis on a 1.0% agarose gel, followed by downward northern blotting.

Statistics and reproducibility

All biological analysis experiments were performed in triplicates.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Related Articles

Type 2 immunity in allergic diseases

Significant advancements have been made in understanding the cellular and molecular mechanisms of type 2 immunity in allergic diseases such as asthma, allergic rhinitis, chronic rhinosinusitis, eosinophilic esophagitis (EoE), food and drug allergies, and atopic dermatitis (AD). Type 2 immunity has evolved to protect against parasitic diseases and toxins, plays a role in the expulsion of parasites and larvae from inner tissues to the lumen and outside the body, maintains microbe-rich skin and mucosal epithelial barriers and counterbalances the type 1 immune response and its destructive effects. During the development of a type 2 immune response, an innate immune response initiates starting from epithelial cells and innate lymphoid cells (ILCs), including dendritic cells and macrophages, and translates to adaptive T and B-cell immunity, particularly IgE antibody production. Eosinophils, mast cells and basophils have effects on effector functions. Cytokines from ILC2s and CD4+ helper type 2 (Th2) cells, CD8 + T cells, and NK-T cells, along with myeloid cells, including IL-4, IL-5, IL-9, and IL-13, initiate and sustain allergic inflammation via T cell cells, eosinophils, and ILC2s; promote IgE class switching; and open the epithelial barrier. Epithelial cell activation, alarmin release and barrier dysfunction are key in the development of not only allergic diseases but also many other systemic diseases. Recent biologics targeting the pathways and effector functions of IL4/IL13, IL-5, and IgE have shown promising results for almost all ages, although some patients with severe allergic diseases do not respond to these therapies, highlighting the unmet need for a more detailed and personalized approach.

Iron homeostasis and ferroptosis in muscle diseases and disorders: mechanisms and therapeutic prospects

The muscular system plays a critical role in the human body by governing skeletal movement, cardiovascular function, and the activities of digestive organs. Additionally, muscle tissues serve an endocrine function by secreting myogenic cytokines, thereby regulating metabolism throughout the entire body. Maintaining muscle function requires iron homeostasis. Recent studies suggest that disruptions in iron metabolism and ferroptosis, a form of iron-dependent cell death, are essential contributors to the progression of a wide range of muscle diseases and disorders, including sarcopenia, cardiomyopathy, and amyotrophic lateral sclerosis. Thus, a comprehensive overview of the mechanisms regulating iron metabolism and ferroptosis in these conditions is crucial for identifying potential therapeutic targets and developing new strategies for disease treatment and/or prevention. This review aims to summarize recent advances in understanding the molecular mechanisms underlying ferroptosis in the context of muscle injury, as well as associated muscle diseases and disorders. Moreover, we discuss potential targets within the ferroptosis pathway and possible strategies for managing muscle disorders. Finally, we shed new light on current limitations and future prospects for therapeutic interventions targeting ferroptosis.

Climate change threatens crop diversity at low latitudes

Climate change alters the climatic suitability of croplands, likely shifting the spatial distribution and diversity of global food crop production. Analyses of future potential food crop diversity have been limited to a small number of crops. Here we project geographical shifts in the climatic niches of 30 major food crops under 1.5–4 °C global warming and assess their impact on current crop production and potential food crop diversity across global croplands. We found that in low-latitude regions, 10–31% of current production would shift outside the climatic niche even under 2 °C global warming, increasing to 20–48% under 3 °C warming. Concurrently, potential food crop diversity would decline on 52% (+2 °C) and 56% (+3 °C) of global cropland. However, potential diversity would increase in mid to high latitudes, offering opportunities for climate change adaptation. These results highlight substantial latitudinal differences in the adaptation potential and vulnerability of the global food system under global warming.

Targeting of TAMs: can we be more clever than cancer cells?

With increasing incidence and geography, cancer is one of the leading causes of death, reduced quality of life and disability worldwide. Principal progress in the development of new anticancer therapies, in improving the efficiency of immunotherapeutic tools, and in the personification of conventional therapies needs to consider cancer-specific and patient-specific programming of innate immunity. Intratumoral TAMs and their precursors, resident macrophages and monocytes, are principal regulators of tumor progression and therapy resistance. Our review summarizes the accumulated evidence for the subpopulations of TAMs and their increasing number of biomarkers, indicating their predictive value for the clinical parameters of carcinogenesis and therapy resistance, with a focus on solid cancers of non-infectious etiology. We present the state-of-the-art knowledge about the tumor-supporting functions of TAMs at all stages of tumor progression and highlight biomarkers, recently identified by single-cell and spatial analytical methods, that discriminate between tumor-promoting and tumor-inhibiting TAMs, where both subtypes express a combination of prototype M1 and M2 genes. Our review focuses on novel mechanisms involved in the crosstalk among epigenetic, signaling, transcriptional and metabolic pathways in TAMs. Particular attention has been given to the recently identified link between cancer cell metabolism and the epigenetic programming of TAMs by histone lactylation, which can be responsible for the unlimited protumoral programming of TAMs. Finally, we explain how TAMs interfere with currently used anticancer therapeutics and summarize the most advanced data from clinical trials, which we divide into four categories: inhibition of TAM survival and differentiation, inhibition of monocyte/TAM recruitment into tumors, functional reprogramming of TAMs, and genetic enhancement of macrophages.

Endosymbionts modulate virus effects on aphid-plant interactions

Vector-borne pathogens frequently modify traits of their primary hosts and vectors in ways that influence disease transmission. Such effects can themselves be altered by the presence of other microbial symbionts, yet we currently have limited understanding of these interactions. Here we show that effects of pea enation mosaic virus (PEMV) on interactions between host plants and aphid vectors are modulated by the presence of different aphid endosymbionts. In a series of laboratory assays, we found strong interactive effects of virus infection and endosymbionts on aphid metabolomic profiles, population growth, behavior, and virus transmission during aphid feeding. Furthermore, the strongest effects—and those predicted to favor virus transmission—were most apparent in aphid lines harboring particular endosymbionts. These findings show that virus effects on host-vector interactions can be strongly influenced by other microbial symbionts and suggest a potentially important role for such interactions in disease ecology and evolution.

Responses

Your email address will not be published. Required fields are marked *