Microsecond Molecular Dynamics Simulations of Influenza Neuraminidase Suggest a Mechanism for the Increased Virulence of Stalk-Deletion Mutants

Enhancing NA immunogenicity through novel VLP designs

Current influenza virus vaccines poorly display key neuraminidase (NA) epitopes and do not robustly induce NA-reactive antibodies; instead, they focus on the induction of hemagglutinin (HA)-reactive antibodies. Next-generation influenza vaccines should be optimized in order to activate NA-reactive B cells and to induce a broadly cross-reactive and protective antibody response. We aimed at enhancing the immunogenicity of the NA on vaccines by two strategies: (i) modifying the HA:NA ratio of the vaccine preparation and (ii) exposing epitopes on the lateral surface or beneath the head of the NA by extending the NA stalk. The H1N1 glycoproteins from the influenza virus A/California/04/2009 strain were displayed on human immunodeficiency virus 1 (HIV-1) gag-based virus-like particles (VLP). Using the baculovirus insect cell expression system, we biased the quantity of surface glycoproteins employing two different promoters, the very late baculovirus p10 promoter and the early and late gp64 promoter. This led to a 1:1 to 2:1 HA:NA ratio, which was approximately double or triple the amount of NA as present on the wild-type influenza A virus (HA:NA ratio 3:1 to 5:1). Furthermore, by insertion of 15 amino acids from the A-New York/61/2012 strain (NY12) which prolongates the NA stalk (NA long stalk; NA-LS), we intended to improve the accessibility of the NA. Six different types of VLPs were produced and purified using a platform downstream process based on Capto-Core 700™ followed by Capto-Heparin™ affinity chromatography combined with ultracentrifugation. These VLPs were then tested in a mouse model. Robust titers of antibodies that inhibit the neuraminidase activity were elicited even after vaccination with two low doses (0.3 μg) of the H1N1 VLPs without compromising the anti-HA responses. In conclusion, our results demonstrate the feasibility of the two developed strategies to retain HA immunogenicity and improve NA immunogenicity as a future influenza vaccine candidate.

Structural and functional characterization of avian influenza H9N2 virus neuraminidase with a combination of five novel mutations

The influenza virus neuraminidase (NA) protein is responsible for actively cleaving the sialic acid (SA) bound to the viral hemagglutinin. In the present study, we identified a combination of five novel amino acid substitutions in the NA, conferring increased substrate binding and altered surface characteristics to a low pathogenic avian influenza (LPAI) H9N2 virus strain. The H9N2 strain reported from India, A/Environmental/India/1726265/2017 (H9N2-1726265) showed the combination of amino acid substitutions T149I, R249W, G346A, W403R and G435R, which were in the vicinity of the enzyme active site cavity. The strain A/chicken/India/99321/2009 (H9N2-99321) did not show these substitutions and was used for comparison. Virus elution was studied using turkey red blood cells (tRBCs). NA enzyme kinetics assays were carried out using the MUNANA substrate, which is an SA analogue. Homology modelling and molecular docking were performed to determine alterations in the surface characteristics and substrate binding. H9N2-1726265 showed enhanced elution from tRBCs. Enzyme kinetics revealed a lower KM of H9N2-1726265 (111.5 μM) as compared to H9N2-99321 (135.2 μM), indicating higher substrate binding affinity of H9N2-1726265, due to which the NA enzyme cleaved the SA more efficiently, leading to faster elution. Molecular docking revealed a greater number of binding interactions of H9N2-1726265 to SA as compared to H9N2-99321 corroborating the greater substrate binding affinity. Changes in the surface charge, hydrophobicity, and contour, were observed in H9N2-1726265 NA due to the five substitutions. Thus, the novel combination of five amino acids near the sialic acid binding site of NA, resulted in altered surface characteristics, higher substrate binding affinity, and virus elution.

Loss of amino acids 67-76 in the neuraminidase protein under antibody selection pressure alters the tropism, transmissibility and innate immune response of H9N2 avian influenza virus in chickens

H9N2 virus has become the most widespread subtype of avian influenza in Chinese poultry. Although many studies have been published on this disease, the pathogenesis of the H9N2 virus remains to be fully understood. In our previous work, we identified 44 viral strains with 67-76 amino acid deletions in the neuraminidase protein (NA∆67-76) from trachea and lung tissues after 20 successive generations in vaccinated chickens. Interestingly, these 10 amino acid deletions are located in the stalk of the NA protein, and all mutations were unique to the viruses under the selection pressure of vaccine antibodies. To investigate the effect of NA∆67-76 on the H9N2 virus, the NA∆67-76 deletion mutant (rF/NAΔ67-76) was constructed in the H9N2 virus A/Chicken/Shanghai/F/98 (F/98) to assess the phenotypic changes between the parental and mutant strains. The results showed that the recombinant virus rF/NAΔ67-76 had no significantly effect on the antigenicity of the virus or on the infectivity of the host cells, but it significantly inhibited the release of virions from host cells. In addition, rF/NAΔ67-76 efficiently enhanced the neuraminidase activity and improved the receptor binding ability of the virus, indicating that the influence of receptor binding ability on the rF/NAΔ67-76 virus is much greater than that of neuraminidase activity. Furthermore, this study revealed that rF/NAΔ67-76 reduced the viral replication ability at 6 and 12 h post-infection, but improved it at 24, 48, and 72 h post-infection. Chicken experiments showed that rF/NAΔ67-76 exhibits a much higher tissue tropism for the trachea rather than lung tissue. rF/NAΔ67-76 still had the ability to infect the upper respiratory tract through aerosol, but its cloaca replication capacity was significantly reduced. Both in vivo and in vitro experiments confirmed that rF/NAΔ67-76 could produce a stronger innate immune response after infecting cells and chickens, especially significantly enhancing the transcription levels of TLR3, TLR4, TLR7, TLR21, MDA5, and NLRP3. Altogether, the results of this study propose that antibody selection pressure plays an important role in the evolution of H9N2 avian influenza virus.

Molecular Dynamics Simulations Reveal the Conformational Transition of GH33 Sialidases

Sialidases are increasingly used in the production of sialyloligosaccharides, a significant component of human milk oligosaccharides. Elucidating the catalytic mechanism of sialidases is critical for the rational design of better biocatalysts, thereby facilitating the industrial production of sialyloligosaccharides. Through comparative all-atom molecular dynamics simulations, we investigated the structural dynamics of sialidases in Glycoside Hydrolase family 33 (GH33). Interestingly, several sialidases displayed significant conformational transition and formed a new cleft in the simulations. The new cleft was adjacent to the innate active site of the enzyme, which serves to accommodate the glycosyl acceptor. Furthermore, the residues involved in the specific interactions with the substrate were evolutionarily conserved in the whole GH33 family, highlighting their key roles in the catalysis of GH33 sialidases. Our results enriched the catalytic mechanism of GH33 sialidases, with potential implications in the rational design of sialidases.

Exploration of Computational Aids for Effective Drug Designing and Management of Viral Diseases: A Comprehensive Review

BACKGROUND

Microbial diseases, specifically originating from viruses are the major cause of human mortality all over the world. The current COVID-19 pandemic is a case in point, where the dynamics of the viral-human interactions are still not completely understood, making its treatment a case of trial and error. Scientists are struggling to devise a strategy to contain the pandemic for over a year and this brings to light the lack of understanding of how the virus grows and multiplies in the human body.

METHODS

This paper presents the perspective of the authors on the applicability of computational tools for deep learning and understanding of host-microbe interaction, disease progression and management, drug resistance and immune modulation through in silico methodologies which can aid in effective and selective drug development. The paper has summarized advances in the last five years. The studies published and indexed in leading databases have been included in the review.

RESULTS

Computational systems biology works on an interface of biology and mathematics and intends to unravel the complex mechanisms between the biological systems and the inter and intra species dynamics using computational tools, and high-throughput technologies developed on algorithms, networks and complex connections to simulate cellular biological processes.

CONCLUSION

Computational strategies and modelling integrate and prioritize microbial-host interactions and may predict the conditions in which the fine-tuning attenuates. These microbial-host interactions and working mechanisms are important from the aspect of effective drug designing and fine- tuning the therapeutic interventions.

FPocketWeb: protein pocket hunting in a web browser

Detecting macromolecular (e.g., protein) cavities where small molecules bind is an early step in computer-aided drug discovery. Multiple pocket-detection algorithms have been developed over the past several decades. Among them, fpocket, created by Schmidtke and Le Guilloux, is particularly popular. Like many programs used in computational-biology research, fpocket requires users to download and install an executable file. That file must also be run via a command-line interface, further complicating use. An existing fpocket server application effectively addresses these challenges, but it requires users to upload their possibly proprietary structures to a third-party server. The FPocketWeb web app builds on this prior work. It runs the fpocket3 executable entirely in a web browser without requiring installation. The pocket-finding calculations occur on the user's computer rather than on a remote server. A working version of the open-source FPocketWeb app can be accessed free of charge from http://durrantlab.com/fpocketweb .

Molecular dynamics of the viral life cycle: progress and prospects

Molecular dynamics (MD) simulations across spatiotemporal resolutions are widely applied to study viruses and represent the central technique uniting the field of computational virology. We discuss the progress of MD in elucidating the dynamics of the viral life cycle, including the status of modeling intact extracellular virions and leveraging advanced simulations to mimic active life cycle processes. We further remark on the prospects of MD for continued contributions to the basic science characterization of viruses, especially given the increasing availability of high-quality experimental data and supercomputing power. Overall, integrative computational methods that are closely guided by experiments are unmatched in the level of detail they provide, enabling-now and in the future-new discoveries relevant to thwarting viral infection.

The Pathobiology of H7N3 Low and High Pathogenicity Avian Influenza Viruses from the United States Outbreak in 2020 Differs between Turkeys and Chickens

An outbreak caused by H7N3 low pathogenicity avian influenza virus (LPAIV) occurred in commercial turkey farms in the states of North Carolina (NC) and South Carolina (SC), United States in March of 2020. Subsequently, H7N3 high pathogenicity avian influenza virus (HPAIV) was detected on a turkey farm in SC. The infectivity, transmissibility, and pathogenicity of the H7N3 HPAIV and two LPAIV isolates, including one with a deletion in the neuraminidase (NA) protein stalk, were studied in turkeys and chickens. High infectivity [<2 log₁₀ 50% bird infectious dose (BID₅₀)] and transmission to birds exposed by direct contact were observed with the HPAIV in turkeys. In contrast, the HPAIV dose to infect chickens was higher than for turkeys (3.7 log₁₀ BID₅₀), and no transmission was observed. Similarly, higher infectivity (<2–2.5 log₁₀ BID₅₀) and transmissibility were observed with the H7N3 LPAIVs in turkeys compared to chickens, which required higher virus doses to become infected (5.4–5.7 log₁₀ BID₅₀). The LPAIV with the NA stalk deletion was more infectious in turkeys but did not have enhanced infectivity in chickens. These results show clear differences in the pathobiology of AIVs in turkeys and chickens and corroborate the high susceptibility of turkeys to both LPAIV and HPAIV infections.

Large-Scale Molecular Dynamics Simulations of Cellular Compartments

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.

Integrating molecular modelling methods to advance influenza A virus drug discovery

Since the discovery of the anti-influenza drugs oseltamivir and zanamivir using computer-aided drug design methods, there have been significant applications of molecular modelling methodologies applied to influenza A virus drug discovery, such as molecular dynamics (MD) simulation, molecular docking, and virtual screening (VS). In this review, we provide a brief general introduction to molecular modelling in the context of drug discovery and then focus on the advances and impact of integrating these methods with specific reference to potential influenza A antiviral drug targets.

H9N2 Influenza Virus Infections in Human Cells Require a Balance between Neuraminidase Sialidase Activity and Hemagglutinin Receptor Affinity

H9N2 avian influenza (AI) virus, one of the most prevalent AI viruses, has caused repeated poultry and human infections, posing a huge public health risk. The H9N2 virus has diversified into multiple lineages, with the G1 lineage being the most prevalent worldwide. In this study, we isolated G1 variants carrying an 8-amino-acid deletion in their NA stalk, which is, to our knowledge, the longest deletion found in H9N2 viruses in the field. The NA stalk length was found to modulate G1 virus entry into host cells, with the effects being species specific and dependent on the corresponding HA binding affinity. Our results suggest that, in nature, H9N2 G1 viruses balance their HA and NA functions by the NA stalk length, leading to the possible association of host range and virulence in poultry and mammals during the evolution of G1 lineage viruses.

Some avian influenza (AI) viruses have a deletion of up to 20 to 30 amino acids in their neuraminidase (NA) stalk. This has been associated with changes in virus replication and host range. Currently prevalent H9N2 AI viruses have only a 2- or 3-amino-acid deletion, and such deletions were detected in G1 and Y280 lineage viruses, respectively. The effect of an NA deletion on the H9N2 phenotype has not been fully elucidated. In this study, we isolated G1 mutants that carried an 8-amino-acid deletion in their NA stalk. To systematically analyze the effect of NA stalk length and concomitant (de)glycosylation on G1 replication and host range, we generated G1 viruses that had various NA stalk lengths and that were either glycosylated or not glycosylated. The stalk length was correlated with NA sialidase activity, using low-molecular-weight substrates, and with virus elution efficacy from erythrocytes. G1 virus replication in avian cells and eggs was positively correlated with the NA stalk length but was negatively correlated in human cells and mice. NA stalk length modulated G1 virus entry into host cells, with shorter stalks enabling more efficient G1 entry into human cells. However, with a hemagglutinin (HA) with a higher α2,6-linked sialylglycan affinity, the effect of NA stalk length on G1 virus infection was reversed, with shorter NA stalks reducing virus entry into human cells. These results indicate that a balance between HA binding affinity and NA sialidase activity, modulated by NA stalk length, is required for optimal G1 virus entry into human airway cells.

IMPORTANCE H9N2 avian influenza (AI) virus, one of the most prevalent AI viruses, has caused repeated poultry and human infections, posing a huge public health risk. The H9N2 virus has diversified into multiple lineages, with the G1 lineage being the most prevalent worldwide. In this study, we isolated G1 variants carrying an 8-amino-acid deletion in their NA stalk, which is, to our knowledge, the longest deletion found in H9N2 viruses in the field. The NA stalk length was found to modulate G1 virus entry into host cells, with the effects being species specific and dependent on the corresponding HA binding affinity. Our results suggest that, in nature, H9N2 G1 viruses balance their HA and NA functions by the NA stalk length, leading to the possible association of host range and virulence in poultry and mammals during the evolution of G1 lineage viruses.

EasyAmber: A comprehensive toolbox to automate the molecular dynamics simulation of proteins

Conformational plasticity of the functionally important regions and binding sites in protein/enzyme structures is one of the key factors affecting their function and interaction with substrates/ligands. Molecular dynamics (MD) can address the challenge of accounting for protein flexibility by predicting the time-dependent behavior of a molecular system. It has a potential of becoming a particularly important tool in protein engineering and drug discovery, but requires specialized training and skills, what impedes practical use by many investigators. We have developed the easyAmber - a comprehensive set of programs to automate the molecular dynamics routines implemented in the Amber package. The toolbox can address a wide set of tasks in computational biology struggling to account for protein flexibility. The automated workflow includes a complete set of steps from the initial "static" molecular model to the MD "production run": the full-atom model building, optimization/equilibration of the molecular system, classical/conventional and accelerated molecular dynamics simulations. The easyAmber implements advanced MD protocols, but is highly automated and easy-to-operate to attract a broad audience. The toolbox can be used on a personal desktop station equipped with a compatible gaming GPU-accelerator, as well as help to manage huge workloads on a powerful supercomputer. The software provides an opportunity to operate multiple simulations of different proteins at the same time, thus significantly increasing work efficiency. The easyAmber takes the molecular dynamics to the next level in terms of usability for complex processing of large volumes of data, thus supporting the recent trend away from inefficient "static" approaches in biology toward a deeper understanding of the dynamics in protein structures. The software is freely available for download at https://biokinet.belozersky.msu.ru/easyAmber, no login required.

Effect of the selection pressure of vaccine antibodies on evolution of H9N2 avian influenza virus in chickens

H9N2 avian influenza virus has spread worldwide, and vaccination with an inactivated virus is currently the major prevention method in China. To further understand the effect of the selection pressure from antibodies on the evolution of H9N2 avian influenza virus, F/98 (A/Chicken/Shanghai/F/98), which is the vaccine representative of H9N2 avian influenza virus in East China, was used for serial passaging for 20 generations in chickens with and without vaccination. After plaque purification from trachea and lung tissues, 390 quasispecies were obtained. The second-generation quasispecies under the selection pressure of vaccine antibodies had undergone 100% antigen variation, while after passaging to the fifth generation, only 30-40% of the quasispecies displayed antigen variation when there was no selection pressure of vaccine antibodies, implying that the selection pressure of vaccine antibodies promotes the antigen variation of F/98. We found for the first time that there were three mutation hotspots in the HA genes of the quasispecies under the selection pressure of vaccine antibodies, which were K131R, A168T, and N201D. Moreover, under the selection pressure of vaccine antibodies, 10 amino acids (67-76) of the NA protein of all quasispecies were deleted, and PB2 of the quasispecies had undergone a high-frequency R355K mutation. However, without selection pressure of vaccine antibodies, NP had undergone two high-frequency mutations, namely, V186I and L466I, and a high-frequency mutation of L77I appeared in the NS gene. This result shows that the vaccine antibody selection pressure could control and regulate gene variation of the F/98 virus. Compared to that of the parental virus F/98, the EID50 of the twentieth passaged virus under the selection pressure of vaccine antibodies did not change, while the EID50 of the twentieth passaged virus without selection pressure of vaccine antibodies was significantly enhanced by 794 times. Furthermore, the twentieth passaged virus with selection pressure from vaccine antibodies lost its lethal ability in embryonated chicken eggs, whereas the EID50 of the twentieth passaged virus without selection pressure of vaccine antibodies increased to 6.3 times that of the F/98 strain. All the above results show that the selection pressure of vaccine antibodies promotes the antigen variation of H9N2 avian influenza virus and plays a role in regulating and controlling gene mutation of H9N2 avian influenza virus.

Mesoscale All-Atom Influenza Virus Simulations Suggest New Substrate Binding Mechanism

Influenza virus circulates in human, avian, and swine hosts, causing seasonal epidemic and occasional pandemic outbreaks. Influenza neuraminidase, a viral surface glycoprotein, has two sialic acid binding sites. The catalytic (primary) site, which also binds inhibitors such as oseltamivir carboxylate, is responsible for cleaving the sialic acid linkages that bind viral progeny to the host cell. In contrast, the functional annotation of the secondary site remains unclear. Here, we better characterize these two sites through the development of an all-atom, explicitly solvated, and experimentally based integrative model of the pandemic influenza A H1N1 2009 viral envelope, containing ∼160 million atoms and spanning ∼115 nm in diameter. Molecular dynamics simulations of this crowded subcellular environment, coupled with Markov state model theory, provide a novel framework for studying realistic molecular systems at the mesoscale and allow us to quantify the kinetics of the neuraminidase 150-loop transition between the open and closed states. An analysis of chloride ion occupancy along the neuraminidase surface implies a potential new role for the neuraminidase secondary site, wherein the terminal sialic acid residues of the linkages may bind before transfer to the primary site where enzymatic cleavage occurs. Altogether, our work breaks new ground for molecular simulation in terms of size, complexity, and methodological analyses of the components. It also provides fundamental insights into the understanding of substrate recognition processes for this vital influenza drug target, suggesting a new strategy for the development of anti-influenza therapeutics.

Infectious salmon anaemia virus-molecular biology and pathogenesis of the infection

Aquaculture has a long history in many parts of the world, but it is still young at an industrial scale. Marine fish farming in open nets of a single fish species at high densities compared to their wild compatriots opens a plethora of possible infections. Infectious salmon anaemia (ISA) is an example of disease that surfaced after large-scale farming of Atlantic salmon (Salmo salar) appeared. Here, a review of the molecular biology of the ISA virus (ISAV) with emphasis on its pathogenicity is presented. The avirulent HPR0 variant of ISAV has resisted propagation in cell cultures, which has restricted the ability to perform in vivo experiments with this variant. The transition from avirulent HPR0 to virulent HPRΔ has not been methodically studied under controlled experimental conditions, and the triggers of the transition from avirulent to virulent forms have not been mapped. Genetic segment reassortment, recombination and mutations are important mechanisms in ISAV evolution, and for the development of virulence. In the 25 years since the ISAV was identified, large amounts of sequence data have been collected for epidemiologic and transmission studies, however, the lack of good experimental models for HPR0 make the risk evaluation of the presence of this avirulent, ubiquitous variant uncertain. This review summarizes the current knowledge related to molecular biology and pathogenicity of this important aquatic orthomyxovirus.

Diversity of sialidases found in the human body - A review

Sialidases are enzymes essential for numerous organisms including humans. Hydrolytic sialidases (EC 3.2.1.18), trans-sialidases and anhydrosialidases (intramolecular trans-sialidases, EC 4.2.2.15) are glycoside hydrolase enzymes that cleave the glycosidic linkage and release sialic acid residues from sialyl substrates. The paper summarizes diverse sialidases present in the human body and their potential impact on development of antiviral compounds - inhibitors of viral neuraminidases. It includes a brief overview of catalytic mechanisms of action of sialidases and describes the origin of sialidases in the human body. This is followed by description of the structure and function of sialidase families with a special focus on the GH33 and GH34 families. Various effects of sialidases on human body are also briefly described. Modulation of sialidase activity may be considered a useful tool for effective treatment of various diseases. In some cases, it is desired to completely suppress the activity of sialidases by suitable inhibitors. Specific sialidase inhibitors are useful for the treatment of influenza, epilepsy, Alzheimer's disease, diabetes, different types of cancer, or heart defects. Challenges and future directions are shortly depicted in the final part of the paper.

Influenza A Virus Hemagglutinin-Neuraminidase-Receptor Balance: Preserving Virus Motility

Influenza A viruses (IAVs) occasionally cross the species barrier and adapt to novel host species. This requires readjustment of the functional balance of the sialic acid receptor-binding hemagglutinin (HA) and the receptor-destroying neuraminidase (NA) to the sialoglycan-receptor repertoire of the new host. Novel techniques have revealed mechanistic details of this HA-NA-receptor balance, emphasizing a previously underappreciated crucial role for NA in driving the motility of receptor-associated IAV particles. Motility enables virion penetration of the sialylated mucus layer as well as attachment to, and uptake into, underlying epithelial cells. As IAVs are essentially irreversibly bound in the absence of NA activity, the fine-tuning of the HA-NA-receptor balance rather than the binding avidity of IAV particles per se is an important factor in determining host species tropism.

Building a macro-mixing dual-basin Gō model using the Multistate Bennett Acceptance Ratio

The dual-basin Gō-model is a structural-based coarsegrained model for simulating a conformational transition between two known structures of a protein. Two parameters are required to produce a dual-basin potential mixed using two single-basin potentials, although the determination of mixing parameters is usually not straightforward. Here, we have developed an efficient scheme to determine the mixing parameters using the Multistate Bennett Acceptance Ratio (MBAR) method after short simulations with a set of parameters. In the scheme, MBAR allows us to predict observables at various unsimulated conditions, which are useful to improve the mixing parameters in the next round of iterative simulations. The number of iterations that are necessary for obtaining the converged mixing parameters are significantly reduced in the scheme. We applied the scheme to two proteins, the glutamine binding protein and the ribose binding protein, for showing the effectiveness in the parameter determination. After obtaining the converged parameters, both proteins show frequent conformational transitions between open and closed states, providing the theoretical basis to investigate structure-dynamics-function relationships of the proteins.

PCAViz: An Open-Source Python/JavaScript Toolkit for Visualizing Molecular Dynamics Simulations in the Web Browser

Molecular dynamics (MD) simulations reveal molecular motions at atomic resolution. Recent advances in high-performance computing now enable microsecond-long simulations capable of sampling a wide range of biologically relevant events. But the disk space required to store an MD trajectory increases with simulation length and system size, complicating collaborative sharing and visualization. To overcome these limitations, we created PCAViz, an open-source toolkit for sharing and visualizing MD trajectories via the web browser. PCAViz includes two components: the PCAViz Compressor, which compresses and saves simulation data; and the PCAViz Interpreter, which decompresses the data in users’ browsers and feeds it to any of several browser-based molecular-visualization libraries (e.g., 3Dmol.js, NGL Viewer, etc.). An easy-to-install WordPress plugin enables “plug-and-play” trajectory visualization. PCAViz will appeal to a broad audience of researchers and educators. The source code is available at http://durrantlab.com/pcaviz/, and the WordPress plugin is available via the official WordPress Plugin Directory.

Extending the Stalk Enhances Immunogenicity of the Influenza Virus Neuraminidase

Influenza viruses express two surface glycoproteins, the hemagglutinin (HA) and the neuraminidase (NA). Anti-NA antibodies protect from lethal influenza virus challenge in the mouse model and correlate inversely with virus shedding and symptoms in humans. Consequently, the NA is a promising target for influenza virus vaccine design. Current seasonal vaccines, however, poorly induce anti-NA antibodies, partly because of the immunodominance of the HA over the NA when the two glycoproteins are closely associated. To address this issue, here we investigated whether extending the stalk domain of the NA could render it more immunogenic on virus particles. Two recombinant influenza viruses based on the H1N1 strain A/Puerto Rico/8/1934 (PR8) were rescued with NA stalk domains extended by 15 or 30 amino acids. Formalin-inactivated viruses expressing wild-type NA or the stalk-extended NA variants were used to vaccinate mice. The virus with the 30-amino-acid stalk extension induced significantly higher anti-NA IgG responses (characterized by increased in vitro antibody-dependent cellular cytotoxicity [ADCC] activity) than the wild-type PR8 virus, while anti-HA IgG levels were unaffected. Similarly, extending the stalk domain of the NA of a recent H3N2 virus enhanced the induction of anti-NA IgGs in mice. On the basis of these results, we hypothesize that the subdominance of the NA can be modulated if the protein is modified such that its height surpasses that of the HA on the viral membrane. Extending the stalk domain of NA may help to enhance its immunogenicity in influenza virus vaccines without compromising antibody responses to HA.IMPORTANCE The efficacy of influenza virus vaccines could be improved by enhancing the immunogenicity of the NA protein. One of the reasons for its poor immunogenicity is the immunodominance of the HA over the NA in many seasonal influenza virus vaccines. Here we demonstrate that, in the mouse model, extending the stalk domain of the NA protein can enhance its immunogenicity on virus particles and overcome the immunodominance of the HA without affecting antibody responses to the HA. The antibody repertoire is broadened by the extended NA and includes additional ADCC-active antibodies. Our findings may assist in the efforts toward more effective influenza virus vaccines.

Evidence of exposure and human seroconversion during an outbreak of avian influenza A(H5N1) among poultry in Cameroon

From May 2016 to March 2017, 22 poultry outbreaks of avian influenza A(H5N1) were reported in Cameroon, mainly in poultry farms and live bird markets. No human cases were reported. In this study, we sought to describe the 2016 A(H5N1) outbreak strain and to investigate the risk of infection in exposed individuals. We find that highly pathogenic influenza subtype A(H5N1), clade 2.3.2.1c from Cameroon is closely related phylogenetically and antigenically to strains isolated in central and western Africa at the time. No molecular markers of increased human transmissibility were noted; however, seroconversion was detected in two poultry workers (1.5% of total screened). Therefore, the continued outbreaks of avian influenza in poultry and the risk of zoonotic human infection highlight the crucial need for continued and vigilant influenza surveillance and research in Africa, especially in areas of high poultry trade, such as Cameroon.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Influenza Hemagglutinin and Neuraminidase: Yin⁻Yang Proteins Coevolving to Thwart Immunity

Influenza A virions possess two surface glycoproteins-the hemagglutinin (HA) and neuraminidase (NA)-which exert opposite functions. HA attaches virions to cells by binding to terminal sialic acid residues on glycoproteins/glycolipids to initiate the infectious cycle, while NA cleaves terminal sialic acids, releasing virions to complete the infectious cycle. Antibodies specific for HA or NA can protect experimental animals from IAV pathogenesis and drive antigenic variation in their target epitopes that impairs vaccine effectiveness in humans. Here, we review progress in understanding HA/NA co-evolution as each acquires epistatic mutations to restore viral fitness to mutants selected in the other protein by host innate or adaptive immune pressure. We also discuss recent exciting findings that antibodies to HA can function in vivo by blocking NA enzyme activity to prevent nascent virion release and enhance Fc receptor-based activation of innate immune cells.

Influenza Virus Neuraminidase Structure and Functions

With the constant threat of emergence of a novel influenza virus pandemic, there must be continued evaluation of the molecular mechanisms that contribute to virulence. Although the influenza A virus surface glycoprotein neuraminidase (NA) has been studied mainly in the context of its role in viral release from cells, accumulating evidence suggests it plays an important, multifunctional role in virus infection and fitness. This review investigates the various structural features of NA, linking these with functional outcomes in viral replication. The contribution of evolving NA activity to viral attachment, entry and release of virions from infected cells, and maintenance of functional balance with the viral hemagglutinin are also discussed. Greater insight into the role of this important antiviral drug target is warranted.

A Computational Assay that Explores the Hemagglutinin/Neuraminidase Functional Balance Reveals the Neuraminidase Secondary Site as a Novel Anti-Influenza Target

Studies of pathogen-host specificity, virulence, and transmissibility are critical for basic research as well as for assessing the pandemic potential of emerging infectious diseases. This is especially true for viruses such as influenza, which continue to affect millions of people annually through both seasonal and occasional pandemic events. Although the influenza virus has been fairly well studied for decades, our understanding of host-cell binding and its relation to viral transmissibility and infection is still incomplete. Assessing the binding mechanisms of complex biological systems with atomic-scale detail is challenging given current experimental limitations. Much remains to be learned, for example, about how the terminal residue of influenza-binding host-cell receptors (sialic acid) interacts with the viral surface. Here, we present an integrative structural-modeling and physics-based computational assay that reveals the sialic acid association rate constants (k _on) to three influenza sites: the hemagglutinin (HA), neuraminidase (NA) active, and NA secondary binding sites. We developed a series of highly detailed (atomic-resolution) structural models of fully intact influenza viral envelopes. Brownian dynamics simulations of these systems showed how structural properties, such as stalk height and secondary-site binding, affect sialic acid k _on values. Comparing the k _on values of the three sialic acid binding sites across different viral strains suggests a detailed model of encounter-complex formation and indicates that the secondary NA binding site may play a compensatory role in host-cell receptor binding. Our method elucidates the competition among these sites, all present on the same virion, and provides a new technology for directly studying the functional balance between HA and NA.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

Structural Changes Due to Antagonist Binding in Ligand Binding Pocket of Androgen Receptor Elucidated Through Molecular Dynamics Simulations

When a small molecule binds to the androgen receptor (AR), a conformational change can occur which impacts subsequent binding of co-regulator proteins and DNA. In order to accurately study this mechanism, the scientific community needs a crystal structure of the Wild type AR (WT-AR) ligand binding domain, bound with antagonist. To address this open need, we leveraged molecular docking and molecular dynamics (MD) simulations to construct a structure of the WT-AR ligand binding domain bound with antagonist bicalutamide. The structure of mutant AR (Mut-AR) bound with this same antagonist informed this study. After molecular docking analysis pinpointed the suitable binding orientation of a ligand in AR, the model was further optimized through 1 μs of MD simulations. Using this approach, three molecular systems were studied: (1) WT-AR bound with agonist R1881, (2) WT-AR bound with antagonist bicalutamide, and (3) Mut-AR bound with bicalutamide. Our structures were very similar to the experimentally determined structures of both WT-AR with R1881 and Mut-AR with bicalutamide, demonstrating the trustworthiness of this approach. In our model, when WT-AR is bound with bicalutamide, Val716/Lys720/Gln733, or Met734/Gln738/Glu897 move and thus disturb the positive and negative charge clumps of the AF2 site. This disruption of the AF2 site is key for understanding the impact of antagonist binding on subsequent co-regulator binding. In conclusion, the antagonist induced structural changes in WT-AR detailed in this study will enable further AR research and will facilitate AR targeting drug discovery.

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the "open" hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine-rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.

Structural comparison among the 2013-2017 avian influenza A H5N6 hemagglutinin proteins: A computational study with epidemiological implications

Avian influenza viruses easily spread allowing viral re-assortment to simply occur which in-turn increases the potential for a pandemic. A novel 2013 H5N6 influenza strain was detected among the avian population and was reported to continuously evolve, however, this was never structurally demonstrated. Here, we elucidated the putative structural evolution of the novel H5N6 influenza strain. Throughout this study, we analyzed 2013-2017 H5N6 HA protein models. Model quality was first verified before further analyses and structural comparison was made using superimposition. We found that Leu was inserted at position 129₁ among the 2013-2015 models while Leu was not inserted among the 2016-2017 models. Moreover, presence of Leu at position 129₁ shifts residue E126₁ by 159.6° affecting nearby residues which may explain the difference between the 2013-2015 and 2016-2017 HA structural groups. Similarly, we believe that our results would support the hypothesis that the current H5N6 strain is still continuously evolving.

Adaptive mutations of neuraminidase stalk truncation and deglycosylation confer enhanced pathogenicity of influenza A viruses

It has been noticed that neuraminidase (NA) stalk truncation has arisen from evolutionary adaptation of avian influenza A viruses (IAVs) from wild aquatic birds to domestic poultry. We identified this molecular alteration after the adaptation of a 2009 pandemic H1N1 virus (pH1N1) in BALB/c mice. The mouse-adapted pH1N1 lost its eight consecutive amino acids including one potential N-linked glycosite from the NA stalk region. To explore the relationship of NA stalk truncation or deglycosylation with viral pathogenicity changes, we generated NA stalk mutant viruses on the pH1N1 backbone by reverse genetics. Intriguingly, either NA stalk truncation or deglycosylation changed pH1N1 into a lethal virus to mice by resulting in extensive pathologic transformation in the mouse lungs and systemic infection affecting beyond the respiratory organs in mice. The increased pathogenicity of these NA stalk mutants was also reproduced in ferrets. In further investigation using a human-infecting H7N9 avian IAV strain, NA stalk truncation or deglycosylation enhanced the replication property and pathogenicity of H7N9 NA stalk mutant viruses in the same mouse model. Taken together, our results suggest that NA stalk truncation or deglycosylation can be the pathogenic determinants of seasonal influenza viruses associated with the evolutionary adaptation of IAVs.

Genetic Characterization of Continually Evolving Highly Pathogenic H5N6 Influenza Viruses in China, 2012-2016

H5N6 is a highly pathogenic avian influenza (HPAI) and a zoonotic disease that causes recurring endemics in East Asia. At least 155 H5N6 outbreaks, including 15 human infections, have been reported in China. These repeated outbreaks have increased concern that the H5N6 virus may cross over to humans and cause a pandemic. In February, 2016, peafowls in a breeding farm exhibited a highly contagious disease. Post-mortem examinations, including RT-PCR, and virus isolation, confirmed that the highly pathogenic H5N6 influenza virus was the causative agent, and the strain was named A/Pavo Cristatus/Jiangxi/JA1/2016. In animal experiments, it exhibited high pathogenicity in chickens and an estimated median lethal dose in mice of ~10^4.3 TCID₅₀. A phylogenetic analysis showed that JA1/2016 was clustered in H5 clade 2.3.4.4. FG594-like H5N6 virus from Guangdong Province was the probable predecessor of JA1/2016, and the estimated divergence time was June 2014. Furthermore, we found that H5N6 influenza viruses can be classified into the two following groups: Group 1 and Group 2. Group 2 influenza viruses have not been detected since the end of 2014, whereas Group 1 influenza viruses have continually evolved and reassorted with the “gene pool” circulating in south China, resulting in the rise of novel subtypes of this influenza virus. An increase in the number of its identified hosts, the expanding range of its distribution, and the continual evolution of H5N6 AIVs enhance the risk that an H5N6 virus may spread to other continents and cause a pandemic.