An asymmetric sheath controls flagellar supercoiling and motility in the leptospira spirochete
Abstract
Spirochete bacteria, including important pathogens, exhibit a distinctive means of swimming via undulations of the entire cell. Motility is powered by the rotation of supercoiled 'endoflagella' that wrap around the cell body, confined within the periplasmic space. To investigate the structural basis of flagellar supercoiling, which is critical for motility, we determined the structure of native flagellar filaments from the spirochete Leptospira by integrating high-resolution cryo-electron tomography and X-ray crystallography. We show that these filaments are coated by a highly asymmetric, multi-component sheath layer, contrasting with flagellin-only homopolymers previously observed in exoflagellated bacteria. Distinct sheath proteins localize to the filament inner and outer curvatures to define the supercoiling geometry, explaining a key functional attribute of this spirochete flagellum.
Introduction
The spirochetes are gram-negative bacteria with helically coiled cells, and encompass a diverse phylum that includes the agents of leptospirosis (Leptospira spp.), syphilis (Treponema pallidum) and Lyme disease (Borrelia burgdorferi). Spirochetes have powerful swimming capabilities that enable them to rapidly disseminate through connective tissue, blood, and organs (Wunder et al., 2016b). This swimming capability is enabled by a unique configuration of the flagellum in these organisms (Charon et al., 2012; Li et al., 2000a; Wolgemuth et al., 2006; Figure 1A), which allows them to ‘drill’ through highly viscous media (Li et al., 2000b; Wolgemuth et al., 2006). By housing their flagella in the periplasm, spirochetes protect these filaments from environmental insults, such as immune surveillance (Vernel-Pauillac and Werts, 2018) (a similar function may be achieved by the membranous sheath that coats the exoflagellum in Vibrio cholerae Yoon and Mekalanos, 2008). Spirochetal flagella exhibit supercoiled architectures (Charon et al., 1992; Dombrowski et al., 2009), which have been linked to productive motility (Wolgemuth, 2015), and amongst pathogenic species, to virulence (Ko et al., 2009; McBride et al., 2005; Sultan et al., 2013; Wunder et al., 2016a).
Leptospira spp. includes pathogenic species that cause human and animal disease, as well as free-living saprophytic species (Picardeau, 2017). Unique among known bacteria, including spirochetes, Leptospira possess flagellar filaments that spontaneously curl into tightly supercoiled shapes when isolated, following nearly coplanar paths (Bromley and Charon, 1979; Figure 1—figure supplement 1). The flagellar filament superposes a curved path on the underlying helical body, resulting in overall ‘spiral’ or ‘hook’ cell end shapes depending on the motor rotation direction (Charon and Goldstein, 2002; Figure 1A). Similar to other spirochetes, Leptospira cells travel rapidly and unidirectionally when flagellar motors at opposing ends of the organism rotate in opposite directions (i.e. counterclockwise vs. clockwise when viewing the motor from the cell exterior), but stall when the motors are both rotating in the same direction (Goldstein and Charon, 1988; Wolgemuth, 2015). Leptospira can also adopt a distinct ‘crawling’ mode while bound to surfaces (Tahara et al., 2018).
As with all spirochetes, Leptospira have a complex flagellar composition compared to exoflagellated bacteria (Charon and Goldstein, 2002; Figure 1B; Table 1; Figure 1—figure supplement 2). Flagellar filaments from non-spirochete bacteria generally have a single protein component (Erhardt et al., 2010) (flagellin). In contrast, spirochetal flagellar filaments are comprised of a flagellin homolog (FlaB) and FlaA (Brahamsha and Greenberg, 1988), a second and completely unrelated component. FlaA and FlaB proteins can occur in multiple isoforms within a single organism (Wolgemuth et al., 2006; Table 1). Deletion of FlaA or FlaB in Leptospira leads to distinct phenotypes: deletion of FlaA affects flagellar curvature and diminishes motility and pathogenicity (Lambert et al., 2012), while deletion of FlaB eliminates the filament entirely (Picardeau et al., 2001). Similar observations have been made in other spirochetes (Li et al., 2000a; Motaleb et al., 2000), indicating that structural and functional adaptations in Leptospira flagella may be relevant to spirochetes as a whole. However, Leptospira flagella have unique features not found in other spirochetes. Recently, two novel Leptospira-specific components of the flagellum, FcpA and FcpB (Table 1), were identified whose deletion abolishes the ability of flagella to assume their characteristic supercoiled form and dramatically reduces motility and virulence (Wunder et al., 2016a; Wunder et al., 2018).
To date, the structure and composition of any spirochetal endoflagellar filament has remained elusive. Here, we present structures of flagellar filaments from the saprophytic spirochete Leptospira biflexa, solved by a combination of cryo-tomography, sub-tomogram averaging, X-ray crystallography and molecular docking. Our structures reveal that a conserved FlaB core assembly is enclosed by an asymmetric assembly of sheath proteins with novel folds and functions. Analysis of mutant filaments, which have lost one or more of these sheath elements, indicates that the sheath enforces different core lattice geometries on different sides of the filament, thus promoting a supercoiled flagellar shape. This asymmetric supercoiling mechanism may be relevant to other spirochetes, contributing to their distinctive swimming modes.
Results
FcpA and FcpB X-ray structures
Using X-ray crystallography, we established that the flagellar sheath components FcpA and FcpB adopt unique folds, unrelated to globular domains found in available flagellin structures (Figure 1C–F; Galkin et al., 2008; Samatey et al., 2001; Wang et al., 2017). FcpA from L. biflexa was crystallized in three different space groups (San Martin et al., 2017; Table 2), revealing a 3D fold previously unseen within the PDB (release July 2019) as reported by several structural alignment algorithms (Holm and Laakso, 2016; Madej et al., 2014). FcpB from the pathogenic species, L. interrogans (ortholog of the L. biflexa protein; Table 1; Figure 1—figure supplement 2), was crystallized in the orthorhombic space group P212121 (Table 2). The resulting fold exhibits distant structural homology to proteins of the Mam33 family, including human mitochondrial protein p32 (Jiang et al., 1999) that exhibit affinity for a number of protein partners. Homology to FcpB was also detected with VC1805 (Sheikh et al., 2008), a protein from Vibrio cholerae that can bind complement protein C1q.
FcpA and FcpB are highly distinct in overall size and shape. FcpA assumes an all-helical fold consisting of a twisted stack of eleven α-helices, arranged in a ‘V’-shaped architecture with arms of unequal length (Figure 1C). The long arm is formed by an α-helical hairpin structure (α3–4) emanating from a 4-helix bundle (comprising helices α2-3-4-6), while the shorter arm is mainly composed of two 3-helix bundles (α6-7-9 and α9-10-11) bridging to the longer arm via a helical hairpin (α9–10). In contrast to FcpA, FcpB adopts a wide and flat oblong shape, consisting of an 8-stranded anti-parallel β-sheet with two long α-helices (α1 and α4) packed on one side (Figure 1D). These major structural differences between FcpA and FcpB are evident even at relatively low resolution (15–20 Å), facilitating their identification in the tomographic reconstructions described below.
3D reconstruction of flagellar filaments
The strongly curved supercoiling geometry of intact Leptospira flagellar filaments (Figure 2A) precluded the use of most structure determination methods, including cryo-EM single-particle 3D reconstruction (Figure 1—figure supplement 1, Methods), motivating us to apply a cryo-electron tomography approach. Overlapping series of cubic sub-volumes following curved filament paths were manually selected from 62 reconstructed tomographic volumes from samples of purified, wild-type L. biflexa flagellar filaments. A total of ~10,800 sub-volumes were input into a customized sub-tomogram alignment procedure (Figure 2—figure supplements 1 and 2, Materials and methods). The resulting averaged map achieved an overall (average) resolution of ~10 Å (Figure 2B–E; Figure 2—figure supplement 2) and exhibits a series of concentric protein layers centered around a ~ 2 nm central pore (Figure 2F), as well as a distinct array of globular features decorating the filament surface (Figure 2G). The morphology of the inner two density layers (Figure 2G) closely resembles that of the helically assembled flagellin domains D0 and D1 of filaments from exoflagellated bacteria (Wang et al., 2017; Yonekura et al., 2003). This structural homology is consistent with earlier predictions that the core of the spirochete filament is composed of FlaB (Li et al., 2008; Nauman et al., 1969), a protein that lacks the D2 and D3 domains present in flagellins from many exoflagellated bacteria (Figure 3—figure supplement 1).
Asymmetric sheath composition
The tomographic map identified sheath density features that surrounded the two core layers but were not uniformly distributed. Sheath densities were mainly localized on the outer curvature of the filament, or ‘convex’ side, thus offsetting the filament center of mass away from the central pore when the filament is viewed in cross section (Figure 2C,F, Video 1). This break in helical subunit arrangement in the Leptospira endoflagellum sharply contrasts with symmetric structures which were previously observed for exoflagella from Salmonella (Sasaki et al., 2018; Wang et al., 2017; Yonekura et al., 2003), Campylobacter (Galkin et al., 2008), Pseudomonas aeruginosa and Bacillus subtilis (Wang et al., 2017).
We used two independent docking strategies to fit FcpA and FcpB molecular envelopes into the reconstructed filament map (see Materials and methods) and identified convergent locations for these two proteins that were distributed across the convex side of the sheath region (Figure 3, Video 1). These sheath sites form an array of globular features that closely follows the canonical 11-protofilament helical lattice identified in flagellar filaments of Salmonella and other non-spirochete bacteria (Figure 3—figure supplement 2). For FcpA, six adjacent, curved rows of ‘V’-shaped density profiles form an interlocking array (Figure 4; Figure 3—figure supplement 3). The FcpB sites were resolved as four adjacent rows of flattened density lobes protruding slightly outward from the FcpA portion of the sheath (Figures 3 and 4). The resulting assembly consists of an inner layer of FcpA molecules in close contact with the FlaB core, and an outer layer of FcpB molecules nestled between rows of FcpA (Video 2).
The asymmetric localization of FcpA and FcpB, specifically on the convex side of the sheath, is distinct from previously characterized exoflagella but is consistent with prior flagellum electron-microscopic imaging studies performed with antibody labels (Wunder et al., 2016a; Wunder et al., 2018). The localization of FcpB sites were further validated by comparison with a similar reconstruction of filaments from an fcpB- mutant (Figure 2I), for which the projected cross-section of the fcpB- map was missing four globular features found in the wild-type map (Figure 2H) that coincide with our FcpB docking solutions. The signal-to-noise ratio was too low to perform similar sub-tomogram averaging for filaments from an fcpA- mutant. However, the filaments from imaged fcpA- mutants were significantly smaller in diameter than filaments from wild-type and fcpB- mutants (see the last Results section, below), indicating that filaments from the fcpA- mutant have lost FcpA and FcpB, in agreement with previous findings (Wunder et al., 2018).
Additional sheath features are found on the filament concave side, highly divergent from the ones observed on the convex face. These inner-side elements manifested as two rows of elongated density units, separated by a groove that exposes the core (Figures 2C and 3A; Figure 3—figure supplement 4). Further modeling of the sheath was precluded as atomic models are not available for FlaA1 and FlaA2, the only other known components of the Leptospira flagellar filament. FlaA is also presumed to reside in the spirochete sheath (Li et al., 2008), and comparing the fcpB– map with the wild-type map, a missing mass within the concave sheath region indeed suggests candidate positions for FlaA1 and FlaA2 (Figure 2I, red asterisk).
The current data do not unambiguously rule out the possibility that FcpA and/or FcpB may occupy one or more sites on the inner curvature. However, a quantitative comparison of core-sheath contact sites for all 11 sites related by pseudo-helical symmetry reveals exactly six protofilaments that share the distinctive ‘V’-shaped footprint defined by our FcpA docking studies, while features of corresponding, pseudosymmetry-related sites diverge sharply in the remaining 5 protofilaments (Figure 3—figure supplement 3B). Thus, if any of the latter 5 sites are occupied by FcpA and/or FcpB, they must adopt distinctly different binding interfaces with the core. A comprehensive accounting of these features within the concave sheath region awaits further structural and biochemical data.
Filament core structure
Additional averaging of the map core region revealed homology with flagellar filaments from Salmonella and other exoflagellates (Wang et al., 2017). Due to a strongly preferred filament orientation in the specimen ice layer (Figure 1—figure supplement 1), the potential effect of a ‘missing wedge’ of Fourier-space data characteristic of the cryo-tomography (Carazo et al., 2006) could not be fully eliminated. The resulting map resolution anisotropy (Figure 2—figure supplement 3) interfered with structural analyses of the core region, where globular features were neither evident nor expected based on available models of the flagellar core (Yonekura et al., 2003). By performing 3D auto-correlation analysis within the core region, however, we detected a signal corresponding to the 11-fold helical symmetry operator present in previously determined flagellar structures (Figure 3—figure supplement 5). The 11-fold symmetry operator is mostly congruent with the observed lattice positions of FcpA and FcpB molecules identified in the sheath (Figure 3—figure supplement 2), although there are minor deviations due to filament curvature. After applying 11-fold averaging using the symmetry operator, additional features were resolved in the core region, including some α-helices, and revealed a structural homology with the D0/D1 region of the Salmonella flagellum (Figure 3B,C; Figure 3—figure supplement 6; Video 2).
Following 11-fold averaging, improved resolution in the core region allowed for straightforward fitting of a core atomic model representing a single 52 Å repeat of FlaB (11 subunits), which we obtained by applying estimated pseudo-helical symmetry parameters to a FlaB monomer homology model (backbone-only atoms) derived from the reported flagellin atomic structure (Wang et al., 2017; Yamashita et al., 1998; Yonekura et al., 2003) (see Materials and methods). The resulting model also shows good agreement with the non-symmetrized map when the superposition is viewed from favorable orientations that minimize missing wedge artefacts (Figure 3B,C; Video 3).
Conserved protein-protein interfaces
The resulting filament atomic model exhibits a consistent set of interactions between FcpA and the FlaB core, reflecting the approximate helical symmetry of fitted components (Figure 4A–C). FcpA molecules project both arms toward the core region, positioning the long arm within contact radius of two predicted loop regions in FlaB (Figure 4D). Such interactions bury the same protein surface of the FcpA α3-α4 hairpin that forms crystal lattice contacts in all three FcpA X-ray structures. The FlaB loop regions engaged in these interactions project outwards from the core surface toward the predicted FcpA arm positions. Thus, FlaB loops in Leptospira (both 6–7 residues in length) likely contribute to the functional interface between FcpA and the core.
Modeled FcpB subunits project their central β-sheets radially away from the core and present several elements within contact radius of the FcpA inner sheath layer, including one of the long FcpB α-helices (α4) as well as several loops located at the proximal edge of the β-sheet (Figure 3D). We identified two highly conserved surface patches on FcpB which are in close proximity with nearby FcpA molecules (Figure 4E), suggesting that these sites may mediate direct interactions between FcpA and FcpB. The corresponding FcpA surface interaction regions predicted by our model are also conserved, and largely match a conserved protein:protein interface found in all FcpA crystal forms. Complementary electrostatic potentials also map onto the FcpA:FcpB interface (Figure 1E,F). The modeled FcpB sites are too far away from the core to make direct contact, providing an explanation for why deletion of FcpA leads to concurrent loss of FcpB from the filament (Wunder et al., 2018).
The wild-type filament map provides evidence of contact between FcpA and unassigned elements within the concave sheath region, particularly on one side of the groove (Figure 3A; bottom left corner). Since the unassigned material in the concave sheath region is likely composed of FlaA1 and/or FlaA2, FcpA may interact directly with one or both FlaA isoforms. Indeed, pull-down assays of purified L. biflexa flagellar filaments recently identified interactions between FcpA and FlaA2 (Sasaki et al., 2018).
The sheath promotes filament curvature
Trajectory analysis of aligned filament subtomograms from wild-type, fcpA– and fcpB– samples reveals a progressive straightening in isolated flagellar filaments as sheath components are lost (Figure 5A–C). This observation supports and extends previous reports associating these mutations with loss of filament curvature together with pronounced functional effects (Wunder et al., 2016a; Wunder et al., 2018). Strikingly, straightening observed in the mutant filaments manifests as population shifts between multiple peaks in the observed curvature distributions. Wild-type filament segments almost exclusively fall within a single peak in the measured curvature distribution, centered at the average value (~5 µm−1; 0.2 µm radius of curvature; Figure 5A). In contrast, mutant fcpB– filament segments (which still contain FcpA Wunder et al., 2018) show a mixed population of curvatures, with only ~50% of segments falling within the wild-type peak, while curvature for most of the remaining segments is estimated at 3 µm−1 or less (Figure 5B). Curvature is further reduced in fcpA– filaments, which lack both FcpA and FcpB (Wunder et al., 2018): only a vestigial population of segments exhibit the wild-type curvature, while the estimated curvature of most segments is again below 3 µm−1 (Figure 5C).
The fcpA– filament segments can be divided into two populations, revealing a further correlation between sheath structure and filament curvature. A ‘thick’ fcpA– segment population, whose estimated diameter (~170 Å) was consistent with the presence of a partial sheath, shows a propensity to form segments with curvature of ~3 µm−1 (Figure 5C, blue bars). A comparable number of ‘thin’ fcpA– segments, with an estimated diameter (~120 Å) indicative of an unsheathed core, shows a population shift towards lower curvature (Figure 5C, orange bars). Segment 2D class averages of fcpA– filaments support these observations, revealing straighter, ‘thin’ classes indicative of a ‘bare’ FlaB core, as well as thicker, more curved classes with a visible, partial sheath that may preferentially localize to the inner curvature (Figure 5—figure supplement 1). Thus, structure comparisons of wild-type and mutant filaments establish a trend where (1) loss of FcpB converts a substantial fraction of filament segments from the flat, tightly supercoiled, wild-type form to a substantially less-curved form; (2) loss of FcpA and FcpB converts nearly all filaments to less-curved forms; and (3) loss of additional sheath components further straightens the less-curved forms, on average. We occasionally observed evidence of the same trend while imaging individual filaments, where an abrupt loss of curvature in a filament coincides with a reduction in filament diameter (indicative of sheath ‘shedding’; Figure 5—figure supplements 1 and 2).
Discussion
Here we have identified a novel, asymmetric subunit arrangement in the Leptospira flagellum, previously unobserved either in exoflagellated bacteria or in other spirochetes. We have further demonstrated that this asymmetric architecture is functionally linked to flagellar supercoiling, a key attribute associated with spirochete motility. This means of controlling supercoiling contrasts with factors previously identified in other bacteria, which include changes in buffer composition (salt or pH) or the introduction of torsional strain due to flagellar rotation (Calladine et al., 2013). However, earlier investigations of flagellar filament supercoiling have focused on exoflagellates, whose filaments are more simply constructed (homopolymers vs. the heteropolymer considered here), and whose functional roles differ in important ways from spirochete flagellar filaments.
An important innovation in the current work was the application of a subtomogram averaging procedure to solve structures of the supercoiled form present in wild-type Leptospira flagellar filaments. This approach contrasts with previous structural studies of flagellar filaments, most of which utilized mutant forms that strongly favor a straight, helically symmetric filament conformation, as was required for high-resolution cryo-EM analysis. One prior study reports the structure of a supercoiled spirochete flagellum (Liu et al., 2010); however, in that case the filament was analyzed under the assumption of helical symmetry and was computationally straightened, obscuring any asymmetry. More recently, a supercoiled form of the Salmonella flagellar hook was resolved to near-atomic resolution by single-particle cryo-EM (Kato et al., 2019; Shibata et al., 2019). However, the methods used in these latter works are not well suited for all supercoil morphologies, particularly large-diameter supercoils (as in the current work) that are unfavorably oriented within the specimen ice layer.
The flattened supercoil shape observed in our wild-type Leptospira sample, while consistent with prior observations of these filaments (Wolgemuth et al., 2006), diverges from forms seen in other bacteria, including exoflagellates (Calladine, 1975; Calladine, 1976; Calladine et al., 2013; Fujii et al., 2008; Kamiya et al., 1980; Wolgemuth et al., 2006). Functional exoflagella, including from Salmonella, interconvert between up to 12 distinct supercoil polymorphs, but preferentially adopt an extended, left-handed supercoil form at rest or during smooth swimming (‘Normal’; superhelical pitch and diameter of 2.66 and 0.6 µm, respectively). The left-handed, low-pitch geometry we identify in wild-type Leptospira filaments (Figure 5—figure supplement 3) more closely resembles an alternative exoflagella supercoil type that forms under low pH conditions (‘Coiled’; superhelical pitch and diameter of ~0 and 1.0 µm, respectively). However, the supercoil diameter in the Leptospira filament (0.45 µm) is less than half that of the ‘Coiled’ exoflagellar form. Overall, wild-type Leptospira filaments are much more curved (curvature is ~5 µm−1; Figure 5A) than observed or predicted polymorphic forms from exoflagellates, which range in curvature from 0 ~ 3 µm−1 (Calladine, 1975; Kamiya and Asakura, 1976).
Loss of sheath proteins reduces the curvature of Leptospira filaments, resulting in supercoiling geometries more consistent with the polymorphic forms observed in exoflagella. Despite their loss of curvature compared with wild-type filaments (Figure 5A,C; Wunder et al., 2016a), fcpA– filaments remain markedly more curved than cryo-EM samples of tobacco mosaic virus (TMV) filaments or amyloid fibrils, which lack inherent supercoiling character. For ‘straight’ TMV and amyloid filaments, an average curvature ~0.7 µm−1 was estimated from cryo-EM specimens (attributable to thermal fluctuations, variable fluid flow during sample preparation and other factors) (Rohou and Grigorieff, 2014), compared to values of up to ~3 µm−1 observed in the fcpA– sample (Figure 5C). Despite loss of most or all of the sheath, the ‘thin’ fcpA– population maintains an average curvature value (~2 µm−1) that exceeds that of straight TMV or amyloid filaments by more than a factor of two (Figure 5C). The estimated curvature for these bare Leptospira filament cores, falls within the range observed in exoflagellar polymorphs (0–3 µm−1) (Calladine, 1975). This characteristic, as well as the occasional observation of ‘wavy’, sinusoidal fcpA– filaments of appropriate shape parameters (Figure 5—figure supplement 4), supports the idea that the Leptospira flagellar filament core may retain polymorphic, supercoiled behavior similar to exoflagellates.
Regardless of specific structural and functional core properties, our data indicate that binding of Leptospira sheath proteins (particularly FcpA and FcpB) initiates an allosteric, cooperative transition of the core from a straighter (but potentially polymorphic) form (Figure 5D,E), to the flattened, highly curved supercoil geometry suited for swimming (Figure 5F,G). A potential mechanism involves recognition by the sheath proteins of distinct core lattice geometries that occur on different sides of a supercoiled filament. For Leptospira, these lattice differences are quite substantial: for the given filament geometry (120 Å core diameter and ~0.45 µm superhelix diameter), the axial repeating distance is 3 Å longer on the filament outer curvature than the inner curvature. Therefore, if FcpA preferentially binds to longitudinally expanded protofilaments on the outer curvature of a supercoil, and the core preferentially adopts a shorter longitudinal repeat distance, a helical symmetry mismatch will be generated and a specific supercoiled geometry can be stabilized (Figure 5F,G; Caspar, 1963; Klug, 1967). In this type of mechanism, binding of an additional sheath component (e.g., FcpB) that is specifically matched to the expanded FcpA lattice would be ideally suited to further reinforce the supercoiling geometry, as is observed here (Figure 5A–B).
A noteworthy feature of the flagellar supercoiling mechanism presented in Figure 5 is that the asymmetric sheath would lock the core into a single, nondegenerate shape. In contrast, a supercoiled, homopolymer flagellar filament can potentially interconvert between 11 equivalent conformations with identical shape parameters but different, rotated patterns of expansion and contraction in the 11 protofilaments. Thus, a supercoiled homopolymer flagellum could propagate a helical wave without rotating (Klug, 1967) or equivalently, rotate without propagating a helical wave; an asymmetrically assembled heteropolymer flagellum as we describe for Leptospira would be unable to do so.
Unique attributes of Leptospira endoflagellar filaments, compared with exoflagella, could be relevant to the different mechanical requirements of these two, divergent propulsion systems. Physics modeling studies indicate that curved filament morphology and stiffness are critical for motility in Leptospira and spirochetes in general (Berg et al., 1978; Dombrowski et al., 2009; Kan and Wolgemuth, 2007). Both of these properties can be enforced via the addition of an asymmetric sheath surrounding the filament core. To support the large loads involved in driving whole-body undulations, spirochete endoflagella possess specialized reinforcements, including extra-large motors capable of exerting higher torques than conventional exoflagellar motors (Beeby et al., 2016) and unique hooks exhibiting covalent crosslinking between FlgE subunits (Miller et al., 2016). Moreover, Leptospira, alone among the known spirochetes, has only a single flagellum at each cell end, whereas Borrelia, Treponema, and other species exhibit multiple copies per cell end (Charon et al., 2012). One function of the sheath elements FcpA and FcpB, to date found only in Leptospira, may therefore be to reinforce the endoflagellum to allow for even higher torque loads required in the absence of additional, supporting flagella (Trachtenberg et al., 1986; Trachtenberg et al., 1987).
Although FcpA and FcpB play conspicuous roles in supercoiling, our images of fcpA– filaments indicate that the remaining sheath components are themselves capable of asymmetric binding, and may do so preferentially on the concave side of the filament (Figure 3; Figure 5—figure supplement 1; Table 3). These sheath proteins include FlaA1 and FlaA2, although we are not currently able to localize these in the filament map. Orthologs of Leptospira FlaA proteins are found in all spirochetes (Table 1). Furthermore, loss of FlaA has been shown to modify the diameter (Li et al., 2000a) and superhelical pitch (Lambert et al., 2012) of spirochetal filaments. These observations raise the possibility that asymmetric arrangement of components in the sheath (and/or the core) may be a general feature that contributes to flagellar supercoiling and motility in spirochetes.
Materials and methods
Strains and cell culturing
Request a detailed protocolLeptospira biflexa serovar Patoc strain Patoc I (Paris) wild type, fcpA- and fcpB- mutant cells were cultured in Ellinghausen–McCullough–Johnson–Harris (EMJH) liquid medium until they reached logarithmic phase at 30°C (Wunder et al., 2016a; Wunder et al., 2018).
Periplasmic flagella purification
Request a detailed protocolPurification of periplasmic flagella was performed as described (Wunder et al., 2016a; Wunder et al., 2018).
Recombinant protein crystallization and crystal structure determination
Request a detailed protocolFcpA (Uniprot B0STJ8) from L. biflexa serovar Patoc strain Patoc I (Paris) was expressed, purified and crystallized as previously reported (San Martin et al., 2017). Three different crystal forms were obtained (Table 2), and the structure was solved ab initio with the hexagonal data set using Arcimboldo (Rodríguez et al., 2009). Resulting electron density maps were automatically traceable (Lamzin et al., 2012), and the model was subsequently refined with Buster (Bricogne G. and Bricogne et al., 2017) iterated with manual rebuilding and validation with Coot (Emsley et al., 2010). The other two monoclinic crystal forms were solved by molecular replacement (McCoy et al., 2007) using the hexagonal structure as a search probe, then further refined following similar procedures. The final refined model spans residues 55 to 291, while the first 54 residues toward the N-terminus were disordered.
The gene LIC11848 (Uniprot Q72RA0) coding for FcpB from L. interrogans serovar Copenhageni strain Fiocruz L1-130, was cloned in expression plasmid pQE80 (Qiagen), including a TEV-cleavable 6xHis tag. Overexpression was achieved in E. coli Rosetta-gami 2 (DE3), induction conditions, chromatographic purification and TEV cleavage procedures were performed as for FcpA (San Martin et al., 2017). After cleavage, FcpB started at sequence SQQNSGS (thus excluding the signal peptide), with an N-terminal overhanging tripeptide GSG derived from the plasmid. Single FcpB crystals were grown at 20°C using vapor diffusion technique in hanging drops, mixing 2 µl protein solution (18 mg/mL in 20 mM Tris.HCl pH 8.0, 150 mM NaCl) with 2 µl reservoir solution 0.4M NH4I, 26% (w/v) PEG 3350, 0.05M MES pH 6.5, 21% glycerol. Initial anomalous signal was detectable only to 4.5 Å, improved signal was achieved by quick-soaking crystals in 0.1M ‘magic triangle’ I3C (Beck et al., 2008) for 20 s. Highly redundant data sets were processed using XDS (Kabsch, 2010) and scaled with Aimless (Evans, 2011). The structure was solved by SAD, exploiting the iodine atoms as anomalous scatterers, the iodine sub-structure was solved with ShelxD (Schneider and Sheldrick, 2002) and refined with Sharp (Bricogne et al., 2003). Density modification was performed with Pirate and resulting electron density maps allowed for manual chain tracing and initial model refinement with Buster. A second dataset diffracting X-rays to higher resolution was eventually used for final refinement with Buster, iterated with manual rebuilding and validation with Coot.
Structural analyses were done with the CCP4 suite (Winn et al., 2011) of programs and with Pymol (Schrodinger LLC, 2015).
Cryo-electron microscopy
Request a detailed protocolL. biflexa wild-type flagella. Purified L. biflexa wild-type flagella were mixed with 6x concentrated Protein A conjugated with 10 nm colloidal Au (EMS Aurion, Hatfield, PA). Prior to sample application, the 300 mesh Cu Quantifoil grids (Ted Pella, Inc, Redding, CA) with 1.2 µm diameter holes were plasma discharged with H2O2 for 30 s in a Model 950, Solarus Advanced Plasma System (GATAN). Approximately 3 uL of 2:1 (v/v) flagella: 6x Protein A Gold was applied to each grid. The grids were plunge frozen into liquid ethane in a Mark III Vitrobot (FEI Company, Eindhoven, The Netherlands) following a 2 min incubation time, blot time of 6–7.5 s, and blot offset of −2 mm at 18°C and 100% humidity.
Tilt series were acquired on a 300 kV Polara cryo-TEM (UT Houston, Houston, TX) using SerialEM (Mastronarde, 2005) (University of Colorado, Boulder, CO) equipped with a K2 direct electron detector (GATAN, Inc, Pleasanton, CA). Tilt series were acquired with a defocus of between −2.0 and −4.0 µm at +/- 53° with 3° increments in counting mode at a dose rate of ~8 e-/pixel/second. At each tilt angle, 12 frames were collected with 0.1 s exposure per frame. The total dosage for each tilt series was ~60 e-/Å2. The pixel size at a magnification of 15500x was calculated to be 2.604 Å at a binning of 1. A total of 120 tilt series were acquired.
L. biflexa fcpA- and fcpB- flagella. Purified fcpA- and, fcpB- flagella were prepared similarly to the wild-type flagella and vitrified as described above. Mutant fcpA- or, separately, fcpB- flagella were mixed with 6x concentrated 10 nm colloidal Au BSA Gold Tracer beads (EMS Aurion) prior to plunge freezing. The 200 mesh Cu C-flat (EMS) grids with 1.2 μm hole spacing had a 3 nm thick layer of carbon applied to ensure distribution of gold beads and flagella. The grids were plasma discharged with a 75% argon/25% oxygen mixture for six seconds in a Solarus Advanced Plasma System, Model 950 (GATAN).Tilt series were then acquired on a 200 kV Tecnai F20 cryo-TEM (CCMI, Yale University, New Haven, CT) using SerialEM (Mastronarde, 2005) (Boulder, CO) and imaged with a K2 direct electron detector (GATAN) in counting mode. The Tilt series were acquired over +/- 60° in 3° increments with a defocus of −2.5 μm. At a magnification of 14500x at spot size 7, the dose rate was ~8 e-/pixel/second. The exposure time per tilt angle was 1.5 s with ~2.25 e-/Å2 and was dose fractionated into 7 frames, each receiving 0.2 s of exposure time. The final dose was ~90 e-/Å2 with a pixel size was 2.49 Å. A total of 26 tilt series were acquired.
Tilt series reconstruction and subtomogram averaging
Request a detailed protocolL. biflexa wild-type flagella. The frames in each tilt angle in a tilt series were motion corrected using IMOD alignframes (Mastronarde and Held, 2017). The tilt series were then aligned in IMOD, version 4.9.4_RHEL6-64_CUDA8.0 (Mastronarde and Held, 2017) using local alignment with the fiduciary markers. After alignment, tilt series underwent CTF-estimation and correction using phase flipping in IMOD, followed by gold subtraction (Himes and Zhang, 2018; Mastronarde and Held, 2017). Final CTF-corrected and aligned tilt series binned by 2, resulting in a pixel size of 5.208 Å and then reconstructed in Tomo3D (Agulleiro and Fernandez, 2011) using weighted back projection. IMOD 3dmod was used to trace a path of particle points along the center of each filament. The addModPts program in PEET (Cope et al., 2011) was used to set the repeat spacing between each particle point to 52 Å. Using the IMOD model2point program, the particle model file was converted to a text file format that designated X, Y, and Z coordinates for each filament and assigned an incrementing number for each individual filament for import into RELION version 2.1.b1-gcccuda-2016.10-cc37 (Scheres, 2012).
The RELION (Scheres, 2012) sub-tomogram averaging module (Bharat et al., 2015) was used for initial sub-tomogram alignment. EmClarity (Himes and Zhang, 2018) was for subsequent high-resolution refinement and 3D reconstruction, and an in-house method for smoothing and interpolation of sub-volume coordinates (Huehn et al., 2018; Figure 2—figure supplement 2).
A variation of the RELION pre-processing PYTHON script (Bharat et al., 2015) was developed in house to sort particles by filament to ensure splitting of particles for gold standard FSC calculation. This ensured that two adjacent particles in a filament would not be randomly sorted into separate groups, thereby avoiding potential particle coordinate overlap if two or more adjacent particles occupied the same repeat within a filament during averaging and alignment. To ensure minimal polarity switching during particle alignment, larger cubed segments of filaments (500 Å) were extracted and CTF-corrected using CTFFind4. The resulting alignment parameters were imported into RELION version 2.1 beta for unsupervised sub-tomogram averaging using RELION 3D auto-refinement. No mask was used. For 3D auto-refinement, particle dimensions were set to 490 Å to include 490/52 = ~ 9.4 subunits along the filament, and the particle symmetry was designated as C1 due to the inherent asymmetry of the filament. Helical reconstruction options are not implemented for subtomogram averaging in this version of RELION, so the refinement was performed in single-particle mode. An initial angular search spacing of 15° was used. A local search threshold of 1.8° were used for auto-sampling. For local searches, the angular search range in psi and theta (tilt) were 10° and 15°, respectively. The final averaged volume was reported to be 18 Å after FSC calculation during RELION post-processing. Subtomogram averaging in RELION utilized a total of the best 62 tomograms with 10,851 particles.
The flagellum filament structure resulting from subvolume averaging in RELION had a final resolution of 18 Å, but evidently suffered from severe reconstruction artifacts. Efforts to use the RELION subvolume average as a search template for the emClarity (Himes and Zhang, 2018) particle picking procedure failed to yield useful particle coordinates. In-house scripts were therefore utilized to export the RELION alignment parameters to emClarity. Following coordinate import, CTF estimation and correction was performed in emClarity. The optimal defocus estimate was set to −3.0 µm and the defocus window was 0.5 µm.
For averaging and alignment steps in emClarity, particle dimensions and box sizes, as well as in-built shape masks were specified in the emClarity parameter file. Following initial estimates for a filament diameter of ~240 Å and a repeat spacing of 52 Å, the particle ‘radius’ was set to (145 Å, 145 Å, 26 Å) in the parameter file (param0.m). A rectangular alignment mask with x, y, z dimensions of 200, 200, and 160 Å was applied, encompassing 6–8 repeats. Anisotropic SSNR calculation was activated (flgCones = 1). To avoid FSC overestimation due to mixing of closely packed particles together between half datasets, the fscGoldSplitOnTomos option was also activated.
Iterative averaging and alignment steps were run in emClarity for 3–5 cycles with increasingly restrictive search angles and translational shifts. Particle coordinates and Euler angles from several cycles (1, 3 and 5) were extracted from the emClarity results file for analysis of individual filaments. Areas in the tomograms where filaments had lower SNR with respect to the background solvent had poorer tracking of particle coordinates along the filament trajectory, leading to ‘breaks’ and translational shifts in the x,y,z coordinates or Euler angle misalignments. For each filament, particle coordinates were extracted and graphed using in-house scripts. Particle outliers were deleted, and an in-house script was used to smooth a trajectory through the remaining coordinates, and to fill in gaps in the smoothed trajectory with particles according to a repeat spacing of 52 Å. The particle assessment, excision, and trajectory smoothing procedures were performed several times until the particle coordinates followed a regular, smoothed trajectory. These final coordinates were then re-imported into emClarity for a final 8 cycles of averaging and alignment at a binning of 2 for cycles 0 through 3, and a binning of 1 for cycles 4 through 8. The angular search range (rawAngleSearch=1st,2nd,3rd,4th,5th) can be set for each alignment cycle with the first two values corresponding to the angular range and increment of the out-of-plane search. The 3rd and 4th values correspond to the angular range and increment for the in-plane search. The 5th value was set to 0 to turn off the helical search, which was not effective for our system. At a binning of 2, the rawAngleSearch range for cycles zero through three, in order, were (0,0,7,1,0), (4,1,0,0,0), (0,0,3,0.75,0), and (2,1,0,0,0).For cycles four through eight, at a binning of 1, the rawAngleSearch ranges were (0,0,5,1,0), (3,1,0,0,0), (0,0,3,0.75,0), (2.25,0.75,0,0,0), and (0,0,1.5,0.5,0). The two half-dataset volumes were combined with a B-factor of 0, and the Gold standard FSC at 0.143 was calculated to be 9.83 A with anisotropic resolutions ranging from 8.89 to 16.17.
L. biflexa fcpA- and fcpB- flagella. A similar sub-tomogram averaging procedure as described above was followed for the fcpA- and fcpB- flagella tomograms. Specifically, IMOD was used for motion correction (using the program alignframes) and the tilt series were locally aligned through use of the fiduciary markers. CTF correction was also performed in IMOD through phase flipping, followed by subtraction of the fiduciary markers from the tomograms. The final aligned tomograms were binned by 2, resulting in a pixel size of 4.97 Å. Tomograms were then reconstructed in Tomo3D using weighted back projection, and 3dmod (IMOD) was used to select the filaments. The addModPts program in PEET was used to add points corresponding to the repeat spacing of 52 Å, and the IMOD program model2point was used to prepare the files for import into RELION.
The fcpB- tomograms were imported into RELION using the same scripts as described above, with a box size of 480 Å. A total of 24 tomograms were imported, corresponding to 5917 particles. In the program 3dautorefine, no initial reference model was used, and the symmetry was C1, the particle mask was 470 Å, and the same angular search parameters as for the wild-type structure were used. The resultant fcpB- structure had a resolution of 26.5 Å, and similar to the wild-type reconstruction, appeared to have large missing wedge artifacts. For fcpA-, a total of 25 tomograms were used, yielding 4327 particles; particles were imported into RELION and refined using similar parameters. Due to limited signal quality and/or sample heterogeneity, the fcpA- refinement did not yield a converged structure, with inconsistent particle axial rotations along single filaments. Nevertheless, the resulting fcpA- subunit x, y, z coordinates yielded continuous trajectories for many of the filaments. Using 2D classification, two distinct fcpA- particle diameters were identified, with 2958 particles in a smaller-diameter class (~120 Å) and 323 particles in a larger-diameter class (~170 Å). Particle coordinates from these two classes were used for subsequent trajectory/curvature analysis.
For the fcpB- dataset, the RELION-refined x, y, z coordinates and Euler angles of all 24 tomograms and 5917 particles were imported into emClarity, using the same in-house scripts as described above. CTF correction and estimation was performed, with an initial defocus estimate and window of −5.0 μm and 2.0 μm respectively. Due to reduced signal to noise ratio in the fcpB- data, the reference alignment procedure failed to unambiguously establish the polarity of many of the filaments, which were therefore discarded. After this analysis, 1163 particles corresponding to 14 filaments from 10 tomograms remained.
Initial emClarity parameters designated the particle radius in x, y, and z as 150, 150, and 20 Å, and applied a rectangular alignment mask of 180,180, 250 Å. As in the wild-type dataset, the flgCones and fscGoldSplitonTomos options were activated. Five cycles of averaging and alignment were used, in a manner similar to that described above, with each cycle using a binning of 2. The rawAngleSearch range for the cycles, in order, were (0,0,7,1,0), (5,1,0,0,0), (0,0,5,1,0), (4,1,0,0,0), (0,0,3,0.75,0), and (0.75,3,0,0,0). After the fifth cycle, the half data sets were combined, giving an average FSC of 18.41 Å, with a range from 14.58 Å to 31.67 Å at a B-factor of 0.
Resolution estimation
Request a detailed protocolInitial resolution estimates of final subtomogram average volumes were estimated using emClarity as described (Himes and Zhang, 2018; Figure 2—figure supplement 2). For additional resolution estimates, emClarity source code was modified to output unmasked copies of both refined half-dataset volumes, and these volumes were input to other software packages: local resolution estimates (Figure 2—figure supplement 4) were made using ‘Relion version 3.0.6 (Zivanov et al., 2018), and direction-dependent resolution estimates (Figure 2—figure supplement 5) were made using 3D FSC (Tan et al., 2017). For local resolution estimates, maps were sampled at a rate of 40 Å (Relion option ‘--locres_sampling 40’), with local regions defined by a spherical mask of radius 20 Å softened by a cosine edge function of (40 Å falloff).
Additional masks used in combination with direction-dependent resolution calculations were defined as follows. A curved, quasi-cylindrical mask enclosing the sheathed filament and following the filament path, with a circular cross-section (radius 240 Å), was defined using the ‘shape tube’ command from UCSF Chimera (Pettersen et al., 2004) to generate curved, tubular masking envelopes. To better match the filament cross-section, the intersection of three such masks, slightly offset from each other, was taken to produce a tighter mask with approximately elliptical cross-section. The filament core was defined using the Chimera ‘shape tube’ command with a radius of 120 Å. These masks were truncated axially to a length of three filament repeats (~106 Å) and then softened by convolution with a spherical cosine edge function (fall-off width of 31 Å).
Model building, fitting of protein components into tomographic maps and minimization of the L. biflexa flagellar filament
Request a detailed protocolModels of L. biflexa FlaB1 and FcpB were built by homology modeling with Rosetta63. 10,000 models of FlaB1 were thus generated using the structure of the flagellar filament of B. subtilis (PDB 5WJT) in a straight symmetric context. In the case of FcpB, we used our L. interrogans experimental structure (PDB 6NQZ) as a template, to generate 20,000 models.
Using SITUS (Kovacs et al., 2018; Wriggers, 2012), the averaged volume of the L. biflexa wild-type flagella generated by emClarity was converted to a SITUS readable file using the SITUS map2map program. Using the colores program, single monomers of FcpA and, separately, FcpB were rigidly fit into the ~10 Å resolution 3D volume with Laplacian filtering using an angular step of 5° where quasi-uniform angular spacing was enforced by the pole sparsing method. The highest-scoring fittings thus obtained gave alignments for two protofilaments of FcpA and one protofilament of FcpB within the outer convex sheath region. To identify additional docking locations within the sheath where resolution anisotropy (due to the missing cone of Fourier data) may have reduced the efficacy of the fitting procedure, we performed additional SITUS searches with the same parameters as before but with modified input maps and models. Specifically, we (i) focused searches on specific sheath regions by segmenting the cryo-EM map into a series of overlapping sub-regions and (ii) increased the signal power of the FcpA and FcpB search models by generating symmetry-related ‘triplexes’ of the initially discovered docking hits. Sub-region segmentation of the flagellar sheath in our 3D volume was performed using the UCSF-Chimera (Pettersen et al., 2004) module Segger. To generate a ‘triplex’ of FcpA or FcpB atomic coordinates, the highest-scoring initial hit was axially shifted by one 52 Å helical repeat forward and backwards to generate a composite PDB file containing three FcpA (or FcpB) monomers at consecutive sites along a single protofilament. Top-scoring hits from this second round of SITUS searching identified four additional protofilaments of FcpA, for a total of 6 protofilaments, and three additional protofilaments of FcpB for a total four protofilaments.
The relative positioning of FcpA and FcpB in our sheath model was further supported by an independent computational docking strategy, in which we utilized pre-formed FcpA/FcpB dimers predicted by computational methods. In order to construct an FcpA:FcpB complex, we first modeled L. biflexa FcpB using our L. interrogans experimental structure (PDB 6NQZ) as a template, using standard homology modeling procedures from Rosetta (Das and Baker, 2008). Through this procedure, more than 20,000 models of FcpB were generated, including loops β3β4 and β7β8, absent in the experimental structure because of poor electron density in those regions. The best 100 FcpB models, according to the Rosetta score, were used as input for an exhaustive full docking procedure with our FcpA X-ray model. The best 5000 docked configurations, according to the binding energy score, were rigidly fit into a segmented portion of the electron density map of the flagellar filament using SITUS colores and ranked according to the correlation coefficient. The highest-scoring solutions from this search clustered in two distinct molecular arrangements, one of which yielded a fair approximation of a single FcpA/FcpB heterodimer from our full FcpA/FcpB sheath model.
Curved FlaB core filament building and fitting
Request a detailed protocolTo build a complete atomic model of the curved filament, an 11-mer of modeled FlaB subunits (co-assembled according to (straight) helical symmetry) was first fit into the pseudo-helically symmetrized core volume. The 11-mer was then deformed to follow a curved path that matched the measured curvature of the tomographically averaged wild-type filament, using UCSF Chimera (Pettersen et al., 2004). This core model, corresponding to a single 52 Å repeat of the filament, was then merged with the best SITUS solutions for FcpA (6 protofilaments) and FcpB (4 protofilaments). The resulting core-sheath atomic model, representing a single 52 Å filament repeat, was then replicated 7 times along the same curved path to generate an assembly representing 8 52 Å repeats. To refine this filament model and eliminate minor clashes, all the protein monomers (protomers) placed in density were subjected to iterative cycles of rigid body, side-chain and backbone minimization using positional and conformational constraints, restrained within the cryoET volume map using a customized protocol in Rosetta (Das and Baker, 2008).
Curvature analysis
Request a detailed protocolCurvature was estimated for each reconstructed filament segment by comparing the 3D coordinates of neighboring filament subunits, and the resulting curvature values summed in histograms. The curvature of the purified flagella was determined in the manner of Crenshaw et al. (2000), utilizing the smoothed three-dimensional coordinates of the filaments. Consider consecutive three-dimensional points P1, P2, P3, …, PN, with each point spaced 52 Å apart (corresponding to the repeat spacing of the filament). Let C = Pi–P(i-10) and D = P(i+10)-Pi from i = 11 to i=(N-10); a spacing of 10 is used to minimize noise that might result from minor deviations between consecutive points. Then, curvature (κ) at a particular point can be determined by:
The final curvature value is then found by averaging the curvature at successive points, in the following manner:
The resultant values are multiplied by 1/(pixel size * binning), giving a curvature with units of 1/Å for points 11 to N-10 along the filament. These values were calculated for the 84 filaments of wild-type that were used in the emClarity analysis, the 14 filaments of the fcpB- sample that were used in the emClarity analysis, and 8 filaments of the fcpA- sample. Due to the small size of the dataset and sample heterogeneity, a high-resolution structure of the fcpA- sample could not be determined.
Estimating supercoil parameters from subtomograms
Request a detailed protocolFor each subtomogram average volume (wild-type or fcpB-), the volume was superposed on a 53 Å axially-shifted copy of itself using the UCSF Chimera command ‘Fit in Map’. Parameters for the superhelical rotation axis were then estimated using the Chimera command ‘measure rotation’. The superhelical diameter was estimated as twice the distance from the superhelical axis to the filament central channel, while the superhelical pitch was estimated as the reported volume shift along the superhelical axis multiplied by 360° divided by the reported volume rotation.
Data availability
Atomic models for three crystal forms of FcpA, and for FcpB, have been deposited in the Protein Data Bank under accession numbers 6NQW, 6NQX, 6NQY, and 6NQZ (respectively). Cryo-EM maps have been deposited in the Electron Microscopy Data Bank (EMDB) under accession number EMD-20504. A pseudo-atomic model consisting of main chain FlaB, FcpA and FcpB docked into the wild-type filament reconstruction has been deposited in the Protein Data Bank (PDB) under the accession number 6PWB.
-
RCSB Protein Data BankID 6NQW. Flagellar protein FcpA from Leptospira biflexa - hexagonal form.
-
RCSB Protein Data BankID 6NQX. Flagellar protein FcpA from Leptospira biflexa / primitive monoclinic form.
-
RCSB Protein Data BankID 6NQY. Flagellar protein FcpA from Leptospira biflexa / ab-centered monoclinic form.
-
RCSB Protein Data BankID 6NQZ. Flagellar protein FcpB from Leptospira interrogans.
-
RCSB Protein Data BankID 6PWB. Rigid body fitting of flagellin FlaB, and flagellar coiling proteins, FcpA and FcpB, into a 10 Angstrom structure of the asymmetric flagellar filament purified from Leptospira biflexa Patoc WT cells resolved via subtomogram averaging.
-
Electron Microscopy Data BankID 20504. 10 Angstrom structure of the asymmetric flagellar filament purified from Leptospira biflexa Patoc WT cells resolved via subtomogram averaging.
-
RCSB Protein Data BankID 5WJT. Cryo-EM structure of B. subtilis flagellar filaments N226Y.
-
RCSB Protein Data BankID 3A5X. L-type straight flagellar filament made of full-length flagellin.
References
-
Fast tomographic reconstruction on multicore computersBioinformatics 27:582–583.https://doi.org/10.1093/bioinformatics/btq692
-
A magic triangle for experimental phasing of macromoleculesActa Crystallographica Section D Biological Crystallography 64:1179–1182.https://doi.org/10.1107/S0907444908030266
-
BookLeptospiral MotilityIn: Stanier R. Y, Rogers H. J, Ward J. B, editors. Relations Between Structure and Function in the Prokaryotic Cell. in 28th Symposium of the Society for General Microbiology. Cambridge: Cambridge University Press. pp. 285–294.
-
Biochemical and cytological analysis of the complex periplasmic flagella from Spirochaeta aurantiaJournal of Bacteriology 170:4023–4032.https://doi.org/10.1128/JB.170.9.4023-4032.1988
-
Generation, representation and flow of phase information in structure determination: recent developments in and around SHARP 2.0Acta Crystallographica Section D Biological Crystallography 59:2023–2030.https://doi.org/10.1107/S0907444903017694
-
Axial filament involvement in the motility of Leptospira interrogansJournal of Bacteriology 137:1406–1412.https://doi.org/10.1128/JB.137.3.1406-1412.1979
-
Design requirements for the construction of bacterial flagellaJournal of Theoretical Biology 57:469–489.https://doi.org/10.1016/0022-5193(76)90016-3
-
A "mechanistic" explanation of the multiple helical forms adopted by bacterial flagellar filamentsJournal of Molecular Biology 425:914–928.https://doi.org/10.1016/j.jmb.2012.12.007
-
BookAlgorithm for three-dimensional reconstruction from the imperfect projection data provided by electron microscopyIn: Frank J, editors. Electron Tomography: Methods for Three-Dimensional Visualization of Structures in the Cell. New York: Springer. pp. 217–243.https://doi.org/10.1007/978-0-387-69008-7
-
Assembly and stability of the tobacco mosaic virus particleAdvances in Protein Chemistry 18:37–121.https://doi.org/10.1016/s0065-3233(08)60268-5
-
Morphology and dynamics of protruding spirochete periplasmic flagellaJournal of Bacteriology 174:832–840.https://doi.org/10.1128/JB.174.3.832-840.1992
-
The unique paradigm of spirochete motility and chemotaxisAnnual Review of Microbiology 66:349–370.https://doi.org/10.1146/annurev-micro-092611-150145
-
Cryo-electron tomography for structural characterization of macromolecular complexesCurrent Protocols in Protein Science 13:1713s65.https://doi.org/10.1002/0471140864.ps1713s65
-
Analysis of the three-dimensional trajectories of organisms: estimates of velocity, curvature and torsion from positional informationThe Journal of Experimental Biology 203:961–982.
-
Macromolecular modeling with rosettaAnnual Review of Biochemistry 77:363–382.https://doi.org/10.1146/annurev.biochem.77.062906.171838
-
The elastic basis for the shape of Borrelia burgdorferiBiophysical Journal 96:4409–4417.https://doi.org/10.1016/j.bpj.2009.02.066
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Bacterial nanomachines: the flagellum and type III injectisomeCold Spring Harbor Perspectives in Biology 2:a000299.https://doi.org/10.1101/cshperspect.a000299
-
An introduction to data reduction: space-group determination, scaling and intensity statisticsActa Crystallographica Section D Biological Crystallography 67:282–292.https://doi.org/10.1107/S090744491003982X
-
Polar, Peritrichous, and lateral flagella belong to three distinguishable flagellar familiesJournal of Molecular Biology 379:273–283.https://doi.org/10.1016/j.jmb.2008.04.012
-
Motility of the spirochete LeptospiraCell Motility and the Cytoskeleton 9:101–110.https://doi.org/10.1002/cm.970090202
-
The actin filament twist changes abruptly at boundaries between bare and cofilin-decorated segmentsJournal of Biological Chemistry 293:5377–5383.https://doi.org/10.1074/jbc.AC118.001843
-
XDSActa Crystallographica. Section D, Biological Crystallography 66:125–132.https://doi.org/10.1107/S0907444909047337
-
Helical transformations of Salmonella flagella in vitroJournal of Molecular Biology 106:167–186.https://doi.org/10.1016/0022-2836(76)90306-5
-
The shape and dynamics of the LeptospiraceaeBiophysical Journal 93:54–61.https://doi.org/10.1529/biophysj.106.103143
-
Structure of the native supercoiled flagellar hook as a universal jointNature Communications 10:5295.https://doi.org/10.1038/s41467-019-13252-9
-
BookThe design of self-assembling systems of equal unitsIn: Warren K. B, editors. Formation and Fate of Cell Organelles. New York: Academic Press. pp. 1–18.https://doi.org/10.1016/B978-0-12-395580-7.50006-X
-
Leptospira: the dawn of the molecular genetics era for an emerging zoonotic pathogenNature Reviews. Microbiology 7:736–747.https://doi.org/10.1038/nrmicro2208
-
BookCrystallography of biological macromoleculesIn: Arnold E, Himmel D. M, Rossmann M. G, editors. International Tables for Crystallography. Wiley Online Library. pp. 525–528.https://doi.org/10.1107/97809553602060000001
-
The spirochete FlaA periplasmic flagellar sheath protein impacts flagellar helicityJournal of Bacteriology 182:6698–6706.https://doi.org/10.1128/JB.182.23.6698-6706.2000
-
Spirochete periplasmic flagella and motilityJournal of Molecular Microbiology and Biotechnology 2:345–354.
-
MMDB and VAST+: tracking structural similarities between macromolecular complexesNucleic Acids Research 42:D297–D303.https://doi.org/10.1093/nar/gkt1208
-
Automated electron microscope tomography using robust prediction of specimen movementsJournal of Structural Biology 152:36–51.https://doi.org/10.1016/j.jsb.2005.07.007
-
Automated tilt series alignment and tomographic reconstruction in IMODJournal of Structural Biology 197:102–113.https://doi.org/10.1016/j.jsb.2016.07.011
-
LeptospirosisCurrent Opinion in Infectious Diseases 18:376–386.https://doi.org/10.1097/01.qco.0000178824.05715.2c
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
Purification, Ultrastructure, and composition of axial filaments from LeptospiraJournal of Bacteriology 98:264–280.https://doi.org/10.1128/JB.98.1.264-280.1969
-
UCSF chimera--a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Virulence of the zoonotic agent of leptospirosis: still terra incognita?Nature Reviews Microbiology 15:297–307.https://doi.org/10.1038/nrmicro.2017.5
-
Frealix: model-based refinement of helical filament structures from electron micrographsJournal of Structural Biology 186:234–244.https://doi.org/10.1016/j.jsb.2014.03.012
-
Crystallization of FcpA from Leptospira , a novel flagellar protein that is essential for pathogenesisActa Crystallographica Section F Structural Biology Communications 73:123–129.https://doi.org/10.1107/S2053230X17002096
-
RELION: implementation of a bayesian approach to cryo-EM structure determinationJournal of Structural Biology 180:519–530.https://doi.org/10.1016/j.jsb.2012.09.006
-
Substructure solution with SHELXDActa Crystallographica. Section D, Biological Crystallography 58:1772–1779.https://doi.org/10.1107/s0907444902011678
-
The PyMOL Molecular Graphics System, Version 2.1.0The PyMOL Molecular Graphics System, Version 2.1.0.
-
Crystal structure of VC1805, a conserved hypothetical protein from a Vibrio cholerae pathogenicity island, reveals homology to human p32Proteins: Structure, Function, and Bioinformatics 71:1563–1571.https://doi.org/10.1002/prot.21993
-
Torque transmission mechanism of the curved bacterial flagellar hook revealed by cryo-EMNature Structural & Molecular Biology 26:941–945.https://doi.org/10.1038/s41594-019-0301-3
-
Motility is crucial for the infectious life cycle of Borrelia burgdorferiInfection and Immunity 81:2012–2021.https://doi.org/10.1128/IAI.01228-12
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
The flagellar Cytoskeleton of the spirochetesJournal of Molecular Microbiology and Biotechnology 11:221–227.https://doi.org/10.1159/000094056
-
Flagellar motility of the pathogenic spirochetesSeminars in Cell & Developmental Biology 46:104–112.https://doi.org/10.1016/j.semcdb.2015.10.015
-
Conventions and workflows for using situsActa Crystallographica. Section D, Biological Crystallography 68:344–351.https://doi.org/10.1107/S0907444911049791
-
FcpB is a surface filament protein of the endoflagellum required for the motility of the spirochete LeptospiraFrontiers in Cellular and Infection Microbiology 8:130.https://doi.org/10.3389/fcimb.2018.00130
-
Structure and switching of bacterial flagellar filaments studied by X-ray fiber diffractionNature Structural Biology 5:125–132.https://doi.org/10.1038/nsb0298-125
-
Decreased potency of the Vibrio cholerae sheathed flagellum to trigger host innate immunityInfection and Immunity 76:1282–1288.https://doi.org/10.1128/IAI.00736-07
Article and author information
Author details
Funding
National Institutes of Health (R01 GM 110530)
- Kimberley H Gibson
- Megan R Brady
- Zhiguo Shang
- Charles Vaughn Sindelar
National Institutes of Health (U01 AI 088752)
- Elsio A Wunder
- Albert Ko
National Institutes of Health (R01 TW009504)
- Elsio A Wunder
- Albert Ko
National Institutes of Health (R01 AI052473)
- Elsio A Wunder
- Albert Ko
National Institutes of Health (R0 AI121207)
- Elsio A Wunder
- Albert Ko
Agencia Nacional de Investigación e Innovación (FCE_3_2016_1_126797)
- Felipe Trajtenberg
- Fabiana San Martin
- Ariel Mechaly
- Alejandro Buschiazzo
Agencia Nacional de Investigación e Innovación (ALI_1_2014_1_4982)
- Felipe Trajtenberg
- Fabiana San Martin
- Ariel Mechaly
- Alejandro Buschiazzo
Agencia Nacional de Investigación e Innovación (FCE_3_2016_1_126797)
- Mathieu Picardeau
Agence Nationale de la Recherche (ANR-18-CE15-0027-1)
- Mathieu Picardeau
Agence Nationale de la Recherche (ANR-08-MIE-018)
- Mathieu Picardeau
Pasteur International Joint Research Unit (Integrative Microbiology of Zoonotic Agents (IMiZA))
- Mathieu Picardeau
Institut Pasteur
- Mathieu Picardeau
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We gratefully acknowledge Garrett Debs for help with the RELION subtomogram averaging pre-processing step, including his customizations of published python scripts for Relion sub-tomogram averaging. We thank the IT Department from Institut Pasteur for protein modeling computations on the TARS cluster and the staff in the Yale High- Performance Computing facility for their maintenance of these facilities. We thank Prof. Felix Rey for providing insightful discussions about the manuscript.
Copyright
© 2020, Gibson et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,401
- views
-
- 283
- downloads
-
- 27
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Microbiology and Infectious Disease
Phage-derived peptidoglycan hydrolases (i.e. lysins) are considered promising alternatives to conventional antibiotics due to their direct peptidoglycan degradation activity and low risk of resistance development. The discovery of these enzymes is often hampered by the limited availability of phage genomes. Herein, we report a new strategy to mine active peptidoglycan hydrolases from bacterial proteomes by lysin-derived antimicrobial peptide-primed screening. As a proof-of-concept, five peptidoglycan hydrolases from the Acinetobacter baumannii proteome (PHAb7-PHAb11) were identified using PlyF307 lysin-derived peptide as a template. Among them, PHAb10 and PHAb11 showed potent bactericidal activity against multiple pathogens even after treatment at 100°C for 1 hr, while the other three were thermosensitive. We solved the crystal structures of PHAb8, PHAb10, and PHAb11 and unveiled that hyper-thermostable PHAb10 underwent a unique folding-refolding thermodynamic scheme mediated by a dimer-monomer transition, while thermosensitive PHAb8 formed a monomer. Two mouse models of bacterial infection further demonstrated the safety and efficacy of PHAb10. In conclusion, our antimicrobial peptide-primed strategy provides new clues for the discovery of promising antimicrobial drugs.
-
- Ecology
- Microbiology and Infectious Disease
Predicting how species diversity changes along environmental gradients is an enduring problem in ecology. In microbes, current theories tend to invoke energy availability and enzyme kinetics as the main drivers of temperature-richness relationships. Here, we derive a general empirically-grounded theory that can explain this phenomenon by linking microbial species richness in competitive communities to variation in the temperature-dependence of their interaction and growth rates. Specifically, the shape of the microbial community temperature-richness relationship depends on how rapidly the strength of effective competition between species pairs changes with temperature relative to the variance of their growth rates. Furthermore, it predicts that a thermal specialist-generalist tradeoff in growth rates alters coexistence by shifting this balance, causing richness to peak at relatively higher temperatures. Finally, we show that the observed patterns of variation in thermal performance curves of metabolic traits across extant bacterial taxa is indeed sufficient to generate the variety of community-level temperature-richness responses observed in the real world. Our results provide a new and general mechanism that can help explain temperature-diversity gradients in microbial communities, and provide a quantitative framework for interlinking variation in the thermal physiology of microbial species to their community-level diversity.