Impact of Dengue Virus Serotype 2 Strain Diversity on Serological

Oct 22, 2018 - Dengue is a mosquito-borne disease caused by four dengue virus serotypes (DENV1–4) that are loosely categorized by sequence ...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of South Dakota

Article

Impact of dengue virus serotype 2 strain diversity on serological immune responses to dengue Sarah Keasey, Jessica Smith, Stefan Fernandez, Anna Durbin, Bryan Zhao, and Robert G Ulrich ACS Infect. Dis., Just Accepted Manuscript • DOI: 10.1021/acsinfecdis.8b00185 • Publication Date (Web): 22 Oct 2018 Downloaded from http://pubs.acs.org on October 23, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Impact of dengue virus serotype 2 strain diversity on serological immune responses to dengue

Sarah L. Keasey1,2, †, Jessica L. Smith1, †, Stefan Fernandez3, Anna P. Durbin4, Bryan M. Zhao1, and Robert G. Ulrich1*

1Molecular

and Translational Sciences Division, U.S. Army Medical Research Institute of

Infectious Diseases, 1425 Porter St, Frederick, Maryland 21702, USA 2Department

of Biological Sciences, University of Maryland Baltimore County, 1000 Hilltop

Circle, Baltimore, Maryland 21250, USA 3Armed

Forces Research Institute of Medical Sciences, Bangkok, 10400, Thailand

4Center

for Immunization Research, Johns Hopkins Bloomberg School of Public Health, 624

North Broadway, Room 251, Baltimore, Maryland 21205, USA

†Both

authors contributed equally.

*correspondence: [email protected]; 301 619-4232

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Dengue is a mosquito-borne disease caused by four dengue virus serotypes (DENV1-4) that are loosely categorized by sequence commonalities and antibody recognition profiles. The highly variable envelope protein (E) that is prominently displayed on the surface of DENV is an essential component of vaccines currently under development, yet the impact of using single strains to represent each serotype in tetravalent vaccines has not been adequately studied. We synthesized chimeric E by replacing highly variable residues from vaccine strain PUO-218 with those from 16 DENV2 lineages spanning 60 years of antigen evolution. Examining sera from human and rhesus macaques challenged with single strains of DENV2, antibody-E interactions were markedly inhibited or enhanced by residues mainly focused within a 480Å2 footprint displayed on the E backbone. The striking impact of E diversity on polyclonal immune responses suggests that frequent antigen updates may be necessary for vaccines to counter shifts in circulating strains.

Keywords. Dengue virus, dengue virus envelope, dengue vaccine, antibody response, dengue strain variability, protein microarray

ACS Paragon Plus Environment

Page 2 of 42

Page 3 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

The clinical symptoms of dengue range from mild febrile illness to severe hemorrhagic fever and shock syndrome with capillary leakage1. Four serotypes of the mosquito-borne dengue virus (DENV1-4) circulate in most tropical or semi-tropical regions of the world, and this extensive distribution leads to greater than 500,000 severe cases and 22,000 deaths annually2, 3. DENV are members of the Flavivirus genus of single-stranded RNA viruses, which also includes Zika, yellow fever, and West Nile viruses. Long-lasting immunity to homologous DENV infections is mediated by antibody and T-cell responses to primary infections by one of the DENV serotypes4, while immunity against other serotypes (heterotypic immunity) may be shortlived5-8 or lower in magnitude9, resulting in continuous cycles of infection within susceptible human populations. A further concern is that secondary infection with a heterologous serotype can in some cases lead to increased risk of severe dengue disease, possibly because antibodies to one serotype may enhance infections with heterologous serotypes by promoting viral entry and infection through Fc receptor-expressing cells (antibody-dependent enhancement or ADE)10. Despite similar disease presentations, DENV1-4 differ11 at the proteomic level by 25-40%, and by 3-6% within serotypes12, 13, while nucleic acid substitution rates are generally consistent across serotypes. The circulation of new viral isolates is driven by selection pressures14, 15 that are imposed during the natural cycling of DENV between mosquito vectors and human hosts. The emerging viral isolates often harbor structural variants of the envelope (E) protein7, 16, 17, the dominant surface antigen of DENV. The E and membrane (M) structural proteins of DENV assemble into 180 heterodimers that form a sheath around the nucleocapsid core. The E monomer folds into three domains (EDI, EDII, and EDIII) comprising the ectoprotein and two alpha helices of a membrane-proximal stem. Three parallel E dimers are complexed with the M

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

protein (prM), which contains a transmembrane anchor peptide that is cleaved during maturation of the virus18-20. The extensive diversity of virus strains moving in and out of human populations presents a significant challenge to the development of effective dengue vaccines that must first incorporate antigenic components from all serotypes. There are several competing tetravalent dengue vaccines (TDV) that are approaching licensure for immunization and distribution21. The first TDV approved for use in a limited number of DENV-endemic countries is a live-attenuated virus with prM and E proteins from one prototype strain of each DENV serotype inserted into the remaining proteome of the yellow fever virus 17D vaccine (CYD-TDV). A limited, overall efficacy rate of 56.5-60.8% was reported in the first advanced clinical trials of CYD-TDV, especially for previously seronegative subjects, and the lowest efficacy was observed with the DENV2 serotype (18.7-55.2%)22-24. Whereas new sequence variants of DENV2 are reported almost every year, the monoclonal prM and E proteins of DENV2 that were used for CYD-TDV originated from a virus (PUO-218) that was isolated in 1980 from a patient in Thailand. Although the basic E protein structure is highly conserved among DENV, amino acid residues that are permissive to mutations are likely to play a role in driving adaptation to selection pressures imposed by host immune responses. Here we examine the relationship between serological immune responses to DENV2 strains representative of two separate lineages and antibody interactions with naturally occurring variants of the DENV2 E protein. Chimeric E proteins were constructed by grafting solvent-exposed amino acid residues from strains that represent major branch points in the evolution of DENV2 onto the PUO-218 backbone, and these recombinant antigens were then used to evaluate human and rhesus macaque antibody responses to disparate DENV2 strains.

ACS Paragon Plus Environment

Page 4 of 42

Page 5 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Results E protein evolution of DENV2 The molecular evolution of DENV2 was examined with 266 non-redundant E proteins encoded by fully-sequenced strains that circulated in 36 countries between 1944 and 2010 (Supplementary File 1). The minimal amino acid sequence identity by pair-wise comparison was 91.92%, and similar to levels observed across DENV2 genotypes12, 25. Bayesian methods26-30 were used to determine the specific lineages that emerged from the most recent common DENV2 ancestor (Fig. 1), which was estimated to have originated in the late 1800s (ca 1897) at a time that was approximately 47 years before the oldest (ca 1944) documented strain of DENV2 (New Guinea C (NGC)) was isolated in Papua, New Guinea (Accession No. AF038403)31. Because the early genomic records were sparse, our estimate of the most recent common ancestor is likely to represent an upper boundary to the date of DENV2 emergence. Higher order phylogenetic relationships of amino acid sequences clustered E proteins into specific clades that mapped to discrete geographical origins and previously described lineages (Fig. 1; Supplementary File 2) of (i) Asian I, (ii) Asian II, (iii) Cosmopolitan, (iv) American, (v) Asian/American, and (vi) sylvatic genotypes32, 33. The clustering of lineages indicated robust support (posterior probabilities >0.9) for the evolutionary divergence of these groups and the constraints of DENV2 E amino acid diversity based on geographical distributions34. For example, a sylvatic clade encompassing isolates primarily from West Africa and Malaysia clustered independently of all other lineages, and diverged much earlier (Fig. 1). Further, Asian/American strains principally originated from the Caribbean region and nearby Latin American countries, but also from geographically distant Southeast Asia (Fig. 1). While distinct strains formed location-specific clusters within the

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Asian/American clade32, the geographical diversity among DENV2 strains may be due to rapid spread of mosquito vectors by urbanization and increased global travel34.

High surface entropy of the DENV2 E antigen Because antibodies resulting from immune responses to DENV interact with amino acid residues that are exposed to solvent on the surface of E proteins35-37, it follows that conserved positions are less likely to drive viral evolution, while selection pressures imposed by host antibodies may lead to changes in residues that are permissive to mutations. We evaluated the evolutionary stability of solvent-exposed positions by using entropy (H) as a measure of amino acid variation for each E residue over time (Fig. 2a). Only 22% of residues were completely conserved across all E proteins of DENV2 (Fig. 2a), suggesting a great amount of tolerance to mutations within each of the three E protein domains (EDI, II, and III). Approximately 30% of E monomer residues (146 of 495 total residues) are exposed38 (Fig. 2b) on the surface of the mature DENV2 virus39 (PDB 3J27; Supplementary Fig. 1), including the conserved glycosylation sites at positions N67 and N153. The remaining side chains are buried in contact surfaces between monomers or adjacent units of symmetry (Supplementary Fig. 1). Further, amino acid residues positioned as buried contacts between subunits (Fig. 2b, c) are highly conserved (34% H > 0). In contrast, 44% of surface exposed residues displayed H > 0, indicating that variability is focused on solvent accessible residues of E proteins (Fig. 2b, c). Notably, DENV2 harbors the least evolutionary stability in E protein surfaces among serotypes, with 32 solvent-accessible residues that are highly variable compared to ≤ 23 for each of the other 3 serotypes (Supplementary Fig. 2).

ACS Paragon Plus Environment

Page 6 of 42

Page 7 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Antibody recognition of E proteins across 60 years of DENV2 evolution Modeling the cumulative distribution function of the exponentially distributed number of strains (Scl) and total entropy (Hsecl) per clade (Supplementary Fig. 3), thresholds for each determinant were established as the value that represented 90% of the population (Scl = 30 and Hsecl = 0.2). Sixteen clades that represented major branch points in the evolution of five DENV2 lineages were identified by using values of Scl 0.2 (Fig. 1; Supplementary Fig. 3). A small cluster of E proteins from seven American strains that represented the sixth DENV2 lineage did not exceed Hsecl > 0.2, and thus did not meet a sufficient level of diversity. However, one of the sixteen selected strains harbored all amino acid substitutions present in the American sequences, except residue 71, which is A in Cosmopolitan strains, D in American, and E in all others. A representative isolate was selected from each of the sixteen DENV2 clades (Fig. 1; Table 1), and the variable E surface residues of the isolates (designated here as variants) were grafted onto the backbone of the PUO-218 E protein by designing synthetic genes of the E soluble ectodomain. The DENV2 proteins were produced by in vitro expression (Supplementary Fig. 4), including all essential components for protein glycosylation and translation from DNA templates. To evaluate the contribution of buried residues, the E variants 2C and 5C were synthesized with surfaces of PUO-218 on backbones that exhibited the lowest (variant 2) and highest (variant 5) total sequence identity to PUO-218 (Table 1). The recombinant proteins from E variants, the prototype strains of DENV1-4, and PUO-218 were isolated by affinity chromatography, and printed in microarrays for analysis of antibody recognition, using previously described methods40-43. Residues Gly104, Gly106, and Leu107 of the fusion-loop peptide of E domain II are highly conserved among the flaviviruses, and site-specific mutations of these residues prevent binding of monoclonal antibody 4G2 to bind to the E glycoprotein monomer44, 45. The

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

printed E proteins were equally recognized by 4G2 and a rabbit polyclonal antibody that was raised against the fusion-loop peptide that harbors the 4G2 epitope (Supplementary Fig. 4). We reasoned that an examination of polyclonal B-cell responses to monoclonal E antigens should allow us to describe a global perspective of adaptive immune responses to dengue that may complement results obtained with monoclonal antibodies35-37. Serological immune responses were studied with a non-human primate model that is used to evaluate efficacy of dengue vaccines40, 46. Asian rhesus macaques (Macaca mulatta) were subcutaneously challenged with representative strains of each DENV serotype, and blood sera were collected for analysis. The DENV2 challenge strain (S16803) was an Asian I lineage that differs from the full-length E polypeptide of PUO-218 by a single conservative substitution of aspartic for glutamic acid (PUO-218) at the surface accessible residue of position 383 (Supplementary File 1). The E protein of DENV4 was recognized ( 3.5), were removed and pixel counts from replicate spots were averaged to obtain mean fluorescence intensities (MFI). Relative binding was calculated with equation 2: (2)

Relative binding = ( x / xi )100,

where x is MFI origination from microarrayed antigens, and i is the infecting virus species. Percent signal change was calculated using equation 3: (3)

% signal change =

( ) × 100, 𝑦 ― 𝑦𝑗 𝑦𝑗

where y is the MFI originating from DENV2 E variant proteins and j is the E variant exhibiting the highest sequence identity to the infecting strain. Statistically significant antibody binding was

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

determined by using Student’s two-tailed t-tests, performed with the Data Analysis add-in for Microsoft Excel. Histograms and exponential curve fitting were performed using OriginPro v9.0 (Origin Lab Corporation, Northampton, MA). For threshold determinations of number of strains and entropy per clade, we selected the value that represented 90% of the distribution based on modeling the cumulative distribution function with equation 4: (4)

0.9 = 1- 𝑒( ― 𝜆𝑥),

where λ = 4.84 (Hsecl) and 0.06206 (Scl) and x is the selected threshold value.

DENV2 E proteins and neutralizing antibodies. A cell-based infection assay was used to evaluate the relationship between virus-neutralizing antibodies and the recombinant E protein variants. To produce virus for infection assays, Vero cells (ATCC® CCL-81TM) grown at 37oC, 5% CO2 in Dulbecco's Modification of Eagle's Medium (DMEM) (Mediatech, Inc, Manassas, VA) containing 10% HycloneTM (GE Healthcare Life Sciences, Logan, UT) fetal bovine serum (FBS) were seeded in T-75 flasks 24 h before infection. DENV2 NGC (BEI Resources NR-84) was diluted in DMEM containing 2% FBS and Vero cells were infected with 0.14 TCID50 units per cell (~MOI of 0.1). Cell culture supernatants harvested 8 days post-infection were centrifuged (1200 x g, 10 min) and passed through a 0.2 µm HT Tuffryn® Acrodisc® filter (Pall Life Sciences, Ann Arbor, MI) prior to storage (-80oC). For infection assays, Vero cells (3.8 x 104/well) were seeded in 24-well plates using DMEM with 10% FBS 24 h before infection. Viral titration experiments were performed in infection media (DMEM containing 2% heatinactivated FBS) to determine the amount of virus (2.5 µL) to infect up to 15% of the cells. A previously described flow cytometry-based infectivity assay71 was used with modifications to examine neutralizing antibody interactions with recombinant variant E proteins. Sera from

ACS Paragon Plus Environment

Page 20 of 42

Page 21 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

human DENV2 challenge subjects were heat-inactivated (55oC, 30 min) and 50% neutralizing (NT50) antibody titers were calculated as the sera dilution factor required to neutralize 50% of the virus. To determine the inhibitory effect of soluble E proteins on antibody neutralization of cell infection, titrated dilutions of the purified recombinant E (2-200 nM) were incubated (37oC, 1 h) with diluted sera (2x NT50) from the DENV2 challenge subjects. Equal volumes of infection media containing no virus (Vero uninfected control) or virus diluted in infection media were incubated (37oC, 1 h) with the sera and E mixtures. Cell culture media were replaced with fresh media containing penicillin-streptomycin (1%), and the Vero cells were cultured for 1.5 h (37oC, 5% CO2,). The supernatants were aspirated and the cells were washed once before incubating (48 h, 37oC) with fresh infection media containing penicillin-streptomycin (1%). DENV2infected cells were washed (2 x) with PBS (Mediatech, Inc, Manassas, VA) and trypsinized to detach from the plastic surfaces. PBS containing 10% FBS was added, and the cells were resuspended and transferred to polystyrene tubes. The cells were centrifuged (500 x g, 5 min), washed in PBS (1 x), treated with 1x BD FACS™ Permeabilizing Solution 2 (22oC, 10 min), washed in PBS (1 x), blocked with 1x PBS/BSA (5%) (22oC, 20 min), and washed again (1 x PBS). The cells were incubated (22oC, 20 min) with anti-dengue virus complex antibody, clone D3-2H2-9-21 (MilliporeSigma, Burlington, MA) diluted (10 µg/mL) in PBS/1% BSA, and washed with PBS (1 x). To detect the primary antibody, the cells were incubated with Alexa Fluor 488-conjugated goat anti-mouse IgG (H+L) (Life Technologies) diluted 1:1000 in PBS (22oC, 20 min), centrifuged, and stored in PBS (4oC) containing 2% formaldehyde (ThermoFisher Scientific). Flow cytometry data were acquired on a BD FACSCalibur™ instrument, using BD CellQuestTM Pro software v 5.2.1, and analyzed using the Super-Enhanced Dmax Subtraction (SED) algorithm in FlowJo v10.3, with uninfected Vero cell histograms used

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 42

as controls in population comparisons. Infected cell populations were quantile normalized by group (pre- and post-challenge), and the percentage of virus neutralized by post-challenge antibodies was calculated as follows with equation 5: (5)

% neutralization =

(

𝑦𝑝𝑟𝑒



) × 100,

𝑦𝑝𝑜𝑠𝑡

𝑦𝑝𝑟𝑒

where y is the percentage of infected cells by virus pre-incubated with sera from each time point. Neutralizing antibody inhibition by recombinant E variants was calculated as the change in the % of virus neutralized, relative to the level of neutralization achieved in the presence of recombinant PUO-218 envelope protein.

Dengue challenge sera. Sixteen healthy, flavivirus naïve rhesus macaques (M. mulatta) were subcutaneously injected with 105 PFU of either DENV1 (West Pac 74), DENV2 (S16803), DENV3 (CH53489), or DENV4 (341750) (n=4 per challenge group), derived from low passage, near wild-type virus isolates. Sera were collected prior to and 30 days post-infection (DPI)40. Sera from human primary DENV2 infections (n = 10) were collected as part of a DENV human challenge model originally developed by the Laboratory of Infectious Diseases at the U. S. National Institutes of Health47, 48. Sera from participants that had no prior history or serological evidence of flavivirus infection were collected prior to and 28 days post-challenge with 103 PFU of rDEN2∆30. rDEN2∆30 induced viremia in all participants by 5 DPI.

Ethics Statement This research was conducted in compliance with the Animal Welfare Act and other federal statutes and regulations relating to animals and experiments involving animals, and adhered principles stated in the Guide for the Care and Use of Laboratory Animals, National Research

ACS Paragon Plus Environment

Page 23 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Council, 1996, under facilities fully accredited by the Association for the Assessment and Accreditation of Laboratory Animal Care International, by an approved protocol through the Institutional Animal Care and Use Committee (IACUC) of the Walter Reed Army Institute of Research (WRAIR). Animal care procedures complied with the regulations detailed under the Animal Welfare Act and in the Guide for the Care and Use of Laboratory Animals, in accordance with the recommendations of the Weatherall report, “The use of non-human primates in research”. Animals were kept at a constant temperature (22 °C ± 2 °C) and relative humidity (50%), with environmental enrichment and 12 h of artificial light per day. If the animal was in pain on the basis of their observations of animal behavior, analgesics were subsequently administered, through a single i.m. injection of 5 mg/kg flunixine (Finadyne®, Schering Plough) in the days after interventions. Research on human subjects was conducted in full compliance with DoD, NIH, federal, and state statutes and regulations relating to the protection of human subjects and adheres to principles identified in the Belmont Report (1979). Sera from adult human (18-50 years of age) primary DENV2 (rDEN2∆30) challenges were collected under approval from the Institutional Review Boards of the National Institutes of Health (NIH), National Institute of Allergy and Infectious Diseases (NIAID), Regulatory Compliance and Human Subjects Protection Branch (RCHSPB), and also the Western Institutional Review Board® (Olympia, Washington). All subjects participated willingly in the study as evidenced by signing the informed consent document.

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Ancillary Information Supporting Information: Supplemental Figures S1-S5 and Supplemental Files 1 and 2.

Corresponding Author Information: [email protected]; 301 619-4232

Present/Current Author Addresses: 1Molecular and Translational Sciences Division, U.S. Army Medical Research Institute of Infectious Diseases, Frederick, MD 21702; 2Department of Biological Sciences, University of Maryland Baltimore County, Baltimore, MD 21250; 3Armed

Forces Research Institute of Medical Sciences, Bangkok, 10400, Thailand; 4Center for

Immunization Research, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland 21205, USA

Author Contributions: SLK and JLS contributed equally as first authors. Conceived study (RGU, SF); prepared and reviewed manuscript (SLK, JLS, SF, APD, BMZ, RGU); performed laboratory experiments (SLK, JLS); protein modeling (SLK, BMZ); data analysis (SLK, JLS, SF, APD, BMZ, RGU).

Acknowledgment: Funding was provided to RGU and SF by the Military Infectious Disease Research Program (https://midrp.amedd.army.mil/), project S0453_14_AF. This research was supported in part by the National Academies of Sciences, Engineering, and Medicine through a National Research Council (NRC) research associateship awarded to JLS, and by an appointment of SLK and BZ to the Postgraduate Research Participation Program administered by Oak Ridge Institute for Science and Education through an interagency agreement with the U.S. Department

ACS Paragon Plus Environment

Page 24 of 42

Page 25 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

of Energy. Opinions, interpretations, conclusions, and recommendations are those of the authors and are not necessarily endorsed by any Agency or branch of the U.S. government.

Abbreviations Used: ADE, antibody-dependent enhancement; CYD-TDV, recombinant chimeric yellow fever-17D-dengue virus, live attenuated, tetravalent dengue vaccine; DENV1-4, dengue virus serotypes 1-4; DENV2-NGC, dengue virus serotype 2 strain New Guinea C; E, envelope protein; EDI, EDII, and EDIII, envelope protein domains I, II, and III; ESS, effective sample size; HPD, highest posterior density; Hsecl, total entropy per clade; IMAC, immobilized metal affinity chromatography; IVT, in vitro translation; MCC, maximum clade credibility tree; MCMC, Markov Chain Monte Carlo; MFI, mean fluorescence intensities; NHP, nonhuman primate; NT50, 50% neutralizing antibody titers; PFU, plaque forming unit; prM, pre-membrane; PRNT, virus plaque reduction neutralization titer; PUO-218, dengue virus serotype 2 vaccine strain; Scl, distribution function of the exponentially distributed number of strains; TVD, tetravalent dengue vaccines; TM, transmembrane; ViPR, Virus Pathogen Database and Analysis Resource

References [1] Pang, T., Cardosa, M. J., and Guzman, M. G. (2007) Of cascades and perfect storms: the immunopathogenesis of dengue haemorrhagic fever-dengue shock syndrome (DHF/DSS), Immunology and cell biology 85, 43-45. DOI: 10.1038/sj.icb.7100008. [2] Gubler, D. J. (2012) The economic burden of dengue, The American journal of tropical medicine and hygiene 86, 743-744. DOI: 10.4269/ajtmh.2012.12-0157. [3] Bhatt, S., Gething, P. W., Brady, O. J., Messina, J. P., Farlow, A. W., Moyes, C. L., Drake, J. M., Brownstein, J. S., Hoen, A. G., Sankoh, O., Myers, M. F., George, D. B., Jaenisch, T., Wint, G. R., Simmons, C. P., Scott, T. W., Farrar, J. J., and Hay, S. I. (2013) The global distribution and burden of dengue, Nature 496, 504-507. DOI: 10.1038/nature12060. [4] Rothman, A. L. (2004) Dengue: defining protective versus pathologic immunity, The Journal of clinical investigation 113, 946-951. DOI: 10.1172/JCI21512.

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[5] Lanciotti, R. S., Gubler, D. J., and Trent, D. W. (1997) Molecular evolution and phylogeny of dengue-4 viruses, The Journal of general virology 78 ( Pt 9), 2279-2284. DOI: 10.1099/0022-1317-78-9-2279. [6] Shrestha, B., Brien, J. D., Sukupolvi-Petty, S., Austin, S. K., Edeling, M. A., Kim, T., O'Brien, K. M., Nelson, C. A., Johnson, S., Fremont, D. H., and Diamond, M. S. (2010) The development of therapeutic antibodies that neutralize homologous and heterologous genotypes of dengue virus type 1, PLoS pathogens 6, e1000823. DOI: 10.1371/journal.ppat.1000823. [7] Wahala, W. M., Donaldson, E. F., de Alwis, R., Accavitti-Loper, M. A., Baric, R. S., and de Silva, A. M. (2010) Natural strain variation and antibody neutralization of dengue serotype 3 viruses, PLoS pathogens 6, e1000821. DOI: 10.1371/journal.ppat.1000821. [8] Brien, J. D., Austin, S. K., Sukupolvi-Petty, S., O'Brien, K. M., Johnson, S., Fremont, D. H., and Diamond, M. S. (2010) Genotype-specific neutralization and protection by antibodies against dengue virus type 3, Journal of virology 84, 10630-10643. DOI: 10.1128/JVI.01190-10. [9] Katzelnick, L. C., Fonville, J. M., Gromowski, G. D., Arriaga, J. B., Green, A., James, S. L., Lau, L., Montoya, M., Wang, C., VanBlargan, L. A., Russell, C. A., Thu, H. M., Pierson, T. C., Buchy, P., Aaskov, J. G., Munoz-Jordan, J. L., Vasilakis, N., Gibbons, R. V., Tesh, R. B., Osterhaus, A. D., Fouchier, R. A., Durbin, A., Simmons, C. P., Holmes, E. C., Harris, E., Whitehead, S. S., and Smith, D. J. (2015) Dengue viruses cluster antigenically but not as discrete serotypes, Science 349, 1338-1343. DOI: 10.1126/science.aac5017. [10] Halstead, S. B., and O'Rourke, E. J. (1977) Dengue viruses and mononuclear phagocytes. I. Infection enhancement by non-neutralizing antibody, The Journal of experimental medicine 146, 201-217. DOI: 10.1084/jem.146.1.201. [11] Sukupolvi-Petty, S., Austin, S. K., Purtha, W. E., Oliphant, T., Nybakken, G. E., Schlesinger, J. J., Roehrig, J. T., Gromowski, G. D., Barrett, A. D., Fremont, D. H., and Diamond, M. S. (2007) Type- and subcomplex-specific neutralizing antibodies against domain III of dengue virus type 2 envelope protein recognize adjacent epitopes, Journal of virology 81, 12816-12826. DOI: 10.1128/JVI.00432-07. [12] Holmes, E. C., and Twiddy, S. S. (2003) The origin, emergence and evolutionary genetics of dengue virus, Infection, genetics and evolution 3, 19-28. DOI: 10.1016/S15671348(03)00004-2. [13] Twiddy, S. S., Farrar, J. J., Vinh Chau, N., Wills, B., Gould, E. A., Gritsun, T., Lloyd, G., and Holmes, E. C. (2002) Phylogenetic relationships and differential selection pressures among genotypes of dengue-2 virus, Virology 298, 63-72. DOI: 10.1006/viro.2002.1447. [14] Anez, G., Morales-Betoulle, M. E., and Rios, M. (2011) Circulation of different lineages of dengue virus type 2 in Central America, their evolutionary time-scale and selection pressure analysis, PloS one 6, e27459. DOI: 10.1371/journal.pone.0027459. [15] OhAinle, M., Balmaseda, A., Macalalad, A. R., Tellez, Y., Zody, M. C., Saborio, S., Nunez, A., Lennon, N. J., Birren, B. W., Gordon, A., Henn, M. R., and Harris, E. (2011) Dynamics of dengue disease severity determined by the interplay between viral genetics and serotype-specific immunity, Science translational medicine 3, 114ra128. DOI: 10.1126/scitranslmed.3003084. [16] de Alwis, R., Smith, S. A., Olivarez, N. P., Messer, W. B., Huynh, J. P., Wahala, W. M., White, L. J., Diamond, M. S., Baric, R. S., Crowe, J. E., Jr., and de Silva, A. M. (2012) Identification of human neutralizing antibodies that bind to complex epitopes on dengue

ACS Paragon Plus Environment

Page 26 of 42

Page 27 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

virions, Proceedings of the National Academy of Sciences of the United States of America 109, 7439-7444. DOI: 10.1073/pnas.1200566109. [17] Chao, D. Y., King, C. C., Wang, W. K., Chen, W. J., Wu, H. L., and Chang, G. J. (2005) Strategically examining the full-genome of dengue virus type 3 in clinical isolates reveals its mutation spectra, Virology journal 2, 72. DOI: 10.1186/1743-422X-2-72. [18] Modis, Y., Ogata, S., Clements, D., and Harrison, S. C. (2003) A ligand-binding pocket in the dengue virus envelope glycoprotein, Proceedings of the National Academy of Sciences of the United States of America 100, 6986-6991. DOI: 10.1073/pnas.0832193100. [19] Modis, Y., Ogata, S., Clements, D., and Harrison, S. C. (2004) Structure of the dengue virus envelope protein after membrane fusion, Nature 427, 313-319. DOI: 10.1038/nature02165. [20] Zhang, W., Chipman, P. R., Corver, J., Johnson, P. R., Zhang, Y., Mukhopadhyay, S., Baker, T. S., Strauss, J. H., Rossmann, M. G., and Kuhn, R. J. (2003) Visualization of membrane protein domains by cryo-electron microscopy of dengue virus, Nature structural biology 10, 907-912. DOI: 10.1038/nsb990. [21] Durbin, A. P., and Whitehead, S. S. (2010) Dengue vaccine candidates in development, Current topics in microbiology and immunology 338, 129-143. DOI: 10.1007/978-3-64202215-9_10. [22] Sabchareon, A., Wallace, D., Sirivichayakul, C., Limkittikul, K., Chanthavanich, P., Suvannadabba, S., Jiwariyavej, V., Dulyachai, W., Pengsaa, K., Wartel, T. A., Moureau, A., Saville, M., Bouckenooghe, A., Viviani, S., Tornieporth, N. G., and Lang, J. (2012) Protective efficacy of the recombinant, live-attenuated, CYD tetravalent dengue vaccine in Thai schoolchildren: a randomised, controlled phase 2b trial, Lancet. DOI: 10.1016/S0140-6736(12)61428-7. [23] Villar, L., Dayan, G. H., Arredondo-Garcia, J. L., Rivera, D. M., Cunha, R., Deseda, C., Reynales, H., Costa, M. S., Morales-Ramirez, J. O., Carrasquilla, G., Rey, L. C., Dietze, R., Luz, K., Rivas, E., Miranda Montoya, M. C., Cortes Supelano, M., Zambrano, B., Langevin, E., Boaz, M., Tornieporth, N., Saville, M., Noriega, F., and Group, C. Y. D. S. (2015) Efficacy of a tetravalent dengue vaccine in children in Latin America, The New England journal of medicine 372, 113-123. DOI: 10.1056/NEJMoa1411037. [24] Capeding, M. R., Tran, N. H., Hadinegoro, S. R., Ismail, H. I., Chotpitayasunondh, T., Chua, M. N., Luong, C. Q., Rusmil, K., Wirawan, D. N., Nallusamy, R., Pitisuttithum, P., Thisyakorn, U., Yoon, I. K., van der Vliet, D., Langevin, E., Laot, T., Hutagalung, Y., Frago, C., Boaz, M., Wartel, T. A., Tornieporth, N. G., Saville, M., Bouckenooghe, A., and Group, C. Y. D. S. (2014) Clinical efficacy and safety of a novel tetravalent dengue vaccine in healthy children in Asia: a phase 3, randomised, observer-masked, placebocontrolled trial, Lancet 384, 1358-1365. DOI: 10.1016/S0140-6736(14)61060-6. [25] Rico-Hesse, R. (1990) Molecular evolution and distribution of dengue viruses type 1 and 2 in nature, Virology 174, 479-493. DOI: 10.1016/0042-6822(90)90102-W. [26] Drummond, A. J., Rambaut A, Xie W. (2002-2010) BEAUTI-Bayesian evolutionary analysis utility, 2013. URL: http://beast.bio.ed.ac.uk. [27] Drummond, A. J., Suchard, M. A., Xie, D., and Rambaut, A. (2012) Bayesian phylogenetics with BEAUti and the BEAST 1.7, Molecular biology and evolution 29, 1969-1973. DOI: 10.1093/molbev/mss075.

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[28] Rambaut, A. (2006-2012) Tree Figure Drawing Tool version 1.4, 2013. URL: http://tree.bio.ed.ac.uk/software/figtree/. [29] Rambaut, A., Drummond A. J. (2003-2009) Tracer: MCMC Trace Analysis Tool v1.5.0, 2013. URL: http://beat.bio.ed.ac.uk/. [30] Rambaut, A., Drummond A. J. (2002-2010) TreeAnnotator v1.7.5, 2013. URL: http://beat.bio.ed.ac.uk/. [31] Gruenberg, A., Woo, W. S., Biedrzycka, A., and Wright, P. J. (1988) Partial nucleotide sequence and deduced amino acid sequence of the structural proteins of dengue virus type 2, New Guinea C and PUO-218 strains, The Journal of general virology 69 ( Pt 6), 1391-1398. DOI: 10.1099/0022-1317-69-6-1391. [32] Twiddy, S. S., Woelk, C. H., and Holmes, E. C. (2002) Phylogenetic evidence for adaptive evolution of dengue viruses in nature, The Journal of general virology 83, 1679-1689. DOI: 10.1099/0022-1317-83-7-1679. [33] Salda, L. T., Parquet, M. D., Matias, R. R., Natividad, F. F., Kobayashi, N., and Morita, K. (2005) Molecular epidemiology of dengue 2 viruses in the Philippines: genotype shift and local evolution, The American journal of tropical medicine and hygiene 73, 796-802. DOI: 10.4269/ajtmh.2005.73.796. [34] Sim, S., and Hibberd, M. L. (2016) Genomic approaches for understanding dengue: insights from the virus, vector, and host, Genome biology 17, 38. DOI: 10.1186/s13059-0160907-2. [35] Austin, S. K., Dowd, K. A., Shrestha, B., Nelson, C. A., Edeling, M. A., Johnson, S., Pierson, T. C., Diamond, M. S., and Fremont, D. H. (2012) Structural basis of differential neutralization of DENV-1 genotypes by an antibody that recognizes a cryptic epitope, PLoS pathogens 8, e1002930. DOI: 10.1371/journal.ppat.1002930. [36] Cockburn, J. J., Navarro Sanchez, M. E., Goncalvez, A. P., Zaitseva, E., Stura, E. A., Kikuti, C. M., Duquerroy, S., Dussart, P., Chernomordik, L. V., Lai, C. J., and Rey, F. A. (2012) Structural insights into the neutralization mechanism of a higher primate antibody against dengue virus, The EMBO journal 31, 767-779. DOI: 10.1038/emboj.2011.439. [37] Fibriansah, G., Tan, J. L., Smith, S. A., de Alwis, A. R., Ng, T. S., Kostyuchenko, V. A., Ibarra, K. D., Wang, J., Harris, E., de Silva, A., Crowe, J. E., Jr., and Lok, S. M. (2014) A potent anti-dengue human antibody preferentially recognizes the conformation of E protein monomers assembled on the virus surface, EMBO molecular medicine 6, 358371. DOI: 10.1002/emmm.201303404. [38] Krissinel, E., and Henrick, K. (2007) Inference of macromolecular assemblies from crystalline state, Journal of molecular biology 372, 774-797. DOI: 10.1016/j.jmb.2007.05.022. [39] Zhang, X., Ge, P., Yu, X., Brannan, J. M., Bi, G., Zhang, Q., Schein, S., and Zhou, Z. H. (2013) Cryo-EM structure of the mature dengue virus at 3.5-A resolution, Nature structural & molecular biology 20, 105-110. DOI: 10.1038/nsmb.2463. [40] Fernandez, S., Cisney, E. D., Tikhonov, A. P., Schweitzer, B., Putnak, R. J., Simmons, M., and Ulrich, R. G. (2011) Antibody recognition of the dengue virus proteome and implications for development of vaccines, Clinical and vaccine immunology : CVI 18, 523-532. DOI: 10.1128/CVI.00016-11. [41] Kamata, T., Natesan, M., Warfield, K., Aman, M. J., and Ulrich, R. G. (2014) Determination of specific antibody responses to the six species of ebola and Marburg

ACS Paragon Plus Environment

Page 28 of 42

Page 29 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

viruses by multiplexed protein microarrays, Clinical and vaccine immunology : CVI 21, 1605-1612. DOI: 10.1128/CVI.00484-14. [42] Natesan, M., Jensen, S. M., Keasey, S. L., Kamata, T., Kuehne, A. I., Stonier, S. W., Lutwama, J. J., Lobel, L., Dye, J. M., and Ulrich, R. G. (2016) Human Survivors of Disease Outbreaks Caused by Ebola or Marburg Virus Exhibit Cross-Reactive and LongLived Antibody Responses, Clinical and vaccine immunology : CVI 23, 717-724. DOI: 10.1128/CVI.00107-16. [43] Keasey, S., Pugh, C., Tikhonov, A., Chen, G., Schweitzer, B., Nalca, A., and Ulrich, R. G. (2010) Proteomic basis of the antibody response to monkeypox virus infection examined in cynomolgus macaques and a comparison to human smallpox vaccination, PloS one 5, e15547. DOI: 10.1371/journal.pone.0015547. [44] Crill, W. D., and Chang, G. J. (2004) Localization and characterization of flavivirus envelope glycoprotein cross-reactive epitopes, Journal of virology 78, 13975-13986. DOI: 10.1128/JVI.78.24.13975-13986.2004. [45] Metz, S. W., Gallichotte, E. N., Brackbill, A., Premkumar, L., Miley, M. J., Baric, R., and de Silva, A. M. (2017) In Vitro Assembly and Stabilization of Dengue and Zika Virus Envelope Protein Homo-Dimers, Sci Rep 7, 4524. DOI: 10.1038/s41598-017-04767-6. [46] Onlamoon, N., Noisakran, S., Hsiao, H. M., Duncan, A., Villinger, F., Ansari, A. A., and Perng, G. C. (2010) Dengue virus-induced hemorrhage in a nonhuman primate model, Blood 115, 1823-1834. DOI: 10.1182/blood-2009-09-242990. [47] Larsen, C. P., Whitehead, S. S., and Durbin, A. P. (2015) Dengue human infection models to advance dengue vaccine development, Vaccine 33, 7075-7082. DOI: 10.1016/j.vaccine.2015.09.052. [48] Kirkpatrick, B. D., Whitehead, S. S., Pierce, K. K., Tibery, C. M., Grier, P. L., Hynes, N. A., Larsson, C. J., Sabundayo, B. P., Talaat, K. R., Janiak, A., Carmolli, M. P., Luke, C. J., Diehl, S. A., and Durbin, A. P. (2016) The live attenuated dengue vaccine TV003 elicits complete protection against dengue in a human challenge model, Science translational medicine 8, 330ra336. DOI: 10.1126/scitranslmed.aaf1517. [49] Butrapet, S., Childers, T., Moss, K. J., Erb, S. M., Luy, B. E., Calvert, A. E., Blair, C. D., Roehrig, J. T., and Huang, C. Y. (2011) Amino acid changes within the E protein hinge region that affect dengue virus type 2 infectivity and fusion, Virology 413, 118-127. DOI: 10.1016/j.virol.2011.01.030. [50] Rogers, D. M., Kent, M. S., and Rempe, S. B. (2015) Molecular basis of endosomalmembrane association for the dengue virus envelope protein, Biochimica et biophysica acta 1848, 1041-1052. DOI: 10.1016/j.bbamem.2014.12.018. [51] Guirakhoo, F., Arroyo, J., Pugachev, K. V., Miller, C., Zhang, Z. X., Weltzin, R., Georgakopoulos, K., Catalan, J., Ocran, S., Soike, K., Ratterree, M., and Monath, T. P. (2001) Construction, safety, and immunogenicity in nonhuman primates of a chimeric yellow fever-dengue virus tetravalent vaccine, Journal of virology 75, 7290-7304. DOI: 10.1128/JVI.75.16.7290-7304.2001. [52] Sin Leo, Y., Wilder-Smith, A., Archuleta, S., Shek, L., Chong, C. Y., Nam Leong, H., Yong Low, C., Oh, M. L., Bouckenooghe, A., Crevat, D., and Wartel, A. (2012) Immunogenicity and safety of recombinant tetravalent dengue vaccine (CYD-TDV) in individuals aged 2-45 y: Phase II randomized controlled trial in Singapore, Hum Vaccin Immunother 8. DOI: 10.4161/hv.21224.

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[53] Capeding, R. Z., Luna, I. A., Bomasang, E., Lupisan, S., Lang, J., Forrat, R., Wartel, A., and Crevat, D. (2011) Live-attenuated, tetravalent dengue vaccine in children, adolescents and adults in a dengue endemic country: randomized controlled phase I trial in the Philippines, Vaccine 29, 3863-3872. DOI: 10.1016/j.vaccine.2011.03.057. [54] Poo, J., Galan, F., Forrat, R., Zambrano, B., Lang, J., and Dayan, G. H. (2010) Liveattenuated Tetravalent Dengue Vaccine in Dengue-naive Children, Adolescents, and Adults in Mexico City: Randomized Controlled Phase 1 Trial of Safety and Immunogenicity, Pediatr Infect Dis J. DOI: 10.1097/INF.0b013e3181fe05af. [55] Guy, B., Briand, O., Lang, J., Saville, M., and Jackson, N. (2015) Development of the Sanofi Pasteur tetravalent dengue vaccine: One more step forward, Vaccine 33, 71007111. DOI: 10.1016/j.vaccine.2015.09.108. [56] Guirakhoo, F., Weltzin, R., Chambers, T. J., Zhang, Z. X., Soike, K., Ratterree, M., Arroyo, J., Georgakopoulos, K., Catalan, J., and Monath, T. P. (2000) Recombinant chimeric yellow fever-dengue type 2 virus is immunogenic and protective in nonhuman primates, Journal of virology 74, 5477-5485. DOI: 10.1128/JVI.74.12.5477-5485.2000. [57] Forshey, B. M., Reiner, R. C., Olkowski, S., Morrison, A. C., Espinoza, A., Long, K. C., Vilcarromero, S., Casanova, W., Wearing, H. J., Halsey, E. S., Kochel, T. J., Scott, T. W., and Stoddard, S. T. (2016) Incomplete Protection against Dengue Virus Type 2 Reinfection in Peru, PLoS neglected tropical diseases 10, e0004398. DOI: 10.1371/journal.pntd.0004398. [58] Parameswaran, P., Charlebois, P., Tellez, Y., Nunez, A., Ryan, E. M., Malboeuf, C. M., Levin, J. Z., Lennon, N. J., Balmaseda, A., Harris, E., and Henn, M. R. (2012) Genomewide patterns of intrahuman dengue virus diversity reveal associations with viral phylogenetic clade and interhost diversity, Journal of virology 86, 8546-8558. DOI: 10.1128/JVI.00736-12. [59] Romano, C. M., Lauck, M., Salvador, F. S., Lima, C. R., Villas-Boas, L. S., Araujo, E. S., Levi, J. E., Pannuti, C. S., O'Connor, D., and Kallas, E. G. (2013) Inter- and intra-host viral diversity in a large seasonal DENV2 outbreak, PloS one 8, e70318. DOI: 10.1371/journal.pone.0070318. [60] Srikiatkhachorn, A., and Yoon, I. K. (2016) Immune correlates for dengue vaccine development, Expert Rev Vaccines 15, 455-465. DOI: 10.1586/14760584.2016.1116949. [61] Barban, V., Munoz-Jordan, J. L., Santiago, G. A., Mantel, N., Girerd, Y., Gulia, S., Claude, J. B., and Lang, J. (2012) Broad neutralization of wild-type dengue virus isolates following immunization in monkeys with a tetravalent dengue vaccine based on chimeric yellow fever 17D/dengue viruses, Virology 429, 91-98. DOI: 10.1016/j.virol.2012.03.007. [62] Wang, C., Katzelnick, L. C., Montoya, M., Hue, K. D., Simmons, C. P., and Harris, E. (2016) Evolutionarily Successful Asian 1 Dengue Virus 2 Lineages Contain One Substitution in Envelope That Increases Sensitivity to Polyclonal Antibody Neutralization, The Journal of infectious diseases 213, 975-984. DOI: 10.1093/infdis/jiv536. [63] VanBlargan, L. A., Mukherjee, S., Dowd, K. A., Durbin, A. P., Whitehead, S. S., and Pierson, T. C. (2013) The type-specific neutralizing antibody response elicited by a dengue vaccine candidate is focused on two amino acids of the envelope protein, PLoS pathogens 9, e1003761. DOI: 10.1371/journal.ppat.1003761.

ACS Paragon Plus Environment

Page 30 of 42

Page 31 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

[64] Keasey, S. L., Pugh, C. L., Jensen, S. M., Smith, J. L., Hontz, R. D., Durbin, A. P., Dudley, D. M., O'Connor, D. H., and Ulrich, R. G. (2017) Antibody Responses to Zika Virus Infections in Environments of Flavivirus Endemicity, Clinical and vaccine immunology : CVI 24. DOI: 10.1128/CVI.00036-17. [65] World Health Organization. (2016) Dengue vaccine: WHO position paper - July 2016, Releve epidemiologique hebdomadaire 91, 349-364. URL: http://www.who.int/wer/2016/wer9130.pdf?ua=1. [66] Pickett, B. E., Sadat, E. L., Zhang, Y., Noronha, J. M., Squires, R. B., Hunt, V., Liu, M., Kumar, S., Zaremba, S., Gu, Z., Zhou, L., Larson, C. N., Dietrich, J., Klem, E. B., and Scheuermann, R. H. (2012) ViPR: an open bioinformatics database and analysis resource for virology research, Nucleic acids research 40, D593-598. DOI: 10.1093/nar/gkr859. [67] Abascal, F., Zardoya, R., and Posada, D. (2005) ProtTest: selection of best-fit models of protein evolution, Bioinformatics 21, 2104-2105. DOI: 10.1093/bioinformatics/bti263. [68] Sievers, F., Wilm, A., Dineen, D., Gibson, T. J., Karplus, K., Li, W., Lopez, R., McWilliam, H., Remmert, M., Soding, J., Thompson, J. D., and Higgins, D. G. (2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega, Molecular systems biology 7, 539. DOI: 10.1038/msb.2011.75. [69] Goujon, M., McWilliam, H., Li, W., Valentin, F., Squizzato, S., Paern, J., and Lopez, R. (2010) A new bioinformatics analysis tools framework at EMBL-EBI, Nucleic acids research 38, W695-699. DOI: 10.1093/nar/gkq313. [70] Hall, T. A. (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT, Nucleic Acids Symposium Series 41, 95-98. [71] de Alwis, R., and de Silva, A. M. (2014) Measuring antibody neutralization of dengue virus (DENV) using a flow cytometry-based technique, Methods Mol Biol 1138, 27-39. DOI: 10.1007/978-1-4939-0348-1_3.

Figure Legends Figure 1. Molecular phylogeny of DENV2 envelope protein sequences. The evolutionary history of 266 unique global DENV2 envelope (E) protein sequences isolated between 1944 and 2010 was inferred using the maximum clade credibility tree obtained using a Bayesian method26, 27.

Clade branches are colored according to geographical area of isolation, which correspond to

the colored country regions in the continent maps. For clarity, triangles represent large clades of geographically related isolates. Triangle size is proportional to the number of isolates contained within the clade. The horizontal length of the triangle reflects the age of the clade, and begins at

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the estimated MRCA of the clade and extends to the date of isolation of the most recently isolated strain. Probability values >0.9 are represented by (*), >0.95 by (**), and >0.99 by (***), inside the nodes. Each of the six DENV2 genotypes are represented in the tree and cluster distinctly, as indicated by the bold black vertical bars to the right of the tree. Numbered clades represent major branch points in the evolution of DENV2 E that exhibit elevated diversity in solvent accessible E amino acids displayed in macromolecular complexes on the virion surface.

Figure 2. Solvent accessible residues of the envelope monomer displayed on the surface of DENV2 exhibit increased diversity in amino acid sequences. Surface exposed envelope (E) residues that exhibited amino acid variability in a multiple sequence alignment of 266 unique E amino acid sequences representing isolates originating between 1944 and 2010 from diverse geographic regions and genotypes were mapped on to E protein structures. (a) Amino acid variability was calculated as Shannon entropy (H) for each column of the E sequence alignment. E residues are colored by E protein domain (red, DI; yellow, DII; royal blue, DIII; cyan, transmembrane (TM) domain). (b) Proportions of total and variable amino acids of the fulllength E polypeptide that are solvent accessible or non-solvent accessible. (c) Entropy values shown in a for each E amino acid were mapped on to an E monomer (space filling) in the structural model of the DENV2 E-M-M-E heterotetramer (PDB 3J2P) viewed from the exterior of the virus (top), side (middle), and interior (bottom). Residue entropy is gradient colored from low (white) to high (non-solvent accessible, red; solvent accessible, blue) according to the key shown at the right. For clarity, only ribbon structures of E DI-III (colored as in a) are shown in exterior and interior views, while E-TM domain (cyan) and M (orange) are included in the side view of the complex.

ACS Paragon Plus Environment

Page 32 of 42

Page 33 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Figure 3. Specificity of antibodies from primary challenges with DENV. Relative binding of convalescent serum antibodies from rhesus macaques and human primary DENV infections to microarrayed DENV envelope (E) proteins (DENV1-Hawaii, DENV2-NGC, DENV3-H87, DENV4-H241). Antibody specificity for E from representative DENV1-4 strains by sera from (a) rhesus macaques challenged independently with DENV1, DENV3, DENV4, or (b) rhesus macaques challenged with Asian I lineage and humans challenged with American lineage DENV2. Values are shown for antibody binding to E relative to the challenge serotype.

Figure 4. Polyclonal antibody recognition of solvent-accessible residues from current and historical DENV2 strains displayed on the PUO-218 envelope backbone. Serum antibodies from rhesus macaques challenged with the Asian-lineage strain S16803 (a, b) or humans challenged with the American-lineage strain Tonga/1974 (c, d). DENV2 lineages of residues displayed: Asian I, AsI; sylvatic, Syl; Asian II, AsII; Cosmopolitan, Cosmo; Asian/American, As/Am. In a and c, antibody binding is shown relative to the envelope (E) variant exhibiting the highest sequence identity to the strain used for challenge: V (PUO-218) for Asian-lineage antibodies; average of variants 12 and 13 for American-lineage antibodies. Antibody interactions (a-e) are shown as decreased (hashed), unchanged (open), or increased (shaded) bars. E proteins that exhibited changes in antibody recognition are indicated (* p < 0.05, ** p < 0.001, or *** p < 0.00001). E surface residues that contributed to differential antibody binding are indicated in b and d. Columns are labeled with residues of the challenge strains, and amino acid substitutions representing each DENV2 clade are indicated within rows. Amino acids are colored by charge (blue, positive; red, negative; unshaded, uncharged). (e) Antibody recognition of variants 2C and

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

5C that display surface residues of PUO-218 on backbones that exhibited the lowest (variant 2) and highest (variant 5) total sequence identity to PUO-218. Non-solvent accessible residues Y59F and T236M of 2C, as well as G228E of 5C, alter antibody binding compared to complementary variants 2 and 5.

Figure 5. Antibodies primarily target solvent-accessible residues within a 480Å2 region on the surface of the DENV2 envelope protein monomer. Surface exposed residues of the E monomer that contributed to differential recognition of envelope (E) protein variants by convalescent antibodies from non-human primates and humans after primary infection with DENV2 were mapped on to E protein structures. E monomers are colored by domain (EDI, red; EDII, yellow; EDIII, blue). (a) The soluble ectodomain structure of the DENV2 E dimer as viewed from the outside of the virus. Solvent accessible residues that were targeted by polyclonal antibodies are shown as green spheres (ribbon) or are green surfaces (space filling monomer), whereas buried residues that presumably affect E surfaces are shown in orange (ribbon). The residues form a hexagonal footprint (dashed line and inset) of ~480Å2. (b) An E protein raft consisting of three E dimers (PDB 3J27). Surface accessible residues associated with polyclonal antibody interactions are shown as green spheres. (c) Potential antibody recognition orientations based on surface accessible residues identified in our study: binding that spans surfaces of a single E dimer (green shading) or bridges neighboring dimers (blue shading).

ACS Paragon Plus Environment

Page 34 of 42

Page 35 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Table 1. DENV2 sequences selected for construction of chimeric envelope proteins. Cladea/ Variant

Accession Date of No. Strain Location isolation Genotype AF038403 NGC Papua New Guinea 1944 Asian II V D00345d PUO-218 Thailand 1980 Asian I 1 EF105385 Dak Ar D20761 Senegal 1974 Sylvatic 2 EF105382 Dak Ar 2039 Burkina Faso 1980 Sylvatic 3 EF105390 Dak Ar 141070 Senegal 1999 Sylvatic 4 HQ891024 V5056 Taiwan 2008 Asian II 5 FJ024454 V1675 Vietnam 2007 Asian I 6 EU482676 V739 Vietnam 2006 Asian I 7 FJ196852 GD01/01 China 2001 Cosmopolitan 8 JN851128 0232Y06 Singapore 2006 Cosmopolitan 9 GU131843 V3502 Burkina Faso 1986 Cosmopolitan 10 EU569707 V1368 USA/PRe 1995 As/Amf 11 EU677144 V1427 USA/PR 1999 As/Am 12 EU482720 V588 USA/PR 2006 As/Am 13 EU482726 V595 USA/PR 2006 As/Am 14 FJ639822 V2262 Venezuela 2006 As/Am 15 FJ182012 V1597 Columbia 2005 As/Am 16 FJ639698 V2021 Cambodia 2002 As/Am aClade number corresponds to numbering in Figure 1. bHse corresponds to the total entropy of surface exposed residues per clade cl cSurface exposed dPUO-218 included in CYD-TDV ePuerto Rico fAsian/American

ACS Paragon Plus Environment

Hseclb

0.42434 0.34658 0.69315 0.42434 0.21217 0.22528 0.34658 0.34658 0.21217 0.34658 0.34658 0.20016 0.20016 0.23468 0.34658 0.34658

% Sequence Identity Total SEc Residues Residues 100 94.34 93.94 94.93 97.98 98.59 98.59 97.78 98.18 98.18 97.98 98.18 98.59 97.17 97.78 97.76 97.78

100 90.41 87.67 90.28 95.21 97.95 97.95 95.89 96.58 95.89 95.89 97.26 97.26 97.26 96.58 96.5 95.89

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1.

ACS Paragon Plus Environment

Page 36 of 42

Page 37 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Figure 2.

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3.

ACS Paragon Plus Environment

Page 38 of 42

Page 39 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Figure 4.

ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5.

ACS Paragon Plus Environment

Page 40 of 42

Page 41 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

For Table of Contents Use Only

ACS Paragon Plus Environment

ACS Infectious Diseases

Dengue virus 2 envelope surface variants

Protein microarray

Detect antibody bound to proteins

Y Y

Y Y

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41

Page 42 of 42

Digital antibody binding data

ACS Paragon Plus Environment