Insight into the oseltamivir resistance R292K mutation in H5N1 influenza virus: a molecular docking and molecular dynamics approach

Characterization of the In Vitro and In Vivo Efficacy of Baloxavir Marboxil against H5 Highly Pathogenic Avian Influenza Virus Infection

Human infections caused by the H5 highly pathogenic avian influenza virus (HPAIV) sporadically threaten public health. The susceptibility of HPAIVs to baloxavir acid (BXA), a new class of inhibitors for the influenza virus cap-dependent endonuclease, has been confirmed in vitro, but it has not yet been fully characterized. Here, the efficacy of BXA against HPAIVs, including recent H5N8 variants, was assessed in vitro. The antiviral efficacy of baloxavir marboxil (BXM) in H5N1 virus-infected mice was also investigated. BXA exhibited similar in vitro activities against H5N1, H5N6, and H5N8 variants tested in comparison with seasonal and other zoonotic strains. Compared with oseltamivir phosphate (OSP), BXM monotherapy in mice infected with the H5N1 HPAIV clinical isolate, the A/Hong Kong/483/1997 strain, also caused a significant reduction in viral titers in the lungs, brains, and kidneys, thereby preventing acute lung inflammation and reducing mortality. Furthermore, compared with BXM or OSP monotherapy, combination treatments with BXM and OSP using a 48-h delayed treatment model showed a more potent effect on viral replication in the organs, accompanied by improved survival. In conclusion, BXM has a potent antiviral efficacy against H5 HPAIV infections.

Emergence of Novel Reassortant H1N1 Avian Influenza Viruses in Korean Wild Ducks in 2018 and 2019

Influenza A virus subtype H1N1 has caused global pandemics like the “Spanish flu” in 1918 and the 2009 H1N1 pandemic several times. H1N1 remains in circulation and survives in multiple animal sources, including wild birds. Surveillance during the winter of 2018–2019 in Korea revealed two H1N1 isolates in samples collected from wild bird feces: KNU18-64 (A/Greater white-fronted goose/South Korea/KNU18-64/2018(H1N1)) and WKU19-4 (A/wild bird/South Korea/WKU19-4/2019(H1N1)). Phylogenetic analysis indicated that M gene of KNU18-64(H1N1) isolate resembles that of the Alaskan avian influenza virus, whereas WKU19-4(H1N1) appears to be closer to the Mongolian virus. Molecular characterization revealed that they harbor the amino acid sequence PSIQRS↓GLF and are low-pathogenicity influenza viruses. In particular, the two isolates harbored three different mutation sites, indicating that they have different virulence characteristics. The mutations in the PB1-F2 and PA protein of WKU19-4(H1N1) indicate increasing polymerase activity. These results corroborate the kinetic growth data for WKU19-4 in MDCK cells: a dramatic increase in the viral titer after 12 h post-inoculation compared with that in the control group H1N1 (CA/04/09(pdm09)). The KNU18-64(H1N1) isolate carries mutations indicating an increase in mammal adaptation; this characterization was confirmed by the animal study in mice. The KNU18-64(H1N1) group showed the presence of viruses in the lungs at days 3 and 6 post-infection, with titers of 2.71 ± 0.16 and 3.71 ± 0.25 log10(TCID50/mL), respectively, whereas the virus was only detected in the WKU19-4(H1N1) group at day 6 post-infection, with a lower titer of 2.75 ± 0.51 log10(TCID50/mL). The present study supports the theory that there is a relationship between Korea and America with regard to reassortment to produce novel viral strains. Therefore, there is a need for increased surveillance of influenza virus circulation in free-flying and wild land-based birds in Korea, particularly with regard to Alaskan and Asian strains.

E-pharmacophore hypothesis strategy to discover potent inhibitor for influenza treatment

The surface protein of Influenza virus, Neuraminidase (NA), is believed to play a critical role in the release of new viral particle and thus spreads infection. It has been recognized as a valid drug target for anti-influenza therapy. Despite the number of available approved drugs for the influenza infection treatment, the emergence of resistant variants with novel mutations are the foremost challenges for the currently used NA inhibitors. Thus, the current investigation was carried out to ascertain potent inhibitors using computational strategies such as e-pharmacophore based virtual screening and docking approach. A three-dimensional e-pharmacophore hypothesis was generated based on the chemical features of complexes of the drugs and NA protein using PHASE module of Schrödinger suite. The generated hypothesis consisted of one hydrogen bond acceptor (A), two hydrogen bond donors (D), one negatively charged group (N) and one aromatic ring (R), ADDNR. The hypothesis was further evaluated for its integrity using enrichment analysis and used to filter out molecules with similar pharmacophoric features from approved, investigational and experimental subsets of DrugBank and ZINC database. In addition, ligand filtration was performed to curb down the molecules to an efficient collection of hit molecules by using Lipinski “rule of five and ADME analysis by using Qikprop module. Overall, the results from our analysis suggest that compound lisinopril and formoterol could serve as potent antiviral compounds for the treatment of influenza A virus infection. It is worth mentioning that the results correlate well with literature evidences.

Hyphenated 3D-QSAR statistical model-drug repurposing analysis for the identification of potent neuraminidase inhibitor

The Influenza A virus is one of the principle causes of respiratory illness in human. The surface glycoprotein of the influenza virus, neuraminidase (NA), has a vital role in the release of new viral particle and spreads infection in the respiratory tract. It has been long recognized as a valid drug target for influenza A virus infection. Oseltamivir is used as a standard drug of choice for the treatment of influenza. However, the emergence of mutants with novel mutations has increased the resistance to potent NA inhibitor. In the present investigation, we have employed computer-assisted combinatorial techniques in the screening of 8621 molecules from Drug Bank to find potent NA inhibitors. A three-dimensional pharmacophore model was generated from the previously reported 28 carbocylic influenza NA inhibitors along with oseltamivir using PHASE module of Schrödinger Suite. The model generated consists of one hydrogen bond acceptor (A), one hydrogen bond donors (D), one hydrophobic group (H), and one positively charged group (P), ADHP. The hypothesis was further validated for its integrity and significance using enrichment analysis. Subsequently, an atom-based 3D-QSAR model was built using the common pharmacophore hypothesis (CPH). The developed 3D-QSAR model was found to be statistically significant with R² value of 0.9866 and Q² value of 0.7629. Further screening was accomplished using three-stage docking process using the Glide algorithm. The resultant lead molecules were examined for its drug-like properties using the Qikprop algorithm. Finally, the calculated pIC₅₀ values of the lead compounds were validated by the AutoQSAR algorithm. Overall, the results from our analysis highlights that lisinopril (DB00722) is predicted to bind better with NA than currently approved drug. In addition, it has the best match in binding geometry conformations with the existing NA inhibitor. Note that the antiviral activity of lisinopril is reported in the literature. However, our paper is the first report on lisinopril activity against influenza A virus infection. These results are envisioned to help design the novel NA inhibitors with an increased antiviral efficacy.

The inhibitory performance of flavonoid cyanidin-3-sambubiocide against H274Y mutation in H1N1 influenza virus

Oseltamivir (Tamiflu) is the most accepted antiviral drug that targets the neuraminidase (NA) protein to inhibit the viral release from the host cell. Few H1N1 influenza strains with the H274Y mutation creates drug resistance to oseltamivir. In this study, we report that flavonoid cyanidin-3-sambubiocide (C3S) compound acts as a potential inhibitor against H274Y mutation. The drug resistance mechanism and inhibitory activity of C3S and oseltamivir against wild-type (WT) and H274Y mutant-type (MT) have been studied and compared based on the results of molecular docking, molecular dynamics, and quantum chemical methods. Oseltamivir has been found less binding affinity with MT. C3S has more binding affinity with WT and MT proteins. From the dynamical study, the 150th loop of the MT protein has found more deformation than WT. A single H274Y mutation induces the conformational changes in the 150th loop which leads to produce more resistance to oseltamivir. The 150th cavity is more attractive target for C3S to stop the conformational changes in the MT, than 430th cavity of NA protein. The C3S is stabilized with MT by more number of hydrogen bonds than oseltamivir. The electrostatic interaction energy shows a stronger C3S binding with MT and this compound may be more effective against oseltamivir-resistant virus strains.

Exploring the Origin of Differential Binding Affinities of Human Tubulin Isotypes αβII, αβIII and αβIV for DAMA-Colchicine Using Homology Modelling, Molecular Docking and Molecular Dynamics Simulations

Tubulin isotypes are found to play an important role in regulating microtubule dynamics. The isotype composition is also thought to contribute in the development of drug resistance as tubulin isotypes show differential binding affinities for various anti-cancer agents. Tubulin isotypes αβII, αβIII and αβIV show differential binding affinity for colchicine. However, the origin of differential binding affinity is not well understood at the molecular level. Here, we investigate the origin of differential binding affinity of a colchicine analogue N-deacetyl-N-(2-mercaptoacetyl)-colchicine (DAMA-colchicine) for human αβII, αβIII and αβIV isotypes, employing sequence analysis, homology modeling, molecular docking, molecular dynamics simulation and MM-GBSA binding free energy calculations. The sequence analysis study shows that the residue compositions are different in the colchicine binding pocket of αβII and αβIII, whereas no such difference is present in αβIV tubulin isotypes. Further, the molecular docking and molecular dynamics simulations results show that residue differences present at the colchicine binding pocket weaken the bonding interactions and the correct binding of DAMA-colchicine at the interface of αβII and αβIII tubulin isotypes. Post molecular dynamics simulation analysis suggests that these residue variations affect the structure and dynamics of αβII and αβIII tubulin isotypes, which in turn affect the binding of DAMA-colchicine. Further, the binding free-energy calculation shows that αβIV tubulin isotype has the highest binding free-energy and αβIII has the lowest binding free-energy for DAMA-colchicine. The order of binding free-energy for DAMA-colchicine is αβIV ≃ αβII >> αβIII. Thus, our computational approaches provide an insight into the effect of residue variations on differential binding of αβII, αβIII and αβIV tubulin isotypes with DAMA-colchicine and may help to design new analogues with higher binding affinities for tubulin isotypes.

The significance of naturally occurring neuraminidase quasispecies of H5N1 avian influenza virus on resistance to oseltamivir: a point of concern

Viral adaptability and survival arise due to the presence of quasispecies populations that are able to escape the immune response or produce drug-resistant variants. However, the presence of H5N1 virus with natural mutations acquired without any drug selection pressure poses a great threat. Cloacal samples collected from the 2004-2005 epidemics in Thailand from Asian open-billed storks revealed one major and several minor quasispecies populations with mutations on the oseltamivir (OTV)-binding site of the neuraminidase gene (NA) without prior exposure to a drug. Therefore, this study investigated the binding between the NA-containing novel mutations and OTV drug using molecular dynamic simulations and plaque inhibition assay. The results revealed that the mutant populations, S236F mutant, S236F/C278Y mutant, A250V/V266A/P271H/G285S mutant and C278Y mutant, had a lower binding affinity with OTV as compared with the WT virus due to rearrangement of amino acid residues and increased flexibility in the 150-loop. This result was further emphasized through the IC50 values obtained for the major population and WT virus, 104.74 nM and 18.30 nM, respectively. Taken together, these data suggest that H5N1 viruses isolated from wild birds have already acquired OTV-resistant point mutations without any exposure to a drug.

Understanding the cross-resistance of oseltamivir to H1N1 and H5N1 influenza A neuraminidase mutations using multidimensional computational analyses

This study embarks on a comprehensive description of the conformational contributions to resistance of neuraminidase (N1) in H1N1 and H5N1 to oseltamivir, using comparative multiple molecular dynamic simulations. The available data with regard to elucidation of the mechanism of resistance as a result of mutations in H1N1 and H5N1 neuraminidases is not well established. Enhanced post-dynamic analysis, such as principal component analysis, solvent accessible surface area, free binding energy calculations, and radius of gyration were performed to gain a precise insight into the binding mode and origin of resistance of oseltamivir in H1N1 and H5N1 mutants. Three significant features reflecting resistance in the presence of mutations H274Y and I222K, of the protein complexed with the inhibitor are: reduced flexibility of the α-carbon backbone; an improved ΔE_ele of ~15 (kcal/mol) for H1N1 coupled with an increase in ΔG_sol (~13 kcal/mol) from wild-type to mutation; a low binding affinity in comparison with the wild-type of ~2 (kcal/mol) and ~7 (kcal/mol) with respect to each mutation for the H5N1 systems; and reduced hydrophobicity of the overall surface structure due to an impaired hydrogen bonding network. We believe the results of this study will ultimately provide a useful insight into the structural landscape of neuraminidase-associated binding of oseltamivir. Furthermore, the results can be used in the design and development of potent inhibitors of neuraminidases.

Single H5N1 influenza A neuraminidase mutation develops resistance to oseltamivir due to distorted conformational and drug binding landscape: multiple molecular dynamics analyses

Clinical studies showed that a single mutation, I117V, develops severe resistance to oseltamivir, the first orally active influenza A neuraminidase inhibitor, in highly pathogenic H5N1 influenza A viruses.

Identification of novel compounds against an R294K substitution of influenza A (H7N9) virus using ensemble based drug virtual screening

Influenza virus H7N9 foremost emerged in China in 2013 and killed hundreds of people in Asia since they possessed all mutations that enable them to resist to all existing influenza drugs, resulting in high mortality to human. In the effort to identify novel inhibitors combat resistant strains of influenza virus H7N9; we performed virtual screening targeting the Neuraminidase (NA) protein against natural compounds of traditional Chinese medicine database (TCM) and ZINC natural products. Compounds expressed high binding affinity to the target protein was then evaluated for molecular properties to determine drug-like molecules. 4 compounds showed their binding energy less than -11 Kcal/mol were selected for molecular dynamics (MD) simulation to capture intermolecular interactions of ligand-protein complexes. The molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method was utilized to estimate binding free energy of the complex. In term of stability, NA-7181 (IUPAC namely {9-Hydroxy-10-[3-(trifluoromrthyl) cyclohexyl]-4.8-diazatricyclo [6.4.0.02,6]dodec-4-yl}(perhydro-1H-inden-5-yl)formaldehyde) achieved stable conformation after 20 ns and 27 ns for ligand and protein root mean square deviation, respectively. In term of binding free energy, 7181 gave the negative value of -30.031 (KJ/mol) indicating the compound obtained a favourable state in the active site of the protein.

Highly efficient iridium-catalyzed asymmetric hydrogenation of β-acylamino nitroolefins

A highly Ir-catalyzed enantioselective hydrogenation of β-acylamino nitroolefins is first reported, which provides straightforward access to chiral β-amino nitroalkanes in excellent enantioselectivities (up to >99.9% ee).