Structure of the Low pH Conformation of Chandipura Virus G Reveals Important Features in the Evolution of the Vesiculovirus Glycoprotein

Chandipura virus (CHAV, belonging to vesiculovirus genus), is an emerging pathogen associated with deadly encephalitis, principally among children, in India. CHAV single envelope glycoprotein G is involved in two successive steps of virus entry: receptor recognition that allows endocytosis of the virion and fusion of viral and cellular membranes that allows release of the viral nucleocapsid into the cytoplasm for the subsequent steps of infection. As for other rhabdovirus, fusion is triggered at low pH and catalyzed by a conformational change of G from its pre-fusion toward its post-fusion conformation. We have determined the crystalline structure of CHAV-G in its post-fusion conformation. An overall comparison of this structure with the previously reported post-fusion conformation of vesicular stomatitis virus G reveals the selective pressure of the humoral immune response: The major antigenic site, located in the most exposed domain of the pre-fusion conformation of the glycoprotein, is also the most structurally divergent region of G. The structure also highlights the remarkable plasticity of G in terms of local structures particularly in some hinge regions and reveals that the vesiculovirus can be classified in two groups based on the nature of their pH sensitive switches which trigger the conformational change of G.

Published in the journal: . PLoS Pathog 11(3): e32767. doi:10.1371/journal.ppat.1004756
Category: Research Article
doi: 10.1371/journal.ppat.1004756


Chandipura virus (CHAV, belonging to vesiculovirus genus), is an emerging pathogen associated with deadly encephalitis, principally among children, in India. CHAV single envelope glycoprotein G is involved in two successive steps of virus entry: receptor recognition that allows endocytosis of the virion and fusion of viral and cellular membranes that allows release of the viral nucleocapsid into the cytoplasm for the subsequent steps of infection. As for other rhabdovirus, fusion is triggered at low pH and catalyzed by a conformational change of G from its pre-fusion toward its post-fusion conformation. We have determined the crystalline structure of CHAV-G in its post-fusion conformation. An overall comparison of this structure with the previously reported post-fusion conformation of vesicular stomatitis virus G reveals the selective pressure of the humoral immune response: The major antigenic site, located in the most exposed domain of the pre-fusion conformation of the glycoprotein, is also the most structurally divergent region of G. The structure also highlights the remarkable plasticity of G in terms of local structures particularly in some hinge regions and reveals that the vesiculovirus can be classified in two groups based on the nature of their pH sensitive switches which trigger the conformational change of G.


Chandipura virus (CHAV; Family Rhabdoviridae, Genus Vesiculovirus) is an emerging human pathogen associated with deadly encephalitis, principally among children in the tropical areas of India. CHAV can be isolated from natural populations of Phlebotomine sandflies [13] and, in the recent years, has caused several outbreaks with high mortality rate [4]. Despite its importance in public health, there are no structural data available on CHAV proteins. In fact, most of the functional and structural information on CHAV has been inferred from studies performed on the prototype vesiculovirus, vesicular stomatitis virus serotype Indiana (VSVIND). This inference is reasonable due to the conservation of the amino acid (aa) sequences of the proteins of the two viruses [5]; for instance, aa sequences of glycoproteins G from both viruses share 40% identity and around 65% similarity [6].

The entry of rhabdoviruses into the host cell is mediated by glycoprotein G which constitutes the spikes that protrude at the viral surface [7]. First, G mediates viral attachment via interaction with a cellular receptor that, in the case of VSVIND, has been demonstrated to be the LDL receptor and its family members [8]. Then, after endocytosis [9,10], upon acidification of the inner environment of the endosome, G undergoes a low-pH-induced structural transition from its trimeric pre-fusion form (PRE) towards its trimeric post-fusion conformation (POST) that drives fusion between the viral and endosomal membranes [11,12]. During this conformational change, internal hydrophobic motifs (the so-called fusion loops) interact with one or both participating membranes [13], resulting in their destabilization and merger. Remarkably, low-pH-induced structural transition of rhabdoviral G is reversible [11,14]. In fact, there is a pH dependent thermodynamic equilibrium between different states of G. At neutral pH, G exists as a population of mostly monomeric species with a minority of PRE trimers [12,15] whereas, at low pH, the G distribution is shifted toward the trimeric POST state [15,16]. This is the main difference between rhabdoviral G and other viral fusion glycoproteins activated at low pH for which the PRE conformation is metastable [17,18].

Mature glycoprotein G is about 500 aa long (495 for VSVIND and 509 for CHAV). G has a huge ectodomain (452 aa for CHAV-G) anchored in the lipid bilayer by a single transmembrane domain (TMD) which is located upstream a short C-terminal intra-viral segment. The crystal structures of both PRE and POST states of VSVIND-G ectodomain (VSVIND-Gth, obtained by limited proteolysis of G at the surface of virions of the Indiana Mudd-Summers strain with thermolysin) have been previously reported [19,20]. This revealed that the POST structure of VSV-G has the same fold as the presumptive post-fusion state of glycoprotein gB of herpesviruses [21] and baculovirus gp64 [22]. This defined a new class (class III) of fusion proteins [23].

G ectodomain folds into three distinct structural domains [24] (S1 Fig.). A central domain (CD) is made of a lateral β-sheet rich region and a central helix which is involved in trimerization of the molecule in both PRE and POST states. A pleckstrin homology domain (PHD) is inserted in a loop of CD and the fusion domain (FD) is itself inserted in a loop of PHD. The FD is made of an extended β-sheet structure at the tip of which are located the two hydrophobic fusion loops. These rigid domains are connected by segments (R1–4) which, together with the C-terminal segment (R5) connecting CD to TMD, undergo major refolding events during the structural change.

During the conformational change from PRE to POST conformation, both FD and the C-terminal part of the ectodomain rotate around CD [20,24] (S1 Fig.). As a consequence, both the fusion loops and the TM domain move from one end of the molecule to the other. The FD moves relatively to PHD due to refolding of segments R2 and R3. There is also a major refolding event in the two segments connecting CD to PHD (R1 and R4) leading to lengthening of the CD central helix F. Similarly, the C-terminal segment of the ectodomain both completely relocates and refolds into α-helix H, which positions itself, in an antiparallel manner, in the grooves of the trimeric core made by the three helices F. This results in the formation of a six-helix bundle that buries acidic amino acids, which constitute pH sensitive molecular switches [19,25]. The complex conformational change has been demonstrated to involve monomeric intermediates [15].

In this work, we have characterized the fusion and biochemical properties of CHAV-G and compared them to those of VSV-G. We also determined the structure of a soluble form of CHAV-G ectodomain (CHAV-Gth, generated by thermolysin limited-proteolysis of recombinant VSV particles in which the G gene was replaced by that of CHAV[26]) in the POST conformation at 3.6 Å resolution. An overall comparison of this structure with VSV-Gth POST, shows a high structural similarity as expected from the comparison of aa sequences. Among the three domains of G, the PHDs appear to be the most divergent and the largest differences are confined to the secondary structure of the major antigenic site of rhabdoviruses. Finally, local differences indicate that CHAV, together with a related group of vesiculoviruses, has evolved alternate structural solutions in hinge regions between PHD and FD and distinct pH sensitive switches.

Results and Discussion

Fusion properties of CHAV glycoprotein

We compared the fusion activity of both CHAV-G and VSV-G in a cell-cell fusion assay. For this, BSR cells were co-transfected with a plasmid encoding G (either VSV-G or CHAV-G) and rabies virus (RABV) phosphoprotein fused to the green fluorescent protein (P-GFP). P-GFP is unable to passively diffuse into the nucleus [27]. This allows an easy visualization and counting of multinucleated cells (syncytia) formed as a result of the fusion activity of the glycoprotein [25,28]. Transfected cells were incubated for 10 minutes with medium adjusted to pH values from 5 to 7, which was then replaced by medium buffered at pH 7.4. The cells were then further incubated for 1h at 37°C before fixation.

We could observe formation of massive syncytia between pH 5.0 and 6.0 for cells expressing either CHAV-G or VSVIND-G. At pH 6.3, syncytia were smaller for both glycoproteins. At pH 6.5, no syncytia were detected for CHAV-G whereas small syncytia (less than 20 nuclei) were still observed for VSVIND-G up to pH 7.0 (Fig. 1).

Fusion activity of CHAV-G and VSV-G analyzed in a cell-cell fusion assay.
Fig. 1. Fusion activity of CHAV-G and VSV-G analyzed in a cell-cell fusion assay.
A) BSR cells were transfected with plasmids expressing the full-length G proteins and RABV P-GFP to facilitate the observation of syncytia. At 24 h post-transfection, the cells were incubated for 10 min with cell medium adjusted to the indicated pH, and then replaced by fresh medium at pH 7.4. The cells were kept at 37°C for 1 h before fixation. Nuclei were stained with DAPI. All images are representative examples from at least three independent experiments. The scale is the same for all micrographs. B) Percentage of fluorescent cells having formed a syncytia (containing more than 3 nuclei) as a function of pH for both VSVIND-G and CHAV-G. C) Distribution of the size of syncytia as a function of pH for both VSVIND-G and CHAV-G (at each pH, the examined field corresponded to about 400 transfected fluorescent cells at pH7).

Purification and biochemical properties of CHAV-Gth

We purified CHAV-G ectodomain from a recombinant VSV in which the G gene was replaced for that of CHAV (VSV/CHAV-G) [29]. This chimeric virus is non-pathogenic and can be produced at high titers in BSR cells. At pH 6.0, thermolysin produces a single cut in CHAV-G between residues 419 and 420 (the cleavage site was determined by mass spectrometry). After proteolysis, the CHAV-G ectodomain (CHAV-Gth, residues 1–419) was purified by successive steps of ultracentrifugation and chromatography as previously described [26].

In the case of VSV-G, the pre-fusion trimer, although detected at the viral surface [12], is not stable in solution and only monomers of G are detected at high pH [11,15]. On the other hand, at low pH, the post-fusion trimer is stable [11,15]. Monomers and post-fusion trimers of G have distinct sedimentation coefficient and can be easily separated in a sucrose gradient [11,25]. Therefore, we analyzed the quaternary structure of CHAV-Gth on linear 5–20% (w/v) sucrose gradients at different pHs and compared it with that of VSVIND-Gth.

The overall behavior of CHAV-Gth in the sucrose gradients was similar to that of VSVIND-Gth (Fig. 2A and 2B). At pH values from 7.0 to 8.0, Gth was found in fractions 9–12 corresponding to monomeric species, while at pH 5.5 and 6.0 only the trimeric post-fusion form was found in lower fractions of the gradient (fractions 4–8). However, the pK of the transition (defined as the pH value at which monomeric and trimeric G are found in equal amounts) was slightly higher for CHAV-Gth (pH 6.7) than for VSV-Gth (pH 6.5).

Biochemical properties of CHAV-G<sub>th</sub> and VSV-G<sub>th</sub>.
Fig. 2. Biochemical properties of CHAV-Gth and VSV-Gth.
Oligomerization assay of CHAV-Gth (A) and VSV-Gth (B). Proteins were preincubated in 50 mM TRIS or MES buffers adjusted to the indicated pH. The solution was loaded on a linear 5–20% sucrose gradient at the same pH and spun at 35000 rpm for 16 h and 4°C. After centrifugation, 12 fractions were collected and analyzed by 10% SDS-PAGE gels stained with Coomassie blue. The dashed line indicates the gradient boundary between monomeric pre- and trimeric post-fusion forms. C) Membrane flotation assay of CHAV-Gth incubated with liposomes over a range of pH values from 5.7 to 8.5. For each experiment, top (t) and bottom (b) fractions were collected and analyzed on 10% SDS PAGE stained with Coomassie Blue. The right panel shows the experiments performed at pH 8.5 and the reversibility of the low-pH-induced association of Gth with liposomes after incubation at pH 5.7 followed by a re-incubation at pH 8.5.

We also investigated the membrane interaction properties of CHAV-Gth by liposome flotation assay. CHAV-Gth/liposome interaction was assessed over a broad pH-range, from pH 5.7 to pH 8.5 (Fig. 2C). From pH 5.7 to 6.6, all the glycoprotein was located in the upper layer indicating a massive association of CHAV-Gth with liposomes. At pH 7.0, the quantity of CHAV-Gth recovered in the bottom fraction increased considerably and only half of the glycoprotein was still detected in association with membranes. From pH 7.3 to 8.5, CHAV-Gth was exclusively found in the bottom fraction indicating a poor membrane interacting activity. Finally, when CHAV-Gth-liposome complexes formed at pH 5.7 were brought to pH 8.5, membrane flotation experiments revealed a complete dissociation of glycoproteins from liposomes.

Taken together, our results indicate that the fusion properties and the associated low-pH-induced conformational change are very similar for CHAV-G and VSV-G. Particularly, as VSV-G [11,15], CHAV-G pre-fusion trimer is not stable in solution and low-pH-induced membrane association of CHAV-G is reversible.

Structure of the post-fusion conformation of CHAV-Gth

We screened several crystallization conditions for CHAV-Gth in presence of n-dodecyl β-maltoside (DDM) to avoid rosette formation via aggregation of the protein through its fusion loops. In these experiments, we could obtain CHAV-Gth crystals at pH 4.6 and determine the structure of the protein at 3.6 Å resolution by molecular replacement. Statistics of data collection and refinement are given in Table 1. The comparatively low resolution of the structure and starting phases obtained from molecular replacement with a 40% identity model (only 25% between PH domains) did not preclude an accurate final model in this case. This is due to a large solvent content and three molecules in the asymmetric unit. This yields over 20,000 unique reflections with which to refine three times 3,250 atoms per protomer with tight non-crystallographic symmetry and stereochemical restraints and conservative temperature factor models (see Material and Methods section for more details). The structure's accuracy is attested by excellent final statistics, particularly agreement between model and data and stereochemical parameters including the most unbiased indicators (Table 1, R-free and Ramachandran plot, respectively). Final temperature factors are quite reasonable for this resolution except in the distal parts of fusion domains, in accordance with these parts being the most flexible tip of an elongated molecule (S2 Fig.). The parts omitted in the initial molecular replacement model showed density and could be built de novo after the first few rounds of refinement and rebuilding, and the final composite simulated annealing omit map confirms unambiguous placing of main chain and most side chains except in the distal parts of FD (S3 Fig.). Accordingly, in the following sections, we discuss only those parts of the model, particularly the hydrogen bond networks, for which density is unambiguous.

Tab. 1. Data collection and refinement statistics.
Data collection and refinement statistics.
Statistics for the highest-resolution shell are shown in parentheses.

The single trimer in the asymmetric unit has the three protomers in a classic post-fusion conformation. The segments R1 to R5 connecting the different domains were traced in the experimental electron density map. The three chains were built from residues 1 to 415, 1 to 416 and 1 to 413, for chains A, B and C, respectively; with internal breaks in segments 26–33 and 27–33 of chain A and C. Those breaks reflect the intrinsic flexibility of segment R1.

The overall structure of CHAV-Gth is very similar to that of VSVIND-Gth in its post-fusion form (Fig. 3A and B) [19]. The structural comparison between CHAV-G and VSVIND-G reveals that the three domains are more conserved than would be expected from the sequence divergence [30,31] (Table 2, Fig. 3C). CDs and the upper part of FDs (excluding the poorly defined extremity of FD) main chains are almost superimposable within experimental error (Table 2, S4 Fig.).

Overall structures of CHAV-G<sub>th</sub> and VSV<sub>IND</sub>-G<sub>th</sub> post-fusion forms.
Fig. 3. Overall structures of CHAV-Gth and VSVIND-Gth post-fusion forms.
(A) Structure of CHAV-Gth and VSVIND-Gth post-fusion trimers and (B) their respective protomers. The proteins are colored by domains and segments as defined in Table 2. The fusion loops at the tip of FD are in pink. The arrows indicate the most noticeable differences between both proteins (located in PHDs and segments R3). The black double headed arrows indicate the spacing between the tips of FDs (between Cα of residue 73 and Cα of residue 117 of the neighboring protomer). The corresponding distance is 15.7Å for VSV- G and 26.2 Å for CHAV-G. (C) Sequence alignment of CHAV-G and VSVIND-G. Conserved residues are in blue boxes whereas similar ones are in red (overall aa sequence conservation is 41%). The black dashes show gaps. Elements of secondary structure are indicated above the sequence for the two structures. They are colored as the domain they belong to, with the color code of Table 1. VSVIND-Gth helices and sheets are named as in [20]. Segments R1 to R5 that refold during the transition are also indicated with the color code of Table 1. Brown triangles indicate the glycosylation sites. Residues located in the fusion loops are labeled by pink stars. The asterisks indicate the last residue of the molecule for which electron density was still interpretable.

Tab. 2. Domains and segments nomenclature used in the text.
Domains and segments nomenclature used in the text.
Root mean square deviation (RMSD) and number of superimposable alpha carbons (Cα) are reported between indicated domains. As the distal (fusion loop-containing) part of FD is the most flexible part in all structures, we defined a “FD upper part” which was used to calculate the RMSD between VSV and CHAV-Gth. We also defined a PHD core corresponding to the 79 residues (out of 90) in CHAV-G which are in structurally equivalent positions in VSV-G. The RMSD between the two (PHD core + FD upper part) is not significantly different from the value obtained for the PHD core alone. This indicates that the relative orientation of the two domains is the same in both proteins.

PHD is the most divergent domain of vesiculovirus glycoproteins

The PHDs are more divergent than FDs and CDs. Only 79 residues of the PHD (called PHD core) out of 90 in CHAV-G are in structurally equivalent positions with VSV-G (Table 2). The largest differences are confined to helices D and E which are missing in CHAV-G (see the orange arrow on Fig. 3B). A comparison of the sequence of both proteins at these positions reveals that CHAV-G bears an insertion of 4 aa residues at the level of helix E resulting instead in the formation of a small two stranded β-sheet (Fig. 4). Also, in the equivalent position to helix D, CHAV-G has a deletion of two residues leading to a loop conformation (Fig. 4).

Ribbon representation of CHAV-G and VSV-G PHDs.
Fig. 4. Ribbon representation of CHAV-G and VSV-G PHDs.
The segments having different structures between the two domains and corresponding to Rhabdoviruses major antigenic site are highlighted in green and purple. Helices and β-sheets are labeled in capital and lower case letters, respectively, according to the previous VSV-G nomenclature. The sequence alignment below the structures corresponds to the two segments highlighted in the ribbon representation. The position of the first and last aa residues of the two sequences is indicated on the structures. Conserved residues between CHAV-G and VSV-G are highlighted in blue.

VSV-G segments containing both helix D and helix E are the most exposed part of the glycoprotein in its PRE conformation at the viral surface and, in fact, this region constitutes the major antigenic site of several rhabdovirus glycoproteins [3235]. This suggests that the structural divergence between CHAV-G and VSVIND-G in this region is the result of the selective pressure from the humoral response of the host immune system.

Distinct structural solutions in the hinge region connecting PHD and FD

A second distinct feature between CHAV-G and VSVIND-G lies in the hinge R2 and R3 connecting both PHD and FD (green arrow on Fig. 3B). Although the relative orientation of the two domains is the same in both proteins (Table 2), there are differences in the secondary structure of these hinges.

In the POST conformation of VSVIND-G, R2 is mainly a large loop and extends until the first residues of the β-strand c of FD, whereas R3 forms helix C. The contact between the two segments essentially involves hydrophobic residues (Fig. 5).

Ribbon representation of CHAV-G and VSV<sub>IND</sub>-G R2-R3 hinge regions.
Fig. 5. Ribbon representation of CHAV-G and VSVIND-G R2-R3 hinge regions.
Protein domains are colored according to Table 2. The main elements of the secondary structure of VSVIND-G and their equivalents in CHAV-G are indicated with capital (for α-helices) and lower case (for β-strands) letters. Residues located at the interface between R2, R3 are in sticks representations. Note their polar character in the case of CHAV-G. The sequence alignment of segments R2 and R3 is displayed below the structures. The position of the first and last aa residues of the two sequences is indicated on the structures. Conserved residues between CHAV-G and VSV-G are highlighted in blue.

For CHAV-G, the amino-terminal part of R2 adopts a β-strand conformation, which is followed by a small loop and a small helix; CHAV-G R3 has a structure which is similar to VSV-G R2, a long loop ending in a β-strand. The interaction between the two segments is mediated by charged and polar residues (Fig. 5).

Despite those differences in R2 and R3 regions, the relative positions of FD and PHD are identical for both glycoproteins (Table 2). This suggests that the hinge regions need only to be flexible to allow the rotation of the FD.

Alignment of G aa sequences for 12 vesiculoviruses (S5 Fig.) reveals that the sequences of R2 and R3 segments are rather divergent among vesiculoviruses. This suggests that the structural constraints (i.e. flexibility and relative position of PHD and FD in the POST conformation) can be accommodated by very diverse amino acid sequences.

Different pH sensitive molecular switches

Four acidic residues of VSVIND-G (D268, D274, D393 and D395), which are brought close together in the POST six-helix bundle (Fig. 6A and 6B), have been demonstrated to play the role of pH-sensitive switches for the transition back from the POST toward the PRE conformation of VSV-G [25]. The acidic character of the residues corresponding to D274, D393 and D395 of VSV-G is conserved among the vesiculovirus genus (S5 Fig.). However, this is not the case for the residue corresponding to D268 which is a major pH-sensitive switch of VSVIND-G [25].

pH sensitive molecular switches in the post-fusion structures of VSV-G<sub>th</sub> and CHAV-G<sub>th</sub>.
Fig. 6. pH sensitive molecular switches in the post-fusion structures of VSV-Gth and CHAV-Gth.
pH sensitive molecular switches located in the 6-helix bundle of the post-fusion form of VSVIND-G (A-B) and CHAV-G (D-E). A and D are lateral views of two protomers interacting by the central helices F. B and E are top views of the trimeric interfaces for both proteins. The acidic character of VSV-G residues D274, D393 and D395 is conserved in CHAV-G whereas D268, a major pH sensor in VSV-G [25], is replaced by an alanine. Instead, in CHAV-G, another pH sensitive switch is formed by the triad H209, E234 and D269. C: Sequence alignments of 12 vesiculoviruses in regions flanking the molecular switches D268 (single arrow) and D269/D234 (double headed arrow) from VSVIND-G and CHAV-G, respectively. Vesicular Stomatitis Indiana Virus (IND, GenBank: AAA48370.1), Maraba Virus (MAR, GenBank: AEI52254.1), Cocal Virus (COC, GenBank: ACB47437.1), Vesicular Stomatitis Alagoas Virus (ALA, GenBank: ACB47442.1), Carajas Virus (CAR, Patent: JP 2010503660-A 2 04-FEB-2010), Vesicular Stomatitis New Jersey Virus (NJE, NCBI: YP_009047084.1), Chandipura Virus (CHA, NCBI: YP_007641380.1), Piry Virus (PIR, Swiss-Prot: Q85213.1), Isfahan Virus (ISF, NCBI: YP_007641385.1), Jurona Virus (JUR, GenBank: AEG25348.1), Perinet Virus (PER, GenBank: AEG25354.1), Malpais Spring Virus (MAL, GenBank: AGI04017.1). Conserved residues through the whole genus are in shade of blue, through the CHAV group are in shade of green and through the IND groups are in shades of brown.

Specifically, in CHAV-G, the residue at this position is an alanine (A271). The absence of this pH-sensitive switch seems to be compensated by CHAV-G’s residues D269 and E234 which form an intra-protomer hydrogen bond with H209 in their vicinity (Fig. 6D and 6E). This clustering is absent for the corresponding, and non-protonable, residues of VSVIND-G POST. Clustering of protonable residues is known to affect the pKa of their side chain. Indeed, when the PropKa program [36] was used to compute the pKa of the side chains of the cluster based on CHAV-G POST structure, it predicted a pKa of about 10.5 for the carboxylate of residue D269. The only other residue with such a pKa shift in CHAV-G POST is residue E399 which corresponds to VSVIND-G D395 which has also been demonstrated to contribute to the pH sensitivity of VSVIND-G POST. Similarly, VSVIND-G D268 was calculated to have a pKa of about 14.5 in agreement with its established role as a major pH sensitive switch in VSVIND-G POST. This further suggests the likely role of the CHAV-G POST cluster D269, E234 and H209 as a pH sensitive switch. The deprotonation of these residues at high pH is certainly destabilizing the POST conformation of CHAV-G.

Remarkably, alignment of vesiculoviruses G aa sequences (Fig. 6) reveals that the vesiculoviruses can be classified in two groups. Either they have an aspartate residue in the position corresponding to VSV-G D268, or they have a histidine in position corresponding to CHAV-G H209 and an acidic residue in position corresponding to CHAV-G E234 and D269 (Figs. 6C and S5). Unsurprisingly, this segregation which relies only on vesiculovirus glycoprotein local sequences is in agreement with global phylogenetic studies on the genus [37,38].

This different switch of CHAV results in a tighter interaction between the base of the six-helix bundle and the PH domain (Fig. 7). As a consequence, the PHD has not exactly the same orientation relative to CD. This is why the FDs are slightly spread apart in CHAV-G (Fig. 3B).

Packing of PHD domain against the central domain in CHAV-G<sub>th</sub> and VSV<sub>IND</sub>-G<sub>th</sub>.
Fig. 7. Packing of PHD domain against the central domain in CHAV-Gth and VSVIND-Gth.
In A and B are represented the top views of the protomers. The upper black lines are crossing the structures at the level of the Cα of two pairs of structurally equivalent residues in both glycoproteins which are S285 (located in helix F) and T10 for CHAV-G, and T282 and Q10 for VSV-G. The lower black lines show the border between PHD (orange) and CD (red). In C is shown the interactions at the base of helices F and H of CHAV-Gth involved in the close packing of PHD relative to CD. Residues implicated in the interaction and their immediate neighbors are in stick representation. In D are represented the equivalent residues of VSVIND-G.

Final conclusions

The comparison of VSVIND-G and CHAV-G structures reveal that the domains do not evolve at the same rate and that the major actor on G evolution is the immune system. As a consequence, the most exposed domain in the PRE conformation (i.e. the PHD), which is the target of neutralizing antibodies, is also by far the most divergent domain among vesiculoviruses. This domain, located at the top of the PRE conformation, is also probably involved in receptor recognition. Therefore, it is probable that the humoral immune response is at the origin of receptor diversification among the rhabdovirus family.

The principal finding of this work is that, vesiculovirus G have diverged quite far in some respects but retained important functional and structural features. These features can be thus pointed at as essential for the protein’s function. As an example, the relative orientation of CD and PHD is different in CHAV-G and VSV-G POST protomers, but the relative orientation of PHD and FD remains the same despite the sequence and local structure divergence of segments R2 and R3 and the largely different sets of interactions between FD. This highlights the mechanistic importance of maintaining PHD-FD orientation. This also shows that the most important feature for segments R2 and R3 is their flexibility, but not their secondary structure. As a consequence, the structural constraints which exert on these segments are weak and R2 and R3 can tolerate quite different aa sequences. As another example, although pH sensitive switches are absolutely required for the correct function of G, those switches can be located in alternate places, as in the CHAV-G D269-E234-H209 cluster, the counterpart of which does not exist in VSV-G (Figs. 6C and S5).

In conclusion, the vesiculovirus glycoprotein, despite experiencing strong functional constraints, exhibits a remarkable plasticity in terms of sequence and even local structures. This is reminiscent of the remarkable structural conservation of the class II architecture despite the lack of any sequence similarity [39,40]. This suggests that evolvability is a key feature for viral fusion glycoproteins: They display a 3D structure that is robust in the face of amino acid sequence variation because of being malleable in local conformations, yet still maintaining the overall trimeric post-fusion conformation that is essential for their function.

Materials and Methods

Production and purification of CHAV-Gth

CHAV-Gth was purified as previously described [26]. Briefly, a recombinant VSV harboring the CHAV-G gene [29] was produced at high titers in BSR cells, clones of BHK-21 (baby hamster kidney, ATCC CCL-10). Purified virions were submitted to proteolysis with thermolysin at pH 6.0 for 1.5h. This treatment generates a soluble fragment of G comprising the residues 1–419 (Gth) as assessed by mass spectrometry. After digestion, the virus was pelleted and the viral pellet was suspended in 50 mM HEPES pH 7.5 containing 2% glycerol. The suspension containing Gth was clarified by ultracentrifugation at 25000 rpm (SW35 rotor) for 1h on a 20% w/v sucrose cushion. Gth purification was achieved on a DEAE trisacryl column (GE healthcare), followed by a size exclusion step using a Superdex S200 10/30 column (GE healthcare). Fractions containing Gth were concentrated up to 10 mg/ml and stored at −80°C until use.

Cell-cell fusion assay

BSR cells plated on glass coverslips at 70% confluence were co-transfected with pCAGGS plasmids encoding VSV-G or CHAV-G, and P-GFP plasmid encoding the phosphoprotein of rabies virus fused to GFP [27]. Twenty four hours after transfection, cells were incubated with fusion buffer (DMEM-10mM MES or HEPES) at various pHs (from 5.0 to 7.0) for 10 minutes at 37°C. Cells were then washed once and incubated with DMEM-10mM HEPES-NaOH buffered at pH 7.4, 0.1% BSA at 37°C for 1 hour. Cells were fixed with 4% paraformaldehyde in 1× PBS for 15 min. Cells nuclei were stained with DAPI and syncytia formation was analyzed with a Zeiss Axiovert 200 fluorescence microscope at 20x magnification.

Oligomerization assay

The oligomeric state of the Gth protein was determined by sucrose gradient centrifugation. 100 μg of VSVIND-Gth or CHAV-Gth were incubated in 200 μl of buffer containing either 50 mM MES-NaOH (for pH 5.5 to 6.7) or 50 mM Tris-HCl (for pH 7 to 8.5) plus 100 mM NaCl at 37°C for 30 minutes. The protein solution was overlaid on a 5-to-20% w/v continuous sucrose gradient in 100 mM NaCl, 50 mM Tris-HCl or 50mM MES-NaOH at the required pHs. Following centrifugation at 35,000 rpm for 16 hours in an SW41 rotor, 12 equal fractions were collected manually from the bottom of the gradient. When required, the fractions were dissolved and reconcentrated to the original volume using ultrafiltration devices (30 kDa molecular weight cut off) to get rid of the excess of sucrose. 20 μl of each fraction was analyzed on 10% SDS-PAGE gels stained with Coomassie blue.

Membrane flotation assay

Membrane flotation assays were performed as previously described [15]. Briefly, 500 μg of fresh liposomes (phosphatidylcholine: phosphatidylethanolamine: gangliosides, 7:3:1) were incubated with 100 μg of purified Gth at the desired pH during 30 minutes at 37°C (final volume 600 μL). The lipid-protein mix was mixed with a 80% w/v sucrose solution so that final concentration of sucrose was around 65%. The mixture was overlaid with 2 mL of 50% w/v sucrose solution and 1.5 mL of 10% w/v sucrose solution, both solutions at the desired pH. Cushions were spun for 16 hr at 20°C in a SW55 Beckman rotor at 45,000 r.p.m. 20 μL of the top and bottom fraction were analyzed by SDS-PAGE. Gth in the top fraction was considered as liposome-bond, Gth in the bottom fraction was considered as not interacting with the membranes. For the reversibility experiment, lipid-protein mix was re-incubated at pH 8.5 after a first stage at pH 5.7.

Crystallization of CHAV-Gth and structure determination

Crystals were obtained at 293 K by the hanging drop vapour diffusion method. Drops were prepared by mixing 1 μl of CHAV-Gth (4 mg/ml) supplemented with 0.2% DDM with 1 μl of the reservoir solution (12% PEG 3350, 0.1 M sodium acetate pH 4.6) and equilibrated against 500 μl of reservoir solution.

Crystals were harvested in a cryoprotecting solution containing reservoir solution supplemented with 30% glycerol and flash-cooled by plunging in liquid nitrogen. Diffraction data were collected at 100 K at PROXIMA 1 and ID-29 beam lines in SOLEIL and ESRF synchrotrons (France). Data for the best crystal were integrated and scaled using XDS/XSCALE [41]. The crystals belong to space group C2 with cell parameters consistent with 3 molecules in the asymmetric unit and 60% solvent content [42].

The structure was solved by molecular replacement (MR) with the PHENIX package [43]. The VSV-G post-fusion protomer (PDB 2CMZ, chain A) was first trimmed of non-conserved sidechains using the default procedure in phenix.sculptor. The three search models were the CD, FD and PHD further trimmed of carbohydrates and of a couple of residues upstream and downstream of refolding segments R1 to R5. The nine expected domains (corresponding to 3 protomers) in the asymmetric unit were placed by phenix.phaser. These initial phases were improved by solvent flattening and NCS averaging. Initial maps were clear enough to rebuild CDs and remove those parts of PHDs that did not fit the electron density maps. The model was iteratively rebuiltin COOT [44] and refined with phenix.refine. The refinement protocol was tailored to reduce the parameter to data ratio as much as possible: Tight non-crystallographic symmetry restraints were enforced between the three protomers in the asymmetric unit. B-factors were parameterized initially by TLS domains, allowing subsequent individual B-factor refinement with the very tight restraints on neighboring atoms implemented in recent versions of phenix.refine. This refinement protocol may allow very high B-factors in distal parts of elongated molecules due to the TLS component. The CDs, PHDs and the upper parts of FDs were carefully rebuilt and segments R1 to R5 added in as density became clear enough to model them. No Ramachandran restraints were applied in refinement at any stage to allow orthogonal validation of stereochemistry. All the structure figures were prepared using the program PyMol (PyMOL Molecular Graphics System. DeLano Scientific LLC, San Carlos, CA, USA.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5


1. Dhanda V, Rodrigues FM, Ghosh SN (1970) Isolation of Chandipura virus from sandflies in Aurangabad. Indian J Med Res 58: 179–180. 5528233

2. Mavale MS, Fulmali PV, Geevarghese G, Arankalle VA, Ghodke YS, et al. (2006) Venereal transmission of Chandipura virus by Phlebotomus papatasi (Scopoli). Am J Trop Med Hyg 75: 1151–1152. 17172384

3. Sudeep AB, Bondre VP, Gurav YK, Gokhale MD, Sapkal GN, et al. (2014) Isolation of Chandipura virus (Vesiculovirus: Rhabdoviridae) from Sergentomyia species of sandflies from Nagpur, Maharashtra, India. Indian J Med Res 139: 769–772. 25027088

4. Menghani S, Chikhale R, Raval A, Wadibhasme P, Khedekar P (2012) Chandipura Virus: an emerging tropical pathogen. Acta Trop 124: 1–14. doi: 10.1016/j.actatropica.2012.06.001 22721825

5. Cherian SS, Gunjikar RS, Banerjee A, Kumar S, Arankalle VA (2012) Whole genomes of Chandipura virus isolates and comparative analysis with other rhabdoviruses. PLoS One 7: e30315. doi: 10.1371/journal.pone.0030315 22272333

6. Masters PS, Bhella RS, Butcher M, Patel B, Ghosh HP, et al. (1989) Structure and expression of the glycoprotein gene of Chandipura virus. Virology 171: 285–290. 2741347

7. Albertini AA, Baquero E, Ferlin A, Gaudin Y (2012) Molecular and cellular aspects of rhabdovirus entry. Viruses 4: 117–139. doi: 10.3390/v4010117 22355455

8. Finkelshtein D, Werman A, Novick D, Barak S, Rubinstein M (2013) LDL receptor and its family members serve as the cellular receptors for vesicular stomatitis virus. Proc Natl Acad Sci U S A 110: 7306–7311. doi: 10.1073/pnas.1214441110 23589850

9. Cureton DK, Massol RH, Saffarian S, Kirchhausen TL, Whelan SP (2009) Vesicular stomatitis virus enters cells through vesicles incompletely coated with clathrin that depend upon actin for internalization. PLoS Pathog 5: e1000394. doi: 10.1371/journal.ppat.1000394 19390604

10. Johannsdottir HK, Mancini R, Kartenbeck J, Amato L, Helenius A (2009) Host cell factors and functions involved in vesicular stomatitis virus entry. J Virol 83: 440–453. doi: 10.1128/JVI.01864-08 18971266

11. Doms RW, Keller DS, Helenius A, Balch WE (1987) Role for adenosine triphosphate in regulating the assembly and transport of vesicular stomatitis virus G protein trimers. J Cell Biol 105: 1957–1969. 2824524

12. Libersou S, Albertini AA, Ouldali M, Maury V, Maheu C, et al. (2010) Distinct structural rearrangements of the VSV glycoprotein drive membrane fusion. J Cell Biol 191: 199–210. doi: 10.1083/jcb.201006116 20921141

13. Durrer P, Gaudin Y, Ruigrok RW, Graf R, Brunner J (1995) Photolabeling identifies a putative fusion domain in the envelope glycoprotein of rabies and vesicular stomatitis viruses. J Biol Chem 270: 17575–17581. 7615563

14. Gaudin Y, Tuffereau C, Segretain D, Knossow M, Flamand A (1991) Reversible conformational changes and fusion activity of rabies virus glycoprotein. J Virol 65: 4853–4859. 1870204

15. Albertini AA, Merigoux C, Libersou S, Madiona K, Bressanelli S, et al. (2012) Characterization of Monomeric Intermediates during VSV Glycoprotein Structural Transition. PLoS Pathog 8: e1002556. doi: 10.1371/journal.ppat.1002556 22383886

16. Roche S, Gaudin Y (2002) Characterization of the equilibrium between the native and fusion-inactive conformation of rabies virus glycoprotein indicates that the fusion complex is made of several trimers. Virology 297: 128–135. 12083843

17. Carr CM, Chaudhry C, Kim PS (1997) Influenza hemagglutinin is spring-loaded by a metastable native conformation. Proc Natl Acad Sci U S A 94: 14306–14313. 9405608

18. Gaudin Y (2000) Reversibility in fusion protein conformational changes. The intriguing case of rhabdovirus-induced membrane fusion. Subcell Biochem 34: 379–408. 10808339

19. Roche S, Bressanelli S, Rey FA, Gaudin Y (2006) Crystal structure of the low-pH form of the vesicular stomatitis virus glycoprotein G. Science 313: 187–191. 16840692

20. Roche S, Rey FA, Gaudin Y, Bressanelli S (2007) Structure of the prefusion form of the vesicular stomatitis virus glycoprotein g. Science 315: 843–848. 17289996

21. Heldwein EE, Lou H, Bender FC, Cohen GH, Eisenberg RJ, et al. (2006) Crystal structure of glycoprotein B from herpes simplex virus 1. Science 313: 217–220. 16840698

22. Kadlec J, Loureiro S, Abrescia NG, Stuart DI, Jones IM (2008) The postfusion structure of baculovirus gp64 supports a unified view of viral fusion machines. Nat Struct Mol Biol 15: 1024–1030. doi: 10.1038/nsmb.1484 18776902

23. Backovic M, Jardetzky TS (2009) Class III viral membrane fusion proteins. Curr Opin Struct Biol 19: 189–196. doi: 10.1016/ 19356922

24. Baquero E, Albertini AA, Vachette P, Lepault J, Bressanelli S, et al. (2013) Intermediate conformations during viral fusion glycoprotein structural transition. Curr Opin Virol 3: 143–150. doi: 10.1016/j.coviro.2013.03.006 23562213

25. Ferlin A, Raux H, Baquero E, Lepault J, Gaudin Y (2014) Characterization of pH sensitive molecular switches triggering the structural transition of VSV Glycoprotein from the post- toward the pre-fusion state. J Virol.

26. Baquero E, Buonocore L, Rose JK, Bressanelli S, Gaudin Y, et al. (2012) Crystallization and preliminary X-ray analysis of Chandipura virus glycoprotein G. Acta Crystallogr Sect F Struct Biol Cryst Commun 68: 1094–1097. doi: 10.1107/S1744309112030151 22949203

27. Pasdeloup D, Poisson N, Raux H, Gaudin Y, Ruigrok RW, et al. (2005) Nucleocytoplasmic shuttling of the rabies virus P protein requires a nuclear localization signal and a CRM1-dependent nuclear export signal. Virology 334: 284–293. 15780878

28. Desmezieres E, Maillard AP, Gaudin Y, Tordo N, Perrin P (2003) Differential stability and fusion activity of Lyssavirus glycoprotein trimers. Virus Res 91: 181–187. 12573496

29. Rose NF, Roberts A, Buonocore L, Rose JK (2000) Glycoprotein exchange vectors based on vesicular stomatitis virus allow effective boosting and generation of neutralizing antibodies to a primary isolate of human immunodeficiency virus type 1. J Virol 74: 10903–10910. 11069984

30. Chothia C, Lesk AM (1986) The relation between the divergence of sequence and structure in proteins. EMBO J 5: 823–826. 3709526

31. Wood TC, Pearson WR (1999) Evolution of protein sequences and structures. J Mol Biol 291: 977–995. 10452901

32. Benmansour A, Leblois H, Coulon P, Tuffereau C, Gaudin Y, et al. (1991) Antigenicity of rabies virus glycoprotein. J Virol 65: 4198–4203. 1712859

33. Kongsuwan K, Cybinski DH, Cooper J, Walker PJ (1998) Location of neutralizing epitopes on the G protein of bovine ephemeral fever rhabdovirus. J Gen Virol 79 (Pt 11): 2573–2581. 9820132

34. Prehaud C, Coulon P, LaFay F, Thiers C, Flamand A (1988) Antigenic site II of the rabies virus glycoprotein: structure and role in viral virulence. J Virol 62: 1–7. 2446011

35. Vandepol SB, Lefrancois L, Holland JJ (1986) Sequences of the major antibody binding epitopes of the Indiana serotype of vesicular stomatitis virus. Virology 148: 312–325. 2417417

36. Olsson MHM, Soendergard CR, Rostkowski M, Jensen JH (2011) PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa predictions. Journal of Chemical Theory and Computation 7: 525–537.

37. Marriott AC (2005) Complete genome sequences of Chandipura and Isfahan vesiculoviruses. Arch Virol 150: 671–680. 15614433

38. Vasilakis N, Widen S, Travassos da Rosa AP, Wood TG, Walker PJ, et al. (2013) Malpais spring virus is a new species in the genus vesiculovirus. Virol J 10: 69. doi: 10.1186/1743-422X-10-69 23497016

39. DuBois RM, Vaney MC, Tortorici MA, Kurdi RA, Barba-Spaeth G, et al. (2013) Functional and evolutionary insight from the crystal structure of rubella virus protein E1. Nature 493: 552–556. doi: 10.1038/nature11741 23292515

40. Modis Y (2014) Relating structure to evolution in class II viral membrane fusion proteins. Curr Opin Virol 5: 34–41. doi: 10.1016/j.coviro.2014.01.009 24525225

41. Kabsch W (2010) Xds. Acta Crystallogr D Biol Crystallogr 66: 125–132. doi: 10.1107/S0907444909047337 20124692

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

43. 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

44. Emsley P, Lohkamp B, Scott WG, Cowtan K (2010) Features and development of Coot. Acta Crystallogr D Biol Crystallogr 66: 486–501. doi: 10.1107/S0907444910007493 20383002

45. Terwilliger TC, Grosse-Kunstleve RW, Afonine PV, Moriarty NW, Adams PD, et al. (2008) Iterative-build OMIT maps: map improvement by iterative model building and refinement without model bias. Acta Crystallogr D Biol Crystallogr 64: 515–524. doi: 10.1107/S0907444908004319 18453687

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

Článek vyšel v časopise

PLOS Pathogens

2015 Číslo 3

Nejčtenější v tomto čísle

Tomuto tématu se dále věnují…

Kurzy Doporučená témata Časopisy
Zapomenuté heslo

Nemáte účet?  Registrujte se

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.


Nemáte účet?  Registrujte se

VIRTUÁLNÍ ČEKÁRNA ČR Jste praktický lékař nebo pediatr? Zapojte se! Jste praktik nebo pediatr? Zapojte se!