How the Strain Origin of Zika Virus NS1 Protein Impacts Its Dynamics and Implications to Their Differential Virulence

Differences in the Membrane-Binding Properties of Flaviviral Nonstructural 1 (NS1) Protein: Comparative Simulations of Zika and Dengue Virus NS1 Proteins in Explicit Bilayers

NS1 in flaviviruses is the only nonstructural protein that is secretory and interacts with different cellular components of the host cell membrane. NS1 is localized in the ER as a dimer to facilitate viral replication. Crystal structures of NS1 homologues from zika (ZIKV) and dengue (DENV) viruses have revealed the organization of different domains in NS1 dimers. The β-roll and the connector and intertwined loop regions of wing domains of NS1 have been shown to interact with the membranes. In this study, we have performed multiple molecular dynamics (MD) simulations of ZIKV and DENV NS1 systems in apo and in POPE bilayers with different cholesterol concentrations (0, 20 and 40%). The NS1 protein was placed just above the membrane surface, and for each NS1-membrane system two to three independent simulations with 600 ns production run were performed. At the end of the production runs, ZIKV NS1 inserts deeper inside the membrane compared to the DENV counterpart. Unlike ZIKV NS1, the orientation of DENV NS1 is asymmetric in which one of the chains in the dimer interacts with the membrane while the other is more exposed to the solvent. The β-roll region in ZIKV NS1 penetrates beyond the headgroup region and interacts with the lipid acyl chains while the C-terminal region barely interacts with the headgroup. Specific residues in the intertwined region deeply penetrate inside the membrane. The role of charged and aromatic residues of ZIKV NS1 in strongly interacting with the membrane components is revealed. The presence of cholesterol affects the extent of insertion in the membrane and interaction of individual residues. Overall, membrane-binding properties of ZIKV NS1 significantly differ from its counterpart in DENV. The differences found in the binding and insertion of NS1 can be used to design drugs and novel antibodies that can be flavivirus specific.

Investigating the aggregation perspective of Dengue virus proteome

Dengue viruses are human pathogens that are transmitted through mosquitoes. Apart from the typical symptoms associated with viral fevers, DENV infections are known to cause several neurological complications such as meningitis, encephalitis, intracranial haemorrhage, retinopathies along with the more severe, and sometimes fatal, vascular leakage and dengue shock syndrome. This study was designed to investigate, in detail, the predicted viral protein aggregation prone regions among all serotypes. Further, in order to understand the cross-talk between viral protein aggregation and aggregation of cellular proteins, cross-seeding experiments between the DENV NS1 (1-30), corresponding to the β-roll domain and the diabetes hallmark protein, amylin, were performed. Various techniques such as fluorescence spectroscopy, circular dichroism, atomic force microscopy and immunoblotting have been employed for this. We observe that the DENV proteomes have many predicted APRs and the NS1 (1-30) of DENV1-3, 2K and capsid anchor of DENV2 and DENV4 are capable of forming amyloids, in vitro. Further, the DENV NS1 (1-30), aggregates are also able to cross-seed and enhance amylin aggregation and vice-versa. This knowledge may lead to an opportunity for designing suitable inhibitors of protein aggregation that may be beneficial for viral infections and comorbidities.

Antiviral Activity of Quercetin Hydrate against Zika Virus

Zika virus (ZIKV) has re-emerged in recent decades, leading to outbreaks of Zika fever in Africa, Asia, and Central and South America. Despite its drastic re-emergence and clinical impact, no vaccines or antiviral compounds are available to prevent or control ZIKV infection. This study evaluated the potential antiviral activity of quercetin hydrate against ZIKV infection and demonstrated that this substance inhibits virus particle production in A549 and Vero cells under different treatment conditions. In vitro antiviral activity was long-lasting (still observed 72 h post-infection), suggesting that quercetin hydrate affects multiple rounds of ZIKV replication. Molecular docking indicates that quercetin hydrate can efficiently interact with the specific allosteric binding site cavity of the NS2B-NS3 proteases and NS1-dimer. These results identify quercetin as a potential compound to combat ZIKV infection in vitro.

NS1 from Two Zika Virus Strains Differently Interact with a Membrane: Insights to Understand Their Differential Virulence

Zika virus (ZIKV) from Uganda (UG) expresses a phenotype related to fetal loss, whereas the variant from Brazil (BR) induces microcephaly in neonates. The differential virulence has a direct relation to biomolecular mechanisms that make one strain more aggressive than the other. The nonstructural protein 1 (NS1) is a key viral toxin to comprehend these viral discrepancies because of its versatility in many processes of the virus life cycle. Here, we aim to examine through coarse-grained models and molecular dynamics simulations the protein-membrane interactions for both NS1_ZIKV-UG and NS1_ZIKV-BR dimers. A first evaluation allowed us to establish that the NS1 proteins, in the membrane presence, explore new conformational spaces when compared to systems simulated without a lipid bilayer. These events derive from both differential coupling patterns and discrepant binding affinities to the membrane. The N-terminal domain, intertwined loop, and greasy finger proposed previously as binding membrane regions were also computationally confirmed by us. The anchoring sites have aromatic and ionizable residues that manage the assembly of NS1 toward the membrane, especially for the Ugandan variant. Furthermore, in the presence of the membrane, the difference in the dynamic cross-correlation of residues between the two strains is particularly pronounced in the putative epitope regions. This suggests that the protein-membrane interaction induces changes in the distal part related to putative epitopes. Taken together, these results open up new strategies for the treatment of flaviviruses that would specifically target these dynamic differences.

Quantitative electrostatic force tomography for virus capsids in interaction with an approaching nanoscale probe

Electrostatic interactions are crucial for the assembly, disassembly and stability of proteinaceous viral capsids. Moreover, at the molecular scale, elucidating the organization and structure of the capsid proteins in response to an approaching nanoprobe is a major challenge in biomacromolecular research. Here, we report on a generalized electrostatic model, based on the Poisson-Boltzmann equation, that quantifies the subnanometric electrostatic interactions between an AFM tip and a proteinaceous capsid from molecular snapshots. This allows us to describe the contributions of specific amino acids and atoms to the interaction force. We show validation results in terms of total electrostatic forces with previous semi-empirical generalized models at available length scales (d > 1 nm). Then, we studied the interaction of the Zika capsid with conical and spherical AFM tips in a tomography-type analysis to identify the most important residues and atoms, showing the localized nature of the interaction. This method can be employed for the interpretation of force microscopy experiments in fundamental virological characterization and in diverse nanomedicine applications, where specific regions of the protein cages are aimed to electrostatically interact with molecular sized functionalized inhibitors, or tailoring protein-cage functional properties for nucleic acid delivery.

Structural dynamics of Zika virus NS1 via a reductionist approach reveal the disordered nature of its β-roll domain in isolation

Flavivirus Non-structural 1 (NS1) protein performs multiple functions and it is highly plausible that significant structural and folding dynamics of NS1 might play a role in its multifunctionality. It is important to understand the structural conformations of NS1 and its domains in isolation, possibly highlighting the implications on the overall NS1 protein dynamics. Therefore, we have employed extensively long molecular dynamic (MD) simulations in understanding the dynamics of the three structural domains (i.e., β-roll, wing, and β-ladder) in isolation, as a reductionist approach. We also found that the β-ladder domain is highly flexible, while the β-roll domain is disordered during long simulations. Further, we experimentally validated our findings using CD spectroscopy and confirmed the intrinsically disordered behavior of NS1 β-roll in isolation and lipid mimetic environments. Therefore, we believe this study may have implications for significant dynamics played by NS1 protein, specifically during oligomerization of NS1.

Self-association features of NS1 proteins from different flaviviruses

Flaviviruses comprise a large group of arboviral species that are distributed in several countries of the tropics, neotropics, and some temperate zones. Since they can produce neurological pathologies or vascular damage, there has been intense research seeking better diagnosis and treatments for their infections in the last decades. The flavivirus NS1 protein is a relevant clinical target because it is involved in viral replication, immune evasion, and virulence. Being a key factor in endothelial and tissue-specific modulation, NS1 has been largely studied to understand the molecular mechanisms exploited by the virus to reprogram host cells. A central part of the viral maturation processes is the NS1 oligomerization because many stages rely on these protein-protein assemblies. In the present study, the self-associations of NS1 proteins from Zika, Dengue, and West Nile viruses are examined through constant-pH coarse-grained biophysical simulations. Free energies of interactions were estimated for different oligomeric states and pH conditions. Our results show that these proteins can form both dimers and tetramers under conditions near physiological pH even without the presence of lipids. Moreover, pH plays an important role mainly controlling the regimes where van der Waals interactions govern their association. Finally, despite the similarity at the sequence level, we found that each flavivirus has a well-characteristic protein-protein interaction profile. These specific features can provide new hints for the development of binders both for better diagnostic tools and the formulation of new therapeutic drugs.

Up State of the SARS-COV-2 Spike Homotrimer Favors an Increased Virulence for New Variants

The COVID-19 pandemic has spread worldwide. However, as soon as the first vaccines-the only scientifically verified and efficient therapeutic option thus far-were released, mutations combined into variants of SARS-CoV-2 that are more transmissible and virulent emerged, raising doubts about their efficiency. This study aims to explain possible molecular mechanisms responsible for the increased transmissibility and the increased rate of hospitalizations related to the new variants. A combination of theoretical methods was employed. Constant-pH Monte Carlo simulations were carried out to quantify the stability of several spike trimeric structures at different conformational states and the free energy of interactions between the receptor-binding domain (RBD) and angiotensin-converting enzyme II (ACE2) for the most worrying variants. Electrostatic epitopes were mapped using the PROCEEDpKa method. These analyses showed that the increased virulence is more likely to be due to the improved stability to the S trimer in the opened state, in which the virus can interact with the cellular receptor, ACE2, rather than due to alterations in the complexation RBD-ACE2, since the difference observed in the free energy values was small (although more attractive in general). Conversely, the South African/Beta variant (B.1.351), compared with the SARS-CoV-2 wild type (wt), is much more stable in the opened state with one or two RBDs in the up position than in the closed state with three RBDs in the down position favoring the infection. Such results contribute to understanding the natural history of disease and indicate possible strategies for developing new therapeutic molecules and adjusting the vaccine doses for higher B-cell antibody production.