Three new powerful oseltamivir derivatives for inhibiting the neuraminidase of influenza virus

Owing to its unique function in assisting the release of newly formed virus particles from the surface of an infected cell, neuraminidase, an antigenic glycoprotein enzyme, is a main target for drug design against influenza viruses. The group-1 neuraminidase of influenza virus possesses a 150-cavity, which is adjacent to the active pocket, and which renders conformational change from the ‘open’ form to the ‘closed’ form when the enzyme is binding with a ligand. Using AutoGrow evolutionary algorithm, one very unique fragment is screened out from the fragment databases by exploiting additional interactions with the 150-cavity. Subsequently, three derivatives were constructed by linking the unique fragment to oseltamivir at its three different sites. The three derivatives thus formed show much stronger inhibition power than oseltamivir, and hence may become excellent candidates for developing new and more powerful drugs for treating influenza. Or at the very least, the findings may stimulate new strategy or provide useful insights for working on the target vitally important to the health of human beings.     

Analogs of zanamivir with modified C4-substituents as the inhibitors against the group-1 neuraminidases of influenza viruses

Unlike the group-2 neuraminidase, the group-1 neuraminidase of influenza virus possesses a flexible loop (the 150-loop) and a cavity (the 150-cavity) adjacent to the active site, and renders a conformational change from the ‘open’ form to the ‘closed’ form on binding with substrate (sialo-glycoprotein) or inhibitor (e.g., zanamivir). Zanamivir derivative 8a having an extended (piperazinocarbonyl)propyl substituent at the internal N-position of the guanidino group is designed as a possible inhibitor on the basis of computer docking to the open form of N1 subtype neuraminidase. Indeed, compound 8a exhibits strong neuraminidase inhibition and good anti-influenza activity against H1N1 virus with IC₅₀ = 2.15 μM and EC₅₀ = 0.77 μM, respectively. This study may provide a clue to future design of better group-1 neuraminidase inhibitors.     

Virtual screening pipeline and ligand modelling for H5N1 neuraminidase

The H5N1 virus neuraminidase structure was solved in two different conformations depending on the inhibitor concentration. In the absence of oseltamivir or at a low concentration, the neuraminidase structure assumes an open form that closes at a high oseltamivir concentration due to the shift of the so-called 150-loop near the active site. Although the close conformation is similar to all the other structurally known neuraminidase types, it doesn’t appear to be the most likely physiological condition for N1.To investigate the specific ligand binding properties of the open form, we screened by docking simulation, a large dataset of ligands and compared the results with closed form. The virtual screening procedure was implemented in a docking pipeline that also performs a step-by-step, target specific, filtering approach for data reduction. The selected ligands display binding ability involving multiple sites of interaction including the active site and an adjacent cavity made available by the 150-loop shift. Two ligands are especially interesting and are proposed as substituents to design oseltamivir derivatives specifically suited for the open conformation.     

Binding mechanism of oseltamivir and influenza neuraminidase suggests perspectives for the design of new anti-influenza drugs

Oseltamivir is a widely used influenza virus neuraminidase (NA) inhibitor that prevents the release of new virus particles from host cells. However, oseltamivir-resistant strains have emerged, but effective drugs against them have not yet been developed. Elucidating the binding mechanisms between NA and oseltamivir may provide valuable information for the design of new drugs against NA mutants resistant to oseltamivir. Here, we conducted large-scale (353.4 μs) free-binding molecular dynamics simulations, together with a Markov State Model and an importance-sampling algorithm, to reveal the binding process of oseltamivir and NA. Ten metastable states and five major binding pathways were identified that validated and complemented previously discovered binding pathways, including the hypothesis that oseltamivir can be transferred from the secondary sialic acid binding site to the catalytic site. The discovery of multiple new metastable states, especially the stable bound state containing a water-mediated hydrogen bond between Arg118 and oseltamivir, may provide new insights into the improvement of NA inhibitors. We anticipated the findings presented here will facilitate the development of drugs capable of combating NA mutations.     

Locking the 150-Cavity Open: In Silico Design and Verification of Influenza Neuraminidase Inhibitors

Neuraminidase (NA) of influenza is a key target for virus infection control and the recently discovered open 150-cavity in group-1 NA provides new opportunity for novel inhibitors design. In this study, we used a combination of theoretical methods including fragment docking, molecular linking and molecular dynamics simulations to design ligands that specifically target at the 150-cavity. Through in silico screening of a fragment compound library on the open 150-cavity of NA, a few best scored fragment compounds were selected to link with Zanamivir, one NA-targeting drug. The resultant new ligands may bind both the active site and the 150-cavity of NA simultaneously. Extensive molecular dynamics simulations in explicit solvent were applied to validate the binding between NA and the designed ligands. Moreover, two control systems, a positive control using Zanamivir and a negative control using a low-affinity ligand 3-(p-tolyl) allyl-Neu5Ac2en (ETT, abbreviation reported in the PDB) found in a recent experimental work, were employed to calibrate the simulation method. During the simulations, ETT was observed to detach from NA, on the contrary, both Zanamivir and our designed ligand bind NA firmly. Our study provides a prospective way to design novel inhibitors for controlling the spread of influenza virus.     

Influenza A virus N5 neuraminidase has an extended 150-cavity

There are 9 serotypes of neuraminidase (NA) from influenza A virus (N1 to N9), which are classified into two groups based on primary sequences (groups 1 and 2). The structural hallmark of the two groups is the presence or absence of an extra 150-cavity (formed by the 150-loop) in the active site. Thus far, structures of NAs from 6 out of the 9 serotypes have been solved. Here, we solved the N5 structure, the last unknown structure group 1 serotype with a unique Asn147 residue in its 150-loop, demonstrating that it has an extended 150-cavity that closes upon inhibitor binding.     

On the structure-based design of novel inhibitors of H5N1 influenza A virus neuraminidase (NA)

The structure-based design of novel H5N1 neuraminidase inhibitors is currently a research topic of vital importance owing to both a recent pandemic threat by the worldwide spread of H5N1 avian influenza and the high resistance of H5N1 virus to the most widely used commercial drug, oseltamivir-OTV (Tamiflu). A specific criterion used in this work for determining fully acceptable conformations of potential inhibitors is a previous experimental proposal of exploiting potential benefits for drug design offered by the ‘150-cavity’ adjacent to the NA active site. Using the crystal structure of H5N1 NA (PDB ID: 2hty) as the starting point, in a set of 54 inhibitors previously proposed by modifying the side chains of oseltamivir, 4 inhibitors were identified using two different computational strategies (ArgusLab4.0.1, FlexX-E3.0.1) both to lower the binding free energy (BFE) of oseltamivir and to have partially acceptable conformations. These 4 oseltamivr structure-based analogues were found to adopt the most promising conformations by identifying the guanidinium side chain of Arg156 as a prospective partner for making polar contacts, but none of the modified 4-amino groups of oseltamivir in the 4 favorable conformations was found to make polar contacts with the guanidinium side chain of Arg156. Hence, the structures of two additional inhibitors were designed and shown to further lower the binding free energy of OTV relative to the previous 54 inhibitors. These two novel structures clearly suggest that it may be possible for a new substituent to be developed by functional modifications at position of the 4-amino group of oseltamivir in order to make polar contacts with the guanidinium side chain of Arg156, and thereby enhance the binding of a more potent inhibitor. Several standpoints of vital importance for designing novel structures of potentially more effective H5N1 NA inhibitors are established.     

Analogue inhibitors by modifying oseltamivir based on the crystal neuraminidase structure for treating drug-resistant H5N1 virus

The worldwide spread of H5N1 avian influenza and the increasing reports about its resistance to the existing drugs have made a priority for the development of the new anti-influenza molecules. The crystal structure of H5N1 avian influenza neuraminidase reported recently by Russell et al. [R.J. Russell, L.F. Haire, D.J. Stevens, P.J. Collins, Y. P. Lin, G.M. Blackburn, A.J. Hay, S.J. Gamblin, J.J. Skehel, The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design, Nature 443 (2006) 45-49] have provided new opportunities for drug design in this regard. It is revealed through the structure that the active sites of the group-1 neuraminidases, which contain the N1 subtype, have a very different three-dimensional structure from those of group-2 neuraminidases. The key difference is in the 150-loop cavity adjacent to the conserved active site in neuraminidase. Based on these findings and by modifying oseltamivir, six analog inhibitors were proposed as candidates for developing inhibitors against H5N1 virus, particularly against the oseltamivir-resistant H5N1 virus strain.     

Design,synthesis, and evaluation of carboxyl-modified oseltamivir derivatives with improved lipophilicity as neuraminidase inhibitors

In this study, a series of carboxyl-modified oseltamivir analogs with improved lipophilicity were designed and synthesized, and their inhibitory activities against neuraminidase from influenza A virus H5N1 subtype were evaluated. The results demonstrated that compound 5m exhibited potent inhibitory activity (IC₅₀?=?1.30?±?0.23?μM), and it targeted the recently discovered 430-cavity. Compound 5m (Log D?=??0.12) is more lipophilic than oseltamivir carboxylate (Log D?=??1.69) at pH 7.4, which is potentially propitious to improved membrane permeability and oral drug absorption. Meanwhile, 5m showed high stability in human liver microsomes. The findings of this study can be valuable in identifying neuraminidase inhibitors with optimal lipophilicity and in the exploration of 430-cavity.     

Structural and functional analysis of laninamivir and its octanoate prodrug reveals group specific mechanisms for influenza NA inhibition

The 2009 H1N1 influenza pandemic (pH1N1) led to record sales of neuraminidase (NA) inhibitors, which has contributed significantly to the recent increase in oseltamivir-resistant viruses. Therefore, development and careful evaluation of novel NA inhibitors is of great interest. Recently, a highly potent NA inhibitor, laninamivir, has been approved for use in Japan. Laninamivir is effective using a single inhaled dose via its octanoate prodrug (CS-8958) and has been demonstrated to be effective against oseltamivir-resistant NA in vitro. However, effectiveness of laninamivir octanoate prodrug against oseltamivir-resistant influenza infection in adults has not been demonstrated. NA is classified into 2 groups based upon phylogenetic analysis and it is becoming clear that each group has some distinct structural features. Recently, we found that pH1N1 N1 NA (p09N1) is an atypical group 1 NA with some group 2-like features in its active site (lack of a 150-cavity). Furthermore, it has been reported that certain oseltamivir-resistant substitutions in the NA active site are group 1 specific. In order to comprehensively evaluate the effectiveness of laninamivir, we utilized recombinant N5 (typical group 1), p09N1 (atypical group 1) and N2 from the 1957 pandemic H2N2 (p57N2) (typical group 2) to carry out in vitro inhibition assays. We found that laninamivir and its octanoate prodrug display group specific preferences to different influenza NAs and provide the structural basis of their specific action based upon their novel complex crystal structures. Our results indicate that laninamivir and zanamivir are more effective against group 1 NA with a 150-cavity than group 2 NA with no 150-cavity. Furthermore, we have found that the laninamivir octanoate prodrug has a unique binding mode in p09N1 that is different from that of group 2 p57N2, but with some similarities to NA-oseltamivir binding, which provides additional insight into group specific differences of oseltamivir binding and resistance.     

Computational study on new natural polycyclic compounds of H1N1 influenza virus neuraminidase

A new strain of influenza A (H1N1) virus is a major cause of morbidity and mortality around the world. The neuraminidase of the influenza virus has been the most potential target for the anti-influenza drugs such as oseltamivir and zanamivir. However, the emergence of drug-resistant variants of these drugs makes a pressing need for the development of new neuraminidase inhibitors for controlling illness and transmission. Here a 3D structure model of H1N1 avian influenza virus neuraminidase type 1 (N1) was constructed based on the structure of the template H5N1 avian influenza virus N1. Upon application of virtual screening technique for N1 inhibitors, two novel compounds (ZINC database ID: ZINC02128091, ZINC02098378) were found as the most favorable interaction energy with N1. Docking results showed that the compounds bound not only in the active pocket, but also in a new hydrophobic cave which contains Arg368, Trp399, Ile427, Pro431 and Lys432 of N1. Our result suggested that both of the screened compounds containing the hydrophobic group bring a strong conjugation effect with Arg293, Arg368 Lys432 of N1 by pi-pi interaction. However, the control inhibitors zanamivir and oseltamivir do not have this effect. The details of N1-compound binding structure obtained will be valuable for the development of a new anti-influenza virus agent.     

基于分子对接技术筛选中草药中潜在的H1N1病毒神经氨酸酶抑制剂

甲型H1N1流感病毒是一种高度接触性急性呼吸道传染病,给人类健康造成了极大威胁。作为一个已经被证明的抗流感药物靶点,神经氨酸酶在流感病毒复制和传播中发挥重要作用,并且其活性中心的氨基酸组成高度保守。本研究中,我们使用分子对接技术筛选从中草药小分子数据库中筛选神经氨酸酶潜在的抑制剂,结果得到4个亲和力高于奥司他韦(Osel-tamivir)的化合物,并确定了他们的中草药来源。为从草药中提取,设计以及实验合成新的神经氨酸酶抑制剂提供了一定的依据和指导。     

Search for novel neuraminidase inhibitors: Design, synthesis and interaction of oseltamivir derivatives with model membrane using docking, NMR and DSC methods

As a part of our ongoing program of developing novel influenza virus inhibitors, some new derivatives of oseltamivir were prepared by modifying the amino group with glycyl, acetyl, benzyl and prolyl moieties. The interactions of these derivatives with neuraminidase have been probed by molecular modeling techniques. Further, the interaction of these derivatives with model membranes prepared from DPPC and the effect on the thermotropic behavior and polymorphism of the bilayers have been investigated by multinuclear NMR and DSC methods. Results indicate that the glycyl derivative of oseltamivir has the most profound effects on the membrane, compared to other derivatives and seems to be the most promising derivative for further pharmacological evaluation as a neuraminidase inhibitor.     

Insights from investigating the interaction of oseltamivir (Tamiflu) with neuraminidase of the 2009 H1N1 swine flu virus

The neuraminidase (NA) of influenza virus is the target of anti-flu drugs oseltamivir and zanamivir. Clinical practices showed that oseltamivir was effective to treat the 2009-H1N1 influenza but failed to the 2006-H5N1 avian influenza. To perform an in-depth analysis on such a drug-resistance problem, the 2009-H1N1-NA structure was developed. To compare it with the crystal 2006-H5N1-NA structure as well as the 1918 influenza virus H1N1-NA structure, the multiple sequential and structural alignments were performed. It has been revealed that the hydrophobic residue Try347 in H5N1-NA does not match with the hydrophilic carboxyl group of oseltamivir as in the case of H1N1-NA. This may be the reason why H5N1 avian influenza virus is drug-resistant to oseltamivir. The finding provides useful insights for how to modify the existing drugs, such as oseltamivir and zanamivir, making them not only become more effective against H1N1 virus but also effective against H5N1 virus.     

The Mechanism by which 146-N-Glycan Affects the Active Site of Neuraminidase

One of the most conserved glycosylation sites of neuraminidase (NA) is 146-N-glycan. This site is adjacent to the 150-cavity of NA, which is found within the active site and thought to be a target for rational drug development against the antiviral resistance of influenza. Here, through a total of 2.4 μs molecular dynamics (MD) simulations, we demonstrated that 146-N-glycan can stabilize the conformation of the 150-loop that controls the volume of the 150-cavity. Moreover, with 146-N-glycan, our simulation result was more consistent with crystal structures of NAs than simulations conducted without glycans. Cluster analysis of the MD trajectories showed that 146-N-glycan adopted three distinct conformations: monomer-bridged, dimer-bridged and standing. Of these conformations, the dimer-bridged 146-N-glycan was the most stable one and contributed to stabilization of the 150-loop conformation. Furthermore, our simulation revealed that various standing conformations of 146-N-glycan could block the entrance of the binding pocket. This result was consistent with experimental data and explained the relatively low activity of inhibitors with flexible substituents toward the 150-cavity. Together, our results lead us to hypothesize that rigid and hydrophobic substituents could serve as better inhibitors targeting the 150-cavity.     

Inhibitory potency of flavonoid derivatives on influenza virus neuraminidase

The constant risk of emerging new influenza virus strains that are resistant to established inhibitors like oseltamivir leaves influenza neuraminidase (NA) a prominent target for drug design. The inhibitory activity of several flavonoid derivatives was experimentally tested in comparison to oseltamivir for the NA expressed by the seasonal influenza virus strains A/California/7/09 (A(H1N1)pdm09), A/Perth/16/09 (A(H3N2)), and B/Brisbane/60/08. IC₅₀ values of polyphenols confirmed moderate inhibition in the μM range. Structurally, the amount and site of glycosylation of tested flavonoids have no significant influence on their inhibitory potency. In a pharmacophore-based docking approach the structure–activity relationship was evaluated. Molecular dynamics simulations revealed highly flexible parts of the enzyme and the contribution of salt bridges to the structural stability of NA. The findings of this study elucidate the impact of flavonoids on viral neuraminidase activity and the analysis of their modes of action provide valuable information about the mechanism of NA inhibition.     

Inhibitor selectivity of a new class of oseltamivir analogs against viral neuraminidase over human neuraminidase enzymes

The viral neuraminidase enzyme is an established target for anti-influenza pharmaceuticals. However, viral neuraminidase inhibitors could have off-target effects due to interactions with native human neuraminidase enzymes. We report the activity of a series of known inhibitors of the influenza group-1 neuraminidase enzyme (N1 subtype) against recombinant forms of the human neuraminidase enzymes NEU3 and NEU4. These inhibitors were designed to take advantage of an additional enzyme pocket (known as the 150-cavity) near the catalytic site of certain viral neuraminidase subtypes (N1, N4 and N8). We find that these modified derivatives have minimal activity against the human enzymes, NEU3 and NEU4. Two compounds show moderate activity against NEU3, possibly due to alternative binding modes available to these structures. Our results reinforce that recognition of the glycerol side-chain is distinct between the viral and human NEU enzymes, and provide experimental support for improving the selectivity of viral neuraminidase inhibitors by exploiting the 150-cavity found in certain subtypes of viral neuraminidases.     

A nonsynonymous SNP in human cytosolic sialidase in a small Asian population results in reduced enzyme activity： potential link with severe adverse reactions to oseltamivir

The use of oseltamivir, widely stockpiled as one of the drugs for use in a possible avian influenza pandemic, has been reported to be associated with neuropsychiatric disorders and severe skin reactions, primarily in Japan. Here we identified a nonsynonymous SNP （single nucleotide polymorphism） in dbSNP database, R41Q, near the enzymatic active site of human cytosolic sialidase, a homologue of virus neuraminidase that is the target of oseltamivir. This SNP occurred in 9.29% of Asian population and none of European and African American population. Our structural analyses and Ki measurements using in vitro sialidase assays indicated that this SNP could increase the unintended binding affinity of human sialidase to oseltamivir carboxylate, the active form of oseltamivir, thus reducing sialidase activity. In addition, this SNP itself results in an enzyme with an intrinsically lower sialidase activity, as shown by its increased Km and decreased Vmax values. Theoretically administration of oseltamivir to people with this SNP might further reduce their sialidase activity. We note the similarity between the reported neuropsychiatric side effects ofoseltamivir and the known symptoms of human sialidase-related disorders. We propose that this Asian-enriched sialidase variation caused by the SNP, likely in homozygous form, may be associated with certain severe adverse reactions to oseltamivir.     

Susceptibility of antiviral drugs against 2009 influenza A (H1N1) virus

The recent outbreak of the novel strain of influenza A (H1N1) virus has raised a global concern of the future risk of a pandemic. To understand at the molecular level how this new H1N1 virus can be inhibited by the current anti-influenza drugs and which of these drugs it is likely to already be resistant to, homology modeling and MD simulations have been applied on the H1N1 neuraminidase complexed with oseltamivir, and the M2-channel with adamantanes bound. The H1N1 virus was predicted to be susceptible to oseltamivir, with all important interactions with the binding residues being well conserved. In contrast, adamantanes are not predicted to be able to inhibit the M2 function and have completely lost their binding with the M2 residues. This is mainly due to the fact that the M2 transmembrane of the new H1N1 strain contains the S31N mutation which is known to confer resistance to adamantanes.     

In silico analysis of drug-resistant mutant of neuraminidase (N294S) against oseltamivir

The recent H1N1 influenza pandemic has attracted worldwide attention due to the high infection rate. Oseltamivir is a new class of anti-viral agent approved for the treatment and prevention of influenza infections. The principal target for this drug is a virus surface glycoprotein, neuraminidase (NA), which facilitates the release of nascent virus and thus spreads infection. Until recently, only a low prevalence of neuraminidase inhibitor (NAI) resistance (<1 %) had been detected in circulating viruses. However, there have been reports of significant numbers of A (H1N1) influenza strains with a N294S neuraminidase mutation that was highly resistant to the NAI, oseltamivir. Hence, in the present study, we highlight the effect of point mutation-induced oseltamivir resistance in H1N1 subtype neuraminidases by molecular simulation approach. The docking analysis reveals that mutation (N294S) significantly affects the binding affinity of oseltamivir with mutant type NA. This is mainly due to the decrease in the flexibility of binding site residues and the difference in prevalence of hydrogen bonds in the wild and mutant structures. This study throws light on the possible effects of drug-resistant mutations on the large functionally important collective motions in biological systems.