Effect of influenza H1N1 neuraminidase V116A and I117V mutations on NA activity and sensitivity to NA inhibitors

High pathogenic avian influenza A(H5) viruses of clade 2.3.4.4b in Europe-Why trends of virus evolution are more difficult to predict

Since 2016, A(H5Nx) high pathogenic avian influenza (HPAI) virus of clade 2.3.4.4b has become one of the most serious global threats not only to wild and domestic birds, but also to public health. In recent years, important changes in the ecology, epidemiology, and evolution of this virus have been reported, with an unprecedented global diffusion and variety of affected birds and mammalian species. After the two consecutive and devastating epidemic waves in Europe in 2020-2021 and 2021-2022, with the second one recognized as one of the largest epidemics recorded so far, this clade has begun to circulate endemically in European wild bird populations. This study used the complete genomes of 1,956 European HPAI A(H5Nx) viruses to investigate the virus evolution during this varying epidemiological outline. We investigated the spatiotemporal patterns of A(H5Nx) virus diffusion to/from and within Europe during the 2020-2021 and 2021-2022 epidemic waves, providing evidence of ongoing changes in transmission dynamics and disease epidemiology. We demonstrated the high genetic diversity of the circulating viruses, which have undergone frequent reassortment events, providing for the first time a complete overview and a proposed nomenclature of the multiple genotypes circulating in Europe in 2020-2022. We described the emergence of a new genotype with gull adapted genes, which offered the virus the opportunity to occupy new ecological niches, driving the disease endemicity in the European wild bird population. The high propensity of the virus for reassortment, its jumps to a progressively wider number of host species, including mammals, and the rapid acquisition of adaptive mutations make the trend of virus evolution and spread difficult to predict in this unfailing evolving scenario.

Computational modelling studies of some 1,3-thiazine derivatives as anti-influenza inhibitors targeting H1N1 neuraminidase via 2D-QSAR, 3D-QSAR, molecular docking, and ADMET predictions

Background

Influenza virus disease remains one of the most contagious diseases that aided the deaths of many patients, especially in this COVID-19 pandemic era. Recent discoveries have shown that the high prevalence of influenza and SARS-CoV-2 coinfection can rapidly increase the death rate of patients. Hence, it became necessary to search for more potent inhibitors for influenza disease therapy. The present study utilized some computational modeling concepts such as 2D-QSAR, 3D-QSAR, molecular docking simulation, and ADMET predictions of some 1,3-thiazine derivatives as inhibitors of influenza neuraminidase (NA).

Results

The 2D-QSAR modeling results showed GFA-MLR ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9192, Q2 = 0.8767, R2adj = 0.8991, RMSE = 0.0959, $$R_{{{\text{test}}}}^{2}$$ R test 2  = 0.8943, $$R_{{{\text{pred}}}}^{2}$$ R pred 2  = 0.7745) and GFA-ANN ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9227, Q2 = 0.9212, RMSE = 0.0940, $$R_{{{\text{test}}}}^{2}$$ R test 2  = 0.8831, $$R_{{{\text{pred}}}}^{2}$$ R pred 2  = 0.7763) models with the computed descriptors as ATS7s, SpMax5_Bhv, nHBint6, and TDB9m for predicting the NA inhibitory activities of compounds which have passed the global criteria of accepting QSAR model. The 3D-QSAR modeling was carried out based on the comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA). The CoMFA_ES ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9620, Q2 = 0.643) and CoMSIA_SED ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.8770, Q2 = 0.702) models were found to also have good and reliable predicting ability. The compounds were also virtually screened based on their binding scores via molecular docking simulations with the active site of the NA (H1N1) target receptor which also confirms their resilient potency. Four potential lead compounds (4, 7, 14, and 15) with the relatively high inhibitory rate (> 50%) and docking (> − 6.3 kcal/mol) scores were identified as the possible lead candidates for in silico exploration of improved anti-influenza agents.

Conclusion

The drug-likeness and ADMET predictions of the lead compounds revealed non-violation of Lipinski’s rule and good pharmacokinetic profiles as important guidelines for rational drug design. Hence, the outcome of this research set a course for the in silico design and exploration of novel NA inhibitors with improved potency.

In-silico modelling studies of 5-benzyl-4-thiazolinone derivatives as influenza neuraminidase inhibitors via 2D-QSAR, 3D-QSAR, molecular docking, and ADMET predictions

Influenza virus disease is one of the most infectious diseases responsible for many human deaths, and the high mutability of the virus causes drug resistance effects in recent times. As such, it became necessary to explore more inhibitors that could avert future influenza pandemics. The present research utilized some in-silico modelling concepts such as 2D-QSAR, 3D-QSAR, molecular docking simulation, and ADMET predictions on some 5-benzyl-4-thiazolinone derivatives as influenza neuraminidase (NA) inhibitors. The 2D-QSAR modelling results revealed GFA-MLR (R train 2 =0.8414, Q2 = 0.7680) and GFA-ANN (R train ​ 2 =0.8754, Q2 = 0.8753) models with the most relevant descriptors (MATS3i, SpMax5_Bhe, minsOH and VE3_D) for predicting the inhibitory activities of the molecules which has passed the global criteria of accepting QSAR models. The results of the 3D-QSAR modelling results showed that CoMFA_ES (R train ​ 2 =0.9030, Q2 = 0.5390) and CoMSIA_EA (R train 2 =0.880, Q2 = 0.547) models are having good predicting ability among other developed models. The molecules were virtually screened via molecular docking simulation with the active site of NA protein receptor (pH1N1) which confirms their resilient potency when compared with zanamivir standard drug. Molecule 11 as the most potent molecule formed more H-bond interactions with the key residues such as TRP178, ARG152, ARG292, ARG371, and TYR406 that triggered the catalytic reactions for NA inhibition. Furthermore, six (6) molecules (9, 10, 11, 17, 22, and 31) with relatively high inhibitory activities and docking scores were identified as the possible leads for in-silico exploration of novel NA inhibitors. The drug-likeness and ADMET predictions of the lead molecules revealed non-violation of Lipinski's rule and good pharmacokinetic profiles respectively, which are important guidelines for rational drug design. Hence, the outcome of this study overlaid a solid foundation for the in-silico design and exploration of novel NA inhibitors with improved potency.

Influenza antivirals and animal models

Influenza A and B viruses are among the most prominent human respiratory pathogens. About 3-5 million severe cases of influenza are associated with 300 000-650 000 deaths per year globally. Antivirals effective at reducing morbidity and mortality are part of the first line of defense against influenza. FDA-approved antiviral drugs currently include adamantanes (rimantadine and amantadine), neuraminidase inhibitors (NAI; peramivir, zanamivir, and oseltamivir), and the PA endonuclease inhibitor (baloxavir). Mutations associated with antiviral resistance are common and highlight the need for further improvement and development of novel anti-influenza drugs. A summary is provided for the current knowledge of the approved influenza antivirals and antivirals strategies under evaluation in clinical trials. Preclinical evaluations of novel compounds effective against influenza in different animal models are also discussed.

Physics of self-rolling viruses

Viruses are right at the interface of inanimate matter and life. However, recent experiments [Sakai et al., J. Virol. 92, e01522-17 (2018)0022-538X10.1128/JVI.01522-17] have shown that some influenza strains can actively roll on glycan-covered surfaces. In a previous letter [Ziebert and Kulić, Phys. Rev. Lett. 126, 218101 (2021)0031-900710.1103/PhysRevLett.126.218101] we suggested this to be a form of viral surface metabolism: a collection of spike proteins that attach to and cut the glycans act as a self-organized mechano-chemical motor. Here we study in more depth the physics of the emergent self-rolling states. We give scaling arguments how the motion arises, substantiated by a detailed analytical theory that yields the full torque-angular velocity relation of the self-organized motor. Stochastic Gillespie simulations are used to validate the theory and to quantify stochastic effects like virus detachment and reversals of its direction. Finally, we also cross-check several approximations made previously and show that the proposed mechanism is very robust. All these results point together to the statistical inevitability of viral rolling in the presence of enzymatic activity.

Theoretical insights into the molecular mechanism of I117V mutation in neuraminidase mediated reduction of oseltamivir drug susceptibility in A/H5N1 influenza virus

The substitution of Ile to Val at residue 117 (I117V) of neuraminidase (NA) reduces the susceptibility of the A/H5N1 influenza virus to oseltamivir (OTV). However, the molecular mechanism by which the I117V mutation affects the intermolecular interactions between NA and OTV has not been fully elucidated. In this study, we performed molecular dynamics (MD) simulations to analyze the characteristic conformational changes that contribute to the reduced binding affinity of NA to OTV after the I117V mutation. The results of MD simulations revealed that after the I117V mutation in NA, the changes in the secondary structure around the mutation site had a noticeable effect on the residue interactions in the OTV-binding site. In the case of the WT NA-OTV complex, the positively charged side chain of R118, located in the β-sheet region, frequently interacted with the negatively charged side chain of E119, which is an amino acid residue in the OTV-binding site. This can reduce the electrostatic repulsion of E119 toward D151, which is also a negatively charged residue in the OTV-binding site, so that both E119 and D151 simultaneously form hydrogen bonds with OTV more frequently, which greatly contributes to the binding affinity of NA to OTV. After the I117V mutation in NA, the side chain of R118 interacted with the side chain of E119 less frequently, likely because of the decreased tendency of R118 to form a β-sheet structure. As a result, the electrostatic repulsion of E119 toward D151 is greater than that of the WT case, making it difficult for both E119 and D151 to simultaneously form hydrogen bonds with OTV, which in turn reduces the binding affinity of NA to OTV. Hence, after the I117V mutation in NA, influenza viruses are less susceptible to OTV because of conformational changes in residues of R118, E119, and D151 around the mutation site and in the binding site.

How Influenza's Spike Motor Works

While often believed to be a passive agent that merely exploits its host's metabolism, the influenza virus has recently been shown to actively move across glycan-coated surfaces. This form of enzymatically driven surface motility is currently not well understood and has been loosely linked to burnt-bridge Brownian ratchet mechanisms. Starting from known properties of influenza's spike proteins, we develop a physical model that quantitatively describes the observed motility. It predicts a collectively emerging dynamics of spike proteins and surface-bound ligands that combined with the virus' geometry give rise to a self-organized rolling propulsion. We show that in contrast to a Brownian ratchet, the rotary spike drive is not fluctuation driven but operates optimally as a macroscopic engine in the deterministic regime. The mechanism also applies to relatives of influenza and to man-made analogs like DNA monowheels and should give guidelines for their optimization.

Laninamivir-Interferon Lambda 1 Combination Treatment Promotes Resistance by Influenza A Virus More Rapidly than Laninamivir Alone

Each year, 5% to 20% of the population of the United States becomes infected with influenza A virus. Combination therapy with two or more antiviral agents has been considered a potential treatment option for influenza virus infection. However, the clinical results derived from combination treatment with two or more antiviral drugs have been variable. We examined the effectiveness of cotreatment with two distinct classes of anti-influenza drugs, i.e., neuraminidase (NA) inhibitor, laninamivir, and interferon lambda 1 (IFN-λ1), against the emergence of drug-resistant virus variants in vitro We serially passaged pandemic A/California/04/09 [A(H1N1)pdm09] influenza virus in a human lung epithelial cell line (Calu-3) in the presence or absence of increasing concentrations of laninamivir or laninamivir plus IFN-λ1. Surprisingly, laninamivir used in combination with IFN-λ1 promoted the emergence of the E119G NA mutation five passages earlier than laninamivir alone (passage 2 versus passage 7, respectively). Acquisition of this mutation resulted in significantly reduced sensitivity to the NA inhibitors laninamivir (∼284-fold) and zanamivir (∼1,024-fold) and decreased NA enzyme catalytic activity (∼5-fold) compared to the parental virus. Moreover, the E119G NA mutation emerged together with concomitant hemagglutinin (HA) mutations (T197A and D222G), which were selected more rapidly by combination treatment with laninamivir plus IFN-λ1 (passages 2 and 3, respectively) than by laninamivir alone (passage 10). Our results show that treatment with laninamivir alone or in combination with IFN-λ1 can lead to the emergence of drug-resistant influenza virus variants. The addition of IFN-λ1 in combination with laninamivir may promote acquisition of drug resistance more rapidly than treatment with laninamivir alone.

Impact of long-term parental exposure to Tamiflu metabolites on the development medaka offspring (Oryzias latipes)

Despite the widespread use of the antiviral drug, Tamiflu®, little is known about the long-term toxic effects of drug or its metabolites in an aquatic ecosystem. This study integrated epidemiological and ecotoxicological methods to determine environmentally relevant concentrations of Tamiflu. A model based on the species medaka (Oryzias latipes) was then used to determine the health status and reproductivity of adults exposed to the drug as well as the embryonic development of offspring. The proposed ecotoxicological model was also used to quantitatively and qualitatively evaluate the toxicodynamic parameters related to egg production, hatchability, and development. Our results revealed that at an environmentally relevant exposure, Tamiflu and its metabolites had no adverse effects on growth, survival, or fecundity of adult medaka. Nonetheless, we observed a reduction in hatchability under exposure to 300 μg L^-1 and a reduction in body length under exposure exceeding 90 μg L^-1. Under exposure to 300 μg L^-1, the estimated spawning time to reach 50% of the maximum percentage of cumulative egg production (ET50) far exceeded that of the control group (without exposure to Tamiflu). We also observed a ∼ 3-fold decrease in maximum egg hatching (E_max). Based on an integrated epidemiological and ecotoxicological model, predictions of environmental concentrations of Tamiflu and its metabolites revealed that the influenza subtypes associated with increases in environmental concentrations: A(H3N2) > A(H1N1) > type B (in order of their effects). We also determined that A(H3N2) posed a potential risk to hatchability and development. Note however, the environmental concentrations of Tamiflu and its metabolites in most countries are lower than the effect concentrations derived in this study, indicating no hazards for aquatic environments. We recommend the use of hatchability and embryonic development as indicators in assessing the effects of long-term parental exposure to Tamiflu metabolites.

Influenza A virus hemagglutinin mutations associated with use of neuraminidase inhibitors correlate with decreased inhibition by anti-influenza antibodies

Background

Vaccination and the use of neuraminidase inhibitors (NAIs) are currently the front lines of defense against seasonal influenza. The activity of influenza vaccines and antivirals drugs such as the NAIs can be affected by mutations in the influenza hemagglutinin (HA) protein. Numerous HA substitutions have been identified in nonclinical NAI resistance-selection experiments as well as in clinical specimens from NAI treatment or surveillance studies. These mutations are listed in the prescribing information (package inserts) for FDA-approved NAIs, including oseltamivir, zanamivir, and peramivir.

Methods

NAI treatment-emergent H1 HA mutations were mapped onto the H1N1 HA1 trimeric crystal structure and most of them localized to the HA antigenic sites predicted to be important for anti-influenza immunity. Recombinant A/California/04/09 (H1N1)-like viruses carrying HA V152I, G155E, S162 N, S183P, and D222G mutations were generated. We then evaluated the impact of these mutations on the immune reactivity and replication potential of the recombinant viruses in a human respiratory epithelial cell line, Calu− 3.

Results

We found that the G155E and D222G mutations significantly increased viral titers ~ 13-fold compared to the wild-type virus. The hemagglutination and microneutralization activity of goat and ferret antisera, monoclonal antibodies, and human serum samples raised against pandemic A(H1N1)pdm09 viruses was ~ 100-fold lower against mutants carrying G155E or D222G compared to the wild-type virus.

Conclusions

Although the mechanism by which HA mutations emerge during NAI treatment is uncertain, some NAI treatment-emergent HA mutations correlate with decreased immunity to influenza virus.