#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Structural and Functional Analysis of Murine Polyomavirus Capsid Proteins Establish the Determinants of Ligand Recognition and Pathogenicity


Viruses are obligate intracellular pathogens, and all of them share one crucial step in their life cycle—the attachment to their host cell via cellular receptors, which are usually proteins or carbohydrates. This step is decisive for the selection of target cells and virus entry. In this study, we investigated murine polyomavirus (MuPyV), which attaches to host gangliosides with its major capsid protein, VP1. We have solved the crystal structures of VP1 in complex with previously known interaction partners as well as with the ganglioside GT1a, which we have identified as a novel functional receptor for MuPyV. Earlier studies have shown that different strains with singular amino acid exchanges in the receptor binding pocket of VP1 display altered pathogenicity and viral spread. Our investigations show that, while these exchanges do not abolish binding or significantly alter interaction modes to our investigated carbohydrates, they have subtle effects on glycan affinity. The combination of receptor specificity, abundance, and affinity reveals a much more intricate regulation of pathogenicity than previously believed. Our results exemplify how delicate changes to the receptor binding pocket of MuPyV VP1 are able to drastically alter virus behavior. This system provides a unique example to study how the first step in the life cycle of a virus can dictate its biological properties.


Published in the journal: . PLoS Pathog 11(10): e32767. doi:10.1371/journal.ppat.1005104
Category: Research Article
doi: https://doi.org/10.1371/journal.ppat.1005104

Summary

Viruses are obligate intracellular pathogens, and all of them share one crucial step in their life cycle—the attachment to their host cell via cellular receptors, which are usually proteins or carbohydrates. This step is decisive for the selection of target cells and virus entry. In this study, we investigated murine polyomavirus (MuPyV), which attaches to host gangliosides with its major capsid protein, VP1. We have solved the crystal structures of VP1 in complex with previously known interaction partners as well as with the ganglioside GT1a, which we have identified as a novel functional receptor for MuPyV. Earlier studies have shown that different strains with singular amino acid exchanges in the receptor binding pocket of VP1 display altered pathogenicity and viral spread. Our investigations show that, while these exchanges do not abolish binding or significantly alter interaction modes to our investigated carbohydrates, they have subtle effects on glycan affinity. The combination of receptor specificity, abundance, and affinity reveals a much more intricate regulation of pathogenicity than previously believed. Our results exemplify how delicate changes to the receptor binding pocket of MuPyV VP1 are able to drastically alter virus behavior. This system provides a unique example to study how the first step in the life cycle of a virus can dictate its biological properties.

Introduction

The engagement of one or several host cell receptors is the first step in the infectious cycle of a virus. A large number of viruses, including many human pathogens, depend on carbohydrate recognition for initial attachment to the cell surface. Viral tropism and the internalization pathway are usually determined by the specificity and affinity of the receptor interaction as well as the glycan distribution on different cell surfaces (reviewed in [1]). Many viruses use glycoproteins, glycolipids, or both as receptors for cell entry [2]. Gangliosides are ubiquitous glycolipids on the outer leaflet of mammalian cell membranes that serve as receptors for a number of viruses. They are composed of a membrane-embedded ceramide moiety linked to a complex carbohydrate structure that projects away from the cell. Gangliosides almost always contain α-5-N-acetyl-neuraminic acid (sialic acid, Neu5Ac) that can be attached to the core of the molecule with [α-2,3], [α-2,6], or [α-2,8] linkages (Fig 1). Gangliosides exist on cell surfaces in complex and poorly understood patterns that are cell type-, age-, and tissue-dependent ([3,4], reviewed in [5]).

Fig. 1. Overview and biosynthetic pathway of the four most prominent ganglioside series.
Overview and biosynthetic pathway of the four most prominent ganglioside series.
The glycan parts of important members are shown for each series. The downstream biosynthetic steps are identical for all members of a row, although they may vary in linkage orientation. The six-membered pyranose rings are numbered counterclockwise, starting from the bottom (C1, except for C2 in Neu5Ac), and the ring oxygen is symbolized with a black dot. Neu5Ac moieties are rearranged for clarity, and all linkages are mediated by O2 or O8. Most of the gangliosides (e.g. LM1) can be further modified, e.g. by fucosylation. Linkages involving Neu5Ac are present in the α conformation, all other linkages are in the β conformation. Boxes represent three distinguishable sialoglycotopes that contain linkages found in GT1a (blue, representative for [α-2,8]), GD1a (green, [α-2,3]), and 3’-6’-iso-LD1 (also referred to as DSLNT, orange, [α-2,6]). The naming is according to the corresponding gangliosides; if possible, the Svennerholm shorthand is used [6466] All biosynthetic routes were verified using the KEGG metabolic pathway database [67]. A prototype glycan that exemplifies the different positions of Gal and Neu5Ac moieties is depicted on the lower right. The glycan portions investigated in this study are highlighted by purple boxes.

Murine Polyomavirus (MuPyV) is a double-stranded DNA virus that can induce tumors in newborn animals. It was long known to engage glycan receptors that contain a minimal motif of sialic acid [α-2,3]-linked to galactose [6,7], and more recently gangliosides GD1a and GT1b were identified as MuPyV receptors [8]. Viral attachment is mediated by the major capsid protein, VP1, which forms pentameric capsomers that assemble into the T = 7d icosahedral capsid of the virus [9,10]. Sialylated oligosaccharide receptors are engaged in a shallow groove on top of VP1 formed by loop structures on the protein surface [1113], similar to other polyomaviruses [1].

MuPyV displays striking differences in pathogenicity and spread among three closely related prototype strains upon infection of newborn virus-free mice. The laboratory-derived RA strain [14] shows limited spread and induces few tumors of strictly mesenchymal origin after a long latency period, while the naturally occurring PTA strain [15,16] has disseminated infection and causes multiple tumors of epithelial and mesenchymal origin within a short time. LID [17,18], another laboratory isolate MuPyV strain, spreads most rapidly, causing early death by damaging host tissues, leading to brain hemorrhages and kidney failure [19]. The differences among the three strains have been mapped to amino acid variations at two positions, 91 and 296, within the receptor-binding region of VP1 [2024]. While RA bears a glycine residue at position 91, this residue is replaced with a glutamate in both PTA and LID. An additional valine-to-alanine exchange at position 296 is present in LID (Table 1). The pathogenicity profile of one strain can be introduced into the other strains by mutating these two residues, confirming that these substitutions are necessary and sufficient to generate a specific phenotype [25]. The same substitutions have also been observed for other strains of MuPyV [21,22]. MuPyV found in feral mice has the VP1 sequence of PTA [26], but the virus is controlled by an intact immune system. As studies of viral spread can be conducted in vivo and virus infectivity can be tested in vitro using ganglioside deficient mouse cells, MuPyV represents an attractive and rare model system to define the relationships between receptor binding and viral spread and tropism.

Tab. 1. Description of the investigated MuPyV strains.
Description of the investigated MuPyV strains.

Crystal structures of the low pathogenicity strain RA have shown how this virus engages 3’-sialyllactose, a short, linear trisaccharide terminating in [α-2,3]-linked sialic acid, as well as an oligosaccharide that additionally contains a second, branching [α-2,6]-linked sialic acid [11,12]. These structures also identified the location of residues 91 and 296 in the carbohydrate-binding region, suggesting that they might modulate interactions of VP1 with its receptors in the higher pathogenicity strains PTA and LID. Modelling suggested that a glutamate side chain at position 91 would lead to electrostatic repulsion of the [α-2,6]-branched sialic acid, thereby preventing binding of such branched structures by either LID or PTA. Branched sugars carrying an [α-2,6]-linked sialic acid could thus act as pseudoreceptors that will not facilitate productive infection but hamper the spread of RA within the host, in contrast to PTA and LID [8,12]. In line with this hypothesis, gangliosides GD1a and GT1b, which do not contain an [α-2,6]-branched sialic acid, have been identified as entry receptors for the PTA [8,16] and RA strain [27] of MuPyV. However, the molecular determinants of GD1a or GT1b receptor interactions with PTA and LID are not understood, because all structural information is limited to date to RA MuPyV.

To define the interactions of the three MuPyV strains with receptors on the cell surface, we have solved high-resolution structures of RA, PTA, and LID VP1 pentamers in complex with three ganglioside glycans that represent common motifs found in members of the four most prominent ganglioside biosynthesis series and that feature [α-2,3]-, [α-2,6]-, and [α-2,8]-linked sialic acids (for carbohydrate structures, nomenclature, and annotations see Fig 1). We have also conducted crystallographic soaking experiments at different ligand concentrations to compare the relative affinities of each of the three strains for their interaction partners. We find that expanding the well-characterized Neu5Ac-[α-2,3]-Gal epitope with a linear [α-2,8]-linked sialic acid (as found for example in GT1a vs. GD1a) leads to additional interactions between carbohydrate and VP1 in all three strains. Consequently, we identify ganglioside GT1a as an infectious receptor for all three strains. Moreover, the branching [α-2,6]-linked sialic acid is close to the strain-defining amino acids, but can be accommodated by all strains, in contrast to the earlier model. However, the amino acid exchanges defining each strain have subtle effects on their affinity for the validated receptor GD1a. Our results exemplify the effect of minimal changes in a binding pocket on the receptor binding properties of a virus.

Results

GT1a, GD1a, and GT1b gangliosides are infectious receptors for MuPyV

Previous efforts to identify receptors for MuPyV used immortalized cell lines, such as Vero or C6 glioma cells that were supplemented with candidate gangliosides before infection [8,28]. We utilized a mouse embryo knock-out fibroblast cell line (Gang-/- MEFs) specifically deficient in ganglioside synthesis and completely resistant to MuPyV infection (S1A Fig and [29]) to test the ability of ganglioside receptors to rescue infection by different strains of MuPyV. Gang-/- MEFs were supplemented with individual gangliosides followed by infection with RA, PTA, and LID MuPyV (Fig 2). Importantly, it should be noted that the three MuPyV strains we used do not have the same particle to PFU ratio. The viruses have been normalized to similar MOIs, but they cannot be quantitatively compared to one another. However, each strain has been normalized to its own infection rate of WT MEFs; therefore, infection rates upon supplementation of gangliosides can be compared within a strain. The previously identified ganglioside receptors GD1a and GT1b [8] rescued RA, PTA, and LID infection of Gang-/- MEFs in a dose-responsive manner. We also analyzed the GT1a ganglioside that had not been previously investigated as a candidate infectious receptor for MuPyV. We found that GT1a, a member of the ganglio-series synthesized from GD1a (Fig 1), also rescued RA, PTA, and LID infection in a dose responsive manner (Fig 2). Moreover, GT1a supplementation of Gang-/- MEFs conferred higher levels of RA, PTA, and LID MuPyV infection than the previously identified receptors GD1a and GT1b. Finally, we tested the ability of the gangliosides GD1b and GM1 to rescue MuPyV infection of Gang-/- MEFs. GD1b and GT1b supplementation has previously been shown to restore BK polyomavirus infection of ganglioside deficient cells [30]; however, GD1b restored little to no MuPyV infection of Gang-/- MEFs. GM1 supplementation has previously been shown to restore infection by SV40 [8]; however, GM1 did not rescue MuPyV infection of Gang-/- MEFs. These data confirm that GT1a is an infectious receptor for all strains of MuPyV.

Fig. 2. GT1a, GT1b, and GD1a supplementation rescues MuPyV infection of Gang-/- MEFs.
GT1a, GT1b, and GD1a supplementation rescues MuPyV infection of Gang-/- MEFs.
Ganglioside knock-out (Gang-/-) MEFs were completely resistant to infection of all strains of MuPyV as shown by the absence of T-antigen positive nuclei at 24 hours post infection (DMSO control). GD1a, GT1b, and GT1a ganglioside supplementation of Gang-/- MEFs restored RA (A), PTA (B), and LID MuPyV (C) infection, while GD1b and GM1 supplementation resulted in little to no infection by any virus strain. Infection levels were quantified at both 2 μM and 4 μM ganglioside supplementation (blue and green bars, respectively). Infection levels are normalized to MuPyV infection of WT MEFs, and error bars correspond to standard error.

We also investigated whether MuPyV cell surface binding to infectious or non-infectious ganglioside receptors correlated with infection. To this end, we measured the levels of free (unbound) virus in each ganglioside supplemented sample at 4 hours post infection. We did not detect significant differences in MuPyV cell surface binding to different ganglioside receptors or WT MEFs, indicating that cell surface binding alone does not determine infection (S1B Fig). Instead, a considerable amount of virus binds to Gang-/- MEFs even in the absence of ganglioside supplementation (S1A Fig). MuPyV is also endocytosed in Gang-/- MEFs, which however does not lead to infection [29]. Taken together, these data confirm that gangliosides are not required for cell surface binding. They are, however, required for infection, and GT1a appears to be more efficient than GD1a and GT1b.

Structure of MuPyV VP1 bound to GT1a

In order to define the mode of recognition of GT1a, particularly to the naturally occurring PTA strain of MuPyV, we have soaked VP1 crystals with the glycan portion of GT1a and solved the structure of the complex (Table 2). While the receptor interaction pocket of RA VP1 has been described [1113], no structural information for the pathogenicity-defining amino acids at positions 91 and 296 in the pockets of PTA and LID has been available. PTA and LID both carry a glutamate at position 91, and this side chain is being held in a characteristic position with the carboxyl group facing away from the glycan receptor due to a salt bridge formed with K186 (Fig 3), as previously predicted [12]. The GT1a glycan is a branched structure with a long disialylated arm, which has the sequence Neu5Acb-[α-2,8]-Neu5Aca-[α-2,3]-Gala-[β-1,3]-GalNAc, and a second short arm, which consists of a single Neu5Acd [α-2,3]-linked to Galb (for carbohydrate structures, nomenclature, and moiety indexing see Fig 1). The disialylated arm of GT1a is clearly visible in the crystal structure of PTA VP1; it is well defined by electron density and makes extensive contacts with the protein (Fig 4B–4D). Overall, the GT1a glycan adopts a twisted horseshoe-like shape, with Neu5Aca and Neu5Acb wrapping around the side chains of Y72 and R77 of VP1. Its longer, disialylated arm contains a Neu5Aca-[α-2,3]-Gala sequence that is also present in GD1a and simpler compounds such as 3’-sialyllactose (3SL), and the interactions of this motif with VP1 are essentially identical to those seen in previous structures [1113]. However, our structure visualizes an additional network of contacts made by the terminal [α-2,8]-linked Neu5Acb (Fig 4C and 4D). Its carboxyl group engages Y72 and forms water-mediated hydrogen bonds with Q71, Y72, as well as D85 of the neighboring monomer (D85*). In addition, the N-acetyl nitrogen of Neu5Acb forms a hydrogen bond with the backbone carbonyl of T67, and O8 and O9 in the glycerol chain of the sugar are hydrogen-bonded with the R77 side chain. The carboxyl groups of Neu5Aca and Neu5Acb are about 4 Å apart, and the positively charged side chain of R77 counteracts their negative charges (Fig 4C and 4D). Neu5Aca and Neu5Acb contribute binding interfaces of approximately 160 Å2 and 190 Å2, respectively (calculated using the PISA server [31]). The remaining Gala-GalNAc-Galb stem of GT1a forms fewer contacts with the protein, which include a hydrogen bond between G78 and the Gala O4 hydroxyl group (Fig 4) as well as several van der Waals interactions. Notably, the Cβ and Cγ atoms of E91 are within van-der-Waals range of O6 and C6 of Gala, and the E91 carboxylate group is close to C6 of GalNAc. The total contact surface for this portion of the glycan is 142 Å2.

Fig. 3. The MuPyV binding pocket.
The MuPyV binding pocket.
Top view on the receptor-binding region of PTA, which is shown with E91 and V296 highlighted in salmon. Residues that are known to participate in receptor binding are contributed by the BC and HI loops and are highlighted as stick models. One monomer is shaded in green and the other monomers are alternatingly shaded light and dark grey for better distinction.

Fig. 4. Binding of GT1a to PTA.
Binding of GT1a to PTA.
A The PTA binding pocket and the GT1a conformation upon binding are shown from an angle parallel to the fivefold axis. A scheme of the glycan is shown in the inset, and the sugar rings are filled according to the coloring scheme from Fig 1. B Simulated annealing Fobs-Fcalc omit map (resolution 1.71 Å, calculated at 3.5 σ, carved 2.3 Å around the glycan). C Possible binding interactions of GT1a and PTA. E91 and V296 are highlighted in salmon. Hydrogen bonds are shown in black, the hydrophobic contact mediated by V296 in gold, and the van-der-Waals contacts of E91 are shown in cyan. Waters that mediate key hydrogen bonds are shown as red spheres. Unique interactions mediated by the novel GT1a-like binding motif are shown in red. D Zoomed view of the binding to the two terminal Neu5Ac moieties. The rest of the glycan is omitted for clarity. Residues except Y72 and R77 as well as waters involved in contacts with these two glycan moieties are pale grey and salmon, respectively.

Tab. 2. Data collection and refinement statistics.
Data collection and refinement statistics.

Because the differences in tumorigenicity and host spread among strains have been mapped to the glycan binding pocket of VP1, and because GT1a appears to be particularly efficient in facilitating productive infection, we set out to determine how the three strains engage GT1a. By solving the crystal structures of RA and LID VP1 complexed with GT1a using the identical strategy used for the PTA-GT1a complex, we found that the overall binding mode of GT1a is very similar across the three strains (Fig 5A), with a conserved binding mode of the [α-2,8]-linked Neu5Acb. Although the replacement of glutamate with glycine at position 91 leads to a contact area decrease of 33 Å2 in RA, the orientation of GT1a in this strain is not altered (compare Fig 5B and 5C). Likewise, the substitution of valine with alanine at position 296 in LID removes a hydrophobic contact but does not affect the conformation of GT1a (Fig 5E; S2 Fig).

Fig. 5. Binding modes of GT1a to the different MuPyV strains.
Binding modes of GT1a to the different MuPyV strains.
A Superposition of the GT1a-binding mode of RA (GT1a in sky blue), PTA (dark blue), and LID (pale blue). The Neu5Acb-[α-2,8]-Neu5Aca-[α-2,3]-Gala motif is shown in solid sticks, together with the adjacent GalNAc moiety. All superpositionings were carried out in PyMOL [61] using ‘align’ for the protein chains only. Surface, E91 and V296 are from PTA/GT1a. All ‘align’ rmsd values are below 0.16 Å. B & C Close view of the van-der-Waals contacts introduced by the E91 side chain present in PTA and LID (C), but not in RA (B). Hydrogen bonds are shown in grey, van-der-Waals contacts in cyan. D & E Close view of the hydrophobic contact mediated by V296 in RA and PTA (D), but not by A296 in LID (E). The 4.0 Å hydrophobic contact is not present in the LID strain, whose pocket is opened to the right. Hydrogen bonds are shown in grey, hydrophobic contacts are shown in gold.

The Neu5Aca-Gala-GalNAc linkages in the long arm of GT1a adopt conformations that have been reported in numerous structures (for example [3234]). While the [α-2,3] linkage between Neu5Aca and Gala adopts the conformation that has been reported for DSLNT and 3SL, the branching Neu5Acd-[α-2,3]-Galb linkage adopts a different conformation, which has been reported for structures containing O-4-substituted galactoses (as in [35,36]). While a higher variability is observed for Neu5Ac-[α-2,8]-Neu5Ac linkages (S2E Fig), this linkage adopts torsion angles that are in agreement with other, related structures such as in the structure of human liver fructose-1,6-bisphosphatase in complex with an allosteric inhibitor [37] or in the complex of tetanus toxin with a GT1b analog [38]. The overall structure is in good agreement with a molecular dynamics simulation using an AMBER force field in an aqueous environment [39]. A well-defined set of water molecules mediates bridged hydrogen bonds between the pyranose moieties, especially between Neu5Acb and Neu5Acd (S3 Fig). Due to these steric constraints, the GT1a complexes feature well-defined electron density not only for the binding epitope, but also for the non-binding, branching NeuNAcd in its preferred solution conformation [40], which brings this moiety to about 5 Å near the end of the long arm and gives the glycan the characteristic, horseshoe-like topology that is observable in all complex structures.

Structures of MuPyV VP1 strains bound to other sialylated glycans

As RA, PTA, and LID VP1 all bind GT1a in a highly similar conformation, we hypothesized that the differences in pathogenicity and spread among the three strains might be due to the recognition of additional carbohydrates by only a subset of MuPyV strains. As shown in Fig 1, the many different gangliosides share a relatively small set of common sialoglycotopes. We therefore investigated the ability of all three VP1 proteins to bind other glycan structures that are representative for these epitopes. We solved structures of VP1 bound to the glycan portions of two of these gangliosides: The GD1a glycan is an established infectious receptor and essentially a truncated version of GT1a lacking the [α-2,8]-linked Neu5Acb in the long arm. The human milk hexasaccharide DSLNT is the glycan portion of the lacto-series ganglioside 3’-6’-isoLD1 (Fig 1) [41], which is overexpressed in the central nervous system. In contrast to GT1a and GD1a, DSLNT does not contain an [α-2,3]-linked Neu5Acd as a short arm but instead a branching [α-2,6]-linked Neu5Acc. This structure is similar to a very common epitope on O-linked glycoproteins [4244]. DSLNT was used in previous studies of MuPyV as a model “pseudoreceptor” [12] and was investigated here to help rationalize these earlier data, to facilitate a comparison among strains, and to establish a binding profile for glycans containing an [α-2,6]-linked sialic acid.

GD1a

The previously identified MuPyV receptor GD1a is similar to a truncated GT1a structure, containing only a Neu5Aca-[α-2,3]-Gala motif instead of Neu5Acb-[α-2,8]-Neu5Aca-[α-2,3]-Gal. The disaccharide engages all three strains in a very similar manner (Fig 6A). Neither the longer E91 side chain (in PTA and LID) nor the shorter A296 side chain (in LID) result in an altered conformation of the ligand.

Fig. 6. Differences in receptor binding patterns across glycans.
Differences in receptor binding patterns across glycans.
A Superposition of the binding modes of GD1a to RA (light green), PTA (dark green), and LID (pale green). The sequence of GD1a is shown in the inset. The Neu5Aca-[α-2,3]-Gala motif is shown in solid sticks, together with the adjacent GalNAc moiety. In all figures, GT1a bound to PTA is overlaid as a grey ghost for comparison, with Neu5Acb omitted for clarity. Deviations exceeding the atomic error of the structure and alignment rmsd values are only found in the stem region of the sugar, starting at Galb. All superpositionings were carried out in PyMOL [61] using ‘align’ for the protein chains only. Surface, E91 and V296 are from PTA/GT1a. All ‘align’ rmsd values are below 0.16 Å. B & C Comparison of the DSLNT binding modes to RA (red), PTA (orange), and LID (violet). In PTA-DSLNT, [α-2,6]-branching causes a 15–20° psi angle shift of the GlcNAc moiety compared to GD1a and GT1a, resulting in a 1 Å sideward twist movement of the stem. In RA-DSLNT, combination of this shift with a 15° shift in the phi angle results in a downward movement of GlcNAc and its branching Neu5Acc compared to PTA-DSLNT. In LID, the shift is already observable for Gala and results in the loss of ordered density for GlcNAc. All angles were analyzed with CARP. The sequence of DSLNT is shown in the inset of panel C.

DSLNT

The DSLNT glycan terminates in a Neu5Aca-[α-2,3]-Gala motif, which is the part of the molecule best defined by electron density in all complexes. DSLNT also contains an additional [α-2,6]-linked, branched Neu5Acc residue, which is not present in either GT1a or GD1a. There is weak electron density for Neu5Acc in one of the five binding pockets of the RA strain, but it only engages in few interactions [12]. While PTA and LID do bind DSLNT, the complex structures do not show any electron density for the Neu5Acc, indicating that this sugar is conformationally flexible and does not contribute contacts. When bound to the PTA strain, the stem of DSLNT is moderately rearranged (Fig 6B). In comparison to GalNAc in GT1a, the GlcNAc moiety is slightly tilted away from E91 due to a ~20° rotation of the psi angle in the Gala-[β-1,3]-GlcNAc linkage (Fig 6B and 6C, assessed using CARP) that propagates throughout the sugar. In addition, there is no visible electron density for the GlcNAc O6 that is engaged in the [α-2,6]-branching as well as an increased B-Factor variance within the glycan (S4 Fig). The reason for the sideward twist and the missing electron density for Neu5Acc observed in PTA is probably the electrostatic repulsion between the carboxyl groups of Neu5Acc and E91. While the charge of E91 is compensated by K186 (Fig 3), as was hypothesized before [12], the two carboxylate groups would come within about 2 Å of one another if DSLNT bound to PTA in the same way as observed in RA. This hypothesis is confirmed by a PTA E91Q mutant that rescued binding of Neu5Acc (S5 Fig). In turn, when bound to RA, DSLNT exhibits a stronger conformational rearrangement (Fig 6B and 6C). Due to the missing side chain at position 91, the psi angle rotation between Gala and GlcNAc is accompanied by an additional 15° rotation of the phi angle, bringing the GlcNAc moiety and the attached Neu5Acc closer to the protein surface [12]. In LID, the valine to alanine mutation at position 296 reduces its van-der-Waals radius. This change results in a broader binding pocket compared to the other strains and the loss of a hydrophobic interaction between position 296 in VP1 and C3 of Neu5Aca for all glycans. This gives room for a stark alteration in the binding mode of DSLNT that starts with a slight tilt of Neu5Aca and propagates through the sugar (Fig 6C), ultimately resulting in a prominent sideward shift of the whole glycan stem. The resulting increased conformational freedom of DSLNT is reflected by a lack of electron density in its stem region as well as by an elevated temperature factor variance (S4 Fig). This alteration of the binding mode in LID is likely to be observed for other glycans that are not conformationally restrained by the [α-2,3]-linked Neu5Acd.

Relative affinities of MuPyV strains for sialylated glycans

Since all three MuPyV strains are able to engage the three different glycan structures in a largely identical manner, we reasoned that the differences in pathogenicity and spread might be attributable to subtle differences in affinity, rather than specificity, among the strains. The affinities of RA VP1 for 3’-sialyllactose and DSLNT were estimated to be in the low mM range [11]. Coupled with the high costs of glycans and the high amount required due to their low binding affinity, weak binding poses technical challenges for classical affinity measurements. We therefore utilized a crystallographic approach to quantitatively compare ligand binding. We crystallized all three VP1 pentamers in the same condition, and soaked each with the oligosaccharide portions of GT1a, GD1a, and DSLNT at different concentrations in parallel. X-ray data of all crystals were collected in the same manner, and the data sets were processed using the same protocol and integrated as described previously [45]. All data sets were processed in the same unit cell, scaled, and the bias-reduced difference electron density around the central Neu5Aca-[α-2,3]-Gala motif was quantified for each data set (see the Methods section for details). Our crystallization condition contains a high amount of ammonium sulfate, which competes with the carboxyl group of Neu5Aca and has to be displaced by the carbohydrates. Therefore, our observed binding is weaker than in a physiological setting. However, while not yielding dissociation constants in the traditional sense, this method enables us to compare relative levels of binding across our three different strains and three different glycans.

The GT1a glycan exhibits the strongest binding in all three VP1 variants compared with DSLNT or GD1a (Fig 7A–7C), with no detectable difference between the strains (Fig 7D). This finding is in accord with our ganglioside add-back experiments in cell culture (Fig 2), which consistently showed higher levels of infection mediated by GT1a compared to GD1a. The stronger overall binding of GT1a can be attributed to the additional [α-2,8]-linked sialic acid present in GT1a (Neu5Acb), which contributes several interactions and an increased buried surface area. These contacts seem to outweigh the differences in van der Waals contacts with the side chains of E91 or V296, at least to the extent discernable in our assay.

Fig. 7. Binding of ligands to MuPyV VP1.
Binding of ligands to MuPyV VP1.
The average simulated annealing Fobs-Fcalc electron density for the Neu5Aca-[α-2,3]-Gala in GT1a, GD1a, and DSLNT is plotted against ligand concentration. RA VP1 is colored blue, PTA VP1, is colored red, and LID VP1 is colored green. GT1a is displayed in a dotted line with triangles, GD1a in a dashed line with circles, and DSLNT in a solid line with squares. The error bars correspond to the standard deviation of the mean electron density observed in the five chains of VP1. A Electron density of GT1a, GD1a, and DSLNT in RA VP1. B Same as in A, but for PTA VP1. C Same as in A, but for LID VP1. D Comparison of GT1a-derived electron density in RA VP1, PTA VP1, and LID VP1. E Same as in D, but for GD1a. F Same as in D, but for DSLNT.

GD1a binds less well to all strains compared to GT1a. In addition, there are differences in binding strength among the three strains. PTA and LID VP1 appear to bind GD1a at the same level and better compared with RA (Fig 7A–7C and 7E) because these two strains gain additional interaction surface and van-der-Waals contacts from their E91 side chain. This effect is more pronounced than in GT1a, because in GD1a it cannot be masked by the additional contacts of the [α-2,8]-linked Neu5Acb.

DSLNT displays the lowest overall affinity to all strains, with levels comparable to GD1a in RA for all three strains (Fig 7A–7C) despite the DSLNT conformation being slightly different in each VP1 complex (Fig 7F). Neither the blocking of Neu5Acc binding by E91, nor the increased conformational freedom in LID appears to alter binding affinity. It is possible that Neu5Acc in RA adopts a conformation that might not be favorable and therefore not heavily contribute to affinity, in spite of the added contact surface. Combined with the fact that electron density for Neu5Acc could only be observed in one binding pocket of RA VP1 [12], we believe that this conformation is possible but not probable in solution. Instead, an increased number of conformational options might make up for a loss of binding contacts.

Discussion

Many viruses engage cell-surface glycans to mount an infection, and subtle differences in the recognition of such receptors can be linked with altered tropism and pathogenicity. Examples include the canine parvovirus and feline panleukopenia virus [46,47], the human BK polyomavirus [48], B-lymphotropic polyomavirus [49,50] as well as avian and human influenzaviruses [51,52]. However, MuPyV is a rare example of a virus in which drastic differences in pathogenicity directly correlate with single amino acid substitutions in the viral capsid.

In order to provide a structural basis for understanding the profoundly different pathogenicities of the three MuPyV strains RA, PTA and LID, we have solved structures of their VP1 proteins and characterized their receptor-binding properties. We show that the ganglioside GT1a serves as a MuPyV receptor and promotes infection with higher potency than the previously identified receptors GD1a and GT1b. Structurally, the increased potency of GT1a can be directly explained by a set of additional contacts involving the [α-2,8]-linked Neu5Acb that is only present in this glycan and that gives it a characteristic horseshoe-like shape. It had previously been suggested that the G91E mutation present in PTA and LID abolishes binding to branched glycans containing [α-2,6]-linked Neu5Ac and thus allows the virus to spread more efficiently in the host [8,11]. However, our analyses show that the presence of a glutamate at position 91 still allows binding of the branched oligosaccharides GT1a, GD1a, or DSLNT to all three strains, albeit with subtle differences in binding affinity. While all three strains bind GT1a with comparable affinity, PTA and LID bind GD1a better than RA. The DSLNT glycan binds similarly to all three strains, with the lowest overall affinity. This is again in line with the structures, which show that the branched Neu5Acc of DSLNT does not engage in any specific contacts. The limited contacts between Neu5Acc and RA observed in an earlier structure [12] have to be considered a crystallization artifact as they were only observed in one out of five binding sites, and this visible Neu5Acc moiety was located near a crystal contact.

The ligand binding promiscuity of MuPyV is surprisingly high. Binding mostly requires a ubiquitous minimal Neu5Ac-[α-2,3]-Gal motif, in agreement with earlier findings [6,7]. It therefore seems plausible that the virus also recognizes other glycans bearing this motif, resulting in differences in pathogenicity and spread. Preliminary studies show that glycans with an N-acetyllactose core (Neu5Ac-[α-2,3]-Gal-[β-1,4]-GlcNAc), as found in neolacto gangliosides such as the predominant ganglioside of peripheral nerve cells, LM1 [53,54], can also be bound in a manner similar to DSLNT and with higher flexibility than GT1a or GD1a (S6 Fig).

Based on our structures, certain requirements that contribute to receptor specificity can be established. For example, branches at Gal-O4 within the minimal motif produce clashes and cannot be tolerated. Therefore, although the GD1a glycan possesses two Neu5Ac-[α-2,3]-Gal motifs, it prefers the one on its longer arm for complex formation. For the same reason, glycans such as GM1 or GM2 that only possess such a branched Neu5Ac-[α-2,3]-Gal epitope cannot engage MuPyV productively. In support of this, the GM1 ganglioside is not able to rescue MuPyV infection of Gang-/- MEFs (Fig 2, [29]), although low-level and probably non-specific interactions with cells can be detected (S1B Fig). GT1b possesses a disialylated arm at Galb and is monosialylated at Gala. We predict that GT1b engages VP1 with its monosialylated arm. The second, disialylated arm is likely to be accommodated in such a binding mode, and the [α-2,8]-linked sialic acid might contribute additional contacts. Binding via the monosialylated arm is in line with our findings that supplementation of Gang-/- cells with GT1b rescues infection at a level between GD1a and GT1a. Some gangliosides whose glycan epitopes are capable of engaging VP1 in vitro might not be infectious receptors in vivo, mainly because of steric complications in the context of the cell membrane. For example, while the crystal structure of PTA with the glycan portion of GD3 shows an identical binding mechanism to GT1a (S7 Fig), supplementation of Gang-/- MEFs with GD3 does not restore infectivity [29]. We reason that the glycan stem of GD3 (and of gangliosides with a similar length such as GM3) is too short to allow efficient attachment of the MuPyV capsid to the cell membrane.

The discrepancy in pathogenicity in MuPyV strains that differ only at one single position is stark. In sharp contrast, the differences among receptor binding between the three strains investigated here are subtle, and a correlation of the structural data with the observed pathogenicity profiles remains challenging. One reason for this is that avidity effects in the virus capsid, which can engage many ligands simultaneously, multiply subtle changes in the affinity of capsomers for single glycans. It was shown for influenza viruses that small changes between millimolar binding affinities of single binding sites can result in dramatically altered viral binding properties [52]. As discussed above, we found the main difference between RA and PTA/LID to be a differing affinity for GD1a, which appears to bind better to the latter strains due to the larger E91 side chain. This might facilitate attachment and productive infection by these strains to cells that display GD1a, and may thus give them an advantage over RA. While we could not show differences between the PTA and LID strain in terms of glycan affinity to isolated VP1 pentamers, it is unclear how this translates to avidity effects. As such, it is possible that capsid avidities differ enough to explain the more limited spread of PTA. Although direct correlations cannot be made, it becomes increasingly clear that the virus needs to uphold a delicate equilibrium between efficient infection and release from infected and lysed cells as well as selective affinity for productive receptors. The absence of the RA and LID strains outside the laboratory [26] emphasizes that this equilibrium is affected by minute changes in the receptor binding properties.

The MuPyV receptor pocket can clearly accommodate several related but distinct glycan structures (Figs 1 and 46). These structures also decorate glycoproteins on many cell surfaces. It therefore seems likely that MuPyV can also engage glycans that are not attached to gangliosides. For instance, the glycan stem of GD1α, which is very similar to DSLNT and prominent on glycoproteins [4244], is a likely receptor candidate. The different cell-surface distribution patterns of glycoproteins and gangliosides may likewise influence MuPyV spread [8]. Glycoprotein receptors with unknown identity have in fact been shown to promote non-productive internalization of MuPyV, which in turn elicits innate immune responses by the host [29]. Along these lines, our results suggest that virus particles adhere to and enter ganglioside deficient MEFs to levels that are not significantly lower than for wild-type and ganglioside supplemented Gang-/- cells, although without detectable infection. Although not representative for other cell types, these results suggest that the amount of non-productive “pseudoreceptors” on the MEF cell surface is much higher than anticipated.

Our data demonstrate that varying affinities for different gangliosides are the key determinants of a successful MuPyV infection, in line with earlier reports [68]. Perhaps unexpectedly, we also find that (even non-specific) attachment of the virus to a host cell can lead to successful internalization, but that this does not necessarily lead to an infection. Thus, we propose that the ratio between productive (ganglioside bound) and non-productive (ganglioside and glycoprotein bound) glycotopes on the host cell itself or in its microenvironment helps to determine the productivity of infection through diverging entry routes, and that differential affinities to these receptors dictate this equilibrium. The nature of these diverging routes, their underlying driving forces, and potential biological consequences other than immune stimulation [29] remain unknown–as does the point at which they diverge. We cannot exclude the possibility that the distribution and binding properties of (pseudo-)receptors are of importance mostly for the post-entry stage rather than for events taking place at the cell surface. A better understanding of the distribution patterns and densities of glycans on specific cells is clearly needed to fully appreciate the many aspects of pathogenesis and tropism of MuPyV as well as many other glycan-binding viruses.

Methods

Ganglioside supplementation and quantification of MuPyV infection

WT and Gang-/- MEFs were seeded onto 96-well Costar 3906 imaging plates in Dulbecco's Modified Eagle's Medium supplemented with 10% fetal bovine serum (FBS). WT (B4+/+St8+/+) and Gang-/- MEFs (B4-/-St8-/-) were provided by Thomas Benjamin at Harvard Medical School. Gangliosides were purchased from Matreya LLC and resuspended in DMSO upon arrival, aliquoted, and stored at -20°C until use. Cells were incubated overnight in serum free media prior to infection. For ganglioside supplemented Gang-/- MEFs, cells were starved in serum free media containing the indicated concentration of ganglioside. Gangliosides were then removed, and cells were washed with serum free media to remove any free ganglioside. Cells were then infected with NG59RA, PTA, and LID MuPyV (MOI ~10–30). At 24 hours post infection cells were washed in phosphate buffered saline and fixed with 4% paraformaldehyde at room temperature for 10 minutes. Cells were then permeabilized with 0.1% Triton X-100, blocked in 10% FBS in PBS, and then stained for the viral protein, T-antigen (E1). Samples were then incubated with Alexa Fluor labeled secondary antibodies (546). Plates were imaged with the Molecular Devices ImageXpress Micro XL High-Content Screener. The percent infected was calculated for each well (5 images were taken per well). Three wells were quantified per sample and the average percent infected, standard error, and standard deviation were calculated for each sample. To quantify infection, T-antigen staining was measured per each DAPI labeled nucleus. For image analysis, the DAPI channel on each image was thresholded, and nuclei were counted using ImageJ (Analyze Particles). These particles were marked as “Regions of Interest” (ROI), and then the average pixel intensity of T-antigen staining was measured for each nucleus (ROI). These were then binned into T-antigen positive or T-antigen negative nuclei to create % infected.

VP1 immunofluorescence staining

WT and Gang-/- MEFs were seeded onto glass coverslips in Dulbecco's Modified Eagle's Medium supplemented with 10% (FBS). Cells were incubated overnight in serum free media prior to infection. For ganglioside supplementation, Gang-/- MEFs were starved in serum free media containing the indicated concentration of ganglioside. Gangliosides were then removed and cells were washed with serum free media to remove any free ganglioside. Cells were then infected with NG59RA. At indicated times post infection the cells were fixed with 4% paraformaldehyde at room temperature. Cells were blocked in 10% FBS in PBS and then stained for GD1a using mAb MAB5606 (Millipore). Cells were then permeabilized with 0.1% Triton X-100 and stained for the viral proteins, VP1 (I58 antibody) and T-antigen (E1 antibody). Samples were washed and then incubated with Alexa Fluor labeled secondary antibodies (488, 546, 647). Slides were then mounted using DAPI prolong gold mounting media. Slides were imaged with a Nikon A1R confocal microscope. All images were taken as a 9 to 13 step (.25μm) z-stacks on a laser scanning confocal microscope. Each z-stack was aligned and compressed into a max intensity Z projection image.

Virus binding to ganglioside supplemented Gang-/- MEFs

WT and Gang-/- MEFs were seeded onto a 24 well dish in Dulbecco's Modified Eagle's Medium supplemented with 10% (FBS). Cells were incubated overnight in serum free media prior to infection. For ganglioside supplemented Gang-/- MEFs, cells were starved in serum free media containing the indicated concentration of ganglioside. Gangliosides were then removed and cells were washed with serum free media to remove any free ganglioside. Cells were then infected with either NG59RA, PTA, or LID at an MOI ~10–30 (250 μL/well). At 4 hours post infection 150 μL of virus supernatant was removed and placed into a microcentrifuge tube. This virus supernatant was then used to infect WT MEFs seeded onto a 96-well plate (50 μL/well). The amount of free virus was then quantified as percent of infection of the 96-well reinfection plate. At 24 hours post virus addition the plate was washed in PBS and fixed with 4% PFA at RT for 10 minutes. Cells were then permeabilized with 0.1% Triton X-100, blocked in 10% FBS in PBS, and then stained for the viral protein, T-antigen (E1). Samples were then incubated with Alexa Fluor labeled secondary antibodies (546). Plates were imaged with the Molecular Devices ImageXpress Micro XL High-Content Screener. The percent infected was calculated for each well (5 images were taken per well) as indicated by T-antigen positive nuclei. Three wells were quantified per sample and the average percent infected, standard error, and standard deviation were calculated for each sample. For image analysis, the DAPI channel on each image was thresholded and nuclei were counted using ImageJ (Analyze Particles). These particles were marked as “Regions of Interest” (ROI) and then the average pixel intensity of T-antigen staining was measured for each nucleus (ROI). These were then binned into T-antigen positive or T-antigen negative nuclei to create % infected.

Expression and purification of VP1 pentamers

DNA encoding residues 33–316 of RA (GenBank # M34958.1) or PTA VP1 (GenBank # PSU27812) was cloned into the expression vector pET15b (Novagen) in frame with an N-terminal hexahistidine tag (His-tag) and a thrombin cleavage site. DNA for LID VP1 (GenBank # PSU27813) was generated by site-directed mutagenesis of PTA VP1 residue 296. VP1 pentamers were overexpressed in E. coli (BL21) after IPTG induction, and purified by nickel affinity chromatography. The His-tag was removed by thrombin cleavage on column for 72 hours (leaving the non-native residues GSHM at the N-terminus), followed by size exclusion chromatography on a Superdex-200 column.

Crystallization and crystal soaking

Pure VP1 pentamers were supplemented with 20 mM DTT, concentrated to 7.5–8 mg/mL (RA VP1) or 8.5–9 mg/mL (PTA and LID VP1), and crystallized by sitting-drop vapor diffusion. RA VP1 was crystallized at 20°C against reservoir solutions containing a range of 1.25–1.8 M ammonium sulfate and 1–10% (v/v) isopropanol. PTA and LID were crystallized at 4°C against reservoir solutions containing 0.1 M HEPES pH 7–8.5 and 1.6–1.8 M K-Na phosphate. For complex formation, the crystals were soaked in the reservoir solution supplemented with the glycan. The detailed crystallization and soaking procedures are listed in S1 Table. The GT1a and GD1a glycans were purchased from Elicityl SA (France), and the DSLNT glycan was purchased from Carbosynth (United Kingdom).

For concentration-dependent soaking VP1 pentamers of all three strains were crystallized at 20°C against a mother liquor containing 1.5 M ammonium sulfate and 6% (v/v) isopropanol. These crystals were soaked in drops of mother liquor containing the appropriate concentration of glycan for 30 minutes.

All crystals were cryoprotected by incubation in mother liquor supplemented with the appropriate concentration of glycan and 25% (v/v) glycerol. They were then flash-frozen in liquid nitrogen.

Structure determination and electron density quantification

Data reduction was carried out in XDS [55], and the structure of native RA VP1 was solved in Molrep [56] using a model generated from the previously solved structure of P16 VP1 (PDB code 1VPN [12]). Other structures were solved by molecular replacement using the RA VP1 structure in Phenix [57]. All structures were completed by alternating rounds of manual model building in Coot [58], followed by restrained coordinate and isomorphous B-factor refinement including TLS refinement and five-fold non-crystallographic symmetry restraints in Refmac5 [59]. TLS parameters were obtained from the TLSMD server [60]. All models agree well with the experimental data and have good geometry (Table 2). The PDB accession codes for the structures are listed in Table 2. Structural figures were prepared in PyMOL [61].

Data collected for concentration-dependent soaking experiments was processed as described above. The unit cell parameters for all datasets were treated as equal for all datasets and isomorphous to the dataset “RA Nat” (S2 Table). They were scaled against “RA Nat” in Scaleit [62] and then subjected to B-factor refinement and simulated annealing in Phenix against models of RA, PTA, or LID VP1, which lacked atoms of all solvent molecules in the receptor binding pocket as well as those of tryptophan residues 98 and 227 as controls. The resulting bias-reduced Fobs-Fcalc electron density for Neu5Aca-[α-2,3]-Gala and the two marker tryptophans was calculated as a summation of values of the grid points in a mask generated 1 Å around these groups using the program Mapman [63]. The overall binding of a sugar at different concentrations influences the electron density of the Neu5Aca-[α-2,3]-Gala portion, which is included in GT1a, GD1a, and DSLNT. In contrast, it has no effect on the electron density of the marker tryptophan residues, which do not differ significantly for all data points. For each data point, the average density of the five chains was plotted against ligand concentration and submitted to a non-linear least squares fit using the equation

where Bmax was the highest observed electron density value overall (constrained to 95.03 AU) and B0 the electron density in the binding pocket at 0 mM ligand concentration. Plotting and fitting was done using the program Prism 6 (GraphPad Software, Inc., La Jolla, California, USA).

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9


Zdroje

1. Ströh LJ, Stehle T (2014) Glycan Engagement by Viruses: Receptor Switches and Specificity. Annu Rev Virol 1: 285–306.

2. Matrosovich M, Herrler G, Klenk HD (2015) Sialic Acid Receptors of Viruses. Top Curr Chem. 367:1–28 doi: 10.1007/128_2013_466 23873408

3. Ngamukote S, Yanagisawa M, Ariga T, Ando S, Yu RK (2007) Developmental changes of glycosphingolipids and expression of glycogenes in mouse brains. J Neurochem 103: 2327–2341. 17883393

4. Yu RK, Macala LJ, Taki T, Weinfield HM, Yu FS (1988) Developmental changes in ganglioside composition and synthesis in embryonic rat brain. J Neurochem 50: 1825–1829. 3131485

5. Yu RK, Tsai YT, Ariga T, Yanagisawa M (2011) Structures, biosynthesis, and functions of gangliosides—an overview. J Oleo Sci 60: 537–544. 21937853

6. Cahan LD, Singh R, Paulson JC (1983) Sialyloligosaccharide receptors of binding variants of polyoma virus. Virology 130: 281–289. 6316632

7. Fried H, Cahan LD, Paulson JC (1981) Polyoma virus recognizes specific sialyligosaccharide receptors on host cells. Virology 109: 188–192. 6258307

8. Tsai B, Gilbert JM, Stehle T, Lencer W, Benjamin TL, et al. (2003) Gangliosides are receptors for murine polyoma virus and SV40. EMBO J 22: 4346–4355. 12941687

9. Liddington RC, Yan Y, Moulai J, Sahli R, Benjamin TL, et al. (1991) Structure of simian virus 40 at 3.8-A resolution. Nature 354: 278–284. 1659663

10. Stehle T, Gamblin SJ, Yan Y, Harrison SC (1996) The structure of simian virus 40 refined at 3.1 A resolution. Structure 4: 165–182. 8805523

11. Stehle T, Harrison SC (1996) Crystal structures of murine polyomavirus in complex with straight-chain and branched-chain sialyloligosaccharide receptor fragments. Structure 4: 183–194. 8805524

12. Stehle T, Harrison SC (1997) High-resolution structure of a polyomavirus VP1-oligosaccharide complex: implications for assembly and receptor binding. EMBO J 16: 5139–5148. 9305654

13. Stehle T, Yan Y, Benjamin TL, Harrison SC (1994) Structure of murine polyomavirus complexed with an oligosaccharide receptor fragment. Nature 369: 160–163. 8177322

14. Feunteun J, Sompayrac L, Fluck M, Benjamin T (1976) Localization of gene functions in polyoma virus DNA. Proc Natl Acad Sci U S A 73: 4169–4173. 186787

15. Dawe CJ, Law LW, Dunn TB (1959) Studies of parotid-tumor agent in cultures of leukemic tissues of mice. J Natl Cancer Inst 23: 717–797. 13814555

16. Gross L (1951) "Spontaneous" leukemia developing in C3H mice following inoculation in infancy, with AK-leukemic extracts, or AK-embrvos. Proc Soc Exp Biol Med 76: 27–32. 14816382

17. Main JH, Dawe CJ (1966) Tumor induction in transplanted tooth buds infected with polyoma virus. J Natl Cancer Inst 36: 1121–1136. 4287623

18. Rowe WP, Hartley JW, Estes JD, Huebner RJ (1959) Studies of mouse polyoma virus infection. 1. Procedures for quantitation and detection of virus. J Exp Med 109: 379–391. 13641563

19. Bolen JB, Fisher SE, Chowdhury K, Shan TC, Williams JE, et al. (1985) A determinant of polyomavirus virulence enhances virus growth in cells of renal origin. J Virol 53: 335–339. 2981359

20. Bauer PH, Bronson RT, Fung SC, Freund R, Stehle T, et al. (1995) Genetic and structural analysis of a virulence determinant in polyomavirus VP1. J Virol 69: 7925–7931. 7494305

21. Dawe CJ, Freund R, Mandel G, Ballmer-Hofer K, Talmage DA, et al. (1987) Variations in polyoma virus genotype in relation to tumor induction in mice. Characterization of wild type strains with widely differing tumor profiles. The American journal of pathology 127: 243–261. 2437801

22. Freund R, Calderone A, Dawe CJ, Benjamin TL (1991) Polyomavirus tumor induction in mice: effects of polymorphisms of VP1 and large T antigen. J Virol 65: 335–341. 1845894

23. Freund R, Dawe CJ, Benjamin TL (1988) The middle T proteins of high and low tumor strains of polyomavirus function equivalently in tumor induction. Virology 167: 657–659. 2849243

24. Freund R, Garcea RL, Sahli R, Benjamin TL (1991) A single-amino-acid substitution in polyomavirus VP1 correlates with plaque size and hemagglutination behavior. Journal of virology 65: 350–355. 1845896

25. Bauer PH, Cui C, Liu WR, Stehle T, Harrison SC, et al. (1999) Discrimination between sialic acid-containing receptors and pseudoreceptors regulates polyomavirus spread in the mouse. J Virol 73: 5826–5832. 10364334

26. Carroll J, Dey D, Kreisman L, Velupillai P, Dahl J, et al. (2007) Receptor-binding and oncogenic properties of polyoma viruses isolated from feral mice. PLoS Path 3: e179.

27. Gilbert J, Benjamin T (2004) Uptake pathway of polyomavirus via ganglioside GD1a. J Virol 78: 12259–12267. 15507613

28. O'Hara SD, Stehle T, Garcea R (2014) Glycan receptors of the Polyomaviridae: structure, function, and pathogenesis. Curr Opin Virol 7: 73–78. doi: 10.1016/j.coviro.2014.05.004 24983512

29. You J, O'Hara SD, Velupillai P, Castle S, Levery S, et al. (2015) Ganglioside and non-Ganglioside Receptors for the Mouse Polyomavirus. PLoS Path 11: e1005175.

30. Low JA, Magnuson B, Tsai B, Imperiale MJ (2006) Identification of gangliosides GD1b and GT1b as receptors for BK virus. J Virol 80: 1361–1366. 16415013

31. Krissinel E, Henrick K (2007) Inference of macromolecular assemblies from crystalline state. J Mol Biol 372: 774–797. 17681537

32. Blaum BS, Hannan JP, Herbert AP, Kavanagh D, Uhrin D, et al. (2015) Structural basis for sialic acid-mediated self-recognition by complement factor H. Nat Chem Biol 11: 77–82. doi: 10.1038/nchembio.1696 25402769

33. Chen Y, Tan M, Xia M, Hao N, Zhang XC, et al. (2011) Crystallography of a Lewis-binding norovirus, elucidation of strain-specificity to the polymorphic human histo-blood group antigens. PLoS Path 7: e1002152.

34. Neu U, Hengel H, Blaum BS, Schowalter RM, Macejak D, et al. (2012) Structures of Merkel cell polyomavirus VP1 complexes define a sialic acid binding site required for infection. PLoS Path 8: e1002738.

35. Merritt EA, Kuhn P, Sarfaty S, Erbe JL, Holmes RK, et al. (1998) The 1.25 A resolution refinement of the cholera toxin B-pentamer: evidence of peptide backbone strain at the receptor-binding site. J Mol Biol 282: 1043–1059. 9753553

36. Neu U, Woellner K, Gauglitz G, Stehle T (2008) Structural basis of GM1 ganglioside recognition by simian virus 40. Proc Natl Acad Sci U S A 105: 5219–5224. doi: 10.1073/pnas.0710301105 18353982

37. Cubrilovic D, Haap W, Barylyuk K, Ruf A, Badertscher M, et al. (2014) Determination of protein-ligand binding constants of a cooperatively regulated tetrameric enzyme using electrospray mass spectrometry. ACS Chem Biol 9: 218–226. doi: 10.1021/cb4007002 24128068

38. Fotinou C, Emsley P, Black I, Ando H, Ishida H, et al. (2001) The crystal structure of tetanus toxin Hc fragment complexed with a synthetic GT1b analogue suggests cross-linking between ganglioside receptors and the toxin. J Biol Chem 276: 32274–32281. 11418600

39. Sharmila DJ, Veluraja K (2006) Conformations of higher gangliosides and their binding with cholera toxin—investigation by molecular modeling, molecular mechanics, and molecular dynamics. J Biomol Struct Dyn 23: 641–656. 16615810

40. Zhang Y, Yamamoto S, Yamaguchi T, Kato K (2012) Application of paramagnetic NMR-validated molecular dynamics simulation to the analysis of a conformational ensemble of a branched oligosaccharide. Molecules 17: 6658–6671. doi: 10.3390/molecules17066658 22728360

41. Kuan CT, Chang J, Mansson JE, Li J, Pegram C, et al. (2010) Multiple phenotypic changes in mice after knockout of the B3gnt5 gene, encoding Lc3 synthase—a key enzyme in lacto-neolacto ganglioside synthesis. BMC Dev Biol 10: 114. doi: 10.1186/1471-213X-10-114 21087515

42. Parkkinen J, Rogers GN, Korhonen T, Dahr W, Finne J (1986) Identification of the O-linked sialyloligosaccharides of glycophorin A as the erythrocyte receptors for S-fimbriated Escherichia coli. Infect Immun 54: 37–42. 2875951

43. Schauer R (2000) Achievements and challenges of sialic acid research. Glycoconj J 17: 485–499. 11421344

44. Thomsson KA, Holmen-Larsson JM, Angstrom J, Johansson ME, Xia L, et al. (2012) Detailed O-glycomics of the Muc2 mucin from colon of wild-type, core 1- and core 3-transferase-deficient mice highlights differences compared with human MUC2. Glycobiology 22: 1128–1139. doi: 10.1093/glycob/cws083 22581805

45. Burmeister WP, Guilligay D, Cusack S, Wadell G, Arnberg N (2004) Crystal structure of species D adenovirus fiber knobs and their sialic acid binding sites. J Virol 78: 7727–7736. 15220447

46. Parrish CR (1990) Emergence, natural history, and variation of canine, mink, and feline parvoviruses. Adv Virus Res 38: 403–450. 2171302

47. Parrish CR (1991) Mapping specific functions in the capsid structure of canine parvovirus and feline panleukopenia virus using infectious plasmid clones. Virology 183: 195–205. 1647068

48. Neu U, Allen SA, Blaum BS, Liu Y, Frank M, et al. (2013) A structure-guided mutation in the major capsid protein retargets BK polyomavirus. PLoS Path 9: e1003688.

49. Neu U, Khan ZM, Schuch B, Palma AS, Liu Y, et al. (2013) Structures of B-lymphotropic polyomavirus VP1 in complex with oligosaccharide ligands. PLoS Path 9: e1003714.

50. Kanda T, Furuno A, Yoshiike K (1986) Mutation in the VP-1 gene is responsible for the extended host range of a monkey B-lymphotropic papovavirus mutant capable of growing in T-lymphoblastoid cells. J Virol 59: 531–534. 3488416

51. Xiong X, Martin SR, Haire LF, Wharton SA, Daniels RS, et al. (2013) Receptor binding by an H7N9 influenza virus from humans. Nature 499: 496–499. doi: 10.1038/nature12372 23787694

52. Xiong X, Coombs PJ, Martin SR, Liu J, Xiao H, et al. (2013) Receptor binding by a ferret-transmissible H5 avian influenza virus. Nature 497: 392–396. doi: 10.1038/nature12144 23615615

53. Ogawa-Goto K, Funamoto N, Ohta Y, Abe T, Nagashima K (1992) Myelin gangliosides of human peripheral nervous system: an enrichment of GM1 in the motor nerve myelin isolated from cauda equina. J Neurochem 59: 1844–1849. 1402926

54. Ishikawa Y, Gasa S, Minami R, Makita A (1987) Characterization of neutral glycosphingolipids from fetal human brain: evidence for stage-specific expression of the globo, ganglio, and neolacto series in the central nervous system. J Biochem 101: 1369–1375. 3667553

55. Kabsch W (2010) Integration, scaling, space-group assignment and post-refinement. Acta Crystallogr D Biol Crystallogr 66: 133–144. doi: 10.1107/S0907444909047374 20124693

56. Vagin A, Teplyakov A (2010) Molecular replacement with MOLREP. Acta Crystallogr D Biol Crystallogr 66: 22–25. doi: 10.1107/S0907444909042589 20057045

57. Adams PD, Afonine PV, Bunkoczi G, Chen VB, Davis IW, et al. (2010) PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr 66: 213–221. doi: 10.1107/S0907444909052925 20124702

58. Emsley P, Cowtan K (2004) Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 60: 2126–2132. 15572765

59. Murshudov GN, Vagin AA, Dodson EJ (1997) Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr 53: 240–255. 15299926

60. Painter J, Merritt EA (2006) Optimal description of a protein structure in terms of multiple groups undergoing TLS motion. Acta Crystallogr D Biol Crystallogr 62: 439–450. 16552146

61. The PyMOL Molecular Graphics System, Version 1.5.0.4, Schrödinger, LLC.

62. Howell PL, Smith GD (1992) Identification of heavy-atom derivatives by normal probability methods. J Appl Cryst 25: 81–86.

63. Kleywegt GJ, Jones TA (1996) xdlMAPMAN and xdlDATAMAN—Programs for reformatting, analysis and manipulation of biomacromolecular electron-density maps and reflection data sets. Acta Crystallogr D Biol Crystallogr 52: 826–828. 15299647

64. Kato Y, Kuan CT, Chang J, Kaneko MK, Ayriss J, et al. (2010) GMab-1, a high-affinity anti-3'-isoLM1/3',6'-isoLD1 IgG monoclonal antibody, raised in lacto-series ganglioside-defective knockout mice. Biochem Biophys Res Commun 391: 750–755. doi: 10.1016/j.bbrc.2009.11.132 19944071

65. Svennerholm L (1969). Comprehensive Biochemistry.

66. Taube S, Jiang M, Wobus CE (2010) Glycosphingolipids as receptors for non-enveloped viruses. Viruses 2: 1011–1049. doi: 10.3390/v2041011 21994669

67. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28: 27–30. 10592173

Štítky
Hygiena a epidemiologie Infekční lékařství Laboratoř

Článek vyšel v časopise

PLOS Pathogens


2015 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#