Induced opening of influenza virus neuraminidase N2 150-loop suggests an important role in inhibitor binding

Identification of Novel Neuraminidase Inhibitors as Potential Anti-Influenza Agents: Virtual Screening, Molecular Docking, in vitro Validation and Molecular Dynamic Simulation Studies

The influenza virus causes approximately hundreds of thousands of deaths annually. Coupled with the emergence of drug resistance, there is an urgent need to develop new drugs for the treatment of influenza. Neuraminidase (NA) has long been recognized as a valid drug target for anti-influenza therapy. Herein, in order to identify potential NA inhibitors with novel structures, we employed a structure-based virtual screening strategy to screen a library containing 1.6 million compounds. Based on XP docking score and free energy calculation results, the three compounds E570-1769, K788-4718, and C071-0424 were selected that may have better binding affinity for the NA protein compared to oseltamivir. Amongst, E570-1769 was identified to be the most potential hit. Docking study showed that E570-1769 bound to NA with a binding energy of -10.3 kcal/mol. Moreover, in silico ADME/T studies demonstrated the druggability of E570-1769 was quite well. Furthermore, in vitro assay demonstrated that E570-1769 inhibited the wild-type and H274Y-muatated NAs with IC50 values of 72.6 μM and 229 μM, respectively. Additionally, molecular dynamic (MD) simulation studies were performed to gain a deep insight into the binding modes of E570-1769 in complex with NA. While less potent than oseltamivir, the novel structure of E570-1769 and promising ADME/T properties indicates it as a promising lead for future research.

Discovery of new quinoline derivatives bearing 1-aryl-1,2,3-triazole motif as influenza H1N1 virus neuraminidase inhibitors

Sporadically and periodically, influenza outbreaks threaten global health and the economy. Antigen drift-induced influenza virus mutations hamper antiviral drug development. Thus, a novel antiviral agent is urgently needed to address medication inefficacy issues. Herein, sixteen new quinoline-triazole hybrids 6a-h and 9a-h were prepared and evaluated in vitro against the H1N1 virus. In particular, 6d, 6e, and 9b showed promising H1N1 antiviral activity with selective index (SI) CC50/IC50 values of 15.8, 37, and 29.15. After that, the inhibition rates for various mechanisms of action (virus replication, adsorption, and virucidal activity) were investigated for the most efficient candidates 6d, 6e, and 9b. Additionally, their ability to inhibit neuraminidase was evaluated. With an IC50 value of 0.30 µM, hybrid 6d demonstrated effective and comparable inhibitory activity to Oseltamivir. Ultimately, molecular modeling investigations, encompassing molecular docking and molecular dynamic simulations, were conducted to provide a scientific basis for the observed antiviral results.

Discovery andsynthesis of novel benzoylhydrazone neuraminidase inhibitors

Neuraminidase (NA) serves as a promising target for the exploration and development of anti-influenza drugs. In this work, lead compound 5 was discovered through pharmacophore-based virtual screening and molecular dynamics simulation, and 14 new compounds were obtained by modifying the lead compound 5 based on pharmacophore features. The biological activity test shows that 5n (IC50 = 0.13 μM) has a better inhibitory effect on wild-type NA (H5N1), while 5i (IC50 = 0.44 μM) has a prominent inhibitory effect on mutant NA (H5N1-H274Y), both of them are better than the positive control oseltamivir carboxylate (OSC). The analysis of docking results indicate that the good activities of compounds 5n and 5i may be attributed to the thiophene ring in 5n can stretch into the 150-cavity of NA, whereas the thiophene moiety in 5i can extend to the 430-cavity of NA. The findings of this study may be helpful for the discovery of new NA inhibitors.

A novel N-heterocycles substituted oseltamivir derivatives as potent inhibitors of influenza virus neuraminidase: discovery, synthesis and biological evaluation

Our previous studies have shown that the introduction of structurally diverse benzyl side chains at the C5-NH2 position of oseltamivir to occupy 150-cavity contributes to the binding affinity with neuraminidase and anti-influenza activity. To obtain broad-spectrum neuraminidase inhibitors, we designed and synthesised a series of novel oseltamivir derivatives bearing different N-heterocycles substituents that have been proved to induce opening of the 150-loop of group-2 neuraminidases. Among them, compound 6k bearing 4-((r)-2-methylpyrrolidin-1-yl) benzyl group exhibited antiviral activities similar to or weaker than those of oseltamivir carboxylate against H1N1, H3N2, H5N1, H5N6 and H5N1-H274Y mutant neuraminidases. More encouragingly, 6k displayed nearly 3-fold activity enhancement against H3N2 virus over oseltamivir carboxylate and 2-fold activity enhancement over zanamivir. Molecular docking studies provided insights into the explanation of its broad-spectrum potency against wild-type neuraminidases. Overall, as a promising lead compound, 6k deserves further optimisation by fully considering the ligand induced flexibility of the 150-loop.

New insights into the neuraminidase-mediated hemagglutination activity of influenza A(H3N2) viruses

Influenza virus neuraminidase (NA) can act as a receptor-binding protein, a role commonly attributed to hemagglutinin (HA). In influenza A(H3N2) viruses, three NA amino acid residues have previously been associated with NA-mediated hemagglutination: T148, D151, and more recently, H150. These residues are part of the 150-loop of the NA monomer. Substitutions at 148 and 151 arise from virus propagation in laboratory cell cultures, whereas changes at 150 occurred during virus evolution in the human host. In this study, we examined the effect of natural amino acid polymorphism at position 150 on NA-mediated hemagglutination. Using the A/Puerto Rico/8/34 backbone, we generated a comprehensive panel of recombinant A(H3N2) viruses that have different NAs but shared an HA that displays poor binding to red blood cells (RBCs). None of the tested substitutions at 150 (C, H, L, R, and S) promoted NA-binding. However, we identified two new determinants of NA-binding, Q136K and T439R, that emerged during virus culturing. Similar to T148I, both Q136K and T439R reduced NA enzyme activity by 48-86% and inhibition (14- to 173-fold) by the NA inhibitor zanamivir. NA-binding was observed when a virus preparation contained approximately 10% of NA variants with either T148I or T439R, highlighting the benefit of using deep sequencing in virus characterization. Taken together, our findings provide new insights into the molecular mechanisms underlying the ability of NA to function as a binding protein. Information gained may aid in the design of new and improved NA-targeting antivirals.

Discovery of N-substituted oseltamivir derivatives as novel neuraminidase inhibitors with improved drug resistance profiles and favorable drug-like properties

To yield potent neuraminidase inhibitors with improved drug resistance and favorable drug-like properties, two series of novel oseltamivir derivatives targeting the 150-cavity of neuraminidase were designed, synthesized, and biologically evaluated. Among the synthesized compounds, the most potent compound 43b bearing 3-floro-4-cyclopentenylphenzyl moiety exhibited weaker or slightly improved inhibitory activity against wild-type neuraminidases (NAs) of H1N1, H5N1, and H5N8 compared to oseltamivir carboxylate (OSC). Encouragingly, 43b displayed 62.70- and 5.03-fold more potent activity than OSC against mutant NAs of H5N1-H274Y and H1N1-H274Y, respectively. In cellular antiviral assays, 43b exerted equivalent or more potent activities against H1N1, H5N1, and H5N8 compared to OSC with no significant cytotoxicity up to 200 μM. Notably, 43b displayed potent antiviral efficacy in the embryonated egg model, in which achieved a protective effect against H5N1 and H5N8 similar to OSC. Molecular docking studies were implemented to reveal the binding mode of 43b in the binding pocket. Moreover, 43b possessed improved physicochemical properties and ADMET properties compared to OSC by in silico prediction. Taken together, 43b appeared to be a promising lead compound for further investigation.

Influenza virus and pneumococcal neuraminidases enhance catalysis by similar yet distinct sialic acid-binding strategies

Influenza A viruses and the bacterium Streptococcus pneumoniae (pneumococci) both express neuraminidases that catalyze release of sialic acid residues from oligosaccharides and glycoproteins. Although these respiratory pathogen neuraminidases function in a similar environment, it remains unclear if these enzymes use similar mechanisms for sialic acid cleavage. Here, we compared the enzymatic properties of neuraminidases from two influenza A subtypes (N1 and N2) and the pneumococcal strain TIGR4 (NanA, NanB, and NanC). Insect cell-produced N1 and N2 tetramers exhibited calcium-dependent activities and stabilities that varied with pH. In contrast, E. coli-produced NanA, NanB, and NanC were isolated as calcium insensitive monomers with stabilities that were more resistant to pH changes. Using a synthetic substrate (MUNANA), all neuraminidases showed similar pH optimums (pH 6-7) that were primarily defined by changes in catalytic rate rather than substrate binding affinity. Upon using a multivalent substrate (fetuin sialoglycans), much higher specific activities were observed for pneumococcal neuraminidases that contain an additional lectin domain. In virions, N1 and especially N2 also showed enhanced specific activity toward fetuin that was lost upon the addition of detergent, indicating the sialic acid-binding capacity of neighboring hemagglutinin molecules likely contributes to catalysis of natural multivalent substrates. These results demonstrate that influenza and pneumococcal neuraminidases have evolved similar yet distinct strategies to optimize their catalytic activity.

Structural and inhibitor sensitivity analysis of influenza B-like viral neuraminidases derived from Asiatic toad and spiny eel

Influenza virus neuraminidase (NA) is an important target for antiviral development because it plays a crucial role in releasing newly assembled viruses. Two unique influenza-like virus genomes were recently reported in the Wuhan Asiatic toad and Wuhan spiny eel. Their NA genes appear to be highly divergent from all known influenza NAs, raising key questions as to whether the Asiatic toad influenza-like virus NA (tNA) and spiny eel NA (eNA) have canonical NA activities and structures and whether they show sensitivity to NA inhibitors (NAIs). Here, we found that both tNA and eNA have neuraminidase activities. A detailed structural analysis revealed that tNA and eNA present similar overall structures to currently known NAs, with a conserved calcium binding site. Inhibition assays indicated that tNA is resistant to NAIs, while eNA is still sensitive to NAIs. E119 is conserved in canonical NAs. The P119E substitution in tNA can restore sensitivity to NAIs, and, in contrast, the E119P substitution in eNA decreased its sensitivity to NAIs. The structures of NA-inhibitor complexes further provide a detailed insight into NA-inhibitor interactions at the atomic level. Moreover, tNA and eNA have unique N-glycosylation sites compared with canonical NAs. Collectively, the structural features, NA activities, and sensitivities to NAIs suggest that fish- and amphibian-derived influenza-like viruses may circulate in these vertebrates. More attention should be paid to these influenza-like viruses because their NA molecules may play roles in the emergence of NAI resistance.

Discovery of Novel Boron-Containing N-Substituted Oseltamivir Derivatives as Anti-Influenza A Virus Agents for Overcoming N1-H274Y Oseltamivir-Resistant

To address drug resistance to influenza virus neuraminidase inhibitors (NAIs), a series of novel boron-containing N-substituted oseltamivir derivatives were designed and synthesized to target the 150-cavity of neuraminidase (NA). In NA inhibitory assays, it was found that most of the new compounds exhibited moderate inhibitory potency against the wild-type NAs. Among them, compound 2c bearing 4-(3-boronic acid benzyloxy)benzyl group displayed weaker or slightly improved activities against group-1 NAs (H1N1, H5N1, H5N8 and H5N1-H274Y) compared to that of oseltamivir carboxylate (OSC). Encouragingly, 2c showed 4.6 times greater activity than OSC toward H5N1-H274Y NA. Moreover, 2c exerted equivalent or more potent antiviral activities than OSC against H1N1, H5N1 and H5N8. Additionally, 2c demonstrated low cytotoxicity in vitro and no acute toxicity at the dose of 1000 mg/kg in mice. Molecular docking of 2c was employed to provide a possible explanation for the improved anti-H274Y NA activity, which may be due to the formation of key additional hydrogen bonds with surrounding amino acid residues, such as Arg152, Gln136 and Val149. Taken together, 2c appeared to be a promising lead compound for further optimization.

Iterative Optimization and Structure-Activity Relationship Studies of Oseltamivir Amino Derivatives as Potent and Selective Neuraminidase Inhibitors via Targeting 150-Cavity

With our continuous endeavors in seeking neuraminidase (NA) inhibitors, we reported herein three series of novel oseltamivir amino derivatives with the goal of exploring the druggable chemical space inside the 150-cavity of influenza virus NAs. Among them, around half of the compounds in series C were demonstrated to be better inhibitors against both wild-type and oseltamivir-resistant group-1 NAs than oseltamivir carboxylate (OSC). Notably, compounds 12d, 12e, 15e, and 15i showed more potent or equipotent antiviral activity against H1N1, H5N1, and H5N8 viruses compared to OSC in cellular assays. Furthermore, compounds 12e and 15e exhibited high metabolic stability in human liver microsomes (HLMs) and low inhibitory effect on main cytochrome P450 (CYP) enzymes, as well as low acute/subacute toxicity and certain antiviral efficacy in vivo. Also, pharmacokinetic (PK) and molecular docking studies were performed. Overall, 12e and 15e possess great potential to serve as anti-influenza candidates and are worthy of further investigation.

Emerging HxNy Influenza A Viruses

The continuous emergence and reemergence of diverse subtypes of influenza A viruses, which are known as "HxNy" and are mediated through the reassortment of viral genomes, account for seasonal epidemics, occasional pandemics, and zoonotic outbreaks. We summarize and discuss the characteristics of historic human pandemic HxNy viruses and diverse subtypes of HxNy among wild birds, mammals, and live poultry markets. In addition, we summarize the key molecular features of emerging infectious HxNy influenza viruses from the perspectives of the receptor binding of Hx, the inhibitor-binding specificities and drug-resistance features of Ny, and the matching of the gene segments. Our work enhances our understanding of the potential threats of novel reassortant influenza viruses to public health and provides recommendations for effective prevention, control, and research of this pathogen.

Influenza A Virus Neuraminidase Inhibitors

Depending on the strain, influenza A virus causes animal, zoonotic, pandemic, or seasonal influenza with varying degrees of severity. Two surface glycoprotein spikes, hemagglutinin (HA) and neuraminidase (NA), are the most important influenza A virus antigens. NA plays an important role in the propagation of influenza virus by removing terminal sialic acid from sialyl decoy receptors and thereby facilitating the release of viruses from traps such as in mucus and on infected cells. Some NA inhibitors have become widely used drugs for treatment of influenza. However, attempts to develop effective and safe NA inhibitors that can be used for treatment of anti-NA drugs-resistant influenza viruses have continued. In this chapter, we describe the following updates on influenza A NA inhibitor development: (i) N-acetylneuraminic acid (Neu5Ac)-based derivatives, (ii) covalent NA inhibitors, (iii) sulfo-sialic acid analogs, (iv) N-acetyl-6-sulfo-β-D-glucosaminide-based inhibitors, (v) inhibitors targeting the 150-loop of group 1 NAs, (vi) conjugation inhibitors, (vii) acylhydrazone derivatives, (viii) monoclonal antibodies, (ix) PVP-I, and (x) natural products. Finally, we provide future perspectives on the next-generation anti-NA drugs.

Influenza Neuraminidase Characteristics and Potential as a Vaccine Target

Neuraminidase of influenza A and B viruses plays a critical role in the virus life cycle and is an important target of the host immune system. Here, we highlight the current understanding of influenza neuraminidase structure, function, antigenicity, immunogenicity, and immune protective potential. Neuraminidase inhibiting antibodies have been recognized as correlates of protection against disease caused by natural or experimental influenza A virus infection in humans. In the past years, we have witnessed an increasing interest in the use of influenza neuraminidase to improve the protective potential of currently used influenza vaccines. A number of well-characterized influenza neuraminidase-specific monoclonal antibodies have been described recently, most of which can protect in experimental challenge models by inhibiting the neuraminidase activity or by Fc receptor-dependent mechanisms. The relative instability of the neuraminidase poses a challenge for protein-based antigen design. We critically review the different solutions that have been proposed to solve this problem, ranging from the inclusion of stabilizing heterologous tetramerizing zippers to the introduction of inter-protomer stabilizing mutations. Computationally engineered neuraminidase antigens have been generated that offer broad, within subtype protection in animal challenge models. We also provide an overview of modern vaccine technology platforms that are compatible with the induction of robust neuraminidase-specific immune responses. In the near future, we will likely see the implementation of influenza vaccines that confront the influenza virus with a double punch: targeting both the hemagglutinin and the neuraminidase.

Design, Synthesis, Biological Evaluation and In Silico Studies of Pyrazole-Based NH2-Acyl Oseltamivir Analogues as Potent Neuraminidase Inhibitors

Oseltamivir represents one of the most successful neuraminidase (NA) inhibitors in the current anti-influenza therapy. The 150-cavity of NA was identified as an additional binding pocket, and novel NA inhibitors have been designed to occupy the 150-cavity based on the structure information of oseltamivir carboxylate (OC) in complex with NA. In this study, a series of C-5-NH2-acyl derivatives of OC containing the pyrazole moiety were synthesized. Several derivatives exhibited substantial inhibitory activity against NA. Moreover, in silico ADME evaluation indicated that the derivatives were drug-like with higher oral absorption rates and greater cell permeability than OC. Additionally, molecular docking studies revealed that the derivatives interacted with both the NA enzyme active site and 150-cavity as expected. The results provided useful information for further structural optimization of OC.

Structure-based design of 5'-substituted 1,2,3-triazolylated oseltamivir derivatives as potent influenza neuraminidase inhibitors

Resistant viruses containing mutant neuraminidases (NAs) with diminished drug affinity continue to emerge, and new anti-influenza agents are urgently required. Several potent inhibitors targeting the hydrophobic 150-cavity of viral NAs have been developed by modifying the antiviral drugs, oseltamivir carboxylate (OSC) and zanamivir, with hydrophobic groups. Here, we describe a different strategy for exploring novel and efficient NA inhibitors by targeting the charged amino acid residues around the entrance to the 150-cavity. We synthesized a C5-substituted OSC derivative (1e) with a 4'-phenyl-1,2,3-triazolyl group capable of entering the 150-cavity, and solved the crystal structure of 1e in complex with influenza A virus N5 NA. Using the resulting structural information, we next designed and synthesized two series of OSC derivatives carrying various polar substituents at the triazolyl group of 1e and 2e, with 2e being a 5'-phenyl-1,2,3-triazole regioisomer of 1e. The NA inhibition assays demonstrated that the 2 series (2e-n) generally had superior activity compared with the 1 series (1e-n). Compound 2j, bearing a 3-phenylamino group on the triazole ring, was the most potent inhibitor of all tested NAs including an N2 NA containing the E119V OSC-resistant mutation. Moreover, 2j potently inhibited viral replication in vitro, and molecular docking studies revealed that its phenylamino group can form an additional strong hydrogen bond with residue D151 near the entrance of the 150-cavity. The design method described in this study provides useful insights into the development of novel NA inhibitors. Compound 2j warrants further structural optimization to obtain a candidate for clinical use.

Bioinformatic analysis of subfamily-specific regions in 3D-structures of homologs to study functional diversity and conformational plasticity in protein superfamilies

Local 3D-structural differences in homologous proteins contribute to functional diversity observed in a superfamily, but so far received little attention as bioinformatic analysis was usually carried out at the level of amino acid sequences. We have developed Zebra3D - the first-of-its-kind bioinformatic software for systematic analysis of 3D-alignments of protein families using machine learning. The new tool identifies subfamily-specific regions (SSRs) - patterns of local 3D-structure (i.e. single residues, loops, or secondary structure fragments) that are spatially equivalent within families/subfamilies, but are different among them, and thus can be associated with functional diversity and function-related conformational plasticity. Bioinformatic analysis of protein superfamilies by Zebra3D can be used to study 3D-determinants of catalytic activity and specific accommodation of ligands, help to prepare focused libraries for directed evolution or assist development of chimeric enzymes with novel properties by exchange of equivalent regions between homologs, and to characterize plasticity in binding sites. A companion Mustguseal web-server is available to automatically construct a 3D-alignment of functionally diverse proteins, thus reducing the minimal input required to operate Zebra3D to a single PDB code. The Zebra3D + Mustguseal combined approach provides the opportunity to systematically explore the value of SSRs in superfamilies and to use this information for protein design and drug discovery. The software is available open-access at https://biokinet.belozersky.msu.ru/Zebra3D.

Structure-Based Modification of an Anti-neuraminidase Human Antibody Restores Protection Efficacy against the Drifted Influenza Virus

The immune system produces antibodies to protect the human body from harmful invaders. The monoclonal antibody (MAb) is one kind of effective antivirals. In this study, we isolated an antibody (Z2B3) from an H7N9 influenza virus-infected child. It shows cross-reactivity to both group 1 (N1) and group 2 (N9) neuraminidases (NAs) but is sensitive to N1 NA with a K432E substitution. Structural analysis of the NA-antibody fragment antigen-binding (Fab) complex provides a clue for antibody modification, and the modified antibody restored binding and inhibition to recently drifted N1 NA and regained protection against the variant influenza strain. This finding suggests that antibodies to NA may be a useful therapy and can be in principle edited to defeat drifted influenza virus.

Here, we investigate a monoclonal antibody, Z2B3, isolated from an H7N9-infected patient, that exhibited cross-reactivity to both N9 (group 2) and a broad range of seasonal and avian N1 (group 1) proteins but lost activity to the N1 with the substitution K432E. This substitution exists in 99.25% of seasonal influenza strains after 2013. The NA-Z2B3 complex structures indicated that Z2B3 binds within the conserved active site of the neuraminidase (NA) protein. A salt bridge between D102 in Z2B3 and K432 in NA plays an important role in binding. Structure-based modification of Z2B3 with D102R in heavy chain reversed the salt bridge and restored the binding and inhibition of N1 with E432. Furthermore, Z2B3-D102R can protect mice from A/Serbia/NS-601/2014 H1N1 virus (NA contains E432) infection while the wild-type Z2B3 antibody shows no protection. This study demonstrates that a broadly reactive and protective antibody to NA can be in principle edited to restore binding and inhibition to recently drifted N1 NA and regain protection against the variant influenza strain.

Identification of Chebulinic Acid and Chebulagic Acid as Novel Influenza Viral Neuraminidase Inhibitors

The influenza A virus (IAV) causes seasonal epidemics and occasional but devastating pandemics, which are of a major public health concern. Although several antiviral drugs are currently available, there is an urgent need to develop novel antiviral therapies with different mechanisms of action due to emergence of drug resistance. In this study, two related compounds, chebulagic acid (CHLA) and chebulinic acid (CHLI), were identified as novel inhibitors against IAV replication. A reporter virus-based infection assay demonstrated that CHLA and CHLI exhibit no inhibitory effect on IAV entry or RNA replication during the virus replication cycle. Results of viral release inhibition assay and neuraminidase (NA) inhibition assay indicated that CHLA and CHLI exert their inhibitory effect on the NA-mediated viral release. Moreover, oseltamivir-resistance mutation NA/H274Y of NA is susceptible to CHLA or CHLI, suggesting a different mechanism of action for CHLA and CHLI. In summary, CHLA and CHLI are promising new NA inhibitors that may be further developed as novel antivirals against IAVs.

Interactive molecular dynamics in virtual reality for accurate flexible protein-ligand docking

Simulating drug binding and unbinding is a challenge, as the rugged energy landscapes that separate bound and unbound states require extensive sampling that consumes significant computational resources. Here, we describe the use of interactive molecular dynamics in virtual reality (iMD-VR) as an accurate low-cost strategy for flexible protein-ligand docking. We outline an experimental protocol which enables expert iMD-VR users to guide ligands into and out of the binding pockets of trypsin, neuraminidase, and HIV-1 protease, and recreate their respective crystallographic protein-ligand binding poses within 5-10 minutes. Following a brief training phase, our studies shown that iMD-VR novices were able to generate unbinding and rebinding pathways on similar timescales as iMD-VR experts, with the majority able to recover binding poses within 2.15 Å RMSD of the crystallographic binding pose. These results indicate that iMD-VR affords sufficient control for users to carry out the detailed atomic manipulations required to dock flexible ligands into dynamic enzyme active sites and recover crystallographic poses, offering an interesting new approach for simulating drug docking and generating binding hypotheses.

Source of oseltamivir resistance due to single E119D and double E119D/H274Y mutations in pdm09H1N1 influenza neuraminidase

Influenza epidemics are responsible for an average of 3-5 millions of severe cases and up to 500,000 deaths around the world. One of flu pandemic types is influenza A(H1N1)pdm09 virus (pdm09H1N1). Oseltamivir is the antiviral drug used to treat influenza targeting at neuraminidase (NA) located on the viral surface. Influenza virus undergoes high mutation rates and leads to drug resistance, and thus the development of more efficient drugs is required. In the present study, all-atom molecular dynamics simulations were applied to understand the oseltamivir resistance caused by the single E119D and double E119D/H274Y mutations on NA. The obtained results in terms of binding free energy and intermolecular interactions in the ligand-protein interface showed that the oseltamivir could not be well accommodated in the binding pocket of both NA mutants and the 150-loop moves out from oseltamivir as an "open" state. A greater number of water molecules accessible to the binding pocket could disrupt the oseltamivir binding with NA target as seen be high mobility of oseltamivir at the active site. Additionally, our finding could guide to the design and development of novel NA inhibitor drugs.

Design, synthesis and biological evaluation of oseltamivir derivatives containing pyridyl group as potent inhibitors of neuraminidase for influenza A

Influenza A neuraminidase plays an indispensable role in the process of replication and transmission of influenza, so the neuraminidase inhibition can prevent the reproduction of the viruses therefore achieve the effect of treatment of influenza. However, drug resistance of neuraminidase inhibitors such as oseltamivir highlights the need to develop novel structural neuraminidase inhibitors. Here we explored a series of oseltamivir derivatives bearing pyridyl group. Among them, compound 23b exhibiting potent inhibitory activity against neuraminidase from H5N1 subtype was comparable to oseltamivir carboxylate. Cytopathic effect inhibition assay in MDCK cells indicated that compound 23b exerted powerful inhibitions on influenza viruses. And compound 23b were nontoxic to MDCK cells. Meanwhile, compound 23b showed high stability towards rat liver microsomes, human liver microsomes and human plasma. This research enriched the structural type of neuraminidase inhibitors.

Design, synthesis and biological evaluation of "Multi-Site"-binding influenza virus neuraminidase inhibitors

Encouraged by our earlier discovery of neuraminidase inhibitors targeting 150-cavity or 430-cavity, herein, to yield more potent inhibitors, we designed, synthesized, and biologically evaluated a series of novel oseltamivir derivatives via modification of C-1 and C5-NH2 of oseltamivir by exploiting 150-cavity and/or 430-cavity. Among the synthesized compounds, compound 15e, the most potent N1-selective inhibitor targeting 150-cavity, showed 1.5 and 1.8 times greater activity than oseltamivir carboxylate (OSC) against N1 (H5N1) and N1 (H5N1-H274Y). In cellular assays, 15e also exhibited greater potency than OSC against H5N1 with EC50 of 0.66 μM. In addition, 15e demonstrated low cytotoxicity in vitro and low acute toxicity in mice. Molecular docking studies provided insights into the high potency of 15e against N1 and N1-H274Y mutant NA. Overall, we envisioned that the significant breakthrough in the discovery of potent group-1-specific neuraminidase inhibitors may lead to further investigation of more potent anti-influenza agents.

Novel N-Substituted oseltamivir derivatives as potent influenza neuraminidase inhibitors: Design, synthesis, biological evaluation, ADME prediction and molecular docking studies

The discovery of novel potent neuraminidase (NA) inhibitors remains an attractive approach for treating infectious diseases caused by influenza. In this study, we describe the design and synthesis of novel N-substituted oseltamivir derivatives for probing the 150-cavity which is nascent to the activity site of NA. NA inhibitory studies showed that new derivatives demonstrated the inhibitory activity with IC₅₀ values at nM level against NA of a clinical influenza virus strain. Moreover, the in silico ADME predictions showed that the selected compounds had comparable properties with oseltamivir carboxylate, which demonstrated the druggablity of these derivatives. Furthermore, molecular docking studies showed that the most potent compound 6f and 10i could adopt different modes of binding interaction with NA, which may provide novel solutions for treating oseltamivir-resistant influenza. Based on the research results, we consider that compounds 6f and 10i have the potential for further studies as novel antiviral agents.

Impact of tetramerization on the ligand recognition of N1 influenza neuraminidase via MMGBSA approach

Influenza virus neuraminidase (NA) is a homotetrameric surface protein that, in contrast to other non-influenza NAs, requires a quaternary assembly to exhibit enzymatic activity, suggesting that the oligomeric state significantly impacts the active site of influenza NA. Nevertheless, most structure-based drug design studies have been reported by employing the monomeric state in the closed or open-loop due to the computational cost of employing the tetrameric NA. In this work, we present MD simulations coupled to the MMGBSA approach of avian N1 type NA in its monomeric and tetrameric closed and open-loop state both with and without the inhibitor oseltamivir and its natural substrate, sialic acid. Structural and energetic analyses revealed that the tetrameric state impacts flexibility as well as the map of interactions participating in stabilizing the protein-ligand complexes with respect to the monomeric state. It was observed that the tetrameric state exerts dissimilar effects in binding affinity, characteristic of positive and negative cooperativity for oseltamivir and sialic acid, respectively. Based on our results, to perform a confident structure-based drug design, as well as to evaluate the impact of key mutations through MD simulations, it is important to consider the tetrameric state closed-loop state.

Novel antiviral drug discovery strategies to tackle drug-resistant mutants of influenza virus strains

INTRODUCTION

The emergence of drug-resistant influenza virus strains highlights the need for new antiviral therapeutics to combat future pandemic outbreaks as well as continuing seasonal cycles of influenza. Areas covered: This review summarizes the mechanisms of current FDA-approved anti-influenza drugs and patterns of resistance to those drugs. It also discusses potential novel targets for broad-spectrum antiviral drugs and recent progress in novel drug design to overcome drug resistance in influenza. Expert opinion: Using the available structural information about drug-binding pockets, research is currently underway to identify molecular interactions that can be exploited to generate new antiviral drugs. Despite continued efforts, antivirals targeting viral surface proteins like HA, NA, and M2, are all susceptible to developing resistance. Structural information on the internal viral polymerase complex (PB1, PB2, and PA) provides a new avenue for influenza drug discovery. Host factors, either at the initial step of viral infection or at the later step of nuclear trafficking of viral RNP complex, are being actively pursued to generate novel drugs with new modes of action, without resulting in drug resistance.

DNA-linked inhibitor antibody assay (DIANA) as a new method for screening influenza neuraminidase inhibitors

Influenza neuraminidase is responsible for the escape of new viral particles from the infected cell surface. Several neuraminidase inhibitors are used clinically to treat patients or stockpiled for emergencies. However, the increasing development of viral resistance against approved inhibitors has underscored the need for the development of new antivirals effective against resistant influenza strains. A facile, sensitive, and inexpensive screening method would help achieve this goal. Recently, we described a multiwell plate-based DNA-linked inhibitor antibody assay (DIANA). This highly sensitive method can quantify femtomolar concentrations of enzymes. DIANA also has been applied to high-throughput enzyme inhibitor screening, allowing the evaluation of inhibition constants from a single inhibitor concentration. Here, we report the design, synthesis, and structural characterization of a tamiphosphor derivative linked to a reporter DNA oligonucleotide for the development of a DIANA-type assay to screen potential influenza neuraminidase inhibitors. The neuraminidase is first captured by an immobilized antibody, and the test compound competes for binding to the enzyme with the oligo-linked detection probe, which is then quantified by qPCR. We validated this novel assay by comparing it with the standard fluorometric assay and demonstrated its usefulness for sensitive neuraminidase detection as well as high-throughput screening of potential new neuraminidase inhibitors.

Optimization of N-Substituted Oseltamivir Derivatives as Potent Inhibitors of Group-1 and -2 Influenza A Neuraminidases, Including a Drug-Resistant Variant

On the basis of our earlier discovery of N1-selective inhibitors, the 150-cavity of influenza virus neuraminidases (NAs) could be further exploited to yield more potent oseltamivir derivatives. Among the synthesized compounds, 15b and 15c were exceptionally active against both group-1 and -2 NAs. Especially for 09N1, N2, N6, and N9 subtypes, they showed 6.80-12.47 and 1.20-3.94 times greater activity than oseltamivir carboxylate (OSC). They also showed greater inhibitory activity than OSC toward H274Y and E119V variant. In cellular assays, they exhibited greater potency than OSC toward H5N1, H5N2, H5N6, and H5N8 viruses. 15b demonstrated high metabolic stability, low cytotoxicity in vitro, and low acute toxicity in mice. Computational modeling and molecular dynamics studies provided insights into the role of R group of 15b in improving potency toward group-1 and -2 NAs. We believe the successful exploitation of the 150-cavity of NAs represents an important breakthrough in the development of more potent anti-influenza agents.

Influenza Infection in Humans Induces Broadly Cross-Reactive and Protective Neuraminidase-Reactive Antibodies

Antibodies to the hemagglutinin (HA) and neuraminidase (NA) glycoproteins are the major mediators of protection against influenza virus infection. Here, we report that current influenza vaccines poorly display key NA epitopes and rarely induce NA-reactive B cells. Conversely, influenza virus infection induces NA-reactive B cells at a frequency that approaches (H1N1) or exceeds (H3N2) that of HA-reactive B cells. NA-reactive antibodies display broad binding activity spanning the entire history of influenza A virus circulation in humans, including the original pandemic strains of both H1N1 and H3N2 subtypes. The antibodies robustly inhibit the enzymatic activity of NA, including oseltamivir-resistant variants, and provide robust prophylactic protection, including against avian H5N1 viruses, in vivo. When used therapeutically, NA-reactive antibodies protected mice from lethal influenza virus challenge even 48 hr post infection. These findings strongly suggest that influenza vaccines should be optimized to improve targeting of NA for durable and broad protection against divergent influenza strains.

Exploration of the 150 cavity and the role of serendipity in the discovery of inhibitors of influenza virus A neuraminidase

Influenza pandemics are an ongoing threat for the human population, as the avian influenza viruses H5N1 and H7N9 continue to circulate in the bird population and the chance of avian to human transmission increases. Neuraminidase, a glycoprotein located on the surface of the influenza virus, plays a crucial role in the viral replication process and, hence, has proven to be a useful target enzyme for the treatment of influenza infections. The discovery that certain subtypes of influenza neuraminidase have an additional cavity, the 150 cavity, near the substrate binding site has triggered considerable interest in the design of influenza inhibitors that exploit this feature. Currently available antiviral drugs, neuraminidase inhibitors oseltamivir and zanamivir, were designed using crystal structures predating this discovery by some years. This mini review is aimed at summarizing our group’s efforts, together with related work from other groups, on neuraminidase inhibitors that are designed to exploit both the catalytic site and the 150 cavity. The design of a parent scaffold that yields a potent inhibitor that is active in cell culture assays and retains activity against several neuraminidases from mutant strains is also described. Finally, the role of serendipity in the discovery of a new class of potent neuraminidase inhibitors with a novel spirolactam scaffold is also highlighted.

Screening for Neuraminidase Inhibitor Resistance Markers among Avian Influenza Viruses of the N4, N5, N6, and N8 Neuraminidase Subtypes

Several subtypes of avian influenza viruses (AIVs) are emerging as novel human pathogens, and the frequency of related infections has increased in recent years. Although neuraminidase (NA) inhibitors (NAIs) are the only class of antiviral drugs available for therapeutic intervention for AIV-infected patients, studies on NAI resistance among AIVs have been limited, and markers of resistance are poorly understood. Previously, we identified unique NAI resistance substitutions in AIVs of the N3, N7, and N9 NA subtypes. Here, we report profiles of NA substitutions that confer NAI resistance in AIVs of the N4, N5, N6, and N8 NA subtypes using gene-fragmented random mutagenesis. We generated libraries of mutant influenza viruses using reverse genetics (RG) and selected resistant variants in the presence of the NAIs oseltamivir carboxylate and zanamivir in MDCK cells. In addition, two substitutions, H274Y and R292K (N2 numbering), were introduced into each NA gene for comparison. We identified 37 amino acid substitutions within the NA gene, 16 of which (4 in N4, 4 in N5, 4 in N6, and 4 in N8) conferred resistance to NAIs (oseltamivir carboxylate, zanamivir, or peramivir) as determined using a fluorescence-based NA inhibition assay. Substitutions conferring NAI resistance were mainly categorized as either novel NA subtype specific (G/N147V/I, A246V, and I427L) or previously reported in other subtypes (E119A/D/V, Q136K, E276D, R292K, and R371K). Our results demonstrate that each NA subtype possesses unique NAI resistance markers, and knowledge of these substitutions in AIVs is important in facilitating antiviral susceptibility monitoring of NAI resistance in AIVs.IMPORTANCE The frequency of human infections with avian influenza viruses (AIVs) has increased in recent years. Despite the availability of vaccines, neuraminidase inhibitors (NAIs), as the only available class of drugs for AIVs in humans, have been constantly used for treatment, leading to the inevitable emergence of drug-resistant variants. To screen for substitutions conferring NAI resistance in AIVs of N4, N5, N6, and N8 NA subtypes, random mutations within the target gene were generated, and resistant viruses were selected from mutant libraries in the presence of individual drugs. We identified 16 NA substitutions conferring NAI resistance in the tested AIV subtypes; some are novel and subtype specific, and others have been previously reported in other subtypes. Our findings will contribute to an increased and more comprehensive understanding of the mechanisms of NAI-induced inhibition of influenza virus and help lead to the development of drugs that bind to alternative interaction motifs.

Exploration of binding and inhibition mechanism of a small molecule inhibitor of influenza virus H1N1 hemagglutinin by molecular dynamics simulation

Influenza viruses are a major public health threat worldwide. The influenza hemagglutinin (HA) plays an essential role in the virus life cycle. Due to the high conservation of the HA stem region, it has become an especially attractive target for inhibitors for therapeutics. In this study, molecular simulation was applied to study the mechanism of a small molecule inhibitor (MBX2329) of influenza HA. Behaviors of the small molecule under neutral and acidic conditions were investigated, and an interesting dynamic binding mechanism was found. The results suggested that the binding of the inhibitor with HA under neutral conditions facilitates only its intake, while it interacts with HA under acidic conditions using a different mechanism at a new binding site. After a series of experiments, we believe that binding of the inhibitor can prevent the release of HA1 from HA2, further maintaining the rigidity of the HA2 loop and stabilizing the distance between the long helix and short helices. The investigated residues in the new binding site show high conservation, implying that the new binding pocket has the potential to be an effective drug target. The results of this study will provide a theoretical basis for the mechanism of new influenza virus inhibitors.

Neuraminidase-mediated haemagglutination of recent human influenza A(H3N2) viruses is determined by arginine 150 flanking the neuraminidase catalytic site

Over the last decade, an increasing proportion of circulating human influenza A(H3N2) viruses exhibited haemagglutination activity that was sensitive to neuraminidase inhibitors. This change in haemagglutination as compared to older circulating A(H3N2) viruses prompted an investigation of the underlying molecular basis. Recent human influenza A(H3N2) viruses were found to agglutinate turkey erythrocytes in a manner that could be blocked with either oseltamivir or neuraminidase-specific antisera, indicating that agglutination was driven by neuraminidase, with a low or negligible contribution of haemagglutinin. Using representative virus recombinants it was shown that the haemagglutinin of a recent A(H3N2) virus indeed had decreased activity to agglutinate turkey erythrocytes, while its neuraminidase displayed increased haemagglutinating activity. Viruses with chimeric and mutant neuraminidases were used to identify the amino acid substitution histidine to arginine at position 150 flanking the neuraminidase catalytic site as the determinant of this neuraminidase-mediated haemagglutination. An analysis of publicly available neuraminidase gene sequences showed that viruses with histidine at position 150 were rapidly replaced by viruses with arginine at this position between 2005 and 2008, in agreement with the phenotypic data. As a consequence of neuraminidase-mediated haemagglutination of recent A(H3N2) viruses and poor haemagglutination via haemagglutinin, haemagglutination inhibition assays with A(H3N2) antisera are no longer useful to characterize the antigenic properties of the haemagglutinin of these viruses for vaccine strain selection purposes. Continuous monitoring of the evolution of these viruses and potential consequences for vaccine strain selection remains important.

Influenza A(H1N1)pdm09 virus exhibiting enhanced cross-resistance to oseltamivir and peramivir due to a dual H275Y/G147R substitution, Japan, March 2016

An influenza A(H1N1)pdm09 virus carrying a G147R substitution in combination with an H275Y substitution in the neuraminidase protein, which confers cross-resistance to oseltamivir and peramivir, was detected from an immunocompromised inpatient in Japan, March 2016. This dual H275Y/G147R mutant virus exhibited enhanced cross-resistance to both drugs compared with the single H275Y mutant virus and reduced susceptibility to zanamivir, although it showed normal inhibition by laninamivir.

Potential New H1N1 Neuraminidase Inhibitors from Ferulic Acid and Vanillin: Molecular Modelling, Synthesis and in Vitro Assay

We report the computational and experimental efforts in the design and synthesis of novel neuraminidase (NA) inhibitors from ferulic acid and vanillin. Two proposed ferulic acid analogues, MY7 and MY8 were predicted to inhibit H1N1 NA using molecular docking. From these two analogues, we designed, synthesised and evaluated the biological activities of a series of ferulic acid and vanillin derivatives. The enzymatic H1N1 NA inhibition assay showed MY21 (a vanillin derivative) has the lowest IC50 of 50 μM. In contrast, the virus inhibition assay showed MY15, a ferulic acid derivative has the best activity with the EC50 of ~0.95 μM. Modelling studies further suggest that these predicted activities might be due to the interactions with conserved and essential residues of NA with ΔGbind values comparable to those of oseltamivir and zanamivir, the two commercial NA inhibitors.

Molecular Basis for Differential Patterns of Drug Resistance in Influenza N1 and N2 Neuraminidase

Neuraminidase (NA) inhibitors are used for the prevention and treatment of influenza A virus infections. Two subtypes of NA, N1 and N2, predominate in viruses that infect humans, but differential patterns of drug resistance have emerged in each subtype despite highly homologous active sites. To understand the molecular basis for the selection of these drug resistance mutations, structural and dynamic analyses on complexes of N1 and N2 NA with substrates and inhibitors were performed. Comparison of dynamic substrate and inhibitor envelopes and interactions at the active site revealed how differential patterns of drug resistance have emerged for specific drug resistance mutations, at residues I222, S246, and H274 in N1 and E119 in N2. Our results show that the differences in intermolecular interactions, especially van der Waals contacts, of the inhibitors versus substrates at the NA active site effectively explain the selection of resistance mutations in the two subtypes. Avoiding such contacts that render inhibitors vulnerable to resistance by better mimicking the dynamics and intermolecular interactions of substrates can lead to the development of novel inhibitors that avoid drug resistance in both subtypes.

Resistance to Mutant Group 2 Influenza Virus Neuraminidases of an Oseltamivir-Zanamivir Hybrid Inhibitor

Influenza virus neuraminidase (NA) drug resistance is one of the challenges to preparedness against epidemic and pandemic influenza virus infections. NA N1- and N2-containing influenza viruses are the primary cause of seasonal epidemics and past pandemics. The structural and functional basis underlying drug resistance of the influenza virus N1 NA is well characterized. Yet drug resistance of the N2 strain is not well understood. Here, we confirm that replacement of N2 E119 or I222 results in multidrug resistance, and when the replacements occur together, the sensitivity to NA inhibitors (NAI) is reduced severely. Using crystallographic studies, we showed that E119 replacement results in a loss of hydrogen bonding to oseltamivir and zanamivir, whereas I222 replacement results in a change in the hydrophobic environment that is critical for oseltamivir binding. Moreover, we found that MS-257, a zanamivir-oseltamivir hybrid inhibitor, is less susceptible to drug resistance. The binding mode of MS-257 shows that increased hydrogen bonding interactions between the inhibitor and NA active site anchor the inhibitor within the active site and allow adjustments in response to active-site modifications. Such stability is likely responsible for the observed reduced susceptibility to drug resistance. MS-257 serves as a next-generation anti-influenza virus drug candidate and serves also as a scaffold for further design of NAIs.

IMPORTANCE

Oseltamivir and zanamivir are the two major antiviral drugs available for the treatment of influenza virus infections. However, multidrug-resistant viruses have emerged in clinical cases, which pose a challenge for the development of new drugs. N1 and N2 subtypes exist in the viruses which cause seasonal epidemics and past pandemics. Although N1 drug resistance is well characterized, the molecular mechanisms underlying N2 drug resistance are unknown. A previous report showed that an N2 E119V/I222L dual mutant conferred drug resistance to seasonal influenza virus. Here, we confirm that these substitutions result in multidrug resistance and dramatically reduced sensitivity to NAI. We further elucidate the molecular mechanism underlying N2 drug resistance by solving crystal structures of the N2 E119V and I222L mutants and the dual mutant. Most importantly, we found that a novel oseltamivir-zanamivir hybrid inhibitor, MS-257, remains more effective against drug-resistant N2 and is a promising candidate as a next-generation anti-influenza virus drug.

Design and synthesis of 1,2,3-triazole-containing N-acyl zanamivir analogs as potent neuraminidase inhibitors

The design of potent metabolically stable neuraminidase (NA) inhibitors represents an attractive approach for treating influenza virus infection. In this study, we describe the exploitation of the 150-cavity in the active site of group 1 NA for the design, synthesis, and in vitro evaluation of new triazole-containing N-acyl derivatives related to Zanamivir. Inhibition studies with influenza virus NAs of group 1 (H1N1) and group 2 (H3N2) revealed that several of them are good inhibitors, with IC50 values in the low nanomolar (2.3 nM-31 nM) range. Substituents that form stable van der Waals interaction with the 150-cavity residues play crucial roles in NA inhibition as demonstrated by the potency of 6a (H1N1 IC50 = 2.3 nM, and H3N2 IC50 = 2.9 nM). Docking studies indicated that the cyclohexane-substituted triazole ring extended toward the hydrophobic region in the active site of group 1 NA in open form. The high potency observed for inhibitor 6a may be attributable to the highly favorable hydrophobic interactions in this region.

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

Deletions in the stalk of the influenza neuraminidase (NA) surface protein are associated with increased virulence, but the mechanisms responsible for this enhanced virulence are unclear. Here we use microsecond molecular dynamics simulations to explore the effect of stalk deletion on enzymatic activity, contrasting NA proteins from the A/swine/Shandong/N1/2009 strain both with and without a stalk deletion. By modeling and simulating neuraminidase apo glycoproteins embedded in complex-mixture lipid bilayers, we show that the geometry and dynamics of the neuraminidase enzymatic pocket may differ depending on stalk length, with possible repercussions on the binding of the endogenous sialylated-oligosaccharide receptors. We also use these simulations to predict previously unrecognized druggable “hotspots” on the neuraminidase surface that may prove useful for future efforts aimed at structure-based drug design.

Phylogenetic analysis of human influenza A/H3N2 viruses isolated in 2015 in Germany indicates significant genetic divergence from vaccine strains

Infections by H3N2-type influenza A viruses (IAV) resulted in significant numbers of hospitalization in several countries in 2014-2015, causing disease also in vaccinated individuals and, in some cases, fatal outcomes. In this study, sequence analysis of H3N2 viruses isolated in Germany from 1998 to 2015, including eleven H3N2 isolates collected early in 2015, was performed. Compared to the vaccine strain A/Texas/50/2012 (H3N2), the 2015 strains from Germany showed up to 4.5 % sequence diversity in their HA1 protein, indicating substantial genetic drift. The data further suggest that two distinct phylogroups, 3C.2 and 3C.3, with 1.6-2.3 % and 0.3-2.4 % HA1 nucleotide and amino acid sequence diversity, respectively, co-circulated in Germany in the 2014/2015 season. Distinct glycosylation patterns and amino acid substitutions in the hemagglutinin and neuraminidase proteins were identified, possibly contributing to the unusually high number of H3N2 infections in this season and providing important information for developing vaccines that are effective against both genotypes.

Novel binding patterns between ganoderic acids and neuraminidase: Insights from docking, molecular dynamics and MM/PBSA studies

Recently, ganoderic acids (GAs) give rise to the attractive candidates of novel neuraminidase (NA) inhibitors. However, there is still no evident conclusion about their binding patterns. To this end, docking, molecular dynamics and MM/PBSA methods were combined to study the binding profiles of GAs with the N1 protein and familiar H274Y and N294S mutations (A/Vietnam/1203/04 stain). It was found that the binding affinities of ganoderic acid DM and Z (ΔGbind, -16.83 and -10.99 kcal mol(-1)) are comparable to that of current commercial drug oseltamivir (-23.62 kcal mol(-1)). Electrostatic interaction is the main driving force, and should be one important factor to evaluate the binding quality and rational design of NA inhibitors. The 150-loop residues Asp151 and Arg152 played an important role in the binding processes. Further analysis revealed that ganoderic acid DM is a potential source of anti-influenza ingredient, with novel binding pattern and advantage over oseltamivir. It had steric hindrance on the 150 cavity of N1 protein, and exerted activities across the H274Y and N294S mutations. This work also pointed out how to effectively design dual-site NA inhibitors and reinforce their affinities. These findings should prove valuable for the in-depth understanding of interactions between NA and GAs, and warrant the experimental aspects to design novel anti-influenza drugs.

Applying high-performance computing in drug discovery and molecular simulation

In recent decades, high-performance computing (HPC) technologies and supercomputers in China have significantly advanced, resulting in remarkable achievements. Computational drug discovery and design, which is based on HPC and combines pharmaceutical chemistry and computational biology, has become a critical approach in drug research and development and is financially supported by the Chinese government. This approach has yielded a series of new algorithms in drug design, as well as new software and databases. This review mainly focuses on the application of HPC to the fields of drug discovery and molecular simulation at the Chinese Academy of Sciences, including virtual drug screening, molecular dynamics simulation, and protein folding. In addition, the potential future application of HPC in precision medicine is briefly discussed.

Synthesis and evaluation of influenza A viral neuraminidase candidate inhibitors based on a bicyclo[3.1.0]hexane scaffold

We describe the synthesis of constrained oseltamivir analogues designed to mimic the proposed boat conformation of the enzymatic transition state.

Discovery of Influenza A virus neuraminidase inhibitors using support vector machine and Naïve Bayesian models

Neuraminidase (NA) is a critical enzyme in the life cycle of influenza virus, which is known as a successful paradigm in the design of anti-influenza agents. However, to date there are no classification models for the virtual screening of NA inhibitors. In this work, we built support vector machine and Naïve Bayesian models of NA inhibitors and non-inhibitors, with different ratios of active-to-inactive compounds in the training set and different molecular descriptors. Four models with sensitivity or Matthews correlation coefficients greater than 0.9 were chosen to predict the NA inhibitory activities of 15,600 compounds in our in-house database. We combined the results of four optimal models and selected 60 representative compounds to assess their NA inhibitory profiles in vitro. Nine NA inhibitors were identified, five of which were oseltamivir derivatives with large C-5 substituents exhibiting potent inhibition against H1N1 NA with IC50 values in the range of 12.9-185.0 nM, and against H3N2 NA with IC50 values between 18.9 and 366.1 nM. The other four active compounds belonged to novel scaffolds, with IC50 values ranging 39.5-63.8 μM against H1N1 NA and 44.5-114.1 μM against H3N2 NA. This is the first time that classification models of NA inhibitors and non-inhibitors are built and their prediction results validated experimentally using in vitro assays.