Stability of trimeric DENV envelope protein at low and neutral pH: An insight from MD study

Homotypic antibodies target novel E glycoprotein domains after natural DENV 3 infection/vaccination

The envelope (E) glycoprotein is the primary target of type-specific (TS) neutralizing antibodies (nAbs) after infection with any of the four distinct dengue virus serotypes (DENV1-4). nAbs can be elicited to distinct structural E domains (EDs) I, II, or III. However, the relative contribution of these domain-specific antibodies is unclear. To identify the primary DENV3 nAb targets in sera after natural infection or vaccination, chimeric DENV1 recombinant encoding DENV3 EDI, EDII, or EDIII were generated. DENV3 EDII is the principal target of TS polyclonal nAb responses and encodes two or more neutralizing epitopes. In contrast, some were individuals vaccinated with a DENV3 monovalent vaccine-elicited serum TS nAbs targeting each ED in a subject-dependent fashion, with an emphasis on EDI and EDIII. Vaccine responses were also sensitive to DENV3 genotypic variation. This DENV1/3 panel allows the measurement of serum ED TS nAbs, revealing differences in TS nAb immunity after natural infection or vaccination.

Unraveling key interactions and the mechanism of demethylation during hAGT-mediated DNA repair via simulations

Alkylating agents pose the biggest threat to the genomic integrity of cells by damaging DNA bases through regular alkylation. Such damages are repaired by several automated types of machinery inside the cell. O6-alkylguanine-DNA alkyltransferase (AGT) is an enzyme that performs the direct repair of an alkylated guanine base by transferring the alkyl group to a cysteine residue. In the present study, using extensive MD simulations and hybrid QM/MM calculations, we have investigated the key interactions between the DNA lesion and the hAGT enzyme and elucidated the mechanisms of the demethylation of the guanine base. Our simulation shows that the DNA lesion is electrostatically stabilized by the enzyme and the Arg135 of hAGT enzyme provides the main driving force to flip the damaged base into the enzyme. The QM/MM calculations show demethylation of the damaged base as a three-step process in a thermodynamically feasible and irreversible manner. Our calculations show that the final product forms via Tyr114 in a facile way in contrast to the previously proposed Lys-mediated route.

Can the local electric field be a descriptor of catalytic activity? A case study on chorismate mutase

The current theoretical perception of enzymatic activity is highly reliant on the determination of the activation energy of the reactions, which is often calculated using computationally demanding quantum mechanical calculations. With the ever-increasing use of bioengineering techniques that produce too many variants of the same enzyme, a fast and accurate way to study the relative efficiency of enzymes is currently in high demand. Here, we propose the local electric field (LEF) of the enzyme along the reaction axis as a descriptor for the enzymatic activity using the example of chorismate mutase in its native form and several variants (R90A, R90G, and R90K/C88S). The study shows a direct correlation between the calculated enzymatic EF and the enzymatic activity for all the complexes. MD simulations of the Michaelis complex and the transition state analog (TSA) show a stabilizing force on the TSA due to the enzymatic EF. QM/MM and QM-only DFT calculations in the presence of an external electric field (EEF) oriented along the reaction axis show that the electric field can interact with the dipole moment of the TS, thereby stabilizing it and thus lowering the activation energy.

Structural and biochemical insights into flavivirus proteins

Flaviviruses are the fastest spreading arthropod-borne viruses that cause severe symptoms such as hepatitis, hemorrhagic fever, encephalitis, and congenital deformities. Nearly 40 % of the entire human population is at risk of flavivirus epidemics. Yet, effective vaccination is restricted only to a few flaviviruses such as yellow fever and Japanese encephalitis viruses, and most recently for select cases of dengue virus infections. Despite the global spread of dengue virus, and emergence of new threats such as Zika virus and a new genotype of Japanese encephalitis virus, insights into flavivirus targets for potentially broad-spectrum vaccination are limited. In this review article, we highlight biochemical and structural differences in flavivirus proteins critical for virus assembly and host interactions. A comparative sequence analysis of pH-responsive properties of viral structural proteins identifies trends in conservation of complementary acidic-basic character between interacting viral structural proteins. This is highly relevant to the understanding of pH-sensitive differences in virus assembly in organelles such as neutral ER and acidic Golgi. Surface residues in viral interfaces identified by structural approaches are shown to demonstrate partial conservation, further reinforcing virus-specificity in assembly and interactions with host proteins. A comparative analysis of epitope conservation in emerging flaviviruses identifies therapeutic antibody candidates that have potential as broad spectrum anti-virals, thus providing a path towards development of vaccines.

Does Antibody Stabilize the Ligand Binding in GP120 of HIV-1 Envelope Protein? Evidence from MD Simulation

CD4-mimetic HIV-1 entry inhibitors are small sized molecules which imitate similar conformational flexibility, in gp120, to the CD4 receptor. However, the mechanism of the conformational flexibility instigated by these small sized inhibitors is little known. Likewise, the effect of the antibody on the function of these inhibitors is also less studied. In this study, we present a thorough inspection of the mechanism of the conformational flexibility induced by a CD4-mimetic inhibitor, NBD-557, using Molecular Dynamics Simulations and free energy calculations. Our result shows the functional importance of Asn425 in substrate induced conformational dynamics in gp120. The MD simulations of Asn425Gly mutant provide a less dynamic gp120 in the presence of NBD-557 without incapacitating the binding enthalpy of NBD-557. The MD simulations of complexes with the antibody clearly show the enhanced affinity of NBD-557 due to the presence of the antibody, which is in good agreement with experimental Isothermal Titration Calorimetry results (Biochemistry2006, 45, 10973-10980).

Computational modelling of flavivirus dynamics: The ins and outs

Enveloped viruses such as the flaviviruses represent a significant burden to human health around the world, with hundreds of millions of people each year affected by dengue alone. In an effort to improve our understanding of the molecular basis for the infective mechanisms of these viruses, extensive computational modelling approaches have been applied to elucidate their conformational dynamics. Multiscale protocols have been developed to simulate flavivirus envelopes in close accordance with biophysical data, in particular derived from cryo-electron microscopy, enabling high-resolution refinement of their structures and elucidation of the conformational changes associated with adaptation both to host environments and to immunological factors such as antibodies. Likewise, integrative modelling efforts combining data from biophysical experiments and from genome sequencing with chemical modification are providing unparalleled insights into the architecture of the previously unresolved nucleocapsid complex. Collectively, this work provides the basis for the future rational design of new antiviral therapeutics and vaccine development strategies targeting enveloped viruses.

Investigating the stability of dengue virus envelope protein dimer using well-tempered metadynamics simulations

We explored the stability of the dengue virus envelope (E) protein dimer since it is widely assumed that the E protein dimer is stabilized by drug ligands or antibodies in an acidic environment, neutralizing the virus's ability to fuse with human cells. During this process, a large conformational change of the E protein dimer is required. We performed Molecular Dynamics simulations to mimic the conformational change and stability of the dimer in neutral and acidic conditions with the well-tempered metadynamics method. Furthermore, as a few neutralizing antibodies discovered from dengue patients were reported, we used the same simulation method to examine the influence of a selected antibody on the dimer stability in both neutral and acidic conditions. We also investigated the antibody's influence on a point-mutated E protein that had been reported to interrupt the protein-antibody interaction and result in more than 95% loss of the antibody's binding ability. Our simulation results are highly consistent with the experimental conclusion that binding of the antibody to the E protein dimer neutralizes the virus, especially in a low pH condition, while the mutation of W101A or N153A significantly reduces the antibody's ability in stabilizing the E protein dimer. We demonstrate that well-tempered metadynamics can be used to accurately explore the antibody's interaction on large protein complexes such as the E protein dimer, and the computational approach in this work is promising in future antibody development.

Single-molecule studies of flavivirus envelope dynamics: Experiment and computation

Flaviviruses are simple enveloped viruses exhibiting complex structural and functional heterogeneities. Decades of research have provided crucial basic insights, antiviral medication and moderately successful gene therapy trials. The most infectious particle is, however, not always the most abundant one in a population, questioning the utility of classic ensemble-averaging virology approaches. Indeed, viral replication is often not particularly efficient, prone to errors or containing parallel routes. Here, we review different single-molecule sensitive fluorescence methods that are employed to investigate flaviviruses. In particular, we review how (i) time-resolved Förster resonance energy transfer (trFRET) was applied to probe dengue envelope conformations; (ii) FRET-fluorescence correlation spectroscopy to investigate dengue envelope intrinsic dynamics and (iii) single particle tracking to follow the path of dengue viruses in cells. We also discuss how such methods may be supported by molecular dynamics (MD) simulations over a range of spatio-temporal scales, to provide complementary data on the structure and dynamics of flaviviral systems. We describe recent improvements in multiscale MD approaches that allowed the simulation of dengue particle envelopes in near-atomic resolution. We hope this review is an incentive for setting up and applying similar single-molecule studies and combine them with MD simulations to investigate structural dynamics of entire flavivirus particles over the nanosecond-to-millisecond time-scale and follow viruses during infection in cells over milliseconds to minutes.

Targeting domain-III hinging of dengue envelope (DENV-2) protein by MD simulations, docking and free energy calculations

The entry of the dengue virus is mediated by the conformational change in the envelope protein due to change in the endosomal pH. The structural study reveals that domain-III of the dengue envelope protein (DENV) shows the largest shift in its position during the entry of the virus. Therefore, targeting the hinge region of the domain-III may block the conformational changes in the DENV. In the present work, we have targeted the domain I/III hinge region using four known ligands used for the dengue envelope protein (serotype-2) and have intended to explore the specificity of one ligand R1 (5-(3-chlorophenyl)-N-(2-phenyl-2H-benzo[d][1,2,3]triazol-6-yl)furan-2-carboxamide) that succeeded the dengue inhibition by the molecular dynamics (MD) simulations in conjunction of the molecular docking and the binding free energy calculations. The residue interactions map shows Lys 296 of domain-III of the DENV-2, which might be responsible for binding small molecules between domain I/III interface, as an important residue conserved in all the dengue serotypes.

Structure-based pKa prediction provides a thermodynamic basis for the role of histidines in pH-induced conformational transitions in dengue virus

pH-induced conformational changes in dengue virus (DENV) are critical to its ability to infect host cells. The envelope protein heterodimers that make up the viral envelope shift from a dimer to a trimer conformation at low-pH during membrane fusion. Previous studies have suggested that the ionization of histidine residues at low-pH is central to this pH-induced conformational change. We sought out to use molecular modeling with structure-based pKa prediction to provide a quantitative basis for the role of histidines in pH-induced conformational changes and identify which histidine residues were primarily responsible for this transition. We combined existing crystallographic and cryo-electron microscopy data to construct templates of the dimer and trimer conformations for the mature and immature virus. We then generated homology models for the four DENV serotypes and carried out structure-based pKa prediction using Rosetta. Our results showed that the pKa values of a subset of conserved histidines in DENV successfully capture the thermodynamics necessary to drive pH-induced conformational changes during fusion. Here, we identified the structural determinants underlying these pKa values and compare our findings with previous experimental results.

•

Structure-based pKa prediction for histidines in dengue virus was carried out.

•

Conserved histidine residues drive pH-dependent conformation changes.

•

The mature dimer form of the envelop protein is destabilized at low-pH.

Structure of a dengue virus envelope protein late-stage fusion intermediate

The final stages of dengue virus fusion are thought to occur when the membrane-proximal stem drives the transmembrane anchor of the viral envelope protein (E) toward the fusion loop, buried in the target cell membrane. Crystal structures of E have lacked this essential stem region. We expressed and crystallized soluble mutant forms of the dengue virus envelope protein (sE) that include portions of the juxtamembrane stem. Their structures represent late-stage fusion intermediates. The proximal part of the stem has both intra- and intermolecular interactions, so the chain "zips up" along the trimer seam. The penultimate interaction we detected involves the conserved residue F402, which has hydrophobic contacts with a conserved surface on domain II. These interactions do not require any larger-scale changes in trimer packing. The techniques for expression and crystallization of sE containing stem reported here may allow further characterization of the final stages of flavivirus fusion.