Purification, crystallisation, data collection and structure determination of FluPol_B with bound CTD peptide

Influenza B/Memphis/13/2003 polymerase, wild type or with the PA K135A mutation to eliminate endonuclease activity, was expressed and purified as described previously [22].

For crystals enabling high resolution visualisation of CTD binding in site 2B, FluPol_B PA mutant K135A at 9 mg ml⁻¹ (35 μM) was mixed with 40 μM of nucleotides 1–13 vRNA 5’ end (5′-pAGUAGUAACAAGA-3′) and 1.8 mM 28-mer CTD peptide (YSPTpSPS)₄ in a buffer containing 50 mM HEPES pH 7.5, 500 mM NaCl, 5% glycerol, 2 mM TCEP. Hanging drops for crystallisation were set up at 20 °C. Rod-shaped crystals growing up to 700 μm in length appeared one week after set-up in mother liquor containing 100 mM tri-sodium citrate and 13% PEG 3350 with a drop ratio of 0.5 μl + 2 μl protein to well solution. Crystals were cryo-protected with additional 20% glycerol and 1.8 mM CTD peptide in mother liquor and flash-frozen in liquid nitrogen. Data were collected on ESRF beamline ID29 and integrated with an ellipsoidal mask using AUTOPROC/STARANISO to an anisotropic resolution of 2.42–2.95 Å. The structure was solved using molecular replacement with PHASER [34] using PDB:5FMZ as model [35]. The model was iteratively corrected and refined using COOT [36] and REFMAC5 [37] and quality-controlled using MOLPROBITY [38]. See Table 1 for data collection and refinement statistics.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 1. Crystallographic data collection and refinement statistics.

https://doi.org/10.1371/journal.ppat.1010328.t001

For crystals enabling simultaneous visualisation of CTD binding in sites 2A and 2B, wild type FluPol_B at 11.7 mg ml⁻¹ (45 μM) was mixed with 52 μM of nucleotides 1–14 (5′-pAGUAGUAACAAGAG-3′) and 1–18 (5′-UAUACCUCUGCUUCUGCU-3′) of respectively the vRNA 5’ and 3’ ends, 104 μM of capped 13-mer (5’-m⁷GpppAAUCUAUAAUAGC-3’), and 380 μM CTD 28-mer Ser5P peptide, in buffer containing 50 mM HEPES pH 7.5, 500 mM NaCl, 5% glycerol, 2 mM TCEP. Hanging drops for crystallisation were set up at 4° C. Diamond-shaped crystals growing up to 200 μm in size appeared in two to three weeks in mother liquor containing 200 mM di-ammonium phosphate, and 100 mM sodium acetate at pH 4.4, with a drop ratio of 1 μl + 2 μl protein to well solution. The drops were soaked with 840 μM CTD peptide for 17 days. Crystals were cryo-protected with an additional 30% glycerol and 885 μM peptide in mother liquor and flash-frozen in liquid nitrogen. Data were collected on ESRF beamline ID30A1 (MASSIF) and processed and refined as described above, using PDB:5MSG as model for molecular replacement. See Table 1.

Structure determination of FluPol_A(H7N9) core with bound CTD peptide

The core of influenza A/Zhejiang/DTID-ZJU01/2013(H7N9) polymerase comprising PA 201–716, PB1 full-length, PB2 1–127 was expressed and purified from insect cells as described previously [18]. A/H7N9 polymerase core at a concentration of 9 mg/ml was co-crystallised with 60 μm of a 12-mer of the vRNA 5’ end (5’-pAGUAGUAACAAG) in sitting drops at 4° C in conditions of 0.1 M Tris pH 7.0, 13% PEG 8K, 0.2 M MgCl₂, 0.1 M guanidine hydrochloride with drop mixing ratios of 1:2 (protein:well). Crystals grew typically within 4–5 days and diffracted to around 3.5 Å resolution. A four-repeat pS5 CTD mimicking peptide (Tyr-Ser-Pro-Thr-pSer-Pro-Ser)₄ was soaked into existing crystals at a concentration of ~2 mM over a period of 24 h. Data were collected on ESRF beamline ID29 and processed and refined as described above, using previously described apo-H7N9 core structure ([18, 39], PDB:6TU5) as model for molecular replacement. See Table 1 for data collection and refinement statistics. The PDB numbers for the new protein structures are provided in Table 1: 7Z42, 7Z43 and 7Z4O.

Cells

HEK-293T cells (purchased at ATCC, CRL-3216) were grown in complete Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with 10% fetal bovine serum (FBS) and 100 U/mL penicillin and 100 μg/mL streptomycin. MDCK cells (provided by the National Influenza Center, Paris, France) were grown in Modified Eagle’s medium (MEM, Gibco) supplemented with 5% FBS, 100 U/mL penicillin and 100 μg/mL streptomycin.

Plasmids

The reverse genetics plasmids derived from the IAV A/WSN/1933 (WSN) [40] and the IBV B/Brisbane/60/2008 [41] were kindly provided by Pr. G. Brownlee (Sir William Dunn School of Pathology, Oxford, UK) and Pr. D. Perez (College of Veterinary Medicine, University of Georgia), respectively. For polymerase activity assays, a pPolI-Firefly plasmid encoding the Firefly luciferase sequence in negative polarity flanked by the 5’ and 3’ non-coding regions of either the IAV or IBV NS segment was used. The pRenilla-TK plasmid (Promega) was used as an internal control. The WSN pcDNA3.1-PB2, -PB1, -PA plasmids [32] and B/Memphis/13/2003 (Memphis) pcDNA3.1-PB2, -PB1, -PA and -NP plasmids [42] were described previously. The WSN-NP Open Reading Frame (ORF) was subcloned into the pCI-plasmid. The WSN and Memphis pCI-PB2-G1 and pCI-PA-G1 plasmids used for Gaussia luciferase complementation assays were constructed as described previously [43]. The RPB1 ORF was obtained from the Human ORFeome resource (hORFeome v3.1), fused to G2 at the N-terminus by PCR and subcloned into pcDNA3.1 (G2-RPB1). The CTD repeats 4 to 51 were deleted from the G2-RPB1 construct by PCR (G2-RPB1ΔCTD). The pCI-RPB2-G2 construct was kindly provided by Dr. B. Delmas (INRAE, Jouy-en-Josas). The wild-type full-length CTD sequences was fused at the C-terminus to an SV40 NLS by PCR using the G2-RPB1 plasmid as a template. The resulting amplicon was cloned in frame downstream the G2 sequence into the pCI vector (G2-CTD). A sequence in which each CTD serine 5 residue was replaced by an alanine was ordered as synthetic gene (GenScript) and subcloned in place of the wild-type CTD sequence into the G2-CTD construct (G2-CTD-S5A). The pCI-G2-NUP62 plasmid was described previously [44]. Mutations were introduced by an adapted QuickChange site-directed mutagenesis (Agilent Technologies) protocol [45]. Primers and plasmid sequences are available upon request.

Protein complementation and minigenome assays

HEK-293T cells were seeded in 96-well white opaque plates (Greiner Bio-One) the day before transfection. For the split-luciferase complementation assays, cells were co-transfected in technical triplicates with 25 ng plasmid encoding the polymerase subunits PB2, PB1 and PA (either PB2-G1 or PA-G1, respectively) and 100 ng of the G2-tagged targets (CTD, RPB1 or RPB2, respectively) using polyethyleneimine (PEI-max, #24765–1 Polysciences Inc). When indicated, the CDK7 inhibitor BS-181-HC (Tocris Bioscience) was added 24 hours post-transfection (hpt) at a final concentration of 20 μM for 1 h. DMSO 0.2% was used as a control. Cells were lysed 20–24 hpt in Renilla lysis buffer (Promega) for 45 min at room temperature under steady shaking (650 rpm) and the Gaussia princeps luciferase enzymatic activity was measured on a Centro XS LB960 microplate luminometer (Berthold Technologies, reading time 10 s after injection of 50 μl Renilla luciferase reagent (Promega)). The Normalized Lumines- cence Ratios (NLRs) were calculated as follows: the luminescence activity (Relative Light Units or RLU) measured in cells co-transfected with the plasmids encoding the PB2/PA-G1 and G2-CTD/RPB1/RPB2 fusion proteins was divided by the sum of the luminescence activities (RLU) measured in control samples co-transfected with either the G2 and PB2/PA-G1 plasmids or the G1 and G2-CTD/RPB1/RPB2 plasmids. For the minigenome assays, cells were co-transfected in technical triplicates with 25 ng of each pcDNA3.1 PB2, PB1, PA, in conjunction with 50, 10 and 5 ng of the pCI-NP, pPolI-Firefly and pTK-Renilla plasmids, respectively. Luciferase activities were measured 20–24 hpt using the the Dual-Glo Luciferase Assay system (Promega) according to the manufacturer’s instructions.

Antibodies and immunoblots

Total cell lysates were prepared in RIPA cell lysis buffer as described before [46]. Immunoblot membranes were incubated with primary antibodies directed against CTD-pS5 (Active Motif, 3EB), CTD-pS2 (Active Motif, 3E10), Gaussia princeps luciferase (New England Biolabs, #E8023) or Tubulin (Sigma-Aldrich, B-5-1-2), and subsequently with the according HRP-tagged secondary antibodies (Jackson Immunoresearch). Membranes were revealed with the ECL2 substrate according to the manufacturer’s instructions (Pierce) and chemiluminescence signals were acquired using the ChemiDoc imaging system (Bio-Rad) and analysed with ImageLab (Bio-Rad).

Production and characterisation of recombinant viruses

The recombinant viruses were produced by transfection of a co-culture of HEK-293T and MDCK cells as described previously [40,41]. The reverse genetics supernatants were titrated on MDCK cells in a standard plaque assay as described before [47]. Plaque diameters were measured upon staining with crystal violet using Fiji [48].

In vitro endonuclease and transcription activity assays

RNA for the activity assays was produced in vitro with T7 polymerase. Recombinant polymerases used corresponding to A/little yellow-shouldered bat/Guatemala/060/2010 and B/Memphis/13/2003 were purified as previously described [42]. 23 nt RNA (5’-GAAUCUAUACAUAAAGACCAGGC-3’) was capped with vaccinia capping enzyme and 2’-O-methyltransferase (NEB) and radiolabelled with [α ³²P]-GTP. For the endonuclease assay 25 nM FluPol_A and 50 nM FluPol_B were incubated with 1.2-fold molar excess of v5’ RNA (FluPol_A: 5’-pAGUAGAAACAAGGC-3’, FluPol_B: 5’-pAGUAGUAACAAGAG-3’) and the capped RNA in reaction buffer containing 50 mM HEPES pH 7.5, 150 mM NaCl, 5 mM MgCl₂, 2 mM tris(2-carboxyethyl)phosphine (TCEP). Transcription reactions were performed with 50 nM FluPol_A or FluPol_B in the reaction buffer, supplemented with 150 nM v3’ template RNA (FluPol_A: 5’-AGUUUGCCUGCUUCUGCU-3’, FluPol_B: 5’-UAUACCUCUGCUUCUGCU-3’) and 250 μM NTP mix (ThermoFisher). 50 μM CTD peptides were added at concentrations corresponding to at least a 10-fold excess over the K_D of the lowest measured affinity for a two-repeat peptide. Two- and four-repeat phosphoserine 5 (pS5) CTD peptides were purchased from Covalab and six-repeat pS5 CTD peptide was synthesised at the Chemical Biology Core Facility at EMBL Heidelberg. Reactions were incubated at 30 °C for 30 min and quenched with RNA loading dye (formamide, 8 M urea, 0.1% SDS, 0.01% bromophenol blue, 0.01% xylene cyanol), supplemented with 50 mM EDTA and boiled at 95 °C. The reaction products were separated on 20% denaturing acrylamide gel (containing 8 M urea) in Tris-Borate-EDTA (TBE) buffer, exposed on a Storage Phosphor screen and recorded with a Typhoon reader. DECADE marker was used as ladder.

CTD sequences and alignment

The RPB1 CTD domain amino acid sequences of Homo sapiens (NP_000928.1), Sus scrofa (XP_020923484.1), Equus caballus (XP_014584045.2), Canis lupus (XP_038521325.1), and Mus musculus (NP_001277997.1) host species were obtained from the RefSeq database [49]. As the predicted RefSeq sequences available for the Gallus gallus (XP_040551262) and Anas platyrhyncos (XM_038172734) RPB1 subunits were only partial, we designed a targeted protein sequence assembly strategy data based on RNA-seq and/or WGS SRA public data available for these two species. To obtain the Gallus gallus RPB1 complete sequence (1969 aa), we first aligned Illumina RNA-seq short reads (ERR2664216) on the human RefSeq curated protein sequence (NP_000928) using DIAMOND algorithm [50], and then used the aligned reads for subsequent Trinity transcript assembly (“–longreads XP_040551262” option to use the partial sequence as a guide) followed by Transdecoder for the ORF prediction [51]. The Anas platyrhyncos RPB1 complete sequence (1970 aa) was obtained by aligning Illumina RNA-seq short reads (SRR10176883) and PACBIO long reads (SRR8718129, SRR8718130) on the JACEUL010000271.1 genomic scaffold by using respectively HISAT2 [52] and minimap2 [53] followed by Stringtie2 [54] with the “–mix” option to allow hybrid de novo gene assembly. The RPB1 coding sequences from Gallus gallus and Anas platyrhyncos are held on the Zenodo repository: https://doi.org/10.5281/zenodo.6467097. The CTD sequences were aligned with SnapGene 6.0 and visualised by Espript 3.0 [55].

Cocrystal structures reveal distinct CTD binding sites in FluPol_B and FluPol_A

Previous structural studies using a four repeat CTD pS5 peptide mimic (YSPTpSPS)₄ [32] revealed two distinct CTD binding sites on FluPol_B, denoted site 1B and site 2B. Site 1B, exclusively on the PA subunit and in which the pS5 phosphate is bound by PA basic residues K631 and R634, is essentially the same as site 1A for bat influenza A polymerase [32] and is thereafter named site 1AB. Site 2B, which extends across the PA-PB2 interface, is unique to FluPol_B and distinct from site 2A for FluPol_A, which is again exclusive to the PA subunit [32]. However, the original crystal structure data for FluPol_B were of insufficient resolution to be able to construct a model for the CTD peptide in site 2B, nor even to define its directionality. To overcome this limitation, we co-crystallised the four repeat pS5 peptide with influenza B/Memphis/13/2003 polymerase in a different P2₁ crystal form, previously used to obtain a structure with the 5’ end of the vRNA [35], and measured anisotropic diffraction data to a resolution of 2.42–2.95 Å (Table 1). The resultant map, which contains two heterotrimers in the asymmetric unit, showed clear electron density in site 2B for both trimers (S1A Fig), into which an unambiguous model for the CTD peptide could be built (Figs 1A and S2A). No significant differences were observed between the two heterotrimers except for small displacements in the poorly defined endonuclease and cap-binding domains, nor in the mode of CTD binding to site 2B, although the density for one of the binding-sites was slightly better due to B-factor differences. Only very weak density for the CTD peptide is observed in site 1AB, perhaps because of competition with a phosphate bound at the position of the phosphoserine. To reconfirm that sites 1B and 2B could be occupied simultaneously, we re-crystallised full promoter-bound FluPol_B with the CTD peptide in the original P3₂21 crystal form, but this time with a capped primer and at lower pH. Under these conditions, the extremity of the vRNA 3’ end is in the RNA synthesis active site [56]. Anisotropic diffraction data to a resolution of 3.12–3.56 Å was measured and the resultant map showed clear electron density for the CTD peptide bound in both sites 1AB and 2B (Fig 1B and Table 1), as reported previously for this crystal form [32] but with slightly improved resolution. Unexpectedly, the CTD peptides bound in site 1AB and site 2B are orientated such that they cannot be linked by the shortest path, as this would be between both N-termini, which are ~ 17 Å apart, whereas the straight-line distance between the C-ter of site 1AB and N-ter of site 2B is ~ 36 (44) Å. These distances suggest that a minimum of 6, probably 7, heptad repeats would be required to occupy both sites contiguously (Fig 1B, dotted red line). This contrasts with the situation in FluPol_A, where the peptide directionality in sites 1AB and 2A allow them to be linked by the shortest path, implying that four heptad repeats is sufficient to occupy both sites (Fig 1C) [32].

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 1. Structural analysis of CTD binding to influenza B polymerase.

A. Overall view of the crystal structure of influenza B polymerase, with bound vRNA 5’ hook (violet) and CTD peptide mimic (slate blue spheres) in site 2B. Ribbon diagram of the polymerase with PA (green), PB1 (cyan), PB2-N (red), PB2-cap-binding (orange), PB2-midlink (magenta), PB2-627 (deep salmon). B. Overall view of the crystal structure of influenza B polymerase with bound promoter (violet and yellow), capped primer (blue) and CTD peptide mimic (slate blue spheres) bound in sites 1AB and 2B. The polymerase is coloured as in (A). The N and C-termini of the two CTD fragments are marked and the red dotted line shows the shortest connection between them with directionality indicated by the arrow. C. Overall view of the crystal structure of bat influenza A polymerase with bound promoter and CTD peptide mimic (slate blue spheres) bound in sites 1AB and 2A ([32], PDB: 5M3H). The colour code is as in (A). The N and C-termini of the two CTD fragments are marked and the red dotted line shows the shortest connection between them with directionality indicated by the arrow. D. Details of the interaction between key residues of the influenza B polymerase PA subunit (green), PB2-N (red) and PB2-627 (deep-salmon) with the CTD peptide (slate blue sticks) in site 2B. Three CTD repeats denoted a (orange), b (cyan) and c (dark green) are involved in this interaction. Hydrogen bonds are indicated as dotted green lines. E. Sequence alignment of the CTD binding sites in the A/WSN/33 (A0A2Z5U3X0) and B/Memphis/13/2003 (Q5V8X3) polymerase subunits PA and PB2. Protein sequences were obtained from UniProt (https://www.uniprot.org/) and aligned with SnapGene 6.0. Key residues for CTD binding are indicated in bold. Identical, similar and non-similar residues are highlighted in purple, light blue and orange, respectively. Grey boxes indicate residues that form a loop. Residues submitted to mutagenesis in this study are indicated with their numbers above (FluPol_A) and below (FluPol_B) the alignment, respectively.

https://doi.org/10.1371/journal.ppat.1010328.g001

Three repeats (designated repeats a, b and c) of the CTD peptide (i.e. Y1aS2aP3aT4apS5aP6aS7a-Y1bS2bP3bT4bpS5bP6bS7b-Y1cS2cP3cT4cpS5cP6c) are visible in site 2B in both structures, including two well-defined phosphoserines (in bold). The N-terminal part of the CTD peptide (Y1a-S2b) forms a compact structure comprising two successive proline turns stabilised by four intra-peptide hydrogen bonds, with P3a stacked on Y1b and P6a stacked on PA/Y597 (Fig 1D). PB2 R134 partially stacks against the other side of the Y1b sidechain, whose hydroxyl group hydrogen bonds to the main-chain of PB2 I135. The phosphate of pS5a forms a strong salt-bridge with PA R608 as well as hydrogen bonding with S7a. FluPol_B-specific PA R608 is in a four-residue insertion (606-GDRV-609) compared to FluPol_A, with hydrogen bond interactions from PA D607 and N611 positioning the side-chain of PA Y597 under the CTD peptide. This configuration of residues seems specifically designed to accommodate the compactly folded CTD peptide. Interestingly, recently identified FluPol_B-like polymerases from fish and amphibians [57] also possess the four-residue insertion in PA. However, only in the Wuhan spiny-eel influenza virus polymerase, which is remarkably similar to human FluPol_B, are all the functional residues Y597, D607, R608, N611 conserved [6,57]. The rest of the CTD peptide (P3b-T4c) has an extended conformation and lies across the PB2 627-domain (Fig 1D). To create the CTD binding surface requires concerted side chain reorientations of PB2 W553, M572 and W575 (Figs 1D and S2A), allowing P6b to pack on W553 and Y1c on M572 and L561, with its hydroxyl group hydrogen bonding to D571. PB2 K556 forms a salt bridge with pS5b.

Most functional studies on CTD are performed with human or avian influenza A polymerase, whereas CTD binding has only been structurally characterised for bat A/little yellow-shouldered bat/Guatemala/060/2010(H17N10) polymerase [32] and C/Johannesburg/1/1966 [33]. Although sequence alignments and mutational studies strongly suggest that the mode of CTD binding is conserved for all IAV-like polymerase [32], we attempted to confirm this by determining the structure of a CTD mimicking peptide bound to influenza A/Zhejiang/DTID-ZJU01/2013(H7N9) polymerase. Previously, we have reported crystals of the A/H7N9 core (PA 201–716, PB1 full-length, PB2 1–127) in the apo-state, which forms symmetrical dimers as described elsewhere [18,58,59]. Here, we soaked the four-repeat pS5 CTD peptide mimic into co-crystals of H7N9 core with the vRNA 5’ hook. The crystals diffracted to a maximum resolution of 3.41 Å (Table 1) and again contain symmetrical dimers of the polymerase core (S2B Fig). We observed clear electron density, not only for the 5’ hook, but also for the CTD peptide bound in site 2A (S1B and S2C Figs), essentially identically bound as previously seen for bat influenza A/H17N10 polymerase (S2D Fig). However, there was no CTD peptide bound in site 1AB. The most likely explanation for this is that in the symmetrical dimeric form of influenza A (core only, or full trimer), both polymerases are in the so-called ‘dislocated’ conformation [18] with an open active site. In particular, PA regions 425–452 and 586–618 are rotated by ~20°, compared to the active, monomeric promoter bound state (e.g. A/H3N2 polymerase structure, [58], PDB:6RR7). This particularly affects the position of key site 1AB binding site residues Y445, E449 and F612 (S2D Fig), thus preventing CTD at this site, while not affecting binding to site 2A.

Alignments of PB2 and PA amino acid sequences from the representative influenza A and B viruses used in our study are shown in S3 and S4 Figs, respectively. The PB2 and PA sequence are very conserved, with > 94% identity among the A/WSN/1933(H1N1), A/NT/60/1968(H3N2) and A/Zhejiang/DTID-ZJU01/2013(H7N9) viruses (details provided in S3 Fig), and > 99% identity between the B/Memphis/13/2003 and B/Brisbane/60/2008 viruses. Key residues for CTD binding (highlighted in S3 and S4 Figs) are strictly conserved between the two representative type B viruses on one hand, and among the H1N1, H3N2 and H7N9 type A viruses, on the other hand. Most of them are identical or similar in the A/little yellow shouldered bat/Guatemala/060/2010 (H17N10) virus compared to the other type A viruses. When compared across type B and type A virus sequences, key residues for CTD binding appear mostly conserved in site 1AB, non-conserved in site 2A and partially conserved in site 2B (Figs 1E and S5).

The FluPol-CTD interaction can be monitored using a cell-based luciferase complementation assay

To confirm the structural findings of this study and investigate the distinctive features of CTD binding sites in FluPol_B and FluPol_A in the cellular context, we set up a CTD-binding assay using the Gaussia princeps luciferase trans-complementing fragments (G1 and G2) [43]. The full-length CTD was fused to G2 and the SV40 nuclear localization signal (G2-CTD). PB2 or PA were fused to G1 at their C-terminus (PB2-G1 and PA-G1, respectively, schematically represented in red in Fig 2A and 2B) as FluPol was shown to retain activity when tagged at these sites [44]. Upon co-expression of G2-CTD and the three polymerase subunits (including PB2-G1 or PA-G1), a luminescence signal resulting from the FluPol-CTD interaction was measured, which was generally higher for FluPol_B compared to FluPol_A (Fig 2A and 2B). The interaction signal decreased when PB1 was omitted, in agreement with previous reports that the FluPol-CTD interaction depends on FluPol assembly [28], and it was independent of FluPol catalytic activity (Fig 2A and 2B, PB1 D444A-D445A mutant [60]). When key CTD-contacting residues of FluPol_A were mutated (PA K289A and R638A [32]), the interaction signal was significantly decreased compared to PA wt (Fig 2C and 2D). To test whether the FluPol-CTD binding assay reflects the dependency on the phosphorylation of the CTD S5 moiety, all S5 residues of the CTD were mutated to alanine (schematically represented in Fig 2E), which prevented S5 phosphorylation as documented by western blot (Fig 2F and 2G, bottom). Although the wt and S5A CTD showed similar steady-state levels of expression, the binding of the S5A CTD to FluPol_B and FluPol_A was significantly decreased compared to WT CTD. Consistently, pharmacological inhibition of CDK7, which represents the major kinase for CTD S5 phosphorylation [61], specifically reduced FluPol_A/B binding to the CTD but not to the FluPol_A interaction partner NUP62 [44] (S6A and S6B Fig). Overall, these data demonstrate that the FluPol-CTD interaction can be accurately monitored in cells using our split-luciferase complementation assay conditions.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 2. Gaussia princeps luciferase-based FluPol—CTD binding assay.

A-B. G2-CTD was expressed by transient transfection in HEK-293T cells together with PB2, PB1 and PA of FluPol_B (B/Memphis/13/2003, blue bars) or FluPol_A (A/WSN/1933, grey bars). Either PB2 (A, hatched bars) or PA (B, filled bars) were C-terminally tagged in frame with G1. As controls, the wild-type (wt) PB1 was replaced by the catalytic inactive PB1 D444A D445A mutant (i) or was omitted (-). Luciferase activities were measured in cell lysates at 24 hpt. Normalised luciferase ratios (NLRs) were calculated as described in the Materials and Methods section. The data shown are the mean ± SD of at least three independent experiments performed in technical triplicates. **p ≤ 0.002; ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test). C-D. The CTD binding of FluPol_A mutants PA K289A and R638A was investigated. HEK-293T cells were transfected as described in (A) and (B), respectively. Relative Light Units (RLUs) are expressed as percentages relative to the FluPol_A PA wt. The data shown are the mean ± SD of three independent experiments performed in technical triplicates. **p ≤ 0.002; ***p ≤ 0.001 (one-way ANOVA; Dunnett’s multiple comparisons test). The dotted line labelled “Ctrl” indicates the background signal i.e. the sum of the luminescence activities measured in control samples co-transfected with either the FluPol-G1 and G2 plasmids, or the G1 and G2-CTD plasmids. E. Schematic representation of the CTD constructs used in (F) and (G): the wild-type G2-CTD (wt, top) and the G2-CTD in which all serine 5 residues were replaced with an alanine (S5A, bottom). F-G. The interaction of the wt or the S5A mutated CTD to FluPol_B (F) or FluPol_A (G) was investigated by transient transfection in HEK-293T cells as described in (A-B). The data shown are the mean ± SD of four independent experiments performed in technical triplicates. *p ≤ 0.033; ***p ≤ 0.001 (two-way ANOVA; Sidak’s multiple comparisons test). In parallel, cell lysates were analysed by western blot using antibodies specific for the pS5 or pS2 CTD, G.princeps luciferase (Gluc) and tubulin. The slow- and fast-migrating bands detected with the pS2 CTD antibody, which likely correspond to the hyperphosphorylated and hypophosphorylated forms of the CTD, respectively, are indicated by a star and a triangle, respectively. The smeared signal in the G2-S5A-CTD samples likely corresponds to the detection of phosphorylation intermediates.

https://doi.org/10.1371/journal.ppat.1010328.g002

Structure-driven mutagenesis confirms FluPol_B and FluPol_A have distinct CTD binding modes on PA

To systematically assess in vivo FluPol-CTD binding, we mutated key residues forming the CTD binding sites in the FluPol_B and/or FluPol_A co-crystal structures and measured the impact of these mutations on CTD-binding using the split-luciferase complementation assay described above (Fig 2). In parallel, we investigated polymerase activity in a minireplicon assay using a viral-like Firefly luciferase reporter RNA, and we rescued recombinant mutant IBVs and IAVs and measured plaque diameters on reverse genetic supernatants as a read-out for viral growth capacity. The A/WSN/1933 residue nature and numbering is used in the text and figures, except when indicated.

The structure and key residues of the CTD binding site 1AB are conserved between FluPol_B and FluPol_A (Figs 1E and 3A). We mutated the pS5 interacting residues PA K635 and R638 to alanines. The mutations did not affect PA accumulation levels (Fig 3B) but significantly decreased in vivo binding to the full-length CTD for both FluPol_B and FluPol_A (Figs 3C and S7A), which is in line with biochemical data obtained in vitro with CTD-mimicking peptides [32]. Consistently, the corresponding recombinant mutant IBVs and IAVs were attenuated or could not be rescued (Fig 3D), and the mutant FluPols activity was generally reduced (Fig 3E). The CTD-binding and minireplicon data were consistent for the PA-K635A-R638A double mutant. However, the single mutants had a stronger impact on FluPol_B than FluPol_A in the CTD-binding assay whereas the opposite trend was observed in the minigenome assay (Fig 3C and 3E). The PA-R638A mutation, which decreased FluPol_B CTD-binding (Fig 3C) and reduced viral growth (Fig 3D), even increased FluPol_B activity in the minigenome assay (Fig 3E). We and others have previously documented discrepancies between FluPol activity as measured in a minigenome assay and viral growth capacity (e.g. [62,63]). The most likely explanation for our data is that each of the PA-K635A and PA-R638A mutations on its own has no or a mild effect on CTD-binding when FluPol_B is incorporated into a transiently expressed vRNP.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 3. FluPol_B and FluPol_A CTD-binding mode at site 1AB.

A. Superposition of the similar CTD binding in sites 1AB on the PA subunit for influenza B (B/Memphis/13/2003, green) and bat influenza A (A/little yellow-shouldered bat/Guatemala/060/2010(H17N10), light grey) polymerases with the CTD peptide as a thin tube (respectively slate blue and light grey). Key conserved residues are indicated in their respective colours, as well as the FluPol_B-specific insertion (PA 608 loop) that is important for part of site 2B. See sequence alignment in Fig 1E. B. HEK-293T cells were transfected with the indicated FluPol_A (A/WSN/1933) and FluPol_B (B/Memphis/13/2003) site 1AB mutants, which were C-terminally tagged with the G1 fragment. Cells were lysed at 24 hpt and analysed by western blot using antibodies specific for G.princeps luciferase (PA-G1) and tubulin. The residue numbering corresponds to FluPol_A (A/WSN/1933). C. In vivo CTD binding of the indicated mutants of FluPol_A (A/WSN/1933, grey bars) and FluPol_B (B/Memphis/13/2003, blue bars). The G2-tagged CTD was expressed by transient transfection in HEK-293T cells together with PB2, PB1 and PA-G1. RLUs are expressed as percentages relative to wt FluPol_A/B. The data shown are mean ± SD of four independent experiments performed in technical triplicates. The dotted line labelled “Ctrl” indicates the background signal i.e. the sum of the luminescence activities measured in control samples co-transfected with either the FluPol-G1 and G2 plasmids or the G1 and G2-CTD plasmids. ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test). D. Characterisation of recombinant IAV (A/WSN/1933, grey dots) and IBV (B/Brisbane/60/2008, blue dots) viruses. Recombinant viruses with the indicated mutations were generated by reverse genetics as described in the Materials and Methods section. Reverse genetic supernatants were titrated on MDCK cells, stained at 72 hpi and plaque diameters were determined using the Fiji software. Each dot represents the diameter of a viral plaque relative to the mean plaque size of IAV wt or IBV wt recombinant virus. (#) not measurable pinhead-sized plaque diameter; (✞) no viral rescue. E. Polymerase activity of CTD-binding site 1AB mutants. FluPol_A (A/WSN/1933, grey bars) or FluPol_B (B/Memphis/13/2003, blue bars) was reconstituted in HEK-293T cells by transient transfection of PB2, PB1, PA, NP and a model RNA encoding the Firefly luciferase flanked by the 5’ and 3’ non-coding regions of the IAV or IBV NS segments, respectively. As an internal control, a RNA-Polymerase II promotor driven Renilla plasmid was used. Luminescence was measured at 24 hpt as described in the Materials and Methods section. Firefly activity was normalised to Renilla activity and is shown as percentages relative to wt FluPol_A/B. The data shown are the mean ± SD of three independent experiments performed in technical duplicates. The dotted lines indicate the background signals in cells transfected with FluPol_A or FluPol_B plasmids minus the PB2_A or PB2_B plasmid, respectively. ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test).

https://doi.org/10.1371/journal.ppat.1010328.g003

CTD binding site 2A differs substantially between FluPol_B and FluPol_A (Figs 1E and 4A). We introduced mutations at residues PA K289, R454 and S420, which are critical for CTD binding to FluPol_A and are not conserved in FluPol_B, and we deleted the PA 550 loop, which buttresses the CTD in FluPol_A and is considerably shortened in FluPol_B. These modifications did not affect PA accumulation levels (Fig 4B) and specifically decreased in vivo CTD binding of FluPol_A but not FluPol_B (Figs 4C and S7B). Consistently, the PA R454A and S420E mutations in the IBV background (IBV numbering: K450A and K416E) did not impair viral growth (Fig 4D) nor did PA K450A affect FluPol_B polymerase activity (Fig 4E). The PA S420E and PA 550 loop deletion impaired FluPol_B activity (Fig 4E), indicating that they hinder a function of the polymerase besides CTD binding.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 4. FluPol_B and FluPol_A CTD-binding mode at site 2A.

A. Superposition of CTD peptide (slate blue tube) bound at site 2A of the PA subunit of FluPol_A (A/Zhejiang/DTID-ZJU01/2013(H7N9), green) with the equivalent region of FluPol_B (B/Memphis/13/2003, wheat), showing similarities and differences in CTD interacting residues. See sequence alignment in Fig 1E. B-E. Protein expression (B), in vivo CTD binding (C), characterisation of recombinant IAV and IBV viruses (D) and polymerase activity (E) of CTD-binding site 2A mutants. Experiments were performed as described in Fig 3B–3E for FluPol_B and FluPol_A site 1AB mutations. C. ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test). D. (#) not measurable pinhead-sized plaque diameter; (✞) no viral rescue, (n.d.) not determined. E. ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test).

https://doi.org/10.1371/journal.ppat.1010328.g004

The PB2 627 domain is involved in CTD binding for FluPol_B but not FluPol_A

The key CTD binding residues and the 3D structure of site 2B are often conserved between FluPol_B and FluPol_A (Figs 1E and 5A). However, CTD binding at site 2B has never been observed in vitro with FluPol_A, and the inserted PA 608 loop (IBV numbering) which buttresses the CTD at the junction between the PA-Cter and PB2 627 domains in FluPol_B is absent in FluPol_A. The PA R608A mutation significantly decreased FluPol_B CTD binding in our cell-based complementation assay (Fig 5B middle panel, no counterpart residue in FluPol_A). Consistently, the corresponding recombinant mutant IBV could not be rescued upon reverse genetics (S8 Fig), and the PA R608 mutant FluPol_B showed reduced polymerase activity (Fig 5B, right panel). We then mutated to alanines the residues W552 and R555 (W553 and R556 according to IBV numbering), which are located on the PB2 627 domain, make contact with the CTD pS5 in the FluPol_B co-crystal and are conserved between IBVs and IAVs (Figs 1E and 5A). The mutations did not affect PB2 accumulation levels (Fig 5C) and either decreased (R555A) or increased (W552A) CTD binding of FluPol_B, whereas they had no effect on FluPol_A CTD binding (Figs 5D and S7C). We speculate that the W552A mutation (W553A according to FluPol_B numbering) could make it easier for the FluPol_B PB2 627 domain to adopt a CTD-bound conformation. To rule out any CTD binding activity on the FluPol_A PB2 627 domain, we deleted the whole domain as described before [64]. The deletion strongly and specifically decreased CTD binding to FluPol_B but not to FluPol_A (Fig 5E). Nevertheless, single amino acid substitutions at residues PB2 W552 and R555 impaired viral growth and polymerase activity of FluPol_B as well as FluPol_A, however with weaker effects on FluPol_A (Figs 1F and 5G). Given the multiple functions attributed to the PB2 627 domain [64], the most likely interpretation of our data is that residues on the PB2 627 domain contribute to the CTD recruitment exclusively for IBVs while they have overlapping CTD-unrelated functions for IBVs and IAVs.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 5. FluPol_A and FluPol_B CTD-binding mode at site 2B.

A. Superposition of CTD peptide (slate blue tube) bound in site 2B of FluPol_B (B/Memphis/13/2003, PA green, PB2-N red, PB2-627 deep salmon) with the equivalent region of FluPol_A (A/NT/60/1968 (H3N2), [58], PDB: 6RR7, PA light green, PB2 pink), showing similarities and differences in CTD interacting residues. See sequence alignment in Fig 1E. B. Protein expression (left), in vivo CTD binding (middle), and polymerase activity (right) of FluPol_B PA R608A. Experiments were performed as described in Fig 3. The data shown are mean ± SD of three independent experiments performed in technical triplicates. ***p ≤ 0.001 (unpaired t test).C-D. Protein expression (C) and in vivo CTD binding (D) of CTD-binding site 2B mutants. Experiments were performed as described in Fig 3B and 3C for FluPol_B and FluPol_A site 1AB mutations. D. **p ≤ 0.002, ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test). E. In vivo CTD binding of PB2 Δ627 domain deletion mutants was investigated. G2-CTD was expressed by transient transfection in HEK-293T cells together with PB2-G1, PB1 and PA of FluPol_A (A/WSN/1933, grey bars) or FluPol_B (B/Memphis/13/2003, blue bars). Luciferase activities were measured in cell lysates at 24 hpt. Normalised luciferase ratios (NLRs) were calculated as described in the Materials and Methods section. ***p ≤ 0.001 (two-way ANOVA; Sidak’s multiple comparisons test). Cell lysates were analysed in parallel by western blot with antibodies specific for the pS5 CTD, G. princeps luciferase (PB2-G1) and tubulin. F-G. Characterisation of recombinant IAV and IBV viruses (F) and polymerase activity (G) of CTD-binding site 2B mutants. Experiments were performed as described in Fig 3D-E for FluPol_B and FluPol_A site 1AB mutations. F. (✞) no viral rescue. G. ***p ≤ 0.001 (two-way ANOVA; Dunnett’s multiple comparisons test).

https://doi.org/10.1371/journal.ppat.1010328.g005

We asked whether this major difference between FluPol_B and FluPol_A CTD binding modes results in different levels of transcriptional activation by CTD mimicking peptides in vitro. A model has been proposed for FluPol_C in which the CTD stabilizes a transcription-competent conformation by binding at the interface of PB1, P3 (PA equivalent), and the flexible PB2 C-terminus [33]. Our observations suggest that the same model could apply to FluPol_B and not to FluPol_A. Therefore, we tested the impact of pS5 CTD mimicking peptides of varying lengths (two, four, or six YSPTSPS repeats) on FluPol_B and FluPol_A in vitro transcriptional activity (Fig 6). The FluPol_B in vitro endonuclease activity (Fig 6A, lane 4) and elongation activity (Fig 6A, lane 8) were increased in the presence of the six-repeat pS5 CTD mimicking peptide compared to the mock control (Fig 6A, lanes 1 and 5), and a similar trend was observed with FluPol_A (Fig 6B, lanes 4 and 8 compared to lanes 1 and 5, respectively). These data complement previous reports that CTD pS5 binding facilitates FluPol_A and FluPol_C transcriptional activity [33] and strengthen the hypothesis that the CTD stabilises FluPol in a transcription-competent conformation [20]. However, our finding that the FluPol_A PB2 627 domain has no CTD binding activity (Fig 5E) indicates that FluPol_A has evolved a divergent mechanism by which the CTD stabilizes the FluPol transcriptase. It also questions whether bridging of the PB2 627 and PA-Cter domains per se is needed for transcriptional activation of FluPols.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 6. Effect of CTD pS5 peptides on in vitro endonuclease and transcription activity of FluPol_B and FluPol_A.

A-B. CTD pS5 peptides of different lengths (two, four and six repeats) were added to (A) FluPol_B (B/Memphis/13/2003) or (B) FluPol_A (bat influenza A (A/little yellow-shouldered bat/Guatemala/060/2010(H17N10)) in vitro activity reactions as described in Materials and Methods. The left four lanes show endonuclease and right four lanes transcription reactions. Lanes 1 and 5: no CTD pS5 peptides. Lanes 2–6, 3–7 and 4–8: addition of two-, four- and six-repeat CTD pS5 peptides, respectively. The different RNA species are labeled. In the FluPol_B samples, a secondary transcription product larger in size than those expected from the size of the template is detected, most likely resulting from stable hybridization of template and primary product as previously described [18]. Quantification of the reaction products of four independent experiments is shown below (FluPol_B in blue and FluPol_A in grey, respectively). The products of the reactions are normalised to the total RNA amount for each reaction and are presented as fractions of the activity of the reaction without peptide. *p ≤ 0.033, ***p ≤ 0.001 (one-way ANOVA; Dunnett’s multiple comparisons test).

https://doi.org/10.1371/journal.ppat.1010328.g006

FluPol_B and FluPol_A bind to the host RNAP II independently of the CTD

The RNAP II transcriptional machinery is highly conserved across eukaryotes, and the CTD in particular shows almost no sequence differences among vertebrate (mammalian and avian) host species susceptible to IAV or IBV infection (S9 Fig). It is therefore unlikely that differences in the CTD amino acid sequence drove the evolution of divergent IBV and IAV CTD-binding modes. We investigated whether FluPol_A/B can interact with the two major RNAP II subunits (RPB1, RPB2) independently of the CTD, using the split-gaussia luciferase complementation assay. The G2-RPB1 and RPB2-G2 fusion proteins were co-expressed with G1-tagged FluPol (PA-G1) by transient transfection as described above. Both combinations resulted in robust and comparable interaction signals (Fig 7A). Interestingly, in the presence of a truncated RPB1 deleted from the CTD (RPB1ΔCTD), a stable interaction signal with FluPol_A/B could still be measured, most likely corresponding to a direct interaction between FluPol_A/B and the core domain of RPB1. RNPA II can form clusters in the nucleus [65] and bridging of the G2-RPB1ΔCTD to G1-tagged FluPol through a wild-type RNAP II complex cannot be formally excluded. However, it is unlikely that an active G. princeps luciferase can be reconstituted in these conditions, given the large size of the complexes involved (RNAPII, > 500 kD; FluPol, > 240 kD). Moreover, as the G2-RPB1ΔCTD protein accumulates at substantially higher levels than the endogenous RPB1 protein (Fig 7A), it is likely to compete with the endogenous RPB1 protein for assembly into RNPAII complexes and their subsequent import into the nucleus, and therefore RNAP II nuclear clusters are likely to contain a majority of G2-tagged RNPA II complexes. Mutations in site 1AB which reduced CTD binding to background levels (Fig 3C) had only weak effects on RPB1 and RPB1ΔCTD binding (PA K631A R634 in Fig 7B, PA K635A R638 in Fig 7C). The same was observed with site 2A mutations (PA K450A in Fig 7B, PA R454A in Fig 7C) and site 2B mutations and (PB2 R555A Fig 7B, PB2 K556A in Fig 7C). These findings, taken together with the relatively low affinity of FluPol for pS5 CTD peptides [32], suggest that the CTD is not the only interface between FluPol and the host RNAP II, and it may not be essential to connect FluPol to the RNAP II but rather to coordinate FluPol cap-snatching.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Fig 7. FluPol_B and FluPol_A binding to RPB1, RPB1ΔCTD and RPB2.

A. Binding of FluPol_A (A/WSN/1933, grey bars) and FluPol_B (B/Memphis/13/2003, blue bars) to RPB1, RPB1ΔCTD and RPB2 was evaluated. RPB1, RPB1ΔCTD and RPB2 were tagged with G2 and expressed by transient transfection in HEK-293T cells together with PB2, PB1 and PA-G1. Normalised luciferase activities (NLRs) were calculated as described in the Materials and Methods section. The data shown are the mean ± SD of five independent experiments performed in technical triplicates. *p ≤ 0.033, **p ≤ 0.002 (two-way ANOVA; Dunnett’s multiple comparisons test). Cell lysates were analysed in parallel by western blot with antibodies specific for the pS5 CTD, G. princeps luciferase and tubulin. Lanes 1 to 3 and 4 to 6 correspond to the co-transfection of FluPol_A and FluPol_B, respectively, with the G2-RPB1, -RPB1ΔCTD and -RPB2 plasmids as indicated. The band detected in lanes 2–3 and 5–6 with the pS5 CTD antibody correspond to the endogenous RPB1 protein. B-C. Binding of (B) FluPol_B (B/Memphis/13/2003) and (C) FluPol_A (A/WSN/1933) mutants in site 2B (PB2 K556/R555A), site 1 (PA K631A R634A/ PA K635A R638A) and site 2A (PA K450A / PA R454A) to RPB1, RPB1ΔCTD and RPB2 was evaluated as described in (A).

https://doi.org/10.1371/journal.ppat.1010328.g007