|Home | About | Journals | Submit | Contact Us | Français|
Human bocavirus (HBoV) was recently discovered and classified in the Bocavirus genus (family Parvoviridae, subfamily Parvovirinae) on the basis of genomic similarity to bovine parvovirus and canine minute virus. HBoV has been implicated in respiratory tract infections and gastroenteric disease in children worldwide, yet despite numerous epidemiological reports, there has been limited biochemical and molecular characterization of the virus. Reported here is the three-dimensional structure of recombinant HBoV capsids, assembled from viral protein 2 (VP2), at 7.9-Å resolution as determined by cryo-electron microscopy and image reconstruction. A pseudo-atomic model of HBoV VP2 was derived from sequence alignment analysis and knowledge of the crystal structure of human parvovirus B19 (genus Erythrovirus). Comparison of the HBoV capsid structure to that of parvoviruses from five separate genera demonstrates strong conservation of a β-barrel core domain and an α-helix, from which emanate several loops of various lengths and conformations, yielding a unique surface topology that differs from the three already described for this family. The highly conserved core is consistent with observations for other single-stranded DNA viruses, and variable surface loops have been shown to confer the host-specific tropism and the diverse antigenic properties of this family.
Human bocavirus (HBoV), a newly discovered member of the family Parvoviridae, was originally isolated in randomly selected nasopharyngeal aspirates (5). Since this initial discovery, HBoV has also been detected worldwide, predominantly in children under the age of 2 years with respiratory infections, in serum, urine, and fecal samples (40). Symptomatic children commonly exhibit acute diseases of the upper and lower respiratory tracts (7, 36, 44, 56) and, possibly, gastroenteritis (31, 56) though the link to gastroenteritis outbreaks has been questioned (12). It is still unclear if HBoV is the sole etiologic agent of respiratory disease as higher rates of coinfections with other respiratory pathogens such as human rhinovirus and Streptococcus spp. are often observed (4). However, Allander et al. recently reported (4) that HBoV was found in 19% of children with acute wheezing, thereby making it the fourth most common virus, after rhinoviruses, enteroviruses, and respiratory syncytial virus, detected in children exhibiting this symptom. These findings suggest that, at high viral load, HBoV could be an etiologic agent of respiratory tract disease (4). HBoV infection is common in the first few years of life, and clinical research suggests it may follow the primary period for acquisition of human parvovirus B19 (B19) though there is no antigenic cross-reactivity between B19 and HBoV (28, 30). By age 5, most people have circulating antibodies against HBoV, as is also true for other respiratory viruses such as respiratory syncytial virus, rhinoviruses, and human metapneumovirus (17). HBoV has also been identified in adults, with ~63% of samples tested being seropositive, showing a positive correlation with age and a slight positive bias toward women (14).
The Parvoviridae is a family of small, nonenveloped viruses that package a single-stranded DNA (ssDNA) genome of ~5,000 bases. These viruses are subdivided into two subfamilies: Parvovirinae and Densovirinae (Table (Table1).1). The Parvovirinae are further subdivided into five genera, all of whose members infect vertebrates. The Densovirinae (four genera) infect only invertebrates. Phylogenetic analysis places HBoV in the recently classified Bocavirus genus (Table (Table1).1). In addition to HBoV, numerous parvoviruses circulate among the human population. Among these are the following: several dependoviruses; adeno-associated virus (AAV) serotypes AAV1 to AAV3, AAV5, and AAV9; the Erythrovirus B19; and the newly discovered human parvovirus genotypes 4 (Parv4) and 5 (Parv5) (23, 27, 50). Of these, only B19 had been implicated in disease until the discovery of HBoV and Parv4, which has been isolated from patients who present symptoms of acute HIV infection (50).
The HBoV genome, like that of all members of the Bocavirus genus, contains three open reading frames (ORFs). The first ORF, at the 5′ end, encodes NS1, a nonstructural protein. The next ORF, unique to the bocaviruses, encodes NP1, a second nonstructural protein. The third ORF, at the 3′ end, encodes the two structural capsid viral proteins (VPs), VP1 and VP2. The HBoV VPs share 42% and 43% amino acid sequence identity with the corresponding VPs of bovine parvovirus and canine minute virus, respectively (5). More recently, two additional HBoV-like viruses, HBoV-2 and HBoV-3, were identified in stool samples from children (8, 31). The genome organization of these viruses is identical to that of HBoV, with the NS1, NP1, and VP proteins of HBoV-2 and HBoV-3 being, respectively, ~80 and 90%, ~70 and 80%, and ~80 and 80% identical to the respective proteins in HBoV (8, 31).
Parvovirus genomes are packaged into a T=1 icosahedral capsid that is assembled from 60 copies of a combination of up to six types of capsid VPs (VP1 to VP6), all of which share a C terminus. VP1 is always a minor component, typically comprising about five copies per capsid, whereas the smallest VP is always the major component. The unique N-terminal region of VP1 (VP1u) contains a conserved phospholipase A2 (PLA2) motif within the first 131 amino acids that is essential for infection (49, 61). Interestingly, Aleutian mink disease virus (AMDV), the only member of the Amdovirus genus, is the only exception in that this motif is absent, which suggests that this virus employs a different mechanism to escape the endosome during infection (54).
The X-ray crystal structures of several parvoviruses show that all VPs contain a conserved, eight-stranded β-barrel motif (βB to βI) that forms the core of the capsid (15). There is also a conserved α-helix (αA) observed in all parvovirus structures determined to date. The bulk of the VP consists of elaborate loops between the strands that form the surface of the capsid. For example, the GH loop between the βG and βH strands is ~230 residues. The composition and topology of these loops encode several important functions, including tissue tropism, pathogenicity, and the antigenic response directed against each parvovirus during infection (2).
A number of parvoviruses have been studied by cryo-electron microscopy (cryo-EM) and three-dimensional (3D) image reconstruction in concert with and complementary to X-ray crystallographic studies (reviewed in reference 15). Reported here is the 3D structure of a recombinant HBoV capsid solved to 7.9-Å resolution using cryo-EM. The capsid of HBoV was compared to that of representative members of the Parvoviridae (Table (Table1)1) with known atomic structures (AAV2, minute virus of mice [MVM], B19, and Galleria mellonella densovirus [GmDNV]) or pseudo-atomic models built into cryo-EM reconstructed density (AMDV) to identify similarities and differences. The capsid topology of the newly emerging HBoV incorporates a combination of surface structural features seen in other members of the Parvovirinae and is closest to that of B19, the only other structurally characterized parvovirus that is pathogenic to humans. A pseudo-atomic model of the HBoV VP2, built into the reconstructed density, identified conserved core secondary structure elements, which are known to be important for parvovirus capsid assembly, and variable surface loops, which likely govern host specific interactions.
The putative HBoV VP2 gene was constructed by de novo synthesis using the published HBoV ST1 sequence (NCBI accession number DQ000495; start codon, nucleotide [nt] 3373), and cloned into a pUC vector (Epoch Biolabs, Inc., Missouri City, TX). The vector was cut with EcoRV and HindIII to release the insert and recloned into a Fastbac1 plasmid (Invitrogen). The sequence and orientation of the resulting plasmid were confirmed, and the plasmid was used to produce a recombinant baculovirus using the standard Bac-to-Bac technology (Invitrogen).
Sf9 insect cells (ATCC) were maintained in Grace's medium (Invitrogen) with 10% fetal calf serum (FCS) and antibiotics; cells were infected with recombinant baculovirus, harvested 4 or 7 days postinfection into the medium, and collected by centrifugation. Production of baculovirus was confirmed by a combination of immunofluorescence of infected cells with mouse antibody against baculovirus (eBioscience) and the typical cytopathic effect of baculovirus-infected cells. For the immunofluorescence testing, slides were prepared by cytocentrifugation (200 × g for 8 min in a Shandon Cytospin 4) or by spotting a concentrated cell suspension onto each of eight spots of a fluorescent antibody slide (Bellco Glass, Inc.). The slides were allowed to air dry and were fixed in 1:1 acetone-methanol at −20°C for 10 min. The monoclonal antibody was diluted 1:50 in phosphate-buffered saline (PBS) plus 10% FCS, and 25 μl of the solution was incubated with each cell spot (~200 cells) for 1 h at 37°C. After cells were washed with PBS, they were incubated with 1:100 anti-mouse fluorescein isothiocyanate (FITC)-conjugated antibody in PBS plus 10% FCS and 1:200 Evans Blue (Sigma) for 1 h at 37°C. After another washing step, the cells were examined by UV microscopy.
Sf9 cells were grown in suspension at 27°C in Sf-900 II SFM medium (Gibco/Invitrogen Corporation) and infected at a multiplicity of infection of 5 PFU per cell. Resultant HBoV virus-like particles (VLPs) were released from infected cells by three freeze-thaw cycles in lysis buffer (50 mM Tris-HCl, pH 8.0, 100 mM NaCl, 1 mM EDTA, 0.2% Triton X-100), with the addition of benzonase (Merck KGaA) after the second cycle. The sample was clarified by centrifugation at 12,100 × g for 15 min at 4°C. The cell lysate was pelleted through a 20% (wt/vol) sucrose cushion (in 25 mM Tris-HCl, pH 8.0, 100 mM NaCl, 0.2% Triton X-100, 1 mM EDTA) by ultracentrifugation at 149,000 × g for 3 h at 4°C. The pellet from the cushion was resuspended in the same buffer overnight at 4°C. The sample was subjected to multiple low-speed spins at 10,000 × g in order to remove insoluble material. The clarified sample was purified by separation within a 10 to 40% (wt/vol) sucrose step gradient spun at 151,000 × g for 3 h at 4°C. Equilibrium dialysis was performed against 20 mM Tris-HCl, pH 7.5, 150 mM NaCl, and 2 mM MgCl2 at 4°C. The final concentration of the sample was 10 mg ml−1 calculated from optical density measurements at 280 nm, with an extinction coefficient of 1.7 M−1 cm−1. The purity and integrity of the viral capsids were monitored by SDS-PAGE and negative-stain electron microscopy, respectively.
Small (3.5 μl) aliquots of purified VLPs (~5 mg ml−1) were applied to a continuous carbon support film that had been glow discharged for 15 s in an Emitech K350 glow-discharge unit and subsequently stained with 1% aqueous uranyl acetate. Micrographs were recorded on a 2,048- by 2,048-pixel Gatan charge-coupled device (CCD) camera in an FEI Sphera microscope at a nominal magnification of ×40,000 at an accelerating voltage of 200 keV.
Small (3.5 μl) aliquots of purified VLPs (~10 mg ml−1) were vitrified via standard rapid freeze-plunging procedures (1, 20). Samples were applied to Quantifoil holey grids that had been glow discharged for ~15 s in an Emitech K350 glow-discharge unit. Grids were then blotted for ~5 s, plunged into liquid ethane, and transferred into a precooled FEI Polara multispecimen holder, which maintained the specimen at liquid nitrogen temperature. Micrographs were recorded on Kodak SO-163 electron image film at 200 keV in an FEI Polara microscope under minimum-dose conditions (~9 e/Å2) at a nominal magnification of ×39,000.
Sixteen micrographs exhibiting minimal astigmatism and specimen drift, recorded at underfocus settings ranging between 0.95 and 3.25 μm, were digitized at 7-μm intervals (representing 1.795-Å pixels) on a Zeiss SCAI scanner. The software package RobEM (http://cryoEM.ucsd.edu/programs.shtm) was used to extract 3,754 individual particle images (each 1872 pixels), preprocess the images, and estimate the defocus level of each micrograph (9). A random-model procedure (59) was used to generate an initial reconstruction at ~25-Å resolution from 150 particle images. This was then used as a starting model to initiate full orientation and origin determinations and refinement of the entire set of images using AUTO3DEM (60). Corrections to compensate for the effects of phase reversals in the contrast-transfer functions of the images were performed as previously described (10, 62), but amplitude corrections were not applied. A final 3D map at an estimated resolution limit of 7.9 Å, based on the 0.5 threshold of the Fourier shell correlation (FSC0.5) criterion (55) (data not shown), was reconstructed from all particle images. To verify the new structure, we computed two, completely independent, ab initio reconstructions from separate sets of micrographs consisting of 1,770 and 1,839 particle images. The resultant reconstructions showed excellent correspondence in all coarse and fine features, as evidenced by good agreement (FSC0.5) at all spatial frequencies out to 1/9.8 Å−1. An inverse temperature factor of 1/100 Å2 was applied to the final 3D reconstruction (26) to enhance fine details and to aid in the analysis and interpretation of the capsid structure. A Gaussian function was applied to attenuate the Fourier data smoothly to zero between 7.24 and 6.75 Å. The absolute handedness of the reconstruction was set to be consistent with surface features seen in parvovirus crystal structures (see below).
Density maps were generated from the crystal structure Cα coordinates of AAV2 (Protein Data Bank [PDB] code 1LP3) (58), B19 (PDB code 1S58) (33), GmDNV (PDB code 1DNV) (51), and MVM (PDB code 1Z1C) (38) and from the Cα pseudo-atomic model of AMDV (42). The coordinates for each VP2 subunit were used to generate a complete icosahedral capsid using the oligomer generator in VIPERdb (13). The CCP4 software suite was used to calculate structure factors (7.9-Å resolution cutoff; B-factor of 100 Å2) and electron density maps, using the SFALL and FFT (fast Fourier transform) routines, respectively (18). The calculated maps were rendered with RobEM (http://cryoEM.ucsd.edu/programs.shtm) and Chimera (48).
VP2 sequences of representative members of the six Parvoviridae genera (taken from the PDB accession files listed above, AMDV [NCBI accession number ABG20974], and HBoV [NCBI accession number ABN46897]) were uploaded to the Biology Workbench (52). A multiple sequence alignment was conducted using ClustalW (53) with default values and with B19 as the reference sequence. An unrooted phylogenetic tree was then created using the PHYLIP algorithm (22).
A homology model was built for the HBoV VP2 (amino acids 37 to 542) using the Swiss-PdbViewer Deep View, version 3.7, software program (25) with the VP2 sequence (NCBI accession number ABN46897) and the B19 VP2 coordinates (PDB code 1S58) supplied as a template. B19 VP2 was chosen as a template for building the HBoV VP2 pseudo-atomic model because the comparison of density maps generated for the members of the different Parvoviridae genera (see above) identified B19 as being the most structurally similar to HBoV even though the AAV2 VP2 sequence had a slightly higher sequence identity.
The HBoV VP2 homology model was used to generate the coordinates for a complete, 60-subunit capsid through icosahedral matrix multiplication in VIPERdb (13). The 60-mer was then visualized within the HBoV cryo-EM density map in COOT (21). The core β-barrel (β-strands B to I) and helix αA elements fitted into density contoured at a 3.3-σ threshold without further adjustment. The loops between β-strands D and E (DE loop) and between β-strands H and I (HI loop) of a reference VP2 monomer were manually adjusted to fit into density contoured at a threshold of 1.0 σ guided by their distinct structural features. The remaining surface loops were adjusted to be within the reconstructed envelope at a density threshold of 0.5 σ. Further model building included addition of amino acids at the N terminus of the VP2 monomer (residues 28 to 36) to satisfy density (0.5-σ threshold) inside the capsid under the icosahedral 5-fold axes. Following these adjustments, the geometry of the model was regularized under idealized geometric constraints in COOT (21). The noncrystallographic symmetry (NCS) function in COOT was used to generate the remaining 59 icosahedral symmetry-related VP2 monomers to verify their fit within the cryo-EM density map. The pseudo-atomic capsid model coordinates were used to generate a density map for comparison with the HBoV cryo-EM density in MAPMAN (37) using the program's similarity function (http://xray.bmc.uu.se/usf/mapman_man.html).
The structures of VP2 from representative members of six Parvoviridae genera (Table (Table1)1) were aligned by superposition of the Cα positions of their atomic coordinates (AAV2, B19, GmDNV, and MVM) or pseudo-atomic coordinates built into cryo-reconstructed densities (AMDV and HBoV) using the MatcherMaker tool (43) in Chimera (48).
Self-assembled VP2 VLPs expressed and purified from insect cells were checked for composition and purity by SDS-PAGE (Fig. (Fig.1A)1A) and for integrity by negative-stain electron microscopy (Fig. (Fig.1B)1B) prior to sample vitrification and cryo-EM (Fig. (Fig.1C).1C). A 3D reconstruction of the HBoV capsid, at an estimated resolution limit of 7.9 Å, was computed from 3,754 individual particle images (Fig. (Fig.2)2) and compared to the capsid structures of representative members of the Parvovirinae (B19, AAV2, AMDV, and MVM) and Densovirinae (GmDNV) subfamilies (Fig. (Fig.33).
The HBoV capsid exhibits three characteristic features common to other vertebrate parvoviruses. The first of these is a dimple-like depression at each icosahedral 2-fold axis (Fig. (Fig.2A2A and and3).3). The second is a large, trimeric protrusion that surrounds each 3-fold axis or is located at the 3-fold axis (Fig. (Fig.2A2A and and3).3). The third is a channel at each 5-fold axis, whose outermost opening is formed by a small, pentameric structure encircled by a wide, canyon-like region (Fig. (Fig.2A2A and and3).3). While the dimple is also observed among the invertebrate parvoviruses, these members lack the 3-fold protrusions and canyon around the 5-fold channel (e.g., GmDNV in Fig. Fig.33).
The external diameter of the HBoV capsid ranges from ~215 Å at the lowest points of the dimple and canyon to ~280 Å at the top of the protrusion. For comparison, the maximum diameters of AMDV and GmDNV (the largest and smallest capsids of the representative parvoviruses being compared) are 310 and 270 Å, respectively, and the minimum diameter calculated from the capsid density maps generated from the atomic or pseudo-atomic coordinates of all the viruses being compared is ~220 Å (Fig. (Fig.3).3). The spherical, relatively smooth topology of the invertebrate GmDNV clearly distinguishes it from the others, especially those with more prominent surface features such as AMDV, MVM, and AAV2. HBoV appears to be morphologically most similar to B19 since both viruses display features intermediate to the above extremes.
Radial density projections of the six representative Parvoviridae capsids reveal that structural similarity among them is greatest at low radii and most divergent at the capsid surface (Fig. (Fig.4).4). Of note, all the vertebrate parvoviruses have strikingly similar structures at radii of <98 Å, suggesting that they share a core. Again, GmDNV stands out among the parvoviruses as having a structure that is most distinct from the others at all radii due to minor positional differences in conserved core secondary structural elements and significant differences in the surface topology of its major capsid VP (see below).
The HBoV capsid surface topology was observed to be most similar to that of B19 (Fig. (Fig.3)3) even though the AAV2 VP2 protein sequence is closer to that of HBoV in percentage identity (Table (Table2).2). Preliminary docking of the crystal structure of the VP2 subunit for both AAV2 and B19 into the HBoV cryo-reconstruction showed that the β-barrel and αA helix in both viruses fitted into the HBoV density (with a MAPMAN  correlation coefficient of 0.6), but for AAV2, the loop regions that form each protrusion extended outside the HBoV envelope (data not shown), consistent with HBoV's smoother surface (Fig. (Fig.3).3). Thus, the B19 crystal structure coordinates (VP2 residues 19 to 554) (33) were used to construct the initial homology model of HBoV VP2 (residues 37 to 542) for docking into the cryo-EM density map (Fig. 5A and B).
A map of HBoV at subnanometer resolution enabled precise, rigid-body fitting of the conserved secondary structural elements (βA, β-barrel core, and αA helix) as well as the DE and HI loops, which required only small manual adjustments, followed by bond geometry refinement (Fig. 5B to D). The remaining interstrand loops and the C terminus of the VP2, including variable regions (VRs) (as defined in reference 24), were adjusted to conform to the HBoV density envelope (Fig. 5A and B). The initial homology model of HBoV started at residue 37, which corresponds to B19 residue 19, the first N-terminal residue reported in the crystal structure (33). However, a clearly defined region of density at the base of the channel adjacent to residue 37 of the initial HBoV model remained uninterpreted following the fitting procedure. This density (Fig. 2C and D, asterisk) was modeled as residues 28 to 36 contributed by five VP2 monomers (Fig. (Fig.5E).5E). These nine residues were modeled as a small loop that projects into the additional density and extends toward the channel, as seen in B19 and MVM (3, 32) (Fig. 5A and E, dark blue loop). This interpretation positions residue 28 of the HBoV model toward the remaining uninterpreted density that lies close to the entrance of the channel and appears to form a “plug-like” structure inside the capsid (Fig. (Fig.2C2C and and5E,5E, black arrow). The correlation coefficient between the map of the 60-subunit capsid (generated from the final HBoV VP2 pseudo-atomic model; residues 28 to 542) and the cryo-EM density map is 0.71. This represents an improvement over 0.6 for the initial B19 model (see above) as a result of model adjustment guided by the reconstructed map. A correlation coefficient of <1.0 is consistent with the inability to precisely place the variable surface loops at the current resolution.
A similar density plug at the base of the 5-fold channel was also observed in the cryo-reconstruction of B19 VLPs comprised solely of VP2 (32). This plug appeared as negative difference density when the B19 VLP cryo-reconstruction was subtracted from reconstructions of wild-type full (VP1, VP2, and DNA) and empty (VP1 and VP2) B19 capsid structures (32). The location of the first ordered N-terminal residue in the B19 VP2 VLP crystal structure (starting with residue 18) overlapped this negative density, suggesting that the remaining, unordered residues lie inside the VLP capsid (32). In the B19 report, no effort was made to model the missing residues into the negative difference density. However, positive difference density was observed in the B19 full and empty capsid minus VLP density map subtractions, which projected into the channel toward the capsid surface, emerging between the DE loop β-ribbons onto the canyon at the capsid surface. This density was modeled as B19 VP2 residues 1 to 17 (32). The differences in locations of the N-terminal residues in B19 VLPs versus capsids with or without DNA and VP1 are likely dictated by the presence of the DNA and/or VP1. The VP1u motif was not observed in B19 capsids (32), most likely because only a few copies are present per particle, and their contributions to the density map would be averaged out by the icosahedral symmetry imposed during structure determination. The B19 channel appears to be closed or tightly constricted at the capsid surface (Fig. (Fig.3)3) (32, 33). Thus, disposition of the N-terminal residues of the B19 capsid proteins on the outer surface is consistent with previous reports that placed VP1u on the B19 surface (6, 29, 34), and this differs from other parvoviruses in which externalization is proposed to occur via the channel. In HBoV the channel is constricted at its base but open at the top (Fig. 2C and D, and and3).3). Thus, it is difficult to assess whether the HBoV VP2 and/or VP1 N termini are also already located on the capsid surface. Our current model of HBoV suggests that the N-terminal residues (Fig. (Fig.5E)5E) may protrude through the channel in virions. Structural studies of capsids that contain VP1 or both VP1 and DNA are needed to verify this prediction.
Sequence alignments and an unrooted phylogenetic tree were generated to quantitatively compare the HBoV VP2 sequence with the corresponding VP sequences in five other parvovirus genera (Table (Table1)1) and to help guide homology modeling based on the HBoV cryo-EM density map. The pairwise sequence identities ranged between 12 and 23%, with GmDNV, the sole invertebrate member in the comparison, being the most divergent from HBoV (Table (Table2).2). Phylogenetic analysis based on a multiple sequence alignment of the six viruses being compared placed them almost equidistant from each other in an unrooted tree (data not shown), suggesting high divergence.
Structural superimposition of the Cα positions of the atomic (B19, AAV2, and MVM) and pseudo-atomic (AMDV) coordinates of the Parvovirinae VP2s with the pseudo-atomic model generated for HBoV (residues 28 to 542) showed that, despite the high (~50%) sequence variability, the eight-stranded β-barrel (BIDG-CHEF) core, βA, and αA helix are conserved in all the viruses (Fig. (Fig.6A).6A). Most VRs in the VP2 structures of the Parvovirinae being compared to HBoV are located within the loops between the β-strands and occur at and around the 3-fold axis (15). In the HBoV pseudo-atomic model, the VRs are also predicted to occur around this axis (Fig. (Fig.5A).5A). However, at the current resolution of the HBoV density map, the exact positions of the loop structures are unknown and thus preclude a comparison with those of the other viruses. As with the Parvovirinae members, the core β-strands are structurally very similar between HBoV and GmDNV though the exact positions of αA differ slightly (Fig. (Fig.6A6A).
Based on the overall surface morphologies of their capsids, parvoviruses have been structurally assigned to three groups (47). Groups I and III comprise members of the Parvovirinae subfamily, all of which share three features as described earlier: dimples at 2-fold axes, protrusions at or surrounding 3-fold axes, and channels surrounded by canyons at the 5-fold axes (Fig. (Fig.3).3). The primary difference between these two groups occurs in the vicinity of each 3-fold axis. Group I capsids have a single, relatively flat, pinwheel-shaped protrusion, whereas group III capsids have three distinct protrusions. Current structural data place members of the Parvovirus genus in group I and members of the Erythrovirus, Dependovirus, and Amdovirus genera in group III (Table (Table1)1) (47). A third topology (group II), seen in GmDNV and other members of the Densovirinae subfamily (11, 51), is characterized by capsids which are relatively spherical and featureless compared to the vertebrate parvoviruses (Table (Table1)1) (47). This morphology, exemplified by GmDNV, results from the reduced length of the surface loops compared to those of members of the Parvovirinae (51).
HBoV possesses features common to groups I and III and also others that are unique. One distinguishing feature in HBoV arises from a difference in surface topology at the 2-fold axis, where the HBoV dimple is narrower and shallower than that in the other parvoviruses (Fig. (Fig.3).3). However, like B19, the HBoV dimple is flanked by wider “walls” between the 2- and 5-fold depressions than are present in the other viruses, and this may reflect a conformational loop difference in the modeled VR-IX located close to the 2-fold axis (Fig. (Fig.5A).5A). Residues that form the walls of the parvovirus dimple have been implicated in receptor attachment for MVM (41) and as determinants of tissue tropism and pathogenicity for other members of the Parvovirus genus (reviewed in reference 38). There is no evidence, however, that this region plays an analogous role in B19 or HBoV.
Other characteristic parvovirus features at the 2-fold region are conserved in HBoV. Helix αA, observed in all parvovirus structures determined to date (reviewed in reference 15), is also present in HBoV (Fig. (Fig.5C5C and and6A).6A). This helix forms stabilizing interactions at the icosahedral 2-fold interface (24). Even though there are fewer interactions at the 2-fold interface than at the 3-fold interface, they are proposed to play a crucial role in parvovirus capsid assembly (24, 57). Conservation of this structural feature in the HBoV VP2 thus suggests that it plays a similar, central role in capsid assembly.
HBoV topology at the 3-fold axis is intermediate between that of groups I and III. The group III viruses, AAV2 and AMDV, have large protrusions that arise from the EF and GH loops. The protrusions in HBoV are less pronounced than in these two parvoviruses, which may result from shorter sequences and/or conformational differences in the HBoV EF and GH loops. This results in a less pronounced depression in this region (Fig. (Fig.3).3). The group I viruses, e.g., MVM, also lack depressions at their 3-fold axes (Fig. (Fig.3).3). Residues near the 3-fold axes in AAV2 and B19 have been implicated as receptor attachment sites (16, 35, 39, 45, 46). The protrusions surrounding these axes are responsible for the antigenic reactivity of members of the Parvovirinae subfamily (reviewed in reference 2). Given this antigenic role of the group I and III surface protrusions, the lack of similar structures in group II viruses is consistent with the fact that invertebrate viruses are not subject to adaptive pressures of the host immune response (51).
As mentioned above, the capsid structures of members of the Parvovirinae have a conserved channel at the icosahedral 5-fold axis, surrounded by a canyon (Fig. (Fig.3).3). The channel is a β-cylinder formed by five symmetry-related β-ribbons inserted between β-strands D and E with a loop at the top (referred to as the DE loop). In HBoV the DE loops form finger-like protrusions that surround the channel (Fig. 5A and D and and6B)6B) and project away from the 5-fold axis in a manner similar to the conformation observed in AAV2, resulting in an open channel (Fig. (Fig.3).3). In contrast, for the other Parvovirinae structures, the DE loops cluster closer to the 5-fold axis and constrict the channel (Fig. (Fig.33 and and6B).6B). This channel in parvoviruses is postulated to serve several functions, for example, as the site of (i) VP2 externalization in members of the Parvovirus genus (e.g., MVM, in which VP3 is generated by a maturation cleavage event), (ii) VP1u externalization (except in B19 and GmDNV) for its PLA2 function, and (iii) genomic DNA packaging. However, because the diameter of the channel varies between ~5 and 25 Å at its narrowest and widest points, respectively, it is likely too small to accommodate extrusion of the ~130 amino acids in the PLA2 domain. It has thus been proposed that flexibility and structural rearrangement of the DE loop (and possibly the HI loop; see below) are required to facilitate these functions. Conservation of this channel in HBoV suggests that it performs analogous functions. Notably, HBoV appears to have the largest surface opening of the parvovirus structures determined to date (Fig. (Fig.33).
The canyon depression in HBoV is less pronounced than the canyons seen in all other parvovirus structures except for GmDNV, which lacks that feature (Fig. (Fig.3).3). In HBoV the HI loop lies in the canyon and is structurally most similar to the HI loop of AAV2 and, although shorter in HBoV, generates a similar surface topology for the two viruses (Fig. (Fig.33 [red loop features on the canyon floor], 5D, and 6B). The HI loop conformation differs significantly in B19 because it projects radially from the capsid surface and lies very close to the DE loop to form a continuous, raised surface in the canyon (Fig. (Fig.3).3). This B19 HI loop conformation is similar to that proposed to occur in AAV2 following receptor attachment (39). The HI loops in MVM and AMDV, which also lie on the canyon floor, are structurally similar to each other but differ from the HI loop of HBoV (Fig. (Fig.6B).6B). In GmDNV the HI loop lies very close to the DE loop, and together they “fill in” the canyon floor (Fig. (Fig.33 and and6B).6B). The HI loop, which has been proposed to undergo structural rearrangements in AAV2 to facilitate VP1 externalization upon receptor binding, has also been reported to play a role in capsid assembly (19, 39).
The structure of the HBoV capsid presented here represents the first for a member of the newly established Bocavirus genus of the Parvoviridae. The subnanometer resolution obtained in the cryo-reconstruction clearly reveals that the newly identified HBoV assembles a capsid that retains the highly conserved structural core adopted by all other parvoviruses and does so despite low sequence similarity between its major capsid protein and the capsid proteins of the other viruses that have existed for many years. Conservation of core structural features in the major coat protein, including the eight-stranded β-barrel, the αA helix, the DE loop, and the HI loop, is consistent with their contributing essential functions in capsid assembly and stability. Conservation of the β-cylinder channel formed by the DE loop strongly suggests that this channel is essential for productive HBoV infection as it likely serves as a portal for VP1u externalization and DNA packaging, as reported for other parvoviruses.
Overall, HBoV shares most of its capsid surface features with the two other human parvoviruses (B19 and AAV2) in this comparative study, and these features result in a unique capsid morphology whose functional domains are yet to be verified. The resemblance of B19 and AAV2 capsids at the 3-fold axes suggests that this region confers similar, receptor attachment capability. On the other hand, the structural difference of the HBoV 5-fold axis from that of B19 (reported to have its VP1/VP2 localized on the capsid surface) and its similarity to that of AAV2 (proposed to externalize its VP1u for its PLA2 function via the channel) suggest that the HBoV capsid is required to undergo conformational transitions similar to the AAV2 capsid during cellular trafficking. Variable HBoV VP2 surface loops, as modeled based on known parvovirus structures (15) and constrained by the reconstructed density envelope, likely reflect adaptations related to host cell recognition and evasion of the adaptive immune response. The structure of HBoV provides an additional framework for the molecular annotation of parvovirus VP structures, which have evolved to facilitate infectivity through the utilization of structural elements and capsid surface features that fulfill myriad functions.
We thank Jinghua Tang for helpful discussions and for verifying the absolute scale of the HBoV density map by quantitative comparison to the B19 crystal structure.
This work was supported in part by NSF grant MCB-0718948 (M.A.-M. and R.M.), by NIH grants R21 AI81072341 (M.A.-M. and R.M.), R37 GM-33050, and 1S10 RR020016 (T.S.B.), by NIH NRSA fellowship F32A1078624 (K.N.P.), and by support from the University of California, San Diego (UCSD), and the Agouron Foundation (T.S.B.) to establish cryo-EM facilities at UCSD.
Published ahead of print on 7 April 2010.