A structural mechanism for bacterial autotransporter glycosylation by a dodecameric heptosyltransferase family
Abstract
A large group of bacterial virulence autotransporters including AIDA-I from diffusely adhering E. coli (DAEC) and TibA from enterotoxigenic E. coli (ETEC) require hyperglycosylation for functioning. Here we demonstrate that TibC from ETEC harbors a heptosyltransferase activity on TibA and AIDA-I, defining a large family of bacterial autotransporter heptosyltransferases (BAHTs). The crystal structure of TibC reveals a characteristic ring-shape dodecamer. The protomer features an N-terminal β-barrel, a catalytic domain, a β-hairpin thumb, and a unique iron-finger motif. The iron-finger motif contributes to back-to-back dimerization; six dimers form the ring through β-hairpin thumb-mediated hand-in-hand contact. The structure of ADP-D-glycero-β-D-manno-heptose (ADP-D,D-heptose)-bound TibC reveals a sugar transfer mechanism and also the ligand stereoselectivity determinant. Electron-cryomicroscopy analyses uncover a TibC–TibA dodecamer/hexamer assembly with two enzyme molecules binding to one TibA substrate. The complex structure also highlights a high efficient hyperglycosylation of six autotransporter substrates simultaneously by the dodecamer enzyme complex.
https://doi.org/10.7554/eLife.03714.001eLife digest
Bacteria release proteins known as virulence factors to help them infect host cells. Many bacteria are surrounded by two membranes, so virulence factors must be able to pass through both of these membranes.
Autotransporters are a group of virulence factors that pass through the inner membrane and anchor themselves in the outer membrane; this allows part of the autotransporter to project from the surface of the bacterial cell and stick to the surface of the cell that the bacterium is about to infect. Many autotransporters are coated with sugar molecules and this increases their ability to adhere to cells. Enzymes called glycosyltransferases ensure that this sugar coating process takes place.
Autotransporters contain two sections: the passenger domain and the beta domain. The passenger domain is important for virulence, while the beta domain forms a pore in the outer membrane that the passenger domain passes through to reach the outer surface of the bacterium. TibA and AIDA-1 are autotransporters associated with two types of bacteria that infect the intestines and cause diarrhea. Similar autotransporters are found in a wide range of bacteria, but the precise details of how these autotransporters are coated with sugar molecules are not fully understood.
Yao et al. now show that a glycosyltransferase called TibC, which is found in many bacteria, adds large numbers of sugar molecules to the passenger domains of both TibA and AIDA-1. To learn more about this process Yao et al. used X-ray diffraction to work out the structure of TibC. Strikingly, this revealed that TibC proteins come together to form a large circular structure that contains two rings, each made of six TibC proteins. The integrity of this structure is maintained by the presence of iron atoms, which also gives TibC a characteristic brown colour.
Yao et al. also studied what happens when TibC binds to TibA; a technique called electron cryo-microscopy revealed that six TibA molecules are distributed along the inner surface of the circular TibC structure, with each TibA protein binding to two TibC proteins. This arrangement allows for the efficient transfer of sugar molecules from the glycosyltransferase to the autotransporter.
https://doi.org/10.7554/eLife.03714.002Introduction
Protein glycosylation is one of the most abundant post-translational modifications in all domains of life (Spiro, 2002). Recent studies have appreciated protein glycosylation in bacteria, which is often associated with pathogen virulence and immune modulation (Szymanski and Wren, 2005; Abu-Qarn et al., 2008; Nothaft and Szymanski, 2010). For instance, flagella glycosylation is found in many bacteria species including Campylobacter jejuni, Helicobacter pylori, Clostridium spp. and Pseudomonas aeruginosa, contributing to bacterial locomotion or virulence (Schirm et al., 2003; Nothaft and Szymanski, 2010). Recently, we and others have shown that secreted effectors from enteropathogenic Escherichia coli (EPEC) and related enteric pathogens harbor an arginine N-acetylglucosamine transferase activity that modifies host death-domain proteins and is essential for bacterial colonization in infected mice (Li et al., 2013; Pearson et al., 2013).
The autotransporter secretion pathway (also known as the type V secretion system) in bacteria delivers autotransporters onto the bacterial surface. Autotransporters, representing the largest family of bacterial virulence factors, share a similar structural organization containing a signal peptide followed by a functional passenger domain and a C-terminal β-barrel translocation domain. Autotransporters play a critical role in diverse aspects of bacterial physiology including proteolytic digestion of host proteins, biofilm formation, adhesion and invasion of host cells, and intracellular motility (Henderson et al., 2004; Lazar Adler et al., 2011; Wells et al., 2007). A subfamily of autotransporters, including AIDA-I (adhesin involved in diffuse adherence) from diffusely adhering E. coli (DAEC) 2787 (Benz and Schmidt, 1989) and TibA from enterotoxigenic E. coli (ETEC) H10407 (Elsinghorst and Weitz, 1994), are glycosylated in their passenger domains (Lindenthal and Elsinghorst, 1999; Benz and Schmidt, 2001; Sherlock et al., 2006), which functions in bacterial autoaggregation and adhesion to host cells (Benz and Schmidt, 1989; Sherlock et al., 2004; Charbonneau and Mourez, 2007). AIDA-I and TibA-like autotransporters are present in diverse bacterial species. Previous studies carried out in the ectopic system (Benz and Schmidt, 2001; Moormann et al., 2002) indicate a role for AAH (autotransporter adhesin heptosyltransferase) in DAEC and its ETEC homologue TibC in AIDA-I and TibA glycosylation. However, AAH and TibC harbor no sequence homology to known glycosyltransferases and there has been no biochemical evidence demonstrating their glycosyltransferase activity.
In a separate parallel study (Lu et al., 2014) we showed that AAH is a bona fide heptosyltransferase belonging to a large bacterial autotransporter heptosyltransferase (BAHT) family. Here we determine the crystal structures of TibC heptosyltransferase, both alone and in complex with ADP-D-glycero-β-D-manno-heptose (ADP-D,D-heptose). The structure shows a symmetric ring-shape dodecamer. The protomer features a β-hairpin thumb and an iron-finger motif, both required for the dodecamer assembly. The ligand-bound structure reveals a sugar transfer mechanism and determinants for the ligand stereoselectivity for TibC and AAH. Furthermore, electron cryomicroscopy (cryo-EM) analyses of a TibC–TibA dodecamer/hexamer enzyme–substrate complex reveal the structural basis for high efficient autotransporter hyperheptosylation by the TibC dodecamer.
Results
TibC catalyzes AIDA-I/TibA heptosylation and confers AIDA-I-mediated bacterial adhesion to host cells
We first confirmed the autotransporter heptosyltransferase activity of TibC and its function in mediating bacterial adhesion to host cells (Lu et al., 2014). Co-expression of TibC or AAH together with AIDA-I in E. coli BL21 cells resulted in a tight adhesion of the bacteria to HeLa cells (Figure 1A) and AIDA-I glycosylation (Figure 1B). Co-expression of TibC also induced TibA glycosylation in the heterologous system (Figure 1C). We further identified a minimal passenger domain fragment of TibA (TibA305-350) that could be efficiently glycosylated by TibC (Figure 1D,E). TibA305-350 tagged with Flag epitopes at both termini was then fused C-terminal to the Small Ubiquitin-like Modifier (SUMO). The SUMO-Flag-TibA305-350-Flag fusion protein showed a slower migration on the SDS-PAGE gel when co-expressed with TibC (Figure 1F). Electrospray ionization (ESI) mass spectrometry analysis gave a mass of 20,758 Da that matched the theoretical mass (20,757 Da) of SUMO-Flag-TibA305-350-Flag (Figure 1G). Co-expression of TibC resulted in mass increase of 576, 768, or 960 Da, corresponding to the addition of 3, 4, and 5 heptoses, respectively (Figure 1G). The modified fusion protein was further digested with the endoproteinase Asp-N for liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. Asp-N digestion is expected to yield two major peptide fragments, DK-305SASKVIQNSGGAVITNTSAAVSGTN329 and 330DNGSFSIAGGSAVNMLLENGG350. We were able to detect abundant signals for both unmodified and modified DK-305SASKVIQNSGGAVITNTSAAVSGTN329 in full scan MS analyses and only weak signals of the modified 330DNGSFSIAGGSAVNMLLENGG350 peptide, probably due to the hydrophobic nature. Mass measurements indicated that both peptides had multiple heptose modifications; a mixed population of modified peptides, containing different numbers of heptose conjugations, was observed. Similarly to that performed with AAH modification of AIDA-I (Lu et al., 2014), serine residues (Ser-313 and Ser-322 within DK-305SASKVIQNSGGAVITNTSAAVSGTN329), but not threonine residues, were identified to be the modification sites by electron transfer dissociation (ETD) mass spectrometry (Figure 1H). Supporting the heptosyltransferase activity of TibC, recombinant TibC was found to directly modify synthetic peptides derived from the passenger domains of either AIDA-I or TibA in which ADP-D,D-heptose, but not the anomer ADP-L,D-heptose, could serve as the sugar ligand (Figure 1I).
AAH and TibC homologues, defined as the BAHT family (Lu et al., 2014), are widely present in pathogenic bacteria species including Citrobacter rodentium, Salmonella enterica serovar Urbana, Shigella sp. D9, Laribacter hongkongensis, Cronobacter sakazakii, and several Burkholderia species (Figure 1—figure supplement 1). These proteins generally bear more than 50% sequence homology to TibC/AAH and contain residues critical for the heptosyltransferase activity (see below). In the genomic locus, genes encoding BAHT are often followed by an autotransporter or a putative autotransporter, consistent with the notion that BAHT modifies its autotransporter partner for functioning (Lu et al., 2014). Consistent with the modification of AIDA-I by TibC, the passenger domain of AIDA-I shows sequence similarities to that of TibA and contains a number of heptosylation motifs identified in TibA previously (Lu et al., 2014) (Figure 1—figure supplement 2).
Crystal and cryo-EM structures of TibC dodecamer
To reveal the mechanism of the BAHT family, we attempted to solve the crystal structure. Among several BAHTs analyzed, TibC behaved the best in recombinant expression and homogeneity. Purified His6-TibC was crystallized, but the good-looking crystals did not diffract due to internal disorders. TibC obtained from glutathione-S-transferase (GST)-fusion expression was also crystallized, which only diffracted to 8 Å despite extensive optimizations. A number of point mutations designed to reduce the surface entropy were screened and a TibC variant with three mutated epitopes (E83A/E84A, K400A/K401A, and Q215A/E216A) produced the highest quality crystal. This mutant exhibited identical biochemical and biological functions to wild-type TibC (Figure 1A–C). A 2.9 Å structure was solved by single wavelength dispersion phasing using anomalous signals from selenomethionine (SeMet) and the naturally bound ferric ion (Table 1).
Strikingly, one asymmetric unit in the solved structure contains 12 TibC molecules, which are assembled into a large ring structure (Figure 2A). The overall shape of the dodecamer ring resembles a circular garland with an external diameter of 145 Å and a height of 72 Å. Parts of the 12 protomers are lined side by side at the central plate to form the stem of the garland with an inner diameter of ∼110 Å (Figure 2A). The inner surface of the stem is decorated with uniformly distributed 12 ferric ions, forming a characteristic iron belt along the ring (Figure 2A). The presence of ferric ions is expected from the brownish color of recombinant TibC protein in solution. Six symmetric TibC protomers project upwards from the central plate while the other six project downwards, generating a gear-like shape at both sides. Two adjacent TibC protomers projecting oppositely contact each other in a back-to-back manner around a dyad axis; six back-to-back dimers are further linked together in a hand-in-hand fashion around a sixfold axis to form the garland-like ring (Figure 2B). Thus, the dodecamer assembly can be viewed as a hexameric oligomerization of dimers. The back-to-back interface is extensive and tilts across the central plate by 45° (Figure 2A,B). The hand-in-hand connection creates 12 tilted grooves extending to the interior of the garland from both up and down sides.
To investigate whether the dodecamer observed in the crystal structure results from crystal packing, single particle cryo-EM analyses of highly purified TibC protein were performed (Figure 2—figure supplement 1). A 3D reconstruction of 11.5 Å was obtained from the cryo-EM images, which showed a remarkably similar dodecamer architecture to that determined by X-ray crystallography (Figure 2C). This suggests that the dodecamer ring architecture is the predominant form of TibC protein in solution and likely represents a functional physiological state. Consistently, purified TibC eluted as an oligomer of ∼500 kDa on the gel filtration column (see below) and the molecular weight was determined to be 578 kDa by analytic ultracentrifugation (Figure 2—figure supplement 2).
Structure of TibC protomer and the glycosyltransferase domain
All the 12 TibC protomers adopt an identical structure containing an N-terminal β-barrel domain, a core catalytic domain, a β-hairpin thumb, and an iron-finger motif which are arranged into a palm-like shape (Figure 3A). The β-hairpin thumb and the iron-finger motif are insertions into the catalytic domain. The flattened β-barrel (residues 1–97) is structured from two four-stranded sheets that pack against each other. The barrel is stacked against the catalytic domain from the side and makes extensive contacts with β9, β12, and the loop linking α3 and β12; several hydrophobic residues from both domains are in close proximity at the interface (Figure 3—figure supplement 1A). The β-barrel is well separated from the catalytic domain and is unlikely to play a direct role in sugar transfer. Notably, deletion of the N-terminal 80 or 95 residues abolished TibC heptosylation of TibA (Figure 3—figure supplement 1B), indicating a possible function of the β-barrel domain in either substrate recognition or maintaining structural integrity of the enzyme.
The catalytic domain contains two lobes (residues 98–195 and 218–406, respectively) linked by a long loop (residues 196–217). The N-lobe, containing β9–12, α1–4 and surrounding loops, adopts a compact Rossmann-like globular fold. Different from the classical Rossmann fold that has α-helices flanking the central β-sheet, the β-sheet in the N-lobe is flanked by α-helices (α1–4) at one side but the β-barrel at the other side (Figure 3A). The C-lobe (α5–11 and β13–18) features a central β-sheet (β13 and β16–18) sandwiched by α7/α8 at one side and α6/α11 at the other side, adopting a typical Rossmann fold.
The catalytic cores of nearly all known glycosyltransferases, despite the extraordinary sequence divergence, converge onto two general folds—namely, GT-A and GT-B (Unligil and Rini, 2000; Breton et al., 2006; Lairson et al., 2008; Breton et al., 2012). The typical GT-A fold consists of two α/β/α domains with a continuous central β-sheet. The GT-A glycosyltransferase features a D×D motif that binds to the nucleotide phosphate in the sugar donor via a metal ion (Bourne and Henrissat, 2001). The GT-B fold has two Rossmann-like domains separated by a cleft, in which a sugar donor is positioned for nucleophilic attack by the substrate. The GT-B glycosyltransferase does not require metal ion for catalyzing sugar transfer (Wrabl and Grishin, 2001; Hu and Walker, 2002). The catalytic domain of TibC structurally more resembles the GT-B fold such as WaaC and MurG, two bacterial glycosyltransferases involved in cell wall synthesis (Hu et al., 2003; Grizot et al., 2006), with its own distinctions (Figure 3B). One major difference is that the N-lobe of the TibC catalytic domain is not a typical Rossmann fold; its central β-sheet is four-stranded rather than six-stranded for a Rossmann fold and also misses the flanking α-helices at one side (Figure 3B).
The β-hairpin thumb and hand-in-hand contact
The TibC catalytic core is accessorized by two unique structural modules, the β-hairpin thumb and the iron-finger motif, both of which are important for the dodecamer assembly. The β-hairpin thumb consists of β15, β16, and a four-residue hairpin loop projecting away from the globular C-lobe (Figure 3A). The β-hairpin thumb mediates the hand-in-hand contact between adjacent dimers (Figure 2B). Phe-265 and Val-266 at the tip of the thumb, together with a nearby Phe-368 from the catalytic domain, contact the same set of residues from the adjacent protomer to form a hydrophobic cluster (Figure 4A). To validate the structural observations, the three residues were each subjected to mutational analyses. TibC F368D mutant was largely insoluble. Both F265D and V266D mutants eluted as a dimer from the gel filtration column (Figure 4B) and failed to catalyze TibA heptosylation in E. coli (Figure 4C). These results highlight the importance of the β-hairpin thumb in mediating the hand-in-hand association for assembly of an active TibC dodecamer. The data also explain the unique presence of two successive hydrophobic residues on TibC β-hairpin loop, which differs from other β-hairpin loops that are usually dominated by hydrophilic residues (Sibanda and Thornton, 1985). The β-hairpin thumb is highly conserved in the BAHT family (Figure 1—figure supplement 1), supporting the view that an intact dodecamer is critical for this family of heptosyltransferase to modify its partner autotransporter.
A unique iron-finger motif and dodecamer assembly
The other insertion module (residues -368) features two knuckles with long exposed loops capping and embellishing the C-lobe (Figure 3A). Within the module, a ferric ion is tetrahedrally coordinated by four conserved cysteine residues (Cys-339, 342, 358, and 370 in TibC) (Figure 3A and Figure 1—figure supplement 1), analogous to the zinc-finger motif. This iron-finger motif lies ∼20 Å from the catalytic cleft (see below) and its architecture belongs to the 1Fe-0S cluster. Differing from the most common zinc-finger motif, the 1Fe-0S iron-finger motif is only observed once in rubredoxin protein found in sulfur-metabolizing bacteria and archaea (Lovenberg and Sobel, 1965; Adman et al., 1975). In rubredoxin, the 1Fe-0S cluster acts as an electron transfer carrier (Sieker et al., 1994). The iron-finger motif in TibC, however, does not seem to function in electron transfer according to our structural and functional analyses. Notably, mutation of any of the four iron-coordinating cysteine residues resulted in a colorless TibC protein (Figure 5A) due to the loss of iron binding. The four cysteine mutants eluted from the gel filtration column as heterogeneous interconverting oligomers with their sizes smaller than a dodecamer (Figure 5B). Thus, the iron finger is required for maintaining structural integrity of the TibC dodecamer. Consistently, the four iron binding-deficient mutants all failed to catalyze TibA heptosylation in the E. coli system (Figure 5C).
Structurally, the iron-finger motif contributes to the back-to-back dimerization. A portion of the iron-finger motif, together with several surface residues from the catalytic domain of one protomer, makes extensive contacts with the β-barrel as well as the catalytic domain of the other protomer in the back-to-back dimer (Figure 5D). A symmetric interface is formed via the same kind of interactions between the two protomers. Extensive polar and hydrophobic contacts occur at the interface (Figure 5D), burying ∼1467 Å2 solvent-accessible surface area of TibC (∼9% of the total surface areas). Dimerization via this large interface provides building blocks for TibC dodecamer assembly.
Structure of the TibC dodecamer in complex with ADP-D,D-heptose
We also crystallized the TibC dodecamer in complex with its natural ligand ADP-D-glycero-β-D-manno-heptose (ADP-D,D-heptose). This was accomplished by soaking the sugar ligand into crystals of the catalytically inactive TibC D110A mutant. Structural determination was achieved by molecular replacement, which revealed extra omit density in the TibC catalytic domain. Placing the sugar into the omit density led to a concomitant loss of Fo–Fc density after the refinement. ADP-D,D-heptose was convincingly modeled into each protomer and a final model of 3.88 Å was obtained (Table 1). The 12 sugar ligands form an array of two parallel circles along the inner surface of the dodecamer ring and face towards the hollow center (Figure 6A). ADP-D,D-heptose is located at the cleft between the two lobes of the catalytic domain (Figure 6B and Figure 3B). The adenine heterocycle fits snugly into a hydrophilic pocket and its N1 atom bears a hydrogen bond contact with Arg-286 of TibC (Figure 6B). The hydrogen bond contact, though weak, is incompatible with a guanine due to its saturated N1, providing a structural basis for the specificity of TibC in using ADP-activated sugar donor. Moreover, the heptose hydroxyl O4 forms a hydrogen bond with Trp-305 indolic N1 (Figure 6B). The β-phosphate has strong interactions with Thr-226 and Lys-230 in the ligand-binding pocket of TibC. Lys-230 is also close to the heptose hydroxyl O5; its positive charge could potentially alleviate the negative charge of β-phosphate, thereby facilitating nucleophilic attack of the heptose C1 by the substrate acceptor.
Superimposition of the TibC sugar binding site with those of the GT-B glycosyltransferases WaaC and MurG identified Asp-110 as the candidate catalytic base that activates the substrate acceptor. Asp-110 is located close to the anomeric hydroxyl of heptose in the TibC–ADP-D,D-heptose complex structure (Figure 6B). All four sugar-binding residues (Arg-286, Trp-305, Thr-226, and Lys-230) as well as the catalytic base (Asp-110) are highly conserved in the BAHT family (Figure 1—figure supplement 1). Alanine substitution of Asp-110, Arg-286, Trp-305, or Lys-230 in TibC all abolished or largely diminished TibA/AIDA-I heptosylation in E. coli (Figure 6—figure supplement 1A); these point mutants also failed to support AIDA-I-mediated bacterial adhesion to HeLa cells (Figure 6—figure supplement 1B). The results also confirm the functional importance of the heptosyltransferase activity of TibC and the BAHT family.
Sugar ligand stereoselectivity
AAH accepts both ADP-D,D-heptose and ADP-L,D-heptose as the sugar donor despite a possible slight preference for the latter (Figure 6C and Figure 6—figure supplement 2). In contrast, TibC exhibited a high stereoselectivity for ADP-D,D-heptose. To reveal the underlying structural basis, we modeled the structure of sugar ligand-bound AAH using the TibC complex as the template and examined TibC/AAH residues within 5 Å from the sugar ligand (Figure 6D). This led to the identification of residue 300 in TibC (294 in AAH) as a possible structural determinant. In TibC, residue 300 is a proline and its main-chain carbonyl oxygen forms two hydrogen bonds with D,D-heptose hydroxyl O4 and O6, respectively (Figure 6E). Notably, the latter interaction is incompatible with L,D-heptose due to the chiral inversion of C6 and the resulting spatial distance. For AAH, the side-chain hydroxyl of Ser-294 can have a polar interaction with hydroxyl O6 in D,D-heptose or hydroxyl O7 in L,D-heptose (Figure 6E), thereby capable of tolerating both sugar ligands. Consistent with these structural indications, TibC P300S mutant could use both ADP-D,D-heptose and ADP-L,D-heptose while AAH S294P mutant instead became selective for ADP-D,D-heptose (Figure 6C and Figure 6—figure supplement 2A). Thus, Pro-300-mediated polar interaction is a critical structural determinant for TibC sugar ligand stereoselectivity.
Cryo-EM analysis of TibC–TibA complex reveals two conformational states
We further succeeded in purifying a TibC–TibA enzyme/substrate complex by co-expression of the sugar binding-deficient TibC K230A mutant with TibA55-350 in E. coli. The complex eluted from the gel filtration column as a homogenous fraction with a size larger than the TibC dodecamer. The molecular weight of the giant enzyme/substrate assembly was determined to be 750 kDa by analytical ultracentrifugation (Figure 2—figure supplement 2), matching well the composition of a TibC dodecamer and six TibA molecules. Cryo-EM micrographs revealed a well-defined shape of the TibC–TibA dodecamer/hexamer complex at an intermediate resolution (Figure 7—figure supplement 1A), allowing for single particle structural analysis (Figure 7—figure supplements 1, 2). Averaging and refinement without classification from the whole dataset of 53,303 particles produced an electron density map of 9.7 Å in which the TibC12–TibA6 stoichiometry became immediately evident (Figure 7A). In this low-resolution model, density agreeing well with the crystal structure of the TibC dodecamer could be clearly identified, but the TibA moiety is slightly blurred and segmented, indicating a conformational heterogeneity. The reconstructed maps could be further classified into two classes based on different features of the TibA density (Figure 7—figure supplement 2). Two maps of higher resolution at 8.2 Å and 8.9 Å were then reconstructed by averaging 35,300 and 18,003 sorted particles, respectively (Figure 7—figure supplement 2). As a result, a higher contrast in the TibA structure was observed between the two models despite a similar overall architecture for the entire complex (Figure 7A). These analyses indicate two major conformational states of the TibC–TibA enzyme/substrate complex in the absence of the sugar donor.
Resting state of the TibC–TibA dodecamer/hexamer complex
The 8.9 Å EM map fits remarkably well with crystal structures of TibA55-350 (PDB ID code: 4Q1Q) (Lu et al., 2014) and the TibC dodecamer (Figure 7B), allowing for reliable atomic model docking. The cross-correlation coefficient (CCC) for TibA in the final model is 0.76 (only a small N-terminal portion of TibA55-350 is out of the density). The reconstruction revealed that six TibA55-350 molecules are distributed along the inner surface of the TibC dodecamer ring and share its sixfold axis (Figure 7C). Importantly, each TibA55-350 binds to two TibC protomers of a back-to-back dimer, suggesting that two enzyme molecules are responsible for catalyzing full heptosylation of one TibA substrate. Ser-176 is the closest sugar acceptor situated at the entrance of the catalytic cleft of one TibC protomer (TibC-A), but the distance to the catalytic base exceeds 20 Å (Figure 7C and Figure 7—figure supplement 3A). A similar distance was observed between the other TibC protomer (TibC-B) and its possible substrate acceptor (Figure 7—figure supplement 3A). Thus, the 8.9 Å TibC–TibA dodecamer/hexamer complex model represents a sugar transfer incompetent state, designated hereafter as the resting state. Moreover, the spiral shaft of TibA β-helix adopts a trend roughly parallel to the sixfold axis of the dodecamer which, together with the patterned solenoid-like surface distribution of heptosylation sites on TibA (Lu et al., 2014), indicates a possible screw propelling-like mechanism for processive heptosylation of TibA by the dodecamer ring.
Active state of the TibC–TibA dodecamer/hexamer complex
The static X-ray structures of TibC dodecamer and TibA55-350 were also docked into the 8.2 Å EM map, but the TibA moiety had a few degrees of bending with a protrusion into the catalytic center of one TibC protomer (Figure 7A,B and Figure 7—figure supplement 3B). The extra protruding density could not be appropriately fit despite extensive attempts at possible solutions. The heterogeneity observed in the cryo-EM image likely reflects the dynamic catalytic motions; the protruding density was tentatively interpreted as the TibA loop bearing the sugar acceptor Ser-176. This interpretation was supported by molecular dynamics (MD) simulation of a resting state TibC2–TibA55-350 trimer. In the 400 ps steered MD simulation, pulling Ser-176 in TibA55-350 towards the heptose led to a proportional increase in the CCC (Figure 7D, E). A subsequent 20 ns unbiased MD simulation produced an ensemble of TibA models that matched the EM density with a best overall CCC value of 0.87 (Figure 7B). The final model of the entire complex was well accommodated by the experimental map (Figure 7B and Figure 7—figure supplement 3B). The protruding loop was assigned to residues 170–180 of TibA, resulting from disruption of a local β-strand within these residues. The loop placed Ser-176 into a position ready for attacking the sugar donor in one TibC protomer (TibC-A) (Figure 7C and Figure 7—figure supplement 3C). This ‘protrusion’ state of TibA was designated as the catalytically competent active state. Notably, the other TibC protomer (TibC-B) in the back-to-back dimer and its possible substrate acceptor appeared to be in the resting state (Figure 7—figure supplement 3C,D), suggesting that only one of the two enzyme molecules in the dimer can be productive in transferring the sugar onto TibA at a time. It is worth mentioning here that the EM dataset presumably covers a continuous series of catalytic stages; the two discrete classifications employed in our reconstruction may not give a complete picture of all catalytic motions of the TibC–TibA enzyme/substrate complex. The active state structure, which was reconstructed from two-thirds of the total particles, probably contains the majority of representative conformational stages in the catalytic process and therefore is used to represent the overall structure of the TibC–TibA dodecamer/hexamer complex (Figure 7C).
Discussion
The large ring-shape platform of TibC dodecamer for efficient heptosylation of TibA
TibC is distinct from known glycosyltransferases in that it adopts an unusual giant dodecamer ring architecture in catalyzing simultaneous glycosylation of six molecules of TibA autotransporter. The selection of a double-layered ring with a large hollow channel is advantageous and provides a favorable environment for cooperative catalysis in a structurally economical manner. Assembly of the dodecamer into a ring can reduce the solvent-accessible surface area and renders a stable enzymatic machinery complex beneficial to efficient catalysis. Moreover, the spatial arrangement of the back-to-back TibC dimer, binding to the different glycosylation sites in one TibA, suggests a divalent catalysis is also probably contributing to the high catalytic efficiency required for hyperglycosylation. Thus, the dodecamer/hexamer enzyme/substrate complex with six symmetric units may represent a best compromise between structural economy and catalytic cooperativity/efficiency.
The TibC complex, compared with other macromolecular assemblies with large interior volumes such as ferritin and multimeric chaperones/proteases, is more permeable with 12 titled grooves extending from outside to the interior chamber. These grooves are optimally imbedded and close to the catalytic center, allowing for easy access of ADP-heptose from outside into the TibC particle. The catalytic center, on the other hand, is structurally best accessible from the interior, where TibA is located and can readily deliver the to-be-glycosylated serine by loop protrusion. We further speculate that the energy used to drive the rotation and shuttle motion of the TibA substrate is partially derived from breakdown of the ADP-heptose glycosidic bond occurring during catalysis. Further experiments such as obtaining a higher resolution cryo-EM structure of the TibC12–TibA6 complex, which may capture more catalytic stages/events, are necessary to validate the proposed model for TibC-catalyzed processive heptosylation of TibA.
The unique iron-finger motif
Another most unusual feature of the TibC protomer structure is the presence of the 1Fe-0S iron-finger motif. The iron-finger motif does not fulfill the previously defined function of serving as an electron carrier in rubredoxin protein. Biochemical and structural studies have underscored the singular importance of this motif in mediating the back-to-back dimer formation and TibA substrate recognition. Interestingly, the iron-finger motif, adopting a topological fold similar to the zinc-finger motif, follows a right-handed helical path in binding to TibA β-helix, which is analogous to DNA binding by zinc-finger motifs. In addition, the recognition specificity and the interaction strength of this motif will have been precisely tuned during macromolecular evolution so that specific substrate recognition, sliding of the TibA β-helix along the reaction channel, and efficient protrusion of the TibA loop can all occur in a catalytically productive manner. It will be interesting to elucidate the detailed mechanism by which iron-finger motif-mediated substrate recognition is effectively coupled with its promoting of the TibA loop protrusion.
Materials and methods
Constructs and recombinant proteins
Request a detailed protocolThe cDNAs encoding TibC and TibA were PCR amplified from genomic DNA of ETEC strain H10407. The full-length tibC was cloned into pGEX-6p-2 vector (GE Healthcare, UK) for recombinant expression in E. coli. To co-express TibA55-350 and TibC K230A for capturing the TibC–TibA55-350 enzyme–substrate complex, the tibA fragment was cloned into the pGEX-6p-2 vector and a catalytically inactive mutant of TibC was inserted into the pACYCDuet vector. AAH and AIDA-I constructs have been described previously (Lu et al., 2014). All constructs were generated by standard molecular cloning. Point mutations were generated by the QuickChang Site-Directed Mutagenesis Kit (Stratagene, La Jolla, CA). All the plasmids were verified by DNA sequencing.
All recombinant proteins and protein complexes were expressed in E. coli BL21 (DE3) Gold strain (Agilent Technologies, Santa Clara, CA) with 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) for overnight induction at 22°C in Luria–Bertani (LB) medium. TibC was purified by glutathione affinity chromatography and the GST tag was removed by PreScission protease digestion. TibC was further purified by anion exchange and gel filtration chromatography using HiTrap Q HP and Superdex 200 Hi-load columns (GE Healthcare), respectively. The TibC–TibA55-350 complex was initially purified using one-step glutathione affinity chromatography. The GST moiety on TibA was removed by PreScission protease cleavage. The complex was further separated from free TibA55-350 by gel filtration chromatography using the Superdex 200 Hi-Load column. SeMet-labeled TibC and TibA were expressed in the methionine-auxotrophic E. coli B834 (DE3) strain (Novagen, Germany) using SeMet-supplemented M9 medium. SeMet protein was purified following the same scheme as that used for the native protein.
Crystallization and data collection
Request a detailed protocolAll the crystallization experiments were carried out at 20°C using the hanging drop vapor diffusion method. The purified TibC was concentrated to 20 mg/ml in a buffer containing 10 mM Tris–HCl (pH 7.6), 100 mM NaCl, and 2 mM DTT. Crystals of both native and SeMet-TibC protein were grown in 8% (wt/vol) PEG 8000, 120 mM magnesium acetate, and 100 mM MES buffer (pH 5.5) for two days. Crystals were cryoprotected with the mother liquid containing 19% (vol/vol) ethylene glycol and 1% (vol/vol) DMSO and flash-frozen with liquid nitrogen for data collection. The TibC D110A mutant crystals were grown at the same condition. Before freezing, the crystals were soaked with 5 mM ADP-D,D-heptose in the mother liquid for 1 hr.
Structure determination
Request a detailed protocolAll the diffraction data were collected at BL-17U of Shanghai Synchrotron Radiation Facility (SSRF) at 100 K and processed with the HKL 2000 suite by the routine procedure (Otwinowski and Minor, 1997). The structure of TibC dodecamer was solved by the single wavelength anomalous dispersion (SAD) method. In the TibC data, a total of 36 heavy atoms were determined using the program SHELXD (Schneider and Sheldrick, 2002). The identified heavy atom sites were then refined and the initial phases were calculated using the program autoSHARP (Global Phasing Ltd). The phase was then refined and extended to 2.9 Å by NCS averaging, solvent flattening and histogram matching using the DM module in CCP4 (Dodson et al., 1997). After density modification, the map was sufficiently interpretable for model building. The model was automatically built using the Buccaneer software (Cowtan, 2006), which generated ∼70% complete model. The remaining parts were manually built and adjusted in Coot (Emsley et al., 2010). The model was refined in Refmac (Murshudov et al., 1997). The TibC–ADP-D,D-heptose complex structure was solved by molecular replacement in Phaser MR (McCoy et al., 2007) using the apo-TibC structure as the template. The final model was adjusted in Coot and refined in Refmac. The structural pictures were generated in Pymol (www.pymol.org).
Gel filtration chromatography analysis of the protein size
Request a detailed protocolPurified recombinant TibC (WT or indicated mutants) were loaded onto a HiLoad 16/600 Superdex 200 column (GE Healthcare) in a buffer containing 20 mM Tris–HCl (pH 7.6) and 20 mM NaCl on an ÄKTA Purifier (GE Healthcare). Fractions corresponding to the specified elution volumes were collected and loaded onto SDS-PAGE gels followed by Coomassie blue staining.
In vivo and in vitro heptosyltransferase enzymatic assays
Request a detailed protocolIn vivo heptosyltransferase activity of TibC was monitored by co-expressing pACYCDuet-Flag-TibA (the Flag tag was inserted at the junction between the signal peptide and the passenger domain) and pET21a-TibC (WT and indicated mutants). 1/10 of 100 μl of total lysates from 1 mL of E. coli culture was used for immunoblotting and glycoprotein detection analyses. Expression levels of TibA and TibC were detected by anti-Flag M2 mouse monoclonal antibody (Sigma-Aldrich, St. Louis, MO) and in-house generated rabbit anti-TibC serum, respectively. Glycosylated proteins were detected with the ECL glycoprotein detection module (GE Healthcare, product code RPN2190) according to the manufacturer’s protocol. For in vitro analysis of TibC and AAH heptosyltransferase activities, 10 μg of the TibA-derived synthetic peptide (GGVQHVSSGGSATSSTINSGG, synthesized by Scilight Biotechnology LLC, China) were incubated with 2 μg of purified TibC or AAH protein (WT or indicated mutants) for 4 hr at 30°C in a 20 μl reaction containing 10 mM Tris–HCl (pH 7.5), 50 mM NaCl, 10 mM MgCl2, 1 mM DTT, and 0.5 mM ADP-L,D-heptose or ADP-D,D-heptose ligands (Zamyatina et al., 2003). The reaction mixtures were subjected to mass spectrometry analysis.
Analytic ultracentrifugation
Request a detailed protocolAnalytic ultracentrifugation was carried out with an XL-I analytical ultracentrifuge (Beckman-Coulter, Indianapolis, IN) with An-60 Tirotor at 4°C. The TibC and TibC–TibA complex proteins were prepared at 0.8 mg/ml in a buffer containing 20 mM Tris–HCl (pH 8.0), 150 mM NaCl, and 2 mM DTT. All data were collected at a speed of 25,000/32,000 rpm. 3 mM path length charcoal-filled Epon centerpieces were used for the TibC–TibA complex due to the sample quality limitation and a 12 mM path length aluminum centerpiece was used for the TibC sample. Samples were monitored real time using interference optics at intervals of 4 min. Buffer alone was used as the reference. The molecular weights were calculated by the SEDFIT program (Schuck, 2000).
Cryo-EM data acquisition, image processing, and reconstruction
Request a detailed protocol10 μg/ml bacitracin (Sigma) was added to the purified protein to obtain mono-dispersed particles and make the orientation distribution more anisotropic. An aliquot of 3.5 μl of TibC or TibC–TibA55-350 complex (∼1 mg/ml) was absorbed onto a glow-discharged Quantifoil holey carbon grid (R2.1, 300 mesh; Jena Biosciences, Germany). The grid was blotted in an FEI Vitrobot Mark IV (FEI Company, Hilsboro, OR) using 4 s blotting time and blotting force 2 with 100% humidity at 283 K, and then plunged into the chilled liquid ethane. The grids were immediately transferred into a Titan Krios electron microscope (FEI Company) operated at 300 kV. For the TibC dodecamer, approximately 3000 images were semi-automatically collected following the Leginon procedure (Carragher et al., 2000) using a 4 K × 4 K Gatan UltraScan4000 CCD (model 895, Gatan Inc, Pleasanton, CA) at a nominal magnification of 75,000×. For the TibC–TibA55-350 complex, ∼1800 cryo-EM images were recorded using a 2K × 2K Gatan Tridiem CCD camera (Gatan Inc) at a nominal magnification of 81,000× with the energy filter turned off. The corresponding pixel size on the specimen level is 1.196 and 1.778 Å/pixel for TibC and TibC–TibA55-350 complex datasets, respectively. The range of nominal defocus is 2.3–3.6 μm for TibC and 2.7–4.1 μm for TibC–TibA55-350 data. The electron dose is 20 e−/Å·s for each CCD frame. For both TibC and TibC–TibA55-350 micrographs, the particles were semi-automatically boxed using e2boxer.py in EMAN2 package (Tang et al., 2007) followed by interactively screening. The particles with the wrong shape and size were discarded. Approximately 10,200 and 56,000 particle images were cropped out and normalized (mean = 0, standard deviation = 1) for TibC and TibC–TibA55-350 complex, respectively. After determining the defocus and B-factor, the particles were phase-flipped for contrast transfer function (CTF) correction using e2ctf.py in EMAN2 (Tang et al., 2007).
For the TibC dataset, two-dimensional multivariate statistical analysis-based reference-free classification was performed in EMAN2 to further eliminate particles belonging to the class averages with blurry appearance. Particles contained in the averages, which were redundant in the ring-like view, were also removed to avoid introducing orientation isotropy. A total of 8546 particles from the TibC data were used for subsequent data processing. The initial model was built using e2initialmodel.py with C1 symmetry applied and was iteratively refined against the phase-flipped particles until no further improvement could be obtained. During the projection-matching refinement process, a Fourier ring correlation (FRC) was chosen to compare the similarity between the class averages and re-projections of the model, and C6 symmetry was enforced as indicated by the reference-free class averages. 20% of the worst particles in each class were discarded in each round of iteration. The information over a 20–100 Å resolution range and SSNR weighting were considered in the comparison in several early iterations. In the last four rounds of iterations, the information was expanded to the full resolution range and the step angle was assigned to 2.0°. The map obtained in the last iteration was then subjected to B-factor correction with a damping factor of −210 Å2 to enhance the high resolution information using the program embfactor (Fernandez et al., 2008).
For the TibC–TibA55-350 complex data, the processing procedure before reconstruction was similar to that described for the TibC data above. A total of 53,303 particles were finally employed for 3D reconstruction. The starting model used was a 60 Å low-pass filtered TibC dodecamer reconstruction (C6 symmetry applied). After several rounds of refinement running, the phase-flipped particle images were imported into RELION (Scheres, 2012) for 3D classification. After 30 rounds of iterations with three-class classification, two classes that contained 35,300 members in total and showed almost identical reconstruction in the TibA density were merged together for subsequent 3D refinement. Finally, two clusters that corresponded to the ‘active’ and ‘resting’ conformations of the complex were acquired. After the 3D refinement converged, the particles corresponding to the two clusters were sorted for subsequent single-model projection-matching refinement in EMAN2. Final maps for the two conformations were independently reconstructed through projection-matching refinement with the step angle of 1.68° for the active conformation and 3.6° for the resting conformation. Resolution was estimated by the gold standard FSC using the 0.143 cutoff value (Rosenthal and Henderson, 2003; Scheres and Chen, 2012). An additional projection-matching reconstruction test for potential model bias was performed with the resting volume as the initial reference against the active state particle set and vice versa. The structural pictures of cryo-EM reconstruction were drawn in UCSF chimera (Pettersen et al., 2004).
Molecular modeling of AAH structure
Request a detailed protocolHomology modeling of AAH was performed in the program MODELLER (version 9v7) (Marti-Renom et al., 2000; Bredenberg and Nilsson, 2001) using the TibC structure as the template. The quality of the modeled structure was ensured by the high sequence identity between AAH and TibC (68%) (Figure 1—figure supplement 1). The top DPOE scored model was chosen and validated with PROCHECK (Laskowski et al., 1993). Structural models of AAH S294P and TibC P300S mutants were generated directly from the corresponding wild-type structures. Residues within a radius of 5 Å from the mutated residue were refined using the Protein Local Optimization Program (PLOP) (Jacobson et al., 2002; Zhu et al., 2007).
Refinement of TibC12–TibA6 complex structures
Request a detailed protocolMolecular dynamics flexible fitting (MDFF) (Trabuco et al., 2009) was used to refine the complex structure by incorporating the EM density map as an external potential into MD sampling. All MD simulations were performed in the program NAMD 2.9 (Phillips et al., 2005) using the CHARMM27 force field (MacKerell et al., 1998) including the CMAP correction (Mackerell et al., 2004). The MDFF simulations were carried out at T = 300 K with a scaling factor ζ = 1 kcal for 10 ns. The snapshot with the top CCC was selected for further unbiased MD simulation. The value of CCC was calculated using the program VMD (version 1.9.1) (Humphrey et al., 1996).
Molecular dynamics simulation
Request a detailed protocolThe CHARMM force field files (topology and parameter) for ADP-D,D-heptose were automatically generated using the server SwissParam (Zoete et al., 2011), except that the partial atomic charges of ADP-D,D-heptose were manually assigned using the analogs in the existing CHARMM force field (Best et al., 2012). The parameters of ferric ions with four coordinating cysteine residues were modified according to the available parameters developed for the zinc-finger motif (Eriksson et al., 1995; Bredenberg and Nilsson, 2001).
To sample the intermediate state of loop protrusion in the TibC–TibA complex, steered molecular dynamics (SMD) simulation (Isralewitz et al., 2001) was carried out with one unit of TibC12–TibA6 complex. By fixing the position of the C1 atom of ADP-D,D-heptose, a force was applied to the OG atom of Ser-176 along a vector connecting the two atoms. The constant pulling velocity was set to 0.06 Å/ps with a spring constant of 10 kcal/mol/Å2. The trajectories of TibA were extracted from a total of 400 ps SMD simulation and then fitted into the cryo-EM density of TibA, which was extracted from the total map using ‘Volume Eraser’ in CHIMERA (Pettersen et al., 2004). The correlation between the CCC values and the distances of loop protrusion was analyzed. The last snapshot of the TibC–TibA subunit from SMD with minimum distance between C1 and OG was chosen as the initial model of the active state conformation, and a 20 ns unbiased MD simulation was then carried out for further refinement. The MD snapshot with the top CCC score was used to generate the final active state TibC–TibA dodecamer/hexamer complex structure.
The MDFF-refined TibC12–TibA6 complex structure was solvated using TIP3P water molecules and then neutralized. The system was minimized for 10,000 integration steps and equilibrated for 100 ps with a time step of 1 fs, and the temperature was gradually increased from 25 K to 300 K. Following this, an unbiased MD simulation was performed under a constant temperature of 300 K and a constant pressure of 1 atm using the Nosé–Hoover Langevin piston method (Feller et al., 1995) with the integration time step of 2 fs. The cutoff distances for electrostatic and van der Waals calculations were set at 10 Å. The long-range electrostatic forces were computed using the particle mesh Ewald method (Darden et al., 1993) with a grid spacing of 1 Å. All covalent bonds involving hydrogen atoms were constrained with the SHAKE algorithm.
Accession numbers
Request a detailed protocolCrystal structural data for apo-TibC and TibC/ADP-heptose complex are deposited in the Protein Data Bank (PDB) under the accession numbers 4RAP and 4RB4, respectively. Cryo-EM structure data are deposited in Electron Microscopy Data Bank (EMD) with accession codes EMD-2755 (apo-TibC), EMD-2756 (averaged TibC12–TibA6 complex), EMD-2757 (active state TibC12–TibA6 complex), and EMD-2758 (resting state TibC12–TibA6 complex).
Data availability
-
Crystal structure of bacterial iron-containing dodecameric glycosyltransferase TibC from enterotoxigenic E. coli H10407Publicly available at RCSB Protein Data Bank.
-
Cryo-electron microscopy of TibC dodecamerPublicly available at Electron Microscopy Data Bank.
-
Cryo-electron microscopy of TibC12-TibA6 octadecamer in averaged conformationPublicly available at Electron Microscopy Data Bank.
-
Cryo-electron microscopy of TibC12-TibA6 octadecamer in active statePublicly available at Electron Microscopy Data Bank.
-
Cryo-electron microscopy of TibC12-TibA6 octadecamer in resting statePublicly available at Electron Microscopy Data Bank.
-
Crystal structure of TibC-catalyzed hyper-glycosylated TibA55-350 fragmentPublicly available at RCSB Protein Data Bank.
-
E. coli heptosyltransferase WaaC with ADP-2-deoxy-2-fluoro heptosePublicly available at RCSB Protein Data Bank.
-
Crystal structure of MurG:GlcNAc complexPublicly available at RCSB Protein Data Bank.
References
-
Not just for Eukarya anymore: protein glycosylation in Bacteria and ArchaeaCurrent Opinion in Structural Biology 18:544–550.https://doi.org/10.1016/j.sbi.2008.06.010
-
NH–S hydrogen bonds in Peptococcus aerogenes ferredoxin, Clostridium pasteurianum rubredoxin, and Chromatium high potential iron proteinProceedings of the National Academy of Sciences of USA 72:4854–4858.https://doi.org/10.1073/pnas.72.12.4854
-
Cloning and expression of an adhesin (AIDA-I) involved in diffuse adherence of enteropathogenic Escherichia coliInfection and Immunity 57:1506–1511.
-
Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone phi, psi and side-chain chi(1) and chi(2) dihedral anglesJournal of Chemical Theory and Computation 8:3257–3273.https://doi.org/10.1021/ct300400x
-
Glycoside hydrolases and glycosyltransferases: families and functional modulesCurrent Opinion in Structural Biology 11:593–600.https://doi.org/10.1016/S0959-440X(00)00253-0
-
Modeling zinc sulfhydryl bonds in zinc fingersInternational Journal of Quantum Chemistry 83:230–244.https://doi.org/10.1002/qua.1214
-
Recent structures, evolution and mechanisms of glycosyltransferasesCurrent Opinion in Structural Biology 22:540–549.https://doi.org/10.1016/j.sbi.2012.06.007
-
Structures and mechanisms of glycosyltransferasesGlycobiology 16:29R–37R.https://doi.org/10.1093/glycob/cwj016
-
Leginon: an automated system for acquisition of images from vitreous ice specimensJournal of Structural Biology 132:33–45.https://doi.org/10.1006/jsbi.2000.4314
-
Functional organization of the autotransporter adhesin involved in diffuse adherenceJournal of Bacteriology 189:9020–9029.https://doi.org/10.1128/JB.01238-07
-
The Buccaneer software for automated model building. 1. Tracing protein chainsActa Crystallographica Section D, Biological Crystallography 62:1002–1011.https://doi.org/10.1107/S0907444906022116
-
Particle mesh Ewald: an N • log (N) method for Ewald sums in large systemsJournal of Chemical Physics 98:10089.https://doi.org/10.1063/1.464397
-
Epithelial cell invasion and adherence directed by the enterotoxigenic Escherichia coli tib locus is associated with a 104-kilodalton outer membrane proteinInfection and Immunity 62:3463–3471.
-
Features and development of CootActa Crystallographica Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Constant pressure molecular dynamics simulation: the Langevin piston methodJournal of Chemical Physics 103:4613.https://doi.org/10.1063/1.470648
-
Sharpening high resolution information in single particle electron cryomicroscopyJournal of Structural Biology 164:170–175.https://doi.org/10.1016/j.jsb.2008.05.010
-
Structure of the Escherichia coli heptosyltransferase WaaC: binary complexes with ADP and ADP-2-deoxy-2-fluoro heptoseJournal of Molecular Biology 363:383–394.https://doi.org/10.1016/j.jmb.2006.07.057
-
Type V protein secretion pathway: the autotransporter storyMicrobiology and Molecular Biology Reviews 68:692–744.https://doi.org/10.1128/MMBR.68.4.692-744.2004
-
Crystal structure of the MurG: UDP-GlcNAc complex reveals common structural principles of a superfamily of glycosyltransferasesProceedings of the National Academy of Sciences of USA 100:845–849.https://doi.org/10.1073/pnas.0235749100
-
Remarkable structural similarities between diverse glycosyltransferasesChemistry & Biology 9:1287–1296.https://doi.org/10.1016/S1074-5521(02)00295-8
-
VMD: visual molecular dynamicsJournal of Molecular Graphics 14:33–38.https://doi.org/10.1016/0263-7855(96)00018-5
-
Steered molecular dynamics and mechanical functions of proteinsCurrent Opinion in Structural Biology 11:224–230.https://doi.org/10.1016/S0959-440X(00)00194-9
-
Force field validation using protein side chain predictionJournal of Physical Chemistry B 106:11673–11680.https://doi.org/10.1021/jp021564n
-
Glycosyltransferases: structures, functions, and mechanismsAnnual Review of Biochemistry 77:521–555.https://doi.org/10.1146/annurev.biochem.76.061005.092322
-
PROCHECK: a program to check the stereochemical quality of protein structuresJournal of Applied Crystallography 26:283–291.https://doi.org/10.1107/S0021889892009944
-
Identification of a glycoprotein produced by enterotoxigenic Escherichia coliInfection and Immunity 67:4084–4091.
-
Rubredoxin: a new electron transfer protein from Clostridium pasteurianumProceedings of the National Academy of Sciences of USA 54:193–199.https://doi.org/10.1073/pnas.54.1.193
-
All-atom empirical potential for molecular modeling and dynamics studies of proteinsJournal of Physical Chemistry B 102:3586–3616.https://doi.org/10.1021/jp973084f
-
Comparative protein structure modeling of genes and genomesAnnual Review of Biophysics and Biomolecular Structure 29:291–325.https://doi.org/10.1146/annurev.biophys.29.1.291
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
Refinement of macromolecular structures by the maximum-likelihood methodActa Crystallographica Section D, Biological Crystallography 53:240–255.https://doi.org/10.1107/S0907444996012255
-
Protein glycosylation in bacteria: sweeter than everNature Reviews Microbiology 8:765–778.https://doi.org/10.1038/nrmicro2383
-
BookProcessing of x-ray diffraction data collected in oscillation modeOtwinowski Z, Minor W, editors. New York: Academic Press.
-
UCSF Chimera–a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Scalable molecular dynamics with NAMDJournal of Computational Chemistry 26:1781–1802.https://doi.org/10.1002/jcc.20289
-
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
-
Prevention of overfitting in cryo-EM structure determinationNature Methods 9:853–854.https://doi.org/10.1038/nmeth.2115
-
Substructure solution with SHELXDActa Crystallographica Section D, Biological Crystallography 58:1772–1779.https://doi.org/10.1107/S0907444902011678
-
Glycosylation of the self-recognizing Escherichia coli Ag43 autotransporter proteinJournal of Bacteriology 188:1798–1807.https://doi.org/10.1128/JB.188.5.1798-1807.2006
-
Rubredoxin in crystalline stateMethods in Enzymology 243:203–216.https://doi.org/10.1016/0076-6879(94)43016-0
-
Protein glycosylation in bacterial mucosal pathogensNature Reviews Microbiology 3:225–237.https://doi.org/10.1038/nrmicro1100
-
EMAN2: an extensible image processing suite for electron microscopyJournal of Structural Biology 157:38–46.https://doi.org/10.1016/j.jsb.2006.05.009
-
Glycosyltransferase structure and mechanismCurrent Opinion in Structural Biology 10:510–517.https://doi.org/10.1016/S0959-440X(00)00124-X
-
Autotransporter proteins: novel targets at the bacterial cell surfaceFEMS Microbiology Letters 274:163–172.https://doi.org/10.1111/j.1574-6968.2007.00833.x
-
Homology between O-linked GlcNAc transferases and proteins of the glycogen phosphorylase superfamilyJournal of Molecular Biology 314:365–374.https://doi.org/10.1006/jmbi.2001.5151
-
Improved methods for side chain and loop predictions via the protein local optimization program: variable dielectric model for implicitly improving the treatment of polarization effectsJournal of Chemical Theory and Computation 3:2108–2119.https://doi.org/10.1021/ct700166f
-
SwissParam: a fast force field generation tool for small organic moleculesJournal of Computational Chemistry 32:2359–2368.https://doi.org/10.1002/jcc.21816
Article and author information
Author details
Funding
Howard Hughes Medical Institute (55007431)
- Feng Shao
Ministry of Science and Technology of the People's Republic of China (2012CB518700 and 2014CB849602)
- Qing Yao
- Qiuhe Lu
- Yue Xu
- Feng Shao
National Natural Science Foundation of China (Program for Distinguished Young Scholars 31225002)
- Feng Shao
Chinese Academy of Sciences (XDB08020202)
- Feng Shao
National Natural Science Foundation of China (Program for Distinguished Young Scholars 31230018)
- Ping Zhu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the staff at Shanghai Synchrotron Radiation Facility (SSRF) for their assistance in data collection. We thank W Chu at Tsinghua University and the China National Center for Protein Sciences Beijing for performing analytical ultracentrifugation. We thank Dr. Rong Meng and the mass spectrometry facility of National Center for Protein Sciences at Peking University for assistance with ETD mass spectrometry analyses. Computational support was provided by the National Supercomputer Center in Tianjin, and the calculations were performed on TianHe-1(A). We are also grateful to all members of the Shao laboratory for discussions and technical assistance. The research was supported in part by an International Early Career Scientist grant from the Howard Hughes Medical Institute to FS. This work was also supported by the National Basic Research Program of China 973 Programs (2012CB518700 and 2014CB849602), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB08020202) to FS, and the China National Science Foundation Program for Distinguished Young Scholars to FS (31225002) and PZ (31230018).
Copyright
© 2014, Yao 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
-
- 3,744
- views
-
- 413
- downloads
-
- 29
- 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
Many voltage-gated potassium (Kv) channels display a time-dependent phenomenon called C-type inactivation, whereby prolonged activation by voltage leads to the inhibition of ionic conduction, a process that involves a conformational change at the selectivity filter toward a non-conductive state. Recently, a high-resolution structure of a strongly inactivated triple-mutant channel kv1.2-kv2.1-3m revealed a novel conformation of the selectivity filter that is dilated at its outer end, distinct from the well-characterized conductive state. While the experimental structure was interpreted as the elusive non-conductive state, our molecular dynamics simulations and electrophysiological measurements show that the dilated filter of kv1.2-kv2.1-3m is conductive and, as such, cannot completely account for the inactivation of the channel observed in the structural experiments. The simulation shows that an additional conformational change, implicating isoleucine residues at position 398 along the pore lining segment S6, is required to effectively block ion conduction. The I398 residues from the four subunits act as a state-dependent hydrophobic gate located immediately beneath the selectivity filter. By mutating I398 to Asparagine, ion permeation can be resumed in the kv1.2-kv2.1-3m channel, which was not a reversion C-type inactivation, since AgTxII fails to block the ionic permeation of kv1.2-kv2.1-3m_I398N. As a critical piece of the C-type inactivation machinery, this structural feature is the potential target of a broad class of QA blockers and negatively charged activators thus opening new research directions towards the development of drugs that specifically modulate gating-states of Kv channels.
-
- Structural Biology and Molecular Biophysics
Kainate receptors are key modulators of synaptic transmission and plasticity in the central nervous system. Different kainate receptor isoforms with distinct spatiotemporal expressions have been identified in the brain. The GluK1-1 splice variant receptors, which are abundant in the adult brain, have an extra fifteen amino acids inserted in the amino-terminal domain (ATD) of the receptor resulting from alternative splicing of exon 9. However, the functional implications of this post-transcriptional modification are not yet clear. We employed a multi-pronged approach using cryogenic electron microscopy, electrophysiology, and other biophysical and biochemical tools to understand the structural and functional impact of this splice insert in the extracellular domain of GluK1 receptors. Our study reveals that the splice insert alters the key gating properties of GluK1 receptors and their modulation by the cognate auxiliary Neuropilin and tolloid-like (Neto) proteins 1 and 2. Mutational analysis identified the role of crucial splice residues that influence receptor properties and their modulation. Furthermore, the cryoEM structure of the variant shows that the presence of exon 9 in GluK1 does not affect the receptor architecture or domain arrangement in the desensitized state. Our study thus provides the first detailed structural and functional characterization of GluK1-1a receptors, highlighting the role of the splice insert in modulating receptor properties and their modulation.