A molecular modelling approach to understand the effect of co-evolutionary mutations (V344M, I354L) identified in the PB2 subunit of influenza A 2009 pandemic H1N1 virus on m7GTP ligand binding

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.

Inter- Versus Intra-Host Sequence Diversity of pH1N1 and Associated Clinical Outcomes

The diversity of RNA viruses dictates their evolution in a particular host, community or environment. Here, we reported within- and between-host pH1N1virus diversity at consensus and sub-consensus levels over a three-year period (2015-2017) and its implications on disease severity. A total of 90 nasal samples positive for the pH1N1 virus were deep-sequenced and analyzed to detect low-frequency variants (LFVs) and haplotypes. Parallel evolution of LFVs was seen in the hemagglutinin (HA) gene across three scales: among patients (33%), across years (22%), and at global scale. Remarkably, investigating the emergence of LFVs at the consensus level demonstrated that within-host virus evolution recapitulates evolutionary dynamics seen at the global scale. Analysis of virus diversity at the HA haplotype level revealed the clustering of low-frequency haplotypes from early 2015 with dominant strains of 2016, indicating rapid haplotype evolution. Haplotype sharing was also noticed in all years, strongly suggesting haplotype transmission among patients infected during a specific influenza season. Finally, more than half of patients with severe symptoms harbored a larger number of haplotypes, mostly in patients under the age of five. Therefore, patient age, haplotype diversity, and the presence of certain LFVs should be considered when interpreting illness severity. In addition to its importance in understanding virus evolution, sub-consensus virus diversity together with whole genome sequencing is essential to explain variabilities in clinical outcomes that cannot be explained by either analysis alone.

Computational analysis of the effect of polymerase acidic (PA) gene mutation F35L in the 2009 pandemic influenza A (H1N1) virus on binding aspects of mononucleotides in the endonuclease domain

An F35L mutation in the N-terminal domain of the polymerase acidic protein (PA-Nter), which contains the active site of the endonuclease, has been reported to result in higher polymerase activity in mouse-adapted strains of the 2009 pandemic influenza A H1N1 virus. We modeled wild and mutant complexes of uridine 5'-monophosphate (UMP) as the endonuclease substrate and performed molecular dynamics simulations. The results demonstrated that the F35L mutation could result in a changed orientation of a helix containing active site residues and improve the ligand affinity in the mutant strain. This study suggests a molecular mechanism of enhanced polymerase activity.

Hydrogen Bond Variations of Influenza A Viruses During Adaptation in Human

Many host specific mutations have been detected in influenza A viruses (IAVs). However, their effects on hydrogen bond (H-bond) variations have rarely been investigated. In this study, 60 host specific sites were identified in the internal proteins of avian and human IAVs, 27 of which contained mutations with effects on H-bonds. Besides, 30 group specific sites were detected in HA and NA. Twenty-six of 36 mutations existing at these group specific sites caused H-bond loss or formation in at least one subtype. The number of mutations in isolations of 2009 pandemic H1N1, human-infecting H5N1 and H7N9 varied. The combinations of mutations and H-bond changes in these three subtypes of IAVs were also different. In addition, the mutations in isolations of H5N1 distributed more scattered than those in 2009 pandemic H1N1 and H7N9. Eight wave specific mutations in isolations of the fifth H7N9 wave were also identified. Three of them, R140K in HA, Y170H in NA, and R340K in PB2, were capable of resulting in H-bond loss. As mentioned above, these host or group or wave specific H-bond variations provide us with a new field of vision for understanding the changes of structural features in the human adaptation of IAVs.