Bacteriophage BPP-1 infects and kills Bordetella species that cause whooping cough. Its diversity-generating retroelement (DGR) provides a naturally occurring phage-display system, but engineering efforts are hampered without atomic structures. Here, we report a cryo electron microscopy structure of the BPP-1 head at 3.5 Å resolution. Our atomic model shows two of the three protein folds representing major viral lineages: jellyroll for its cement protein (CP) and HK97-like (‘Johnson’) for its major capsid protein (MCP). Strikingly, the fold topology of MCP is permuted non-circularly from the Johnson fold topology previously seen in viral and cellular proteins. We illustrate that the new topology is likely the only feasible alternative of the old topology. β-sheet augmentation and electrostatic interactions contribute to the formation of non-covalent chainmail in BPP-1, unlike covalent inter-protein linkages of the HK97 chainmail. Despite these complex interactions, the termini of both CP and MCP are ideally positioned for DGR-based phage-display engineering.https://doi.org/10.7554/eLife.01299.001
Whooping cough is a respiratory illness caused by bacteria in the Bordetella genus. Among the general public, Bordetella species have become a hot topic in recent years due to the re-emergence of whooping cough in the United States and elsewhere. Scientists, meanwhile have become interested in a virus called BPP-1 that can kill the Bordetella species.
BPP-1 is a double-stranded DNA virus, and such viruses have long been of interest to scientists because they are the most abundant organisms on Earth. These viruses are also noteworthy because their shells (also known as capsids) are capable of withstanding the very high pressures (up to about 40 atmospheres) that are created by packing so much DNA into the very small volume inside the capsid.
BPP-1 is of particular interest because it is capable of making large-scale changes to its own DNA in order to adapt to changes in its hosts and environment. Of all the organism that do not contain nuclei within their cells (collectively known as prokaryotes), BPP-1 is the only one that is capable of making such changes to its DNA. However, efforts to exploit the properties of BPP-1 for bioengineering applications have been hampered because its detailed structure is not known. Now Zhang et al. have used cryo electron microscopy to study the structure of BPP-1 at the atomic level.
Most viruses belong to one of three major lineages, with each lineage having a characteristic fold in its capsid proteins. Zhang et al. found that BPP-1 contains two of these folds, which suggests that it is a hybrid of two of these lineages. This is the first time that such a structure has been observed. Moreover, Zhang et al. found that one of the folds has an unusual topology that has not been seen before. The atomic structure reveals how double-stranded DNA viruses use a variety of non-covalent interactions and a type of protein ‘chainmail’ to form a highly stable capsid that is capable of withstanding very high pressures.
In addition to enabling applications in bioengineering, the new structure might also provide insights into the evolution of prokaryotes.https://doi.org/10.7554/eLife.01299.002
Capsid proteins of non-enveloped viruses fall, so far, into three major structural classes: the β-jellyroll, the HK97 fold, and the fold of dsRNA-virus shell proteins (Bamford et al., 2005; Oksanen et al., 2012). The RNA bacteriophage MS2 subunits have a fourth structure, not yet found in eukaryotic viruses (Valegard et al., 1990). The HK97 fold is present not only in a large number of dsDNA bacteriophage capsids, including the T-phages, lambdoid phages, etc, but also in the major capsid protein of eukaryotic herpesviruses (Bamford et al., 2005). In HK97 itself, the intersubunit contacts are reinforced by a post-assembly covalent linkage—an isopeptide bond between adjacent gp5 subunits, so positioned that the entire capsid is topologically interlinked in a ‘chainmail’ arrangement (Duda, 1998; Wikoff et al., 2000). In other cases, such as bacteriophage λ, the non-covalent interactions between the subunits are reinforced by an additional ‘cement’ protein, which binds on the outer surface of the capsid at positions close to those of the isopeptide bonds in HK97 (Lander et al., 2008). Heads of bacteriophage λ defective in this cement protein break down during DNA packaging (Sternberg and Weisberg, 1977; Fuller et al., 2007; Lander et al., 2008).
Bacteriophage BPP-1 is a short-tailed, dsDNA virus and a member of the Podoviridae family. It infects and kills Bordetella species that cause whooping cough in humans and respiratory diseases in other mammals. It has a T = 7l icosahedral capsid, ∼670 Å in diameter. A 7 Å resolution cryoEM structure of BPP-1 showed that its capsid protein has an HK97-like fold, but the shell has an additional protein component (Dai et al., 2010). Lacking information at the time about their genetic identities, these proteins were named according to their structural roles: a major capsid protein (MCP) and a ‘cement protein’ (CP) that decorates the shell.
From a biotechnology standpoint, BPP-1 has emerged as an attractive phage-display system thanks to its unique and well characterized diversity-generating retro-element (DGR). As the only known source of massive DNA sequence deviation in prokaryotes, DGRs use a unique reverse transcriptase-based mechanism to introduce targeted diversity into protein-coding DNA sequences to accelerate the evolution of adaptive traits (Liu et al., 2002; Guo et al., 2008). As such, BPP-1 is a naturally occurring diversity-generating system and can be engineered to display foreign proteins with adaptive heterologous sequences. An atomic description of the BPP-1 head will enhance bioengineering efforts, reveal non-covalent molecular interactions conducive of stable bacteriophage capsid formation, and clarify evolutionary relationships of BPP-1 MCP and CP with proteins in other viruses.
In this study, we report the three-dimensional (3D) structure of the BPP-1 head at ∼3.5 Å resolution determined by single-particle cryoEM and derived an atomic model. Both our structural and structure-based mutagenesis studies reveal major novelties in the BPP-1 MCP and CP structures and their interactions. We also show that the C-termini of both MCP and CP are ideally positioned to display DGR-diversified peptide libraries for protein engineering applications.
Genomic analysis indicated that the 42.5 kb BPP-1 genome encodes up to 50 viral proteins (Liu et al., 2004). To identify the genes coding for MCP and CP, we first carried out SDS-PAGE analysis and showed that the two most abundant protein bands have molecular masses of 36.3 kD and 15.2 kD (Figure 1A). Mass spectrometry analysis confirmed that the bands correspond to Bbp17 and Bbp16, respectively (Figure 1B). The theoretical molecular masses of Bbp17 and Bbp16 (i.e., 36417.2D and 14458.3D) match their apparent molecular masses by SDS-PAGE, suggesting that the two capsid proteins in the mature bacteriophage particles are not cleaved (Figure 1A). Finally, to verify that genes bbp16 and bbp17 encode components essential for forming infectious phage particles, we constructed in-frame deletions of the bbp16 and bbp17 genes. As expected, no infectious phage particles were produced from BPP-1Δbrt lysogens carrying deletions in bbp16 or bbp17, confirming that both MCP (Bbp17) and CP (Bbp16) are required for phage production. For the sake of clarity and ease of comparison with structural homologs of other viruses, we will continue to use CP and MCP to refer to Bbp16 and Bbp17 of BPP-1, respectively. As shown below, these assignments are directly verified through atomic modeling in which side-chain structures of amino acid residues in sequences match side-chain densities visualized in the cryoEM map (Figure 2).
We reconstructed the 3D structure of the BPP-1 icosahedral head by single-particle cryoEM (Figure 2, Video 1). The quality of MCP and CP densities was further improved by additional averaging of seven CP or six hexameric MCP subunits (except for the pentameric MCP) in the asymmetric unit (Zhang et al., 2010a). Based on the reference-based Fourier shell correlation criterion (FSC = 0.143), the resolution of the capsid is 3.58 Å (Figure 2C) (Rosenthal and Henderson, 2003). Consistently, the R-factors of atomic models are better than 0.5 at the zone of 1/3.5 Å (for capsid) or 1/3.4 Å (for averaged CP and MCP) (Figure 2C; Table 1), which correspond to a Fourier shell correlation (FSC) coefficient greater than 0.143 (Wolf et al., 2010). Therefore, the resolution of the capsid and averaged densities were estimated to be ∼3.5 Å and ∼3.4 Å, respectively. This assessment of resolution is also consistent with clearly visible side chain densities of bulky amino acid residues, such as arginine and phenylalanine (Figures 2D and 4E–F; Videos 2–4).
Arranged on a T = 7 icosahedral shell, MCP and CP each have seven copies or conformers in each asymmetric unit. MCP alone forms a complete icosahedral shell without any noticeable gaps. CP forms dimers bound to the underlying MCP shell at local and icosahedral twofold axes, making the maximum diameter of the BPP-1 (670 Å) bigger than that of the HK97 capsid (659 Å). When CP dimers are computationally removed, the overall size and architecture of the BPP-1 capsid is almost identical to that of the HK97 capsid, and the backbone model of the HK97 capsid is nearly super-imposable with the BPP-1 density map with minor mismatches. Based on the averaged density maps of CP and MCP, we built initial atomic models of CP and MCP with Coot (Emsley and Cowtan, 2004), and refined them with Phenix (Adams et al., 2010). Models of the seven MCP and CP monomers within the asymmetric unit were subsequently obtained by adjusting the initial MCP and CP models with Coot and refined with Phenix, with overall R/R-free factors of 0.26/0.27 at 3.5 Å resolution (Table 1) (‘Materials and methods’).
We traced 139 (Met1 to Ile139) of the 140 amino acid residues of the cement protein Bbp16 in the cryoEM density map. CP contains two β-sheets each consisting of four anti-parallel strands, and a ∼36-Å long extension (Thr125 to Ala140) with the C-terminus exposed on the external surface (Figure 3; Video 2). The backbones of the seven CP subunits in the asymmetric unit, including their C-terminal extensions, are nearly identical, with an RMSD of 0.4–0.66 Å when superimposed (Figure 3D).
The topological organization of strands in the two β sheets (i.e., BIDG and CHEF; Figure 3B) of CP is identical to that of the jellyroll motif. This structural motif was first discovered in spherical RNA viruses (Harrison et al., 1978; Abad-Zapatero et al., 1980; Hogle et al., 1985; Rossmann et al., 1985; Chelvanayagam et al., 1992), and subsequently found widely in other DNA viruses, such as ϕX174, bacteriophage PRD1 and human adenovirus (McKenna et al., 1992; Abrescia et al., 2004; Zubieta et al., 2005; Liu et al., 2010; Krupovic and Bamford, 2011). Unlike its role as the major capsid protein in the above-mentioned viruses, the jellyroll motif in BPP-1 forms an auxiliary protein to stabilize the viral capsid made by proteins of another fold. Notably, in some dsRNA viruses, the jellyroll motif has been adopted as a domain in a stabilizer/adaptor protein that forms trimers located at an intermediate layer of viral capsids (Grimes et al., 1998; Mathieu et al., 2001; Liemann et al., 2002; Zhang et al., 2010a). In BPP-1, CP subunits form a dimer through β-sheet augmentation between the F-strands of two CHEF sheets, forming an 8 stranded β-sheet visible on the external surface (Figure 3C). Although the BIDG sheet of the jellyroll motif faces but do not interact with the underneath MCP shell, instead, as described in detail below, CP interacts with the MCP shell mainly through its N- and C-termini as well as a linking loop between its β-strands B and C.
We traced 327 amino acid residues (Ser5 to Val331) of the total 331 residues of the major capsid protein Bbp17 in the cryoEM density map (Figure 4; Video 3). Although little sequence identity (CLUSTALW score = 5) was detected between BPP-1 MCP and HK97 gp5 proteins with ClustalW (Thompson et al., 1994), the overall architecture of the MCP resembles that of HK97 gp5 (Wikoff et al., 2000) (Figure 4—figure supplement 1A), with an axial (A) domain, a peripheral (P) domain and an extended loop (E-loop) (Figure 4C). Both the subunit organization and domain orientations of the seven MCPs in the BPP-1 capsid are also identical to those of the corresponding gp5 subunits in the HK97 capsid. The A-domain contains a 5-stranded β sheet flanked by two helices on one side and a C-terminal loop on the other. The P-domain contains a characteristic long, 7-turn helix sandwiched by a three-stranded β sheet and the N-terminal loop. The E-loop contains a single amino acid sequence segment folded into an extended hairpin loop projecting from a 2-stranded β sheet. Notably, both N- and C-termini of MCP are exposed on the external surface of the capsid, thus are accessible for tethering peptides in phage-display applications.
The structures of the seven MCP conformers in the asymmetric unit are nearly identical with some minor differences (Figure 4—figure supplement 1B), and RMSDs of their backbones (residues from 31–60 to 81–331) are 0.55–0.98 Å among six hexameric MCPs and 1.2–1.8 Å between the pentameric and the hexameric MCPs. The minor differences among the different MCP conformers include: (1) Each N-terminal portion (Ser5–Glu30) of the seven MCPs adopts a slightly different conformation, and backbone RMSDs of this segment are 1.4–5.0 Å among seven MCP copies; (2) The E-loops of the seven MCPs have backbone RMSDs of 0.8–8.3 Å, larger than other portions of MCPs, possibly resulting from different local shell curvatures at the regions of the E-loops (Figure 4E); (3) Although the C-terminal proximal loop (Ala295–Val309) of every hexameric MCP is well resolved, the corresponding segment in the pentameric MCP is not resolved, suggesting disordered conformation of this segment in pentons possibly due to steric hindrance at the central channel of the MCP pentamer (Figure 4—figure supplement 2).
Despite the similarities in the architectural appearances of BPP-1 MCP and HK97 gp5, their atomic structures differ in four significant ways. Firstly, the conformations of the N-terminal loops are different (Figure 4—figure supplements 1, 2C). In BPP-1 MCP, the N-terminal peptide extends radially to interact with a CP by augmenting its CHEF β-sheet. In HK97 gp5, the N-terminus extends circumferentially and interacts with three adjacent gp5 subunits, augmenting a β-sheet in one of them. Secondly, a C-terminal loop (Lys293–Val309) of the pentameric MCP in BPP-1 is disordered while the corresponding loop (Arg294–Thr304) in the HK97 pentameric gp5 is ordered (Figure 4—figure supplement 2B,C) (Wikoff et al., 2000). Thirdly, the electrostatic potential properties of the two proteins are strikingly different. Finally, and perhaps most interestingly, the folding topologies of BPP-1 MCP and HK97 gp5 are different, which, as described below, provides the first opportunity to explore possible alternative ways to build the highly conserved HK97-like folds (Figure 5, Figure 5—figure supplements 1 and 2; Video 4).
As remarked above, despite of the similar architectural appearance (or fold) of BPP-1 Bbp17 and HK97 gp5, the peptide traces of these two proteins follow different folding topologies (Figure 5A–B, Figure 5—figure supplements 1 and 2). A careful comparison of the two structures has led us to identify three structural elements of the canonical HK97-like, or ‘Johnson’ fold: N-, β- and α-elements (marked cyan, purple and green in Figures 4A–D and 5A–C). These three structural elements join together through a central 5-stranded (F, E, G, K, L in Figure 5A–C, in an up-up-down-up-up topology) β-sheet, flanked by two short helices (Figure 5A–C, Figure 5—figure supplement 2). The β-element (purple) consists of two anti-parallel (i.e., G, K; down-up) β-strands forming a hairpin located at the middle of the Johnson fold. The α-element (green) contains two parallel (i.e., F, E; up-up) β-strands and the two short helices. The N-element (cyan) contributes the last (i.e., L; up) β-strand to the Johnson fold and also contains the two other characteristic secondary structures of the Johnson fold: a long, kinked ‘spine’ α-helix and the extended loop (E-loop).
In BPP-1 MCP, the N-, β-, and α-elements consist of residues 1–168, 169–241, and 242–331, respectively (Figure 4A). The orders of these three structural elements in BPP-1 and HK97 are N-β-α (‘BPP topology’, Figures 4A–C and 5A), N-α-β (‘HK97 topology’, Figure 5B), respectively. To verify this, we swapped the α and β structural elements in our de novo BPP-1 MCP model to create an ‘interchanged model’ that matches the HK97 topology. This ‘interchanged model’ was then refined with Phenix for five cycles. Despite this refinement effort, the final ‘interchanged model’ does not agree with our experimental cryoEM density (Figure 4—figure supplement 3). In particular, many side chains in the ‘interchanged model’ do not match those resolved in the cryoEM density map (Figure 4—figure supplement 3B,C). This verification provides additional support to the BPP topology established by our de novo modeling approach.
Interchange of two of the structural elements can lead to proteins with different topologies of the Johnson fold (Figure 5C). From a pure mathematical point of view, permuting three structural elements can give rise to a total of six different topologies (i.e., 3! = 6), considering that sequence polarity (N- to C-termini) of the structural elements must be preserved. However, the N-element (cyan) cannot participate in this permutation due to the remote disposition (∼42 Å) of its N-terminus with respect to the C-termini of other structural elements, and thus the total possible topologies is reduced to only 2 (i.e., 2! = 2) and the permutation is non-circular (Figure 5—figure supplements 1A–B, 2G–H). Remarkably, the newly discovered BPP topology presented here is therefore the second topology predicted from the above rule (Figure 5A), the first being the topology discovered in HK97 gp5 (Wikoff et al., 2000) (Figure 5B) and subsequently seen in many other viruses, as well as in bacteria and archaea cells (Akita et al., 2007; Sutter et al., 2008). Insertions into any of the three structural elements can give rise to more complex architectures without increasing the number of folding topologies of the Johnson fold, as would be expected for larger viral particles such as the herpesvirus capsid.
To test whether a functional protein can be produced from BPP-1 MCP by permuting it to the HK97 topology without introducing other changes, we genetically interchanged the primary order of the β- and α-elements in the BPP-1 MCP gene bbp17 (Figure 5—figure supplement 1C), and made two slightly different constructs (Figure 5—figure supplement 1C) (i.e., PM1 and PM2, ‘Materials and methods’). These two engineered proteins were expressed successfully, as indicated by Western blot analysis (Figure 5D). However, no plaque formation was detected for the PM1 and PM2 lysates on RB50 cells transformed with either the wt bbp17, the PM1 or the PM2 construct (Figure 5E), suggesting that the PM1 and PM2 gene products are not functional.
The BPP-1 capsid has a chainmail structure similar to that of HK97 capsid (Figure 6A–B; Video 5). Instead of covalent bonding as in HK97, the BPP-1 chainmail is stabilized by non-covalent interactions, such as: (1) strong electrostatic interactions between adjacent MCPs (Figure 6), and (2) additional CP-MCP interactions contributed by CPs (Figure 7).
Electrostatic interactions between MCPs of BPP-1 are stronger than those between gp5 subunits in the HK97 capsid (Figure 6D–H, Figure 6—figure supplements 1) for two reasons. Firstly, in the MCP hexamer, the electrostatic interactions between adjacent MCPs that appears to be stronger than those between HK97 hexameric gp5 proteins (Figure 6—figure supplement 1A–B). Secondly, the interface between the MCP E-loop and an adjacent MCP protein at each local threefold region contains complementary electrostatic interactions (Figure 6D–E). Specifically, the E-loop is mainly positive charged and its MCP interface is mainly negative charged, and the interaction at this interface appears to be stronger than that of HK97 gp5 in the same interface. Interestingly, the electrostatic properties of BPP-1 MCP and Hk97 gp5 are opposite in this interface (Figure 6—figure supplements 1A–B). The stronger electrostatic interaction between the E-loops and the local threefold region of the BPP-1 MCP may be the reason for the absence (Figure 8) of the salt bridge found between HK97 gp5 proteins at the local threefold interface that is critical for the assembly, stability and maturation of the HK97 capsid (Gertsman et al., 2010).
In addition, the BPP-1 head is further stabilized by interactions between CP and MCP. Complementary electrostatic potential surfaces were also identified at the interface between CP and MCP (Figure 7A–C). Each CP interacts extensively with five MCP subunits underneath, mainly involving three loops: (a) N-terminal loop, (2) C-terminal loop, and (3) the linking loop between strands C and D of the jellyroll motif (Figure 7D–F). Firstly, the N-terminal loop of CP interacts with two MCPs: the N-terminal loop of first MCP (green) and the E-loop of the second MCP (yellow) (Figure 7D–E), Secondly, the C-terminal loop of CP interacts with two MCPs: a β-strand of the β-element of the first MCP (orange) and both N- and C-terminal loops of a second MCP (green) (Figure 7D–E). Thirdly, the linking loop of CP interacts with three MCPs: two β-strands of the β-element of the first MCP (blue), the N-terminal loop of a second MCP (green) and a short loop (Asn48–Glu52) of the third MCP (orange) (Figure 7D–E). Fourthly, β-strand C of CP jellyroll interacts with the N-terminal loop of an MCP, augmenting the CHEF β-sheet of CP (Figure 9C). Finally, two CP monomers form a dimer through anti-parallel interaction between their F-strands of jellyrolls, forming an augmented, 10-stranded β sheet and further stabilizing the capsid (Figures 3C, 7D and 9C).
We found that BPP-1 CP and MCP adopt two of the three characteristic folds of major virus lineages, jellyroll, and HK-97/Johnson. Moreover, the topology of the Johnson fold in BPP-1 MCP is non-circularly permuted as compared to the previously known topology of the Johnson fold. Until now, the Johnson fold has been observed at atomic detail in the Siphoviridae and Myoviridae families of the three tailed-bacteriophage families, exemplified by gp5 of HK97 (Siphoviridae) (Wikoff et al., 2000) and gp24 of T4 (Myoviridae) (Fokine et al., 2005). Our atomic model derived from the 3.5 Å cryoEM map of BPP-1 represents the first atomic model for a virus in the third tailed-bacteriophage family, the Podoviridae. Indeed, BPP-1 represents the first among all viruses shown to contain both the jellyroll and HK97-like folds. In the virosphere, the conventional wisdom is that viruses have evolved from a few distinctive lineages separately. Three major lineages have been proposed based on three highly distinctive structures: jellyroll-like (i.e., jellyroll lineage, include single- and double-jellyroll), HK97-like (i.e., HK97 lineage) and BTV-like (i.e., BTV lineage), and members of each lineage contain MCPs with conserved folds originating from a common ancestor (Rossmann and Johnson, 1989; Bamford et al., 2005; Abrescia et al., 2012). Our result suggests that BPP-1 is likely a hybrid virus that has adopted the characteristic structures from both the jellyroll and HK97 lineages.
The head of bacteriophage Epsilon15 also contains two capsid proteins (gp7 and gp10) located at similar positions as CP (Bbp16) and MCP (Bbp17) in BPP-1. In particular, BPP-1 Bbp17 (MCP) shares ∼46% sequence identity with gp7 of Epsilon15 phage, indicating that the two proteins are likely to have nearly identical structures at least at the level of protein backbones. Indeed, although initially misinterpreted to have a different topology (Figure 5—figure supplement 2H–I) due to limited resolution (Jiang et al., 2008), the model of Epsilon15 gp7 was recently updated on the basis of an improved 4.5 Å resolution map (Baker et al., 2013) to bear the same BPP topology as we initially reported (Zhang et al., 2012) and detailed further in this study.
As described in the results above, the BPP topology and HK97 topology represent the only two feasible topologies of the Johnson fold (Figure 5). To date, all other known atomic structures of proteins with the Johnson fold adopt the HK97 topology, including gp5 of HK97 in Siphoviridae (Wikoff et al., 2000), gp24 of T4 in Myoviridae (Fokine et al., 2005), a 39kD spherical particle-forming protein in archaea (Akita et al., 2007) and an encapsulin protein in bacteria (Sutter et al., 2008) (Figure 5—figure supplement 2A–E). In these two topologies, the primary positions of their β- and α-elements, are swapped in a non-circularly permutated fashion (Figure 5A–C, Figure 5—figure supplement 1), which rarely occurs in nature (Vogel and Morea, 2006). This, combined with the fact that no recognizable sequence similarity can be identified between the BPP-1 MCP and other HK97-like MCPs (CLUSTALW scores: ∼5), obscures the origin of the BPP topology.
From an engineering stand point, it is fortunate that both termini of MCP and the C-terminus of CP are exposed on the external surface of the capsid, ideally positioned for peptide display. Indeed, we observed infectious phage particle formation when either CP or MCP was fused to a peptide of 22 amino acids at their C-terminus (HG and JFM, unpublished observation). As BPP-1 DGRs diversify DNA sequences in vivo, this suggests they may be engineered to display heterologous sequences on the phage capsid by fusion to the termini of CP or MCP. For DGR-based phage-display, the diversity of the DNA library is not limited by the efficiency of DNA transformation, and both the library construction and optimization can occur entirely inside bacterial cells without the need for in vitro manipulation. Moreover, the selected phage particles can be easily re-amplified for iterative selection since the tail, not the engineered head, is involved in phage-receptor recognition. The atomic model of the BPP-1 head presented in this study provides a road map forward to harness the great potentials afforded by the unique properties of BPP-1 for phage-display engineering.
500 ml of LB medium was inoculated with a single colony of BPP-1 lysogen. The cultures were incubated at 37°C on a rotary shaker until log phase when phage production was induced by adding mitomycin C to a final concentration of 2 mg/l. After 3 hr of induction, CHCl3 was added with shaking to facilitate cell lysis and phage release. Cellular debris was removed by centrifugation at 5000×g. Phage particles were then precipitated using 10% PEG8000/500 mM NaCl, pelleted by centrifugation, and resuspended in 50 mM Tris-HCl, 250 mM NaCl, 1 mM MgCl2, pH7.5. The resuspended phage particles were further purified by 15–45% (buffered with 50 mM Tris-HCl, 250 mM NaCl, 1 mM MgCl2, pH7.5) sucrose gradient ultracentrifugation at 35,000 rpm for 90 min using an SW41 rotor in a Beckman L8-80M ultracentrifuge. The phage band was visualized by illuminating the gradient tube and carefully collected by side puncture. Purified phage was buffer-exchanged and concentrated in 50 mM Tris-HCl, 250 mM NaCl, 1 mM MgCl2, pH7.5 by ultra-filtration with a MW cutoff of 100 kD (Millipore, Billerica, MA). The final concentration of BPP-1 was ∼1014 pfu/ml.
To identify capsid proteins using SDS−PAGE and mass spectrometry, purified BPP-1 particles were incubated in 1× SDS sample buffer (0.0625 M Tris-HCl, 1.25% SDS, 5% glycerol and 0.02% bromophenol blue, pH 6.8) at 95°C for 5 min before loading onto the gel. SDS−PAGE was performed using a discontinuous gel (4% polyacrylamide stacking region and 15% polyacrylamide resolving gel) prepared with a discontinuous gel system (Bio-Rad Laboratories, Hercules, CA). Protein bands were visualized by Coomassie Blue staining (Bio-Rad Laboratories). The two major bands (Figure 1A) were manually excised, digested in-gel with sequencing-grade modified trypsin (Promega, Madison, WI), and analyzed using Matrix-assisted laser desorption/ionisation-time of flight mass spectrometry (MALDI-TOF MS). Scaffold (version.3.00.04, Proteome Software Inc., Portland, OR) was used to validate MS/MS-based peptide and protein identifications.
To determine the atomic structure, low dose images of liquid-nitrogen-cooled, frozen hydrated BPP-1 were recorded on Kodak SO-163 film on an FEI Titan Krios cryo electron microscope operated at 300 kV with dose of ∼25e−1/Å2 on specimen. The nominal magnification of images is 59,000×, which was previously calibrated to be 57,660× using tobacco mosaic virus as a standard. Imaging condition was optimized by using parallel illumination and by minimizing beam tilt with a Coma-free alignment procedure.
895 films were recorded and scanned using Nikon Coolscan9000 scanners under step-size of 6.35 μm, corresponding to 1.1 Å/pixel on the specimen scale. Under-defocus value of these films were determined using CTFFIND (Mindell and Grigorieff, 2003), and based on the cross correlation coefficients from CTFFIND3, 340 films with defocus values of −0.8∼−2.27 μm were selected for further processing. Particles were selected by an in-house automatic procedure and followed by visual screening to remove particles near edge of films, and total 43,156 particles with an image size of 880 × 880 were boxed. Data processing and 3D reconstruction were accomplished with an integrative approach using Frealign (Grigorieff, 2007) and eLite3D (Zhang et al., 2010b) on high performance computers including graphics processing units (GPU) as previously described (Zhang et al., 2010a). For global search, 2× binned images were used to save time. At the end of each cycle of image processing, the effective resolution was determined and refinement for the next cycle would be carried out by including data up to this resolution. In the last refinement cycle, data up to a spatial frequency of 1/3.7 Å−1 was included for the refinement and 39,549 particles were used for the final reconstruction. To further improve signal/noise ratio, local averaging of seven CP or six hexameric MCP subunits (the pentameric MCP is slightly different from the hexameric MCP subunits and was not included in the averaging) in the asymmetric unit was performed as previously described (Zhang et al., 2010a).
The effective resolutions were estimated based on the 0.5 R-factors between the cryoEM density map [equivalent to FSC≥0.143 (Rosenthal and Henderson, 2003; Wolf et al., 2010) and the final atomic model calculated by Phenix (Adams et al., 2010). The calculated R-factors reached ∼0.5 at 3.5 Å, 3.4 Å for the capsid density and the averaged CP and MCP densities, respectively (Table 1). These estimations are consistent with the structural features present in the maps (Figure 2, Videos 2–4). The capsid and averaged maps were filtered to 1/(3.4 Å) and 1/(3.5 Å) spatial frequency, respectively; and sharpened using a reverse B-factor of −200 Å2 (for capsid) or −250 Å2 (for averaged densities), which were estimated through a trial-and-error procedure by optimizing side-chain densities and noise level simultaneously. Visualization and segmentation of density maps were done with UCSF Chimera (Pettersen et al., 2004).
Based on the averaged density maps of CP and MCP, we first built initial Cα and full atom models for CP and MCP with Coot (Emsley and Cowtan, 2004) without referring to any existing models of other proteins. For MCP, the N- and C-termini were distinguished based on the ‘Christmas tree’ polarity of α helices and confirmed by landmark, bulky side chain densities. For CP, because there is no helix to be used to reveal the N- to C-terminal polarity, we determined the N-terminus using the side chain densities of some landmark amino acids, such as Phe27, Tyr 102, Tyr106, Tyr107 and Tyr133. The initial full atom models were regularized by constraining both Ramachandran geometry and secondary structures in Coot (Emsley and Cowtan, 2004) but without including hydrogen atoms.
These initial full atom models were iteratively refined using structural information of both amplitude and phase (from Fourier transform of the cryoEM maps) in the following three distinct (two automatic and one manual) steps. The first step (automatic) was performed in Phenix (Adams et al., 2010) using Ramachandran restrain. The second step (automatic) is to regularize the new model also using Phenix. To do so, hydrogens were added to all atoms of the model from the last refinement, and followed by regularization and removal of hydrogens. The latest models were refined iteratively until no further improvement was apparent based on both Ramachandran geometry and R-factors. Then in the third, manual step, amino acid residues with invalid Ramachandran backbone geometries were identified and their backbone Psi–Phi angels were manually corrected in Coot. This process of automatic- and manual model refinement steps was iterated until no further improvement on both Ramachandran geometry and R-factors was evident.
Atomic model of an asymmetric unit including seven MCPs and seven CPs was subsequently obtained by adjusting the refined MCP and CP models according to the density maps of different conformers. This model of the asymmetric unit was refined with Phenix under the constraints of Ramachandran geometry, secondary structures, and icosahedral symmetry (Adams et al., 2010). Clashes at the molecular boundaries across different asymmetric units in the entire capsid were minimized by including icosahedral symmetry constraints in this refinement. The final atomic model of the full capsid was obtained after 14 cycles of refinement.
Because our BPP-1 MCP chain trace differs from that of the HK97 gp5 in the order of the α and β structural elements (‘Results’), we made an extra cautious effort to verify our trace by swapping the α and β structural elements in our model to create an interchanged model that matches the HK97 topology. This interchanged model was then refined with Phenix (Adams et al., 2010) for five cycles. Most of the side chains of the refined interchanged model do not match those in the cryoEM density, further confirming our de novo model (Figure 4—figure supplement 3).
To engineer BPP-1 MCP to match the topology of HK97 gp5, we swapped the order of the β- and α-elements in BPP-1 MCP (Figure 5—figure supplement 1C). In construct PM1, the β-element (peptide 169–241) was cut from wt bbp17 gene and then pasted to the C-terminal end of α-element. Construct PM2 was obtained the same way except that the cut sites were shifted by three amino acid residues on both sides, resulting in a six residue-longer β-element (peptide 166–244). In both the constructs, the N-terminal end of the α-element was pasted to the C-terminal end of the N-element.
Plasmids expressing either wt bbp17, PM1 or PM2 genes B. bronchiseptica were transformed into RB50 cells transformed by electroporation. These transformed RB50 cells were grown on Bordet–Gengou agar containing 15% sheep blood, 25 μg/ml streptomycin and 25 μg/ml chloramphenicol. The protein expressing levels of the wt bbp17, PM1 and PM2 genes (each tagged with 6xhistidine at the C-terminus) in these RB50 cells were determined by Western blot. Single colonies were inoculated into Luria–Bertani media containing 25 μg/ml streptomycin, 25 μg/ml chloramphenicol and 10 mM nicotinic acid and grown at 37°C overnight. Cells in the amount of 1 ml at OD600 = 1.0 were pelleted and then grown in 2.5 ml of Stainer Scholte media with 25 μg/ml streptomycin and 25 μg/ml chloramphenicol to induce expression of the bbp17, PM1 or PM2 gene. Cells expressing each construct were grown for 3 and 6 hr. Equal amounts of the cells (by OD600) were harvested and lysed by boiling in 1 × SDS-PAGE loading buffer, and subsequently analyzed with SDS-PAGE. Western blot was done with a mouse monoclonal antibody against 6xhistidine as the primary antibody and a horse reddish peroxidase-conjugated goat anti-mouse antibody as the secondary antibody with an Amersham kit.
To obtain phage lysates for plaque assays, we generated a lysogen (BPP-1ΔbrtΔbbp17) that has the bbp17 gene deleted. Then, we transformed plasmids expressing either the wt bbp17 gene, the PM1 or PM2 genes into the BPP-1ΔbrtΔbbp17 lysogens. These lysogen cells were first grown at 37°C for 3 hr in Stainer Scholte medium containing 25 μg/ml streptomycin and 25 μg/ml chloramphenicol to induce to the Bvg+ phase, leading to the expression of the wt bbp17 and the two mutants PM1 and PM2. Mitomycin C (0.2 μg/ml) was then added to the cell cultures to induce phage production. After 3 hr, chloroform was added to the cultures, followed by vortexing and centrifugation to remove cellular debris. The resultant supernatants were collected for plaque assays on B. bronchiseptica RB50 cells transformed with plasmids expressing either wt bbp17, PM1 or PM2 gene. As above, these transformed RB50 cells were grown on Bordet–Gengou agar containing 15% sheep blood, 25 μg/ml streptomycin and 25 μg/ml chloramphenicol.
CryoEM density maps, atomic models and surface charge properties were visualized with Chimera (Pettersen et al., 2004).
The cryoEM density map and the atomic model of BPP-1 have been deposited to databanks with accession numbers EMD (5764, 5765, 5766) and PDB (3J4U), respectively.
Structure unifies the viral universeAnnu Rev Biochem 81:795–822.https://doi.org/10.1146/annurev-biochem-060910-095130
Validated near-atomic resolution structure of bacteriophage epsilon15 derived from cryo-EM and modelingProc Natl Acad Sci USA 110:12301–12306.https://doi.org/10.1073/pnas.1309947110
What does structure tell us about virus evolution?Curr Opin Struct Biol 15:655–663.https://doi.org/10.1016/j.sbi.2005.10.012
Three-dimensional structure of tropism-switching Bordetella bacteriophageProc Natl Acad Sci U S A 107:4347–4352.https://doi.org/10.1073/pnas.0915008107
Virus universe: can it be constructed from a limited number of viral architecturesIn: G Witzany, editors. Viruses: essential agents of life. Netherlands: Springer. pp. 83–105.https://doi.org/10.1007/978-94-007-4899-6_5
UCSF Chimera–a visualization system for exploratory research and analysisJ Comput Chem 25:1605–1612.https://doi.org/10.1002/jcc.20084
Icosahedral RNA virus structureAnnu Rev Biochem 58:533–573.https://doi.org/10.1146/annurev.bi.58.070189.002533
Structural basis of enzyme encapsulation into a bacterial nanocompartmentNat Struct Mol Biol 15:939–947.https://doi.org/10.1038/nsmb.1473
Stephen C HarrisonReviewing Editor; Harvard Medical School, United States
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for sending your work entitled “A new topology of the HK97-like fold revealed in Bordetellabacteriophage: non-covalent chainmail secured by jellyrolls” for consideration at eLife. Your article has been favorably evaluated by a Senior editor and 3 reviewers, two of whom, Stephen Harrison and NikoGrigorieff, have agreed to reveal their identity.
Here is a combined summary of the reviews, with the five key areas that need attention before the paper can be accepted:
The structure of the Bordetella phage BPP-1, determined by cryoEM at a resolution of about 3.5 Å, contains three broadly interesting findings. The first of these is a description of the major capsid protein chain trace, showing an interchange (in the sequence) of structural elements, with respect to their order HK97. The second is the structure of the “cement” protein, which is present in a number of other bacteriophages (e.g., lambda) for which chain-trace resolution structures have not been available. The third is the overall structure of the capsid, and information about the positions of the N- and C-termini of the two proteins, because of the potential application to a phage display platform. A revision that convinces us even more completely of the chain trace and that pays attention to several additional issues should be appropriate for publication in eLife.
Chain trace: the difference between the chain trace for this virus (which is the same as the “revised” trace from Houston for phage epsilon15) and HK97 is curious, although the authors (whose work is, as usual for this group, technically superb) present good evidence for it. Some additional tests and illustrations are requested in time (4) below. Because the Houston group’s trace of epsilon15 was initially wrong (they've now revised it to a version that agrees with the BPP-1 trace), the present authors need to be completely convincing. (The authors have to convince us completely that they now have it right, even if HK97 and T4 and other viruses have it the other way.)
The following revisions should make the paper suitable for publication. Overall, it is too long and redundant, and the illustrations are often unclear. In addition to avoiding repetition (they can easily shorten the paper to no more than 2/3 its current length, even with the additions suggested below), the authors need to pay attention to the following five points.
1) The authors misunderstand and misstate the history of discovery and understanding of viral capsid subunit structure.
2) They need to provide more details about data processing.
3) They confuse structural observations with energetic ones, so their speculations about “withstanding high pressures” (moreover, is 40 atm really that high, on a molecular scale?) are misplaced.
4) They need to present some further validation of the difference between BPP-1 and HK97 and to describe the chain trace and the overall structure more clearly. The current description makes it impossible for anyone except the initiated to understand the interesting relationship between BPP-1 MCP and the HK97 coat protein, unless the reader realizes that all they need do is look at Figure 6. Moreover, the description of structural “blocks” is misleading, as the two strands in the central sheet and the long beta hairpin are parts of distinct structural modules.
5) They should tone down evolutionary speculations. (Francis Crick once said something to the effect that we should always imagine evolutionary scenarios and talk about them, to sharpen our thinking, but we shouldn't publish them, unless there is good evidence.)
1) History. It was never assumed that “all viruses would use the jellyroll fold”. There were, of course, lumpers and splitters – at least one reviewer was basically agnostic (see Crick quote). But the RNA phage structure (1990) showed right away that you could make a perfectly nice, T=3 capsid with a different fold. Amusingly enough, that fold has only been found in RNA phage. Anyhow, the dsRNA viruses and HK97 show two more designs, although most of the dsRNA viruses have another protein (e.g., VP6 of rotavirus) that does have a jellyroll. So while it remains puzzling why jellyrolls show up in so many non-enveloped viruses, playing the “my fold is more curious/novel/interesting than your fold” game is silly.
To rectify the misunderstanding, and to avoid unnecessary “fold-ology” (folds are hardly the most interesting thing about viruses), the following two paragraphs are a proposed revision and condensation of the relevant sections. The revision takes into account two points from item (3) – the misunderstanding about pressure and the proper use of the “chainmail” analogy.
“Capsid proteins of non-enveloped viruses fall, so far, into three major structural classes: the β-jellyroll, the HK97 fold, and the fold of dsRNA-virus shell proteins. The RNA phage subunits have a fourth structure, not yet found in eukaryotic viruses. The HK97 fold is present not only in a large number of dsDNA bacteriophage capsids, including the T-phages, lambdoid phages, etc., but also in the major capsid protein of eukaryotic herpesviruses. In HK97 itself, the intersubunit contacts are reinforced by a post-assembly covalent linkage - an isopeptide bond between adjacent gp5 subunits, so positioned that the entire capsid is topologically interlinked in a “chainmail” arrangement. In other cases, such as lambda, the non-covalent interactions between the subunits are reinforced by an additional “cement” protein, which binds on the outer surface of the capsid at positions close to those of the isopeptide bonds in HK97. Heads of lambda phage defective in this cement protein break down during DNA packaging.”
“BPP-2 is a short-tailed, dsDNA virus, a member of the Podoviridae family. It infects and kills Bordetella species that cause whooping cough in humans and respiratory diseases in other mammals. It has a T=7 icsoaehedral capsid, ∼670 Å in diameter. A 7 Å resolution cryoEM structure of BPP-1 showed that its capsid protein has an HK97-like fold, but further suggested that the shell has an additional protein component. Lacking information at the time about their genetic identities, these proteins were named according to their structural roles: a major capsid protein (MCP) and a “cement protein” (CP) that decorates the shell.”
2) Details about data processing. The authors should include, in the section of “CryoEM image, data processing and resolution assessment”, answers to the following questions: i) How was the magnification calibrated? ii) At what resolution was the refinement done? iii) Was pixel binning used? iv) What filtering of the final map was done (low-pass, B-factor)? v) How many images were collected? vi) How many particles were selected? vii) How many particles were used in the final reconstruction? viii) How were “good” particles selected? ix) How was the CTF determined? x) What defocus range was used? xi) How big a box was used for particle windowing?
One reviewer was not particularly impressed by an R-factor comparison between calculated and observed maps, especially if 0.5 is the cut-off. Refinement and free R do not mean quite the same thing here as in x-ray refinements. The structure factors against which one refines are already correlated, so saving out a test set doesn't mean much, and the non-capsid symmetry (ncs) averaging correlates them further. Visual comparison of the maps is obviously the most useful criterion, and that is indeed the criterion on which the authors have largely relied. Use of the Scheres“gold standard” might be helpful to assess the resolution more correctly. Even if the authors choose not to do that calculation, full FSC and R-factor curves should be shown because their shapes sometimes indicate problems with the refinement, e.g., over-fitting.
3) Structure and energetics. It turns out that 40atm is not so difficult a pressure to withstand. 40 atm is 4 MPa. To expand the shell by ΔR = 1 Å, probably a reasonably tolerable range for the intermolecular forces holding the shell together (the circumference would change by 6 Å, distributed over more than six intersubunit contacts as you go around, so any intersubunit interface would stretch by less than 1Å, even assuming that the subunits themselves were utterly inextensible), the DNA would do work PΔV ∼ 4 x 10^-18 Joules or 2400 kJ/mole (taking 300 Å as the radius). That seems a lot, but it is only about 5.5 kJ/mole of subunit (there are 420), or roughly the strength of an extra H-bond (i.e., 1.5 kcal/mole). Not particularly impressive. So the chainmail of HK97 is probably not about withstanding “enormous forces”, but rather about withstanding the slings and arrows of the slime through which the poor phage particle must pass to get from one bacterial host to another. That doesn't make the alternative strategies of cement protein and isopeptide bond less interesting – but it is important to avoid hyperbole based on intuitive macroscopic analogies.
“Chainmail” requires topological interlinks. You can't have “noncovalent chainmail”, without stretching the armor analogy farther than a reasonable metaphor should go.
4) Chain trace comparisons.
First, it is essential to be 100% sure that the chain trace is correct. The authors should show clear examples of how an alternative trace violates the density. In particular, as the altered connections are at the surface of the particle, where local disorder, etc., could obscure correct links, they should show that any plausible sequence register for the HK97-like trace is not possible. One simple way to do this would be to examine a series of possible alternative sequence registers for residues following the first of the “non-HK97-like” connections and show that the pattern of hydrophilic and hydrophobic residues is completely different. In other words, that if you look at the possible sequence registers for the long beta-hairpin that runs parallel to the long helix (probably just 4-5 alternative starting and ending points), do any of them look plausible enough to test by the refinement criteria that Zhang et al. use for their trace? A suitable figure showing various sequence registers might be all that is needed to make this point.
The authors should also show density, as in Figure 2D, for two segments that are affected by the interchange, to show that there is no ambiguity in how the interchanged sequences fit density. Density extending from the new junctions and including side chains consistent with the sequence in the reordered form would be the best way to do this.
Second, assuming the trace to be correct, the description is confusing and doesn't follow the way we generally describe domains or modules in proteins. First, there are only two modules here that can reasonably be called “domains” (although this is a funny protein, so I'm not sure that “domain” is the right word). The “blocks” that the authors name here are not proper folding units, as the long β-hairpin in their β-block is really part of the P-domain (to use the terminology in Figure 4) and the β-strands are part of the A-domain. Thus, Figure 6B and D are misleading. They are formally correct, but they do no correspond to how we imagine a protein might fold. The authors should stick to FigureA and C and not try to force the question with three “building blocks” as in Figure 6E, the unnecessary calculation of possible “permutations”, etc.
Third, the right word for what has happened here is an “interchange” of the order in the polypeptide chain of a long, two part β-element (effectively a β-hairpin with a right-angle bend, which allows it to participate in two distinct sheets). This element is at the C-terminus in HK97, and it has been pasted in front of strand E (L is barely a strand in HK97) in BPP-1. “Non-circular permutation”, while perhaps a mathematically correct use of the word “permutation”, will confuse phage biologists, for whom “circular permutation” has an important reference in the way T-even phage genomes relate from particle to particle.
5) Evolution and structure comparisons. We have no idea how viral evolution relates to host organism, species jumping, etc. It is not a great idea to speculate about “popular views” of how viruses “originated”. While anyone who contributes solid data to the published literature buys a license for a few paragraphs of speculation, and the authors are entitled, should they wish, in the paragraph beginning “One popular view argues that viruses originated from ancient cells…”, it might be wise to reduce or eliminate it.
Three other points in connection with structure comparisons.i) Is the hand laevo or dextro? All T=7 capsids examined so far are laevo; if BPP-1 turned out to be dextro, this would be another distinctive feature (what is epsilon 15?). ii) When considering subunit folds of different lineages, many authors differentiate between the single jellyroll fold observed in the first ssRNA virus structures determined and the double jellyroll found in dsDNA viruses infecting all three domains of life. When multiple jellyrolls from RNA viruses are in a single polypeptide chain like the large subunit of cowpea mosaic virus or all 3 subunits of tobacco ring spot virus they clearly appear as “beads on string”, while the double jelly roll structure integrates the two domains into a contiguous protein. iii) The authors state that BPP1 probably has a different assembly and maturation pathway compared to HK97. In order to make this statement meaningful they need to discuss the scaffolding protein of BPP1 and any intermediates in maturation that may have been characterized. HK97 has a covalently associated scaffold that is proteolytically removed during maturation. Is there any sequence of events known for BPP1? If not, that statement should be softened.https://doi.org/10.7554/eLife.01299.024
- Z Hong Zhou
- Jeff F Miller
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Pavel Afonine and Paul Adams for advice on using Phenix to refine the atomic models. We acknowledge the use of the cryoEM facility in the Electron Imaging Center for Nanomachines by NIH (1S10RR23057 to ZHZ) and CNSI at UCLA.
- Stephen C Harrison, Reviewing Editor, Harvard Medical School, United States
- Received: July 27, 2013
- Accepted: October 27, 2013
- Version of Record published: December 17, 2013 (version 1)
© 2013, Zhang 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.