eLife

eLife. 6: e30483

Cryo-EM structure of the bifunctional secretin complex of Thermus thermophilus

Edoardo D'Imprimahttps://orcid.org/0000-0002-9830-7929¹, Ralf Salzer²†‡, Ramachandra M Bhaskarahttps://orcid.org/0000-0002-7742-0391³†, Ricardo Sánchez⁴, Ilona Rose², Lennart Kirchner², Gerhard Hummer³⁵, Werner Kühlbrandthttp://orcid.org/0000-0002-2013-4810¹, Janet Vonckhttp://orcid.org/0000-0001-5659-8863¹, Beate Averhoff²

1. Department of Structural BiologyMax Planck Institute of BiophysicsFrankfurtGermany

2. Molecular Microbiology and BioenergeticsInstitute of Molecular Biosciences, Goethe University FrankfurtFrankfurtGermany

3. Department of Theoretical BiophysicsMax Planck Institute of BiophysicsFrankfurtGermany

4. Sofja Kovalevskaja GroupMax Planck Institute of BiophysicsFrankfurtGermany

5. Institute of BiophysicsGoethe University FrankfurtFrankfurtGermany

Wesley I Sundquist, Reviewing Editor

‡. Structural Studies Division, Medical Research Council – Laboratory of Molecular Biology, Cambridge, United Kingdom.

†. These authors contributed equally to this work.

© 2017, D'Imprima et al

DOI: 10.7554/eLife.30483

Abstract

Secretins form multimeric channels across the outer membrane of Gram-negative bacteria that mediate the import or export of substrates and/or extrusion of type IV pili. The secretin complex of Thermus thermophilus is an oligomer of the 757-residue PilQ protein, essential for DNA uptake and pilus extrusion. Here, we present the cryo-EM structure of this bifunctional complex at a resolution of ~7 Å using a new reconstruction protocol. Thirteen protomers form a large periplasmic domain of six stacked rings and a secretin domain in the outer membrane. A homology model of the PilQ protein was fitted into the cryo-EM map. A crown-like structure outside the outer membrane capping the secretin was found not to be part of PilQ. Mutations in the secretin domain disrupted the crown and abolished DNA uptake, suggesting a central role of the crown in natural transformation.

Introduction

Natural transformation is a major mode of horizontal gene transfer (Blokesch, 2017; Mell and Redfield, 2014), by which bacteria take up DNA directly from their environment. In this way, bacteria gain novel genetic information, for example metabolic traits, pathogenicity determinants and resistance genes as a driving force for bacterial adaptation and evolution. By enabling pathogenic bacteria to adapt to the human host, natural transformation is clinically relevant.

Natural transformation is a complex process mediated by multi-component transport machineries consisting of so-called competence proteins, which are highly conserved throughout the bacterial world (Averhoff, 2009; Chen et al., 2005; Koomey, 1998). In Gram-negative bacteria this machinery spans the entire periplasm and connects the inner (IM) and outer membranes (OM), mediating DNA binding on the cell surface and subsequent translocation through the periplasm into the cell.

Many conserved proteins of the bacterial DNA uptake machineries are similar to components of protein secretion and type IV pilus biogenesis systems and often play a dual role (Hobbs and Mattick, 1993). Members of the secretin protein family are conserved key components of natural transformation systems in Gram-negative bacteria, and also play important roles in protein secretion, type IV pilus extrusion and the assembly and extrusion of filamentous bacteriophages (Ayers et al., 2010; Korotkov et al., 2011). The OM-embedded C-terminal secretin domain likely provides an aperture for DNA and protein translocation through the OM, connecting the periplasm to the extracellular environment. The members of the secretin family form homo-oligomers consisting of 12 to 19 subunits surrounding a central channel (Bayan et al., 2006; Costa et al., 2015; Jain et al., 2011). The size of the complex, its oligomeric state and the structures of its periplasmic N-terminal domains vary considerably between organisms. Koo et al. (Koo et al., 2016) showed the overall structure of the secretin domain of PilQ from Pseudomonas aeruginosa at intermediate resolution. Moreover, atomic models from cryo-EM structures of type 2 (T2SS) (Yan et al., 2017) and type 3 (T3SS) secretion systems (Worrall et al., 2016) revealed the conserved architecture of the secretin domain as a double β-barrel (‘barrel-in-barrel’ structure) forming the inner gate and the external scaffold. The N-terminal domains of secretins consist of copies of stacked rings which form a periplasmic channel that connects to IM-associated proteins, forming a multi-component secretion system. Structural information on the N domains is scarce, and the N0 domain is usually missing or disordered (Berry et al., 2012; Koo et al., 2016; Worrall et al., 2016; Yan et al., 2017).

In previous studies we discovered a unique secretin (PilQ) complex in the thermophilic bacterium Thermus thermophilus, which is essential for natural transformation and extrusion of type IV pili (Friedrich et al., 2002; Schwarzenlander et al., 2009). We showed that T. thermophilus PilQ contains a thermostable C-terminal secretin domain and an exceptionally long N-terminal tail of six stacked rings (N0-N5) (Burkhardt et al., 2011; 2012; Salzer et al., 2016). Furthermore, we determined the in situ structure of the entire T4PS machinery by cryo-electron tomography (cryo-ET) in the piliated and unpiliated state, revealing conformational changes of the N-terminal ring-forming domains (Gold et al., 2015).

We now report a 7 Å resolution cryo-EM structure of the 13-fold symmetric PilQ complex revealing unique molecular architecture. We used a novel image processing tool, REcenter Particles (REP), and state-of-the-art remote homology detection and modeling methods to determine the structure of the large and flexible PilQ complex. A crown-like density outside the OM was unaccounted for by the molecular model of PilQ and represents an unidentified DNA translocator component essential for the function of the secretin in natural transformation.

Results

Cryo-EM structure of PilQ

We used carbon-back-coated holey carbon films for cryo-EM sample preparation (Koo et al., 2016; Low et al., 2014; Reichow et al., 2010; Schraidt and Marlovits, 2011; Tosi et al., 2014; Yan et al., 2017) to ensure an even distribution of PilQ particles. By this procedure, a large majority of the particle images are side views (Figure 1). Top views were too rare to determine the stoichiometry of the complex reliably by two-dimensional (2D) classification. Top views were subjected to multivariate statistical analysis for an unbiased assessment of particle symmetry (see Materials and methods). This procedure indicated unambiguous 13-fold rotational symmetry of the PilQ complex (Figure 1—figure supplement 1).

Figure 1.View larger versionFigure 1. Electron micrograph of PilQ complex and selected 2D class averages.Typical micrograph of purified WT PilQ recorded at 1.8 μm defocus. Individual PilQ particles are clearly recognizable, either as frequent side-views (blue boxes) or less frequent top-views (green boxes). Insets show a selection of 2D class averages. Class I to IV were obtained by averaging 1096, 935, 728, and 691 particles respectively. Class II (2x magnification) shows the salient features of the PilQ complex with Gate 1 and 2 (blue lines), rings N5 to N0 (red lines) forming the tail (dark blue bracket) and the cap (pink bracket) containing the secretin domain (yellow bracket). Classes III and IV indicate a rocking motion of the cap and N5 ring relative to the long axis of PilQ, with a hinge at the N4 ring. Green arrow, detergent micelle (DDM).

The particle images on the grid and the 2D class averages (Figure 1) show the typical rod shape of T. thermophilus PilQ (Burkhardt et al., 2011; Salzer et al., 2016). The complex has a C-terminal secretin domain, which inserts into the OM, and six N-terminal domains that form six stacked rings N0 to N5, which were found to extend deeply into the periplasm (Burkhardt et al., 2011; Gold et al., 2015). Rings N1 to N5 consist of predicted –βαββα– folds whereas the N0 ring has an unusual –ααβαββα– fold (Burkhardt et al., 2012; Salzer et al., 2016).

Although the secretin domain resists denaturation even at high temperatures in the presence of SDS (Burkhardt et al., 2011), the 2D class averages show clearly that the overall complex is flexible. The distances between cap and tail and between the rings vary and the cap module shows a distinct rocking movement around an apparent hinge at the interface between rings N4 and N5 (Figure 1I,II,III,IV). Consequently, initial 3D reconstructions had a resolution of only ~20 Å (Figure 2F). To improve the map we split the structure into two modules at the N4 ring and reconstructed both parts separately by residual signal subtraction (RSS) (Bai et al., 2015). This resulted in a tail module, comprising rings N0-N4 and a cap module containing the secretin domain along with the N5 ring. 3D classification and auto-refinement improved the resolution of both maps to ~15 Å. Residual signal subtraction improved the map especially in the secretin domain, but the alignment accuracy was poor, due to the low signal of the subtracted particles.

Figure 2.View larger versionFigure 2. Effects of recentering in 2D and 3D reconstruction.(A) Typical 2D class average of 600 PilQ particles. (B–E) Reference-free 2D classification of the same particle subset after residual signal subtraction (RSS) and recentering, which enhances high-resolution features (red arrows) in the cap (B), N4 ring (C), N2N3 ring pair (D), and the N0N1 ring pair (E). (F–H) The improvement upon applying a combination of RSS and recentering is evident after 3D refinement. (F) Standard reconstruction yields a map of ~20 Å resolution. (G) Reconstructing the cap module (pink) and tail module (blue) separately using the RSS procedure improves the resolution to 15 Å. (H) Combination of RSS and recentering after dividing the PilQ complex into four modules doubled the resolution to 7 Å. Images were realigned, classified and refined independently for the cap and N5 ring (1763 particles, orange), N4 ring (3888 particles, yellow), N2N3 ring pair (1585 particles, light green) and N0N1 ring pair (2685 particles, blue). Secondary structure features were resolved in all component maps, revealing the single-polypeptide chain connections between the N0N1 and N2N3 ring pairs. (I) The recentering procedure does not introduce overfitting, as independently aligned and refined modules (N0N1 and N2N3) display consistent features. Features of the N2 ring along with N1-N2 linker region from both modules (blue and green) show maximal overlap. The position of the outer membrane (OM) is indicated by a transparent grey bar.

To improve the map further, we wrote the program REcenter Particles (REP) (Sanchez, 2017) based on the approach of Ilca et al. (2015) and Hou et al., 2017 (see the Appendix for a detailed description of the program). By this procedure the residual signal can be moved to any desired position in the box to maximize alignment precision, and the windows can be cropped to a smaller box size.

We applied the REP procedure to create four new particle data sets, one for each of the rigid parts of the PilQ structure, i.e. the secretin domain with the N5 ring, the N4 ring, and the N2N3 and N0N1 ring pairs (Figure 2B–E,H). The new stacks of re-centered particles were reprocessed using featureless spheres and cylinders as initial references for the secretin-N5 ring and ring modules. This improved the accuracy of the rotational alignment from ~12° to 2°, which doubled the resolution and revealed a wealth of new, interpretable map features (Figure 2 and Figure 2—figure supplement 1). The final component maps have respective resolutions of 7.6 Å, 6.7 Å, 7.6 Å, 7.0 Å for the secretin-N5 ring, N4 ring, and the N2N3 and N0N1 ring pairs (Figure 2—figure supplement 1).

The PilQ model

The higher-resolution component maps provided us with precise information on the position of each domain in the PilQ complex. By remote homology detection (see Materials and methods) we were able to find and align homologous 3D structures for all PilQ domains (Table 1; Figure 3—figure supplement 1). This provided us with starting templates and allowed us to build atomic models of the individual PilQ domains. C₁₃ symmetry was imposed on the ring domains and the symmetrized models were fitted to their map regions, followed by optimization with a flexible fitting procedure (see Materials and methods; Table 1). The linkers connecting the domains were then modeled and refined in a final round of molecular dynamics flexible fitting (MDFF) (Trabuco et al., 2008), resulting in a full-length atomic model of the PilQ complex (Figure 3B–D). The C-terminal secretin domain of PilQ was modeled on the basis of the recently determined high-resolution cryo-EM structures of type 2 (Yan et al., 2017) and type 3 secretion systems (Worrall et al., 2016) and fills a substantial proportion of the secretin-N5 region. The remainder of the PilQ complex is formed by six rings that consist primarily of the N-domains. The rings of domains N1 to N5 show the conserved -βαββα- motif, typical of the N-terminal part of the secretion systems (Video 1). Domains N0 and N1 of T. thermophilus PilQ differ from other known N0 and N1 domain structures (Berry et al., 2012; Spreter et al., 2009) by long insertions at positions where there are loops in the template structures (Figure 3—figure supplement 1). The accuracy of this part is therefore limited by the quality of the sequence alignment for the central scaffold domains N0 and N1 (Figure 3—figure supplement 1). Overall the independently reconstructed ring modules matched the corresponding C₁₃ atomic models reasonably well, as shown by computation of local cross-correlation (See Materials and methods; Figure 3E). The initial best fitted ring models are structurally close to the MDFF refined models (with small RMSD, Table 2).

Figure 3.View larger versionFigure 3. Fit of the PilQ model.(A) Independent reconstructions of the four modules after recentering combined into a single high-resolution map (grey), as indicated by the low-resolution envelope of the entire complex (transparent light blue). (B) Atomic model of C₁₃ symmetric PilQ complex (rainbow colored by chain) fitted into maps of individual modules after MDFF refinement (Trabuco et al., 2008). (C) Cross section of (B), showing the position of gates 1 and 2 relative to the outer membrane (grey bar) and the overall fit of individual domains into each ring. The density of the crown region outside the outer membrane is not accounted for by the PilQ sequence. (D) Top and side views indicating the helical arrangement of PilQ protomers (one protomer red) around the central symmetry axis. (E) Per voxel local cross-correlation values (CC) of C₁₃ symmetric model-derived maps and REP-reconstructed EM maps (outlines) were computed and used to color (red-white-blue) the atomic models (shown as spheres) before MDFF to highlight model/map correspondence.

See full tableTable 1. Modelling and fitting of PilQ domains.Sequences of PilQ domains aligned to their respective templates were used to create atomic models of each domain. Models with the best scores (lower DOPE scores [Sali and Blundell, 1993]) were fitted into the corresponding density maps and screened for best fit orientations (highest cross correlation coefficient).

See full tableTable 2. Statistics of PilQ domain fits.Homology models of the different PilQ domains were fitted into the respective ring module maps. Each domain was initially placed within the density maps in multiple orientations (n = 10,000), spanning the entire Euler angular space, and subjected to steepest descent local optimization. The resulting fits were clustered and ranked by their CC-scores with the corresponding density maps fits. The table lists the top-5 fits for each domain in the Euler-angle search, along with their CC-scores and the number of independent optimization runs ending in this cluster. The top-ranked fits were assessed for compatibility of the N-to-C-terminal domain orientations (bottom-to-top). The compatible orientations with the highest CC-value was selected as initial fit. Subsequent refinement using MDFF flexible fitting (Trabuco et al., 2008) further improved the fits with minimal structural change.

File type mp4 not supported by EPUB. Download File from Web

View larger versionVideo 1. Cryo-EM structure of the PilQ complex.Overall architecture of the PilQ complex. Sagittal sections along the central symmetry axis show the gates; top views of clipped planes indicate the detailed arrangement of all domains and both gates in PilQ from top to bottom.

File type mp4 not supported by EPUB. Download File from Web

View larger versionVideo 2. Atomic model of PilQ protomer.Sagittal and longitudinal sections through the single PilQ domains show the quality of the fit of the homology models of the PilQ protomers. The overall arrangement of one PilQ protomer (chain B, cyan) after flexible fitting and linkers built and the adjacent protomers corresponding to chains A and C (blue and green respectively) show how the they are positioned relative to each other.

Densities corresponding to linker regions between domains N0 and N1, N2 and N3, and N5 and secretin were resolved within the four independent, recentered EM maps (Figure 2I,H). This made it possible to model the linker segments in a straightforward manner after an individual fit of the protomer domains (Figure 3—figure supplement 1). Densities for loops connecting N3 to N4 and N4 to N5 were not resolved. These linker segments were therefore modeled based on the arrangement of rings and their connections as seen in the homologous V. cholerae GspD (Video 2) (see Materials and methods).

Molecular architecture of the PilQ complex

The PilQ complex is 350 Å long and 125 Å wide (Figure 3A). Thirteen PilQ protomers assemble to form a cylinder with a central channel that connects the outer membrane space with the periplasm. The DDM detergent micelle at the top of the secretin domain (Figure 1I–IV) identifies the membrane-embedded region of PilQ (Figure 3—figure supplement 1). The conserved secretin domain has two distinct features: an outer β-barrel cage consisting of four long β-strands per protomer, and two β-hairpins that constitute gate 1 (Figure 3; Figure 3—figure supplement 2), as observed in the other secretin domains (Koo et al., 2016; Worrall et al., 2016; Yan et al., 2017) (Figure 4—figure supplement 1).

Sequence alignment to secretins of known 3D structure indicated that PilQ lacks the lip domain that forms a cap on the extracellular side of the OM in some species (Yan et al., 2017). This feature is present in V. cholerae GspD, while it is absent in other characterized secretion systems like S. enterica InvG, E. coli K12 GspD and K. oxytoca PulD (Worrall et al., 2016) (Figure 4). Surprisingly, the PilQ complex has an additional density outside the OM, which we refer to as the ‘crown’ (Figure 3, Figure 5A, Video 1). The crown consists of 13 distinct spikes, ~40 Å high, which make contact with the inner surface of the outer β-barrel of the transmembrane part of the secretin.

Figure 4.View larger versionFigure 4. Comparison of different secretion systems with T. thermophilus PilQ.Cryo-EM maps of seven different secretion systems aligned against the OM (grey bar) indicate two major classes. (A) E. coli R388 T4SS, S. enterica serovar typhimurium T3SS (consisting of proteins InvG, PrgK and PrgH), E. coli K12 T2SS and P. aeruginosa PilQ do not have any density outside the OM (grey bar). (B) T. thermophilus PilQ, N. meningitidis PilQ and V. cholerae GspD display a conical density outside the OM. Red arrows indicate poorly resolved regions of the N-terminal domains. (C) The structure of the lip domain (-αβ-βα-; cyan) is resolved in V. cholerae GspD (left) and forms the feature outside the OM in (B), but it is absent in PilQ from T. thermophilus (right). (D) Multiple sequence alignment of secretin domains from well-characterized complexes indicates that the lip region (cyan) and the C-terminal helical segments (olive) are absent in all PilQ homologues. We conclude that the features outside the OM (B) are not formed by the PilQ protomer.

Figure 5.View larger versionFigure 5. PilQ mutants K728A and R730A destabilize the crown module.(A) EM map of PilQ complex with fitted model (rainbow-colored from N to C terminus). Residues R730 and K728 (dark blue) are located just below the detergent belt (grey bar), close to the TM region of PilQ in the modeled homo-oligomer. Cross section through the cap region showing the crown module (red) above the OM surface. (B and C) Class averages of negatively stained PilQ-K728A/R730A (B) and wt-PilQ (C) obtained by averaging 374 and 3730 particles, respectively. The crown is absent in PilQ-KR, which shows weak, fuzzy densities instead, indicating that the crown is disordered.

The top of the secretin domain has a ~70 Å aperture, which leads to Gate 1 (Figure 3—figure supplement 2A). Gate 1 is located immediately below the membrane-embedded portion of the β-barrel and defines a ~30 Å aperture. The aperture of Gate 1 is surrounded by multiple lever-like features that emerge from the inner barrel of the secretin domain (Figure 3—figure supplement 2B, Video 1). The N5 ring below Gate 1 has an inner diameter of ~70 Å. The N4 ring is situated below the N5 ring; 2D class averages indicate considerable flexibility in this region (Figure 1I-IV). Ring N4 is separated by ~20 Å from the N2N3 pair. The connection between this ring pair and N4 is not resolved. Given that the linker between domains N3 and N4 is 17 residues long, it is possible that consecutive rings are offset by one or two domains along the circumference of the ring (Figure 3—figure supplement 3) in either direction. The end-to-end distances measured for these alternative topologies suggest a range of 30–60 Å. A 17-residue linker could easily span this distance. Therefore, we modelled this linker segment in two alternative topologies (Figure 3—figure supplement 3).

The N2 and N3 rings form a tight pair (Figure 2D,H). Rings N1 and N0 likewise interact tightly (Figure 2E,H). The N1 ring harbors Gate 2, which has a 30 Å aperture similar to Gate 1 (Figure 3—figure supplement 2B). Gate 2 features a diaphragm-like arrangement with smaller twisted lever elements, which are best seen in top views (Figure 3—figure supplement 2B; Video 1). The extra loops in domain N1 contain several glycine residues (Figure 3—figure supplement 2B, Figure 3—figure supplement 1), which could confer the flexibility that may be required for opening and closing the aperture.

K728 and R730 residues stabilize the crown region of PilQ

The crown region of the PilQ complex outside the outer membrane is not accounted for by the PilQ sequence (Figure 3). Our structural model identifies the region of the PilQ sequence that interacts with the crown (Figure 5A). We exchanged two positively charged residues (K728 and R730) at the interface between the detergent belt and the crown module to alanine by site-directed mutagenesis (see Materials and methods). These PilQ variants were expressed in a pilQ negative T. thermophilus mutant (ΔpilQ::bleo).

EM analysis of the PilQ K728A/R730A (PilQ-KR) variant revealed that it forms PilQ complexes lacking the crown region (Figure 5B–C). To analyze the impact of these residues on PilQ complex stability or assembly, we harvested a T. thermophilus ΔpilQ mutant (ΔpilQ::bleo, negative control), a ΔpilQ mutant complemented with the pilQwt gene (positive control) and a pilQ deletion mutant complemented with pilQ-KR in stationary phase and analyzed the PilQ complexes after boiling the membrane extracts for 30 min in SDS sample buffer (Figure 6A). No PilQ complexes were detected in the ΔpilQ::bleo (negative control) mutant, while the pilQ mutant complemented with pilQwt (positive control) and pilQ-K728A/R730A produced assembled PilQ complexes (Figure 6B). We conclude that K728 and R730 are not essential for PilQ oligomerization or thermostability (Figure 6B), but necessary for crown assembly.

Figure 6.View larger versionFigure 6. Stability of the PilQ-KR variants and twitching motility assay.(A, B) Stability of purified PilQ wildtype and PilQ-KR as indicated by SDS-PAGE (A) and Western blot analysis using PilQ antibodies (1:13000) (B) of membrane extracts boiled for 30 min in SDS sample buffer. (C) Contrast-inverted images of culture plates to quantify the twitching motility of the T. thermophilus strains (Scale bar = 5 mm).

Next, we investigated the overall structural arrangement of PilQ variants. The PilQ complex formed by PilQ-KR was purified from T. thermophilus cells and analyzed by negative-stain EM (see Materials and methods). A total of 2911 particles were manually picked with the boxer module in EMAN (Ludtke et al., 1999) and processed for 2D reference-free classification. This procedure identified 1066 particles (~36%) without crowns (Figure 5B). In these particles the density close to the OM is fuzzy (Figure 5B), while the remaining complex appears unchanged. This observation suggests that the crown module in this PilQ variant has bound but is not properly folded. Interestingly no perceivable differences between the pilQ-K728A/R730A and WT-PilQ complex were observed in the remaining 64% of the particles (Figure 5C). These observations illustrate that residues K728 and R730 are important for PilQ crown assembly or folding.

K728 and R730 play a role in functionality of PilQ in twitching motility and adherence and are essential for natural transformation

To characterize the effect of PilQ-K728A/R730A mutation on the function of the PilQ complex, we quantified the twitching motility of T. thermophilus. The T. thermophilus strain containing WT PilQ-complex has a twitching zone of 2.2 cm after 3 days of incubation. This value was set to 100%. The ΔpilQ::bleo (negative control) was not motile (Figure 6C). The pilQ-KR double mutant showed a twitching zone of 1.4 ± 0.3 cm, which was smaller than that of the WT strain (64% of the wild-type). Next, we measured the adherence of bacterial cells to plastic surfaces, which is mediated by type-IV pili (see Materials and methods). T. thermophilus HB27 was used as positive control and showed a ratio (570 nm/600 nm) of 3.9 ± 1.5 × 10⁻². The negative control (ΔpilQ::bleo) did not adhere at all. The pilQ-KR showed a ratio of 1.9 ± 1.3 × 10⁻². The pilQ-KR mutant therefore exhibited only 49% of the wild-type adherence. These findings suggest that the crown plays a role in type IV pilus-mediated adhesion and twitching motility.

To analyze the role of K728 and R730 in natural transformation, we quantified the natural transformation frequency of the pilQ-mutant complemented with PilQ-KR. The PilQ wildtype cells and the pilQ mutant were used as controls (see Materials and methods). The transformation frequency of the wild-type strain was 2.7 × 10⁻³ ± 4.9 x 10⁻⁴, while ΔpilQ::bleo (negative control) was not transformable. The transformability of the pilQ-KR mutant was reduced to 6% (frequency = 1.6 × 10⁻⁴± 5.4 × 10⁻⁵). This suggests that K728 and R730 play an essential, functional role in DNA uptake and natural transformation.

Discussion

The recentering protocol improves the resolution of flexible particles

Achieving sub-nanometer resolution for this complex required a novel tool for image processing. By applying our new recentering procedure after residual signal subtraction (Bai et al., 2015), we were able to align and refine each rigid module of PilQ complex independently. It also enabled us to visualize the densities of the connecting loops between some of the N-terminal ring domains (Figure 2I).

The overall alignment improved substantially (Figure 2 and Figure 2—figure supplement 1), compensating for the small number of available particles. The improvement is easily explained as follows: A small alignment error at the center of the box becomes more severe with increasing distance from the center of rotation. Moving the particles from the edge of the box therefore results in a significant enhancement in the accuracy of local particle alignment. This feature of our recentering procedure should be especially useful for large, flexible complexes with multiple components such as PilQ and other secretion systems (Worrall et al., 2016; Yan et al., 2017), respiratory supercomplexes (Sousa et al., 2016), spliceosomes (Nguyen et al., 2015), apoptosomes (Pang et al., 2015), or stressosomes (Marles-Wright et al., 2008). In case of our PilQ complex, this procedure improved the map resolution to ~7 Å, whereas the conventional global alignment procedure limits the resolution to ~20 Å (Figure 2F–H).

The ease of implementing the recentering procedure and interfacing it with other cryo-EM software, combined with a substantial reduction in computational cost make it an attractive feature for processing cryo-EM images of large, flexible particles. Reducing the particle box size economizes on computation, particularly when using GPUs (Figure 2—figure supplement 1). Recentering requires only two inputs: the coordinates of the new center and the alignment parameters (see Appendix 1 for details).

Architecture of the PilQ complex

The structure of the T. thermophilus PilQ complex differs from that of other secretion systems (Figure 4). PilQ is evolutionarily related to other secretins and partly resembles T2SS (Figure 4—figure supplement 1). Both complexes have an OM-embedded homo-oligomeric secretin domain with a barrel-within-barrel structure (Figure 4 and Figure 4—figure supplement 1).

Other known secretion systems feature a single gate in the OM (Gate 1) (Costa et al., 2015; Korotkov et al., 2011), which is a conserved feature amongst all secretins. PilQ from T. thermophilus has an additional gate (Gate 2) on the periplasmic side close to the platform apparatus, which connects PilQ to a still not yet identified channel in the inner membrane (Gold et al., 2015). Gate 2 is composed of long loops between the N0 and N1 ring (Figure 3—figure supplement 2). The unusual number of six N-terminal ring domains makes PilQ one of the largest secretion systems to have been characterized so far (Figure 4). The six rings enable the complex to span most of the wide periplasmic space of T. thermophilus (Burkhardt et al., 2011; Burkhardt et al., 2012). Different orientations of the individual domains and inter-protomer interactions account for slight changes in the diameter and height of individual rings. The N1 and N4 domains appear to be slightly tilted towards the inside, compared to the N2N3 pair (Figure 3). This inward tilting reduces the inner diameter of rings N1 and N4. In combination with the N1 inward tilt and extended loops connecting the β1-α1 and β2-β3 of the N-terminal ring, the N1 domain fills additional density within the channel-forming Gate 2 (Video 1, Figure 3—figure supplement 2). The N0 ring is well resolved in the map, whereas in other secretins it was not visible, apparently due to its flexibility (Koo et al., 2016; Worrall et al., 2016; Yan et al., 2017).

Role of the crown module

Multiple sequence alignment of secretin domains (Figure 4D) reveals a cap/lip segment in some species that forms a protrusion outside the OM (Koo et al., 2016; Worrall et al., 2016; Yan et al., 2017). T. thermophilus PilQ does not contain this sequence; however, our EM map shows a large crown module outside the OM (Video 1, Figure 3), which cannot be assigned to any part of the PilQ sequence (Figure 5A) and thus must be another protein component of the complex. The crown module outside the OM is also visible in sub-tomogram averages of the piliated and the non-piliated complex in situ (Gold et al., 2015) (Figure 7). The crown module appears to be difficult to separate from the PilQ secretin domain and resists denaturation with SDS at high temperature (Burkhardt et al., 2011), which indicates a very tight interaction. Negative stain EM showed that the PilQ-K728A/R730A mutation affects the structure of the PilQ complex and interferes with the crown folding or assembly (Figure 5B–C). This mutation has a minor effect on motility and adherence of bacterial cells. However, the transformation efficiency of these cells is severely impaired. Our previous finding that PilQ is essential for DNA binding (Schwarzenlander et al., 2009) leads us to conclude that the crown structure is either directly or indirectly implicated in DNA binding and uptake. We speculate that it may be the extracellular switch that converts the pilus extrusion machinery into a DNA uptake system. Experiments to identify the crown module are currently in progress.

Figure 7.View larger versionFigure 7. Comparison of T. thermophilus T4PS complex in situ to detergent-solubilized PilQ.Cross section of sub-tomogram averages of piliated (light blue) and unpiliated T4PS (light yellow) (Gold et al., 2015) with superposed 7 Å cryo-EM structure of PilQ complex.

PilQ forms the core of the bifunctional T4PS/DNA uptake system and is suggested to interact with a number of periplasmic and outer and inner membrane anchored proteins (Friedrich et al., 2002; Rose et al., 2011; Rumszauer et al., 2006; Salzer et al., 2014; Schwarzenlander et al., 2009; Tsai and Tainer, 2016). Our structure provides a firm base for mechanistic investigations of the complex apparatus that enables natural transformation.

Materials and methods

Expression and purification of PilQ complex and variants

Thermus thermophilus ΔpilQ::bleo mutant cells producing his-tagged wildtype PilQ or his-tagged K728A/R730A PilQ variants encoded from the E.coli/Thermus shuttle vector pDM12 were grown in TM⁺ complex (Burkhardt et al., 2012). Cells were harvested and 40 g cell pellet was resuspended in buffer 1 (40 mM Tris/HCl, 300 mM NaCl, 1.0 M glycine betaine, 5 mM MgCl₂, 5 mM CaCl₂, pH 8.5). The cells were disrupted by ‘French Press’ (3 × 1000 Psi). Membranes were prepared by ultracentrifugation (260,000 g, 4°C, 45 min), and then resuspended in buffer 1. PilQ complexes were solubilized using 3% n-dodecyl-ß-D-maltoside (DDM). After solubilization the supernatant was diluted to a final concentration of 1.5% DDM and incubated with 30 ml of ‘TALON Metal Affinity Resin’ (Clontech Laboratories, Inc., St-Germain-en-Laye, France) for 1 hr at 4°C. The column was washed with buffer 2 (40 mM Tris/HCl, 300 mM NaCl, 5 mM MgCl₂, 5 mM CaCl₂, 0.05% DDM, 12.5 mM imidazole, pH 8.5). Elution was performed with buffer 2 containing 150 mM imidazole. The elution fraction was diluted 1:2 with buffer 3 (40 mM Tris/HCl, 5 mM MgCl₂, 5 mM CaCl₂, 0.05% DDM, pH 8.5) and subjected to anion exchange chromatography (Q-Sepharose), as described recently (Salzer et al., 2016).

E. coli cultures were grown in LB medium using a similar protocol with appropriate kanamycin, ampicillin, and/or bleomycin were added. T. thermophilus strains were cultivated at 68°C in TM⁺ medium to stationary growth phase and cells were harvested. Membrane extracts containing PilQ complex were prepared as described (Salzer et al., 2016). The membrane extracts were boiled for 30 min in SDS sample buffer and separated by 3–12% SDS-PAGE. Proteins were transferred to a nitrocellulose membrane. Western blot analysis was performed with polyclonal PilQ antibodies using a dilution of 1:13000. K728A and R730A variants of PilQ were obtained through site directed mutagenesis. pDM12-pilQ-6his plasmid (Burkhardt et al., 2011) along with primers, SDM_KR-for GATCGGCGAGCTGTTCGCCCAGGCCACGAACGAGAGCACCGACAAGG and SDM_KR-rev AGGGGATGTCCATGAGGAGGGGCACC were used. The plasmids were sequenced to confirm the mutations and transformed into T. thermophilus ΔpilQ::bleo mutant. The resulting pilQ gene (pilQK728A/R730A) was expressed from the plasmid pDM12 under the bc1 promoter (Salzer et al., 2014) in a ΔpilQ::bleo mutant.

Twitching motility, adherence and transformation assay

T. thermophilus strains were grown for 3 days on minimal medium containing 1% BSA. The plates were subsequently stained with Coomassie Blue. T. thermophilus cells were detached from the agar and the twitching zones. After decantation of the cell suspension, colorless regions are visible, where the cells attach. This region (diameter in cm) is defined as the twitching zone. This is made visible as a dark zone by contrast inversion.

Adhesion of bacterial cells, which is mediated by type-IV pili, was quantified using the microtiter adhesion assay (Burkhardt et al., 2012). Cells were inoculated to an optical density of 0.05 (at 600 nm) and incubated for 3 days at 55°C, 64°C, or 72°C. After incubation, the optical density (600 nm) was measured again, followed by a washing step to remove non-sessile cells. Remaining adherent cells were stained with 0.1% crystal violet solution. Excess staining solution was removed by washing three times with water. Crystal violet from the adhering cells was then dissolved in ethanol, and its absorbance was measured at 570 nm. Adherence was quantified as a ratio of absorbance at 570 nm to 600 nm.

Transformation efficiency of WT and mutant cells was analyzed by performing transformations at 68°C on TM⁺ agar medium using 5 μg of genomic DNA of a spontaneous streptomycin-resistant T. thermophilus HB27 mutant (Friedrich et al., 2001). The transformation frequency was calculated as number of transformants per living count. All transformation assays were performed in triplicate.

Negative stain EM of PilQ-KR mutants

We used 2% ammonium molybdate for negative staining of a PilQ-KR mutant sample diluted to a final concentration of 0.01 mg/ml (Salzer et al., 2016). Micrographs were recorded in an FEI Tecnai G2 Spirit operated at 120 kV, at a nominal magnification of 29,000, yielding a pixel size at the specimen level of 4.2 Å. 2D classification was carried out using Xmipp 3.1 using the classification algorithm CL2D (Sorzano et al., 2010).

Cryo-EM data collection

Three µl of wt PilQ sample (concentration = 1.5 mg/ml) was applied to freshly glow-discharged (at 15 mA for 25 s) Quantifoil R1.3/1.2 holey carbon grids (Quantifoil Micro Tools, Germany) back-coated with a thin carbon layer. The grids were vitrified in an FEI Vitrobot Mark IV plunge-freezer at 70% humidity and 10° C after blotting for 8–10 s. Cryo-EM images were collected in a JEOL 3200 FSC electron microscope operating at 300 kV, after coma free alignment, equipped with an in-column energy filter at a slit width of 18 eV. Images were recorded manually at a nominal magnification of 20,000x, yielding a calibrated pixel size of 1.63 Å, on a K2 direct electron detector (Gatan) operating in counting mode. Dose-fractionated 9 s movies of 45 frames were recorded with an electron dose of 0.75 e⁻/Å²/frame at a defocus of 1.2–3.4 μm.

Cryo-EM image processing

A total of 932 micrographs was collected. Whole-image drift correction of each movie was performed using MotionCorr (Li et al., 2013). A second round of whole-image drift correction was performed using Unblur (Grant and Grigorieff, 2015) and the aligned movies were summed and exposure-filtered using the program Summovie (Grant and Grigorieff, 2015). The CTF was determined using CTFFIND4 (Rohou and Grigorieff, 2015) in the RELION workflow (Scheres, 2012). A dataset containing 11,457 particle images manually picked with EMAN boxer (Ludtke et al., 1999) (324 pixels x 324 pixels) was subjected to 2D reference-free classification in RELION (Scheres, 2012) to check the quality of the particle images. To assess the symmetry of the particles, 363 top-views were separately picked and extracted with a smaller box size (168 × 168). Top views were subjected to reference-free classification, but the 2D class averages did not display reliable features because of the low number of particles and the attractiveness effect of the maximum likelihood approach (Sorzano et al., 2010). To assess the symmetry of 2D class averages, the eigenimages were evaluated after multivariate statistical analysis (MSA) (Böttcher et al., 2015; Dube et al., 1993) in IMAGIC (van Heel et al., 1996). From a class average obtained in RELION, a clear top view was selected and 50 randomly rotated copies were produced with the ROTATE-RANDOMLY command, realigned translationally, and subjected to a single round of MSA with the MSA-RUN command. This procedure was repeated with 10 other top views to assess the reproducibility. The eigenimages clearly indicated 13-fold cyclic symmetry (Figure 1—figure supplement 1).

Next a consensus 3D refinement was performed with all particles. As a starting reference a featureless cylinder was generated in Xmipp 3.1 (Sorzano et al., 2010) with the module xmipp_transform_mask with --create_mask --mask cylinder options. After the consensus refinement, the resulting map was used as a starting reference for 3D classification with local angular search. The same map was also used as a template for obtaining a soft-edge shaped mask with the relion_mask_create module in RELION. Several runs of 3D classification were conducted, varying the number of classes in order to assess the consistency of the results. Different symmetries (C₉ to C₁₈) were tested, but only C₁₃ provided consistent reconstructions in terms of interpretable features. Local angular search was performed with 1.875° ± 10° and then 0.9° ± 3° respectively. A dataset of 2646 particle images was obtained after combining the classes displaying the best interpretable features ,which was then subjected to the 3D auto-refinement procedure in RELION with only local angular search using the best alignment parameters obtained at the end of 3D classification. This led to two maps with a resolution of ~20 Å. To account for the observed flexibility of the PilQ complex, the map was divided into two parts: the cap, consisting of the secretin domain and the N5 ring, and the tail containing rings N0-N4. A residual signal subtraction (RSS) procedure (Bai et al., 2015) was applied to both parts, subtracting the signal for the other part from the aligned images, yielding two new datasets of particle images for cap and tail regions (Figure 2,G). Each dataset was then subjected again to a global alignment through a consensus refinement, 3D classification and auto-refinement with only local alignment as described above for the complete complex. Although some new features were partially resolved, the global alignment of the subtracted particles was poor. To resolve this problem, we used our in-house program REP (see the appendix for a detailed explanation) to re-center the residual signal of the particles within a smaller box size. This allowed us to subdivide the PilQ complex into four modules incorporating all the particles of the initial dataset with different box sizes: secretin-N5 domains (160 × 160), N4 ring (120 × 120), N2N3 ring pairs (150 × 150) and N0N1 ring pairs (160 × 160). These new datasets were reprocessed with FREALIGN (Grigorieff, 2007), first in a global alignment cycle (mode 3) and followed by three iterations of local alignments (mode 1) using featureless spheres and cylinders as initial references for the cap and N-terminal ring modules. Subsequently each module was subjected to 3D classification in mode 1, refining the angles every three iterations and performing the classification with 3 to 9 classes in order to test reproducibility. Alternatively, RELION 1.4 (Scheres, 2012) was used to accurately estimate rotational alignment accuracy as function of the box size. After post-processing, maps of the cap, N4 ring, N2N3 and N0N1 ring pairs had a resolution of 7.6, 6.7, 7.0, and 7.6 Å, respectively. Before visualization all density maps were corrected for the modulation transfer function (MTF) of the K2 direct detector. Maps were sharpened by applying a B-factor of −250 Å².

PilQ modeling

To model the structure of PilQ we first identified homologous proteins with known structure by performing a BLAST search (Altschul et al., 1990) of the PilQ Uniprot sequence (Q72IW4) against the protein data bank (PDB). Hits were filtered by query coverage of 70% and an E-value cut-off of 3.0. The resulting template PDB structures for each of the PilQ domains were pruned by structural resolution and realigned with their respective templates using the 3D-Coffee algorithm (Taly et al., 2011) to maximize alignment quality and coverage. N1 and secretin domains were realigned with their respective template sequences using HHPred (Söding et al., 2005). These alignments were then used for generating homology models for all the individual domains using Modeller (Sali and Blundell, 1993; Webb and Sali, 2014). The modelling procedure is summarized in Table 1.

Fitting homology models to the cryo-EM density map

From a series of domain deletion mutants of PilQ, we deduced that individual consecutive domains form discrete stacked ring-like density segments and are arranged with the N-terminal N0 domain at the bottom end and the C-terminal secretin domain at the top end (Salzer et al., 2016). This analysis enabled us to assign all domains in the low-resolution map of the complete PilQ assembly. Further reprocessing of EM images by recentering subdomains provided us with four better-resolved density maps corresponding to (a) secretin domain and N5 domain, (b) N4 domain, (c) N3 and N2 domains, and (d) N1 and N0 domains, respectively. The homology models of individual domains were then fitted into masked ring-like densities, first as individual protomers using the fit-in-map tool of Chimera (Pettersen et al., 2004). The best-fit protomer orientations of individual domains were identified from the 15 top scoring orientations with high cross-correlation coefficients. The best models were analyzed to ensure correct arrangements of domain boundaries for each domain within the ring structures, that is, N terminus towards the base (tail) and C terminus towards the top (cap). Connecting linker segments between domains N0N1, N2N3 and, N5-secretin were modeled using the Modeller loop modeling protocol with DOPE scoring (Sali and Blundell, 1993) at higher precision. C₁₃ symmetry-related chains (A to M) corresponding to modeled protomer fragments were added to make C₁₃-symmetric rings of PilQ domains/fragments, followed by an exhaustive 6D rigid body search to fit into the corresponding localized density regions with an angular sampling of 20 degrees using Colores (Lasker et al., 2010). To account for the flexibility and conformational dynamics of the entire PilQ oligomeric complex, we optimized the 13-mer fragment rigid-body fits using a molecular dynamics flexible fitting (MDFF) procedure (Trabuco et al., 2008). MDFF adds biasing forces proportional to the observed EM density gradients. MDFF runs were set up using VMD (Humphrey et al., 1996) and used the CHARMM36m force field (Huang et al., 2017). Simulations were carried out using NAMD (Phillips et al., 2005) in vacuum at a temperature of 300K for 500 ps. Forces corresponding to EM density gradients were coupled using an MDFF scaling factor ξ of 0.3. The MDFF runs were followed by short energy minimization runs for 10000 steps with ξ of 10. Restraints to preserve secondary structure, chirality, and trans-peptide bond geometry were employed. The loop segments within the N0 and N1 domains are longer and more disordered than in the top-scoring template structures. To model these flexible regions, ten loop conformations were modeled and filtered to identify the model exhibiting maximum overlap with the corresponding density regions in the N0N1-masked map. Four successive MDFF runs were performed for the N0N1 domains, each for 500,000 steps, using increasing values of the scaling factor, ξ = 0.3, 0.5, 1.0, and 3.0, respectively. All MDFF runs were monitored by measuring the evolution of the backbone RMSD and the cross-correlation coefficient changes with the target density maps. To optimize the atomic model of the fully assembled complex, a complete high-resolution EM map for the entire PilQ complex was built by stacking the individually processed ring density maps, guided by a low resolution template map of the whole PilQ complex. The four optimized C₁₃ symmetric models were fitted into corresponding density regions, followed by modeling of missing loop segments between N1-N2 rings, N3-N4 rings, and N4-N5 rings using Modeller (Sali and Blundell, 1993; Webb and Sali, 2014) to preserve the domain connectivity as in V. cholerae GspD. To prevent overfitting, the C₁₃ symmetry-related oligomer was used for flexible fitting of the tridecameric PilQ model into the complete EM map.

Validation of initial fits

We first generated individual domain maps of the entire PilQ protein. The initial homology domain models for all PilQ domains were then subjected to extensive sampling (n = 10000) of complete Euler angular space (ɸ ∈ [0, 2π]; θ ∈ [0, π]; Ψ ∈ [0, 2π]) to obtain initial orientation ensemble for placing the respective domains within their corresponding EM maps. Following local steepest descent optimization, we obtained several unique solutions for each domain. These unique fits were characterized by their fitted-orientation, local CC and the number of runs converging into these solutions. The top best fit solutions were ranked by their local cross correlation values computed with the corresponding domain density maps using colores from the Situs package (Lasker et al., 2010). The previously obtained initial best fit orientations for each PilQ domain were re-evaluated by comparing them with top scoring solutions before and after MDFF refinement (Figure 3—figure supplements 4, 5).

Model/map correspondence

The initial best fitted C₁₃ homology models before refinement were used to generate EM maps at the same corresponding nominal resolution of the recentered EM maps. Subsequently the model-generated maps were used to compute per voxel local cross-correlation (using the vop localCorrelation module of UCSF Chimera) against the experimental ring-module EM maps to show model/map correspondence. This showed overall good agreement between the initial C₁₃ homology models and the EM maps. Regions that were not resolved in the EM maps, corresponding to linker segments and domain termini connecting individual domains, were subsequently refined using MDFF.

Data deposition

The EM maps have been deposited in the EMDB with accession codes EMD-3985 for the full complex and EMD-3995, EMD-3996, EMD-3997 and EMD-3998 for the cap, N4, N2-N3, and N0-N1 domains, respectively. The Recentering program can be downloaded from https://github.com/rkms86/REP (a copy is archived at https://github.com/elifesciences-publications/REP).

Funding Information

This paper was supported by the following grants:

http://dx.doi.org/10.13039/501100005687Goethe-Universität Frankfurt am Main to Beate Averhoff.

http://dx.doi.org/10.13039/501100004189Max-Planck-Gesellschaft to Gerhard Hummer, Werner Kühlbrandt.

http://dx.doi.org/10.13039/501100001659Deutsche Forschungsgemeinschaft Av9/6-2 to Beate Averhoff.

Acknowledgements

We thank Özkan Yildiz and Juan Castillo for computer support and Deryck J Mills and Simone Prinz for EM support. We are grateful to Martin Wilkes for discussions about sample preparation and data acquisition, to Katharina van Pee and Martin Centola for discussions about sample stabilization, and to Ahmad Reza Mehdipour for discussions on modelling and fitting. We also thank Pilar Cossio-Tejada for critical comments on the manuscript. We thank Bernd Ludwig for providing the pDM12 vector. This project was funded by the Max Planck Society and Goethe University, Frankfurt. Part of the work was supported by a grant of the German Research Foundation (DFG) Av9/6-2 to BA.

Additional information

Competing interests

Reviewing Editor, eLife.

No competing interests declared.

Author contributions

Conceptualization, Formal analysis, Writing—original draft, Writing—review and editing, Preparation of cryo-EM samples, Devising data acquisition strategy, Acquisition of EM data, Conceptualization of REP.

Writing—original draft, performed cloning, purification and biochemical analysis, Provided samples for structure determination with help from IR and support from BA.

Writing—original draft, Writing—review and editing, Performed modelling, Fitting and refinement of structural model with inputs from GH.

Wrote the REP program, Wrote the Appendix to the manuscript.

Provided samples for EM with support from BA.

Investigation.

Supervision, Funding acquisition, Writing—review and editing.

Funding acquisition, Writing—review and editing.

Conceptualization, Formal analysis, Supervision, Funding acquisition, Writing—review and editing.

Conceptualization, Funding acquisition, Formal analysis, Supervision, Writing—review and editing.

Supplementary Material

Open In Web Browser

Appendix 1

Description of the Recentering program

DOI: 10.7554/eLife.30483.022

Accurate determination of the Euler angles of each projection of a protein complex is a prerequisite for high-resolution 3D reconstruction. If the complex has different conformations or it is inherently flexible, the orientations of individual domains in a particle can be different. In that case, reconstruction of subdomains by residual signal subtraction (RSS) can be used. However, if the center of the subdomain is far away from the center of the particle box, the alignment becomes increasingly inaccurate with the distance from the center of rotation. A program for recentering (REP) was written to move the center of mass of the subdomain to coincide with the center of rotation (Sanchez, 2017). This makes it possible to reduce the box size of the subdomain images, which maximizes local alignment and reduces computational cost.

The center of a protein complex with N subdomains can be defined as follows.

where Cp is the center of the protein complex with total mass Mp and ciis the center of the i-th subdomain with mass mi. For every projection k the center Cp(k) can be written in terms of deviations of the mean positions of the N subdomains

where Cp(k) is the center of a protein complex for the (k-th)-projection; Δci(k)is a shift of the i-th subdomain from its center ci. This shift reduces the accuracy of the registration (particle alignment) and eventually the resolution of the final 3D reconstruction. To solve this problem, (Ilca et al., 2015) focused the registration step on the i-th subdomain of the protein complex independently to find its particular shift and orientation. The proposed method extracts the projections of the i-th subdomain from the original stack (using for example the residual signal subtraction method of Bai et al., [Bai et al., 2015]) and stores them in a new stack with a new box size centered around the i-th subdomain. Reprojection of the center of the i-th subdomain, from the complete reconstructed complex (with a standard projection matching method) on to the (k-th)-projection, can be obtained by recentering the i-th subdomain.

where ci(k) is the center of the i-th subdomain; Rp(k)is the orientation of the protein complex, Ĉp(k)is its estimated center, for the (k-th)-projection and ci is the estimated center of the i-th subunit. The proj2D[] is an operator that projects 3D points to the XY plane.

The above procedure is implemented in a stand-alone program REP which reads STAR file formats as input that contain the alignment information derived from the consensus refinement, an MRCS file associated with the STAR file, and the subdomain centers from reconstruction ci. REP generates a new STAR file with its corresponding MRCS file as output. The coordinates of the new center can be obtained from UCSF Chimera (Pettersen et al., 2004) or IMOD (Kremer et al., 1996). The alignment parameters, in a STAR file format, are automatically passed on to RELION (Scheres, 2012) or can be generated quickly using various cryo-EM software packages like EMAN2 (Tang et al., 2007), Frealign (Grigorieff, 2007) and Imagic (van Heel et al., 1996). The standalone program REP is written in C++ and can be downloaded from https://github.com/rkms86/REP (Sanchez, 2017). A copy is archived at https://github.com/elifesciences-publications/REP.

Contributor Information

Wesley I Sundquist, Reviewing Editor, University of Utah School of MedicineUnited States.

Articles from eLife are provided here courtesy of eLife Sciences Publications, Ltd

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ, authors. 1990;Basic local alignment search tool. Journal of Molecular Biology. 215:403–410. [Cross Ref] [PubMed]

Averhoff B, author. 2009;Shuffling genes around in hot environments: the unique DNA transporter of Thermus thermophilus. FEMS Microbiology Reviews. 33:611–626. [Cross Ref] [PubMed]

Ayers M, Howell PL, Burrows LL, authors. 2010;Architecture of the type II secretion and type IV pilus machineries. Future Microbiology. 5:1203–1218. [Cross Ref] [PubMed]

Bai XC, Rajendra E, Yang G, Shi Y, Scheres SH, authors. 2015;Sampling the conformational space of the catalytic subunit of human γ-secretase. eLife. 4:e11182. [Cross Ref] [PubMed]

Bayan N, Guilvout I, Pugsley AP, authors. 2006;Secretins take shape. Molecular Microbiology. 60:1–4. [Cross Ref] [PubMed]

Berry JL, Phelan MM, Collins RF, Adomavicius T, Tønjum T, Frye SA, Bird L, Owens R, Ford RC, Lian LY, Derrick JP, authors. 2012;Structure and assembly of a trans-periplasmic channel for type IV pili in Neisseria meningitidis. PLoS Pathogens. 8:e1002923. [Cross Ref] [PubMed]

Blokesch M, author. 2017;In and out-contribution of natural transformation to the shuffling of large genomic regions. Current Opinion in Microbiology. 38:22–29. [Cross Ref] [PubMed]

Burkhardt J, Vonck J, Averhoff B, authors. 2011;Structure and function of PilQ, a secretin of the DNA transporter from the thermophilic bacterium Thermus thermophilus HB27. Journal of Biological Chemistry. 286:9977–9984. [Cross Ref] [PubMed]

Burkhardt J, Vonck J, Langer JD, Salzer R, Averhoff B, authors. 2012;Unusual N-terminal ααβαββα fold of PilQ from Thermus thermophilus mediates ring formation and is essential for piliation. Journal of Biological Chemistry. 287:8484–8494. [Cross Ref] [PubMed]

Böttcher B, Prazak V, Rasmussen A, Black SS, Rasmussen T, authors. 2015;The Structure of YnaI Implies Structural and Mechanistic Conservation in the MscS Family of Mechanosensitive Channels. Structure. 23:1705–1714. [Cross Ref] [PubMed]

Chen I, Christie PJ, Dubnau D, authors. 2005;The ins and outs of DNA transfer in bacteria. Science. 310:1456–1460. [Cross Ref] [PubMed]

Collins RF, Davidsen L, Derrick JP, Ford RC, Tønjum T, authors. 2001;Analysis of the PilQ secretin from Neisseria meningitidis by transmission electron microscopy reveals a dodecameric quaternary structure. Journal of Bacteriology. 183:3825–3832. [Cross Ref] [PubMed]

Costa TR, Felisberto-Rodrigues C, Meir A, Prevost MS, Redzej A, Trokter M, Waksman G, authors. 2015;Secretion systems in Gram-negative bacteria: structural and mechanistic insights. Nature Reviews Microbiology. 13:343–359. [Cross Ref] [PubMed]

Dube P, Tavares P, Lurz R, van Heel M, authors. 1993;The portal protein of bacteriophage SPP1: a DNA pump with 13-fold symmetry. The EMBO Journal. 12:1303–1309. [PubMed]

Friedrich A, Hartsch T, Averhoff B, authors. 2001;Natural transformation in mesophilic and thermophilic bacteria: identification and characterization of novel, closely related competence genes in Acinetobacter sp. strain BD413 and Thermus thermophilus HB27. Applied and Environmental Microbiology. 67:3140–3148. [Cross Ref] [PubMed]

Friedrich A, Prust C, Hartsch T, Henne A, Averhoff B, authors. 2002;Molecular analyses of the natural transformation machinery and identification of pilus structures in the extremely thermophilic bacterium Thermus thermophilus strain HB27. Applied and Environmental Microbiology. 68:745–755. [Cross Ref] [PubMed]

Gold VA, Salzer R, Averhoff B, Kühlbrandt W, authors. 2015;Structure of a type IV pilus machinery in the open and closed state. eLife. 4:e07380. [Cross Ref] [PubMed]

Grant T, Grigorieff N, authors. 2015;Measuring the optimal exposure for single particle cryo-EM using a 2.6 Å reconstruction of rotavirus VP6. eLife. 4:e06980. [Cross Ref] [PubMed]

Grigorieff N, author. 2007;FREALIGN: high-resolution refinement of single particle structures. Journal of Structural Biology. 157:117–125. [Cross Ref] [PubMed]

Hobbs M, Mattick JS, authors. 1993;Common components in the assembly of type 4 fimbriae, DNA transfer systems, filamentous phage and protein-secretion apparatus: a general system for the formation of surface-associated protein complexes. Molecular Microbiology. 10:233–243. [Cross Ref] [PubMed]

Hou L, Wan W, Lee K-H, Hwang J-N, Okopal G, Pitton J, authors. 2017;Robust Human Tracking Based on DPM Constrained Multiple-Kernel from a Moving Camera. Journal of Signal Processing Systems. 86:27–39. [Cross Ref]

Huang J, Rauscher S, Nawrocki G, Ran T, Feig M, de Groot BL, Grubmüller H, MacKerell AD, authors. 2017;CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nature Methods. 14:71–73. [Cross Ref] [PubMed]

Humphrey W, Dalke A, Schulten K, authors. 1996;VMD: visual molecular dynamics. Journal of Molecular Graphics. 14:33–38. [Cross Ref] [PubMed]

Ilca SL, Kotecha A, Sun X, Poranen MM, Stuart DI, Huiskonen JT, authors. 2015;Localized reconstruction of subunits from electron cryomicroscopy images of macromolecular complexes. Nature Communications. 6:8843. [Cross Ref] [PubMed]

Jain S, Mościcka KB, Bos MP, Pachulec E, Stuart MC, Keegstra W, Boekema EJ, van der Does C, authors. 2011;Structural characterization of outer membrane components of the type IV pili system in pathogenic Neisseria. PLoS One. 6:e16624. [Cross Ref] [PubMed]

Koo J, Lamers RP, Rubinstein JL, Burrows LL, Howell PL, authors. 2016;Structure of the Pseudomonas aeruginosa Type IVa Pilus Secretin at 7.4 Å. Structure. 24:1778–1787. [Cross Ref] [PubMed]

Koomey M, author. 1998;Competence for natural transformation in Neisseria gonorrhoeae: a model system for studies of horizontal gene transfer. APMIS. 84:56–61. [Cross Ref] [PubMed]

Korotkov KV, Gonen T, Hol WG, authors. 2011;Secretins: dynamic channels for protein transport across membranes. Trends in Biochemical Sciences. 36:433–443. [Cross Ref] [PubMed]

Kremer JR, Mastronarde DN, McIntosh JR, authors. 1996;Computer visualization of three-dimensional image data using IMOD. Journal of Structural Biology. 116:71–76. [Cross Ref] [PubMed]

Lasker K, Sali A, Wolfson HJ, authors. 2010;Determining macromolecular assembly structures by molecular docking and fitting into an electron density map. Proteins: Structure, Function, and Bioinformatics. 78:3205–3211. [Cross Ref] [PubMed]

Li X, Mooney P, Zheng S, Booth CR, Braunfeld MB, Gubbens S, Agard DA, Cheng Y, authors. 2013;Electron counting and beam-induced motion correction enable near-atomic-resolution single-particle cryo-EM. Nature Methods. 10:584–590. [Cross Ref] [PubMed]

Low HH, Gubellini F, Rivera-Calzada A, Braun N, Connery S, Dujeancourt A, Lu F, Redzej A, Fronzes R, Orlova EV, Waksman G, authors. 2014;Structure of a type IV secretion system. Nature. 508:550–553. [Cross Ref] [PubMed]

Ludtke SJ, Baldwin PR, Chiu W, authors. 1999;EMAN: semiautomated software for high-resolution single-particle reconstructions. Journal of Structural Biology. 128:82–97. [Cross Ref] [PubMed]

Marles-Wright J, Grant T, Delumeau O, van Duinen G, Firbank SJ, Lewis PJ, Murray JW, Newman JA, Quin MB, Race PR, Rohou A, Tichelaar W, van Heel M, Lewis RJ, authors. 2008;Molecular architecture of the "stressosome," a signal integration and transduction hub. Science. 322:92–96. [Cross Ref] [PubMed]

Mell JC, Redfield RJ, authors. 2014;Natural competence and the evolution of DNA uptake specificity. Journal of Bacteriology. 196:1471–1483. [Cross Ref] [PubMed]

Nguyen TH, Galej WP, Bai XC, Savva CG, Newman AJ, Scheres SH, Nagai K, authors. 2015;The architecture of the spliceosomal U4/U6.U5 tri-snRNP. Nature. 523:47–52. [Cross Ref] [PubMed]

Pang Y, Bai XC, Yan C, Hao Q, Chen Z, Wang JW, Scheres SH, Shi Y, authors. 2015;Structure of the apoptosome: mechanistic insights into activation of an initiator caspase from Drosophila. Genes & Development. 29:277–287. [Cross Ref] [PubMed]

Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE, authors. 2004;UCSF Chimera--a visualization system for exploratory research and analysis. Journal of Computational Chemistry. 25:1605–1612. [Cross Ref] [PubMed]

Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kalé L, Schulten K, authors. 2005;Scalable molecular dynamics with NAMD. Journal of Computational Chemistry. 26:1781–1802. [Cross Ref] [PubMed]

Reichow SL, Korotkov KV, Hol WG, Gonen T, authors. 2010;Structure of the cholera toxin secretion channel in its closed state. Nature Structural & Molecular Biology. 17:1226–1232. [Cross Ref] [PubMed]

Rohou A, Grigorieff N, authors. 2015;CTFFIND4: Fast and accurate defocus estimation from electron micrographs. Journal of Structural Biology. 192:216–221. [Cross Ref] [PubMed]

Rose I, Biuković G, Aderhold P, Müller V, Grüber G, Averhoff B, authors. 2011;Identification and characterization of a unique, zinc-containing transport ATPase essential for natural transformation in Thermus thermophilus HB27. Extremophiles. 15:191–202. [Cross Ref] [PubMed]

Rumszauer J, Schwarzenlander C, Averhoff B, authors. 2006;Identification, subcellular localization and functional interactions of PilMNOWQ and PilA4 involved in transformation competency and pilus biogenesis in the thermophilic bacterium Thermus thermophilus HB27. FEBS Journal. 273:3261–3272. [Cross Ref] [PubMed]

Sali A, Blundell TL, authors. 1993;Comparative protein modelling by satisfaction of spatial restraints. Journal of Molecular Biology. 234:779–815. [Cross Ref] [PubMed]

Salzer R, D'Imprima E, Gold VA, Rose I, Drechsler M, Vonck J, Averhoff B, authors. 2016;Topology and Structure/Function Correlation of Ring- and Gate-forming Domains in the Dynamic Secretin Complex of Thermus thermophilus. Journal of Biological Chemistry. 291:14448–14456. [Cross Ref] [PubMed]

Salzer R, Herzberg M, Nies DH, Joos F, Rathmann B, Thielmann Y, Averhoff B, authors. 2014;Zinc and ATP binding of the hexameric AAA-ATPase PilF from Thermus thermophilus: role in complex stability, piliation, adhesion, twitching motility, and natural transformation. The Journal of Biological Chemistry. 289:30343–30354. [Cross Ref] [PubMed]

Sanchez R, author. 2017. Recentering and subboxing of particles (REP)Github. 1.0https://github.com/rkms86/REP

Scheres SH, author. 2012;RELION: implementation of a Bayesian approach to cryo-EM structure determination. Journal of Structural Biology. 180:519–530. [Cross Ref] [PubMed]

Schraidt O, Marlovits TC, authors. 2011;Three-dimensional model of Salmonella's needle complex at subnanometer resolution. Science. 331:1192–1195. [Cross Ref] [PubMed]

Schwarzenlander C, Haase W, Averhoff B, authors. 2009;The role of single subunits of the DNA transport machinery of Thermus thermophilus HB27 in DNA binding and transport. Environmental Microbiology. 11:801–808. [Cross Ref] [PubMed]

Sorzano CO, Bilbao-Castro JR, Shkolnisky Y, Alcorlo M, Melero R, Caffarena-Fernández G, Li M, Xu G, Marabini R, Carazo JM, authors. 2010;A clustering approach to multireference alignment of single-particle projections in electron microscopy. Journal of Structural Biology. 171:197–206. [Cross Ref] [PubMed]

Sousa JS, Mills DJ, Vonck J, Kühlbrandt W, authors. 2016;Functional asymmetry and electron flow in the bovine respirasome. eLife. 5:e21290. [Cross Ref] [PubMed]

Spreter T, Yip CK, Sanowar S, André I, Kimbrough TG, Vuckovic M, Pfuetzner RA, Deng W, Yu AC, Finlay BB, Baker D, Miller SI, Strynadka NC, authors. 2009;A conserved structural motif mediates formation of the periplasmic rings in the type III secretion system. Nature Structural & Molecular Biology. 16:468–476. [Cross Ref] [PubMed]

Söding J, Biegert A, Lupas AN, authors. 2005;The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Research. 33:W244–W248. [Cross Ref] [PubMed]

Taly JF, Magis C, Bussotti G, Chang JM, Di Tommaso P, Erb I, Espinosa-Carrasco J, Kemena C, Notredame C, authors. 2011;Using the T-Coffee package to build multiple sequence alignments of protein, RNA, DNA sequences and 3D structures. Nature Protocols. 6:1669–1682. [Cross Ref] [PubMed]

Tang G, Peng L, Baldwin PR, Mann DS, Jiang W, Rees I, Ludtke SJ, authors. 2007;EMAN2: an extensible image processing suite for electron microscopy. Journal of Structural Biology. 157:38–46. [Cross Ref] [PubMed]

Tosi T, Estrozi LF, Job V, Guilvout I, Pugsley AP, Schoehn G, Dessen A, authors. 2014;Structural similarity of secretins from type II and type III secretion systems. Structure. 22:1348–1355. [Cross Ref] [PubMed]

Trabuco LG, Villa E, Mitra K, Frank J, Schulten K, authors. 2008;Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics. Structure. 16:673–683. [Cross Ref] [PubMed]

Tsai CL, Tainer JA, authors. 2016;The ATPase motor turns for type IV Pilus assembly. Structure. 24:1857–1859. [Cross Ref] [PubMed]

van Heel M, Harauz G, Orlova EV, Schmidt R, Schatz M, authors. 1996;A new generation of the IMAGIC image processing system. Journal of Structural Biology. 116:17–24. [Cross Ref] [PubMed]

Webb B, Sali A, authors. 2014;Comparative protein structure modeling using MODELLER. Current Protocols in Bioinformatics. 47:1–32. [Cross Ref] [PubMed]

Worrall LJ, Hong C, Vuckovic M, Deng W, Bergeron JR, Majewski DD, Huang RK, Spreter T, Finlay BB, Yu Z, Strynadka NC, authors. 2016;Near-atomic-resolution cryo-EM analysis of the Salmonella T3S injectisome basal body. Nature. 540:597–601. [Cross Ref] [PubMed]

Yan Z, Yin M, Xu D, Zhu Y, Li X, authors. 2017;Structural insights into the secretin translocation channel in the type II secretion system. Nature Structural & Molecular Biology. 24:177–183. [Cross Ref] [PubMed]

eLife. 6: e30483  » Sub-article

Decision letter

Wesley I Sundquist, Reviewing Editor

DOI: 10.7554/eLife.30483.036

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

Thank you for submitting your article "Cryo-EM structure of the bifunctional secretin complex of Thermus thermophilus" for consideration by eLife. Your article has been favorably evaluated by Richard Aldrich (Senior Editor) and three reviewers, one of whom, Wesley I Sundquist (Reviewer #1), is a member of our Board of Reviewing Editors.

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

The manuscript by D'Imprima et al. describes a cryoEM structure of a PilQ type IV pilus secretin complex from Thermus thermophilus at ~7 Å resolution. The resolution of the reconstruction was greatly improved using an in-house image processing tool, REP, that allowed the complex to be segmented into individual rigid components that were each processed separately. The T. thermophilus secretin complex has an especially large number of periplasmic domains, N0-N5, all of which are resolved in the EM map. Homology models of individual domains were generated and fit into the EM map. Linkers between several of the domain rings are well-resolved in the structure. A "crown" of density is apparent on the extracellular side of the complex that is not accounted for by the PilQ sequence. This density was not apparent in a PilQ variant with substitutions in two basic residues at the interface with the crown; this variant had a mild motility phenotype and a severe transformation phenotype suggesting that these residues and perhaps the crown itself are involved in DNA uptake. The secretin structure is of high quality and the results related to the crown are intriguing. Overall the work appears well done and paper is clearly and concisely written, but there are still some issues that need to be addressed before the manuscript will be acceptable for publication.

Major issues that need to be addressed:

1) A significant shortcoming is that the authors have not identified the protein that comprises the crown. This should be feasible by mass spec/proteomics approaches and seems necessary to complete the study, particularly as the information could substantially augment the existing structure and/or functional interpretation if the missing component proved to be homologous to another protein of known structure or function (e.g., a DNA binding protein). At a minimum, if this is not possible, then the authors should explain the lengths they have taken to identify this protein and why they have not been successful.

2) Data quality and assessments of the modeling/fitting. Firstly, the authors should explain the reason(s) that a high resolution structure could not be obtained. Secondly, docking at 7A resolution can be quite inaccurate and needs the maximum possible validation. Specifically, the authors need to provide more information regarding the accuracy of the homology models and how well they fit into the cryoEM maps (e.g., please provide percent similarities for the alignments used for homology modeling, statistical analyses of the fitting of domains into the maps, close-up/extracted images of protomers with more than one view for each of the domain fits so that readers can assess their quality, and apply a consistency test by generating a 7A map from the model and comparing it with the experimental map). For example, all published high resolution secretin structures show a continuous beta-stranded belt close to the outer membrane, but in the structure described here, the EM-density is quite discontinuous, which largely contrasts the proposed atomic model (Figure 2H, Video 1, Figure 3).

Finally, the authors should also show how much the homology models diverged from the final models in cases where they were flexibly fitted into the EM map.

3) The C13 symmetry assignment needs to be solidified. Why were the eigenimages in Figure 1—figure supplement 1 performed on only a very small subset of particles (n=27)? What was the basis for selecting 27 out of 363 top-views (subsection “Cryo-EM image processing”, first paragraph). How does the (relevant) eigenimage appear if no random rotation is performed? It is expected that the second eigenimage will only show the 13-fold arrangement if all particles in top views are also not tilted. Can this be excluded? It is surprising that the 13-fold arrangement in the eigenimage (F) is observed at a level that is outside the ring. Could it be that the alternating pattern is a result of noise, tilted views, etc.? The authors should perform an eigenimage analysis on tightly ring-masked particles only, excluding all potential noise outside and inside the ring. Increasing the number of particles would also be of benefit.

4) The authors describe a novel recentering (REP) approach. Yet, how much does the result differ from Ilca et al. 2015? Did the authors do a comparative study?

Contributor Information

Wesley I Sundquist, Reviewing Editor, University of Utah School of MedicineUnited States.

eLife. 6: e30483  » Sub-article

Author response

DOI: 10.7554/eLife.30483.037

Major issues that need to be addressed:

1) A significant shortcoming is that the authors have not identified the protein that comprises the crown. This should be feasible by mass spec/proteomics approaches and seems necessary to complete the study, particularly as the information could substantially augment the existing structure and/or functional interpretation if the missing component proved to be homologous to another protein of known structure or function (e.g., a DNA binding protein). At a minimum, if this is not possible, then the authors should explain the lengths they have taken to identify this protein and why they have not been successful.

We fully agree with the reviewer that the protein that comprises the crown is of major interest. We found that the crown structure is essential for natural transformation and so we put a lot of effort into the identification of the crown protein. We started with a mass spectrometry approach of the purified PilQ complex. This led to the detection of more than 150 proteins in the purified PilQ preparations. Since the PilQ purification was performed with membrane extracts several well-known membrane associated/integrated proteins were detected, such as transporters, ATPases, the pyruvate DH complex, a cytochrome c oxidase, solute binding proteins, porins, lipoproteins and others. However, these studies did not reveal a good candidate for the crown protein.

The mass spectrometry analysis also led to the detection of 50 uncharacterized or hypothetical proteins. Several of them were closely associated with pilin-like genes or known genes of the natural transformation machinery and we speculate that one of these genes might encode the crown protein. We selected five of these uncharacterized/hypothetical genes, linked to pilins or genes of the transformation machinery and generated deletion mutants. However, none of these mutants were defective in natural transformation, which indicates that these proteins might not be directly involved in crown-mediated natural transformation.

We analyzed all proteins from T. thermophilus having sequence signatures resembling the V. cholerae GspD lip domain. We found one hit with a hypothetical membrane protein. This was also detected as part of a mass-spectrometry analysis of PilQ complex preparations. However, deletion of this candidate gene did not affect transformation, which indicates that this protein is not involved directly in DNA binding or natural transformation mediated by the crown.

2) Data quality and assessments of the modeling/fitting. Firstly, the authors should explain the reason(s) that a high resolution structure could not be obtained.

The organization of seven individual PilQ domains into several rigid ring structures connected by flexible linkers segments makes this complex highly dynamic. The increased flexibility enhances inter-ring motions, resulting in several dynamic states. The stabilization of these states depends on several external factors such as buffer conditions/temperature and imaging parameters. Our approach of reconstructing of individual rigid modules improved the resolution dramatically from 20 to 7 Å, and an even bigger improvement could not be expected.

Secondly, docking at 7A resolution can be quite inaccurate and needs the maximum possible validation. Specifically, the authors need to provide more information regarding the accuracy of the homology models and how well they fit into the cryoEM maps (e.g., please provide percent similarities for the alignments used for homology modeling, statistical analyses of the fitting of domains into the maps, close-up/extracted images of protomers with more than one view for each of the domain fits so that readers can assess their quality, and apply a consistency test by generating a 7A map from the model and comparing it with the experimental map). For example, all published high resolution secretin structures show a continuous beta-stranded belt close to the outer membrane, but in the structure described here, the EM-density is quite discontinuous, which largely contrasts the proposed atomic model (Figure 2H, Video 1, Figure 3).

We are aware of the challenges in docking/fitting homology models into EM maps, both because of the quality of the structural models and because of the resolution of the map. To address this issue, we now (i) list percent of the similarities of the alignments used for homology modelling, (ii) provide detailed statistical analysis of the fitting of domains into the maps, and (iii) show closer views of the models in the maps.

To assess the statistics of the fits, we compared the cross-correlation coefficient (CC) after fitting of PilQ domains to the distribution of the CCs obtained in extensive sampling of orientations of the initial model (n = 10000; complete ensemble of fits). We show that the best fitting orientations in this initial search rank within the top five best solutions from the complete ensemble of fits. This implies that the initial best fit orientations had good coverage. To assess convergence, we checked that repeated optimizations produced identical best-fit orientations. We summarize the results of fitting procedures by showing the cumulative distribution of the entire ensemble of fits for each domain along with the CC-value of its initial chosen best fit before and after refinement using MDFF approach. We have added this information to the Materials and methods and also added two new figures (Figure 3—figure supplement 4, 5) and a new Table 2 summarizing these results. We also provide close-up images of initial and refined fits in multiple orientations for all PilQ domains (Figure 3—figure supplement 6). Further we have also added a video (Video 2) showing the final atomic model fit into the density maps. The movie displays individual domain fits in close-up by sectioning through the various ring modules of the entire PilQ complex.

Modelling details:

Table 1 (modified now) also provides information about the homology modelling procedure of individual domain models. It clearly provides information about individual domain templates, their sequence identities/similarities with the target PilQ domain along with individual DOPE scores/Z-scores.

Model/map validation:

The initial best fitted C-13 homology models before refinement were used to generate model derived EM maps at the same nominal resolution as experimental recentered EM maps. Comparison of model derived and experimental EM maps was performed by computing per voxel local cross-correlation (using the vop localCorrelation module of UCSF Chimera) to show model/map correspondence. This showed an overall good agreement between the initial C-13 symmetric homology models and the EM maps. Regions that were not resolved in original EM maps, i.e. corresponding to linker segments and domain termini connecting individual maps were subsequently refined using MDFF.

Discontinuity of beta-barrel density:

All published high resolution structures of Secretins display a continuous beta-belt density corresponding to the outer beta barrel. Beta sheets are only resolved in to individual strand at resolution higher than 4.7Å and are hard to interpret at lower resolution. The discontinuous density in the beta-barrel region of our map (Figure 2H and Figure 3) is highly dependent on the contour level selection. In Figure 2H and Figure 3, the contour level of the map was chosen to highlight features of the crown and N5 regions, which would otherwise not be clear. The beta barrel-region has a continuous density when visualized at a lower contour level. Therefore, our homology model, which is quite similar to other well characterized secretins (Figure 4), is compatible with the current map. We have made an additional figure (Author response image 1) to highlight the continuity/discontinuity as a function of the contour level.

DOI: 10.7554/eLife.30483.024 Author response image 1. Author response image 1. View larger version

Finally, the authors should also show how much the homology models diverged from the final models in cases where they were flexibly fitted into the EM map.

Structural divergence after MDFF:

MDFF refinement adds forces derived from density gradients in the EM maps to a molecular simulation model. Constraints on the secondary structure and chirality of residues ensure that the model structure is changed minimally to avoid overfitting. The variation of the structure after refinement is therefore a measure of the model divergence required to obtain a satisfactory fit (better/improved CC). We find that the average Cα-RMSD for individual domain models before and after MDFF is ~2.7 Å, indicating that the divergence of domain structures during refinement procedure was minimal and restricted to long loops and termini only. The largest change was observed for the N1 domain during refinement (4.04 Å), which contains several long loops connecting individual secondary structure elements. We have added another column to Table 2 providing this information.

3) The C13 symmetry assignment needs to be solidified. Why were the eigenimages in Figure 1—figure supplement 1 performed on only a very small subset of particles (n=27)? What was the basis for selecting 27 out of 363 top-views (subsection “Cryo-EM image processing”, first paragraph).

Out of 363 top-views 27 where chosen because they seem not to be tilted in respect to the imaging plane, therefore ensuring the interpretability of the eigenimages.

How does the (relevant) eigenimage appear if no random rotation is performed?

The results are consistent. Multiple copies randomly rotated are generated in order to get more robust statistics, see Böttcher et al., Structure 23:1705-14, 2015.

It is expected that the second eigenimage will only show the 13-fold arrangement if all particles in top views are also not tilted. Can this be excluded?

The reviewer is correct, for this reason we selected only the 27 particles that look untilted. The fact that we do detect 13-fold symmetry shows that this assumption was mostly correct.

It is surprising that the 13-fold arrangement in the eigenimage (F) is observed at a level that is outside the ring. Could it be that the alternating pattern is a result of noise, tilted views, etc.? The authors should perform an eigenimage analysis on tightly ring-masked particles only, excluding all potential noise outside and inside the ring. Increasing the number of particles would also be of benefit.

We made sure that the mask which we used to calculate the principal components was not too tight and supply a figure to show this (Author response image 2). The circular mask applied during 2D classification includes some noise surrounding the particle images (A). However, the binary annular mask used during multivariate statistical analysis (B) which in (C) is multiplied by the class average clearly show that it excludes the noise outside the particle ring and just includes the symmetric spikes radiating from the particle top views. The same mask is shown multiplied by the second eigenimage (D) which corresponds to Figure 1—figure supplement 1F.

The fact that ultimately we obtained a three-dimensional structure showing secondary structure validates our assessment of the structure as 13-fold symmetric.

Author response image 2.View larger versionAuthor response image 2. Assessment of PilQ symmetry and oligomeric state.(A) 2D class average derived from 27 top views was not sufficient to assess the symmetry unambiguously. (B) Binary annular mask applied for multivariate statistical analysis. (C) Multiplication of (A) by (B) shows that the applied mask is not too lose and it just includes signal from the ring densities. (D) Second eigenimage after multivariate statistical analysis.

4) The authors describe a novel recentering (REP) approach. Yet, how much does the result differ from Ilca et al. 2015? Did the authors do a comparative study?

Our idea was to provide a method which is simple, flexible, modular and fast.

The differences between REP and the program from Ilca et al. 2015 relates to two core points: what the program actually does (1) and how it operates (2).

(1) Ilca's localized reconstruction tool provides a final refinement step to a specific area of the 3D map. Its general procedure consists of four steps: a initialization step which includes preparation of the necessary metadata for the localized reconstruction and the Residual Signal Substraction (RSS) procedure. This is followed by extraction of the subparticles and a final reconstruction with refinement.

On the other hand, our implementation, REP, is designed as a separate processing step of the whole reconstruction procedure. We only compute the shifts and orientations of the subparticles followed by extraction. This way we can apply the recentering procedure to the data whether this has been already processed with RSS or not. The modularity of REP allows extra flexibility, as the subtraction procedure and final refinement steps are optional and can be performed by any processing software package.

(2) The current implementation of Ilca's localized reconstruction uses Python as programming language, and it requires the installation of the Scipion and Relion-1.4 Python libraries.

We think that this is a drawback, as it forces the user to adopt specific programs which may not even be compatible with newer versions. With this in mind, we decided to implement REP using C++ and without using any additional library. By doing this we ensure that we can always use the state-of-the-art algorithms for signal subtraction, reconstruction and refinement provided by any image processing package. The concept of recentering is not a novelty per-se. The corresponding author of Ilca et al. 2015 was for example already co-author of one of the first examples, to our knowledge, of a similar approach (Briggs et al., J Struct Biol.150:332-9, 2005).

We did not do a comparative study since the program from Ilca et al. 2015 makes use of the same concept of particle recentering. The main novelty of our approach is that REP is a stand-alone program completely independent from any other reconstruction software.

[Back]

Figure 1. Electron micrograph of PilQ complex and selected 2D class averages.

Typical micrograph of purified WT PilQ recorded at 1.8 μm defocus. Individual PilQ particles are clearly recognizable, either as frequent side-views (blue boxes) or less frequent top-views (green boxes). Insets show a selection of 2D class averages. Class I to IV were obtained by averaging 1096, 935, 728, and 691 particles respectively. Class II (2x magnification) shows the salient features of the PilQ complex with Gate 1 and 2 (blue lines), rings N5 to N0 (red lines) forming the tail (dark blue bracket) and the cap (pink bracket) containing the secretin domain (yellow bracket). Classes III and IV indicate a rocking motion of the cap and N5 ring relative to the long axis of PilQ, with a hinge at the N4 ring. Green arrow, detergent micelle (DDM).

Figure 1—figure supplement 1.View larger versionFigure 1—figure supplement 1. Assessment of PilQ symmetry and oligomeric state.(A) 2D class average of 27 top views, which was not sufficiently clear for an unambiguous assessment of rotational symmetry. (B) The same 2D class average with 13-fold rotation symmetry applied. (C) Applying 12-fold symmetry (Collins et al., 2001) resulted in a featureless ring. (D) Low-pass filtering to 4 Å enhances the contrast of a typical top view of the complex. Randomly rotated copies (n = 50) of the top view in (D) were subjected to multivariate statistical analysis. An annular mask excluded the noise and limited the analysis to the features of the ring. (E) The first eigenimage, showing the average of all images and the annular mask. (F) The second eigenimage shows 13 distinct positive (blue arrowheads) and negative features, indicating 13-fold rotational symmetry.

[Back]

[Back]

Figure 1—figure supplement 1. Assessment of PilQ symmetry and oligomeric state.

(A) 2D class average of 27 top views, which was not sufficiently clear for an unambiguous assessment of rotational symmetry. (B) The same 2D class average with 13-fold rotation symmetry applied. (C) Applying 12-fold symmetry (Collins et al., 2001) resulted in a featureless ring. (D) Low-pass filtering to 4 Å enhances the contrast of a typical top view of the complex. Randomly rotated copies (n = 50) of the top view in (D) were subjected to multivariate statistical analysis. An annular mask excluded the noise and limited the analysis to the features of the ring. (E) The first eigenimage, showing the average of all images and the annular mask. (F) The second eigenimage shows 13 distinct positive (blue arrowheads) and negative features, indicating 13-fold rotational symmetry.

[Back]

[Back]

Figure 2. Effects of recentering in 2D and 3D reconstruction.

(A) Typical 2D class average of 600 PilQ particles. (B–E) Reference-free 2D classification of the same particle subset after residual signal subtraction (RSS) and recentering, which enhances high-resolution features (red arrows) in the cap (B), N4 ring (C), N2N3 ring pair (D), and the N0N1 ring pair (E). (F–H) The improvement upon applying a combination of RSS and recentering is evident after 3D refinement. (F) Standard reconstruction yields a map of ~20 Å resolution. (G) Reconstructing the cap module (pink) and tail module (blue) separately using the RSS procedure improves the resolution to 15 Å. (H) Combination of RSS and recentering after dividing the PilQ complex into four modules doubled the resolution to 7 Å. Images were realigned, classified and refined independently for the cap and N5 ring (1763 particles, orange), N4 ring (3888 particles, yellow), N2N3 ring pair (1585 particles, light green) and N0N1 ring pair (2685 particles, blue). Secondary structure features were resolved in all component maps, revealing the single-polypeptide chain connections between the N0N1 and N2N3 ring pairs. (I) The recentering procedure does not introduce overfitting, as independently aligned and refined modules (N0N1 and N2N3) display consistent features. Features of the N2 ring along with N1-N2 linker region from both modules (blue and green) show maximal overlap. The position of the outer membrane (OM) is indicated by a transparent grey bar.

Figure 2—figure supplement 1.View larger versionFigure 2—figure supplement 1. FSC curves and refinement statistics.FSC curves (top panel) of independently aligned and refined secretin-N5 rings (orange), N4 rings (yellow), the N2N3 (green) and N0N1 ring pairs (blue) after applying a combination of RSS and recentering. The dashed line marks the 0.143 FSC threshold. The lower panel shows 3D reconstruction statistics.

[Back]

[Back]

Figure 2—figure supplement 1. FSC curves and refinement statistics.

FSC curves (top panel) of independently aligned and refined secretin-N5 rings (orange), N4 rings (yellow), the N2N3 (green) and N0N1 ring pairs (blue) after applying a combination of RSS and recentering. The dashed line marks the 0.143 FSC threshold. The lower panel shows 3D reconstruction statistics.

[Back]

[Back]

Figure 3. Fit of the PilQ model.

(A) Independent reconstructions of the four modules after recentering combined into a single high-resolution map (grey), as indicated by the low-resolution envelope of the entire complex (transparent light blue). (B) Atomic model of C₁₃ symmetric PilQ complex (rainbow colored by chain) fitted into maps of individual modules after MDFF refinement (Trabuco et al., 2008). (C) Cross section of (B), showing the position of gates 1 and 2 relative to the outer membrane (grey bar) and the overall fit of individual domains into each ring. The density of the crown region outside the outer membrane is not accounted for by the PilQ sequence. (D) Top and side views indicating the helical arrangement of PilQ protomers (one protomer red) around the central symmetry axis. (E) Per voxel local cross-correlation values (CC) of C₁₃ symmetric model-derived maps and REP-reconstructed EM maps (outlines) were computed and used to color (red-white-blue) the atomic models (shown as spheres) before MDFF to highlight model/map correspondence.

Figure 3—figure supplement 1.View larger versionFigure 3—figure supplement 1. Sequence alignments of PilQ domains.Homology models were based on sequence alignment of individual PilQ domains with their respective best hit PDB structure templates (Table 1). The predicted secondary structure and DSSP-assigned secondary structure for the templates are indicated by coloured bars (alpha helices: red, beta strands: green).

Figure 3—figure supplement 2.View larger versionFigure 3—figure supplement 2. Gates of PilQ complex.(A) Cross section along the central axis of the PilQ map (transparent grey). The position of Gates 1 and 2 in the PilQ complex in the sub-tomogram average of the non-piliated, closed state (EMD-3021; Gold et al., 2015) (center) is indicated. The crown module (red) and Gate 1 (orange box) are shown on the left. Gate two is formed by the N1N0 ring pair in the blue box (right). The fitted homology model is rainbow-colored by residue. (B) Top and bottom views of clipped density maps of Gate 1 (left) and Gate 2 (right). The lower panels show the map density with fitted PilQ domains. The Gate one aperture consists of bifurcated density segments that are visible in the cross sections of side views (Figure 3C). This feature is conserved in the gates of other secretins (see Figure 4—figure supplement 1). Note that the apertures of both gates consist of loop segments in the secretin domain (Gate 1) and N1 domain (Gate 2).

Figure 3—figure supplement 3.View larger versionFigure 3—figure supplement 3. Alternative linker models connecting N3 and N4 rings.PilQ domains N3 and N4 are connected by a 17-residue linker. There is no clear map density between N3 and N4, but the inter-ring distance constrains the possible connectivity between the C-terminus of N3 and the N-terminus of N4. Domains are most likely offset by one or two along the ring, resulting in alternative models with either right-handed (green dashed lines) or left-handed helicity (red dashed lines). The connection is closest in the model with a domain offset of one and right-handed helicity (green linker segment).

Figure 3—figure supplement 4.View larger versionFigure 3—figure supplement 4. Assessment of the PilQ domain fits.(A) Location of individual PilQ domains along a single protomer showing the best fit models before (grey) and after MDFF refinement (cyan). (B) Individual PilQ protomer domain models were used in exhaustive Euler angle sampling to obtain multiple fits (red, n = 10000) ranked by their cross-correlation (CC) values. The red lines indicate the mean CC value along with the standard error of the mean (green). The best initial fits (grey lines) span the protomer domain density and have the highest CC values. Upon MDFF refinement of the initial fits CC values increase (cyan lines).

Figure 3—figure supplement 5.View larger versionFigure 3—figure supplement 5. Model refinement.Initial best fits (grey) of PilQ domains into corresponding ring-modules density maps after 6D search. The localization information of the individual domains within the ring modules was obtained from domain deletion mutations described previously (Salzer et al., 2016). The initial fits were used to build C₁₃ symmetric atomic models of the individual rings followed by refinement using a Molecular Dynamics Flexible Fitting procedure (MDFF) to obtain the final model (cyan). The refinement procedure improved the overall fit into the individual maps with minor or no changes to the central domain scaffolds, thereby confirming that the folds derived from template structures are preserved in PilQ. The N- and C-termini of the individual domains of the N0N1 fragment undergo local rearrangement to fit better into the density map during the MDFF runs, despite using a weak coupling of MD forces to the map derived density gradient (g = 0.3).

[Back]

[Back]

Figure 3—figure supplement 1. Sequence alignments of PilQ domains.

Homology models were based on sequence alignment of individual PilQ domains with their respective best hit PDB structure templates (Table 1). The predicted secondary structure and DSSP-assigned secondary structure for the templates are indicated by coloured bars (alpha helices: red, beta strands: green).

[Back]

[Back]

Figure 3—figure supplement 2. Gates of PilQ complex.

(A) Cross section along the central axis of the PilQ map (transparent grey). The position of Gates 1 and 2 in the PilQ complex in the sub-tomogram average of the non-piliated, closed state (EMD-3021; Gold et al., 2015) (center) is indicated. The crown module (red) and Gate 1 (orange box) are shown on the left. Gate two is formed by the N1N0 ring pair in the blue box (right). The fitted homology model is rainbow-colored by residue. (B) Top and bottom views of clipped density maps of Gate 1 (left) and Gate 2 (right). The lower panels show the map density with fitted PilQ domains. The Gate one aperture consists of bifurcated density segments that are visible in the cross sections of side views (Figure 3C). This feature is conserved in the gates of other secretins (see Figure 4—figure supplement 1). Note that the apertures of both gates consist of loop segments in the secretin domain (Gate 1) and N1 domain (Gate 2).

[Back]

[Back]

Figure 3—figure supplement 3. Alternative linker models connecting N3 and N4 rings.

PilQ domains N3 and N4 are connected by a 17-residue linker. There is no clear map density between N3 and N4, but the inter-ring distance constrains the possible connectivity between the C-terminus of N3 and the N-terminus of N4. Domains are most likely offset by one or two along the ring, resulting in alternative models with either right-handed (green dashed lines) or left-handed helicity (red dashed lines). The connection is closest in the model with a domain offset of one and right-handed helicity (green linker segment).

[Back]

[Back]

Figure 3—figure supplement 4. Assessment of the PilQ domain fits.

(A) Location of individual PilQ domains along a single protomer showing the best fit models before (grey) and after MDFF refinement (cyan). (B) Individual PilQ protomer domain models were used in exhaustive Euler angle sampling to obtain multiple fits (red, n = 10000) ranked by their cross-correlation (CC) values. The red lines indicate the mean CC value along with the standard error of the mean (green). The best initial fits (grey lines) span the protomer domain density and have the highest CC values. Upon MDFF refinement of the initial fits CC values increase (cyan lines).

[Back]

[Back]

Figure 3—figure supplement 5. Model refinement.

Initial best fits (grey) of PilQ domains into corresponding ring-modules density maps after 6D search. The localization information of the individual domains within the ring modules was obtained from domain deletion mutations described previously (Salzer et al., 2016). The initial fits were used to build C₁₃ symmetric atomic models of the individual rings followed by refinement using a Molecular Dynamics Flexible Fitting procedure (MDFF) to obtain the final model (cyan). The refinement procedure improved the overall fit into the individual maps with minor or no changes to the central domain scaffolds, thereby confirming that the folds derived from template structures are preserved in PilQ. The N- and C-termini of the individual domains of the N0N1 fragment undergo local rearrangement to fit better into the density map during the MDFF runs, despite using a weak coupling of MD forces to the map derived density gradient (g = 0.3).

[Back]

[Back]

Table 1. Modelling and fitting of PilQ domains.

Sequences of PilQ domains aligned to their respective templates were used to create atomic models of each domain. Models with the best scores (lower DOPE scores [Sali and Blundell, 1993]) were fitted into the corresponding density maps and screened for best fit orientations (highest cross correlation coefficient).

DOI: 10.7554/eLife.30483.012

PilQ
Domain
Size
(aa)
Consensus
Secondary structural motif
Template protein (PDB code)
Sequence identity/similarity
Resolution
(Å)/Method
DOPE score
(z-score)

N0
116
βαββαββα
N0 domain of PilQ from N. meningitidis
4ARO_A
11.45/
20.61
NA/NMR structure
−9457.03 (0.230)

N1
107
αββββα
Periplasmic domain of secretin EscC from enteropathogenic E. coli 3GR5_A
N1 domain model from PilQ from N. meningitidis
4AV2_A
8.77/
14.91
3gr5: 2.05/X-ray
4av2: NA/NMR, Model and cryo-EM
−7861.37
(0.351)

N2
57
βαββα
Extra membrane domain of secritin HOFQ from A. actinomycetemcomitans
2Y3M_A
9.72/
20.83
2.3/X-ray
−5616.16
(−1.684)

N3
61
βαββα
N domain of secretin XcpQ from P. aeruginosa
4E9J_A
17.28/
33.33
2.03/X-ray
−5858.94
(−1.331)

N4
56
βαββα
N domain of secretin XcpQ from P. aeruginosa
4E9J_A
14.28/
29.87
2.03/X-ray
−5261.98
(−1.392)

N5
82
βαββα
N1 domain model from PilQ from N. meningitidis
4AV2_A
18.88/
27.77
NA/NMR, Model and cryo-EM
−6353.61
(−1.415)

Secretin
230
βββββββββααβ
Secretin domain of T2SS GspD from V. cholerae
5WQ8_A
8.89/
14.18
3.26/cryo-EM
−17222.63
(0.777)

[Back]

[Back]

Table 2. Statistics of PilQ domain fits.

Homology models of the different PilQ domains were fitted into the respective ring module maps. Each domain was initially placed within the density maps in multiple orientations (n = 10,000), spanning the entire Euler angular space, and subjected to steepest descent local optimization. The resulting fits were clustered and ranked by their CC-scores with the corresponding density maps fits. The table lists the top-5 fits for each domain in the Euler-angle search, along with their CC-scores and the number of independent optimization runs ending in this cluster. The top-ranked fits were assessed for compatibility of the N-to-C-terminal domain orientations (bottom-to-top). The compatible orientations with the highest CC-value was selected as initial fit. Subsequent refinement using MDFF flexible fitting (Trabuco et al., 2008) further improved the fits with minimal structural change.

DOI: 10.7554/eLife.30483.013

PilQ domain fits
Total no. of fits
No. of unique fits
Top 5 solutions
Initial chosen
protomer fit.
CC (Cluster #)
Fit after MDFF.
CC
(RMSD, Å)
CC for C₁₃
oligomeric ring
model
CC after MDFF
run with masked
density map

Rank #
CC with map
Runs

N0
10000
152
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.7028
0.6715
0.6696
0.6681
0.6670
576
327
167
33
174
0.7028
(Rank 1)
0.7863
(3.210)
0.70
0.85

N1
10000
115
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.6673
0.6639
0.6583
0.6570
0.6551
530
682
336
110
211
0.6673
(Rank 1)
0.7682
(4.044)
0.70
0.85

N2
10000
48
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.8381
0.8270
0.8270
0.8256
0.8248
784
815
364
475
953
0.8381
(Rank 1)
0.9139
(1.883)
0.84
0.91

N3
10000
32
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.8279
0.8253
0.8249
0.8237
0.8222
473
457
889
829
408
0.8249 (Rank 3)
0.9026
(3.376)
0.84
0.91

N4
10000
32
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.7982
0.7931
0.7793
0.7521
0.7479
1106
757
483
434
511
0.7982 (Rank 1)
0.8539
(2.166)
0.77
0.93

N5
10000
65
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.7394
0.7371
0.7274
0.7218
0.7203
206
292
364
559
329
0.7203 (Rank 5)
0.7562
(2.378)
0.77
0.82

Secretin
9527
86
Rank 1
Rank 2
Rank 3
Rank 4
Rank 5
0.7634
0.7063
0.7013
0.6958
0.6899
672
156
285
104
370
0.7634 (Rank 1)
0.8314
(1.872)
0.77
0.82

[Back]

[Back]

Video 1. Cryo-EM structure of the PilQ complex.

Overall architecture of the PilQ complex. Sagittal sections along the central symmetry axis show the gates; top views of clipped planes indicate the detailed arrangement of all domains and both gates in PilQ from top to bottom.

[Back]

[Back]

Video 2. Atomic model of PilQ protomer.

Sagittal and longitudinal sections through the single PilQ domains show the quality of the fit of the homology models of the PilQ protomers. The overall arrangement of one PilQ protomer (chain B, cyan) after flexible fitting and linkers built and the adjacent protomers corresponding to chains A and C (blue and green respectively) show how the they are positioned relative to each other.

[Back]

[Back]

Figure 4. Comparison of different secretion systems with T. thermophilus PilQ.

Cryo-EM maps of seven different secretion systems aligned against the OM (grey bar) indicate two major classes. (A) E. coli R388 T4SS, S. enterica serovar typhimurium T3SS (consisting of proteins InvG, PrgK and PrgH), E. coli K12 T2SS and P. aeruginosa PilQ do not have any density outside the OM (grey bar). (B) T. thermophilus PilQ, N. meningitidis PilQ and V. cholerae GspD display a conical density outside the OM. Red arrows indicate poorly resolved regions of the N-terminal domains. (C) The structure of the lip domain (-αβ-βα-; cyan) is resolved in V. cholerae GspD (left) and forms the feature outside the OM in (B), but it is absent in PilQ from T. thermophilus (right). (D) Multiple sequence alignment of secretin domains from well-characterized complexes indicates that the lip region (cyan) and the C-terminal helical segments (olive) are absent in all PilQ homologues. We conclude that the features outside the OM (B) are not formed by the PilQ protomer.

Figure 4—figure supplement 1.View larger versionFigure 4—figure supplement 1. Comparison of secretin domains of different secretion systems.Cross sections through the cryo-EM maps indicate structurally conserved features in the secretin domain. All maps are low-pass filtered to the same nominal resolution (7 Å).

[Back]

[Back]

Figure 4—figure supplement 1. Comparison of secretin domains of different secretion systems.

Cross sections through the cryo-EM maps indicate structurally conserved features in the secretin domain. All maps are low-pass filtered to the same nominal resolution (7 Å).

[Back]

[Back]

Figure 5. PilQ mutants K728A and R730A destabilize the crown module.

(A) EM map of PilQ complex with fitted model (rainbow-colored from N to C terminus). Residues R730 and K728 (dark blue) are located just below the detergent belt (grey bar), close to the TM region of PilQ in the modeled homo-oligomer. Cross section through the cap region showing the crown module (red) above the OM surface. (B and C) Class averages of negatively stained PilQ-K728A/R730A (B) and wt-PilQ (C) obtained by averaging 374 and 3730 particles, respectively. The crown is absent in PilQ-KR, which shows weak, fuzzy densities instead, indicating that the crown is disordered.

[Back]

[Back]

Figure 6. Stability of the PilQ-KR variants and twitching motility assay.

(A, B) Stability of purified PilQ wildtype and PilQ-KR as indicated by SDS-PAGE (A) and Western blot analysis using PilQ antibodies (1:13000) (B) of membrane extracts boiled for 30 min in SDS sample buffer. (C) Contrast-inverted images of culture plates to quantify the twitching motility of the T. thermophilus strains (Scale bar = 5 mm).

[Back]

[Back]

Figure 7. Comparison of T. thermophilus T4PS complex in situ to detergent-solubilized PilQ.

Cross section of sub-tomogram averages of piliated (light blue) and unpiliated T4PS (light yellow) (Gold et al., 2015) with superposed 7 Å cryo-EM structure of PilQ complex.

[Back]

[Back]

Author response image 1.10.7554/eLife.30483.024Author response image 1.

[Back]

[Back]

Author response image 2. Assessment of PilQ symmetry and oligomeric state.

(A) 2D class average derived from 27 top views was not sufficient to assess the symmetry unambiguously. (B) Binary annular mask applied for multivariate statistical analysis. (C) Multiplication of (A) by (B) shows that the applied mask is not too lose and it just includes signal from the ring densities. (D) Second eigenimage after multivariate statistical analysis.

[Back]