Structure of Mycobacterium tuberculosis cytochrome bcc in complex with Q203 and TB47, two anti-TB drug candidates
Abstract
Pathogenic mycobacteria pose a sustained threat to global human health. Recently, cytochrome bcc complexes have gained interest as targets for antibiotic drug development. However, there is currently no structural information for the cytochrome bcc complex from these pathogenic mycobacteria. Here, we report the structures of Mycobacterium tuberculosis cytochrome bcc alone (2.68 Å resolution) and in complex with clinical drug candidates Q203 (2.67 Å resolution) and TB47 (2.93 Å resolution) determined by single-particle cryo-electron microscopy. M. tuberculosis cytochrome bcc forms a dimeric assembly with endogenous menaquinone/menaquinol bound at the quinone/quinol-binding pockets. We observe Q203 and TB47 bound at the quinol-binding site and stabilized by hydrogen bonds with the side chains of QcrBThr313 and QcrBGlu314, residues that are conserved across pathogenic mycobacteria. These high-resolution images provide a basis for the design of new mycobacterial cytochrome bcc inhibitors that could be developed into broad-spectrum drugs to treat mycobacterial infections.
Editor's evaluation
The authors describe the cryo EM structures of cytochrome bcc from mycobacteria in complex with clinical drug candidates Q203 and TB47, which are being trialled for the treatment of tuberculosis. The cryo EM structures were obtained by expressing and purifying a hybrid super complex; the bcc part from M. tuberculosis and the cytochrome oxidase part from M. smegmatis. The cryo EM structures show how the Q203 and TB47 inhibitors bind and block the binding of quinone at the Qo site, thus inhibiting electron-transfer. The complex structures should facilitate the structure-based development of Q203 and TB47 drug compounds against tuberculosis.
https://doi.org/10.7554/eLife.69418.sa0Introduction
Mycobacteria, which belong to the phylum Actinobacteria, have coevolved with humans over thousands of years (Chisholm et al., 2016). Approximately 200 species of mycobacteria have been identified that have diverse lifestyles, morphologies, and metabolic pathways (Tortoli et al., 2017). Mycobacteria can broadly be grouped into two categories: tuberculosis-causing mycobacteria and non-tuberculous mycobacteria (NTM). Mycobacterium leprae is often represented in a distinct genetic clade owing to its genetic and phenotypic differences compared to other mycobacterial species (Cole et al., 2001). Mycobacteria can be further classified into fast-growing and slow-growing species or species complexes; these assignments are according to the physiological, phenotypic and phylogenetic characteristics (Rastogi et al., 2001). Although nearly all mycobacteria are saprophytes or non-pathogenic to humans, a few species cause diseases, resulting in pulmonary and extra-pulmonary infections that can affect nearly all organs. Infections, which are caused by strict or opportunistic pathogenic mycobacteria (Figure 1), pose a sustained threat to human health. Of these, tuberculosis (TB), caused by Mycobacterium tuberculosis (M. tuberculosis), is the most serious, leading to ~1.2 million fatalities per year (World Health Organization, 2019). Infections involving other pathogenic mycobacteria, for example, Mycobacterium abscessus and Mycobacterium avium complex, are on the rise with some outnumbering those caused by M. tuberculosis in countries including the United States (Donohue, 2018; Johansen et al., 2020). These infections are notoriously difficult to treat due to intrinsic or emerging resistance to many common antibiotics, thus exacerbating the challenge to find suitable drug targets.
Oxidative phosphorylation (OXPHOS) has gained interest as a target space for antibiotic drug development (Cook et al., 2017; Hards et al., 2020). In OXPHOS, the protein complexes of the electron transport chain (ETC) establish a proton motive force (PMF) across a biomembrane that drives the synthesis of adenosine triphosphate (ATP) by ATP synthase (Mitchell, 1961). Maintenance of PMF and ATP homeostasis is required for the survival of both replicative and non-replicative (often referred to as dormant) mycobacteria, and its dissipation leads to a rapid loss of cell viability and cell death (Koul et al., 2008; Rao et al., 2008). The reliance on the PMF and ATP homeostasis thus highlights the importance of the mycobacterial proton-pumping cytochrome bcc-aa3 supercomplex, which consists of a bcc menaquinol reductase (complex III, CIII) and an aa3 oxidase (complex IV, CIV) that are tightly associated (Gong et al., 2018; Kim et al., 2015; Megehee et al., 2006). Several studies support the attractiveness of cytochrome bcc-aa3 for mycobacterial drug development (de Jager et al., 2020; Liu et al., 2019; Lu et al., 2019; Pethe et al., 2013; Scherr et al., 2018). Given the strict sequence conservation of this complex (Figure 1), broad-spectrum activity of a therapeutic within the pathogenic mycobacteria is likely (Lee et al., 2020b). Interestingly, all the cytochrome bcc-aa3 inhibitors published to date appear to target the QcrB subunit (Figure 1) of the cytochrome bcc complex and are likely bound to the menaquinol-binding (Qp) site of the QcrB subunit (Lee et al., 2020b). The most advanced of these are Q203 and TB47, which have been shown to clear infections due to M. tuberculosis (de Jager et al., 2020; Lu et al., 2019; Pethe et al., 2013) and Mycobacterium ulcerans (Liu et al., 2019; Scherr et al., 2018). Q203 has recently completed phase II clinical trials for TB treatment (ClinicalTrials.gov number, NCT03563599; de Jager et al., 2020). TB47 has also been evaluated in a preclinical study (http://www.newtbdrugs.org/pipeline/clinical).
To progress an understanding of cytochrome bcc and its interaction with new drug leads, here we have determined the atomic resolution cryo-electron microscopy (cryo-EM) structures of M. tuberculosis cytochrome bcc alone (2.68 Å resolution) and in complex with Q203 (2.67 Å resolution) and TB47 (2.93 Å resolution). These high-resolution structures will greatly accelerate efforts towards structure-guided drug discovery.
Results and discussion
Structure of M. tuberculosis cytochrome bcc
Considering that the hybrid supercomplex consisting of M. tuberculosis CIII and Mycobacterium smegmatis CIV can be stabilized as a functional assembly (Kim et al., 2015), we first chose to express and purify this complex to a high level of homogeneity (Figure 2—figure supplement 1). We confirmed this complex to be active with a turnover number of 23.3 ± 2.4 e-s–1 (mean ± SD, n = 4), in agreement with the previous study that shows M. tuberculosis CIII can functionally complement native M. smegmatis CIII and maintain the growth of M. smegmatis (Kim et al., 2015). The structure of this hybrid supercomplex was determined by cryo-EM to an overall resolution of 2.68 Å, allowing that the components to be clearly assigned in the density map (Figure 2—figure supplements 2 and 3, Table 1). The hybrid supercomplex including the M. tuberculosis CIII dimer forms a pseudo twofold symmetrical compact rod, but with a slight curvature in the membrane plane, as previously observed in the M. smegmatis CIII2CIV2 supercomplex (Gong et al., 2018; Wiseman et al., 2018; Figure 2A and B). The hybrid supercomplex is structurally similar to the supercomplex isolated from M. smegmatis (Gong et al., 2018; Wiseman et al., 2018; Figure 2—figure supplement 4), except that subunits LpqE and Unk (probably MSMEG_0987) (Wiseman et al., 2018) were not observed here. The absence of these two subunits may be due to their depletion during purification of the supercomplex or their map density signal was averaged to background noise during structural determination. As expected, the topology of M. tuberculosis cytochrome bcc and M. smegmatis CIV in the hybrid supercomplex is also similar to that of the equivalent subunits in the M. smegmatis CIII2CIV2 supercomplex (Figure 2C and D, Figure 2—figure supplements 4 and 5). As such, there is no notable non-native contacts that resulted from the hybrid assembly in the hybrid supercomplex. This is attributable to the high structural similarity between the M. tuberculosis cytochrome bcc and M. smegmatis cytochrome bcc. In M. tuberculosis cytochrome bcc, QcrA has three transmembrane helices (TMHs) and has a ‘U’-shaped structure (Figures 2D and 3A). The N-terminal region with TMH1/2 and the TMH3 make up the two arms of the ‘U’ structure. These arms are linked by the region located near the cytoplasmic side. Attached to QcrATMH3 is the C-terminal domain, which faces the periplasmic side of the membrane and holds the [2Fe-2S] cluster in place. QcrB has eight TMHs (Figures 2D and 3B). Four of these are responsible for burying two functionally important heme b cofactors (high potential heme bH and low potential heme bL). Heme bL and heme bH are located between TMH I/II and TMH III/IV, respectively. Heme bL is near the periplasmic side and heme bH near the cytoplasmic side. QcrC is a transmembrane protein with a C-terminal TMH located between QcrBTMH5 and QcrBTMH7 (Figures 2D and 3C). The N-terminal periplasmic portion of QcrC can be divided into two heme-containing cytochrome c domains designated D1 and D2. The D1 domain protrudes out of the core of CIII, whereas the D2 domain interacts extensively with QcrA and QcrB. In the M. smegmatis CIII2CIV2 supercomplex, the cytochrome cc head domain of QcrC adopts an ‘open’ or a ‘closed’ conformation (Gong et al., 2018; Wiseman et al., 2018). However, in this structure it is only the closed conformation (Figure 2—figure supplement 4). Considering the high topology similarity of the subunits and the arrangement of prosthetic groups compared to previous M. smegmatis CIII2CIV2 supercomplex, the mechanism of action of the hybrid supercomplex including M. tuberculosis cytochrome bcc is expected to be the same as that for the M. smegmatis CIII2CIV2 supercomplex (Gong et al., 2018).
Quinone and quinone-binding pockets of M. tuberculosis cytochrome bcc
Quinone-binding sites are essential to the function of respiratory chain complexes and thus are good targets for drug discovery (Harikishore et al., 2021; Lee et al., 2020b). Structurally diverse quinones such as ubiquinone and menaquinone (MK) bind in the mitochondrial ETC (Zhang et al., 1998) and mycobacterial ETC (Gong et al., 2018; Wiseman et al., 2018), respectively, indicating that there are structural differences in the quinone-binding sites of the different species. Structural differences in the ubiquinone-binding site are also observed when the mitochondrial and Escherichia coli respiratory chain complex IIs are compared (Huang et al., 2021). Therefore, it is feasible that species-specific quinone-binding-site inhibitors could be designed. In the present structure, two canonical quinone-binding sites are identified (Figure 3D and E): (i) the quinol oxidation site (QP site) and (ii) the quinone reduction site (QN site). In addition, other quinone-binding sites with quinone bound are also identified as observed in our previous study (Gong et al., 2018; Figure 3—figure supplement 1). These quinones are remote from the canonical quinone-binding sites and into the membrane space, suggesting that they have a structural rather than functional role. The QP site responsible for menaquinol (MKH2) oxidation is near heme bL, whereas the QN site responsible for MK reduction is close to heme bH (Figure 3D and E, Figure 3—figure supplement 1). The QP site is at the center of an inverted triangle structure surrounded by helices (Figure 3D). One MK molecule was identified at this site with its naphthoquinone ring surrounded by hydrophobic residues, QcrBPhe158, QcrBTyr161, QcrBLeu180, QcrBIle183, QcrBMet310, and QcrBMet342. Its hydrophobic tail, which contains multiple isoprenoid groups, wraps around QcrBTMH6 down to its cytoplasmic end and in so doing interacts with QcrBMet187, QcrBAla339, QcrBLeu344, and QcrBVal347. The edge-to-edge distance from MK to heme bL is 15 Å, exceeding the 14 Å limit for efficient physiological electron transfer (Page et al., 1999). Furthermore, there are no observed hydrogen bonds to the carbonyl groups of MK, which are needed to help deprotonate MKH2. In addition, structural superposition with the inhibitor-bound bc1 complex also shows that the MK binds deeper into the QP pocket in the bc1 complex (Birth et al., 2014; Solmaz and Hunte, 2008) than in the currently described complex (Figure 3—figure supplement 2). Hence, the endogenous electron donor MKH2 should bind deeper inside the pocket and close to polar residues such as QcrBTyr161, QcrBThr313, and QcrBAsp314 (Figure 3D) to facilitate electron transfer. In addition, MK could also be modeled to fit deeper inside the pocket (Figure 3—figure supplement 3). Thus, we speculate that what is observed here is the oxidized product as it leaves the QP site. It is worth noting that all the reported inhibitors including Q203 (Pethe et al., 2013) and TB47 (Lu et al., 2019) are suggested to interact with this QP site. In addition, the QN site is formed largely by QcrBTMH1, QcrBTMH4, QcrBTMH5 and one loop region of QcrB (Figure 3E). The head group of MK is bound in this pocket interacting with QcrBPhe39, QcrBGlu49, QcrBLeu225, QcrBLeu232, QcrBTrp236, and QcrBPhe262, and its long hydrophobic tail extends along QcrBTMH1 towards the periplasmic side. MK/MKH2 are part of the Q-cycle hypothesis and essential for electron transfer in the cytochrome bcc complex (Gong et al., 2018). Given the crossspecies activity of this complex (Lee et al., 2020b) and high homology of the QcrB subunits across mycobacterial pathogens (Figure 1), these data open the way for the discovery of broad-spectrum mycobacterial agents based on rational, structure-based inhibitor design principles.
Q203 interactions in M. tuberculosis cytochrome bcc binding pocket
Q203 has recently been subjected to a phase II clinical study for M. tuberculosis treatment (de Jager et al., 2020). This compound has also been shown to be strongly bactericidal against M. ulcerans (Scherr et al., 2018). It is suggested to be an inhibitor that competes with endogenous substrate binding (QP site) of the cytochrome bcc complex (Pethe et al., 2013), but this hypothesis is yet to be verified by direct experimental evidence. To obtain atomic information on the mode of binding of Q203 to cytochrome bcc, we have determined the structure of the hybrid supercomplex as described above in the presence of Q203 by cryo-EM to an overall resolution of 2.67 Å (Figure 4—figure supplements 1 and 2, Table 2). We observe that close to the Qp-binding pocket within the membrane of each QcrB of cytochrome bcc there is density for Q203 (Figure 4A and B). All of the Q203 molecules fill each QcrB subunit binding deeply into the QP pocket and with identical binding modes. The key interactions that anchor Q203 are (i) a hydrogen bond between the hydroxyl oxygen of the side chain of QcrBThr313 and the amine in the carboxamide linker of Q203 (3.0 Å), (ii) a halogen bond between the chlorine atom of the heterocyclic group and an ordered water molecule that simultaneously forms a hydrogen bond with the hydroxyl oxygen of the side chain of QcrBTyr164 (Figure 4—figure supplement 3), (iii) a hydrogen bond between the side chain of QcrBGlu314 and the nitrogen atom in the imidazopyridine ring (3.0 Å), and (iv) a hydrogen bond between the side chain of QcrAHis375 and the nitrogen atom in the imidazopyridine ring (2.98 Å) (Figure 4C). In addition, the carbon atoms of Q203 make hydrophobic interactions with QcrBGly175, QcrBAla179, QcrBLeu180, QcrBThr184, QcrBSer304, QcrBPro306, QcrBMet310, QcrBAla317, and QcrBMet342. These extensive interactions are in agreement with the fact that the activity of the supercomplex is inhibited by Q203 in vitro according to the DMNQH2/oxygen oxidoreductase activity assay. After addition of Q203, the turnover number of the hybrid supercomplex reduces to 5.8 ± 2.4 e-s–1 from 23.3 ± 2.4 e-s–1 (mean ± SD, n = 4; Figure 4—figure supplement 4). In addition, functional studies have shown that substitution of QcrBThr313 to alanine confers Q203 resistance (Pethe et al., 2013). The binding of Q203 involves residues from both QcrA and QcrB. Due to the need to form stabilizing interactions between subunits, resistance may be more difficult to achieve here than if the binding site is within only one subunit. Consistently, the mapping of reported mutations in Q203-resistant M. tuberculosis reveals that they are positioned directly where Q203 binds in this structure (Lupien et al., 2020; Figure 4—figure supplement 5).
TB47-binding mode of M. tuberculosis cytochrome bcc
TB47, also currently being evaluated in preclinical studies, has been suggested to target the QcrB of cytochromes bcc from M. tuberculosis (Lu et al., 2019) and M. ulcerans (Liu et al., 2019). The 2.93 Å cryo-EM map shows density for TB47 and confirms that it binds in the same location as Q203 (Figure 5A and B, Figure 5—figure supplements 1 and 2, Table 3). Three hydrogen bond interactions are observed involving the side chains of QcrBThr313, QcrBGlu314, and QcrAHis375. Similar interactions are also observed when Q203 binds (Figure 5C). QcrBTyr161, QcrBLeu171, QcrBGly175, QcrBAla179, QcrBLeu180, QcrBThr184, QcrBMet187, QcrBLeu194, QcrBSer304, QcrBGly305, and QcrBMet342 also contribute to TB47 binding, largely through hydrophobic interactions (Figure 5C). Unlike when Q203 binds (Figure 4—figure supplement 3), there is no interaction between QcrBTyr164 and TB47 (Figure 5—figure supplement 3). However, after addition of TB47, the values of turnover number for the hybrid supercomplex reduced to 5.1 ± 2.9 e-s–1 (mean ± SD, n = 4; Figure 5—figure supplement 4), a similar value that was observed when Q203 binds. Thus, the absence of this does not greatly diminish inhibition. A mutation in TB47-resistant M. smegmatis (M. tuberculosis: H195Y) is close to the Qp-binding site (Lu et al., 2019; Figure 5—figure supplement 5). As a result of the change in shape, it would appear to contribute to steric interference with the binding of TB47, thus accounting for the observed resistance.
Specificity of Q203 and TB47 for mycobacterial cytochrome bcc complex
The basis for the high specificity of Q203 and TB47 toward the Qp site of mycobacterial cytochromes bcc becomes apparent in the structural comparison between the QcrB subunit of M. tuberculosis and counterparts from other species (Figure 6). The highly conserved residues that are involved in the binding of these two molecules in this region (Figure 6—figure supplement 1) suggest a consistent overall fold and binding site exists in mycobacteria. This is also in agreement with the fact that Q203 and TB47 show antimycobacterial activity across many species (de Jager et al., 2020; Liu et al., 2019; Lu et al., 2019; Pethe et al., 2013; Scherr et al., 2018). It is worth noting that the structures of cytochrome bcc from M. tuberculosis and M. smegmatis have high similarity (Figure 2—figure supplement 5), and no steric hindrance is observed between the Q203 and M. smegmatis cytochrome bcc (Figure 6). This observation indicates that Q203 should have a similar binding mechanism and a similar effect on the activity of cytochrome bcc from M. smegmatis and M. tuberculosis. This is in good agreement with previous antimycobacterial activity data and inhibition data for the bcc complexes from M. smegmatis and M. tuberculosis (Gong et al., 2018; Lu et al., 2018). In contrast, in other prokaryotic, eukaryotic and human Qp-binding pockets, for example, from Saccharomyces cerevisiae (Lange and Hunte, 2002), Rhodobacter sphaeroides (Esser et al., 2008), and human (Guo et al., 2017), many of the observed interactions would not be possible (Figure 6). This suggests that Q203 and TB47 should have low-binding affinity toward its counterpart QcrB in non-mycobacterial bacteria and in eukaryotes. Even if there is some flexibility in the Qp-binding pocket that enables some level of binding, key residues that enable the binding of Q203 and TB47 in the mycobacteria are not present in other bacteria and eukaryotes (Figure 6—figure supplement 1). In addition, the binding free energies of Q203 between wild-type and mutant hybrid supercomplex were calculated and used to further elucidate the importance of some key residues that could play a role in determining the inhibitor specificity. According to a previous functional study (Pethe et al., 2013), structural comparison (Figure 6), and sequence alignment (Figure 6—figure supplement 1), QcrBThr313 and QcrBGlu314 of hybrid supercomplex were selected and mutated to their counterparts in the human supercomplex, namely residues Ala and Tyr, respectively. It is believed that Ala would decrease the affinity of Q203 (Figure 4; Pethe et al., 2013) and the Tyr could lead to steric hindrance of Q203 binding (Figure 6). 75 ns simulations were used for the calculation of binding free energy (Figure 6—figure supplement 2). The different values for the root mean squared deviation (RMSD) suggest that three mutants might have the different effects on the structural conformation of QcrB, and the dual mutant (Thr313Ala and Glu314Tyr) induced significantly structural changes compared to the wild-type complex. The 75 ns simulation was further equally divided into three phases, and the average values were used as the final relative binding free energy (Table 4), which was calculated as the binding free energy of Q203 in the mutant minus that of the wild-type systems. It is suggested that the wild-type system has stronger affinity than the mutant systems (Thr313Ala or Glu314Tyr). The relative free energy of binding calculation for Q203 in the dual mutant (Thr313Ala and Glu314Tyr) indicates that its binding affinity is significantly weaker compared to the wild-type or single-mutant systems. Therefore, the suggestions drawn from the calculation of binding free energy of Q203 provide further evidence for specificity of Q203 inhibition. In summary, these observations correlate with the observed low general antibacterial activity and low cytotoxicity of Q203 and TB47 (Liu et al., 2019; Lu et al., 2019; Pethe et al., 2013; Scherr et al., 2018).
Implications of the Q203 and TB47 inhibitory mechanism
To gain further insights into the mechanism of action of Q203, we compared the structures of M. tuberculosis cytochrome bcc in the presence and absence of Q203 (Figure 7—figure supplement 1A). The structure of apo cytochrome bcc is almost identical with the Q203-bound structure (rmsd of 0.50 Å for all Ca atoms), which suggests that Q203 binding does not significantly affect the overall architecture of cytochrome bcc. A comparison of the Q203-bound and Q203-free cytochrome bcc structures shows residues involved in the binding pocket move outward, thus adapting to the shape of Q203 (Figure 7—figure supplement 1B). Specifically, the side chains of QcrBSer304, QcrBGlu313, QcrBGlu314, and QcrBMet342 undergo significant conformational changes to form hydrogen bonds with Q203. The binding of TB47 to M. tuberculosis cytochrome bcc also induces very similar conformational changes in the Qp-binding pocket to those seen for Q203 (rmsd of 0.454 Å for all Ca atoms) (Figure 7—figure supplement 1B). Differences in binding are due to the different ethyl group and methyl moieties in the head groups of Q203 and TB47. It is also important to note that one endogenous substrate molecule is also bound to the Qp site in the apo structure of cytochrome bcc, which potentially affects the evaluation of the conformational changes upon the binding of Q203 or TB47.
When analyzing the superimposed structures (Figure 7—figure supplement 1), it is apparent that Q203 and TB47 act competitively with the quinol binding as they almost completely occupy the Qp pocket. We therefore conclude that Q203 and TB47 are bona fide analogs of the substrate, and thus ultimately function by hindering the downstream synthesis of ATP (Figure 7). These two compounds are also highly bactericidal against M. ulcerans, almost certainly targeting the Qp-binding site (Liu et al., 2019; Scherr et al., 2018). In summary, the sequences of the QcrB subunits have high homology across pathogenic mycobacteria (Lee et al., 2020b) and the essential residues (QcrBGlu313 and QcrBGlu314) that are involved in hydrogen-bonding interactions with the inhibitors (Pethe et al., 2013; Scherr et al., 2018) are conserved across pathogenic mycobacteria (Figure 7—figure supplement 2).
Conclusions
We have determined the apo- and Q203 and Tb47-bound structures of a hybrid pathogenic M. tuberculosis/M. smegmatis cytochrome bcc complex. The study shows the structural features of M. tuberculosis cytochrome bcc and how it is specifically inhibited by Q203 and TB47. The extensive interactions between Q203 or TB47 and the Qp-binding pocket account for the highly specific binding of these two inhibitors to pathogenic M. tuberculosis cytochrome bcc compared to eukaryotic counterparts. Two conservative residues involved with the formation of hydrogen bonds are observed across the pathogenic mycobacteria. These structures provide a long-sought basis for rational, structure-based inhibitor design to accelerate the development of Q203 and TB47 analogs as drug leads for mycobacterial infections.
Materials and methods
Expression of hybrid supercomplex consisting of M. tuberculosis CIII and M. smegmatis CIV
Request a detailed protocolThe hybrid supercomplex was obtained according to a previous study (Kim et al., 2015) but with some modifications. The M. tuberculosis cytochrome bcc complex is encoded by three putative genes (Rv2194-2196). Genes were amplified from H37Rv genomic DNA by PCR using Phanta Max DNA polymerase (Vazyme), and two-step PCR was used to inset a 10× His tag at the C-terminus of the qcrB (Rv2196). Genes encoding the entire cytochrome bcc complex operon were then cloned into the pVV16 expression vector. The resultant plasmid was transformed into M. smegmatis mc2 51 (Li et al., 2014) cells whose qcrCAB operon encoding M. smegmatis cytochrome bcc had already been knocked out. The cells were cultivated in Luria-Bertani broth (LB) liquid media supplemented with 50 μg/mL hygromycin, 25 μg/mL streptomycin, and 0.1% Tween 80. Cell pellets were harvested by centrifugation when the cells were grown to an optical density (OD600 of 1.0–1.2) at 37°C (220 rpm). Harvested cells were frozen at –80°C until use.
Purification of the hybrid supercomplex
Request a detailed protocolCell pellets were thawed and resuspended in buffer A containing 20 mM 3-(N-morpholino)propanesulfonic acid (MOPS), pH 7.4, 100 mM NaCl, and then lysed by passing through a French Press at 1200 bar three times. Cell debris and non-lysed cells were removed by centrifugation at 14,000 rpm for 10 min at 4°C. The supernatant was collected and ultra-centrifuged at 36,900 rpm and 4°C for 2 hr. The membrane fraction was solubilized by addition of 1% (w/v) lauryl maltose neopentyl glycol (LMNG) in buffer A and incubated for 2 hr at 4°C with slow stirring. The suspension was ultra-centrifuged, and the supernatant was applied to Ni-NTA agarose beads (GE Healthcare) at 4°C. The beads were further washed in buffer A with 50 mM imidazole and 0.004% (w/v) LMNG. The buffer was exchanged to buffer B (20 mM MOPS, pH 7.4, 100 mM NaCl, and 0.1% [w/v] digitonin) and then washed in resin in batch mode. The protein was eluted from the beads with buffer B containing 500 mM imidazole. Protein was then concentrated and loaded onto a Superdex 6 increase (10/300GL, GE Healthcare) column equilibrated in buffer B. Peak fractions were pooled and concentrated to ~8 mg/mL for electron microscopy studies. The protein sample was analyzed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), and the bands were then identified through mass spectrometry.
Activity assays
Request a detailed protocol2,3-Dimethyl-1,4-naphthoquinone (DMNQ, CAS 2197-57-1) was synthesized by WuXi AppTec. In order to obtained reduced DMNQH2, 20 mM DMNQ was ultrasonically dissolved in 1 mL ethanol with 6 mM HCl and reduced with sodium borohydride (NaBH4) in the ice bath. An appropriate amount of 12 N HCl was added into quench unreacted NaBH4 under the protection of argon. Oxidase activity was determined by a method described previously (Bonner et al., 1986; Kusumoto et al., 2000; Safarian et al., 2019). The oxygen consumption was monitored using the Clark-type oxygen electrode (Oxytherm+, Hansatech) in 1 mL reaction buffer (20 mM MOPS, pH 7.4, 100 mM NaCl, 0.004% LMNG) at 25℃, containing 62.5 nM hybrid supercomplex. The reaction was started by addition of 100 μM DMNQH2. The oxygen consumption curve was plotted using GraphPad Prime 8.0 software, from which an estimate of the oxygen-reduction rate (v0) was obtained (corrected for autoxidation). In the inhibition assay, 500 nM Q203 or TB47, which is approximately the value of the median inhibitory concentration (IC50) according to our previous study (Gong et al., 2018), was chosen and incubated with 62.5 nM hybrid supercomplex for 20 min at 25℃. Inhibition curves and oxygen-reduction rates (vi) with the different inhibitors were recorded. The inhibition rate (1 v0/vi) of Q203 and TB47 at 500 nM is reported. This assay was conducted using four groups of parallel experiments.
Cryo sample preparation and data collection
Request a detailed protocol300-mesh Quantifoil R0.6/1.0 grids (Quantifoil, Micro Tools GmbH, Germany) were glow-discharged at H2/O2 atmosphere for 25 s. 3 μL aliquots of protein complex at a concentration of 10 mg/mL were applied to the grid and then blotted for 3 s with force 0 at 8°C and 100% humidity using a Vitrobot IV (Thermo). Images were collected using a Titan Krios 300 keV electron microscope (Thermo), equipped with K3 Summit direct electron detector camera (Gatan). Data were recorded at 29,000× magnification with a calibrated super-resolution pixel size 0.82 Å/pixel. The exposure time was set to 2.4 s with 40 subframes and a total dose of 60 electrons per Å2. All images were automatically recorded using SerialEM with a defocus range from 1.2 μm to 1.8 μm (Mastronarde, 2003). For the datasets of apo, Q203-bound and TB47-bound M. tuberculosis cytochrome bcc, a total of 4141, 3763, and 2968 images were collected, respectively.
Image processing
Request a detailed protocolAll dose-fractioned stacks were motion-corrected and dose-weighted using MotionCorr2 (Zheng et al., 2017) in RELION 3.03 (Zivanov et al., 2019). CTF estimation was conducted using cryoSPARC patch CTF estimation (Punjani et al., 2017). For the dataset of apo hybrid M. tuberculosis cytochrome bcc, 1,208,054 particles were picked automatically using EMD-9610 map as the template and extracted with a box size of rescaled 256 pixels (binned 2). 327,188 particles were selected after two rounds of 2D classification. 100,000 particles were used to perform ab initio reconstruction in four classes, and these four classes were used as 3D volume templates for heterogeneous refinement with all selected particles. 112,804 particles converged into one class with clear signals and then re-extracted with 512 pixels (binned 1). Next, this particle set was used to do homogeneous refinement and local refinement, yielding the final resolution 2.68 Å. For the dataset of Q203-bound and TB47-bound M. tuberculosis cytochrome bcc, the data processing was performed in a similar pipeline, resulting in the final reconstruction resolution at 2.67 Å and 2.93 Å, respectively (detailed parameters shown in supplementary figures).
Model building and validation
Request a detailed protocolThe M. smegmatis respiratory complex CIII2CIV2 (PDB code: 6ADQ) model (Gong et al., 2018) as rigid body was fitted into EM density maps using UCSF Chimera 1.12 (Pettersen et al., 2004). Next, the resultant atomic model was manually modified according to the subunit sequences of M. tuberculosis cytochrome bcc and refined in COOT 0.8.9.1 (Emsley et al., 2010), followed by real-space refinement in PHENIX (Adams et al., 2010). The smile strings of Q203 and TB47 were generated and copied from ChemDraw (Li et al., 2004) and defined in PHENIX elBOW. Q203 and TB47 were manually built into the corresponding EM densities. The local resolution map was calculated in cryoSPARC (Punjani et al., 2017). All reported resolutions were based on the gold-standard FSC 0.143 criteria (Rosenthal and Henderson, 2003).
MK9 docking study
Request a detailed protocolThe binding pose of MK9 with QcrB was generated by molecular docking using the Glide (Schrödinger, LLC). The spatially neighboring subunits QcrA/B were extracted from the hybrid supercomplex and used as the model. The protein structure was processed using the Protein Preparation Wizard module. In this process, hydrogens were added to heavy atoms and bond orders were assigned to each bond. The protonation state of each amino acid was predicted at pH 7.0 using the Epik algorithm (Shelley et al., 2007). All resolved waters were removed from the structure. At the same time, the structure of MK9 was prepared using the LigPrep module, and a low-energy conformation of MK9 was generated. The binding pocket of Q203 in QcrB was defined as the docking site using the Receptor Grid Generation module. The docking site was confined to an enclosed cube with side length of 20 Å, which was centered on the centroid of Q203. MK9 was then docked using the Ligand Docking module. Finally, 10 docking poses were generated for MK9 after post-docking minimization based on the force field of OPLS3 (Harder et al., 2016). These poses were ranked by the scoring function of GlideScore. The pose with the lowest GlideScore was chosen for binding mode analysis.
Binding free energies of Q203 between wild-type and mutant hybrid supercomplex
Request a detailed protocolBinding free energies of Q203 between wild-type and mutant hybrid supercomplex were computed using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method based on all-atom molecular dynamic simulations (Greene et al., 2019; Homeyer and Gohlke, 2012; Homeyer and Gohlke, 2015). This task can be divided into three phases: (i) simulation system construction, (ii) molecular dynamic simulation, and (iii) binding free energy computation.
Considering that QcrB is the core subunit interacting with Q203 in the hybrid supercomplex, the QcrB with Q203 bound was the only region extracted from the Q203-bound hybrid supercomplex for simulation system construction. Membrane Builder in CHARMM-GUI (Jo et al., 2007) was used to build the protein/membrane system solvated in water. PDB file containing the coordination information of Q203-bound QcrB was read while only protein residues and Q203 were retained. The CHARMM General Force Field (CGenFF) was applied to parameterize the Q203 (Vanommeslaeghe and MacKerell, 2012). DPPC lipid molecules were used to explicitly build the lipid bilayers and bulky water layers were added to the top and bottom sides of the membrane. In addition, 0.15 M KCl was added to neutralize the charge of system using the distance-based ion placing method. Next, the protein was inserted in the lipid bilayers and covered by the water layers despite exposure of some residues to the solvents. In order to generate topology and parameter files for subsequent AMBER dynamic simulations, the force fields of ff19SB (Tian et al., 2020) and Lipids17 (Lee et al., 2020a) were respectively applied to the parameterization of protein and lipid molecules. Four Q203-bound QcrB simulation systems were built, including wild-type QcrB, T313A mutant, E314Y mutant, and T313A/E314Y mutant.
Molecular dynamic simulation was carried out using AMBER 2020 (Case et al., 2005). Before dynamic simulation, 5000 steps of energy minimization were performed on waters to remove potential steric clashes between the solute and solvents. Harmonic restraints with a force constant of 10 kcal/mol Å2 and 2.5 kcal/mol Å2 were placed on protein and lipids in the minimization, respectively. Then, two stages of equilibrium simulations were carried out to relax the protein and lipid molecules by gradually decreasing the restraint force. The first stage was further divided into three short simulations each taking 375 ps with the time step set to 1 fs. The constant volume and temperature (NVT) ensemble was applied to this stage of simulation. The second one was a much longer simulation with the constant pressure and temperature (NPT) ensemble, which took 1.5 ns to completely relax the whole system with extremely low force. At this stage, the time step was increased to 2 fs while the SHAKE algorithm was used to fix the bonds with hydrogens (Bailey and Lowe, 2009). Subsequently, a 10 ns NPT simulation without any restraints was performed to further relax the system. Langevin dynamics was used to control the temperature at 303.15 K in NVT and NPT simulations (Grønbech-Jensen and Farago, 2014). In the NPT simulations, additional semi-isotropic pressure coupling and constant surface tension were applied to the simulation of lipid bilayers (Bennun et al., 2007). Once the equilibrium simulation was finished, a 100 ns NPT simulation free of restraints was carried out for each molecular system to produce the final trajectory used for energy calculation. Coordinates were printed every 100 ps so that in total 1000 frames were contained in the final trajectory. The CPPTRAJ module in AMBER was applied to analyze the trajectory (Roe and Cheatham, 2013).
MM/PBSA methods were employed to calculate the binding free energy of Q203 in wild-type and mutant QcrB (Homeyer and Gohlke, 2012). As shown in the equation (ΔG= EMM + Epol + Enp - TS), free energy can be calculated by combining the MM part (EMM) that represents the gas phase energy contribution as well as the solvation free energy components including polar (Epol) and non-polar (Enp) contributions. The gas phase free energy was calculated according to the force field. As for the solvation free energy, the polar part representing electrostatic contribution was calculated using the Poisson–Boltzmann (PB) equation based on implicit solvent model (Honig and Nicholls, 1995). In AMBER, the equation also provides an implicit membrane model to calculate the membrane-mediated electrostatic interactions. Here, geometric multigrid based on iterative solver was selected for PB equation calculation (Harris et al., 2013). On the other hand, the LCPO method was applied to calculate the non-polar part that represents the hydrophobic contribution (Weiser et al., 1999). The Python script MMPBSA.py (Miller et al., 2012) in AMBER was used for the computation based on previously produced simulation trajectories, from which one frame was extracted every 1 ns for the energy calculations. Ions, water, and lipid molecules were all striped from the molecular systems before calculation. For implicit membrane model, the membrane dielectric constant was set to 7.0, and membrane thickness was set to 42 Å. Entropy contributions (TS) could be neglected as the four molecular systems were assumed to have similar entropy changes due to the same protein-ligand structures except for the mutated residues.
Creation of figures
Request a detailed protocolAll the figures were created using UCSF Chimera (Pettersen et al., 2004) or PyMOL (Schrödinger LLC, 2017).
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided. The accession numbers for the 3D cryo-EM density map of apo, Q203-bound and TB47-bound hybrid supercomplex in present study are EMD-30943, EMD-30944 and EMD-30945, respectively. The accession numbers for the coordinates for the apo, Q203-bound and TB47-bound hybrid supercomplex in this study are PDB: 7E1V, PDB: 7E1W and PDB: 7E1X, respectively.
-
Electron Microscopy Data BankID EMD-30943. Cryo-EM structure of apo hybrid respiratory supercomplex consisting of Mycobacterium tuberculosis complexIII and Mycobacterium smegmatis complexIV.
-
Electron Microscopy Data BankID EMD-30944. Cryo-EM structure of hybrid respiratory supercomplex consisting of Mycobacterium tuberculosis complexIII and Mycobacterium smegmatis complexIV in the presence of Q203.
-
Electron Microscopy Data BankID EMD-30945. Cryo-EM structure of hybrid respiratory supercomplex consisting of Mycobacterium tuberculosis complexIII and Mycobacterium smegmatis complexIV in presence of TB47.
-
RCSB Protein Data BankID 7E1V. Cryo-EM structure of apo hybrid respiratory supercomplex consisting of Mycobacterium tuberculosis complexIII and Mycobacterium smegmatis complexIV.
-
RCSB Protein Data BankID 7E1W. Cryo-EM structure of hybrid respiratory supercomplex consisting of Mycobacterium tuberculosis complexIII and Mycobacterium smegmatis complexIV in the presence of Q203.
-
RCSB Protein Data BankID 7E1X. Cryo-EM structure of hybrid respiratory supercomplex consisting of Mycobacterium tuberculosis complexIII and Mycobacterium smegmatis complexIV in presence of TB47.
References
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica. Section D, Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
MILCH SHAKE: An efficient method for constraint dynamics applied to alkanesJournal of Computational Chemistry 30:2485–2493.https://doi.org/10.1002/jcc.21237
-
Simulations of biomembranes and water: Important technical aspectsFluid Phase Equilibria 261:18–25.https://doi.org/10.1016/j.fluid.2007.07.056
-
The Amber biomolecular simulation programsJournal of Computational Chemistry 26:1668–1688.https://doi.org/10.1002/jcc.20290
-
Telacebec (Q203), a new antituberculosis agentThe New England Journal of Medicine 382:1280–1281.https://doi.org/10.1056/NEJMc1913327
-
Features and development of CootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/s0907444910007493
-
Inhibitor-complexed structures of the cytochrome bc1 from the photosynthetic bacterium Rhodobacter sphaeroidesThe Journal of Biological Chemistry 283:2846–2857.https://doi.org/10.1074/jbc.M708608200
-
Heterogeneous dielectric implicit membrane model for the calculation of MMPBSA binding free energiesJournal of Chemical Information and Modeling 59:3041–3056.https://doi.org/10.1021/acs.jcim.9b00363
-
Constant pressure and temperature discrete-time Langevin molecular dynamicsThe Journal of Chemical Physics 141:194108.https://doi.org/10.1063/1.4901303
-
OPLS3: a force field providing broad coverage of drug-like small molecules and proteinsJournal of Chemical Theory and Computation 12:281–296.https://doi.org/10.1021/acs.jctc.5b00864
-
Two for the price of one: Attacking the energetic-metabolic hub of mycobacteria to produce new chemotherapeutic agentsProgress in Biophysics and Molecular Biology 152:35–44.https://doi.org/10.1016/j.pbiomolbio.2019.11.003
-
Targeting the menaquinol binding loop of mycobacterial cytochrome bd oxidaseMolecular Diversity 25:517–524.https://doi.org/10.1007/s11030-020-10034-0
-
Influence of grid spacing in Poisson-Boltzmann equation binding energy estimationJournal of Chemical Theory and Computation 9:3677–3685.https://doi.org/10.1021/ct300765w
-
Free energy calculations by the molecular mechanics Poisson-Boltzmann surface area methodMolecular Informatics 31:114–122.https://doi.org/10.1002/minf.201100135
-
Extension of the free energy workflow FEW towards implicit solvent/implicit membrane MM–PBSA calculationsBiochimica et Biophysica Acta (BBA) - General Subjects 1850:972–982.https://doi.org/10.1016/j.bbagen.2014.10.013
-
Crystallographic investigation of the ubiquinone binding site of respiratory Complex II and its inhibitorsBiochimica et Biophysica Acta. Proteins and Proteomics 1869:140679.https://doi.org/10.1016/j.bbapap.2021.140679
-
Non-tuberculous mycobacteria and the rise of Mycobacterium abscessusNature Reviews. Microbiology 18:392–407.https://doi.org/10.1038/s41579-020-0331-1
-
The Phyre2 web portal for protein modeling, prediction and analysisNature Protocols 10:845–858.https://doi.org/10.1038/nprot.2015.053
-
Diarylquinolines are bactericidal for dormant mycobacteria as a result of disturbed ATP homeostasisThe Journal of Biological Chemistry 283:25273–25280.https://doi.org/10.1074/jbc.M803899200
-
CHARMM-GUI supports the Amber force fieldsThe Journal of Chemical Physics 153:035103.https://doi.org/10.1063/5.0012280
-
Targeting the cytochrome oxidases for drug development in mycobacteriaProgress in Biophysics and Molecular Biology 152:45–54.https://doi.org/10.1016/j.pbiomolbio.2020.02.001
-
Personal experience with four kinds of chemical structure drawing software: review on ChemDraw, ChemWindow, ISIS/Draw, and ChemSketchJournal of Chemical Information and Computer Sciences 44:1886–1890.https://doi.org/10.1021/ci049794h
-
SerialEM: a program for automated tilt series acquisition on Tecnai microscopes using prediction of specimen positionMicroscopy and Microanalysis 9:1182–1183.https://doi.org/10.1017/S1431927603445911
-
MMPBSA.py: An efficient program for end-state free energy calculationsJournal of Chemical Theory and Computation 8:3314–3321.https://doi.org/10.1021/ct300418h
-
Discovery of Q203, a potent clinical candidate for the treatment of tuberculosisNature Medicine 19:1157–1160.https://doi.org/10.1038/nm.3262
-
UCSF Chimera-A visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
The mycobacteria: an introduction to nomenclature and pathogenesisRevue Scientifique et Technique 20:21–54.https://doi.org/10.20506/rst.20.1.1265
-
PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory dataJournal of Chemical Theory and Computation 9:3084–3095.https://doi.org/10.1021/ct400341p
-
SoftwareThe PyMOL molecular graphics system, version 2.0Schrödinger LLC.
-
Epik: a software program for pK(a) prediction and protonation state generation for drug-like moleculesJournal of Computer-Aided Molecular Design 21:681–691.https://doi.org/10.1007/s10822-007-9133-z
-
Structure of complex III with bound cytochrome c in reduced state and definition of a minimal core interface for electron transferThe Journal of Biological Chemistry 283:17542–17549.https://doi.org/10.1074/jbc.M710126200
-
ff19SB: Amino-acid-specific protein backbone parameters trained against quantum mechanics energy surfaces in solutionJournal of Chemical Theory and Computation 16:528–552.https://doi.org/10.1021/acs.jctc.9b00591
-
The new phylogeny of the genus Mycobacterium: The old and the newsFection, Genetics and Evolution 56:19–25.https://doi.org/10.1016/j.meegid.2017.10.013
-
Automation of the CHARMM General Force Field (CGenFF) I: bond perception and atom typingJournal of Chemical Information and Modeling 52:3144–3154.https://doi.org/10.1021/ci300363c
-
Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO)Journal of Computational Chemistry 20:217–230.https://doi.org/10.1002/(SICI)1096-987X(19990130)20:2<217::AID-JCC4>3.0.CO;2-A
-
Structure of a functional obligate complex III2IV2 respiratory supercomplex from Mycobacterium smegmatisNature Structural & Molecular Biology 25:1128–1136.https://doi.org/10.1038/s41594-018-0160-3
-
SoftwareGlobal Tuberculosis ReportWHO.
Article and author information
Author details
Funding
National Key Research and Development Program of China (2017YFC0840300)
- Zihe Rao
National Key Research and Development Program (2020YFA0707500)
- Zihe Rao
National Natural Science Foundation of China (81520108019)
- Zihe Rao
National Natural Science Foundation of China (813300237)
- Zihe Rao
National Natural Science Foundation of China (32100976)
- Hongri Gong
Natural Science Foundation of Tianjin City (20JCQNJC01430)
- Hongri Gong
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Chao Peng of the Mass Spectrometry System at the National Facility for Protein Science in Shanghai (NFPS), Zhangjiang Lab, SARI, China, for data collection and analysis and Prof. Lei Li (College of Life Sciences, Nankai University) for their technical support on Clark-type oxygen electrode and oxygen consumption assay. We thank Prof. Kaixia Mi (CAS Key Laboratory of Pathogenic Microbiology and Immunology, Institute of Microbiology, CAS) for sharing the strain M. smegmatis mc2 51. We would like to thank Prof. Gregory M Cook (School of Biomedical Sciences, University of Otago, New Zealand) and Prof. Xiaoyun Lu (School of Pharmacy, Ji′nan University, China) for kindly providing TB47 for this study. We would also like to thank the Bio-Electron Microscopy Facility of ShanghaiTech University.
Copyright
© 2021, Zhou 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,272
- views
-
- 311
- downloads
-
- 31
- 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
-
- Structural Biology and Molecular Biophysics
Cyclic nucleotide-binding domain (CNBD) ion channels play crucial roles in cellular-signaling and excitability and are regulated by the direct binding of cyclic adenosine- or guanosine-monophosphate (cAMP, cGMP). However, the precise allosteric mechanism governing channel activation upon ligand binding, particularly the energetic changes within domains, remains poorly understood. The prokaryotic CNBD channel SthK offers a valuable model for investigating this allosteric mechanism. In this study, we investigated the conformational dynamics and energetics of the SthK C-terminal region using a combination of steady-state and time-resolved transition metal ion Förster resonance energy transfer (tmFRET) experiments. We engineered donor-acceptor pairs at specific sites within a SthK C-terminal fragment by incorporating a fluorescent noncanonical amino acid donor and metal ion acceptors. Measuring tmFRET with fluorescence lifetimes, we determined intramolecular distance distributions in the absence and presence of cAMP or cGMP. The probability distributions between conformational states without and with ligand were used to calculate the changes in free energy (ΔG) and differences in free energy change (ΔΔG) in the context of a simple four-state model. Our findings reveal that cAMP binding produces large structural changes, with a very favorable ΔΔG. In contrast to cAMP, cGMP behaved as a partial agonist and only weakly promoted the active state. Furthermore, we assessed the impact of protein oligomerization and ionic strength on the structure and energetics of the conformational states. This study demonstrates the effectiveness of time-resolved tmFRET in determining the conformational states and the ligand-dependent energetics of the SthK C-terminal region.
-
- Structural Biology and Molecular Biophysics
Ciliary rootlets are striated bundles of filaments that connect the base of cilia to internal cellular structures. Rootlets are critical for the sensory and motile functions of cilia. However, the mechanisms underlying these functions remain unknown, in part due to a lack of structural information of rootlet organization. In this study, we obtain 3D reconstructions of membrane-associated and purified rootlets from mouse retina using cryo-electron tomography. We show that flexible protrusions on the rootlet surface, which emanate from the cross-striations, connect to intracellular membranes. In purified rootlets, the striations were classified into amorphous (A)-bands, associated with accumulations on the rootlet surface, and discrete (D)-bands corresponding to punctate lines of density that run through the rootlet. These striations connect a flexible network of longitudinal filaments. Subtomogram averaging suggests the filaments consist of two intertwined coiled coils. The rootlet’s filamentous architecture, with frequent membrane-connecting cross-striations, lends itself well for anchoring large membranes in the cell.