Introduction

The replication of genomic material is a key stage of the lifecycle of all viruses. Some of the molecular details of replication are well defined as a result of elegant in vitro analyses. It has also been demonstrated that for many viruses, replication occurs in macromolecular complexes involving both viral and cellular components1,2,3,4,5,6. However, the molecular details of the architecture of replication complexes remains unknown and thus we do not fully understand how proteins and viral genetic material are associated during the replication process or how dynamic rearrangements of higher-order assemblies may regulate RNA synthesis during viral replication.

The genomes of positive-sense RNA viruses act directly as a messenger RNA and hence are translated immediately after infection. For picornaviruses such as poliovirus (PV), rhinovirus and foot-and-mouth disease virus (FMDV), proteolytic processing of the viral polyprotein ultimately releases the nonstructural proteins involved in genome replication. A key protein is the RNA-dependent RNA polymerase, RdRp, termed 3Dpol. This enzyme (in complex with other viral proteins e.g., 2C, 3A, and 3CD and a replication primer protein termed 3B) produces negative-strand genomic intermediates, from which new positive strands are synthesised7,8,9,10.

Biochemical and structural studies have yielded a very good understanding of the stages involved in each step of nucleotide incorporation11. In common with other DNA and RNA polymerases, the structure of the picornavirus 3Dpol resembles a cupped right hand, with three defined subdomains, termed thumb, fingers and palm12,13. The finger subdomains are extended in picornaviral 3Dpol structures, when compared to many other polymerases, and contact the tip of the thumb subdomain. This contact is important for maintaining polymerase structure whilst the geometry of the fingers against the thumb form part of the NTP entry channel. The single-stranded template RNA enters 3Dpol through a channel at the face distal to the active site and takes a right-angle turn at the palm active site when the product-template duplex is formed. This emerges from the large front opening14,15,16,17,18,19. There is, however, no strand separation system identified from polymerase biochemical-structural data alone, but such a mechanism must presumably exist to prevent the formation of a stable RNA duplex.

In contrast to its catalytic role in replication, the non-catalytic roles of 3Dpol are less defined. It forms part of a larger precursor 3CD, and is essential for RNA replication via genomic cyclisation9,20,21. Previous studies have also suggested that 3Dpol (or larger precursors thereof) has an essential cis-acting role in replication (i.e., a structural function which needs to be provided from the replicating genome), which is distinct from a separate trans-acting role (i.e., the enzymatic function of the polymerase that can be supplied by a separate genome)19,22. Dimers and higher-order oligomers of RdRp molecules have been reported in a number of positive-sense RNA viruses23,24,25,26,27. Furthermore, the 3Dpol enzyme of both PV and FMDV has been previously shown to form higher-order fibrillar structures25,28,29,30,31. PV 3Dpol has also been demonstrated to form planar lattices and higher-order structures have been visualised within infected cells25. It has therefore been postulated that the higher-order structures act as a replication scaffold; however, no RNA has been shown within the PV 3Dpol structures and furthermore, these were demonstrated to form in vitro in the absence of RNA28,30,31. In contrast, the FMDV fibrillar structures were only detected in vitro in the presence of RNA. Indeed, our previous studies showed that primer, RNA template and nucleotides were required for fibril formation, thus suggesting that they could be part of an active replication complex29.

Using immunofluorescent confocal and electron microscopy, FMDV has been shown to dysregulate Golgi and ER-derived membranes in a similar way to PV1,5,6,32,33, but distinct membrane-bound replication complexes comprised of viral RNA, structural and nonstructural proteins, and host-cell proteins have yet to be identified for FMDV. Although both viruses are members of the Picornaviridae, the Aphthovirus FMDV has several unique features that distinguish it from other members of the family including an extended 5′ untranslated region and triplication of the primer encoding region.

Picornaviruses are responsible for several important human and animal diseases. For example, FMDV, the causative agent of foot-and-mouth disease, is a significant economic problem, such as occurred in the UK in 2001 (where economic losses to farming and tourism were estimated at ~£15 billion at current value)34,35. Therefore, a more complete understanding of the composition of the replication machinery is important for therapeutic development. Here, we have trapped the products of a replication assay by crosslinking and thus captured higher-order structures in a range of conformations. Using cryo-electron microscopy, we have calculated the structures of 3Dpol fibrils in 9 conformations at 7–9 Å resolution, with some of these occurring within the same fibril. We also probed the importance of 3Dpol-3Dpol interfaces within the fibril by mutagenesis and we use these data to build a model of conformational changes in higher-order assemblies of 3Dpol that may be important in replication.

Results

Wild-type 3Dpol can catalyse incorporation of UTP into nascent RNA

We first sought to demonstrate that the recombinant 3Dpol protein was able to bind a synthetic template RNA and was catalytically active. We produced wild-type (WT) 3Dpol in E. coli as previously described15,29, alongside 3Dpol with mutations in the active site (3Dpol-GNN) as a control. To assess the affinity of RNA binding by 3Dpol, a fluorescence polarisation anisotropy (FPA) assay was performed. Samples of 3Dpol-WT and 3Dpol-GNN were titrated with FITC-labelled 13mer poly-A RNA (Fig. 1a), resulting in Kd values of WT and GNN of 5.83 ± 1.3 and 1.92 ± 0.67 μM, respectively (Fig. 1b). These data showed that both active and inactive enzymes were able to bind RNA in vitro, with the mutant demonstrating significantly higher affinity (p < 0.01).

Fig. 1: RNA binding affinity.
figure 1

Fluorescence polarisation anisotropy (FPA) to determine RNA binding affinities of recombinant 3Dpol-WT and 3Dpol-GNN. a Binding of FITC-labelled 13mer poly-A RNA to recombinant 3Dpol calculated as average anisotropy after 30 min incubation. b Mean Kd values comparing the affinity of 3Dpol-WT and 3Dpol-GNN with 13mer FITC-labelled poly-A RNA (n = 3 ± SEM, analysed by interpolation of a standard curve, two-tailed unpaired t-test, **p < 0.01).

The ability of the protein to incorporate [α-32P] UTP into RNA, was assessed as described in29 utilising a poly-A template and an oligo-d(T) primer, supplemented with a 5:1 ratio of UTP to [α-32P] UTP. After incubation at 30 °C, polymerase activity was assessed by dot-blot and quantified by scintillation counting (Fig. 2a, b). Normalised scintillation counts quantifying the dot-blot samples incubated for 30 min (Fig. 2b) show a 25-fold greater incorporation of [α-32P] UTP for WT compared to the enzymatically inactive 3Dpol-GNN.

Fig. 2: Polymerase activity assay.
figure 2

Polymerase activity assays were used to determine activity of recombinant 3Dpol. a Dot-blot showing that incorporation of [α-32P] UTP is only possible when all assay components are included. b Quantification of levels of incorporated radiolabelled nucleotides after 30 min incubation by scintillation counts per minute normalised to WT (n = 3 ± SEM, analysed by two-tailed unpaired t test, ****p < 0.0001). The GNN mutant indicates background levels.

Formation of higher-order complexes

We have previously demonstrated the formation of higher-order fibrillar complexes of 3Dpol by negative stain transmission electron microscopy (TEM)29. To investigate the assembly of the fibrils, we employed glutaraldehyde crosslinking both to trap higher-order oligomers and stabilise the fibrils. Samples of recombinant 3Dpol were incubated at 30 °C in the presence of glutaraldehyde and the RNA template-primer used for the activity assay described above. At time intervals, the reaction was stopped by the addition of Laemmli buffer. The samples were resolved by gradient SDS-PAGE and 3Dpol-WT or 3Dpol-GNN multimers detected by anti-3Dpol immunoblot (Fig. 3a, c, respectively) and quantified by densitometry (Fig. 3b, d, respectively). An increase in oligomer formation over time was evident, which seemed to correlate predominantly to the addition of dimers. Higher-order oligomer formation appeared to be more rapid for 3Dpol-GNN than WT, however, the differences were not statistically significant. Although higher-order oligomers were present in the absence of RNA, this did not correlate with the formation of fibrils. At later timepoints, complexes were too large to migrate into the gel (Fig. 3a, c).

Fig. 3: Elongation of protein monomers into higher-order oligomers.
figure 3

Polymerase activity assay in the presence of glutaraldehyde, showing the formation of higher-order oligomers by representative Western blot of a polymerase activity assay of a 3Dpol-WT and c 3Dpol-GNN, b and d. Quantification of 3Dpol higher-order oligomer formation of 3Dpol -WT and 3Dpol-GNN, respectively, over time by densitometry. Control reactions were included without primer (-P), without RNA template (-R), or without glutaraldehyde (-G) and analysed after 120 min. Values normalised to the average densitometry values (a.u) at each timepoint across three biological repeats using ImageJ. Approximate molecular weights are indicated, by comparison to Coomassie-stained gels. Comparisons between timepoints are within 3Dpol-WT and 3Dpol-GNN respectively (n = 3, ± SEM, analysed by two-way ANOVA, *p < 0.05, **p < 0.001).

Structural analysis of FMDV 3Dpol fibrils

To investigate higher-order fibril structures in the presence of glutaraldehyde, samples of 3Dpol-WT and 3Dpol-GNN were taken at 30 min. Samples were prepared for TEM imaging using negative stain (Fig. 4a, b). This analysis showed that 3Dpol-WT was able to form higher-order fibrillar structures, similar to those we reported previously29 (Fig. 4a). Larger structures produced by 3Dpol-GNN appeared as protein aggregates (Fig. 4b), similar to those occasionally found in samples of 3Dpol-WT incubated in the absence of RNA, suggesting that the faster oligomerisation of 3Dpol-GNN (shown in Fig. 3) may result in less ordered assembly. The fibrils were susceptible to RNase treatment (Fig. S1), in agreement with previous data29. We have previously reported the X-ray structure of a 3Dpol bearing a mutation in the RNA exit channel19. We found this mutant to also be deficient in fibril formation, consistent with a deficiency in RNA binding (Fig. S1). The availability of fixed, stable 3Dpol-WT fibrils led us to investigate their structure by cryo-EM.

Fig. 4: Structure of 3Dpol fibrils.
figure 4

Negative stain TEM of a sample from a polymerase activity assay in the presence of glutaraldehyde after 30-minute incubation. a 3Dpol-WT forms visible fibrils, b 3Dpol-GNN does not show any distinct fibril structure formation. Scale bar 100 nm. c 2D classification and averaging of fibril section images from cryo-EM data shows variation in fibril diameter. dg Helical reconstruction led to the calculation of density maps for nine classes of fibril at 7.3–9.5 Å resolution, four representative maps are illustrated here coloured according to radius (see key below panel g). d The broadest fibril reconstruction; B1 had a diameter of 22.3 nm, C2 symmetry, a helical twist of 39o and an axial rise per subunit of 27 Å. e Conformation B2 has a more open structure, is also C2 symmetric and has a diameter of 21.7 nm, twist of 40o and rise of 31 Å. f N3 is a narrow (Ø = 21.5 nm) and tightly packed fibril which may also be described by a 1-start helical symmetry with a twist of -158.8o and axial rise of 15 Å. g N9 is a narrow, open fibril structure having a diameter of 20 nm, twist of -158.6° and axial rise of 17.4 Å. h Fitting of the X-ray structure for 3Dpol-WT in complex with template primer RNA and ATP (PDB: 2EC0) into each density map shows that these fibrils are two-start helical lattices comprising protofilaments of 3Dpol dimers. i Protofilament one is coloured violet/thistle, while protofilament two is coloured teal/aquamarine, two shades are used in each protofilament to highlight how 3Dpol monomers form ribbons of dimers.

To this end, 3Dpol-WT fibrils were produced as above, prepared for cryo-EM and imaged on a Thermo-Fisher Titan Krios cryo-electron microscope. Fibril sections were picked and subjected to 2D class-averaging, which appeared to show variation in fibril diameters (Fig. 4c). Measurements from 1D projections of these 2D classes along the fibril axis confirmed this, showing a bimodal distribution of diameters that ranged between ~20 nm and ~23 nm. The modal diameters were 21.3 nm and 22.6 nm, indicating the presence of at least two forms of fibril (Fig. S2). These values correlate well with the size of fibrils previously measured from negative stain TEM images29.

Helical reconstruction with 3D classification was performed, separating the dataset into fifteen self-similar classes. In this analysis, nine classes gave reconstructions that exhibited sharp density. These were individually refined leading to the definition of fibril reconstructions with diameters that reflected the bimodal distribution seen in the 2D class averages. Of these nine conformations, two were defined as ‘broad’, B1 and B2—corresponding to the 2D classes having a modal diameter of 22.6 nm. Seven reconstructions were classed as ‘narrow’; N3-N9, corresponding to the larger group of fibrils having a modal diameter of 21.3 nm. The distinction between broad and narrow fibrils was made based on the helical symmetry of each group (see below). Reconstructions ranged in resolution between 7.3 Å and 9.5 Å (Fig. 4d-g, Fig. S3, Table 1, Movie S1).

Table 1 Table to show imaging and reconstruction statistics for each of the fibril reconstructions.

Each reconstruction was readily interpreted by rigid-body fitting of the crystal structure for FMDV 3Dpol (PDB ID 2EC0) revealing that all the solved fibril conformations are right-handed 2-start helical assemblies, each comprising a lattice of repeating 3Dpol dimers (Fig. 4h, i). Dimers are rotationally symmetric about their interface and each of the two protofilaments is a ribbon of 3Dpol dimers that is therefore antiparallel. The nine conformations that were resolved through our analysis show that variation in fibril diameter is a consequence of the two protofilaments sliding and flexing relative to each other (Movie S1). Among these reconstructions, four distinct forms of fibril were identified. The presence of a bimodal distribution of fibril diameters reflects the presence of two helical symmetries. Broad fibrils were found to be 2-start helices having C2 symmetry about the helix axis. Protofilaments had between 9 and 9.2 dimeric subunits per helix turn. Narrow fibrils were not C2 symmetric and exhibited 2-start helical symmetry with between 8.4 and 8.6 3Dpol dimers per helix turn in each protofilament ribbon. A notable difference between the broad and narrow classes of fibrils was that narrow fibrils could also be described by a left-handed 1-start helix of 3Dpol dimers having a twist of ~-159° and an axial rise of 15–17 Å while the twofold symmetric, broad-fibrils could not (Table 1).

Within the narrow and broad classes, fibrils were found to form either tightly packed or open conformations. N3–N8 were tightly packed fibrils in which flexing of protofilaments resulted in a progressive tightening of packing from N3 to N8, with N8 being the most densely packed of the narrow fibril forms. N9 shows a rather different assembly with a longer axial rise and, although the fibril is narrower, a more open structure with large gaps between protofilaments. Similarly, B1 and B2 showed tight and open forms respectively, with B2 being narrower than B1 and having a longer axial rise per subunit. Three-dimensional classification did not distinguish more subtle variations for broad fibrils (similar to those seen in the tightly packed narrow fibril classes (N3-N8)), possibly a consequence of fewer broad fibrils being present in the dataset. In both N9 and B2 (the open fibril conformations) the spaces between protofilaments contained noisy density that was not accounted for by fitting the atomic model for 3Dpol, we postulate that this may be RNA (see below). It is also important to note that the RNA entry and exit sites are not occluded (except by this putative RNA density).

Movie S1 shows morphing between fibril conformations, this is intended to illustrate the differences in these structures rather than model any conformational changes that may occur. To determine whether the fibril classes solved in our analysis represented discrete structures or continuously varying helices, we traced fibril image sections from each of the reconstructions back to their original micrograph (Fig. 5). We found that the narrow fibrils that contributed images to reconstructions N3–N8 varied between these forms along their length, suggesting that these reconstructions represent snapshots of a continuously flexing assembly. Fibrils were seen to switch conformation between tight and open forms. In most cases, changes in fibril conformation were marked by visible constrictions and the different forms presented as discrete regions of longer fibrils, rather than the continuous variation seen for the narrow/tight classes. Overall, these data suggest that narrow and broad forms were distinct from one another. Fibrils that appear to show both broad and narrow forms along their length were, however, noted and in these images pronounced discontinuities/constrictions were observed.

Fig. 5: Filaments annotated to show conformations.
figure 5

Four micrographs are shown annotated to show the position of fibril image sections that contributed to each of the three reconstructions. This shows that fibrils may comprise regions in all three states, however, while narrow/tight-packed forms are seen to alternate along single fibrils, open forms are seen in discrete sections of longer fibrils. In general broad and narrow conformations appear to be distinct, some instances of apparent transition from one narrow to broad are seen, but are accompanied by discontinuities in the density (a-black arrow), changes from tight to open forms are marked by the presence of constrictions (a, d – white arrows).

Inter- and intradimer interactions

As noted above, an inspection of our reconstructed density maps along with rigid-body fitting of high-resolution coordinates for 3Dpol indicated that fibrils each comprise two protofilaments that are ribbons of dimers. Analysis of contacts and clashes between 3Dpol monomers in our models confirms this. These data are also consistent with fibril formation by dimer addition as shown in Fig. 3. In all conformations, the dimer interfaces and dimer-dimer interfaces are similar, while different interfaces between protofilaments were seen in tight and open forms (Figs. 5a and 6a-e; Movie S2; Table S1). To describe the interactions between monomers within fibrils we use a notation incorporating the ‘palm’, ‘fingers’, and thumb’ subdomain nomenclature, while the relative spatial organisation of monomers within the helix are denoted as follows. Where n represents a monomer and then n’ is the twofold symmetry-related monomer within a single dimeric subunit (Fig. S4). Here, n + 1 represents the next monomer, on the opposite protofilament (and n’ + 1 its partner within a dimer). For conformations N3-N9 the transformation from n to n + 1 and all subsequent increments is defined by the left-handed 1-start helical parameters: a rotation of -158° and an axial rise of 15–17 Å. For the C2 symmetric broad fibrils, which do not show 1-start helical symmetry, the transformation from n to n + 1 (and subsequent odd-numbered monomers) is simply a rotation of −180°. The transformation from n + 1 to n + 2 (and all subsequent even-numbered monomers) is a rotation of -140° and an axial rise of ~27–31 Å. Under this notation, for all fibril conformations, n + 2, 4, 6 etc. are successive monomers on the same protofilament as n, while n + 1, 3, 5 etc. are on the other. The dimer interface is between thumb subdomains, predominantly at the C-terminal helix α16 (residues 453-467) of n and n’ (Fig. 6b). Dimers repeat along each protofilament ribbon with the primary interaction being between the thumb/fingertips of n and the palm sub-subdomain of n + 2 (and between the palm subdomain of n’ and thumb/fingertips of n’+2; Fig. 6c). Interactions between protofilaments in each conformation are also twofold rotationally symmetric owing to the antiparallel nature of each protofilament. The interactions are not, however, conserved between fibril forms. In tight conformations, the contact interface is mainly within a loop at amino acid residues 9-17, between n and n’+9 or n’+7 for B1 and N3–N8, respectively (Fig. 6d). In conformation N8, there is an additional point of contact between n and n’+5 between V4 in strand β1 and K65 in the N-terminal loop of the palm subdomain. In open fibrils, the contacts/clashes between protofilaments are between n and n’+7 for B2, and n and n’+5 for N9. In both cases β1 strands in each monomer come together, possibly forming an extended antiparallel β-sheet (Fig. 6f).

Fig. 6: Molecular interactions within the fibril.
figure 6

a The pseudoatomic model for fibril N3 is shown as cartoons coloured according to the subdomain structure of 3Dpol; palm subdomain (plum), fingers subdomain (powder blue) and thumb subdomain (gold). Three interactions between monomers are highlighted: dimer interface (red circle), inter-dimer interface (blue circle) and inter-protofilament interface (yellow circle). Close-up views of the dimer interface b, inter-dimer interface c, and inter-protofilament interface d are presented for N3 with contacts/clashes highlighted by white dashed lines and sidechains of putative contact residues shown. e The pseudoatomic model of N9 is annotated to highlight the inter-protofilament interaction in this conformation (yellow circle), which is shown in close-up view with contacts/clashes highlighted with white dashed lines and putative hydrogen bonds indicated by blue dashed lines f.

Comparison of our structures for fibrils formed of FMDV 3Dpol with previously published oligomeric forms of PV 3Dpol shows that these molecules pack differently. X-ray crystallography studies have shown that PV 3Dpol molecules form functional dimers via a strong interaction termed ‘interface 1’ between the C-terminal helix in the thumb subdomain, and the back face of the palm subdomain23. Furthermore, PV helical polymerase lattices form much larger (50 nm) fibrils, helical reconstructions of which have shown a 6-start helix of protofilaments formed of interface 1 dimers31. However, it should be noted that the PV 3Dpol was untagged, whereas the FMDV 3Dpol used here was C-terminally his tagged. 470 residues (of 476 excluding the tag) were modelled. Although density extends beyond this, it is not possible to determine the position of the tag.

Effect of 3Dpol mutations on replication

The cryo-EM data suggested that the fibril protofilaments are assembled via two main intermolecular interactions, a thumb-to-thumb interaction involving residues Y394, G395, T396, F398, P436, Y455, R456, S457, Y459, L460, V463, N464, V466, C467, D469, and A470, and a fingers/thumb-to-palm interaction (involving residues A22, P23, V25, L73, R76, H322, Y323, E324, G325, D329, T330, Y346, D347, A415, R416, R417, I420, P445, L449, and F450). It should be noted that these contact residues were identified by the fitting of high-resolution data into an intermediate resolution map, thus no account is taken for potential rearrangement of sidechains upon polymerisation. Helical reconstruction of FMDV 3Dpol fibrils revealed considerable conformational flexibility. Many attempts were made to achieve sufficient resolution as to allow ab initio atomic modelling, these were not however successful.

The third interaction site (between protofilaments) is highly variable. Our fitting experiments led to the definition of contact interfaces for tight forms comprising four residues. The considerable flexing of protofilaments argues for a weak interaction at this site. The open forms show a more intriguing potential interaction—formation of an extended antiparallel β-sheet via hydrogen bonding between residues 4–8. Mutagenesis to target this interface would require disruption of secondary structure. For these reasons, we chose to focus our mutagenesis studies on those interactions that support protofilament formation.

To understand whether the intermolecular interactions above are important for viral replication we employed a replicon system36,37. We hypothesised that if the polymerase intermolecular interactions described here are biologically important, disruption of these interactions would impair function, resulting in reduced replicon replication. We therefore selected several residues from both the thumb-to-thumb interface (Y394, F398, R456, Y459, L460, V463, and C467) and the fingers/thumb-to-palm interface (P23, V25, E324, Y346, D347, I420, P445, L449) for mutation. Selecting residues at random from both interfaces prevented biasing the outcome by focusing the mutagenesis on specific secondary structures or amino acid physical properties; however, it should be noted that all of the residues are conserved between FMDV isolates. Each of the selected residues was individually mutated to alanine in the context of a FMDV replicon in which the structural proteins were replaced by GFP. Measuring GFP expression therefore allows the replication competency of each mutation to be assessed. From each mutant replicon, RNA was generated by T7 transcription and transfected into BHK-21 cells. Replicons with 3Dpol-WT and 3Dpol -GNN were included as controls, the latter for the level of translation from the input RNA. RNA replication was monitored by GFP fluorescence, with replication shown as the number of GFP positive cells at 8 h post-transfection, analysed according to our established protocols37 (Fig. 7a).

Fig. 7: Replication of 3Dpol mutants.
figure 7

Intramolecular contact residues were identified and were mutated to alanine within the 3D-encoding region of ptGFP-expressing subgenomic replicon. a Replicon RNA was transfected into BHK-21 cells and replication determined by measuring the levels of GFP expression in comparison to the level of translation (GNN). Residues that significantly reduced the level of GFP expression are highlighted in pink. b and c Location of the amino acid residues that were mutated in this analysis, those that are coloured orange did not significantly affect replicon efficiency, while those that are coloured dark magenta significantly reduced the levels of GFP expression and were located at the modelled 3Dpol dimerisation interface only (PDB: 2EC0). Statistical analysis by unpaired two-tailed t test, ± SEM, **p < 0.001, ***p < 0.0001.

Of the fifteen amino acids mutated only five: F398A, R456A, Y459A, L460A, and C467A, significantly reduced FMDV replicon efficiency compared to the WT replicon. All of these residues were located within the thumb-to-thumb interaction (coloured dark magenta, Fig. 7b). All the other mutations tested showed levels of replication equivalent to a WT replicon (coloured orange Fig. 7b, c). This included two mutations (V463 A and V466A) and all of the mutations made at the fingers/thumb-to-palm interface. Overall, these data suggest that residues involved in dimer formation are important for FMDV replication, potentially through allowing the formation of intermolecular 3Dpol-3Dpol interactions within the FMDV replication complex.

Unassigned density within 3D reconstructions may represent RNA

Within each of our fibril reconstructions, density within the active site of individual 3Dpol molecules was visible, however, this was very weak in comparison to the protein density. Nonetheless, the data suggest the presence of retained RNA at low occupancy within fibrils (Fig. 8, Movie S3). In the open conformations (reconstructions B2 and N9) there is an additional unassigned density that extends between two 3Dpol molecules, bridging an opening in the fibril between protofilaments. The density is poorly resolved but unambiguously present, we suggest that this may represent RNA entering or extruded from the active sites.

Fig. 8: Unassigned density in fibril reconstructions indicates the presence of RNA.
figure 8

Close up views of 3Dpol reconstructions for fibril B1 a and N3 b zoomed into a single monomer and coloured according to the radius (key as in Fig. 4f), shows the presence of weak unassigned density in the active site of 3Dpol that we attribute to RNA (white arrows). d, e the same view is presented as a transparent isosurface with PDB 2EC0 shown. The atomic model is shown with protein represented as a rainbow coloured ribbon diagram and RNA as sticks coloured mauve, to highlight the assignment of protein and RNA to cryo-EM density. c, f Conformation N9 also shows unassigned density in the active site as well as a finger of density linking two protofilaments (red rectangle) that may be RNA template or product extruded from the active site.

Discussion

The RNA-dependent RNA polymerase, 3Dpol, was one of the first proteins to be characterised both in PV38 and in FMDV39 and is the key enzyme in viral replication. Both PV and FMDV 3Dpol have been shown to form higher-order structures in vitro in the form of fibrils25,28,29,30,31,40 that have been postulated to be part of viral replication complexes. It would be expected that replication complexes would be dynamic structures, and although there are structural data available for the PV fibrils these are at relatively low resolution, did not include any RNA and were thus unable to give insight into any dynamics30,31. Here, we report multiple conformations of the fibrils formed by FMDV 3Dpol. Some assemblies exhibited multiple conformations along their length. Furthermore, we report the presence of additional density in fibril reconstructions that we attribute to RNA. We probed the interfaces by mutagenesis and are using these data to propose a model for the potential role of conformational changes in the fibril during genome replication.

We have here demonstrated that the recombinant 3Dpol protein was able to bind RNA and was catalytically active, this led to the formation of fibrils. In contrast, two mutants that were deficient in replication did not form ordered products. We analysed the structures of the assembled fibril products of a polymerase activity assay after crosslinking. This provided a total of 9 different 3Dpol reconstructions, all of which are at higher resolution and distinct from the previously reported structures of PV fibrils, in terms of both dimensions and the orientation of molecules within the helix. The dimerisation interface of the FMDV 3Dpol molecules is novel and distinct from the dimerisation of PV 3Dpol molecules. The formation of fibrils by dimer addition is also supported by the cross-linking data presented here. Mutations that affected the ability of the polymerase to transcribe replicon RNA were introduced within the dimerisation thumb-thumb interface. This contact region has been implicated in the formation of the PV fibrils and shown to be necessary in the production of virus, despite being distant from the active site41. In PV, a model has been described whereby RNA would travel through the active site between polymerase molecules; however, RNA was not detected25,30. In support of this model, density corresponding to RNA was detected in the FMDV fibril structures presented here. This density appears in the active site of all our reconstructions and is also seen linking individual molecules in the open fibril conformations.

The formation of a classical intracellular, membrane-associated replication complex for FMDV has yet to be demonstrated within infected cells. One could speculate, in the context of a cellular replication complex, that the polymerase molecules could coat newly synthesised negative-sense intermediates and prevent the formation of double-stranded RNA, a potent innate immune response trigger. Avoiding double-stranded RNA formation would also allow for a more energetically-favourable replication process. By coating RNA by the addition of dimers, helical fibrils similar to those described here could be generated. The concept of 3Dpol acting in a similar way to the N protein of negative-sense viruses or the coat protein of helical viruses such as tobacco mosaic virus is intriguing.

The structures described here could alternatively represent trapped conformers of the active replication machinery within the functional complex. The different fibril conformations observed could represent different stages of replication. Although there is some density for RNA in the active site of 3Dpol molecules within the fibrils, a model where RNA is looping out from the body of a dynamic fibril would be consistent with the relatively weak RNA density seen overall and with the observation that the fibrils are sensitive to RNaseA. However, it is also tempting to speculate that the density seen between 3Dpol molecules in the open N9 and B2 conformations correlates with newly synthesised RNA being extruded from the active site.

The 9 fibril structures described here are broadly representative of 4 types: narrow, broad, tight and open. Data for the narrow fibrils were consistent with these being a continuous, flexing assembly, with the broad fibrils a separate, distinct class. Given the breadth of fibril morphology shown here, it is therefore possible that both models described above are correct. Some of the fibrils could form by the simple addition and dissociation of 3Dpol subunits and act to coat newly formed RNA transcripts. Others could form part of a dynamic replication complex. Different pools of 3Dpol molecules involved in replication would be consistent with the models we proposed as a result of complementation experiments19. Correlative light-electron microscopy studies in replicon-containing cells are underway in order to further elucidate the formation and dynamics of complexes within cells and to establish if this and other nonstructural proteins, such as 3CD, 3A and 3B, are able to interact with the FMDV 3Dpol fibrils.

Methods

Expression and purification of his-tagged 3Dpol

WT C-terminal His-tagged 3Dpol was expressed from vector pET-28a (a kind gift from Esteban Domingo, Centro De Biologia Molecular Severo Ochoa, Madrid, Spain)42 in Codon plus (+)BL21 using methods described previously15. For lysis, the cells were suspended in 20 mM Tris-HCl [pH 8] lysis buffer with 20 mg/ml lysozyme before sonication as described previously15. Lysates were clarified and the supernatant filtered through 0.22 µm filters and subjected to nickel affinity chromatography using 1 ml HisTrap HP columns (GE Healthcare). Elution was carried out via a 25–500 mM imidazole gradient in 50 mM Tris [pH 8], 500 mM NaCl buffer. Fractions containing the protein were pooled, dialyzed and stored in aliquots at -20 °C as described previously19. A replication-deficient enzyme (termed GNN which includes the active site mutations D338N and D339N19 were expressed and purified in the same way and the mutations confirmed by mass spectrometry.

3Dpol RNA-binding by fluorescence polarisation anisotropy

The ability of FMDV 3Dpol to bind to RNA was examined using fluorescence polarisation anisotropy (FPA) adapted from published protocols43. FPA measures the tumbling rate of a fluorescently labelled molecule. After excitation by polarised light, the excited fluorophore emits light with the degree of polarisation inversely proportional to the tumbling rate of the fluorescent molecule in the solution. Therefore, the smaller the molecule, the higher the tumbling rate resulting in the depolarisation of the emitted light as the fluorophore reorients in solution during the lifetime of the excitation. FPA was undertaken in a black opaque 384-well format (Perkin Elmer). 20 µl of RNA binding buffer (20 mM Tris-HCl (pH 8), 100 mM NaCl, 0.01 % Triton X-100) was added to each well followed by a serial dilution of purified 3Dpol or mutant 3Dpol-GNN. 20 µl of 20 nM 3′FITC-labelled 13mer poly-A (Dharmacon) (kind gift from Thomas Edwards, University of Leeds) was overlaid into each well and samples incubated at room temperature for 1 h. Control samples that contained no fluorescent RNA to establish background polarisation were included. Polarisation was measured using a Tecan plate reader with a 480-nm excitation filter and 530 nm S and P channel emission filters.

3Dpol Activity Assay

Activity of purified FMDV 3Dpol was assessed by the ability of the protein to incorporate [α-32P] UTP into nascent RNA transcripts, as previously described29,44. The final concentration of 3Dpol that was included in the assays was 18.52 µM. Reactions were stopped by the addition of 1.5 µl 500 mM EDTA at the indicated times. For scintillation, 2 µl of each reaction was blotted onto filter paper (Whatmann) and allowed to dry prior to being washed as described previously44. The dry filter paper was exposed to film for imaging before the individual dots on the filter paper were excised and scintillated with the addition of 3 ml scintillation fluid (Perkin Elmer). Radioactive counts were measured as scintillation counts per minute.

For glutaraldehyde crosslinking, 3Dpol activity assays, substituting [α-32P] UTP with 4 mM UTP were employed. Samples were incubated with 5 mM glutaraldehyde at 30 °C for increasing time periods and reactions stopped by the addition of 30 µL 2x Laemmli buffer. The samples were denatured and analysed by gradient SDS-PAGE (12–4% resolving gels, and 4% stacking gel). Immunoblotting was performed as described previously45 with primary rabbit anti-3D 397 polyclonal antibody (kind gift from Francisco Sobrino, Centro De Biologia Molecular Severo Ochoa, Madrid, Spain) and detected with anti-rabbit HRP (Sigma-Aldrich).

Transmission electron microscopy (TEM)

Carbon-coated copper grids were prepared by glow discharge at 25 mA for 30 seconds. The activity assay master mix was prepared as described for the crosslinking assay. Samples were incubated for 30 min at 30 °C and then immediately 10 µl was adsorbed onto carbon-coated copper EM grids for 1 min. The excess sample was removed and the grids washed twice with 10 µl dH2O for 5 seconds prior to staining with 10 µl 2% uranyl acetate for 1 min, the excess was then removed with filter paper. TEM was analysed on the Jeol 1400 (Jeol USA Inc.) at 30,000 x magnification.

Cryo-electron microscopy (Cryo-EM)

Activity assay master mix was prepared as described above and incubated at 30 °C for 30 min. 10 µl of the master mix was loaded on a C-flat (protochips) CF1.2/1.3-4 C grid (1.2 µm holes) that had been prepared by addition of a continuous carbon support film and glow discharged at 25 mA for 1 min. The grids were blotted for 6 seconds at blot force 10, using an FEI Vitrobot Mark IV in chamber conditions of 4°C and 95% relative humidity. These were screened and initial analysis performed using a JEOL 2200 TEM equipped with a Direct Electron DE20 camera, optimised grids were then imaged on a Thermo-Fisher Titan Krios in the Astbury Biostructure Laboratory using a Falcon III detector in integrating mode operated at 40 frames per second. Data were collected at nominal 75,000 x magnification, giving a sampling frequency of 1.065 Å/pixel. Two-second exposures were collected at a dose rate of 54 e2 /second and a defocus range of −1.5 to 2 microns.

Image processing

Preliminary analysis of a small dataset collected on a JEOL 2200 was performed using Relion 2.0.3. A single 2D class average was used to generate a starting model using SPIDER. Initial estimates of approximate helical parameters were derived by taking measurements from 2D class averages (twist = 40° and rise per subunit of 29 Å). Calculation of an initial helical reconstruction was performed giving rise to an interpretable density map with twist = 42.4° and rise of 30 Å. Inspection of this preliminary map revealed a two-start helical assembly (Fig. S5).

Collection of the substantially larger dataset was undertaken at the Astbury BioStructure Laboratory. 5417 micrographs were recorded from which 266,263 images of fibril sections were selected. Motion correction was performed using Motioncor2 with dose-weighting46. CTF estimation was performed using GCTF47. Helical reconstruction was performed using Relion 3.0 followed by 3.148,49. Following 2D classification the dataset was reduced to 247,623 fibril section images that had contributed to sharp class averages. These 2D class average images were carefully realigned to ensure that in all images the helix was parallel to the y axis, using SPIDER50, they were then projected to one dimension (1D) using Bsoft51. Radial density for each class average was plotted using Prism 8 (GraphPad), and fibril diameters were estimated as the distance between the minima on each side of the helix axis. Based on this 2D classification analysis, a dataset of fibril images that gave sharp class averages was selected. 3D classification was performed with k = 15, leading to the definition of nine classes of fibrils having diameters that reflected the full distribution of diameters identified by measurements of the radial density profiles. Inspection of the resulting reconstructions revealed that maps corresponding to the narrower class of fibril (Ø = 20–21.5 nm) could be described as left-handed 1-start helices with twist values of ~-158.8° and axial rise values of ~15–17 Å. The broader class of fibrils (Ø = 21.7–22.3 nm) were twofold symmetric about the helix axis and were processed as right-handed 2-start helices (twist = ~40°, axial rise = 27–31 Å) with C2 symmetry imposed. 3D automated refinement led to the calculation of density maps for each class at resolutions ranging between 7.6 and 10.9 Å. As the reconstructions comprised an asymmetric unit that was a dimer, there was a twofold symmetry present that could not be accommodated within the helical reconstruction process. This was subsequently imposed on half-maps using UCSF Chimera52. Briefly, two copies of each half-map were opened, one was rotated 180° about the x axis (perpendicular to the helix axis), the fit in map function was then used to align the two maps and the density was averaged. Symmetrized half-maps were then imported into Relion for postprocessing analysis, resolution estimates ranged between 7.3 and 9.5 Å.

Fibril images contributing to each reconstruction were traced back to their original micrographs using rln-parts-trace.py (https://github.com/attamatti/relion_particle_backtrace).

Reconstructed density maps were interpreted by docking of published coordinates for the FMDV 3Dpol protein using UCSF Chimera52 and visualised using UCSF ChimeraX53.

Replicon plasmids

The FMDV replicon plasmid, pRep-ptGFP, along with equivalent constructs containing the replication-defective GNN point mutation to the 3Dpol polymerase, have been described previously37. In this replicon, the capsid protein-coding region has been replaced with a green fluorescent reporter protein. Mutations were introduced into the pRep-ptGFP replicon 3Dpol coding region by overlapping PCR site-directed mutagenesis. The sequences of primers are available on request.

FMDV replicon replication assays

To generate in vitro transcribed RNA for transfection, replicon plasmid DNA linearized with AscI (NEB, Ipswich, Massachusetts, USA) was used to setup T7 in vitro transcription reactions, as previously described37. T7 RNA transcripts were purified using a Zymogen RNA Clean and Concentrator-25 (Zymo Research, Irvine, California USA) according to manufacturer’s instructions and the RNA quality was determined by MOPS-formaldehyde agarose gel electrophoresis.

To assay replicon replication, BHK-21 cells seeded into 24-well tissue cultures plates were allowed to adhere for 16 h before transfection in duplicate with in vitro transcribed replicon transcripts using Lipofectin reagent as previously described37. Fluorescent reporter protein expression was monitored using an IncuCyte Zoom Dual Colour FLR (Essen BioSciences), within a 37 °C humidified incubator, scanning hourly up to 24 h post-transfection. Images were captured and analysed using the associated software for fluorescent protein expression. Control transfections (untransfected and the 3Dpol-GNN transfection to assess the level of input translation) were used to determine fluorescence thresholds and identify positive objects from background fluorescence. The number of positive cells per well was determined from the average of up to nine non-overlapping images per well.

Statistics and Reproducibility

The methods for statistical analysis and sizes of the samples (defined as n) are specified in both the results section and individual figure legends for all of the quantitative data. All values are presented as mean ± SEM with the indicated sample size where applicable. Statistical differences were determined using two-tailed unpaired t test and two-way analysis of variance (ANOVA) where appropriate. Statistics were performed using GraphPad Prism 8.

Reporting Summary

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