Abstract
Free full text
Crystal structure of the Helicobacter pylori vacuolating toxin p55 domain
Abstract
Helicobacter pylori VacA, a pore-forming toxin secreted by an autotransporter pathway, causes multiple alterations in human cells, contributes to the pathogenesis of peptic ulcer disease and gastric cancer, and is a candidate antigen for inclusion in an H. pylori vaccine. Here, we present a 2.4-Å crystal structure of the VacA p55 domain, which has an important role in mediating VacA binding to host cells. The structure is predominantly a right-handed parallel β-helix, a feature that is characteristic of autotransporter passenger domains but unique among known bacterial protein toxins. Notable features of VacA p55 include disruptions in β-sheet contacts that result in five β-helix subdomains and a C-terminal domain that contains a disulfide bond. Analysis of VacA protein sequences from unrelated H. pylori strains, including m1 and m2 forms of VacA, allows us to identify structural features of the VacA surface that may be important for interactions with host receptors. Docking of the p55 structure into a 19-Å cryo-EM map of a VacA dodecamer allows us to propose a model for how VacA monomers assemble into oligomeric structures capable of membrane channel formation.
Helicobacter pylori is a Gram-negative bacterium that chronically infects the stomachs of >50% of the human population. H. pylori infection is a significant risk factor for the development of peptic ulcer disease, gastric adenocarcinoma, and gastric lymphoma (1, 2). One of the important virulence factors expressed by this organism is VacA, a secreted toxin named for its capacity to induce extensive vacuolation in the cytoplasm of mammalian cells (3–5). VacA is translated as a 140-kDa protoxin that undergoes N- and C-terminal cleavage during the secretion process to yield a mature 88-kDa toxin, p88 (6–8). Secretion is thought to occur through a Type Va or autotransporter pathway in which an N-terminal signal sequence directs the protein to the periplasm, and a C-terminal β-barrel domain facilitates transport of the p88 “passenger” domain across the outer membrane (8, 9) (Fig. 1a). VacA is able to intoxicate multiple types of human cells to yield a variety of cellular effects (4, 5). In addition to causing cell vacuolation, VacA can cause depolarization of membrane potential, alteration of mitochondrial membrane permeability, apoptosis, detachment of cells from the basement membrane, activation of mitogen-activated protein kinases, inhibition of antigen presentation, and inhibition of T cell activation and proliferation (4, 5). Although the mechanisms by which these processes occur are not fully understood, many of these effects depend on the capacity of VacA to form anion-selective channels in host-cell membranes (4, 5).
The first step in host-cell intoxication is binding of the toxin to the cell surface. Multiple cell surface receptors for VacA have been identified, including RPTPα, RPTPβ, and various lipids (10–13). VacA can then insert into the plasma membrane to form channels (14) or can undergo endocytosis and eventually localize with late endocytic compartments or mitochondria (15, 16). VacA-induced vacuoles correspond to swollen late endosomes, which are presumed to arise as a consequence of anion flux through VacA channels in the membranes of these compartments (5, 17).
VacA oligomerization is thought to precede membrane-channel formation (18). VacA can assemble into a variety of water-soluble oligomeric structures, including single-layered hexamers and heptamers and double-layered structures consisting of 12 or 14 subunits (19–24). The structures resemble “flowers” in which a central ring is surrounded by peripheral “petals.” Although the highest-resolution images are 19-Å cryo-EM maps of VacA dodecamers (24), atomic-force microscopy and electrophysiological studies suggest that membrane-associated VacA channels are single-layered structures (13, 22).
Two domains of VacA, p33 and p55, have been identified based on partial proteolysis of p88 into fragments of 33 and 55 kDa, respectively (7) (Fig. 1a). When expressed independently and then mixed, p33 and p55 can physically interact and reconstitute vacuolating toxin activity (25–27). The N-terminal p33 domain (residues 1–311) contains a hydrophobic sequence (residues 6–27) involved in pore formation (18, 28), whereas the p55 domain (residues 312–821) contains one or more cell-binding domains (29–31). When expressed intracellularly, the minimum portion of VacA required for cell-vacuolating activity comprises the entire p33 domain and ≈110 aa from the N terminus of p55 (25).
Analysis of H. pylori strains isolated from unrelated humans indicates a high level of genetic diversity among vacA alleles (32–35). Although frequent recombination events have eliminated phylogenetic structure from the 5′ region of vacA alleles (32, 33), phylogenetic analysis of sequences from a vacA midregion (located within p55, Fig. 1a) indicates the existence of two large families of sequences, termed m1 and m2 (34, 35). Within this ≈281-aa midregion (roughly corresponding to amino acids D455 to V735 in the VacA sequence of H. pylori strain 60190), the amino acid sequences of types m1 and m2 VacA proteins are only ≈55% identical. Differences in cell-type specificity have been noted for types m1 and m2 VacA proteins, a phenomenon that may result from binding of m1 and m2 VacA proteins to different cell-surface receptors (36–39). With the rare exception of a few m1/m2 chimeras, there has been little evidence of recombination between m1 and m2 vacA alleles within the vacA midregion, and, therefore, a phylogenetic distinction between m1 and m2 VacA sequences has remained intact (35). This distinction is the basis for a widely used typing scheme for H. pylori isolates (34, 35) and is relevant clinically because m1 strains are associated with gastric cancer (40). Despite a wealth of sequence information, little is known about the structural and functional correlates of VacA sequence diversity. To gain insights into the structural features of VacA that contribute to its cell-binding properties and to better understand the structural correlates of VacA sequence variation, we set out to determine the structure of the VacA p55 domain.
Results and Discussion
Crystal Structure of the VacA p55 Domain.
The crystal structure of the VacA p55 domain was determined to a minimum Bragg spacing of dmin = 2.4 Å by using experimental phases from multiple isomorphous replacement and anomalous scattering [see Methods and supporting information (SI) Table 1]. The p55 structure is predominantly a right-handed parallel β-helix (residues 355–735) but has a small globular domain at the C terminus (residues 736–811) with mixed α/β secondary structure elements (Fig. 1b). The structure resembles a sock in which the C-terminal domain curves from the heel and extends to the tip of the foot. The β-helical “calf” is 65 Å long with widths of 25–31 Å, whereas the C-terminal “foot” is ≈17 × 24 × 43 Å. In the crystal, pairs of p55 molecules meet at their N-terminal strands about a crystallographic 2-fold axis (SI Fig. 5). This packing is consistent with the elution of p55 as a dimer from gel filtration columns (see Methods) and a previous EM study indicating that a p55 domain secreted by H. pylori formed dimers (29).
Comparison of the VacA p55 Structure with Structures of Other Autotransporter Passenger Domains.
Crystal structures have been determined for two other type Va autotransporter passenger domains: pertactin (an adhesin from Bordetella pertussis) and Hbp (a hemoglobin protease from Escherichia coli). Pertactin (41), Hbp (42), and VacA p55 do not share sequence similarities, but all contain a β-helix fold. The β-helix fold is composed of multiple ≈25-aa quasirepeats, each of which forms a single coil of the helix. Conservation of the β-helix fold among autotransporter passenger domains suggests that this structural feature is required for efficient secretion across the outer membrane and folding (43). A C-terminal β-helix cap may be important as a nucleus and/or chaperone for secretion and folding (44). The end of the VacA p55 β-helix is capped by a β-hairpin, similar to that observed in the structure of pertactin (41). There is also a short α-helix in this region of VacA that appears to be unique among β-helix structures (Fig. 1c).
A notable difference between VacA p55 and the two other autotransporter passenger domain structures is the presence of a disulfide in the p55 C-terminal subdomain (Fig. 1c). The low cysteine content observed in autotransporter sequences is consistent with a model in which passenger domains translocate across the outer membrane in an unfolded conformation. Nevertheless, a number of autotransporter sequences have a single, closely spaced pair of cysteines near the C terminus of their passenger domains (45). In VacA, these cysteines are either 11 or 13 aa apart. Mutation of either of these cysteines to serine results in a decrease in toxin secretion but has no effect on the vacuolating activity of the toxin (ref. 45 and M.S.M., unpublished results). Here, we show experimentally that these cysteines form a disulfide and that they are positioned to buttress both the β-helix cap and a conserved C-terminal pocket (discussed below). Another feature of p55 that differs from other autotransporter passenger domains is the presence of multiple kinks that disrupt what would otherwise be continuous β-sheets. We have divided the VacA p55 β-helix into five subdomains to reflect these disruptions (Fig. 2a) and note that the divisions correlate with predicted sites of homologous recombination (as discussed further below).
Sequence Variation Among VacA Proteins.
In an effort to understand how genetic differences between m1 and m2 vacA alleles relate to structural and functional differences in the proteins they encode, we have examined VacA sequence polymorphisms in the context of the p55 structure. We identified sequences for 62 m1 and 27 m2 VacA proteins in GenBank that were complete in the p55 region. In addition, we identified three m1/m2 chimeric VacA proteins for which sites of recombination between m1 and m2 sequences were easily recognizable. We aligned these sequences in various combinations using the program ClustalW (46) and mapped the sequence similarity scores to the three-dimensional m1 p55 structure using ESPript (47). The majority of residues pointing into the interior of the β-helix are either strictly or highly conserved, consistent with the idea that mutation of buried residues would result in a detrimental loss of structure and, therefore, function. The surface-exposed residues are also fairly conserved when exclusively m1 or exclusively m2 sequences are analyzed (SI Fig. 6). However, when m1 and m2 sequences are analyzed together, the surface-exposed residues are highly variable (Fig. 3). There are only two regions of the surface with notable sequence conservation. One is located at the N terminus (Fig. 3a) and will be discussed in the next section with respect to protein oligomerization. The second conserved surface is located at the C terminus of the β-helix in a cavity formed by two loops (residues 668–678 and 730–734) and the disulfide-linked C-terminal domain (Fig. 3b). The strict sequence conservation, the fact that many binding sites are located in clefts or cavities, and the fact that this is the only cavity observed in the VacA p55 structure suggest that this area may represent a receptor-binding site that is shared by m1 and m2 forms of the toxin.
Because multiple surface-exposed regions are highly divergent, it is difficult to identify a single region in the p55 structure that can explain the different binding properties of m1 and m2 VacA proteins (Fig. 3). The biggest difference between m1 and m2 sequences is the presence of a 23-aa insertion in m2 sequences (Fig. 2a, asterisk; and SI Fig. 7). The insertion is notable in that its sequence is an imperfect repeat of the sequence located immediately upstream, and the length approaches the average length of a single β-helix coil (≈25 aa) (SI Fig. 7). The 23-aa m2 insertion could therefore result in an additional β-helix coil similar in structure to the coil that precedes it. Based on an analysis of engineered m1/m2 chimeras, a region important for m1-specific binding to HeLa cells has been mapped to residues 460–496 in the m1 sequence (37). The m2 23-aa insertion is located within this region, but an engineered chimera in which these 23 m2 amino acids were inserted into an m1 sequence did not result in a loss of HeLa cell toxicity (37). Excluding the m2 insertion, there are only 10 amino acids in the 460–496 m1 sequence that consistently differ between the m1 and m2 forms (SI Fig. 7). The position of these residues in the context of the 23-aa m2 insertion is impossible to predict, and, therefore, an understanding of structural features that give rise to VacA cell-type specificity may need to await analysis of an m2 p55 structure.
Analysis of the VacA sequence alignments in the context of the m1 p55 structure reveals that disruptions of β-sheet contacts within the β-helix correlate with sites of homologous recombination events. For example, VacA residue D455 (strain 60190) represents the beginning of the midregion that defines m1 and m2 strains (Fig. 2b and SI Fig. 7) and follows a kink in the structure that disrupts β-sheet contacts between β8 and β11 (Fig. 2a and SI Fig. 7). This structural disruption and the associated sequence transition are illustrated with a change from red to orange (Fig. 2 and SI Fig. 7). Other examples of a correlation between VacA subdomains and recombination events are evident based on analysis of rare naturally occurring m1/m2 chimeras, which have m1 character in the N-terminal portion of the VacA midregion but m2 character in the C-terminal portion. The transition point from m1 to m2 in the chimeric VacA from strain 249a3 is at D493, whereas the transition point for the ch2 and v225 strains of VacA is at N616 (Fig. 2b and SI Fig. 7). As illustrated in Fig. 2a, these sequence transitions correlate with disruptions in β-sheet contacts between β14 and β17 and between β29 and β32, respectively. The final subdomain begins at residue G648. Although we did not identify a sequence with evidence of recombination at this site, the subdomain (blue) has longer loops than other subdomains and does not make β-sheet contacts with β33 from the previous subdomain (Fig. 2a).
In summary, the polymorphisms that differentiate m1 and m2 forms correspond almost entirely to surface-exposed residues. Diversification in these surface-exposed regions of VacA may have been driven by immune selective pressures. Recombination occurs commonly in H. pylori (32, 33), but m1/m2 chimeras are only rarely identified (35). This suggests that H. pylori strains with intact m1 or m2 VacA sequences have favorable functional properties and, thus, have a selective advantage compared with strains containing chimeric m1/m2 sequences (36–38). The rare naturally occurring m1/m2 chimeras shown in Fig. 2b and SI Fig. 7 probably arose as the result of recombination events that maintained favorable structural and functional properties in the encoded β-helix structures.
Assembly of VacA into Oligomeric Structures.
To understand the p55 structure in the context of VacA oligomers, we have docked p55 into the 19-Å cryo-EM map of the wild-type dodecamer (24) using the program COLORES (Fig. 4a) (48). The elongated shape of the β-helix and the curve of the C-terminal foot allow for unambiguous placement of 12 p55 subunits into the petal-like features of the map. The docking suggests that contacts between the p55 components of the two hexameric layers are mediated by three loops and that the N-terminal end of the p55 β-helix extends into the central density of the map (Fig. 4b). The conserved pocket that we propose as a potential common receptor binding site for m1 and m2 forms of VacA is located on the sides of the petals and would be fully accessible to cell-surface receptors when VacA is assembled in either a single- or double-layered oligomeric structure (Fig. 4b).
We hypothesize that a large portion of p33 will adopt and extend the β-helix fold observed in p55. This idea is supported by two observations. First, an extended β-helix structure is consistent with the shape of a p88 monomer obtained by EM (Fig. 4c Inset) (24). Second, the program BetaWrapPro identifies the stretch of amino acids between residues 120 and 249 as a five-coil β-helix and predicts additional β-helix strands in the region between 252 and 288 (49). BetaWrapPro uses profile wrapping for prediction and comparative modeling of β-helices and has been shown to identify the β-helix motif with high sensitivity and selectivity (49).
Based on our prediction that a large portion of p33 will extend the β-helix fold observed in p55, we suggest a model in which VacA oligomerization is mediated by contacts between p33 and the N terminus of p55 from a neighboring subunit (Fig. 4c). Specifically, one would predict the p55 oligomerization surface to contain the β-strands β3 and β6. The surface-exposed residues of β3 and β6 are accessible within the central density of the EM map and are strictly conserved among the 92 m1 and m2 sequences we surveyed, suggesting selective pressure to preserve this surface (Fig. 4b). An important functional role of this region is also supported by biochemical data indicating that VacA residues 1–422 can induce vacuolation in cultured cells when expressed intracellularly, but p33, p55, and VacA residues 1–394 cannot (25, 50). Finally, this model is supported by a yeast two-hybrid experiment in which p33 and residues 313–478 of p55 were shown to interact (27). The oligomerization surface may also contain parts of the 312–354 sequence that was not visible in this structure. Much of this region is strictly conserved among the m1 and m2 VacA sequences (SI Fig. 7).
p33 (residues 1–311) contains a putative α-helical pore-forming domain at its N terminus (28, 51). Comparison of the VacA wild-type and Δ6–27 dodecamer structures by cryo-EM suggests that the N-terminal domain of p33 is located in the central density of the structure (24). Secondary structure prediction analyses suggest that p33 has α-helical structural elements between residues 1 and 71. Short β-strands are then predicted to begin at residue 87 and extend through a region (residues 120–288) that is predicted to have a β-helical structure, based on BetaWrapPro (49). We therefore anticipate that a single subunit of VacA will adopt a roughly symmetric shape within the context of the oligomer and that residues in the N-terminal region of p33 (likely C-terminal to the pore-forming region) would be positioned to mediate oligomerization with the N-terminal region of a neighboring p55 subunit (Fig. 4c).
A protective H. pylori vaccine would potentially reduce the incidence of H. pylori infection and the serious complications of peptic ulcer disease and gastric adenocarcinoma (52), but the high level of diversity among H. pylori strains and the high frequency of genetic recombination among strains are likely to present significant challenges to vaccine development. VacA is a candidate antigen because immunization with VacA confers protective immunity in a mouse model of H. pylori infection (52, 53). We show here that two VacA surfaces are strictly conserved among all surveyed m1 and m2 sequences and propose that these surfaces are under selective pressure to be preserved to mediate receptor-binding and oligomerization functions. Efforts to selectively target these regions could result in progress toward the development of a vaccine that confers protection against multiple strains of H. pylori.
Methods
Expression and Purification of p55.
The p55 domain of VacA from H. pylori strain 60190 (a type m1 form of VacA; GenBank accession no. Q48245) with an N-terminal hexahistidine tag was expressed as described (26) except that E. coli BL21(DE3) cells were used. Harvested cells were resuspended in lysis buffer (50 mM potassium phosphate, 300 mM NaCl, 20 mM imidazole at pH 8.0) and lysed by three passages through a pressure dispersion homogenizer at 20,000 psi. Cell lysates were centrifuged at 48,000 × g for 20 min at 4°C. Native p55 was purified from the supernatant by Ni-affinity, ion exchange, and gel filtration chromatography. The sizing column retention time suggested a molecular mass of ≈130 kDa. We hypothesize that this result reflects the presence of an elongated p55 dimer. Seleno-methionine p55 was produced from E. coli BL834(DE3) in a minimal medium containing 40 mg/liter l-seleno-methionine and purified by using methods similar to those used for the native protein.
Crystallization and Preparation of Heavy-Atom Derivatives.
The native and seleno-methionine p55 were concentrated in a 100 mM NaCl, 50 mM potassium phosphate buffer and crystallized at 21°C by the hanging-drop vapor-diffusion method in which protein and precipitant were mixed in a 1:1 ratio. Crystals grew from starting protein concentrations of 5.5 mg/ml, and reservoirs containing 18–23% PEG1500. Heavy-atom derivatives were prepared by soaking crystals in 50 mM potassium tetrabromoplatinate (IV) (K2PtBr6) for 3 d, 5 mM 1,4-diacetoxymercuri-2,3-dimethoxybutane for 24 h, or 5 mM gold chloride (HAuCl4) for 24 h. Crystals were mounted on cryo loops, sequentially soaked in cryoprotectant solutions containing 25% PEG1500, and 5%, 10%, and 15% glycerol, and flash-cooled in liquid nitrogen.
Structure Determination and Refinement.
X-ray data were collected from single crystals at 100 K on beamline ID-22 at the Advanced Photon Source (Argonne, IL) by using a MAR300 image plate detector. All of the diffraction data were indexed, integrated, scaled, and merged with HKL2000 (SI Table 1) (54). Native and derivative crystals were space group P3221 and had unit cell dimensions of a = b = 57 Å, and c = 257 Å. Phases were determined by MIRAS using a native protein crystal and seleno-methionine, platinum, mercury, and gold derivatives (SI Table 1). Heavy-atom positions were located and refined with the autoSHARP program suite (55). The structure was subjected to iterative rounds of model building in O (56) and refinement in CNS (57) before applying TLS refinement (58) in REFMAC (59). The refined model (Rcryst = 18.1%, Rfree = 22.2%) consists of amino acids 355–366, 377–811, and 206 water molecules. SDS/PAGE and N-terminal sequence analysis of dissolved crystals indicated that proteolysis had occurred at the N terminus of p55 during crystallization, thus explaining the absence of residues 312–354 and 367–376 from the crystal structure. Proteolysis was prevented by the addition of protease inhibitors, but crystals did not form under these conditions.
Acknowledgments
We thank Catherine El-Bez and Jacques Dubochet for providing the cryo-EM map. X-ray data were collected at the Southeast Regional Collaborative Access Team (SER-CAT) 22-ID beamline at the Advanced Photon Source, Argonne National Laboratory (Argonne, IL). Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under contract no. W-31-109-Eng-38. This work was also supported by National Institutes of Health Grant R01 AI39657 (to T.L.C.), the Medical Research Service of the Department of Veterans Affairs (T.L.C.), the Sartain-Lanier Family Foundation and an American Cancer Society Institutional Research Grant IRG-58-009-48 (to D.B.L.), and the T32 GM 08320 Molecular Biophysics Training Program (K.A.G.).
Footnotes
The authors declare no conflict of interest.
Data deposition: The atomic coordinates and structure factors have been deposited in the Protein Data Bank, www.pdb.org (PDB ID code 2QV3).
This article contains supporting information online at www.pnas.org/cgi/content/full/0707447104/DC1.
References
Articles from Proceedings of the National Academy of Sciences of the United States of America are provided here courtesy of National Academy of Sciences
Full text links
Read article at publisher's site: https://doi.org/10.1073/pnas.0707447104
Read article for free, from open access legal sources, via Unpaywall: https://europepmc.org/articles/pmc2042200?pdf=render
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
Smart citations by scite.ai
Explore citation contexts and check if this article has been
supported or disputed.
https://scite.ai/reports/10.1073/pnas.0707447104
Article citations
Synthesis, DFT, and in silico biological evaluation of chalcone bearing pyrazoline ring against Helicobacter pylori receptors.
Heliyon, 10(14):e34540, 14 Jul 2024
Cited by: 0 articles | PMID: 39130476 | PMCID: PMC11315094
Structural Analysis of Membrane-associated Forms of Helicobacter pylori VacA Toxin.
J Mol Biol, 436(4):168432, 30 Dec 2023
Cited by: 1 article | PMID: 38161000 | PMCID: PMC11090358
There Are No Insurmountable Barriers: Passage of the <i>Helicobacter pylori</i> VacA Toxin from Bacterial Cytoplasm to Eukaryotic Cell Organelle.
Membranes (Basel), 14(1):11, 28 Dec 2023
Cited by: 1 article | PMID: 38248700 | PMCID: PMC10821523
Review Free full text in Europe PMC
In Silico Discovery of Anticancer Peptides from Sanghuang.
Int J Mol Sci, 23(22):13682, 08 Nov 2022
Cited by: 1 article | PMID: 36430160 | PMCID: PMC9693127
Phylogenetic Classification and Functional Review of Autotransporters.
Front Immunol, 13:921272, 01 Jul 2022
Cited by: 11 articles | PMID: 35860281 | PMCID: PMC9289746
Review Free full text in Europe PMC
Go to all (105) article citations
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
Genes & Proteins
- (1 citation) UniProt - Q48245
Protein structures in PDBe
-
(1 citation)
PDBe - 2QV3View structure
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Analysis of a beta-helical region in the p55 domain of Helicobacter pylori vacuolating toxin.
BMC Microbiol, 10:60, 23 Feb 2010
Cited by: 11 articles | PMID: 20178613 | PMCID: PMC2836311
Cryo-EM Analysis Reveals Structural Basis of Helicobacter pylori VacA Toxin Oligomerization.
J Mol Biol, 431(10):1956-1965, 05 Apr 2019
Cited by: 13 articles | PMID: 30954575 | PMCID: PMC6625667
Sticky socks: Helicobacter pylori VacA takes shape.
Trends Microbiol, 16(3):89-92, 14 Feb 2008
Cited by: 19 articles | PMID: 18280164
Review
Functional properties of the p33 and p55 domains of the Helicobacter pylori vacuolating cytotoxin.
J Biol Chem, 280(22):21107-21114, 06 Apr 2005
Cited by: 46 articles | PMID: 15817461
Funding
Funders who supported this work.
NIAID NIH HHS (2)
Grant ID: R01 AI39657
Grant ID: R01 AI039657
NIGMS NIH HHS (2)
Grant ID: T32 GM 08320
Grant ID: T32 GM008320