M. Gordon Joyce1,2, Amy Bao1, Man Chen1, Ivelin S. Georgiev1,3,4,5, Li Ou1, Tatsiana Bylund1, Aliaksandr Druz1, Wing-Pui Kong1, Dongjun Peng1, Emily J. Rundlet1, Joseph G. Van Galen1, Shuishu Wang1, Yongping Yang1, Baoshan Zhang1, Gwo-Yu Chuang1, Jason S. McLellan1,6, Barney S. Graham1, John R. Mascola1, Peter D. Kwong1
1
Vaccine Research Center; National Institute of Allergy and Infectious
Diseases; National Institutes of Health; Bethesda, Maryland
2
Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc.;
Bethesda, Maryland
3
Vanderbilt Vaccine Center; Vanderbilt University Medical Center; Nashville,
Tennessee
4
Department of Pathology, Microbiology, and Immunology; Vanderbilt
University Medical Center; Nashville, Tennessee
5
Department of Electrical Engineering and Computer Science; Vanderbilt
University; Nashville, Tennessee
6
Department of Molecular Biosciences; College of Natural Sciences; The
University of Texas at Austin; Austin, Texas
Peter D. Kwong
pdkwong@nih.gov
Background: Respiratory syncytial virus (RSV) subtypes, A and B, co-circulate in annual epidemics and alternate in dominance. We have shown that a subtype A RSV fusion (F) glycoprotein, stabilized in its prefusion conformation by DS-Cav1 mutations, is a promising RSV-vaccine immunogen, capable of boosting RSV-neutralizing titers in healthy adults. In both humans and vaccine-tested animals, neutralizing titers elicited by this subtype A DS-Cav1 immunogen were ~ 2- to 3-fold higher against the homologous subtype A virus than against the heterologous subtype B virus.
Methods: To understand the molecular basis for this subtype difference, we introduced DS-Cav1 mutations into RSV strain B18537 F, determined the trimeric crystal structure, and carried out immunogenicity studies.
Results: The B18537 DS-Cav1 F structure at 2-Å resolution afforded a precise delineation of prefusion F characteristics, including those of antigenic site Ø, a key trimer-apex site. Structural comparison with the subtype A prefusion F indicated 11% of surface residues to be different, with an alpha-carbon root-mean-square deviation (RMSD) of 1.2 Å; antigenic site Ø, however, differed in 23% of its surface residues and had an alpha-carbon RMSD of 2.2 Å. Immunization of vaccine-tested animals with DS-Cav1-stabilized B18537 F induced neutralizing responses ~100-fold higher than with postfusion B18537 F. Notably, elicited responses neutralized RSV subtypes A and B at similar levels and were directed towards both conserved equatorial and diverse apical regions.
Conclusion: We propose that structural differences in apical and equatorial sites–coupled to differently focused immune responses–provide a molecular explanation for observed differences in elicited subtype A and B neutralizing responses.
Keywords: antigenic site; crystal structure; epitope; fusion glycoprotein; immunogenicity; neutralization; RSV subtype; vaccine
Respiratory syncytial virus (RSV) was first identified by Chanock and coworkers in 1957 [1, 2] as the etiological agent responsible for the majority of acute respiratory infections in children. In the United States, RSV is the leading cause of pediatric intensive care hospitalizations [3]. Worldwide, RSV annually causes approximately 3 million acute lower respiratory infections in children and is typified by recurrent infections in infants [4]. Two major RSV subtypes, A and B, circulate with 1-2 epidemic periods per year depending on the climate [5, 6]. Subtype A infections appear to be more prevalent, with disease from both subtypes of similar but variable severity [7, 8]. The RSV fusion (F) glycoprotein is a dominant antigen and a prominent vaccine target [9-11]. Although most subtyping is done by sequencing of the RSV G glycoprotein, molecular phylogeny analysis of RSV F from circulating strains has helped to define RSV evolution and subtype prevalence over time [12-14]. A review by Melero and Moore [15] summarizes the differences between RSV subtype infection and the resulting effect on pathogenesis and immunity. RSV infection tends to display a regular pattern of periodicity of subtype A or B dominance in a given geographic area [16-22] with adult re-infection rates of 5%-10% [23-25]. Analysis of infection and re-infection, from a birth cohort in rural Kenya or over 8 seasons of RSV infection in Japan in a pediatric cohort, shows that 80% of RSV re-infections are due to either subtype or genotype differences [26, 27].
Analysis of human sera indicates that the majority of neutralizing antibodies elicited by natural RSV infection target the prefusion form of the F glycoprotein [28, 29], which is metastable and transitions to a more stable postfusion form. The prefusion conformation can be stabilized by the structure-based design of disulfide (DS) and cavity-filling (Cav1) mutations [30] and further stabilized with additional mutations [31, 32]. These mutations introduce covalent linkages or steric hindrance that prevent structural rearrangements required for the transition from prefusion to postfusion. A DS-Cav1-stablilized subtype A RSV F elicited levels of RSV-neutralizing antibodies in mice and macaques many times the protective threshold [30]; however, the level of elicited neutralization against a subtype B virus was only one-third the level observed against the homologous subtype A strain. In a phase I trial in adult humans [33] immunizations with a DS-Cav1-stabilized subtype A RSV F induced a higher vaccine-induced serum neutralization response against subtype A virus than subtype B virus, similar to that observed in animal models.
Several RSV-targeting monoclonal antibodies that are subtype specific, such as the highly potent monoclonal antibodies RSE20 [34] and 5C4 [35], have been identified. There are substantial differences in subtype neutralization elicited by RSV F stabilized in the prefusion conformation [35, 36], indicating that control of RSV disease may require an improved understanding of F subtype differences. To provide structural and immunogenic bases of F subtype differences, we crystallized and determined the structure of a soluble prefusion RSV F from subtype B strain B18537, stabilized by DS-Cav1 mutations. We compared this structure to prior DS-Cav1 structures [30, 37], and analyzed sequence diversity, B factor diversity, and structural variations at the apical and equatorial antigenic sites. We also assessed the immunogenicity of this subtype B prefusion-stabilized F trimer. The results define RSV F subtype differences at the atomic-level and provide insight into subtype-specific immunogenicity of the prefusion RSV F glycoprotein.
We purchased 293F human embryonic kidney cell lines from Invitrogen (Carlsbad, CA) which were maintained in Freestyle 293 expression medium (Invitrogen). Monoclonal antibodies used in binding studies, AM22, D25, 5C4, motavizumab, and MPE8, were expressed by transient co-transfection of both the heavy and light chain gene-containing plasmids into 293F cells in suspension culture at 37°C. Supernatant was harvested after 6 days, and antibodies were purified from a protein A or protein G affinity column (Qiagen). Antigen-binding fragment (Fab) was produced by endoproteinase Lys-C digestion of purified antibody in 25mM Tris-Cl, 1mM EDTA (pH 8) for 6 hours at 37oC. The reaction was stopped with 1mM Nα-p-tosyl-L-lysine chloromethyl ketone and 0.4mM leupeptin and the mixture passed over a protein A or protein G column to remove the Fc fragment. Flowthrough fractions containing Fab were concentrated and further purified to remove Lys-C by size-exclusion chromatography (16/60 Superdex-200, GE) in phosphate buffered saline (PBS) or a buffer containing 350mM NaCl, 2.5mM Tris pH7.1, 0.02% NaN3.
The RSV F gene constructs from subtype B strains B18537 and B1 were made by synthesis of DNA containing codons 1 to 513 from Uniprot sequence P13843 and O36634, respectively, followed by a T4 fibritin trimerization motif [38], thrombin cleavage site, 6x His-tag, and a StreptagII tag. To stabilize the prefusion state, the F gene was altered to introduce “DS-Cav1” mutations, Ser155Cys, Ser190Phe, Val207Leu, and Ser290Cys. The postfusion RSV B18537 F was made by gene synthesis with codons for fusion peptide residues 137 to 146 deleted from Uniprot sequence P13843 followed by a thrombin cleavage site, 6xHis-tag, and a StreptagII tag.
Soluble RSV F proteins were expressed by transient transfection in Expi293F cells using TrueFect-Max (United BioSystems, MD). The culture supernatant was harvested 5 days post transfection, centrifuged at 10,000g to remove cell debris, and then sterile filtered. The supernatant was buffer exchanged and concentrated using tangential flow filtration. RSV F glycoprotein was purified by nickel- (Roche) and Strep-Tactin-affinity (IBA lifesciences) chromatography. Purification tags were removed by thrombin digestion overnight at 4oC, and the protein was further purified by size-exclusion chromatography.
Robotic screening yielded crystals with the RSV strain B18537 F glycoprotein in 1.2M ammonium sulfate, 0.1M sodium acetate pH5.5, 0.1M lithium sulfate. Crystallization conditions were manually optimized in hanging drops by mixing 1 μL of protein complex with 1 μL of the reservoir solution. Crystals were harvested and flash-cooled in liquid nitrogen using 3.2M ammonium sulfate as a cryoprotectant. Data were collected at a wavelength of 1.00 Å at the Southeast Regional Collaborative Access Team (SER-CAT) beamline ID-22 (Advanced Photon Source, Argonne National Laboratory).
Diffraction data were processed using HKL2000, and crystal unit cell analysis [39] indicated 1 RSV F protomer per asymmetric unit. Molecular replacement searches were carried out with PHASER [40], using as a search model the RSV F subtype A structure (PDB ID:4MMS) [30]. Model building was carried out with Coot [41], and refinement was performed with Refmac5 [42], PHENIX [43], and Buster-TNT [44, 45].
The final refined structures were analyzed with MolProbity [46]. Structural comparison and alignments were carried out using LSQKAB [47, 48]. Surface area analyses and contact residue analysis were carried out using PyMOL [49], Pdbsum [50] and PISA [51]. All structural images were created using PyMOL.
Mice were housed and cared for in a facility accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International at the Vaccine Research Center (VRC), National Institute of Allergy and Infectious Diseases, National Institutes of Health in accordance with local, state, federal, and institute policies. Animal experiments were performed in compliance with the Animal Welfare Act requirements for environment enhancement which were adequate to promote the psychological and physical well-being of small animals. The animal care protocol of this study was approved by the VRC Animal Care and Use Committee.
To assess the immunogenicity of DS-Cav1-stabilized RSV F of B18537, we immunized CB6F1/J mice (10 in each group) by intramuscular injection with 10 μg RSV F combined with 50 μg poly I:C adjuvant at weeks 0 and 3. Serum samples were collected at week 5. RSV F B18537 postfusion was used as a control for a side-by-side comparison.
For RSV neutralization assays, sera were serially diluted and mixed with an equal volume of recombinant mKate-RSV expressing F glycoprotein from strain A2 or B18537 and the Katushka fluorescent protein. The mixtures were incubated at 37ºC for 1 hour, and 50 μL of each was added to 1.5 × 104 HEp-2 cells in 30 μL minimal essential medium in each well of 384-well black optical bottom plates. The plates were incubated for 20-22 hours and read with 588 nm excitation and 635 nm emission (SpectraMax Paradigm, Molecular Devices, CA). The EC50 for each sample was calculated by non-linear regression using GraphPad Prism (GraphPad Software Inc., CA). P values were determined by Student's t-test.
A fortéBio Octet Red384 instrument was used to measure binding kinetics of RSV B18537, B1 and A2 DS-Cav1 F glycoproteins to antibodies AM22, 5C4, D25, MPE8, and motavizumab [34, 35, 52-54]. All assays were performed at 30°C in solid black 96-well plates (Greiner Bio-One), with agitation at 1000 rpm in PBS supplemented with 1% BSA in order to minimize nonspecific interactions [30, 32]. The final volume was 100 μL/well. His-tagged F trimers were captured on Ni-NTA sensors at a capture level between 0.8 and 1 nm with a variability of < 0.1 nm within a row of 8 tips. The biosensor tips were equilibrated for 300 seconds in PBS plus 1% BSA prior to binding measurements in the Fab solutions. Binding reactions proceeded for 300 seconds, followed by 300 seconds dissociation. Dissociation wells were used only once to prevent contamination. Systematic baseline drift was corrected by subtracting the measurements recorded for a sensor loaded with RSV-specific antibodies but incubated with PBS buffer supplemented with 1% BSA. Non-specific binding responses were measured with HIV-specific antibody 2F5 [55], which does not bind RSV F, and subtracted from RSV antibody-response data. Data analysis and curve fitting were done using Octet software version 7.0. Experimental data were fitted with a binding equation describing a 1:1 interaction. Global analyses of the complete data sets assumed binding was reversible (full dissociation) and used nonlinear least squares fitting, allowing a single set of binding parameters to be obtained simultaneously for all concentrations used in each experiment.
Serum antibody binding analysis was carried out in a similar manner to the kinetic analysis described above. RSV B18537 F DS-Cav1 and postfusion forms were immobilized to NTA biosensors through a 6x His-tag. The biosensor probes were then equilibrated for 300 seconds in PBS + 1% BSA buffer. Sera were diluted 50- or 100-fold in PBS + 1% BSA, and binding to the immobilized F proteins was assessed for 300 seconds. To analyze antibody competition for binding F trimers, RSV F-loaded probes were incubated in MPE8 or D25 solutions prior to assays in the diluted sera. Sera depletion was carried out by using 1 μg of DS-Cav1 or postfusion F proteins per 1 μL of animal sera. Parallel correction to subtract non-specific sera binding was carried out by subtracting binding levels of an unloaded probe incubated with the sera.
The RSV F glycoprotein nucleotide sequences used in this study were downloaded from Genbank and the DNA databank of Japan. Accession numbers of utilized sequences are as follows for subtype B strains: JX576729 to JX576762, AF013254 (B1), AY353550 (9320), and for subtype A strains: JQ901447 to JQ901458, JX015479 to JX015499, M74568 (A2), AY911262 (Long), FJ614813 (Line19), NC_001803 (RSS-2), JX482018 to JX482038, JX682715 to JX682823, JX477455 to JX477594.
The evolutionary history was inferred by using the Maximum Likelihood method based on the JTT matrix-based model [56]. The tree with the highest log likelihood (-2466.1477) was calculated. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model, and then selecting the topology with the superior log likelihood value. The analysis involved 352 RSV F sequences. All positions containing gaps or missing data were eliminated. There were 273 positions in the final dataset. Evolutionary analyses were conducted in MEGA6 [57].
Sequence variability for RSV subtypes A and B F protein was analyzed using 2 approaches: Shannon entropy and mutual information [58]. A total of 87 (54 A and 33 B) unique RSV F sequences were used in the analyses. Entropy measures the degree of sequence variability for a given residue position, with higher entropy values corresponding to greater variability for that residue, for example, residue 209 has an entropy score of 0.65127 when both sequences from subtypes A and B are assessed, corresponding to a distribution of 64% Lys and 36% Gln, while residue 276 has an entropy score of 0.59977, corresponding to 29% Asn and 71% Ser. Entropy was computed separately for subtypes A and B, as well as for sequences of both subtypes combined.
Mutual information scores take into account the association between mutations at a given RSV F residue position and the prevalence of a residue within a given subtype. Higher mutual information scores are given for a residue where particular mutations are preferentially found in 1 subtype. For example, residue 209 had a mutual information score of 0.56455, which indicates a substantially different amino acid distribution for that residue in RSV F subtype A versus subtype B. Specifically, in RSV F subtype A sequences, residue 209 is a completely conserved Lys, whereas in subtype B, that residue is 94% Gln and only 6% Lys. On the other hand, residue 276 had a mutual information score of 0.171243, indicating the existence of some, though not as pronounced, differences in amino acid distribution between A and B sequences; residue 276 is 46% Asn and 54% Ser in A sequences, whereas it is 100% Ser in B sequences.
To stabilize the RSV glycoprotein of strain B18537 in its prefusion conformation, we introduced disulfide and cavity-filling mutations (S155C, S190F, V207L, and S290C; termed DS-Cav1 [30]), and appended a trimerization domain from phage T4 (Foldon) [38] to the C terminus of the F ectodomain. The entire coding region was placed into a pAH mammalian expression vector derived from pLEXm [59] and expressed by transient transfection. This DS-Cav1-stabilized RSV F was fully cleaved into the appropriate F2 and F1 polypeptides and eluted from a gel filtration column as a single monodisperse peak at the expected retention volume for a glycosylated trimer (Figure S1). Similar results were obtained for an F glycoprotein from strain B1. The typical yield for these DS-Cav1 subtype B RSV F variants was ~1.5 mg of purified trimer per liter of growth media, similar to that of subtype A strain A2 [30]. The purified RSV B18537 F had similar physical stability to RSV A2 F [32] (Supplementary Figure 1C).
We determined antigenic recognition of RSV F using biolayer interferometry (Octet). The F glycoproteins were immobilized on the biosensor tip, and their affinity to a panel of RSV F reactive antibodies was measured (Table 1). Notably, the recognition of both subtypes A and B was similar for most antibodies, except for 3 site Ø-specific antibodies; recognition by antibodies D25 and AM22 was reduced by > 100-fold for a given subtype B strain while antibody 5C4 recognition was fully ablated. Similar results have been reported between another subtype B strain, B9320, and A2 [60]. B18537 has an Arg residue at position 202 located within antigenic site Ø, which is uncommon in subtype B and possibly explains its recognition by AM22; the B1 strain has the more typical Gln residue at position 202. This Arg202 could play a role in AM22 binding similar to Lys201 of subtype A.
Table 1. Antigenic characterization of RSV DS-Cav1 F glycoproteins.
RSV F Antigen |
Antibody Fab/ Antigenic site |
Binding kinetics |
||
|
kon (1/Ms) |
koff (1/s) |
KD (M) |
|
B18537 |
|
2.99 × 104 |
2.20 × 10-2 |
7.36 × 10-7 |
B1 |
D25/site Ø |
1.54 × 105 |
1.66 × 10-2 |
1.07 × 10-7 |
A2 |
7.76 × 104 |
1.14 × 10-5 |
1.47 × 10-10 |
|
B18537 |
|
1.59 × 104 |
< 1.00 × 10-7 |
< 6.29 × 10-12 |
B1 |
AM22/site Ø |
1.77 × 104 |
2.75 × 10-4 |
1.55 × 10-8 |
A2 |
1.53 × 105 |
1.69 × 10-7 |
1.10 × 10-12 |
|
B18537 |
NB |
NB |
NB |
|
B1 |
5C4/site Ø |
NB |
NB |
NB |
A2 |
2.65 × 105 |
3.52 × 10-3 |
1.33 × 10-8 |
|
B18537 |
8.60 × 104 |
7.41 × 10-5 |
8.62 × 10-10 |
|
B1 |
Motavizumab/site II |
1.34 × 105 |
1.15 × 10-4 |
8.56 × 10-10 |
A2 |
2.10 × 105 |
8.59 × 10-6 |
4.09 × 10-11 |
|
B18537 |
4.79 × 105 |
5.89 × 10-6 |
1.23 × 10-11 |
|
B1 |
MPE8/site III |
2.38 × 105 |
8.30 × 10-5 |
3.49 × 10-10 |
A2 |
8.43 × 104 |
4.04 × 10-5 |
4.79 × 10-10 |
NB = No binding
Crystallization trials of the DS-Cav1-stabilized strain B18537 F
glycoprotein yielded crystals of symmetrical cubes of 50 to 80 μm in size,
at pH 5.5 in space group P4132 with 1 protomer per
asymmetric unit (Table 2). We solved the structure using a starting model
from PDB: 4MMS using molecular replacement with Phaser. After iterative
cycles of model building and TLS refinement with 6 groups, 454 residues of
the RSV F molecule were built in the final refined model. Data collection
and refinement statistics are presented in Table 2.
Table 2. Crystallographic data and refinement statistics.
PDB accession code |
6Q0S |
Growth condition |
1.2M ammonium sulfate, 0.1M lithium sulfate, 0.1M sodium acetate pH 5.5 |
Data collection |
|
Space group |
P 4132 |
Cell constants |
|
a, b, c (Å) α, β, γ (°) |
167.9, 167.9, 167.9 90.0, 90.0, 90.0 |
Wavelength (Å) |
1.00 |
Resolution (Å) |
50.0-1.94 (2.09-2.01, 2.01-1.94) |
R sym |
0.10 (0.648, 0.788) |
I / σI |
11.96 (1.94, 1.12) |
Completeness (%) |
95.9 (96.1, 79.7) |
Redundancy |
4.0 (3.3, 2.4) |
Refinement |
|
Resolution (Å) |
50.0-1.94 |
Unique reflections |
57,616 |
R work / Rfree (%) |
18.74/21.15 |
No. atoms |
|
Protein |
3552 |
Ligand/ion |
5 |
Water |
401 |
B -factors (Å2) |
|
Protein |
46.2 |
Ligand/ion |
78.3 |
Water |
53.2 |
R.m.s. deviations |
|
Bond lengths (Å) |
0.017 |
Bond angles (°) |
1.50 |
Ramachandran |
|
Favored regions (%) |
95.3 |
Allowed regions (%) |
4.5 |
Disallowed regions (%) |
0.2 |
The F glycoprotein formed a homotrimer, each protomer comprising 2 different subunits, F2 and F1 (Figure 1). Residues 26 to 106 of the F2 subunit as well as residues 137 to 509 of the F1 subunit were well defined by the experimental electron density. Each protomer comprised 4 major structural domains (Figure 1A), the membrane proximal α10 helix, the DI domain including the F2 N-terminus, the DII domain (residues 401-463), and the membrane distal DIII domain (residues 49-308) including the F2 C-terminus. Each protomer had major interactions (2250 Å2 involving ~ 75 residues) with neighboring protomers resulting in an interweaved structure, where the DIII domain was located directly above the DII domain of an adjacent protomer. The fusion peptide of each protomer was located in a central cavity, interacting with adjacent protomers. We determined the surface areas of antigenic sites on the prefusion trimer using the previously defined antigenic site residues [35, 61-65] (Figure 1B). Antigenic site III had the largest area at 8,779 Å2, and it overlapped with sites II and IV. Antigenic site V [66], also defined as site VIII [67], overlapped with sites II, III, and IV.
Figure 1. Crystal structure of the fusion glycoprotein of RSV B18537. (A) Ribbon representation of a single protomer of the RSV B18537 F glycoprotein. RSV F domains, antigenic sites, fusion peptide, and subunit termini are labeled. Antigenic site III is discontinuous and overlaps with sites II and IV. The residues of site III overlapping with site II or IV are colored only as site II or IV, respectively. Another known antigenic site, site V, located between sites Ø and II and overlapped with sites Ø, II, and III, is not shown for clarity. (B) Surface representation of the RSV B18537 F trimer with the antigenic sites colored as in (A) and the antigenic surface area (Å2/protomer) indicated. The total trimeric accessible surface area (residues 18 to 509) is 52,225 Å2, calculated using a surface accessibility radius of 1.4 Å.
Most of the DS-Cav1 mutations were clearly visible in electron density maps, as exemplified by the S190F mutation shown in Figure S2. The S155C and S290C mutations formed a prefusion-stabilizing disulfide bond (Figure 2A). The Cav1 mutations S190F (Figure 2A) and V207L (Figure 2B) stabilized intra-protomer interactions. The interweaved nature of the protomers was apparent at the local level where fusion peptide residues 143 and 144 main-chain hydrogen-bonded with domain II strands β16 of an adjacent protomer (Figure 2C), while Ser146 to Ile148, located just after the fusion peptide, interacted with the F2 C-terminal region of the same protomer (Figure 2D).
Figure 2. Atomic-level details of the RSV B18537 F glycoprotein structure . (A) Prefusion stabilizing DS-Cav1 mutations (S155C, S190F, S290C) are shown with van der Waals surface (red dots) with the adjacent secondary structure including the domain III β-propeller fold shown in ribbon format. (B) Antigenic site Ø located at the apex of each protomer is highlighted in red with side chains in stick representation. The V207L mutation, which is part of the stabilizing mutations, is highlighted in dot representation. (C) Inter-protomer hydrogen bonds between the fusion peptide of one protomer and β strands 16, 17, and 20 from an adjacent protomer are indicated by dotted lines. (D) Interaction between the F2 C-terminus and the fusion peptide within a single protomer. The F2 C-terminus is shown in surface representation and the fusion peptide shown with dotted van der Waals surface.
In general the structures of the RSV subtypes A and B F glycoproteins were very similar. The root-mean-square deviation (RMSD) between the structures of B18537 (determined here) and A2 (PDB ID: 4MMU) was 1.18 Å over the Cα atoms of 454 residues present in both structures. Despite 28 residues being different in the mature sequence between the 2 proteins, their overall secondary structures were well conserved (Figure 3). Most of the differences were at the apex and F2 C-terminal region. In both structures, residues 101 to 106 of the F2 C-terminus were located adjacent to the fusion peptide of the same protomer (Figure 3C,D); in the A2 structure, however, this region also interacts with a short helical region from residues 354-358 in an adjacent protomer, resulting in a shift of 8.5 Å of Arg106 between the 2 structures (Figure 3A,D). Both the F2 C-terminus and F1 N-terminus were in the central cavity of the RSV F glycoproteins.
Figure 3. Structural comparison of RSV strain B18537 DS-Cav1 F glycoprotein with A2 DS-Cav1 F (PDB ID: 4MMU) and with B9320 DS-Cav1 F (PDB ID: 5UDE). (A) Structural superposition of a protomer of B18537 F (colored as in Figure 1) with that of A2 (gray) shown in ribbon representation. The 2 structures have a RMSD of 1.18 Å over 454 Cα atoms. There is a shift of 8.5 Å between the F2 C-termini of B18537 and A2 F glycoproteins. The F2 C-terminus Arg106 residues of the 2 strains are shown in stick representation with a dashed line indicating the distance between them. (B) Structural superposition of a protomer of B18537 F with that of B9320 (gray). The 2 structures have a RMSD of 1.06 Å over 444 Cα atoms. (C) Close-up view of the F2 C-terminus of the B18537 F structure and its interactions with the fusion peptide and other parts of the same protomer. Shown in pink is an adjacent protomer. (D) Close-up view of the F2 C-terminus of the A2 F structure and its interactions with an adjacent protomer (cyan). In (C) and (D), yellow dashes indicate hydrogen-bonds.
Comparison between the 2 structures of RSV F subtype B, B18537 and B9320 [37], showed that they had an RMSD of 1.06 Å over 444 Cα atoms present in both structures, even though the 2 proteins are >99% identical in sequence. Similar to that described above between subtypes A and B, most of the differences were at the apical region and termini where the B factors were high (Figure 4), indicating these regions to be structurally flexible. Some of the observed variation was likely to be due in part to the low resolution of the B9320 structure, whose F2 C-terminal region is partially disordered, missing residues 100-106, which are present in the B18537 F structure. One notable difference was at the F1 N-terminal region, whose sequence is conserved between the 2 subtype B strains (detailed below). This region is involved in the interface between protomers.
Figure 4. DS-Cav1 stabilized RSV F glycoproteins have higher flexibility at the apical region. (A) Worm diagrams of the RSV F trimers with thickness proportional to B factor values. The figures were prepared with PyMOL from crystal structures of B18537, B9320 (5UDE), and A2 (4MMU). The positions of antigenic sites Ø and III on the front face are marked with red and pale-yellow contours on the diagram of the B18537 structure. (B) Pairwise RMSD between the 3 structures at 2 antigenic sites with the corresponding sequence identity of the epitope residues. RMSD was calculated based on Cα atoms.
Analysis of 352 RSV F sequences (partial or complete; 109 complete B sequences and 164 complete A sequences) showed that subtype A and subtype B sequences were separated into 2 distinct phylogenetic branches (Figure 5A). From this phylogenetic tree, it was clear that RSV subtypes could be categorized based on F sequences, in addition to the conventional G sequence-based subtyping [68]. Subtypes A and B F glycoproteins differed by ~5%-10% from each other in amino acid sequence, and B18537 F differed from A2 F by 38 residues, including 10 in the 27-residue peptide between furin-cleavage sites (Figure 5B). Within each subtype there was significantly less sequence variation, while there was slightly more diversity among A strains than among B strains. The high-resolution structure of B18537 F allowed for a clear definition of the secondary structure of the RSV F glycoproteins, thereby updating the nomenclature of the loop, helices, and β-strands (Figure 5B). There are several well-characterized antigenic sites on the F glycoprotein [69], site I (residues 387-392) targeted by 2F and 131-2a, site II (residues 254-276) targeted by palivizumab, site III (residues 45, 50-54, 150, 154, 178, 180, 186-187, 261-273, 305-307, 309-312, 344-347, 364, and 377 from one protomer and residues 425, 427-431, 448-449, 456 and 458 from an adjacent protomer) targeted by MPE8, site IV (residues 427-436) targeted by 101F, and site Ø (residues 62-69 and 196-209) targeted by D25, AM22, and 5C4. Sequence entropy analysis indicated differences among the strains to be evenly spaced throughout the sequence. There was only 1 residue with significant entropy in sites I (residue 389), II (residue 276), and III (residue 305), while none was observed in site IV. By contrast, there were 5 residues in the prefusion-specific site Ø with significant entropy (Figure 5B). This scarcity of sequence variation in most antigenic sites except site Ø has been observed in an analysis of RSV clinical isolates [70].
Figure 5. RSV F phylogeny, subtype sequence alignment, secondary structure definition, and sequence entropy. (A) Molecular phylogenetic analysis of RSV F sequences by Maximum Likelihood. The tree was drawn to scale, with branch lengths measured in the number of substitutions per site. Strains B1, B18537, B9320, and A2, are labeled. (B) Sequence alignment of RSV F proteins from B18537, B9320 and A2. Antigenic sites Ø, I, II, III, and IV are marked above the sequences with color-coded lines. Site III (colored pale yellow) is discontinuous and has substantial overlaps with sites II and IV. Residues that differ between the strains are highlighted in purple and those also found in antigenic sites are highlighted in bold. The DS-Cav1 mutations (S155C, S190F, V207L, S290C) are highlighted in green. Secondary structural elements of RSV F as defined from the B18537 structure are shown below the sequence alignment. Sequence entropy calculated from 190 RSV F sequences is plotted above the sequence with a color gradient according to the color key shown at the lower right. Residues between 2 arrows (residues 110-136) are removed by protease digestion and therefore not present in the mature proteins.
Shannon entropy and mutual information analysis of subtypes A and B F sequences indicated that subtype-specific residues accounted for an approximate total of 5700 Å2 of solvent accessible area, which is approximately 11% of the F trimer surface (Figure 6). These surface changes were spread throughout the protein, but there were 2 regions where differences clustered. The first region (residues 42, 326, 380, 384, 466, and 477) was in the DII domain, which is similar in structure in both prefusion and postfusion conformations. However, to date, there have been no neutralizing antibodies mapped to this region. The second region included antigenic site Ø, which has large structural changes between the prefusion and postfusion conformations, and, as noted above, is a focus of potent RSV-neutralizing human antibodies.
Figure 6. Surface antigenic differences between subtypes A and B. Surface antigenicity does not differ significantly within a subtype, but significant differences are observed between subtypes, with antigenic site Ø having approximately 23% surface variation. (A) Sequence entropy mapped on protein surface with a purple color scale as shown. (Left) Residue differences found within subtype B are mapped on the B18537 F structure; (Center) differences within subtype A are mapped on the A2 F structure; (Right) differences between subtypes A and B are mapped on the B18537 F structure. The antigenic sites are outlined with borders colored as in Figure 1. (B) Mutual information comparison of subtype B and subtype A sequences mapped on the B18537 F structure. Total antigenic site surface areas (Å2) as shown in Figure 1 are indicated to the right and colored the same as the antigenic site outline. Differences between subtypes are indicated in black fonts for each antigenic site. (C) Zoom-in view of subtype differences at antigenic site Ø as defined by sequence entropy. (D) Antigenic site Ø differences between subtypes A and B as defined by a mutual information metric. In (C) and (D), transparent trimer surface is shown overlaid on the cartoon representation with side chains shown in sticks for residues with significant difference.
Antigenic site I antibodies 2F, 44F, and 45F are subtype specific, with neutralization dependent on the residue 389 (Figure 5B). This residue is Pro in subtype A strains sensitive to 2F, 44F, and 45F, while subtype A strains with Ser/Leu/His at 389 are resistant; Ser389 is conserved in all subtype B strains, preventing their neutralization by 2F, 44F, or 45F [71, 72]. In antigenic site II, Ser276 is highly common in subtype B including B18537 [14], whereas Asn is found at 276 in A2. Residue 276 is at the edge of the site II epitope with its side chain facing away from the bound antibody and making minimal contacts [61]; this residue is not involved in recognition by site II-specific antibodies. As shown in Table 1, the antigenic recognition of B18537 F by motavizumab (an optimized version of palivizumab) was comparable to that of A2 F. There have been reports of subtype B resistance to palivizumab, but in addition to the N276S change, there is a requirement for a K272E mutation [73] or N268I and N276Y mutations [74]. Lastly, the sequence of antigenic site IV on the DII domain was fully conserved in all strains analyzed.
The prefusion-specific antigenic site Ø, which is located at the apex of the F glycoprotein, differed between B18537 and A2 in 5 of the 22 residues of this antigenic site (Figures 5B and 6), and structural alignment analysis indicated an RMSD of 2.2 Å in this region (Figure 4B). This subtype difference at site Ø was consistent with results from antigenic analysis; 2 of the 3 site Ø-specific antibodies revealed a significant difference in binding affinity (Table 1). D25 had an affinity ~100-fold lower against B18537 F than against A2 F; 5C4 was unable to bind B18537 F, whereas it bound A2 F with a nanomolar affinity. Structural variation at site Ø is partly due to structural flexibility at the apex of RSV F prefusion structures, as revealed by the high B factors at this region for all 3 structures (Figure 4A). The RMSD of site Ø residues is greater between B18537 and B9320 than between B18537 and A2, despite the sequence identity being greater between B18537 and B9320 (Figure 4B). Thus, antigenic differences were minor between strains B18537 and A2 F glycoproteins, except at antigenic site Ø, where residues of higher entropy clustered and substantial differences in affinity were observed for several of the antibodies that target this site.
To evaluate the potential of RSV subtype B F glycoprotein as a vaccine candidate, we immunized mice with prefusion-stabilized B18537 DS-Cav1 F and analyzed the serum immunogenicity (Figure 7). The DS-Cav1-immunized sera had more than 100-fold higher titers than those immunized with B18537 postfusion F, when assessed against either the homologous B18537 virus or against the heterologous A2 virus. In contrast to the immunization results from subtype A DS-Cav1 [30], the prefusion F subtype B-immunized sera neutralized both subtypes A and B at comparable levels (Figure 7B). Assessment of serum antigenicity for binding to prefusion or postfusion probes indicated serum responses to be specific to the conformation of both immunogen and probe (Figure 7C). We further delineated the prefusion-specific responses, by using prefusion-specific antibodies MPE8 (which binds to an equatorial site) and D25 (which binds to an apical site). We observed MPE8 at 1μM to reduce substantially the recognition by DS-Cav1-immunized sera. D25, which binds B18537 F DS-Cav1 at a lower affinity, was able to reduce the recognition by DS-Cav1-immunized sera to a similar level at 5μM, a concentration substantially higher than its KD (Table 1). These results suggested that the immunization-elicited responses to B18537 F DS-Cav1 were likely to be directed equally to both apical and equatorial sites.
Figure 7. Prefusion-stabilized RSV B18537 DS-Cav1 F elicits high neutralization titers against both subtypes A and B viruses. (A) Mice were immunized with either DS-Cav1 or postfusion B18537 F glycoprotein at weeks 0 and 3. Serum samples were taken at week 5. (B) Mouse sera were assayed for neutralization against subtype A and B strains. Each dot represents the serum titer from a single mouse. Average titer of each group is marked with a line and also labeled on the horizontal axis. (C) Sera from mice immunized with RSV B18537 F prefusion DS-Cav1 or postfusion immunogen were assessed for binding to the immobilized RSV B18537 F DS-Cav1, postfusion, DS-Cav1 bound by MPE8, or DS-Cav1 bound by D25. Average response of each group is marked with a red line. Statistical analyses in (B) and (C) were calculated by Mann-Whitney 2-tailed t test, ****: P < 0.001.
Production of a safe and effective vaccine against RSV is a major worldwide health objective, and elicitation of neutralizing antibodies that target the F glycoprotein may be the most likely avenue for achieving this goal (reviewed in [9]). The structure of the prefusion form of the RSV B18537 F glycoprotein provides insights to the antigenic differences and immune responses observed in patients. Of the antigenic sites defined to date, site Ø is the most heavily altered between the 2 major subtypes. These alterations do not change the overall secondary structure of the antigenic site, but they do affect antibody recognition. Given that this site is targeted by potent RSV-neutralizing antibodies that appear to make up the majority of the RSV-neutralizing response in humans, site Ø differences may provide an explanation for the lack of subtype cross-protection. Further understanding of site Ø-specific antibodies and their mechanism of recognizing F glycoproteins from both subtypes may be helpful to further efforts at vaccine design, as the isolation and characterization of antibodies from immunized animals or from infected patients should allow us to understand why infection with one subtype does not fully protect infection from the other subtype.
The RSV-neutralizing response elicited by B18537 prefusion-stabilized F in mice appeared not to be focused entirely on site Ø. Rather, equal responses against both equatorial and apical sites were observed, strongly competed by both MPE8 and D25. This targeting, together with the high conservation of sequence and structure in the equatorial regions, provide an explanation for the similar neutralizing responses elicited by B18537 DS-Cav1 F against both RSV subtypes. It remains to be seen if similar equatorial responses will be elicited by B18537 DS-Cav1 F in humans.
The availability of prefusion-stabilized RSV F glycoproteins from subtypes A and B should enable the assessment of dual subtype immunization mixtures, the sequential immunizations of subtypes, as well as other immunization regimens aimed at generating improved neutralizing titers against both subtypes. In terms of optimizing the vaccine-generated protective response, the high-resolution RSV B18537 structure should allow vaccine designers to incorporate subtype differences into structure-based efforts. The RSV F prefusion structures are well conserved between subtypes, with the observed structural differences resulting from sequence variation as well as the high flexibility at the apical regions (Figure 8). It may be possible to generate an RSV F glycoprotein with a chimeric antigenic site Ø that can elicit improved neutralization breadth against both subtypes A and B strains. A self-assembling nanoparticle displaying DS-Cav1-stabilized RSV F induces 10-fold more potent serum neutralizing responses than the trimeric DS-Cav1 RSV F [75]; that use of a dual-component nanoparticle, such as insect ferritin [76], to present prefusion-stabilized RSV F of both subtypes A and B on a single nanoparticle, could increase the potency of responses against both subtypes.
Figure 8. Schematic of RSV F highlighting structural features and subtype differences in immunogenicity. A composite RSV F trimer comprising a subunit each from the structures of B18537, B9320 (PDB 5UDE), and A2 (PDB 4MMU) shows a shape similar to the individual structures. Each subunit is shown as worm representation with thickness indicating B factor values. On the subunit of B18537 is mapped the RMSD between the structures of B18537 and A2 F in a gradient of red; on that of B9320, the RMSD between the structures of B18537 and B9320 F in a gradient of green; and on that of A2 is mapped the sequence entropy of RSV F subtypes A and B sequences in a gradient of purple. The positions of sites Ø and III on the front face are marked with contours in red and pale-yellow, respectively.
The authors thank C.-H. Shen for assistance with phylogenetic analysis and J. Stuckey for assistance with figures. We appreciate the insightful comments and discussions from the members of the Structural Biology Section, Structural Bioinformatics Section, Laboratory of Virology, and the Vector Core Section at the Vaccine Research Center. The authors thank T. Beaumont and H. Spits for antibodies D25 and AM22, A. Lanzavecchia for antibody MPE8, N. Xia for antibody 5C4, and MedImmune for palivizumab and motavizumab. The authors thank M. Moore and A. Hotard, Department of Pediatrics, Emory University, and Children's Healthcare of Atlanta, Atlanta, Georgia, for engineered RSVs used in neutralization assays and protocols.
Support for this work was provided by the Intramural Research Program of the Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health. Use of sector 22 (SER-CAT) at the Advanced Photon Source was supported by the US Department of Energy, Basic Energy Sciences, Office of Science, under contract number W-31-109-Eng-38.
The authors declare no conflict of interests.
1. Chanock R, Roizman B, Myers R. Recovery from infants with respiratory illness of a virus related to chimpanzee coryza agent (CCA). I. Isolation, properties and characterization. Am J Hyg. 1957;66(3):281-90. PubMed PMID: 13478578.
2. Chanock RM. Recovery of a new type of myxovirus from infants with croup. Ann N Y Acad Sci. 1957;67(8):287-95. PubMed PMID: 13411967.
3. Hall CB. The burgeoning burden of respiratory syncytial virus among children. Infect Disord Drug Targets. 2012;12(2):92-7. PubMed PMID: 22335498.
4. Nair H, Nokes DJ, Gessner BD, Dherani M, Madhi SA, Singleton RJ, O'Brien KL, Roca A, Wright PF, Bruce N, Chandran A, Theodoratou E, Sutanto A, Sedyaningsih ER, Ngama M, Munywoki PK, Kartasasmita C, Simões EAF, Rudan I, Weber MW, Campbell H. Global burden of acute lower respiratory infections due to respiratory syncytial virus in young children: a systematic review and meta-analysis. The Lancet. 2010;375(9725):1545-55.
5. Haynes AK, Manangan AP, Iwane MK, Sturm-Ramirez K, Homaira N, Brooks WA, Luby S, Rahman M, Klena JD, Zhang Y, Yu H, Zhan F, Dueger E, Mansour AM, Azazzy N, McCracken JP, Bryan JP, Lopez MR, Burton DC, Bigogo G, Breiman RF, Feikin DR, Njenga K, Montgomery J, Cohen AL, Moyes J, Pretorius M, Cohen C, Venter M, Chittaganpitch M, Thamthitiwat S, Sawatwong P, Baggett HC, Luber G, Gerber SI. Respiratory syncytial virus circulation in seven countries with Global Disease Detection Regional Centers. J Infect Dis. 2013;208 Suppl 3:S246-54. PubMed PMID: 24265484. doi: 10.1093/infdis/jit515
6. Borchers AT, Chang C, Gershwin ME, Gershwin LJ. Respiratory syncytial virus--a comprehensive review. Clin Rev Allergy Immunol. 2013;45(3):331-79. PubMed PMID: 23575961. doi: 10.1007/s12016-013-8368-9
7. Martinello RA, Chen MD, Weibel C, Kahn JS. Correlation between respiratory syncytial virus genotype and severity of illness. J Infect Dis. 2002;186(6):839-42. PubMed PMID: 12198620. doi: 10.1086/342414
8. Jafri HS, Wu X, Makari D, Henrickson KJ. Distribution of respiratory syncytial virus subtypes A and B among infants presenting to the emergency department with lower respiratory tract infection or apnea. Pediatr Infect Dis J. 2013;32(4):335-40. PubMed PMID: 23337904. doi: 10.1097/INF.0b013e318282603a
9. Ruckwardt TJ, Morabito KM, Graham BS. Immunological Lessons from Respiratory Syncytial Virus Vaccine Development. Immunity. 2019;51(3):429-42. PubMed PMID: 31533056. doi: 10.1016/j.immuni.2019.08.007
10. Swanson KA, Balabanis K, Xie Y, Aggarwal Y, Palomo C, Mas V, Metrick C, Yang H, Shaw CA, Melero JA, Dormitzer PR, Carfi A. A monomeric uncleaved respiratory syncytial virus F antigen retains prefusion-specific neutralizing epitopes. J Virol. 2014;88(20):11802-10. PubMed PMID: 25078705. Pubmed Central PMCID: PMC4178725. doi: 10.1128/JVI.01225-14
11. Brock LG, Liu X, Liang B, Lingemann M, Liu X, Herbert R, Hackenberg AD, Buchholz UJ, Collins PL, Munir S. Murine Pneumonia Virus Expressing the Fusion Glycoprotein of Human Respiratory Syncytial Virus from an Added Gene Is Highly Attenuated and Immunogenic in Rhesus Macaques. J Virol. 2018;92(17). PubMed PMID: 29925656. Pubmed Central PMCID: PMC6096832. doi: 10.1128/JVI.00723-18
12. Chi H, Liu HF, Weng LC, Wang NY, Chiu NC, Lai MJ, Lin YC, Chiu YY, Hsieh WS, Huang LM. Molecular epidemiology and phylodynamics of the human respiratory syncytial virus fusion protein in northern Taiwan. PLoS One. 2013;8(5):e64012. PubMed PMID: 23734183. Pubmed Central PMCID: 3667090. doi: 10.1371/journal.pone.0064012
13. Tan L, Coenjaerts FE, Houspie L, Viveen MC, van Bleek GM, Wiertz EJ, Martin DP, Lemey P. The comparative genomics of human respiratory syncytial virus subgroups A and B: genetic variability and molecular evolutionary dynamics. J Virol. 2013;87(14):8213-26. PubMed PMID: 23698290. Pubmed Central PMCID: 3700225. doi: 10.1128/JVI.03278-12
14. Xia Q, Zhou L, Peng C, Hao R, Ni K, Zang N, Ren L, Deng Y, Xie X, He L, Tian D, Wang L, Huang A, Zhao Y, Zhao X, Fu Z, Tu W, Liu E. Detection of respiratory syncytial virus fusion protein variants between 2009 and 2012 in China. Arch Virol. 2014;159(5):1089-98. PubMed PMID: 24297488. doi: 10.1007/s00705-013-1870-9
15. Melero JA, Moore ML. Influence of respiratory syncytial virus strain differences on pathogenesis and immunity. Curr Top Microbiol Immunol. 2013;372:59-82. PubMed PMID: 24362684. doi: 10.1007/978-3-642-38919-1_3
16. Waris M. Pattern of respiratory syncytial virus epidemics in Finland: two-year cycles with alternating prevalence of groups A and B. J Infect Dis. 1991;163(3):464-9. PubMed PMID: 1995719.
17. Zlateva KT, Vijgen L, Dekeersmaeker N, Naranjo C, Van Ranst M. Subgroup prevalence and genotype circulation patterns of human respiratory syncytial virus in Belgium during ten successive epidemic seasons. J Clin Microbiol. 2007;45(9):3022-30. PubMed PMID: 17609323. Pubmed Central PMCID: 2045289. doi: 10.1128/JCM.00339-07
18. Reiche J, Schweiger B. Genetic variability of group A human respiratory syncytial virus strains circulating in Germany from 1998 to 2007. J Clin Microbiol. 2009;47(6):1800-10. PubMed PMID: 19386848. Pubmed Central PMCID: 2691087. doi: 10.1128/JCM.02286-08
19. Shobugawa Y, Saito R, Sano Y, Zaraket H, Suzuki Y, Kumaki A, Dapat I, Oguma T, Yamaguchi M, Suzuki H. Emerging genotypes of human respiratory syncytial virus subgroup A among patients in Japan. J Clin Microbiol. 2009;47(8):2475-82. PubMed PMID: 19553576. Pubmed Central PMCID: 2725673. doi: 10.1128/JCM.00115-09
20. Viegas M, Mistchenko AS. Molecular epidemiology of human respiratory syncytial virus subgroup A over a six-year period (1999-2004) in Argentina. J Med Virol. 2005;77(2):302-10. PubMed PMID: 16121371. doi: 10.1002/jmv.20457
21. Scott PD, Ochola R, Ngama M, Okiro EA, Nokes DJ, Medley GF, Cane PA. Molecular epidemiology of respiratory syncytial virus in Kilifi district, Kenya. J Med Virol. 2004;74(2):344-54. PubMed PMID: 15332285. doi: 10.1002/jmv.20183
22. Arnott A, Vong S, Mardy S, Chu S, Naughtin M, Sovann L, Buecher C, Beaute J, Rith S, Borand L, Asgari N, Frutos R, Guillard B, Touch S, Deubel V, Buchy P. A study of the genetic variability of human respiratory syncytial virus (HRSV) in Cambodia reveals the existence of a new HRSV group B genotype. J Clin Microbiol. 2011;49(10):3504-13. PubMed PMID: 21865418. Pubmed Central PMCID: 3187327. doi: 10.1128/JCM.01131-11
23. Beem M. Repeated infections with respiratory syncytial virus. J Immunol. 1967;98(6):1115-22. PubMed PMID: 4290804.
24. Falsey AR. Respiratory syncytial virus infection in adults. Semin Respir Crit Care Med. 2007;28(2):171-81. PubMed PMID: 17458771. doi: 10.1055/s-2007-976489
25. Falsey AR, Hennessey PA, Formica MA, Cox C, Walsh EE. Respiratory syncytial virus infection in elderly and high-risk adults. N Engl J Med. 2005;352(17):1749-59. PubMed PMID: 15858184. doi: 10.1056/NEJMoa043951
26. Agoti CN, Mwihuri AG, Sande CJ, Onyango CO, Medley GF, Cane PA, Nokes DJ. Genetic relatedness of infecting and reinfecting respiratory syncytial virus strains identified in a birth cohort from rural Kenya. J Infect Dis. 2012;206(10):1532-41. PubMed PMID: 22966119. Pubmed Central PMCID: 3475639. doi: 10.1093/infdis/jis570
27. Yamaguchi M, Sano Y, Dapat IC, Saito R, Suzuki Y, Kumaki A, Shobugawa Y, Dapat C, Uchiyama M, Suzuki H. High frequency of repeated infections due to emerging genotypes of human respiratory syncytial viruses among children during eight successive epidemic seasons in Japan. J Clin Microbiol. 2011;49(3):1034-40. PubMed PMID: 21177891. Pubmed Central PMCID: 3067727. doi: 10.1128/JCM.02132-10
28. Magro M, Mas V, Chappell K, Vazquez M, Cano O, Luque D, Terron MC, Melero JA, Palomo C. Neutralizing antibodies against the preactive form of respiratory syncytial virus fusion protein offer unique possibilities for clinical intervention. Proc Natl Acad Sci U S A. 2012;109(8):3089-94. PubMed PMID: 22323598. Pubmed Central PMCID: 3286924. doi: 10.1073/pnas.1115941109
29. Ngwuta JO, Chen M, Modjarrad K, Joyce MG, Kanekiyo M, Kumar A, Yassine HM, Moin SM, Killikelly AM, Chuang GY, Druz A, Georgiev IS, Rundlet EJ, Sastry M, Stewart-Jones GB, Yang Y, Zhang B, Nason MC, Capella C, Peeples ME, Ledgerwood JE, McLellan JS, Kwong PD, Graham BS. Prefusion F-specific antibodies determine the magnitude of RSV neutralizing activity in human sera. Sci Transl Med. 2015;7(309):309ra162. PubMed PMID: 26468324. Pubmed Central PMCID: PMC4672383. doi: 10.1126/scitranslmed.aac4241
30. McLellan JS, Chen M, Joyce MG, Sastry M, Stewart-Jones GB, Yang Y, Zhang B, Chen L, Srivatsan S, Zheng A, Zhou T, Graepel KW, Kumar A, Moin S, Boyington JC, Chuang GY, Soto C, Baxa U, Bakker AQ, Spits H, Beaumont T, Zheng Z, Xia N, Ko SY, Todd JP, Rao S, Graham BS, Kwong PD. Structure-based design of a fusion glycoprotein vaccine for respiratory syncytial virus. Science. 2013;342(6158):592-8. PubMed PMID: 24179220. doi: 10.1126/science.1243283
31. Zhang L, Durr E, Galli JD, Cosmi S, Cejas PJ, Luo B, Touch S, Parmet P, Fridman A, Espeseth AS, Bett AJ. Design and characterization of a fusion glycoprotein vaccine for Respiratory Syncytial Virus with improved stability. Vaccine. 2018;36(52):8119-30. PubMed PMID: 30340881. doi: 10.1016/j.vaccine.2018.10.032
32. Joyce MG, Zhang B, Ou L, Chen M, Chuang GY, Druz A, Kong WP, Lai YT, Rundlet EJ, Tsybovsky Y, Yang Y, Georgiev IS, Guttman M, Lees CR, Pancera M, Sastry M, Soto C, Stewart-Jones GBE, Thomas PV, Van Galen JG, Baxa U, Lee KK, Mascola JR, Graham BS, Kwong PD. Iterative structure-based improvement of a fusion-glycoprotein vaccine against RSV. Nat Struct Mol Biol. 2016;23(9):811-20. PubMed PMID: 27478931. Pubmed Central PMCID: PMC5016229. doi: 10.1038/nsmb.3267
33. Crank MC, Ruckwardt TJ, Chen M, Morabito KM, Phung E, Costner PJ, Holman LA, Hickman SP, Berkowitz NM, Gordon IJ, Yamshchikov GV, Gaudinski MR, Kumar A, Chang LA, Moin SM, Hill JP, DiPiazza AT, Schwartz RM, Kueltzo L, Cooper JW, Chen P, Stein JA, Carlton K, Gall JG, Nason MC, Kwong PD, Chen GL, Mascola JR, McLellan JS, Ledgerwood JE, Graham BS, Team VRCS. A proof of concept for structure-based vaccine design targeting RSV in humans. Science. 2019;365(6452):505-9. PubMed PMID: 31371616. doi: 10.1126/science.aav9033
34. Corti D, Bianchi S, Vanzetta F, Minola A, Perez L, Agatic G, Guarino B, Silacci C, Marcandalli J, Marsland BJ, Piralla A, Percivalle E, Sallusto F, Baldanti F, Lanzavecchia A. Cross-neutralization of four paramyxoviruses by a human monoclonal antibody. Nature. 2013;501(7467):439-43. PubMed PMID: 23955151. doi: 10.1038/nature12442
35. McLellan JS, Chen M, Leung S, Graepel KW, Du X, Yang Y, Zhou T, Baxa U, Yasuda E, Beaumont T, Kumar A, Modjarrad K, Zheng Z, Zhao M, Xia N, Kwong PD, Graham BS. Structure of RSV fusion glycoprotein trimer bound to a prefusion-specific neutralizing antibody. Science. 2013;340(6136):1113-7. PubMed PMID: 23618766. Pubmed Central PMCID: PMC4459498. doi: 10.1126/science.1234914
36. Tian D, Battles MB, Moin SM, Chen M, Modjarrad K, Kumar A, Kanekiyo M, Graepel KW, Taher NM, Hotard AL, Moore ML, Zhao M, Zheng ZZ, Xia NS, McLellan JS, Graham BS. Structural basis of respiratory syncytial virus subtype-dependent neutralization by an antibody targeting the fusion glycoprotein. Nat Commun. 2017;8(1):1877. PubMed PMID: 29187732. Pubmed Central PMCID: PMC5707411. doi: 10.1038/s41467-017-01858-w
37. Zhu Q, McLellan JS, Kallewaard NL, Ulbrandt ND, Palaszynski S, Zhang J, Moldt B, Khan A, Svabek C, McAuliffe JM, Wrapp D, Patel NK, Cook KE, Richter BWM, Ryan PC, Yuan AQ, Suzich JA. A highly potent extended half-life antibody as a potential RSV vaccine surrogate for all infants. Sci Transl Med. 2017;9(388). PubMed PMID: 28469033. doi: 10.1126/scitranslmed.aaj1928
38. Frank S, Kammerer RA, Mechling D, Schulthess T, Landwehr R, Bann J, Guo Y, Lustig A, Bächinger HP, Engel J. Stabilization of short collagen-like triple helices by protein engineering. J Mol Biol. 2001;308(5):1081-9.
39. Matthews BW. Solvent content of protein crystals. J Mol Biol. 1968;33(2):491-7. PubMed PMID: 5700707.
40. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ. Phaser crystallographic software. J Appl Crystallogr. 2007;40(4):658-74. doi: 10.1107/S0021889807021206
41. Emsley P, Lohkamp B, Scott WG, Cowtan K. Features and development of Coot. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 4):486-501. PubMed PMID: 20383002. Pubmed Central PMCID: 2852313. doi: 10.1107/S0907444910007493
42. Murshudov GN, Vagin AA, Dodson EJ. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr. 1997;53(Pt 3):240-55. PubMed PMID: 15299926. doi: 10.1107/S0907444996012255
43. Adams PD, Afonine PV, Bunkoczi G, Chen VB, Davis IW, Echols N, Headd JJ, Hung LW, Kapral GJ, Grosse-Kunstleve RW, McCoy AJ, Moriarty NW, Oeffner R, Read RJ, Richardson DC, Richardson JS, Terwilliger TC, Zwart PH. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 2):213-21. PubMed PMID: 20124702. Pubmed Central PMCID: 2815670. doi: 10.1107/S0907444909052925
44. Bricogne G, Blanc E, Brandl M, Flensburg C, Keller P, Paciorek W, Roversi P, Sharff A, Smart OS, Vonrhein C, Womack TO. Buster version 2.10.0. 2011.
45. Smart OS, Womack TO, Flensburg C, Keller P, Paciorek W, Sharff A, Vonrhein C, Bricogne G. Exploiting structure similarity in refinement: automated NCS and target-structure restraints in BUSTER. Acta Crystallogr D Biol Crystallogr. 2012;68(Pt 4):368-80. PubMed PMID: 22505257. Pubmed Central PMCID: 3322596. doi: 10.1107/S0907444911056058
46. Chen VB, Arendall WB, 3rd, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, Murray LW, Richardson JS, Richardson DC. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 1):12-21. PubMed PMID: 20057044. Pubmed Central PMCID: 2803126. doi: 10.1107/S0907444909042073
47. Kabsch W. A solution for the best rotation to relate two sets of vectors. Acta Cryst. 1976;A32:922-3.
48. CCP4. The CCP4 suite: programs for protein crystallography. Acta Crystallogr D Biol Crystallogr. 1994;50(Pt 5):760-3. PubMed PMID: 15299374. doi: 10.1107/S0907444994003112
49. Schrodinger, LLC. The PyMOL Molecular Graphics System, Version 1.8. 2015.
50. Laskowski RA. PDBsum new things. Nucleic Acids Res. 2009;37(Database issue):D355-9. PubMed PMID: 18996896. Pubmed Central PMCID: 2686501. doi: 10.1093/nar/gkn860
51. Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372:774-97.
52. Beaumont T, Bakker AQ, Yasuda E, inventors. RSV specific binding molecule. 2012. Patent Application 12/898,325.
53. Spits H, Beaumont T, inventors. RSV-specific binding molecules and means for producing them. 2010. Patent Application 12/600,950.
54. Wu H, Pfarr DS, Johnson S, Brewah YA, Woods RM, Patel NK, White WI, Young JF, Kiener PA. Development of motavizumab, an ultra-potent antibody for the prevention of respiratory syncytial virus infection in the upper and lower respiratory tract. J Mol Biol. 2007;368(3):652-65. PubMed PMID: 17362988. doi: 10.1016/j.jmb.2007.02.024
55. Ofek G, Tang M, Sambor A, Katinger H, Mascola JR, Wyatt R, Kwong PD. Structure and mechanistic analysis of the anti-human immunodeficiency virus type 1 antibody 2F5 in complex with its gp41 epitope. J Virol. 2004;78(19):10724-37. PubMed PMID: 15367639. Pubmed Central PMCID: 516390. doi: 10.1128/JVI.78.19.10724-10737.2004
56. Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8(3):275-82. PubMed PMID: 1633570.
57. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725-9. PubMed PMID: 24132122. Pubmed Central PMCID: 3840312. doi: 10.1093/molbev/mst197
58. Chuang GY, Acharya P, Schmidt SD, Yang Y, Louder MK, Zhou T, Kwon YD, Pancera M, Bailer RT, Doria-Rose NA, Nussenzweig MC, Mascola JR, Kwong PD, Georgiev IS. Residue-level prediction of HIV-1 antibody epitopes based on neutralization of diverse viral strains. J Virol. 2013;87(18):10047-58. PubMed PMID: 23843642. Pubmed Central PMCID: 3753990. doi: 10.1128/JVI.00984-13
59. Aricescu AR, Lu W, Jones EY. A time- and cost-efficient system for high-level protein production in mammalian cells. Acta Crystallogr D Biol Crystallogr. 2006;62(Pt 10):1243-50. PubMed PMID: 17001101. doi: 10.1107/S0907444906029799
60. Jones HG, Battles MB, Lin CC, Bianchi S, Corti D, McLellan JS. Alternative conformations of a major antigenic site on RSV F. PLoS Pathog. 2019;15(7):e1007944. PubMed PMID: 31306469. Pubmed Central PMCID: PMC6658013 following competing interests: JSM is an inventor on patents entitled "Prefusion RSV F proteins and their use," (US patent No. 9,738,689 and 10,017,543). SB and DC are employees of Vir Biotechnology and hold shares in Vir Biotechnology. doi: 10.1371/journal.ppat.1007944
61. McLellan JS, Chen M, Chang J-S, Yang Y, Kim A, Graham BS, Kwong PD. Structure of a major antigenic site on the respiratory syncytial virus fusion glycoprotein in complex with neutralizing antibody 101F. J Virol. 2010;84(23):12236-44. doi: 10.1128/jvi.01579-10
62. Wen X, Mousa JJ, Bates JT, Lamb RA, Crowe JE, Jr., Jardetzky TS. Structural basis for antibody cross-neutralization of respiratory syncytial virus and human metapneumovirus. Nat Microbiol. 2017;2:16272. PubMed PMID: 28134915. Pubmed Central PMCID: PMC5308794. doi: 10.1038/nmicrobiol.2016.272
63. Anderson LJ, Bingham P, Hierholzer JC. Neutralization of respiratory syncytial virus by individual and mixtures of F and G protein monoclonal antibodies. J Virol. 1988;62(11):4232-8.
64. McLellan JS, Chen M, Kim A, Yang Y, Graham BS, Kwong PD. Structural basis of respiratory syncytial virus neutralization by motavizumab. Nat Struct Mol Biol. 2010;17(2):248-50. PubMed PMID: 20098425. Pubmed Central PMCID: 3050594. doi: 10.1038/nsmb.1723
65. Gilman MS, Castellanos CA, Chen M, Ngwuta JO, Goodwin E, Moin SM, Mas V, Melero JA, Wright PF, Graham BS, McLellan JS, Walker LM. Rapid profiling of RSV antibody repertoires from the memory B cells of naturally infected adult donors. Sci Immunol. 2016;1(6). PubMed PMID: 28111638. Pubmed Central PMCID: PMC5244814. doi: 10.1126/sciimmunol.aaj1879
66. Gilman MSA, Furmanova-Hollenstein P, Pascual G, A BvtW, Langedijk JPM, McLellan JS. Transient opening of trimeric prefusion RSV F proteins. Nat Commun. 2019;10(1):2105. PubMed PMID: 31068578. Pubmed Central PMCID: PMC6506550. doi: 10.1038/s41467-019-09807-5
67. Mousa JJ, Kose N, Matta P, Gilchuk P, Crowe JE, Jr. A novel pre-fusion conformation-specific neutralizing epitope on the respiratory syncytial virus fusion protein. Nat Microbiol. 2017;2:16271. PubMed PMID: 28134924. Pubmed Central PMCID: PMC5463187. doi: 10.1038/nmicrobiol.2016.271
68. Sullender WM. Respiratory syncytial virus genetic and antigenic diversity. Clin Microbiol Rev. 2000;13(1):1-15, table of contents. PubMed PMID: 10627488. Pubmed Central PMCID: 88930.
69. Graham BS. Biological challenges and technological opportunities for respiratory syncytial virus vaccine development. Immunol Rev. 2011;239(1):149-66. PubMed PMID: 21198670. Pubmed Central PMCID: 3023887. doi: 10.1111/j.1600-065X.2010.00972.x [doi]
70. Mas V, Nair H, Campbell H, Melero JA, Williams TC. Antigenic and sequence variability of the human respiratory syncytial virus F glycoprotein compared to related viruses in a comprehensive dataset. Vaccine. 2018;36(45):6660-73. PubMed PMID: 30292456. Pubmed Central PMCID: PMC6203811. doi: 10.1016/j.vaccine.2018.09.056
71. Lopez JA, Bustos R, Orvell C, Berois M, Arbiza J, Garcia-Barreno B, Melero JA. Antigenic structure of human respiratory syncytial virus fusion glycoprotein. J Virol. 1998;72(8):6922-8.
72. Magro M, Andreu D, Gomez-Puertas P, Melero JA, Palomo C. Neutralization of human respiratory syncytial virus infectivity by antibodies and low-molecular-weight compounds targeted against the fusion glycoprotein. J Virol. 2010;84(16):7970-82. doi: 10.1128/jvi.00447-10
73. Adams O, Bonzel L, Kovacevic A, Mayatepek E, Hoehn T, Vogel M. Palivizumab-resistant human respiratory syncytial virus infection in infancy. Clin Infect Dis. 2010;51(2):185-8. PubMed PMID: 20528472. doi: 10.1086/653534
74. Crowe JE, Firestone CY, Crim R, Beeler JA, Coelingh KL, Barbas CF, Burton DR, Chanock RM, Murphy BR. Monoclonal antibody-resistant mutants selected with a respiratory syncytial virus-neutralizing human antibody fab fragment (Fab 19) define a unique epitope on the fusion (F) glycoprotein. Virology. 1998;252(2):373-5. PubMed PMID: 9878616. doi: 10.1006/viro.1998.9462
75. Marcandalli J, Fiala B, Ols S, Perotti M, de van der Schueren W, Snijder J, Hodge E, Benhaim M, Ravichandran R, Carter L, Sheffler W, Brunner L, Lawrenz M, Dubois P, Lanzavecchia A, Sallusto F, Lee KK, Veesler D, Correnti CE, Stewart LJ, Baker D, Lore K, Perez L, King NP. Induction of Potent Neutralizing Antibody Responses by a Designed Protein Nanoparticle Vaccine for Respiratory Syncytial Virus. Cell. 2019;176(6):1420-31 e17. PubMed PMID: 30849373. Pubmed Central PMCID: PMC6424820. doi: 10.1016/j.cell.2019.01.046
76. Hamburger AE, West AP, Jr., Hamburger ZA, Hamburger P, Bjorkman PJ. Crystal structure of a secreted insect ferritin reveals a symmetrical arrangement of heavy and light chains. J Mol Biol. 2005;349(3):558-69. PubMed PMID: 15896348. doi: 10.1016/j.jmb.2005.03.074
Supplementary Figure 1. Purification of RSV strain B18537 DS-Cav1 F glycoprotein. (A) Size-exclusion chromatograph of RSV B18537 F glycoprotein in PBS buffer on a 120 mL Superdex-200 gel filtration column. Elution volume of protein standards with known molecular weights are indicated along the base of the chromatogram. (B) SDS-PAGE of the peak fraction from size-exclusion chromatography. (C) Physical stability of RSV B18537 F glycoprotein in comparison with those of RSV subtype A F. Physical stability was measured by the fractional value of the D25 reactivity retained after exposure to various physical extremes, compared with the initial D25 reactivity, as determined by Octet biolayer interferometry.
Supplementary Figure 2. Electron density of RSV strain B18537 DS-Cav1 F glycoprotein around the Cav1 residue Phe190 shown in stereo image. The 2Fo-Fc electron density map is contoured at 1 σ and shown in gray mesh format. Amino acid residues are shown in stick representation.