1. Structural Biology and Molecular Biophysics

Structural analyses at pseudo atomic resolution of Chikungunya virus and antibodies show mechanisms of neutralization

  1. Siyang Sun
  2. Ye Xiang
  3. Wataru Akahata
  4. Heather Holdaway
  5. Pankaj Pal
  6. Xinzheng Zhang
  7. Michael S Diamond
  8. Gary J Nabel
  9. Michael G Rossmann  Is a corresponding author
  1. Purdue University, United States
  2. National Institute of Allergy and Infectious Diseases, National Institutes of Health, United States
  3. Washington University School of Medicine, United States
https://doi.org/10.7554/eLife.00435
Share this article
Cite this article
  1. Siyang Sun
  2. Ye Xiang
  3. Wataru Akahata
  4. Heather Holdaway
  5. Pankaj Pal
  6. Xinzheng Zhang
  7. Michael S Diamond
  8. Gary J Nabel
  9. Michael G Rossmann
(2013)
Structural analyses at pseudo atomic resolution of Chikungunya virus and antibodies show mechanisms of neutralization
eLife 2:e00435.
https://doi.org/10.7554/eLife.00435
  2. Download BibTeX
  3. Download .RIS

Abstract

A 5.3 Å resolution, cryo-electron microscopy (cryoEM) map of Chikungunya virus-like particles (VLPs) has been interpreted using the previously published crystal structure of the Chikungunya E1-E2 glycoprotein heterodimer. The heterodimer structure was divided into domains to obtain a good fit to the cryoEM density. Differences in the T = 4 quasi-equivalent heterodimer components show their adaptation to different environments. The spikes on the icosahedral 3-fold axes and those in general positions are significantly different, possibly representing different phases during initial generation of fusogenic E1 trimers. CryoEM maps of neutralizing Fab fragments complexed with VLPs have been interpreted using the crystal structures of the Fab fragments and the VLP structure. Based on these analyses the CHK-152 antibody was shown to stabilize the viral surface, hindering the exposure of the fusion-loop, likely neutralizing infection by blocking fusion. The CHK-9, m10 and m242 antibodies surround the receptor-attachment site, probably inhibiting infection by blocking cell attachment.

https://doi.org/10.7554/eLife.00435.001

eLife digest

The Chikungunya virus is carried by mosquitos and can cause a number of diseases in humans including encephalitis, which can be fatal in some cases, and severe arthritis. A recent mutation in the E1 protein of the virus has allowed it to efficiently reproduce in a different species of mosquitos, leading to a Chikungunya epidemic in Réunion Island in 2005 and the subsequent infection of millions of individuals in Africa and Asia. The virus also has the potential to spread to many areas of Europe and the Americas.

Chikungunya virus has a single-stranded RNA genome that codes for four non-structural proteins and five structural proteins. Based on this knowledge it has been possible to develop virus-like particles that can be used to immunise non-human primates against Chikungunya infection by inducing antibody production. However, the development of vaccines for Chikungunya in humans will require a deeper understanding of how these antibodies produced by the vaccine interact with the virus and more detailed information about the structures of the virus and antibodies.

Sun et al. have used two techniques – X-ray crystallography and electron cryo-microscopy – to determine the structure of Chikungunya virus-like particles, and to obtain new insights into the interactions of these particles with four related antibodies. Electron cryo-microscopy was used to figure out the structure of the particles at near atomic resolution, and X-ray crystallography was used to determine the atomic resolution structures of two of the four Fab antibodies that neutralize the Chikungunya virus. Electron cryo-microscopy was also used to probe the complex formed by the interactions between the virus-like particles and the antibodies.

Sun et al. were able to identify the likely viral receptor site that is blocked by three of the antibodies when they neutralize the virus; the fourth antibody is thought to act by immobilizing one of the domains of protein E2, thereby hiding the “fusion loop” that allows the virus to enter and infect human tissue. It is hoped that these findings will contribute to efforts to combat the spread of the Chikungunya virus worldwide.

https://doi.org/10.7554/eLife.00435.002

Introduction

Chikungunya virus (CHIKV) is a mosquito-transmitted viral pathogen that causes fever, myalgia, rash, and severe arthritis in humans (Powers and Logue, 2007; Simon et al., 2008). The first reported human CHIKV infections occurred in East Africa in 1952 (Robinson, 1955). Prior to the 2005 epidemic on Réunion Island, CHIKV was not regarded as a highly prevalent virus. An adaptive mutation in the E1 protein (E1-A226V) that allowed CHIKV to replicate more efficiently in Aedes albopictus is considered to be the primary reason for its recent extensive spread, infecting millions of individuals in Africa and Asia (Tsetsarkin et al., 2007; Kumar et al., 2008; Santhosh et al., 2008). In some CHIKV-infected patients, severe damage in joint tissues can cause debilitating chronic arthritis. In the recent outbreaks, a change in pathogenesis was observed, with some cases progressing to fatal encephalitis. An autochthonous CHIKV outbreak in Italy in 2007 and the presence of the Aedes albopictus vector in many areas of Europe and the Americas have raised concern of further spread of the virus. Currently, there is no vaccine or antiviral agent approved for use in humans.

CHIKV belongs to the alphavirus genus of the Togaviridae family (Kuhn, 2007). Alphaviruses are a group of positive-sense, single-stranded RNA, enveloped viruses transmitted by arthropods (Griffin, 2007). The alphavirus genome encodes four non-structural and five structural proteins. The non-structural proteins are required for virus replication, protein modification, and immune antagonism. The structural proteins (capsid-E3-E2-6K-E1) are synthesized as a polyprotein from a subgenomic promoter, and are cleaved post-translationally into separate proteins by an autoproteinase and signalase. The E1 glycoprotein participates in cell fusion (Lescar et al., 2001), whereas the E2 glycoprotein binds to cellular receptors (Smith et al., 1995; Zhang et al., 2005) and initiates clathrin-dependent endocytosis (Solignat et al., 2009). Virus core assembly is initiated by interactions between the genomic RNA and the nuclear capsid protein (NCP) (Tellinghuisen et al., 1999; Tellinghuisen and Kuhn, 2000; Linger et al., 2004) in the cytoplasm. The E3 protein is essential for the proper folding of p62, the precursor to E2, and the formation of the p62-E1 heterodimer (Mulvey and Brown, 1995; Carleton et al., 1997). Although E3 remains part of mature Semliki Forest virus (SFV) and Venezuelan equine encephalitis virus (VEEV), apparently, it is not a component of mature CHIKV (Simizu et al., 1984). The small 6 kDa protein, 6K, associates with the p62-E1 heterodimer and is transported to the plasma membrane prior to assembly. The 6K protein facilitates particle morphogenesis but is not stoichiometrically incorporated into virions (Gaedigk-Nitschko and Schlesinger, 1990, 1991).

Alphaviruses are icosahedral particles that have T = 4 quasi-icosahedral symmetry (Paredes et al., 1993; Venien-Bryan and Fuller, 1994; Cheng et al., 1995; Zhang et al., 2002, 2011; Kostyuchenko et al., 2011). The ectodomain forms 80 spikes, each consisting of three copies of E1-E2 heterodimers. There are 20 icosahedral “i3” spikes, situated on icosahedral 3-fold axes, and 60 quasi-3-fold “q3” spikes in general positions containing a quasi-3-fold axis consistent with T = 4 symmetry (Figure 1A). Thus, each icosahedral asymmetric unit contains one complete q3 spike and one-third of an i3 spike. E1 of CHIKV is a 439 amino acid protein with an N-linked glycosylation site at residue Asn141. The ectodomain is formed by the 404 N-terminal residues, followed by a 30 residue transmembrane (TM) helix and a five amino acid cytoplasmic domain. The E1 ectodomain consists of three β-barrel domains (Lescar et al., 2001; Gibbons et al., 2004; Li et al., 2010; Voss et al., 2010). Domain I is spatially located between domains II and III with the fusion peptide lying at the distal end of domain II. E2 of CHIKV is 423 amino acids in length and has two N-linked glycans at positions Asn263 and Asn345. The E2 364 residue ectodomain is followed by a 26-residue TM helix and a 33-residue cytoplasmic domain. The E2 ectodomain consists of three distinct immunoglobulin-fold domains with domain A, the putative receptor-binding domain, lying between domains B and C (Li et al., 2010; Voss et al., 2010). In the mature virus, domain B covers the fusion loop in domain II of E1. The virus becomes fusogenic in acid pH (Wahlberg and Garoff, 1992) when sets of three E1 molecules combine to make a fusogenic trimer after expelling the three E2 proteins from the center of the spike. The three fusion loops are then exposed outwards for insertion into the host cell membrane (Gibbons et al., 2004). Similar large conformational changes occur in other enveloped viruses (Kielian and Rey, 2006) although little is known of the structural intermediates through which these changes take place.

Figure 1
Download asset Open asset
Structure of the CHIK VLPs.

(A) Surface-shaded figure of ectodomain (left) and surface-shaded figure of nucleocapsid (right), colored according to the radial distance from the center of the virus. White triangles indicate one icosahedral asymmetric unit. (B) Cross-section of the virus showing density above 1.5 σ also colored according to the radial distance from the center of the virus. (C) Resolution of β-strands in the E1 domain III.

https://doi.org/10.7554/eLife.00435.003

The nucleocapsid core consists of 30 hexamers around each 2-fold axis and 20 pentamers at each 5-fold vertex of the nuclear capsid (Figure 1A), and is separated from the ectodomains of E1 and E2 by a lipid membrane (Figure 1B). The N-terminal 118 amino acids of Sindbis (SINV) (Choi et al., 1991) and SFV (Choi et al., 1997) nuclear capsid contain positively charged residues that associate with the genomic RNA (Zhang et al., 2002). A similar positively charged N-terminal region exists in CHIKV NCP and functions to recognize the viral genome for packaging during NCP assembly (Perera et al., 2001; Hong et al., 2006). The carboxy-terminal region of the NCP has a chymotrypsin-like fold and acts as an auto-catalytic proteinase during processing of the structural polyprotein (Hahn and Strauss, 1990; Choi et al., 1991, 1997). In addition, there is a hydrophobic pocket in the C-terminal domain of the NCP that binds the 33 amino acid cytoplasmic tail of E2, effectively linking the proteins of the inner core with those on the surface of the virion (Lee et al., 1996; Skoging et al., 1996). This association of E2 with the nuclear capsid is essential for virus assembly (Owen and Kuhn, 1997).

Many studies have analyzed the interaction of antibodies with alphaviruses (Meyer and Johnston, 1993; Griffin et al., 2001; Hunt et al., 2010). Most neutralizing antibodies bind to E2, as this glycoprotein is more exposed on the viral surface than E1. The mode of neutralization differs considerably among antibodies, with some blocking attachment to cells and others inhibiting the conformational changes required for fusion (Hernandez et al., 2008). The only published structural studies of alphavirus-antibody complexes were a cryoEM analysis of Ross River virus complexed with an antibody that blocked receptor attachment (Smith et al., 1995) and a cryoEM investigation of SINV complexed with an antibody that inhibited the conformational changes required to form fusogenic E1 trimers (Hernandez et al., 2008).

Expression of the CHIKV structural C-E3-E2-6K-E1 polyprotein gives rise to virus-like particles (VLPs) that have been used to immunize and protect non-human primates against CHIKV infection (Akahata et al., 2010). Non-human primates immunized with these VLPs developed high-titer neutralizing antibodies that protected against viremia after high-dose challenges. Moreover, immunodeficient mice that were passively administered purified IgG from immunized non-human primates were protected against CHIKV challenge (Couderc et al., 2009). Previously, we published the structure of the CHIK VLPs at a resolution of 18 Å using cryoEM reconstruction of single particles (Akahata et al., 2010). Here, we report a 5.3 Å three-dimensional high-resolution reconstruction of the CHIK VLP. We also describe pseudo-atomic resolution structures of CHIK VLPs complexed with Fab fragments from four different neutralizing mouse monoclonal antibodies (MAb CHK-9, CHK-152, m10 and m242). These studies are the first structural investigations of CHIKV complexed with neutralizing MAbs, based on 15 Å resolution cryoEM reconstructions of the complexes, the 5.3 Å resolution structure of the VLPs and 3 Å resolution crystallographic studies of the Fab fragments. These structures suggest that MAb CHK-152 neutralizes infectivity by inhibiting fusion whereas the other MAbs probably neutralize by blocking attachment.

Results and discussion

CryoEM structure of CHIK VLPs

The structure of CHIK VLPs was determined to a resolution of ∼5.3 Å. The trimeric appearance of individual spikes was recognizable immediately, along with the inclined positions of the three E1 molecules surrounding each spike. Further inspection of the E1 and E2 molecules showed β-barrel structures surrounding empty cavities, thus yielding a good indication of the positions of each domain. Although resolution of individual β-strands was possible only in some instances (Figure 1C), identification of individual side chains was not possible. The E3 glycoprotein, which is observed in some alphavirus structures (Zhang et al., 2011), was not apparent in the CHIK VLP cryoEM density map. Densities for the E2 B-domain were weak in all four E1-E2 heterodimers, consistent with the crystal structures of the CHIKV E1-E2 heterodimer (Voss et al., 2010) and the SINV trimeric spike (Li et al., 2010).

Within the CHIK VLP structure, the four quasi-equivalent chymotrypsin-like nucleocapsid proteins were organized as pentamers around the 5-fold vertices and hexamers around the 2-fold icosahedral symmetry axes (Figure 1A). The cryoEM density height of the nucleocapsid proteins was roughly 5% less than that of the glycoproteins. The glycoproteins and the nuclear capsid were separated by a ∼45 Å-wide lipid membrane. The external and internal lipid leaflets had a density of about 5 σ, whereas the protein density of the glycoproteins and the NCPs had a density of about 7 σ. The two lipid leaflets were separated by a ∼15 Å wide region, where the density was less than 1.5 σ, representing the loosely packed aliphatic chains of the lipid molecules (Figure 1B). The membrane was traversed by 60 × 4 pairs of α helices representing the E1 and E2 carboxy-terminal regions (Figure 1B). The quality of the membrane density in the icosahedrally averaged map demonstrates that the symmetry of the glycoproteins and nuclear capsid impose icosahedral symmetry on the flexible membrane by confining it to a limited space.

Refinement and verification of the structure

A number of procedures have been developed for the refinement of atomic parameters when the resolution of the structural data is insufficient to follow the polypeptide chains with certainty or to recognize the identity of amino acids (Chapman, 1995; Volkmann and Hanein, 1999; Tama et al., 2004; Fabiola and Chapman, 2005; Goddard et al., 2007; Topf et al., 2008; Trabuco et al., 2008; Lasker et al., 2010; Schröder et al., 2010; Zhu et al., 2010). These procedures increase the number of refinable parameters beyond those necessary to position and orient a rigid known structure into a relatively low-resolution cryoEM density map, thus reducing the ratio of observed data points (e.g. structure amplitudes) to refinable parameters. However severe restraints are placed on the geometry of the structure (thus providing additional observational data) while improving agreement between observed and calculated data in real or reciprocal space. Verification in reciprocal space can be achieved by tests on data that had been excluded in the refinement process, for instance by calculating Rfree (Brünger, 1992). Verification in real space can be established by observing anticipated structural features that had not been considered in the restrained refinement process.

In the analysis of the CHIKV VLPs described here, the initial interpretation was performed by fitting the structure of the CHIKV E1-E2 heterodimer (Voss et al., 2010) into the cryoEM map of the CHIKV VLPs assuming strict T = 4 symmetry between quasi-equivalent structures. Subsequently, the number of refinable parameters was increased by breaking the structure into domains, based on the assumption that the domains likely are more rigid than the hinges in the polypeptide chain that hold the domains together. In addition, constraints imposed by the assumption of strict T = 4 symmetry were abandoned, thus further multiplying the number of refinable parameters by four. Verification came from the improved correlation of the domain-fitted structure with the amino acid sequences, the improved correlation between density heights of quasi-equivalent parts of the structure with each other, and from the emergence of features such as the position of the N-linked carbohydrate moieties.

Fitting of the CHIKV E1-E2 crystal structure into the cryoEM density map

In contrast to the interpretation of the cryoEM density of VEEV (Zhang et al., 2011), which had a resolution of 4.8 Å, we were able to use the crystallographically determined structure of the glycoproteins to aid in the analysis. The crystal structure of CHIKV E1-E2 (Voss et al., 2010) was fitted into the 5.3 Å cryoEM map using the program EMfit to maximize the average value of the density (sumf) at all atomic positions (Table 1A), to minimize the steric clash between symmetry related structures, and to minimize the percentage of atoms in low density (Rossmann et al., 2001). Initially, the E1-E2 heterodimer crystal structure was fitted as a rigid body generating the four quasi-equivalent structures at positions #1, #2, #3, and #4 (Figure 2A) while refining the best positions of the quasi-symmetry axes (see Materials and Methods), and using the icosahedral symmetry to generate the structure in neighboring icosahedral asymmetric units. The resultant locations of the four quasi-equivalent CHIKV heterodimers were used as starting positions for refinement of individual domains. The average height of the density at all fitted atoms (sumf) improved by ∼12% for the independently fitted domains. The value of sumf was at least 38% lower for the B domains than for any of the other domains. This is consistent with the disorder of the B domain in the crystal structure of the SINV spike (Li et al., 2010) and the large temperature factors for the atoms in the B domain in the crystal structure of the CHIKV heterodimer (Voss et al., 2010). The height of the average value of the cryoEM density taken over all atomic positions within a molecular domain was higher at the quasi-equivalent position #4 than at the other three quasi-equivalent positions (Table 1A). This might be due to crowding around the 5-fold symmetry axes, which reduces the variability in the heterodimer position.

Table 1

Quality of fitting the atomic structural fragments into the cryoEM density

https://doi.org/10.7554/eLife.00435.004
A) Average height of the densities at the atomic positions (sumf) (Rossmann et al., 2001) on fitting the E1E2 heterodimer CHIKV structure into the cryoEM map at four quasi equivalent positions
All atomsT = 4 fittingIndependent domain fitting
#1#2#3#4Average#1#2#3#4Average
E1*I16.115.817.017.716.618.217.018.519.918.5
II15.717.218.816.717.018.018.619.018.518.6
III15.714.215.218.015.718.715.918.220.718.3
E2*A15.518.019.017.517.418.817.819.919.418.8
B10.09.49.69.69.810.911.311.712.411.7
C15.315.419.019.116.821.019.920.222.420.5
D13.316.415.816.115.416.417.216.016.816.5
Average14.915.817.016.816.117.817.318.219.018.1
Main chain atoms onlyT = 4 fittingIndependent domain fitting
#1#2#3#4Average#1#2#3#4Average
E1*I20.518.421.522.920.823.221.622.626.323.4
II17.020.923.020.220.321.723.323.522.822.8
III18.915.317.522.318.522.317.922.626.622.4
E2*A16.521.122.219.419.822.020.323.422.622.1
B10.79.49.910.110.010.811.213.013.412.1
C18.719.527.426.923.127.727.127.430.828.3
D13.418.718.518.517.319.419.819.119.019.3
Average16.517.620.020.018.521.020.221.723.121.5
B) Average height of the densities at the atomic positions (sumf) (Rossmann et al., 2001) on fitting the CHIKV homology model of the capsid protein into the cryoEM map at the four quasi equivalent positions
All atomsT = 4 fittingIndependent domain fitting
#1#2#3#4Average#1#2#3#4Average
D1†14.913.313.916.114.616.114.615.417.315.8
D2†17.017.418.118.117.618.319.117.419.018.6
Average16.015.416.017.216.117.417.116.518.217.3
Main chain atoms onlyT = 4 fittingIndependent domain fitting
#1#2#3#4Average#1#2#3#4Average
D1†18.516.517.818.917.920.416.919.222.719.8
D2†21.321.521.423.221.922.824.220.425.323.2
Average19.919.019.621.119.921.821.019.924.121.6
C) Average height of the densities at the atomic positions (sumf) (Rossmann et al., 2001) on fitting the CHIKV model of the transmembrane protein into the cryoEM map at the four quasi equivalent positions
All atomsT = 4 fittingIndependent domain fitting
#1#2#3#4Average#1#2#3#4Average
E1‡13.714.613.218.615.013.914.414.717.515.1
E2‡14.712.811.817.514.215.114.113.718.915.5
Average14.413.512.317.914.614.514.214.218.215.3
Main chain atoms onlyT = 4 fittingIndependent domain fitting
#1#2#3#4Average#1#2#3#4Average
E1‡16.217.616.024.118.515.717.316.723.918.4
E2‡16.715.013.620.716.517.716.715.823.718.5
Average16.516.014.422.017.516.717.016.223.818.4
  1. *

    Domain definition: E1 I (residues 1-36, 132-168, 273-293), II (residues 37-131, 169-272), III (residues 294-393). E2 A (residues 16-134), B (residues 173-231), C (residues 269-342), D (residues 7-15, 135-172, 232-268).

  2. †

    Domain definition: D1(residues 119 to 183), D2(residues 184 to 267).

  3. ‡

    Domain definition: E1 (residues 394-439), E2(residues 343-423).

Figure 2
Download asset Open asset
Diagrammatic representation of contacts between amino acid residues in adjacent subunits.

(A) Diagrammatic organization of the E1 and E2 subunits according to T = 4 icosahedral symmetry. The white numbers show the sequence in which the four independent subunits were generated. The black capital letters indicate the sequence of generating seven of the icosahedral asymmetric units. Icosahedral symmetry elements are shown as filled triangles and ellipses. Quasi-symmetry elements are shown as red outlined triangles and ellipses. (B) Contacts between the E1 (blue) and E2 (green) molecules. (C) Contacts between NCPs. (B) and (C) show the number of contacts between the indicated molecules. The quasi T = 4 related positions #1, #2, #3 and #4 are indicated, prefixed by their icosahedral symmetry identification A, B, C and D. The center of the “i3” icosahedral spike is indicated by a filled red triangle and the center of the “q3” quasi-3-fold spike is indicated by an outlined red triangle.

https://doi.org/10.7554/eLife.00435.005

Once the individual domains were fitted (Figure 3A), the bond geometry between the carboxy end of one domain and the amino end of the next domain were regularized. The distance between these ends prior to regularization (Table 2) did not exceed 6.9 Å (excluding the flexible domain B), indicating that the quality of fit was reasonable. Comparison of the hinge angles between domains in the four quasi-equivalent positions (Table 3) did not exceed 13°, except for defining the orientation of the B domain in E2. Considering that the length of the longest domain (DII of E1) is about 50 Å, this implies a relative difference of about 6 Å for atoms furthest from the hinge, corresponding to a maximum movement of about three standard deviations between superimposed structures.

Figure 3
Download asset Open asset
Fit of the atomic structure into the cryoEM density.

(A) Fit of the E1E2 heterodimer, (B) Fit of the capsid protein. The quasi equivalent subunit closest to the icosahedral 5-fold axis was chosen for display.

https://doi.org/10.7554/eLife.00435.006
Table 2

Distances in Å between N and C termini of fitted domains before regularization

https://doi.org/10.7554/eLife.00435.007
Quasi equivalent positions
Domains#1#2#3#4
E1I 36 to II 374.65.65.53.5
II 131 to I 1324.84.14.44.1
I 168 to II 1692.86.45.45.1
II 272 to I 2735.84.73.92.1
I 293 to III 2945.46.93.75.0
E2D 15 to A 163.86.35.75.9
A134 to D1353.15.04.42.3
D172 to B1733.48.83.05.7
B231 to D2326.27.75.97.2
NCPD268 to C2695.15.24.14.0
(D1)183 to (D2)1841.22.21.61.6
Table 3

Hinge angle change between domains in different heterodimers

https://doi.org/10.7554/eLife.00435.008
Hinge angle change in E1 between domain I and domain II (in degrees)
Position #1Position #2Position #3Position #4
Crystal structure-14.4-6.67.5-8.3
Position #18.09.36.8
Position #23.12.1
Position #3-3.3
Hinge angle change in E1 between domain I and domain III (in degrees)
Position #1Position #2Position #3Position #4
Crystal structure-7.3-8.0-9.42.1
Position #1-3.510.37.7
Position #213.09.2
Position #37.6
Hinge angle change in E2 between domain A and domain B (in degrees)
Position #1Position #2Position #3Position #4
Crystal structure-8.410.422.2-22.3
Position #1-14.125.5-29.7
Position #226.6-28.6
Position #3-26.3
Hinge angle change in E2 between domain A and domain C (in degrees)
Position #1Position #2Position #3Position #4
Crystal structure-11.5-10.4-7.6-17.1
Position #1-4.8-8.6-7.3
Position #27.28.3
Position #3-14.8

After independent fitting of the domains, the root mean square (rms) distance between equivalent Cα atoms when making pairwise comparisons between E1-E2 heterodimers ranged from 0.5 Å to 2.4 Å (Table 4). This showed that the structures of the four quasi-equivalent parts of the model were more alike to each other than to the crystal structure of the CHIKV heterodimer (Voss et al., 2010). As such similarity is unlikely to occur by accident, this observation helped to confirm the validity of the process used to interpret the cryoEM map. The systematic change of the E1-E2 heterodimers in the VLPs compared to the crystal structure is the result of accommodating the different quasi-equivalent environments. In contrast to the results for the related alphavirus, VEEV (Zhang et al., 2011), the presence of differences larger than the standard deviation between equivalenced Cα atoms between the four quasi-equivalent structural components of the CHIKV VLPs explains why averaging of the quasi equivalent density did not improve the map.

Table 4

RMS deviation in Å between quasi equivalent Cα atoms

https://doi.org/10.7554/eLife.00435.009
E1#1#2#3#4Xtal
#10.001.211.251.081.69
#20.001.281.271.52
#30.001.101.68
#40.001.17
Xtal0.00
E2#1#2#3#4Xtal
#10.002.181.791.901.50
#20.002.342.382.17
#30.002.281.84
#40.001.92
Xtal0.00
E1&E2#1#2#3#4Xtal
#10.001.941.651.682.28
#20.002.162.142.24
#30.001.872.37
#40.001.88
Xtal0.00
Capsid#1#2#3#4Xtal
#10.000.490.730.710.48
#20.001.090.940.71
#30.000.630.56
#40.000.52
Xtal0.00

The quality of the side chain placements was determined by correlating their densities in the four quasi-equivalent positions. The density associated with each residue was taken as the average density of all the non-hydrogen atoms in that residue. A high correlation coefficient (see Materials and Methods for definition) indicates that the variation of cryoEM density along each polypeptide chain was similar. As the average density for each residue included all atoms in the side chain, a high correlation also implied similarity of the side chain densities. These correlation coefficients (Table 5A) varied from 0.44 to 0.30, which was significantly greater than the values (0.40–0.08) for the initial T = 4 rigid body fitting. Furthermore, there was a reasonable correlation of the average crystallographic temperature factors for each amino acid residue of the CHIKV heterodimer crystal structure (Voss et al., 2010) with the height of the cryoEM density at each amino acid position (Table 5A).

Table 5

Correlation between amino acid sequence and cryoEM density

https://doi.org/10.7554/eLife.00435.010
A) Correlation between cryoEM densities in the four quasi equivalent positions of the E1 and E2 glycoproteins.
All atomsT = 4 fittingIndependent domain fitting
1234B1234B
11.000.170.230.080.221.000.400.440.400.21
21.000.280.380.301.000.410.300.28
31.000.400.281.000.360.28
41.000.191.000.15
B1.001.00
Main chain atoms onlyT = 4 fittingIndependent domain fitting
12#34B1234B
11.000.110.200.110.191.000.400.410.410.22
21.000.310.300.301.000.360.320.26
31.000.410.271.000.370.26
41.000.161.000.12
B1.001.00
B) Correlation between cryoEM densities of the four quasi equivalent capsid proteins.
All atomsT = 4 fittingIndependent domain fitting
#1#2#3#4#1#2#3#4
#11.000.410.180.171.000.200.190.30
#21.000.260.251.000.450.42
#31.000.131.000.41
#41.001.00
Main chain atoms onlyT = 4 fittingIndependent domain fitting
#1#2#3#4#1#2#3#4
#11.000.370.160.261.000.220.170.33
#21.000.310.271.000.340.43
#31.000.151.000.31
#41.001.00
C) Correlation between cryoEM densities of the four quasi equivalent E1&E2 TM and endodomain regions.
All atomsT = 4 fittingIndependent domain fitting
#1#2#3#4#1#2#3#4
#11.000.430.170.411.000.370.340.35
#21.000.240.381.000.440.34
#31.000.201.000.46
#41.001.00
Main chain atoms onlyT = 4 fittingIndependent domain fitting
#1#2#3#4#1#2#3#4
#11.000.280.080.271.000.230.110.26
#21.000.170.431.000.400.36
#31.000.141.000.50
#41.001.00

The N-linked carbohydrate sites

The fitting of the independent domains of the CHIKV heterodimer into the cryoEM map was further validated by finding the densities associated with N-linked glycans, after setting to zero all densities at grid points within 3.0 Å of every fitted atom. The remaining densities should correspond to the sugar moieties on E1 and E2 (Pletnev et al., 2001). The shape and size of the carbohydrate moiety for each of the three N-linked glycosylation sites is consistent in each of the four quasi-equivalent positions. The average volume of the carbohydrates at the four E1 (Asn141) sites in DI of E1 is 258 Å3 with an average height of 2.2 σ, with most of the noise density not being greater than 1.0 σ. This would correspond to about one sugar moiety. The densities corresponding to the four quasi-equivalent E2 (Asn263) sites in the E2 connecting ribbon were somewhat fragmented, making it difficult to determine their volumes. However, at lower (15 and 10 Å) resolution, the densities were well defined. The volume of the density and average height associated with the four quasi-equivalent sites at E2 (Asn345) in domain C of E2 were as large as the sites at E1 (Asn141), but because of their proximity to the viral membrane, their boundaries were uncertain. The centers of the carbohydrate moieties were ∼6 Å from the N atom of the glycosylated Asn residues.

Interactions between the E1 and E2 glycoproteins

Contacts between the E1 and E2 glycoproteins were determined (Materials and Methods) for the T = 4 fitted structure and compared with the improved, independently fitted, domain structure (Table 6). The E1 and E2 glycoproteins that form the four heterodimers in the VLPs have more contacts with each other in the independently fitted structure than in the crystallographic heterodimer structure. Particularly striking was the increase in hydrophobic contacts in the independently fitted results, which demonstrated how the heterodimer adapted itself to the quasi-equivalent symmetry of the VLPs compared to the constraints imposed by the crystal lattice.

Table 6

Inter-molecular contacts (See methods for descriptions)

https://doi.org/10.7554/eLife.00435.011
Spike positionTotalHydrophobicPossible H bondsPossible salt bridges
A) Number of atom-to-atom contacts between E1-E2 in heterodimer
i3#123499 (42%)34 (13%)0
q3#2254112 (44%)25 (10%)0
#317959 (33%)33 (18%)0
#4302118 (39%)39 (13%)4 (1%)
Crystal861 (1%)32 (37%)1 (1%)
B) Number of atom-to-atom contacts of E2 with other E2s within a spike
i3#122946 (20%)30 (13%)82 (36%)
q3#2708 (11%)10 (14%)33 (47%)
#34714 (30%)6 (13%)10 (21%)
#474118 (23%)24 (46%)7 (9%)
C) Number of atom-to-atom contacts between glycoprotein spikes
q3 to q3A#3 - D #412342 (34%)21 (17%)7 (5%)
A#4 - D#413052 (40%)18 (14%)1 (1%)
i3 to q3A#2 - B#15827 (47%)7 (12%)0
A#1 - A#28829 (33%)19 (22%)0
A#1 - A#3118 (73%)2 (8%)0
D) Number of atom-to-atom contacts between E1 and E2 within the spikes.
q3A#2 - A#33016 (34%)5 (17%)0
A#3 - A#4208 (40%)3 (14%)0
i3A#4 - A#25625 (47%)11 (12%)0
A#1 - B#1224112 (50%)21 (9%)22 (10%)
E) Number of atom-to-atom contacts between capsid proteins*
A#2 - B#17920 (25%)13 (16%)15 (19%)
A#1 - A#3113 (27%)1 (9%)3 (27%)
A#3 - D#218351 (28%)22 (12%)38 (21%)
A#4 - D#4195 (26%)2 (11%)8 (42%)
A#1 - B#1209 (45%)2 (10%)2 (10%)
F) Number of atom-to-atom contacts between E1 & E2 in TM region
#110060 (60%)3 (3%)0
#28243 (52%)8 (10%)0
#32012 (60%)1 (5%)0
#42824 (86%)00
  1. *

    See Figure 2.

Unlike the T = 4 fitted structure, the independently fitted structure lacked exact symmetry between the i3 and q3 spikes. Nevertheless, the q3 spike had apparent, although not exact, 3-fold symmetry within itself and with its surroundings. A comparison of the i3 and q3 structures showed that the E1 molecules wrap around the central E2 molecules more tightly in the i3 spike as indicated by the increased number of interactions between the E2 molecules in the i3 spike compared to the q3 spike. Thus, the i3 trimers in the VLPs were more compact with a greater number of interactions between its component molecules than the q3 trimer. The greater compactness of the i3 spike was confirmed by determining the radius of gyration of the i3 and the q3 spikes as a function of the radial distance from the viral center (Materials and Methods). At the base of the spike, where the E1 molecules make contacts between spikes, the radius of gyration was up to 5.8 Å less for the i3 spike than for the q3 spike (Table 7). At the top of the spike, the i3 and q3 spikes showed less difference. In addition, the whole of the q3 spike was translated away from the center of the VLPs by ∼1.8 Å. Thus, in comparison with the 20 i3 spikes, the 60 q3 spikes may be poised to release the E2 molecules from the center of the spikes in preparation for the formation of the E1 fusogenic homotrimers early in the infection (Li et al., 2010).

Table 7

Spike radii of gyration

https://doi.org/10.7554/eLife.00435.012
Distance from particle center (Å)Radii of gyration (Å)# of atoms
i3q3SINV* crystali3q3SINV* crystal
25558.361.063.3423462
26553.557.354.9131126118
27541.146.942.6125122101
28530.133.130.1849080
29533.234.233.8858091
30533.331.931.3908982
31530.731.224.5708092
32526.725.94147
33531.029.011
Overall40.542.441.9669669626
  1. *

    Sindbis virus (SINV)

Prior studies (Li et al., 2010) have suggested that the trigger for movement of domain B to uncover the fusion loop on E1 at low pH might be the protonation of the highly conserved His169 and His256 residues in SINV on the β-ribbon that connects the ends of the B domain to the A and C domains of E2 (Li et al., 2010). In the CHIK VLP structure, His170 (equivalent to SINV His 169) does not interact closely with other negatively charged residues that would be required for a pH-sensing switch whereas CHIKV His 256 (equivalent to SINV His256) contacts Glu166.

Fitting of a homology model of CHIKV NCP structure into the cryoEM density map

The crystal structures of SINV (Choi et al., 1991; Tong et al., 1992) and SFV (Choi et al., 1997) NCPs had been determined previously. Of these, the amino acid sequence of the NCP of SFV has greater similarity to CHIKV. Therefore, a homology model for CHIKV NCP (residues 119 to 267) was generated based on the crystal structure of the SFV NCP (Protein Data Bank accession number 1VCP). This homology structure was fitted into the CHIK VLP cryoEM density map by assuming T = 4 symmetry, using the EMfit program (Table 1B). The positions of the quasi-symmetry axes were refined to optimize the fit to the density and to minimize clashes between symmetry related NCP structures. The refinement found that the placement of the quasi-symmetry axes did not change significantly from those optimized while fitting the glycoproteins. The orientation and placement of the homology NCP in the cryoEM density was essentially the same as described for other alphaviruses (Zhang et al., 2002). Furthermore, the slight asymmetric twist of the pentameric and hexameric rings about the 5 and 2-fold axes, respectively (Figures 1A,2B) was the same as that observed for Ross River virus (Cheng et al., 1995).

The lack of a substantial rigid association between the two lobes of the NCP's chymotrypsin-like structure suggests that there could be conformational adjustments between the two lobes. Consequently, the NCP structure was split into two domains (amino acid residues 119 to 183 and 184 to 267) representing the two lobes on either side of the substrate-binding cleft. These domains were then fitted independently at each of the four quasi-equivalent positions. The sumf values (Table 1B) for fitting the NCP into the pentameric capsomer density (position #4, Figure 3B) were slightly greater than for the NCPs in the hexameric capsomer (positions #1, #2, and #3), as was the case for the E1 protein fitting near the 5-fold axes. The correlation of the average residue densities between the four quasi-equivalent NCP positions was roughly the same as for fitting of the glycoproteins (Table 5), indicating the same level of amino acid recognition for the glycoproteins and the NCPs.

Contrary to the organization of the E1-E2 glycoproteins, where there are extensive inter-heterodimer contacts within the trimeric spikes, there were few contacts among the three NCP molecules around the icosahedral 3-fold axes and no contacts among the three NCP molecules around the quasi-3-fold axes. However, the six NCP molecules that form the hexameric capsomers around the icosahedral 2-fold axes made extensive contacts (Figure 3B and Table 6). These interactions were between the loop formed by residues 172 to 183 in one NCP molecule with residues in loops 190 to 206, 239 to 243 and 261 to 262 in a neighboring NCP molecule. Moreover, salt bridges were formed between Lys178 and Glu240. The NCP molecules around the icosahedral 5-fold axes made similar contacts and salt bridges, as expected by their quasi-similar environments.

It was unexpected that the positions of the quasi-2- and 3-fold symmetry elements were matched almost perfectly between the external organization of the glycoproteins into 80 trimers and the internal organization of the NCP into 12 pentamers and 30 hexamers (Figure 2). The nucleoplasmic cores assemble in the cytoplasm (Weiss et al., 1989; Owen and Kuhn, 1996; Soonsawad et al., 2010) and are transported to the plasma membrane where they interact with the endodomain of the E2 protein. This provides an opportunity for the C-terminal regions of E2 to bind to the hydrophobic pocket on the surface of the nucleocapsid (Lee et al., 1996), thus nucleating the E1-E2 glycoproteins to form an icosahedral shell in synchrony with the icosahedral symmetry of the cores.

Model building of the E1-E2 TM region and endodomains

The crystal structure of CHIKV E1-E2 (Voss et al., 2010) consists of residues 1-393 of E1 and 7-342 of E2 but does not include the TM regions and endodomains. Cα models were built of the missing part of the E1-E2 ectodomains, the TM helices, and the E2 endodomain, based on Rosetta structural predictions (Lange and Baker, 2012) and on the cryoEM density. The alignment of the amino acid sequence with the cryoEM density was based on assigning the kink observed in the TM helix of E1 to the likely flexible Thr-Gly-Gly sequence. In addition, the alignment was based on placing the interaction of Tyr 400 into the hydrophobic pocket around Trp 251 in the NCP, corresponding to SINV NCP residue Trp 247 (Lee et al., 1996). The TM regions of E1 and E2 were built independently for each of the four quasi-equivalent regions of the cryoEM map (Figure 3C). The E1 and E2 helices were parallel and situated close together with some hydrophobic contacts between them (Table 6). The T = 4 quasi-symmetry operators that relate the pairs of E1-E2 TM helices to other E1-E2 TM helical pairs were in similar positions to the symmetry operators that relate E1-E2 ectodomains to each other.

Although the lipid membrane is presumably flexible and likely to have different organization of lipid molecules in different particles, the cryoEM density showed a preference for roughly radial density features. Presumably, these densities represent an average of similarly placed lipid molecules, as was observed in the 4 Å resolution X-ray crystallographic map of the dsDNA PRD1 bacteriophage (Cockburn et al., 2004).

Comparison with the VEEV structure

The only other high resolution structure of an alpha virus is that of VEEV (Zhang et al., 2011). One way of comparing the structure of VEEV with the structure of CHIK VLPs is by superimposing their icosahedral symmetry elements. This shows that the overall radius of these two viruses is essentially the same. However, the positions and orientations of the molecular protein components have changed by up to 9.7 Å and 11.7°, respectively (Table 8), which is far greater than the RMS Cα atom differences between quasi equivalent positions of the component molecular components in the CHIK VLPs (maximum RMS distance is 2.4 Å, Table 4). It is also much larger than the RMS difference between Cα atoms on superimposing the equivalent molecular components of VEEV and CHIK VLPs without regard for the position and orientation of the icosahedral symmetry axes (Table 8). Thus the change in molecular positions and orientation with respect to the icosahedral framework is highly significant. Although the capsid protein amino acids are more conserved than any of the other structural proteins (Table 8), by far the largest differences in position and orientation of the molecular components occur for the capsid proteins (Table 8). The greater conservation of the amino acid sequences for the capsid protein is also reflected in a lower RMS differences between Cα atoms between superimposed capsid structures of VEEV and CHIK VLPs (less than 1.2 Å) as opposed to a comparison of the other molecular components (greater than 2.0 Å, Table 8).

Table 8

Comparison between VEEV and CHIK VLPs.

https://doi.org/10.7554/eLife.00435.013
E1-E2TMCapsid
A1A2A3A4A1A2A3A4A1A2A3A4
d in Å2.63.53.52.24.82.32.13.09.67.79.09.7
κ in degrees2.33.63.73.45.48.210.410.09.78.010.011.7
dCα in Å2.72.52.42.12.62.22.52.01.11.21.21.2
num69668070971411510099101149149149149
% identityE1 = 56%, E2 = 34%TM = 24%Capsid = 62%

  1. Differences between the positions of the centers of mass (d) and orientation (κ) relative to the icosahedral symmetry axes of the (E1-E2) heterodiners, the trans membrane (TM) helices and capsid proteins at positions A1, A2, A3 and A4 (see Figure 2). The RMS distances (dCα) are given between the number (num) of equivalent Cα positions of the superimposed molecules A1 to A4. The percentage of identical amino acids in the aligned proteins is given on the last line of the table

The above analysis shows that these two alpha-viruses maintain a fairly constant tertiary structure, notwithstanding moderate amino acid sequence variations. However, the quaternary organization has greater variability, although still retaining the same T = 4 quasi symmetry. Considering that the nucleocapsid core is assembled first and the subsequent association of the glycoproteins is based on the pre-existing cores, it is surprising that the difference of structure between the CHIKV and VEEV cores is less evident in the virus' ectodomains.

CryoEM structures of the CHIK VLPs in complex with neutralizing antibody Fab fragments

Four cryoEM structures were determined to ∼15 Å resolution of CHIK VLPs complexed with the Fab fragments of the neutralizing mouse MAbs CHK-9, CHK-152, m242 and m10 (Figure 4). All Fab fragments bound the VLPs with one Fab fragment per E2 molecule. Except for CHK-152, the concentration of the Fab molecules required for 50% neutralization was 100 or more times greater than the intact IgG (Figure 5). This implies a fundamental difference between how CHK-152 binds to the VLPs compared to the other Fabs. Because CHK-152 MAb and Fab fragment neutralize CHIKV at almost the same low concentration, it is unlikely that this MAb requires bivalent binding to inhibit infection.

Figure 4
Download asset Open asset
CHIKV VLP Fab complexes.

Right: “Difference maps” showing the surface structure of the Fab molecules. The difference maps were produced by setting all density within 1.7 Å of an atom in the VLP structure to zero. Left: “Road maps” showing the projected surfaces of the VLP-Fab complexes for one (triangular) icosahedral asymmetric unit. The coloring is according to the radial distance of the surface from the center of the VLP. The footprints of the Fabs are shown in yellow. The radial projection of the whole Fab molecules onto the surface of the VLP is shown with white contours representing the height of the projected density. The mosaic background shows the amino acids that form the viral surface.

https://doi.org/10.7554/eLife.00435.014
Figure 5
Download asset Open asset
Inhibition of CHIKV by IgG and Fab fragments.

Neutralizing activity of intact MAb and Fab fragments against CHIKV is shown for (A) CHK-9, (B) CHK-152, (C) m10 and (D) m242, as determined by the reduction in the number of focus-forming units (FFU) on Vero cells. MAbs and Fab fragments were mixed with 100 FFU of infectious CHIKV (strain La Reunion 2006 OPY-1) for one hour at 37°C before infecting Vero cells. Each data point is the average of three independent experiments performed in triplicate. Error bars represent standard deviations.

https://doi.org/10.7554/eLife.00435.015

The E1-E2 coordinates from the high-resolution VLP map were placed into the Fab bound cryoEM maps by aligning the icosahedral symmetry axes. The coordinates of the E2 B domain in the VLPs did not agree well with the B domain density in the maps of the VLP complexed with the CHK-9 or m242 Fab fragments, suggesting that the E2 B domain had a different orientation in these Fab-bound structures than in the non-complexed VLPs. To determine the modified position of the E2 B domain in these two complexes, the cryoEM density distributions were modified by setting to zero all density that was within 3 Å of all atoms in E1 and of all atoms in the A, C and D domains of E2. The crystallographic coordinates of the E2-B domain (Voss et al., 2010) were then fitted into the modified cryoEM density distributions, assuming T4 quasi-icosahedral symmetry.

The crystal structures of CHK-9 and m242 Fab fragments were determined to 3.0 Å and 3.1 Å resolution, respectively (Table 9). These Fab structures were fitted into the cryoEM maps of the corresponding VLP-Fab complexes. For the VLP-CHK-152 and VLP-m10 Fab complexes, for which no crystal structures of the corresponding Fab fragments had been determined, both the CHK-9 and m242 Fab crystal structures were tested to determine which Fab yielded the better fit to the cryoEM map of the complex. VLP-CHK-152 and VLP-m10 Fab cryoEM maps were best fitted with the Fab structures of CHK-9 and of m242, respectively (Figure 6). The T = 4 symmetry parameters, determined by fitting the E1-E2 heterodimer into the higher resolution VLP density, were used for fitting the Fab fragments into the cryoEM densities of the appropriate VLP-Fab complexes. The height of the density for the variable domains of the Fab fragments was comparable to that of the CHIKV glycoprotein, suggesting a nearly 100% binding occupancy of the Fab fragments under saturating binding conditions. The density for the constant domains (remote from the VLPs) of the Fab fragments was about 17% lower, consistent with the flexibility of the elbow angle between the constant and variable domains. The footprints of the Fab fragments on the surface of the VLPs were determined with the RIVEM program (Xiao and Rossmann, 2007) (Figures 4,7).

Table 9

Data collection and refinement statistics for the m242 and CHK-9 Fab molecules

https://doi.org/10.7554/eLife.00435.016
Fab m242Fab CHK9
Data collection
X-ray source23-ID-B23-ID-B
Wavelength (Å)1.031.03
Resolution (Å)3.13
Space groupP21212C2
Unit cell (Å)a = 137.0, b = 89.1, c = 94.4a = 87.9, b = 57.7, c = 118.1
Unique reflections21,11535,876
Average redundancy3.75.8
I/σ*17.5 (4.9)25.0 (5.0)
Completeness (%)98.3 (96.4)99.6 (99.9)
Rmerge (%)†10.1 (23.4)5.5 (39.2)
Refinement
Resolution (Å)3.11.8
Rworking (%)‡29.3 (34.5)18.7 (23.3)
Rfree (%)§33.1 (34.3)22.0 (27.7)
rmsd bonds (Å)0.0050.01
rmsd angels (°)0.850.91
# of residues Ramachandran disallowed01
  1. *

    Values in parentheses throughout the table correspond to the outermost resolution shell

  2. †

    Rmerge = Σ| I - <I> | / Σ I, where I is the measured intensity of reflections

  3. ‡

    Rworking = Σ||Fobs| -|Fcalc|| / Σ|Fobs|

  4. §

    Rfree has the same formula as Rworking except that calculation was made with the structure factors from the test set

Figure 6
Download asset Open asset
Fit of the Fab structures (blue) into the cryoEM density “difference” maps (hatched grey surface) calculated as described for Figure 4.

Shown also is how the Fab molecules bind to the E1(red)-E2 (green) heterodimer.

https://doi.org/10.7554/eLife.00435.017
Figure 7
Download asset Open asset
“Road maps” showing footprint of four neutralizing Fabs on the VLP surface at position #3 as defined in Figure 2.

In order to differentiate between amino acids in different quasi 3-fold related subunits, their identity is defined as the amino acid sequence number in E1 + 2000, 3000, and 4000, and in E2 + 2500, 3500, and 4500 for positions #2, 3, and 4, respectively (see Figure 2). The surface is colored according to radial distance from the center of the VLP. The A, B, and D domains of E2 are bounded by a black, white and dashed white line, respectively. Residues in the putative receptor binding site on domain A of E2 are bounded by a yellow dashed line. The footprint of the Fabs onto the VLP surface is outlined in yellow.

https://doi.org/10.7554/eLife.00435.018

The footprint of CHK-9 and m242 Fabs onto the VLP localized primarily onto the domain A and did not include any part of domain B. Indeed, binding of these fragments to the VLPs had nudged domain B ∼10 Å sideways, out of the way of the Fab fragment (Figure 7). In contrast, the footprint of the m10 Fab mapped exclusively on domain B and did not cause any observable movement (Figure 7). The footprint of CHK-152, however, spanned both domains B and A as well as the linker peptides to domain B represented by domain D. Thus, unlike the other Fabs, CHK-152 has the capacity to cross-link domain B with domain A. Alphaviruses become fusogenic at low pH when domain B, held in place by the flexible β-ribbons of domain D, becomes mobile and exposes the fusion loop on domain II of E1. As CHK-152 efficiently blocks viral fusion with host membranes (P Pal and M Diamond, unpublished observations), it may do so by cross-linking the flexible domain B to the rest of the virus surface. An analogous mechanism of inhibition was described for neutralization of West Nile virus, in which the viral envelope proteins were cross-linked by the bound FAB fragments (Kaufmann et al., 2010).

All four MAbs examined here bind to the end of the trimeric spikes that includes the putative receptor attachment region (Davis et al., 1991; Dubuisson and Rice, 1993; Kinney et al., 1993; Klimstra et al., 1998; Bernard et al., 2000; Lee et al., 2002; Wang et al., 2003; Pierro et al., 2007, 2008; Ryman et al., 2007; Li et al., 2010), although only the CHK-9 Fab footprint actually overlaps the site (Figure 7). Nevertheless the proximity of these footprints to the possible receptor-binding site could probably allow neutralization by competitively limiting access of the virus to the cell surface receptor. To block formation of a fusogenic E1 trimer it is necessary to bind a Fab molecule to only one of the three E2 molecules in every spike. However, to block receptor attachment it is necessary to block all three receptor-binding sites on every spike. Such a stoichiometry would require substantially higher concentrations of Fabs to inhibit infectivity, as seen with CHK-9, m10 and m242 Fabs (Figure 5). The difference between MAb and Fab neutralization efficiency is presumably the result of bivalent attachment and cross-linking on a single virion or aggregation between virions. In summary our structural studies suggest that the CHK-9, m10, and m242 Fab fragments neutralize CHIKV infectivity by blocking the cellular binding site on the A domain of E2 and that the CHK152 inhibits fusion by stabilizing domain B of E2, preventing exposure of the fusion loop on E1.

Materials and methods

CHIK VLP production and purification

Request a detailed protocol

CHIK VLPs have been reported as structural and immunologically indistinguishable from the mature infectious virus (Akahata et al., 2010); the use of CHIKV VLPs instead of live infectious virus has allowed us to avoid performing all structural studies under biosafety level 3 conditions. CHIK VLPs were produced and purified as described previously (Akahata et al., 2010). Following purification, the buffer was exchanged to PBS and the VLPs were concentrated to 1 mg/ml.

Antibody production and purification of Fab fragments

Request a detailed protocol

MAbs m242, CHK-9, CHK-152, and m10 were purified from hybridoma superantants by protein A Sepahrose chromatography, and then digested with papain. The resultant Fab fragments were purified by sequential protein A Sepahrose and Superdex 75 16/60 size exclusion chromatography.

Neutralization assays

Request a detailed protocol

Serial dilutions of MAbs or their Fab fragments were incubated with 100 focus-forming units (FFU) of CHIKV (CHIKV La Reunion 2006 OPY-1) for one hour at 37°C. MAb or Fab-virus complexes were added to Vero cells in 96-well plates. After 90 min, cells were overlaid with 1% (w/v) methylcellulose in Modified Eagle Media supplemented with 4% FBS. Plates were harvested 18 hr later, and fixed with 1% PFA in PBS. The plates were incubated sequentially with 500 ng/ml of mouse-human chimeric version of CHK-9 (containing a human γ1 constant region, P Pal and M Diamond, unpublished results) and horseradish peroxidase conjugated goat anti-human IgG in PBS supplemented with 0.1% saponin and 0.1% BSA. CHIKV-infected foci were visualized using TrueBlue peroxidase substrate (KPL), quantitated on an ImmunoSpot 5.0.37 macroanalyzer (Cellular Technologies Ltd), and analyzed using GraphPad Prism software.

CryoEM imaging, processing and three-dimensional reconstruction of CHIK VLPs

Request a detailed protocol

Three microliter aliquots of purified CHIK VLP sample (1 mg/ml) were applied to 400 mesh C-flat grids (1.2 μm hole size) and double blotted inside a 100% humidified FEI vitrobot chamber, using a blotting time of 2 s. The grids were frozen by plunging into liquid ethane. A FEI Titan Krios electron microscope operated at 300 kV was used to collect cryoEM images on Kodak SO-163 films at a magnification of 59,000× and an electron dosage of ∼25 e−/Å2 at the boiling point of liquid nitrogen. The microscope beam was aligned parallel with the optical axis of the microscope using the coma-free alignment technique (Kimoto et al., 2003). The quality of the alignment was checked by taking a CCD image of the carbon film under the same dosage as was used for data collection. The alignment was accepted when the Thon rings on the carbon support were visible beyond (1/4.5) Å−1. A total of 1532 films were developed using full strength D19 (Kodak) solution. Micrographs were digitized with a Nikon Coolscan 9000ED scanner at a step size of 6.35 μm/pixel, which yielded a pixel size of 1.05 Å/pixel on the specimen. The absolute pixel size was determined by varying the pixel size from 1.05 to 1.14 in steps of 0.01 Å to find the “best” fit of the CHIKV E1E2 heterodimer crystal structure into the final cryoEM density. The fit was based on the average density of all the non-hydrogen fitted atoms (sumf). The highest value of sumf was found when the pixel separation was 1.11 Å.

Of the 1532 micrographs, 1012 did not have drift or charging problems and had Thon rings beyond (1/6.0) Å−1 observed through the sample. These were selected for further processing. A total of 52,183 particles were boxed with the EMAN2 program e2boxer (Tang et al., 2007). The selected micrographs were found to be under-focused by between 1.1 μm and 2.9 μm using the EMAN program CTFIT (Ludtke et al., 1999). Initially a “four-fold binned” map was used in which each set of 4 × 4 pixels was replaced by one pixel with the average height of the original sixteen pixels. Thus, the pixels in this map were separated by 4.44 Å, limiting the resolution to ∼13 Å. An initial CHIK VLP model was produced by selecting particles with 5- 3- and 2-fold symmetric projections from which an icosahedral reconstruction was made. This model was iteratively refined using a coarse 5° angular step size to improve the orientation and positioning of each particle until convergence was reached, as indicated by the lack of change of the Fourier shell correlation between successive cycles. At this point, the model had the characteristics of an alphavirus (Cheng et al., 1995). In further iterations, the model was modified progressively by setting to zero both the inner density corresponding to the RNA genome and also the outer noise background density. Orientation and centers for each of the particles were determined by comparing the observed two-dimensional projections with projections of the model generated using a 1° angular step size. Each observed projection was allocated to the three best fits in the calculated projection classes. The images associated with each class were averaged and used to generate an improved model for the next cycle. Only those particles whose orientations and centers had changed their orientation and position by less than 2° and 2 pixels between successive cycles were retained for the next reconstruction. Unstable particles were discarded after each cycle and were not reconsidered.

When the global searches had converged, local fine angular and positional refinements were conducted using 2-fold instead of 4-fold binned images. Particles were discarded when their three top solutions did not converge. After each cycle, the resolution of the map was estimated and used to modify the model by filtering out the Fourier coefficients beyond the current resolution. Further improvement was hindered because it was assumed that the magnification of each micrograph was the same and that the defocus distance of each particle on a given micrograph was the same. The program Frealign (Grigorieff, 2007) was used to refine orientation, position, defocus distance and relative pixel size for each micrograph using the original scanned data without binning. Particles were rejected that were found after refinement to have changes in any of these parameters by more than 2σ. Overall convergence was established when the plot of the Fourier shell correlation vs resolution did not improve as estimated by visual inspection. The final map was reconstructed using 36,236 particles. The resolution of this map was estimated to be 4.6 Å, based on the spacing at which the Fourier shell correlation fell below 0.5 (Van Heel, 1987). To verify the effective resolution of the map (Chang et al., 2012), the orientation of each image was determined with respect to the final map after having been low passed filtered to 7 Å. The resultant orientations of the images were used to calculate the Fourier shell correlation between two equal random groups of all the images. The map resolution was found to be 5.3 Å using the 0.5 Fourier shell correlation criterion. A series of maps were calculated using a B factor applied to the Fourier coefficients of the inverted map to estimate the noise level and backbone density continuity. The best value of B was -40 Å2, which slightly increased the emphasis on the high order terms.

Fitting of CHIKV crystal structure into the EM map

Request a detailed protocol

Initially, the E1-E2 heterodimer crystal structure was fitted as a rigid body close to an icosahedral 3-fold axis (Figure 2) defined by the polar angles psi = 69.1°, fi = -90° (Rossmann and Blow, 1962). This position was called “#1” (Figure 2A). The icosahedral 3-fold axis was used to generate the complete icosahedral “i3” spike. Next, the i3 spike was rotated by the quasi-2-fold axis at psi = 74.6°, fi = -80.9° (an orientation appropriate for T = 4 quasi-symmetry) to produce the “q3” spike, generating positions #2, #3, #4 of the heterodimer (Figure 2A). The position of this quasi-2-fold axis was refined to maximize the average height of the fitted density (sumf) while minimizing the number of clashes between symmetry related E1-E2 heterodimers (Rossmann et al., 2001).

The resultant positions of the four quasi-equivalent CHIKV heterodimers were used as starting positions for further position and orientation refinement of individual domains. E1 was divided into domain I (residues 1-36, 132-168, 273-293), II (residues 37-131, 169-272) and III (residues 294-393). E2 was divided into domain A (residues 16-134), B (residues 173-231), C (residues 269-342), and D the β-ribbon connector (residues 7-15, 135-172, 232-268). These domains were fitted sequentially (I, II, III, A, C, B, and D). After each domain was fitted independently, the map was modified by zeroing out the density around each fitted atom within a radius of 3.0 Å. Thus, when fitting the next domain, atoms placed into zero density acted as a restraint by reducing the overall fitting criterion, “rcrit” (Rossmann et al., 2001). This process substituted for avoiding clashes had all domains been fitted simultaneously. After the individual domains had been fitted, the bond geometry between the carboxy end of one domain and the amino end of the next domain were regularized.

The average density height of the final atomic structure (sumf) was better for the individual domain fitting results than for the original T = 4 rigid body fitting results (Table 1A), validating the fitting strategy. The quality of the independently fitted domain structure was also validated by comparing the fit of the atomic structures at each of the four quasi-equivalent densities. The density of each non-hydrogen atom in the CHIKV model was determined by interpolation using the densities at the eight surrounding grid points. The average density (ρi) of residue i was then set to the average densities of all the atoms associated with this residue. The correlation coefficient (CC) was calculated between the densities ρi(A) and ρi(B) where these are the cryoEM densities of residue i at the quasi-equivalent positions A and B. The sums are taken over i = 1, N where N is the number of residues being matched.

CC=∑[(ρi(A)-<ρi(A)>)(ρi(B)-<ρi(B)>)][∑(ρi(A)-<ρi(A)>)2∑(ρi(B)-<ρi(B)>)2]½

Identifying non-covalent contacts between molecules that form the CHIKV VLPs

Request a detailed protocol

Non-covalent contacts between molecular entities were determined and classified by identifying the number of times a specific atom in one molecule was within a defined distance, Dlim = 3.5 Å, of any atom in a specific neighboring molecule. The contacts of each atom in a given molecule were added to yield a measure of how many close contacts the first molecule had with the second molecule. This count is roughly equivalent to measuring the surface contact area between the two molecules. In addition to measuring the total number of contacts, a count was made of that portion of contacts that represented hydrophobic interactions by counting only carbon-carbon atom distances less than Dlim. Also, the number of potential hydrogen bonds was estimated by counting the number of short (Dlim < 3.5 Å) distances between oxygen and nitrogen atoms. Lastly, a count was made of possible formation of salt bridges by counting the number of occasions that an Arg, Lys or His residue approached to within Dlim of an Asp or Glu residue.

Radius of gyration of the trimeric spikes

Request a detailed protocol

The radius of gyration of the spikes was computed as:

rg=∑(ri2)N,

where the Sum (i = 1 to N) is taken over the N Cα atoms in E1 and E2 of one spike, and ri is the radius of the ith atom measured from the spike axis. The latter is defined by the positions of the centers between the 3-fold (in the i3 spike), or quasi-3-fold (in the q3 spike) related atoms.

Comparison of VEEV and CHIK VLPs

Request a detailed protocol

Alignment of the molecular components between VEEV and CHIK VLPs was performed with the HOMOlogy program (Rossmann and Argos, 1975). The position and orientation of each molecular component was calculated with respect to the equivalenced Cα atoms.

CryoEM imaging, processing and three-dimensional reconstruction of the Fab-CHIK VLP complexes

Request a detailed protocol

CHIK VLP particles were incubated with Fab fragments at 4°C for two hours using a stoichiometric ratio of about four Fab fragments per E2 molecule. Samples were hand-blotted, flash-frozen on holey grids in liquid ethane. Images were recorded at 47,000× magnification with a CM300 field emission gun microscope using electron dose levels of approximately 20 electrons per Å2. All micrographs were digitized at 6.35 μm per pixel using a Nikon scanner. Individual particle images were boxed using the program e2boxer in the EMAN2 package. Subsequently, the boxed images were two-fold binned, resulting in a sampling of the specimen images at 2.69 Å intervals. The CTFIT program from the EMAN package was used to determine the contrast transfer function parameters. A cryoEM reconstruction of CHIK VLP, low pass filtered to 18 Å, was used as an initial model for orientation determination and further refinement. The number of particles incorporated into the final reconstruction was 1728, 1820, and 1599 for Fabs m242, CHK9 and m10, giving final resolutions of 15.6 Å, 14.9 Å and 16.9 Å on the basis of a 0.5 Fourier shell correlation threshold, respectively.

Crystallization and structure determination of the m242 and CHK-9 Fab fragments

Request a detailed protocol

The m242 Fab was crystallized in 2M ammonium sulfate and 0.1 M sodium acetate pH 4.5 and the CHK-9 Fab was crystallized in 25% PEG 3350, 0.1M Tris-Cl pH 8.5 and 0.2M lithium sulfate. Crystals were grown by vapor diffusion in hanging drops at 20°C. Crystals were flash-frozen and data were collected at 100K at the Advanced Photon Source (APS) GM/CA-CAT 23-ID-B beamline. Data were processed with the HKL2000 program (Otwinowski and Minor, 1997). The m242 Fab crystals diffracted to 3.1 Å resolution and had a P21212 space group with cell dimensions a = 137.0, b = 89.1 and c = 94.4 Å. There were two Fab molecules per asymmetric unit. The CHK-9 crystals diffracted to 3.0 Å resolution and had a C2 space group with cell dimensions a = 87.9, b= 57.7 and c = 118.1 Å. There was one Fab molecule per asymmetric unit. The structures of the m242 Fab and CHK-9 Fab crystals were determined by molecular replacement using a Fab crystal structure (Protein Data Bank: 3DGG) as a search model with the program Phaser (Mccoy et al., 2007). The structures were refined using the program Refmac (Murshudov et al., 1997) (Table 9). The rms difference between the Cα atoms of the two molecules in the asymmetric unit of the m242 Fab crystals was 0.7 Å.

Data availability

The following data sets were generated

References

    1. Brünger AT
    (1992)
    Free R value: A novel statistical quantity for assessing the accuracy of crystal structures
    Nature 355:472–475.
    1. Carleton M
    2. Lee H
    3. Mulvey M
    4. Brown DT
    (1997)
    Role of glycoprotein PE2 in formation and maturation of the Sindbis virus spike
    J Virol 71:1558–1566.
    1. Chapman MS
    (1995)
    Restrained real-space macromolecular atomic refinement using a new resolution-dependent electron-density function
    Acta Crystallogr A 51:69–80.
    1. Cheng RH
    2. Kuhn RJ
    3. Olson NH
    4. Rossmann MG
    5. Choi HK
    6. Smith TJ
    et al. (1995)
    Nucleocapsid and glycoprotein organization in an enveloped virus
    Cell 80:621–630.
    1. Choi HK
    2. Lu G
    3. Lee S
    4. Wengler G
    5. Rossmann MG
    (1997)
    Structure of Semliki Forest virus core protein
    Proteins 27:345–359.
    1. Davis NL
    2. Powell N
    3. Greenwald GF
    4. Willis LV
    5. Johnson BJ
    6. Smith JF
    et al. (1991)
    Attenuating mutations in the E2 glycoprotein gene of Venezuelan equine encephalitis virus: construction of single and multiple mutants in a full-length cDNA clone
    Virology 183:20–31.
    1. Dubuisson J
    2. Rice CM
    (1993)
    Sindbis virus attachment: isolation and characterization of mutants with impaired binding to vertebrate cells
    J Virol 67:3363–3374.
    1. Gaedigk-Nitschko K
    2. Schlesinger MJ
    (1990)
    The Sindbis virus 6K protein can be detected in virions and is acylated with fatty acids
    Virology 175:274–281.
    1. Gaedigk-Nitschko K
    2. Schlesinger MJ
    (1991)
    Site-directed mutations in Sindbis virus E2 glycoprotein's cytoplasmic domain and the 6K protein lead to similar defects in virus assembly and budding
    Virology 183:206–214.
    1. Griffin DE
    (2007)
    Fields virology
    Fields virology, ., 5th ed., Philadelphia, Lippincott Williams & WIlkins.
    1. Griffin DE
    2. Ubol S
    3. Despres P
    4. Kimura T
    5. Byrnes A
    (2001)
    Role of antibodies in controlling alphavirus infection of neurons
    Curr Top Microbiol Immunol 260:191–200.
    1. Hahn CS
    2. Strauss JH
    (1990)
    Site-directed mutagenesis of the proposed catalytic amino acids of the Sindbis virus capsid protein autoprotease
    J Virol 64:3069–3073.
    1. Kinney RM
    2. Chang GJ
    3. Tsuchiya KR
    4. Sneider JM
    5. Roehrig JT
    6. Woodward TM
    et al. (1993)
    Attenuation of Venezuelan equine encephalitis virus strain TC-83 is encoded by the 5'-noncoding region and the E2 envelope glycoprotein
    J Virol 67:1269–1277.
    1. Klimstra WB
    2. Ryman KD
    3. Johnston RE
    (1998)
    Adaptation of Sindbis virus to BHK cells selects for use of heparan sulfate as an attachment receptor
    J Virol 72:7357–7366.
    1. Kuhn RJ
    (2007)
    Fields virology
    Fields virology, 5th ed., Philadelphia, Lippincott Williams & Wilkins.
    1. Lange OF
    2. Baker D
    (2012)
    Resolution-adapted recombination of structural features significantly improves sampling in restraint-guided structure calculation
    Proteins 80:884–895.
    1. Lee P
    2. Knight R
    3. Smit JM
    4. Wilschut J
    5. Griffin DE
    (2002)
    A single mutation in the E2 glycoprotein important for neurovirulence influences binding of sindbis virus to neuroblastoma cells
    J Virol 76:6302–6310.
    1. Lee S
    2. Owen KE
    3. Choi HK
    4. Lee H
    5. Lu G
    6. Wengler G
    et al. (1996)
    Identification of a protein binding site on the surface of the alphavirus nucleocapsid and its implication in virus assembly
    Structure 4:531–541.
    1. Lescar J
    2. Roussel A
    3. Wien MW
    4. Navaza J
    5. Fuller SD
    6. Wengler G
    et al. (2001)
    The Fusion glycoprotein shell of Semliki Forest virus: an icosahedral assembly primed for fusogenic activation at endosomal pH
    Cell 105:137–148.
    1. Linger BR
    2. Kunovska L
    3. Kuhn RJ
    4. Golden BL
    (2004)
    Sindbis virus nucleocapsid assembly: RNA folding promotes capsid protein dimerization
    RNA 10:128–138.
    1. Meyer WJ
    2. Johnston RE
    (1993)
    Structural rearrangement of infecting Sindbis virions at the cell surface: mapping of newly accessible epitopes
    J Virol 67:5117–5125.
    1. Mulvey M
    2. Brown DT
    (1995)
    Involvement of the molecular chaperone BiP in maturation of Sindbis virus envelope glycoproteins
    J Virol 69:1621–1627.
    1. Otwinowski Z
    2. Minor W
    (1997)
    Methods in Enzymology
    Methods in Enzymology, New York, Academic Press.
    1. Owen KE
    2. Kuhn RJ
    (1996)
    Identification of a region in the Sindbis virus nucleocapsid protein that is involved in specificity of RNA encapsidation
    J Virol 70:2757–2763.
    1. Paredes AM
    2. Brown DT
    3. Rothnagel R
    4. Chiu W
    5. Schoepp RJ
    6. Johnston RE
    et al. (1993)
    Three-dimensional structure of a membrane-containing virus
    Proc Natl Acad Sci USA 90:9095–9099.
    1. Pletnev SV
    2. Zhang W
    3. Mukhopadhyay S
    4. Fisher BR
    5. Hernandez R
    6. Brown DT
    et al. (2001)
    Locations of carbohydrate sites on alphavirus glycoproteins show that E1 forms an icosahedral scaffold
    Cell 105:127–136.
    1. Robinson MC
    (1955)
    An epidemic of virus disease in Southern Province, Tanganyika territory, in 1952-1953
    Trans R Soc Trop Med Hyg 49:28–32.
    1. Rossmann MG
    2. Argos P
    (1975)
    A comparison of the heme binding pocket in globins and cytochrome b5
    J Biol Chem 250:7525–7532.
    1. Rossmann MG
    2. Blow DM
    (1962)
    The detection of sub-units within the crystallographic asymmetric unit
    Acta Crystallogr 15:24–31.
    1. Simizu B
    2. Yamamoto K
    3. Hashimoto K
    4. Ogata T
    (1984)
    Structural proteins of Chikungunya virus
    J Virol 51:254–258.
    1. Skoging U
    2. Vihinen M
    3. Nilsson L
    4. Liljestrom P
    (1996)
    Aromatic interactions define the binding of the alphavirus spike to its nucleocapsid
    Structure 4:519–529.
    1. Smith TJ
    2. Cheng RH
    3. Olson NH
    4. Peterson P
    5. Chase E
    6. Kuhn RJ
    et al. (1995)
    Putative receptor binding sites on alphaviruses as visualized by cryoelectron microscopy
    Proc Natl Acad Sci USA 92:10648–10652.
    1. Tellinghuisen TL
    2. Hamburger AE
    3. Fisher BR
    4. Ostendorp R
    5. Kuhn RJ
    (1999)
    In vitro assembly of alphavirus cores by using nucleocapsid protein expressed in Escherichia coli
    J Virol 73:5309–5319.
    1. Tellinghuisen TL
    2. Kuhn RJ
    (2000)
    Nucleic acid-dependent cross-linking of the nucleocapsid protein of Sindbis virus
    J Virol 74:4302–4309.
    1. Tong L
    2. Choi HK
    3. Minor W
    4. Rossmann MG
    (1992)
    The structure determination of Sindbis virus core protein using isomorphous replacement and molecular replacement averaging between two crystal forms
    Acta Crystallogr A 48:430–442.
    1. Van Heel M
    (1987)
    Similarity measures between images
    Ultramicroscopy 21:95–100.
    1. Wahlberg JM
    2. Garoff H
    (1992)
    Membrane fusion process of Semliki Forest virus. I: low pH-induced rearrangement in spike protein quaternary structure precedes virus penetration into cells
    J Cell Biol 116:339–348.
    1. Wang E
    2. Brault AC
    3. Powers AM
    4. Kang W
    5. Weaver SC
    (2003)
    Glycosaminoglycan binding properties of natural venezuelan equine encephalitis virus isolates
    J Virol 77:1204–1210.
    1. Weiss B
    2. Nitschko H
    3. Ghattas I
    4. Wright R
    5. Schlesinger S
    (1989)
    Evidence for specificity in the encapsidation of Sindbis virus RNAs
    J Virol 63:5310–5318.
    1. Zhang W
    2. Mukhopadhyay S
    3. Pletnev SV
    4. Baker TS
    5. Kuhn RJ
    6. Rossmann MG
    (2002)
    Placement of the structural proteins in Sindbis virus
    J Virol 76:11645–11658.
    1. Zhu J
    2. Cheng L
    3. Fang Q
    4. Zhou ZH
    5. Honig B
    (2010)
    Building and refining protein models within cryo-electron microscopy density maps based on homology modeling and multiscale structure refinement
    J Mol Biol 397:835–851.

Article and author information

Author details

  1. Siyang Sun

    Department of Biological Sciences, Purdue University, West Lafayette, United States
    Present address
    8211 Scicor Drive, Indianapolis, United States
    Contribution
    SS, Collected cryoEM data and produced cryoEM reconstructions, Determined the crystal structures of the Fab fragments and performed the model fitting, Wrote and critiqued the paper
    Contributed equally with
    Ye Xiang
    Competing interests
    The authors declare that no competing interests exist.

  2. Ye Xiang

    Department of Biological Sciences, Purdue University, West Lafayette, United States
    Present address
    School of Medicine, Tsinghua University, Beijing, China
    Contribution
    YX, Collected cryoEM data and produced cryoEM reconstructions, Wrote and critiqued the paper
    Contributed equally with
    Siyang Sun
    Competing interests
    The authors declare that no competing interests exist.

  3. Wataru Akahata

    Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States
    Contribution
    WA, Produced the CHIK VLPs as well as the m10 and m242 MAbs and Fab fragments
    Competing interests
    The authors declare that no competing interests exist.

  4. Heather Holdaway

    Department of Biological Sciences, Purdue University, West Lafayette, United States
    Present address
    Cleveland Center for Membrane and Structural Biology, Case Western Reserve University, Cleveland, United States
    Contribution
    HH, Helped with cryoEM data collection
    Competing interests
    The authors declare that no competing interests exist.

  5. Pankaj Pal

    Departments of Medicine, Molecular Microbiology, Pathology & Immunology, Washington University School of Medicine, St Louis, United States
    Contribution
    PP, Produced the CHK-9 and CHK-152 MAbs and Fab fragments, and performed all neutralization assays
    Competing interests
    The authors declare that no competing interests exist.

  6. Xinzheng Zhang

    Department of Biological Sciences, Purdue University, West Lafayette, United States
    Contribution
    XZ, Re-determined the resolution of the CHIK VLP cryoEM map and helped in making additional figures suggested by the reviewers
    Competing interests
    The authors declare that no competing interests exist.

  7. Michael S Diamond

    Departments of Medicine, Molecular Microbiology, Pathology & Immunology, Washington University School of Medicine, St Louis, United States
    Contribution
    MSD, Produced the CHK-9 and CHK-152 MAbs and Fab fragments and performed all neutralization assays, Wrote and critiqued the paper
    Competing interests
    The authors declare that no competing interests exist.

  8. Gary J Nabel

    Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, United States
    Contribution
    GJN, Produced the CHIK VLPs as well as the m10 and m242 MAbs and Fab fragments, Wrote and critiqued the paper
    Competing interests
    The authors declare that no competing interests exist.

  9. Michael G Rossmann

    Department of Biological Sciences, Purdue University, West Lafayette, United States
    Contribution
    MGR, Analyzed the fitted structures, wrote and critiqued the paper
    For correspondence
    mr@purdue.edu
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Institutes of Health (AI095366)

  • Siyang Sun
  • Ye Xiang
  • Heather Holdaway
  • Xinzheng Zhang
  • Michael G Rossmann

National Institutes of Health (AI104545)

  • Pankaj Pal
  • Michael S Diamond

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank the referees for pointing out that “a lesson to be drawn from this work is the merit of obtaining a cryoEM structure even at a lower resolution than the crystal structure to reveal the structure of the virus proteins free of potential crystallization artifacts and close to the solution state conformation”.

We thank Syd Johnson (MacroGenics) for the chimeric CHK-9 Mab, the sequence of the CHK-9 Fab fragment, and purifying CHK-9 and CHK-152 IgG. We are grateful to Josh Yoder for preparing the m242 Fab fragment. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. We thank Sheryl Kelly for formatting the manuscript.

Copyright

© 2013, Sun et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 6,427
    views
  • 807
    downloads
  • 131
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Siyang Sun
  2. Ye Xiang
  3. Wataru Akahata
  4. Heather Holdaway
  5. Pankaj Pal
  6. Xinzheng Zhang
  7. Michael S Diamond
  8. Gary J Nabel
  9. Michael G Rossmann
(2013)
Structural analyses at pseudo atomic resolution of Chikungunya virus and antibodies show mechanisms of neutralization
eLife 2:e00435.
https://doi.org/10.7554/eLife.00435

Share this article

https://doi.org/10.7554/eLife.00435

Categories and tags

Research organism

Further reading

    1. Structural Biology and Molecular Biophysics

    Ligand-coupled conformational changes in a cyclic nucleotide-gated ion channel revealed by time-resolved transition metal ion FRET

    Pierce Eggan, Sharona E Gordon, William N Zagotta
    Research Article

    Cyclic nucleotide-binding domain (CNBD) ion channels play crucial roles in cellular-signaling and excitability and are regulated by the direct binding of cyclic adenosine- or guanosine-monophosphate (cAMP, cGMP). However, the precise allosteric mechanism governing channel activation upon ligand binding, particularly the energetic changes within domains, remains poorly understood. The prokaryotic CNBD channel SthK offers a valuable model for investigating this allosteric mechanism. In this study, we investigated the conformational dynamics and energetics of the SthK C-terminal region using a combination of steady-state and time-resolved transition metal ion Förster resonance energy transfer (tmFRET) experiments. We engineered donor-acceptor pairs at specific sites within a SthK C-terminal fragment by incorporating a fluorescent noncanonical amino acid donor and metal ion acceptors. Measuring tmFRET with fluorescence lifetimes, we determined intramolecular distance distributions in the absence and presence of cAMP or cGMP. The probability distributions between conformational states without and with ligand were used to calculate the changes in free energy (ΔG) and differences in free energy change (ΔΔG) in the context of a simple four-state model. Our findings reveal that cAMP binding produces large structural changes, with a very favorable ΔΔG. In contrast to cAMP, cGMP behaved as a partial agonist and only weakly promoted the active state. Furthermore, we assessed the impact of protein oligomerization and ionic strength on the structure and energetics of the conformational states. This study demonstrates the effectiveness of time-resolved tmFRET in determining the conformational states and the ligand-dependent energetics of the SthK C-terminal region.

    1. Structural Biology and Molecular Biophysics

    A cryo-electron tomography study of ciliary rootlet organization

    Chris van Hoorn, Andrew P Carter
    Research Article

    Ciliary rootlets are striated bundles of filaments that connect the base of cilia to internal cellular structures. Rootlets are critical for the sensory and motile functions of cilia. However, the mechanisms underlying these functions remain unknown, in part due to a lack of structural information of rootlet organization. In this study, we obtain 3D reconstructions of membrane-associated and purified rootlets from mouse retina using cryo-electron tomography. We show that flexible protrusions on the rootlet surface, which emanate from the cross-striations, connect to intracellular membranes. In purified rootlets, the striations were classified into amorphous (A)-bands, associated with accumulations on the rootlet surface, and discrete (D)-bands corresponding to punctate lines of density that run through the rootlet. These striations connect a flexible network of longitudinal filaments. Subtomogram averaging suggests the filaments consist of two intertwined coiled coils. The rootlet’s filamentous architecture, with frequent membrane-connecting cross-striations, lends itself well for anchoring large membranes in the cell.

    1. Structural Biology and Molecular Biophysics

    Role of the αC-β4 loop in protein kinase structure and dynamics

    Jian Wu, Nisha A Jonniya ... Susan S Taylor
    Research Article

    Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.