Molecular basis for dengue virus broad cross-neutralization by humanized monoclonal antibody 513

Dengue is a widespread viral disease with 3.6 billion people at risk worldwide. Humanized monoclonal antibody (mAb) 513, currently undergoing clinical trials in Singapore, targets an epitope on the envelope protein domain III exposed at the surface of the viral particle. This antibody potently neutralizes all four dengue virus serotypes in a humanized mouse model that recapitulates human dengue infection, without signs of antibody-mediated enhancement of the disease. The crystal structure of single-chain variable fragment (scFv) 513 bound to the envelope protein domain III from dengue virus serotype 4 was used as a template to explore the molecular origins of the broader cross-reactivity and increased in vivo potency of mAb 513, compared to the parent murine mAb 4E11, using molecular dynamics simulations and network analyses. These two methods are a powerful complement to existing structural and binding data and detail specific interactions that underpin the differential binding of the two antibodies. We found that a Glu at position H55 (GluH55) from the second Complementarity Determining Region of the Heavy chain (CDR-H2) which corresponds to Ala in 4E11, is a major contributor to the enhancement in the interactions of mAb 513 compared to 4E11. Importantly, we also validate the importance of GluH55 using site-directed mutagenesis followed by isothermal titration calorimetry measurements.


Results
Electrostatic surface variations in the epitopes across the four DENV serotypes. An alignment of the amino-acid sequences of DIII domains of DENV1-4 used in this study, spanning amino-acid residues 296-400 from the envelope E protein ectodomain, is shown in Fig. 1A. In agreement with phylogenetic analysis conducted on a large number of DIII sequences from both DENV epidemic and sylvatic strains 23 , DENV1 and DENV3 share the highest amino-acid sequence identity (60-70%) and are both closer to DENV2 than to DENV4. DENV4 apparently diverged earlier with amino-acid sequence identities of only 40-50% with DENV1-3 24 . Thus, as anticipated, the evolution and mutations observed among DENV serotypes correlate directly with the cross-neutralization efficacy of mAbs 4E11, 1A1D-2, 9D12 and 513, all of which bind an epitope centered on the A-strand of DIII 15,18,[25][26][27][28] . Crystallographic and NMR studies of the flavivirus E protein ectodomain and of isolated DIII domains demonstrated that despite significant sequence variations between the four serotypes, the DIII structure is well conserved: an IgC-like fold with seven β-strands labeled from A to G 15,18,[25][26][27][28] . The N-terminal region of DIII as well as three strand-connecting loops BC, DE and FG are the most exposed at the virion surface and engage in contacts with various neutralizing antibodies described so far 23 . Moreover, the molecular surface of the viral particle is not a static structure but "breathes", especially at the physiological temperature of 37 °C, and this dynamic behavior gives antibodies access to epitopes that are partially hidden at lower temperatures 29 . As observed in homology models and experimental crystal structures 14,15 , amino-acid substitutions across the dengue serotypes modulate both the shape of the epitope and its electrostatic properties (Fig. 1B), largely explaining the differences in mAbs 513 and 4E11 binding affinities. Substitutions at residues 361-362 projecting from the DE loop, A-strand (305-312) and G-strand (388-390) are the main contributors to these subtle variations in shape and charge distribution: DENV1, 2 and 3 have Lys at positions 361 and 388, while DENV4 has a polar residue (Thr) (Fig. 1C). Another noticeable alteration of the electrostatic surface is observed at position 307 of the epitope, with the loss of a net positive charge: Lys in DENV1-2 is replaced by Val in DENV3 and by Ser in DENV4. Likewise, position 323 is a positively charged Arg or Lys in DENV2-4 and is replaced by a polar Gln in DENV1. At position 325, a positively charged Lys is present in DENV1 and DENV4, while a Gln is found in Scientific RepoRtS | (2018) 8:8449 | DOI: 10.1038/s41598-018-26800-y DENV2 and a negatively charged Glu in DENV3. The sections below examine in detail the influence of such substitutions on contacts that enhance mAb 513 binding towards the various DENV DIII antigens.
Analysis of scFv513-DIII-DENV4 contacts from the crystal structure of the complex. The scFv513 recombinant protein was produced in E. coli by linking the variable heavy chain region (VH) of mAb 513 to the variable light chain region (VL). We obtained two crystal structures from the Protein Data Bank (with PDB access codes 5AAM and 5AAW) providing eight views of the complex formed between scFv513 and DIII-DENV4 that was briefly described earlier 16 . An overall view of the complex ( Fig. 2A) as well as selected portions of the interactions established between scFv513 and DIII from DENV4 are displayed in Fig. 2, highlighting changes introduced between 513 and 4E11. A close-up view of the Ser26Δ deletion of the elbow in the framework region preceding CDR-H1 is displayed in Fig. 2B. The extra salt bridge between Lys323 (DIII) and Glu55 (CDR-H2) is  Table 2). Residue numbers correspond to the full length DENV E protein.
Residues that belong to β-strands (labeled from A-G) are below orange arrows. Residues marked with red and black dots under the alignment, represent conserved and variable epitope residues respectively. Colored blocks depicting the average surface accessibility of residues in white (accessible) to green (buried) are displayed under the alignment. (B) and (C) Structural comparison of the epitopes from the constructed DIII homology models. The surface charge distribution of the epitopes for the four DENV serotypes are shown in (B) The surface electrostatic potential is represented as a color gradient from blue (positive) to red (negative) between ±1 kT and was calculated using APBS plugin in PyMOL. Respective panels in C display the epitope residues as sticks. Residues labeled in italics are conserved in all four serotypes, while others are mutated across the four serotypes. The region boxed by a rectangle highlights the unique loop orientation adopted by residues 381-386 in DENV3.
shown in Fig. 2C and compared with the equivalent region in the 4E11-DIII-DENV4 complex (Fig. 2D). Overall, the scFv513-DIII-DENV4 complex ( Fig. 2A) is closely superimposable with the scFv4E11-DIII-DENV4 complex 16 and a detailed comparison between their atomic contacts and MD studies described below can account for the higher cross-reactivity displayed by 513 (Table 1). In the following sections, we examine the structural dynamics of these complexes and focus on the regions where 513 was modified and establishes contacts with the DIII antigen that differ from those established by the parent 4E11.
Overall structure of the complexes derived from MD simulations. To reveal the structural details underlying neutralization efficacy and binding variations between 513 and 4E11, MD simulations on atomic models of each scFv fragment bound to DIII from the various serotypes were used. The three complexes between 513 and DIII proteins from DENV1-3 were modeled using the crystal structure of scFv513 bound to DIII from DENV4 as template (PDB code 5AAM). Cartesian coordinates from the experimental crystal structures of 4E11 bound to DIII from the four DENV serotypes were obtained from the PDB (see Table S1) 15 . Each of these eight 3D models was optimized and refined by energy minimization and subjected to MD simulations (see Methods). Overall, the resulting minimized complexes between DIII from each serotype with either 513 and 4E11 can be superimposed with root mean square deviations (RMSD) ranging from 0.7 to 1.0 Å. For DENV serotypes 1 and 4, a larger accessible surface area of DIII becomes buried (depicted as the buried surface area or BSA) upon its interaction with 513 compared to 4E11 (Fig. 3A). In contrast, BSA values are comparable for serotypes 2 and 3 for both mAbs. For both 513 and 4E11, a smaller accessible surface area becomes buried upon interaction with DIII of DENV3 in comparison to other serotypes. This is a consequence of the unique orientation adopted by residues 381-386 in DIII of DENV3 that differs from the other serotypes (Fig. 1C). As a result, this bulging loop restricts the engagement of the antibody in the complex with DIII from DENV3. Per-residue BSA calculations for 513 vs 4E11 (bound to each DENV serotype) indicate that CDR residues are more buried in the case of 513 than 4E11 (Fig. S2). For CDRs of the heavy chain, major differences are observed at residues 25-37 (H1), 50-55 (H2) and 99-108 (H3) and at residues 30-37 (L1) and 52-62 (L2) of the light chain (Fig. S2). Overall, the structures of the DIII from all four DENV serotypes in complex with their respective antibodies are stable during the MD simulations ( Fig. 3B) compared to their free states (Fig. S3). Next, we examined the conformations derived from the MD simulations to provide a mechanism for the observed enhanced binding of 513 to DENV4 compared to 4E11.    (Table 2). Overall the binding of 513 is favored over 4E11 in all cases, especially for DENV4 (Table 2). In the case of DENV1-3, the favorable binding energy derives from gains in electrostatics that are offset by desolvation penalties. Side-chains projecting from β-strands A (residues 305-312) and G (387-391) predominantly define the epitope of all four DENV serotypes. The core of the antigen-antibody interface is composed of hydrophobic residues at positions 308, 312, 387, and 389 as well as residues that undergo only conserved substitutions that maintain their aliphatic character and also the overall shape of the epitope. The formation of conserved contacts between both antibodies and the DIII proteins of all serotypes largely explains the cross-reactivity of 513 and 4E11. In this respect, CDR residues Glu H102 , Tyr L32 , Gln L59 and Trp L60 that interact with main-chain atoms from DIII residues located at the center of the interface allow 513 to bind to all DENV serotypes (Fig. 4A,B). While residues Glu H102 and Tyr L32 are conserved, Gln L59 and Trp L60 are two of the six affinity-enhancing mutations introduced in 513. Moreover, these residues significantly contribute to the total interaction free energies as seen in the next section. Based on the decomposition of the energy terms, it is clear that the increased affinity of 513 for DIII from DENV4 originates from an improved electrostatic stabilization ( Table 2). This is in agreement with the extra salt bridges observed in the crystallographic structure (Table 1). Several substitutions accounting for the observed variations between 513 and 4E11 in binding to DIII from serotype 4 are found at positions 307, 312, 329, 361, 364, 388 and 390 of  by the introduction of six affinity-enhancing mutations: Thr H33 Val (CDR-H1), Ala H55 Glu (CDR-H2) of VH, Arg L31 Lys (CDR-L1), Asn L57 Glu, Glu L59 Gln and Ser L60 Trp (CDR-L2) of VL and Ser26Δ in VH that provided higher surface complementarity (Fig. 2). Remarkably, these substitutions guided by computational modelling enhanced 513 affinity for DIII-DENV4 in the experiments, without affecting the affinity to other serotypes 15 . The impact of these substitutions can be seen in Fig. 5, which shows the free energy contributions of antibodies residues with each of the four DENV serotypes. From the antigen side, a similar energetic analysis is shown in Fig. 6. In addition, these six affinity-enhancing mutations were also subjected to "computational alanine scanning" to evaluate individual side-chain contributions to the overall binding energies. When these 'hot spot' residues were mutated to Alanine in either 513 or 4E11, we observed an overall destabilization of the respective antigen-antibody interface, as measured by the ΔΔG binding energy (Fig. 7). Across all serotypes, residues Glu H55 , Lys L31 and Glu L57 contribute favorably to 513 binding. Using this "computational mutation scanning analysis", Glu H55 from CDR-H2 in 513 appears to be a major contributor towards the enhancement of the binding energy, as analyzed below in greater detail. The Ala55Glu mutation in CDR-H2 dramatically improves the affinity of 513 for DENV4. In the crystal structure, the carboxylic side-chain of Glu H55 forms a salt-bridge with the side-chain of Lys323 (Fig. 2C), a basic residue conserved in serotypes 3 and 4 and substituted by Arg in DENV2 and Gln in DENV1. These three side-chains from the antigen provide comparable favorable electrostatic interactions with Glu H55 , with Gln unsurprisingly being the least favorable (Fig. 6). These interactions were recapitulated in the MD simulations of the respective DENV serotypes in complex with either 4E11 or 513. Upon binding, no major structural changes were observed in the chi1, chi2 and chi3 distribution plot of Glu H55 (Fig. S4). Contributions of individual residues from the CDRs to the total interaction energy show that Glu H55 from 513 contributes significantly (up to −6 kcal/mol) except in the case of DENV1 and DENV3 (Fig. 5). Likewise, residues at position 323 in DENV contribute significantly to the interactions with 513 (Fig. 6). Structural examination of DIII from the four serotypes in complex with 513 reveals that stable salt bridges and hydrogen bonds mediate the interactions between side chains of positively charged (Lys/Arg) residues at position 323 in DENV 2, 3 and 4 with Glu H55 and Asp H32 (with over 90% occupancy). In the case of DENV1, the carbonyl atom of Asp H32 forms a stable H-bond with the Gln323 side chain, with 85% occupancy in both 513 and 4E11. The introduction of a negative charge at Glu H55 in CDR-H2 in the 513 paratope is also likely to favor binding, because of the proximity of Lys310 in DENV1-4 (Figs 1A and 2C). Lys310 has been classified as a conserved epitope residue for 4E11 epitopes 15,18,[25][26][27][28] . It belongs to a bulging loop with conserved shape and charge complementarity across all serotypes (Fig. 1B,C). In the crystal structures, the key interactions between the side chains of Asp H31 (conserved in both 4E11 and 513) and Lys310 of DIII in all serotypes is strictly conserved in the complexes for each serotype 15,18,[25][26][27][28] . The energetics of interactions computed from the MD simulations show that Lys310 contributes an additional 1-4 kcal/mol towards complex formation for 513 compared to 4E11 (Fig. 5A). Interestingly, MD simulations suggest that this is a consequence of the formation of a salt bridge between the side chains of Glu H55 and Lys310, with over 75% occupancy, across all four serotypes. This interaction is obviously missing in 4E11 because of the absence of charge as well as shape complementarity due to the presence of Ala H55 at this position. However, in the aforementioned complexes, Lys310 forms a stable H-bond interaction (over 80% occupancy) with Asp H53 in all DENV serotypes. This interaction is also observed in the case of 513 with similar occupancy.
Experimental validation of residue H55 as a key determinant of the binding affinity. In order to validate experimentally the crucial role of residue H55 for the interaction that was predicted computationally (Fig. 5) and also by a simple visual inspection of the structures (Fig. 2), we targeted this residue from the paratope in both mAb 513 and the parent mAb 4E11. Hence, we introduced the mutation GluH55Ala in scFv513 and conversely the AlaH55Glu mutation in scFv4E11. Using ITC, we then experimentally measured the dissociation constants Kd of these two single mutant antibodies for DIII from DENV2 and compared them with the wild-type (WT) scFv proteins. The reaction parameters (equilibrium dissociation constant K d , enthalphy variation ΔH and stoichiometry N) are summarized in Table 3 for each mutant and WT. As shown in Fig. S5, (Fig. S5). Taken together, this experimentally validates position H55 at the antibody paratope as a hot spot for the binding energetics and shows that the other five modifications introduce in 4E11 to obtain 513 also contribute favorably to the binding energy.
Substitutions in CDR-L2. Interestingly, despite the well-documented dominance of the VH region for determining antigen recognition and affinity, we found that substitutions in the CDR-L1 and CDR-L2 of 4E11 played an important role in affinity improvement displayed by 513. The Asn57Glu mutation (CDR-L2) provides a higher charge complementarity for 513 compared to 4E11. In the crystal structure, Glu L57 (Asn in 4E11) forms a salt bridge with Lys305 in DENV4, a polymorphic residue across the various DENV serotypes. While serotypes 2 and 4 have Lys at this position, serotypes 1 and 3 have small polar side chains: Ser and Thr respectively. The effect of this mutation can be clearly observed as a significant change in the energy contribution (from positive to negative) of Lys305 in serotypes 2 and 4 ( Fig. 6). This originates from a strong and long-lived (over 90% occupancy) H-bond between Glu L57 and the conserved Lys305 in serotypes 2 and 4 (Table S2). This salt bridge contributes to the improved binding of 513 with serotypes 2 and 4. In contrast, Ser305 in serotype 1 is unable to engage in any H-bond interaction and hence makes no contribution (Fig. 6). On the other hand, the side-chain oxygen atom of Thr305 (DENV3) H-bonds with Arg L54 (NH1), however with a moderate strength (distance 3 Å) and lifetime (70% occupancy) in both 4E11 and 513 complexes (Table S2). Another substitution in CDR-L2 that significantly contributes towards binding energy is Glu59Gln (4E11 to 513). The effect of this mutation was examined in the MD simulations, where a stable H-bond is observed between the side chains of Gln L59 and Lys307 in serotypes 1 and 2 with 513 (Table S2). The Lys307 side-chain (DENV1 and 2) also forms an additional H-bond with the backbone atoms of Glu H102 , when in complex with 4E11 and 513. Position 307 is another polymorphic site with a positively charged Lys in serotypes 1 and 2 and Val and Ser in serotypes 3 and 4 respectively. Due to a loss of a positive charge in serotypes 3 and 4 at this position, no favorable/stabilizing interaction is observed in the complex with 4E11. However, the side chain of Ser307 in serotype 4 is stabilized by a strong H-bond (over 95% occupancy) with the Glu H102 side chain (Table S2). Another key observation noted in the 513-DIII-DENV4 complex is that Gln L59 interacts with Glu327 with more than 90% H-bond occupancy (Table S2, Fig. S6). No such interaction is observed in the complex with 4E11 due to the presence of Glu L59 at this position (Fig. S6). Gln L59 in the 513-DIII-DENV4 complex contributes ~ −1.5 kcal/mol towards the total binding energy, as opposed to Glu L59 in the 4E11-DIII-DENV4-complex, which has an unfavorable total energy contribution of ~ 1 kcal/mol (Fig. 5). The per-residue energy decomposition suggests that Glu L59 in DENV4-4E11 causes electrostatic repulsions between negative charge (E ele = 14 kcal/mol for Glu), which is removed by mutation to an uncharged Gln L59 residue in 513 complex and is therefore electrostatically favorable (E ele = −5.5 kcal/mol for Gln). Thus, the Glu59Gln substitution improves the binding of 513 for DENV4.

Identification of hotspot residues in 513 using orthogonal inter-residue interaction approach.
Aside from binding energetics and SASA calculations, we employed an orthogonal inter-residue interactions approach 30,31 to investigate the contributions of ScFv513 residues on DIII binding. This approach enables a quantitative analysis of interactions by computing a numerical network score for each residue based on its interatomic interactions with other residues in its neighborhood 30,31 . The numerical network score for a residue can be correlated to its mutability; the higher the network score, the lower the mutability. Extrapolating from this concept, it is expected that the network scores of paratope residues of an antibody would increase in the antigen-bound state, due to the new pairwise interactions in the epitope-paratope interface. Consequently, residues experiencing a large increase in the network score are likely to be critical hotspots. Next, we computed the differences in the network score of ScFv513 residues in the bound (PDB: 5AAM) versus unbound states. As expected, residue positions undergoing network score changes are restricted to the six CDR loops (Fig. 8). Consistent with the predictions made by the MD simulations, Ala55Glu (CDR-H2) experiences a large increase in the network score upon binding DIII. The network score of affinity enhancing substitutions Thr33Val (CDR-H1), Arg31Lys (CDR-L1) and Ser60Trp (CDR-L2) 19 also increase (Fig. 8). Unexpectedly, Asn57Glu and Glu59Gln (CDR-L2), which led to ~4-5 fold increase in binding affinity towards DENV4, did not experience network score increase. Residue positions other than the affinity enhancing substitutions that showed large increase in network score (e.g. Tyr H34 , Trp H101 , Try L32 ) (Fig. 8) are known hotspot residues, whose modification leads to significant drop in DIII binding 19 .   15 . To characterize these differences in binding affinity of 513 within the DENV4 serotypes observed experimentally, we analyzed the structure and dynamics of the two complexes with either 513 or 4E11. The sequences of DIII from these two DENV4 isolates (New-Caledonia and Philippines) are 97% identical and differ by only three mutations, namely V335I, Y360N, T365I respectively (Fig. 9A). Among these, only residue 365 is near the epitope. The conformations of the two proteins in complex with 513 display an r.m.s.d. of ~1 Å, indicating some differences in antibody positioning: the CDR-L1 is shifted with respect to the end of the A-strand, leading to a loss of charge-charge interactions between Glu311 and Lys L31 . The AG loop and G strand, that form the core of the epitope, are also affected. Polar contact between Thr388 and Asn L34 are also lost. In DENV4 N.Cal all three residues 335, 360 and 365 belong to flexible loops, while in DENV4 Philippines , Thr365 is a part of a β-strand. During the MD simulations, residues 363-365 belong to a disordered region in DENV4 N.Cal , while in DENV4-Philippines a stable β-strand is observed throughout the trajectory (Fig. 9B,C). Comparison of energy contributions of DIII residues from the two strains bound to 513 shows that residues 362-365 of DENV4 N.Cal exhibit favorable binding energies in both electrostatics and van der Waals components (Fig. 10). In the DENV4 Philippines -513 complex, the Thr365 side chain forms an intra-molecular H-bond with the Phe357 backbone (in the adjacent loop), thus restraining the secondary structure in this region (Fig. 10). However, in the DENV4 N.Cal -513 complex, Phe357 is re-orientated due to the longer side chain of Ile335 (Val335 in DENV4 Philippines ). Due to this altered conformation, Phe357 no longer interacts with Ile365 (Fig. 10C) making the corresponding regions more dynamic. However, this structural flexibility allows for spatial adjustment (induced fit) for recognition by the antibody. A comparison of the buried surface area (BSA) of DIII DENV4 Philipines bound to 513 compared to 4E11  indicates a ≅3-fold increase (BSA 4E11:DENV4Philippines = 318 Å and BSA 513:DENV4Philippines = 900 Å). Differences in energy contributions of residues from DENV4 Philippines in complex with 4E11 or 513 suggest that the designed substitutions Thr33Val, Ala55Glu in VH and Glu59Gln, Ser60Trp in VL improved affinity towards DENV4 Philippines (Fig. 10). Moreover, the A-strand of DENV4 Philippines (residues 305-310) has higher shape complementarity in complex with 513 and led to the formation of new H-bonds between Lys305 and Glu L57 and between Phe306 and Arg L54 . Both these interactions are well conserved for over 95% of the equilibrated trajectory. Additional contacts between Glu327 with Gln L59 and Trp L60 are also conserved with over 75% H-bond occupancy (Table S2). In contrast, these favorable interactions are not observed in the case of 4E11 (Fig. 11), which has Asn L57 , Glu L59 and Ser L60 respectively at those positions. Residue His390 also engages in a stabilizing strong H-bond interaction (95% occupancy) with Tyr L32 of 513, an interaction that is absent in the case of 4E11. These improved as well as additional contacts account for the net increase in the affinity of 513 for DENV4 Philippines strain compared to 4E11.

Discussion
MD simulations suggest that the broad recognition by 513 of the four DENV serotypes results from a combination of a larger surface area buried upon complexation compared to 4E11 and specific interactions with conserved functional epitope core residues consisting of Lys310, a polymorphic site 323, Leu387 and Leu/Ile389. This is in agreement with previous experimental observations that were made for the interactions of mAb 4E11 with DENV1-4 32 . While hydrophobic interactions such as those mediated by Leu387 and Leu/Ile 389 remain similar between 4E11 and 513, the conserved charged residue Lys310, a residue deemed important for antibody-antigen recognition across all four serotypes, contributes ~1-2 kcal/mol more to the interactions with 513 (Fig. 6).
Of the six mutations introduced in 4E11 to obtain 513 (CDR-L1: Arg31Lys, Thr33Val; CDR-L2: Asn57Glu, Glu59Gln, Ser60Trp; CDR-H2: Ala55Glu) (Fig. 5), it is clear that Arg (the original residue of 4E11) and Gln are largely preferred at positions 31 (CDR-L1) and 59 (CDR-L2) respectively. Mutations at Glu L57 and Trp L60 also lead to substantial increases in affinity for 513 33,34 . Glu H55 has a very large contribution in improving the interactions of 513 with the four serotypes. Indeed, computational Ala mutation analysis (Fig. 7) shows that residue at position H55 is the most important. It also has a large network score (Fig. 8) calculated from the orthogonal inter-residue interaction approach. In addition to the contacts observed in the crystal structure, MD simulations suggest that Glu H55 in 513 forms a stable interaction with the strictly conserved charged Lys310 in the DENV serotypes. The charge complementarity is further maintained in this region as the adjacent conserved acidic residue, Glu311, forms a stable interaction with Lys L31 in the 513 (CDR-L1) paratope. In addition, Glu H55 recognizes and stabilizes residue 323, which is a polymorphic residue on the DENV epitope, with significantly higher binding energy. In contrast to the observations made in the crystal structure (wherein Lys H31 is observed to form a H-bond with Lys323), MD simulations for all four DENV serotypes bound to 513 indicate that Asp H31 and Glu H55 stabilize this polymorphic site in concert. The Ala55Glu mutation increases charge complementarity across the paratope-epitope in this region, while the extended side chain of Glu endows broader specificity with increased contributing surface area. Importantly, the energetic contribution of Glu H55 was confirmed experimentally by site-directed mutagenesis followed by affinity measurement using ITC: Glu55Ala mutation in 513 showed similar affinity as the 4E11 WT, whereas Ala55Glu mutation in 4E11 showed enhanced affinity towards DENV2 DIII ( Fig. S5 and Table 3).
The second major contributor towards the binding efficacy of 513 identified from our work is the Asn57Glu (CDR-L2) mutation, which specifically elicits favorable response due to charge complementarity for DENV2 and 4. The other two most important contributing residues are Tyr H34 and Tyr L32 (Figs 5 and 8). These two paratope residues, conjointly with Asp H32 and Lys L31 , engage in stable H-bond interactions with the A-strand (311, 312), B-strand (323) and G strand (390) residues (Table S2) forming strong interactions at the periphery of the interface. Moreover, Gln and Trp residues at positions 59 and 60 (CDR-L2) form stable H-bonds with Glu327 specific for DENV4-513 complex that further enhances binding energetics. This bond is not observed in other DENV serotypes due to conformational variability.
In general, 513 is endowed with cross-neutralizing and enhancing activity for DENV1-4 serotypes, however significant differences in the binding efficacies were observed within DENV4 strains (see Table 2 of ref. 14 ) such as DENV4 N.Cal and DENV4 Phillipines . From our MD simulations, subtle sequence variations affect the local conformation of the functional epitope, altering shape complementarity and hence antibody binding. However, even though the binding constants are significantly different, the adjacent regions around the mutations largely preserve charge complementarity and the paratope-epitope interface.
In conclusion, a combination of physico-chemical and structural changes brought about by engineering six mutations in 4E11 result in the generation of mAb 513, which potently inhibits all four DENV serotypes. Combining MD simulations and residue-network analyses can therefore be a powerful complement to crystallographic and binding data and illuminate detailed interactions that underpin the differential binding of antibodies to homologous viral antigens. This has implications for the design of therapeutic cross-reactive antibodies against emerging viruses and also for the design of vaccines by "reverse vaccinology".

Materials and Methods
Homology Modeling. The 35 . Table S2 shows the list of PDB structures used to model the virus isolates as well as the mAbs and the complex structures. Point mutations in various virus serotypes were incorporated in the models and missing regions, if any, were modeled. The scFv513 and 4E11 antibody structures were extracted from the crystal structures 4UDZ and 3UZQ, followed by modeling the complexes with the respective serotypes. The files for electrostatic calculations were generated using PDB2PQR 36 server and APBS 37 plug-in in PyMol 38 was used for electrostatic surface illustrations. Figure 11. Analysis of energetics of DENV4 New Caledonia and DENV4 Philippine isolates (right panel). Each plot depicts the free energy contributions of residues from the DENV4 New Caledonia (left panels) and DENV4 Philippine (right panels), when bound mAb 4E11 and mAb 513. For each isolate, a negative value depicts favorable contribution to energetics when bound to mAb 513, whereas a positive value depicts favorable contribution to energetics when bound to mAb 4E11.
Scientific RepoRtS | (2018) 8:8449 | DOI:10.1038/s41598-018-26800-y Molecular Dynamics. The modeled complexes were subjected to molecular dynamics simulations for further refinement. The ff99SB force field as implemented in the program AMBER 39 was used. In all the complexes, hydrogen atoms were added using the Xleap module of Amber12 39 to prepare the systems for MD simulations. Each system was neutralized by adding sufficient number of Na + and Clcounter ions. All the 8 complexes were then solvated in an octahedral box with TIP3P 40 water molecules that extended 10 Å from any protein atom. The short range non-bonded van der Waals interactions were truncated at 9 Å while the long range electrostatics were approximated by the particle mesh Ewald 41 method. The covalent bonds involving hydrogen atoms were constrained using SHAKE 42 . The Sander module was used for minimization of the complexes with 250 steps of steepest decent algorithm, followed by 8000 steps of conjugate gradient algorithm. Initially, the antibody atoms, solvent water molecules and counter ions were relaxed, by keeping the virus residues restrained. This was followed by unrestrained energy minimization to remove any steric clashes. The systems were subject to over 250 ps of heating from 50 to 300 K with weak restraints on the heavy atoms, followed by gradual reduction of the restraints over the next 250 ps, until the restraints were reduced to 0. For next 2 ns the system was equilibrated at 300 K under 1 atm constant pressure. After equilibration production runs of 50 ns with 2 fs time step were carried out on each complex and the atomic coordinates were saved every 10 ps. Analyses of the simulation trajectories were performed using the ptraj module in Amber. The conformations generated were oriented to a common frame to remove overall translation/rotations. The temporal evolution of the RMSD of the Cα atoms of various sampled conformations suggested that the systems had equilibrated within 10 ns. The free energy binding was calculated for each complex using MMPBSA 43,44 methodology. For the binding energy calculations, a total of 50 structures were extracted at regular intervals from last 30 ns of the equilibrated trajectory. Water molecules and ions were stripped from all extracted snapshots and born implicit solvent model of 5 (igb = 5) was used to compute the free energies of binding. The effective binding energies were decomposed into contributions of individual residues using MMGBSA energy decomposition scheme. To investigate the impact of mutations in the antibodies on serotype recognition, selected hot-spot residues from the antibodies in each complex were computationally mutated to alanine. The free energy changes (ΔΔG bind ) were defined as ΔG bind (original residue) -ΔG bind (alanine mutant) for the alanine mutation. A negative ΔΔG bind value for that specific residue from the antibody indicated that the existing residue has a more favorable binding towards the respective DENV serotype. The solvent accessible surface areas were calculated for every 50 th frame in the equilibrated trajectory. Using the average values obtained above, BSA values were then calculated as the sum of the solvent accessible surface area of individual molecules minus the solvent accessible surface area of their complex. Calculations were performed using Naccess V2.1.1 program 45,46 with a probe radius of 1.4 Å. The simulated trajectory was viewed in VMD 47 and figures were generated using PyMol 38 .