Dual-Action Heteromultivalent Glycopolymers Stringently Block and Arrest Influenza A Virus Infection In Vitro and Ex Vivo.

Here, we demonstrate the concerted inhibition of different influenza A virus (IAV) strains using a low-molecular-weight dual-action linear polymer. The 6'-sialyllactose and zanamivir conjugates of linear polyglycerol are optimized for simultaneous targeting of hemagglutinin and neuraminidase on the IAV surface. Independent of IAV subtypes, hemagglutination inhibition data suggest better adsorption of the heteromultivalent polymer than homomultivalent analogs onto the virus surface. Cryo-TEM images imply heteromultivalent compound-mediated virus aggregation. The optimized polymeric nanomaterial inhibits >99.9% propagation of various IAV strains 24 h postinfection in vitro at low nM concentrations and is up to 10000× more effective than the commercial zanamivir drug. In a human lung ex vivo multicyclic infection setup, the heteromultivalent polymer outperforms the commercial drug zanamivir and homomultivalent analogs or their physical mixtures. This study authenticates the translational potential of the dual-action targeting approach using small polymers for broad and high antiviral efficacy.

N1 neuraminidase of H5N1 avian influenza A virus complexed with sialic acid and zanamivir - A study by molecular docking and molecular dynamics simulation.

Development of antiviral drugs is an urgent need to control and prevent the presently circulating H5N1 avian influenza virus which is affects the human respiratory tract. The complex crystal structure of N1-N-acetylneuranamic acid (sialic acid, SIA) is not available as complex and hence SIA and zanamivir (ZMR) are docked into the binding site of N1 neuraminidase. Based on the analysis, the initial complex structures have been simulated for 120 ns to get insight into the binding modes and interaction between protein-ligand complex systems. NAMD pair interaction energy and MM-PBSA binding free energy are calculated and show that there are two possible binding modes (BM1 and BM2) for N1-SIA and a single binding mode (BM1) for and N1-ZMR complex structures respectively. BM1 of N1-SIA is the most preferred binding mode. On contrary to the currently available drugs in which the chair conformation is distorted, in both the binding modes of N1-SIA, the binding pocket of N1 neuraminidase is able to accommodate SIA in 2C5 chair conformation which is the preferred conformation of SIA in solution state. In N1-ZMR complex, ZMR is bind in a distorted chair conformation. The neuraminidase binding pocket is also able to accommodate galactose of SIAα(2→3)GAL and SIAα(2→6)GAL. RMSD, RMSF and hydrogen bonding analyses have been carried out to identify the conformational flexibility and structural stability of each complex system. All the analyses show that SIA can be used as an inhibitor for N1 neuraminidase of H5N1 influenza viral infection. Communicated by Ramaswamy H. Sarma.

Zanamivir-Cholesterol Conjugate: A Long-Acting Neuraminidase Inhibitor with Potent Efficacy against Drug-Resistant Influenza Viruses.

Antiviral therapy of influenza virus infections depends heavily on two viral neuraminidase (NA) inhibitors, oseltamivir (OSV) and zanamivir (ZNV). The efficacy of OSV is challenged by the development of viral resistance, while the clinical use of ZNV is limited by its poor pharmacokinetic profile and requirement for twice-daily intranasal administration. We have developed a novel NA inhibitor by conjugating ZNV to cholesterol. The ZNV-cholesterol conjugate showed markedly improved antiviral efficacy and plasma half-life compared with ZNV. Single-dose administration of the conjugate protected the mice from lethal challenges with wild-type or mutant H1N1 influenza viruses bearing an OSV-resistant H275Y-substitution. Mechanistic studies showed that the conjugate targeted the cell membrane and entered the host cells, thereby inhibiting the NA function and the assembly of progeny virions. The ZNV-cholesterol conjugate represents a potential new treatment for influenza infections with sustained effect. Cholesterol conjugation may be an effective strategy for improving the pharmacokinetics and efficacy of other small-molecule therapeutics.

In Pursuit of Synthetic Efficiency: Convergent Approaches.

The field of total synthesis has reached a stage in which emphasis has been increasingly focused on synthetic efficiency rather than merely achieving the synthesis of a target molecule. The pursuit of synthetic efficiency, typically represented by step count and overall yield, is a rich source of inspiration and motivation for synthetic chemists to invent innovative strategies and methods. Among them, convergent strategy has been well recognized as an effective approach to improve efficiency. This strategy generally involves coupling of fragments with similar complexity to furnish the target molecule via subsequent cyclization or late-stage functionalization. Thus, methodologies that enable effective connection of fragments are critical to devising a convergent plan. In our laboratory, convergent strategy has served as a long-standing principle for pursuing efficient synthesis during the course of planning and implementing synthetic projects. In this Account, we summarize our endeavors in the convergent synthesis of natural products over the last ten years. We show how we identify reasonable bond disconnections and employ enabling synthetic methodologies to maximize convergency, leading to the efficient syntheses of over two-dozen highly complex molecules from eight disparate families.In detail, we categorize our work into three parts based on the diverse reaction types for fragment assembly. First, we demonstrate the application of a powerful single-electron reducing agent, SmI2, in a late-stage cyclization step, forging the polycyclic skeletons of structurally fascinating Galbulimima alkaloids and Leucosceptrum sesterterpenoids. Next, we showcase how three different types of cycloaddition reactions can simultaneously construct two challenging C-C bonds in a single step, providing concise entries to three distinct families, namely, spiroquinazoline alkaloids, gracilamine, and kaurane diterpenoids. In the third part, we describe convergent assembly of ent-kaurane diterpenoids, gelsedine-type alkaloids, and several drug molecules via employing some bifunctional synthons. To access highly oxidized ent-kaurane diterpenoids, we introduce the hallmark bicyclo[3.2.1]octane ring system at an early stage, and then execute coupling and cyclization by means of a Hoppe's homoaldol reaction and a Mukaiyama-Michael-type addition, respectively. Furthermore, we showcase how the orchestrated combination of an asymmetric Michael addition, a tandem oxidation-aldol reaction and a pinacol rearrangement can dramatically improve the efficiency in synthesizing gelsedine-type alkaloids, with nary a protecting group. Finally, to address the supply issue of several drugs, including anti-influenza drug zanamivir and antitumor agent Et-743, we exploit scalable and practical approaches to provide advantages over current routes in terms of cost, ease of execution, and efficiency.

A universal dual mechanism immunotherapy for the treatment of influenza virus infections.

Seasonal influenza epidemics lead to 3-5 million severe infections and 290,000-650,000 annual global deaths. With deaths from the 1918 influenza pandemic estimated at >50,000,000 and future pandemics anticipated, the need for a potent influenza treatment is critical. In this study, we design and synthesize a bifunctional small molecule by conjugating the neuraminidase inhibitor, zanamivir, with the highly immunogenic hapten, dinitrophenyl (DNP), which specifically targets the surface of free virus and viral-infected cells. We show that this leads to simultaneous inhibition of virus release, and immune-mediated elimination of both free virus and virus-infected cells. Intranasal or intraperitoneal administration of a single dose of drug to mice infected with 100x MLD50 virus is shown to eradicate advanced infections from representative strains of both influenza A and B viruses. Since treatments of severe infections remain effective up to three days post lethal inoculation, our approach may successfully treat infections refractory to current therapies.

Five Novel Non-Sialic Acid-Like Scaffolds Inhibit In Vitro H1N1 and H5N2 Neuraminidase Activity of Influenza a Virus.

Neuraminidase (NA) of influenza viruses enables the virus to access the cell membrane. It degrades the sialic acid contained in extracellular mucin. Later, it is responsible for releasing newly formed virions from the membrane of infected cells. Both processes become key functions within the viral cycle. Therefore, it is a therapeutic target for research of the new antiviral agents. Structure-activity relationships studies have revealed which are the important functional groups for the receptor-ligand interaction. Influenza virus type A NA activity was inhibited by five scaffolds without structural resemblance to sialic acid. Intending small organic compound repositioning along with drug repurposing, this study combined in silico simulations of ligand docking into the known binding site of NA, along with in vitro bioassays. The five proposed scaffolds are N-acetylphenylalanylmethionine, propanoic 3-[(2,5-dimethylphenyl) carbamoyl]-2-(piperazin-1-yl) acid, 3-(propylaminosulfonyl)-4-chlorobenzoic acid, ascorbic acid (vitamin C), and 4-(dipropylsulfamoyl) benzoic acid (probenecid). Their half maximal inhibitory concentration (IC50) was determined through fluorometry. An acidic reagent 2'-O-(4-methylumbelliferyl)-α-dN-acetylneuraminic acid (MUNANA) was used as substrate for viruses of human influenza H1N1 or avian influenza H5N2. Inhibition was observed in millimolar ranges in a concentration-dependent manner. The IC50 values of the five proposed scaffolds ranged from 6.4 to 73 mM. The values reflect a significant affinity difference with respect to the reference drug zanamivir (p < 0.001). Two compounds (N-acetyl dipeptide and 4-substituted benzoic acid) clearly showed competitive mechanisms, whereas ascorbic acid reflected non-competitive kinetics. The five small organic molecules constitute five different scaffolds with moderate NA affinities. They are proposed as lead compounds for developing new NA inhibitors which are not analogous to sialic acid.

Topology-Matching Design of an Influenza-Neutralizing Spiky Nanoparticle-Based Inhibitor with a Dual Mode of Action.

In this study, we demonstrate the concept of "topology-matching design" for virus inhibitors. With the current knowledge of influenza A virus (IAV), we designed a nanoparticle-based inhibitor (nano-inhibitor) that has a matched nanotopology to IAV virions and shows heteromultivalent inhibitory effects on hemagglutinin and neuraminidase. The synthesized nano-inhibitor can neutralize the viral particle extracellularly and block its attachment and entry to the host cells. The virus replication was significantly reduced by 6 orders of magnitude in the presence of the reverse designed nano-inhibitors. Even when used 24 hours after the infection, more than 99.999 % inhibition is still achieved, which indicates such a nano-inhibitor might be a potent antiviral for the treatment of influenza infection.

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease.

BACKGROUND: The COVID-19 has now been declared a global pandemic by the World Health Organization. There is an emergent need to search for possible medications. METHOD: Utilization of the available sequence information, homology modeling, and in slico docking a number of available medications might prove to be effective in inhibiting the SARS-CoV-2 two main drug targets, the spike glycoprotein, and the 3CL protease. RESULTS: Several compounds were determined from the in silico docking models that might prove to be effective inhibitors for SARS-CoV-2. Several antiviral medications: Zanamivir, Indinavir, Saquinavir, and Remdesivir show potential as and 3CLPRO main proteinase inhibitors and as a treatment for COVID-19. CONCLUSION: Zanamivir, Indinavir, Saquinavir, and Remdesivir are among the exciting hits on the 3CLPRO main proteinase. It is also exciting to uncover that Flavin Adenine Dinucleotide (FAD) Adeflavin, B2 deficiency medicine, and Coenzyme A, a coenzyme, may also be potentially used for the treatment of SARS-CoV-2 infections. The use of these off-label medications may be beneficial in the treatment of the COVID-19.

[The model of prediction of the inhibition of neuraminidases of influenza A and B based on a reduced set of energy terms]. / Model' predskazaniia ingibirovaniia neiraminidaz virusa grippa A i B na baze ogranichennogo nabora énergeticheskikh vkladov.

The overall model for prediction of IC50 values for inhibitors of neuraminidase influenza virus A and B has been created. It combines data about IC50 values of complexes of 40 variants of neuraminidases of influenza A (7 serotypes) and B and three known inhibitors (oseltamivir, zanamivir, peramivir). The model also uses only data of enthalpy contributions to the potential energy of inhibitor/protein and substrate (MUNANA)/protein complexes. The calculation procedures are ported to use software with support of GPU accelerators, that significant decrease the computation time. The corresponding correlation coefficient (R²) for pIC50 prediction was within 0.45-0.58, the SEM values of around 0.7 (the range of used pIC50 data set is from 4.55 to 10.22).

A Small-Molecule Compound Has Anti-influenza A Virus Activity by Acting as a ''PB2 Inhibitor".

With regular influenza epidemics and the prevalence of drug-resistant influenza virus strains, it is extremely crucial to develop effective and low-toxicity anti-influenza A virus drugs that act on conserved sites of novel targets. Here, we found a new anti-influenza virus compound, 1,3-dihydroxy-6-benzo[ c]chromene (D715-2441), from a library of 8026 small-molecule compounds by cell-based MTT assay and explored the underlying mechanisms. Our results revealed that D715-2441 possessed antiviral activities against multiple subtypes of influenza A viruses (IAVs) strains, including H1N1, H5N1, H7N9, H3N2, the clinical isolate 690 (H3), and oseltamivir-resistant strains with the H274Y NA mutation, and suppressed the early steps in the virus replication cycle. Further mechanistic studies indicated that D715-2441 clearly inhibited viral polymerase activity and directly influenced the location of the PB2 protein. Moreover, binding affinity analyses confirmed that D715-2441 bound specifically to the PB2cap protein. Further, protein sequence alignment and a computer-aided molecular docking indicated that highly conserved amino acid residues in the cap-binding pocket of PB2cap were possible binding sites for D715-2441, which indicates that D715-2441 might be employed as a cap-binding competitor. Moreover, the combination of D715-2441 and zanamivir possessed a remarkable synergistic antiviral effect, with an FICI value of 0.40. In conclusion, these results strongly suggest that D715-2441 has potential as a promising candidate against IAV infection. More importantly, our work offers novel options for the strategic development of PB2cap inhibitors of IAV.

Molecular dynamics simulations of viral neuraminidase inhibitors with the human neuraminidase enzymes: Insights into isoenzyme selectivity.

Inhibitors of viral neuraminidase enzymes have been previously developed as therapeutics. Humans can express multiple forms of neuraminidase enzymes (NEU1, NEU2, NEU3, NEU4) that share a similar active site and enzymatic mechanism with their viral counterparts. Using a panel of purified human neuraminidase enzymes, we tested the inhibitory activity of 2-deoxy-2,3-dehydro-N-acetylneuraminic acid (DANA), zanamivir, oseltamivir, and peramivir against each of the human isoenzymes. We find that, with the exceptions of DANA and zanamivir, these compounds show generally poor activity against the human neuraminidase enzymes. To provide insight into the interactions of viral inhibitors with human neuraminidases, we conducted molecular dynamics simulations using homology models based on coordinates reported for NEU2. Simulations revealed that an organized water is displaced by zanamivir in binding to NEU2 and NEU3 and confirmed the critical importance of engaging the binding pocket of the C7-C9 glycerol sidechain. Our results suggest that compounds designed to target the human neuraminidases should provide more selective tools for interrogating these enzymes. Furthermore, they emphasize a need for additional structural data to enable structure-based drug design in these systems.

Inhibition of sialidase activity and cellular invasion by the bacterial vaginosis pathogen Gardnerella vaginalis.

Bacterial vaginosis is a genital tract infection, thought to be caused by transformation of a lactobacillus-rich flora to a dysbiotic microbiota enriched in mixed anaerobes. The most prominent of these is Gardnerella vaginalis (GV), an anaerobic pathogen that produces sialidase enzyme to cleave terminal sialic acid residues from human glycans. Notably, high sialidase activity is associated with preterm birth and low birthweight. We explored the potential of the sialidase inhibitor Zanamavir against GV whole cell sialidase activity using methyl-umbelliferyl neuraminic acid (MU-NANA) cleavage assays, with Zanamavir causing a 30% reduction in whole cell GV sialidase activity (p < 0.05). Furthermore, cellular invasion assays using HeLa cervical epithelial cells, infected with GV, demonstrated that Zanamivir elicited a 50% reduction in cell association and invasion (p < 0.05). Our data thus highlight that pharmacological sialidase inhibitors are able to modify BV-associated sialidase activity and influence host-pathogen interactions and may represent novel therapeutic adjuncts.

Acylguanidine derivatives of zanamivir and oseltamivir: Potential orally available prodrugs against influenza viruses.

Zanamivir (ZA) and guanidino-oseltamivir carboxylic acid (GOC) are very potent inhibitors against influenza neuraminidase (NA). The guanidinium moiety plays an important role in NA binding; however, its polar cationic nature also hinders the use of ZA and GOC from oral administration. In this study, we investigated the use of ZA and GOC acylguanidine derivatives as possible orally available prodrugs. The acylguanidine derivatives were prepared by coupling with either n-octanoic acid or (S)-naproxen. The lipophilic acyl substituents were verified to improve cell permeability, and may also improve the bioavailability of acylguanidine compounds. In comparison, the acylguanidines bearing linear octanoyl chain showed better NA inhibitory activity and higher hydrolysis rate than the corresponding derivatives having bulky branched naproxen moiety. Our molecular docking experiments revealed that the straight octanoyl chain could extend to the 150-cavity and 430-cavity of NA to gain extra hydrophobic interactions. Mice receiving the ZA octanoylguanidine derivative survived from influenza infection better than those treated with ZA, whereas the GOC octanoylguanidine derivative could be orally administrated to treat mice with efficacy equal to oseltamivir. Our present study demonstrates that incorporation of appropriate lipophilic acyl substituents to the polar guanidine group of ZA and GOC is a feasible approach to develop oral drugs for influenza therapy.

In silico prediction of drug resistance due to S247R mutation of Influenza H1N1 neuraminidase protein.

We present here in silico studies on antiviral drug resistance due to a novel mutation of influenza A/H1N1 neuraminidase (NA) protein. Influenza A/H1N1 virus was responsible for a recent pandemic and is currently circulating among the seasonal influenza strains. M2 and NA are the two major viral proteins related to pathogenesis in humans and have been targeted for drug designing. Among them, NA is preferred because the ligand-binding site of NA is highly conserved between different strains of influenza virus. Different mutations of the NA active site residues leading to drug resistance or susceptibility of the virus were studied earlier. We report here a novel mutation (S247R) in the NA protein that was sequenced earlier from the nasopharyngeal swab from Sri Lanka and Thailand in the year 2009 and 2011, respectively. Another mutation (S247N) was already known to confer resistance to oseltamivir. We did a comparative study of these two mutations vis-a-vis the drug-sensitive wild type NA to understand the mechanism of drug resistance of S247N and to predict the probability of the novel S247R mutation to become resistant to the currently available drugs, oseltamivir and zanamivir. We performed molecular docking- and molecular dynamics-based analysis of both the mutant proteins and showed that mutation of S247R affects drug binding to the protein by positional displacement due to altered active site cavity architecture, which in turn reduces the affinity of the drug molecules to the NA active site. Our analysis shows that S247R may have high probability of being resistant.

Different Inhibitory Potencies of Oseltamivir Carboxylate, Zanamivir, and Several Tannins on Bacterial and Viral Neuraminidases as Assessed in a Cell-Free Fluorescence-Based Enzyme Inhibition Assay.

Neuraminidaseis a key enzyme in the life cycle of influenza viruses and is present in some bacterial pathogens. We here assess the inhibitory potency of plant tannins versus clinically used inhibitors on both a viral and a bacterial model neuraminidase by applying the 2'-(4-methylumbelliferyl)-α-d-N-acetylneuraminic acid (MUNANA)-based activity assay. A range of flavan-3-ols, ellagitannins and chemically defined proanthocyanidin fractions was evaluated in comparison to oseltamivir carboxylate and zanamivir for their inhibitory activities against viral influenza A (H1N1) and bacterial Vibrio cholerae neuraminidase (VCNA). Compared to the positive controls, all tested polyphenols displayed a weak inhibition of the viral enzyme but similar or even higher potency on the bacterial neuraminidase. Structure-activity relationship analyses revealed the presence of galloyl groups and the hydroxylation pattern of the flavan skeleton to be crucial for inhibitory activity. The combination of zanamivir and EPs® 7630 (root extract of Pelargonium sidoides) showed synergistic inhibitory effects on the bacterial neuraminidase. Co-crystal structures of VCNA with oseltamivir carboxylate and zanamivir provided insight into bacterial versus viral enzyme-inhibitor interactions. The current data clearly indicate that inhibitor potency strongly depends on the biological origin of the enzyme and that results are not readily transferable. The therapeutic relevance of our findings is briefly discussed.

Multiscale Simulation of Receptor-Drug Association Kinetics: Application to Neuraminidase Inhibitors.

A detailed understanding of the drug-receptor association process is of fundamental importance for drug design. Due to the long time scales of typical binding kinetics, the atomistic simulation of the ligand traveling from bulk solution into the binding site is still computationally challenging. In this work, we apply a multiscale approach of combined Molecular Dynamics (MD) and Brownian Dynamics (BD) simulations to investigate association pathway ensembles for the two prominent H1N1 neuraminidase inhibitors oseltamivir and zanamivir. Including knowledge of the approximate binding site location allows for the selective confinement of detailed but expensive MD simulations and application of less demanding BD simulations for the diffusion controlled part of the association pathway. We evaluate a binding criterion based on the residence time of the inhibitor in the binding pocket and compare it to geometric criteria that require prior knowledge about the binding mechanism. The method ranks the association rates of both inhibitors in qualitative agreement with experiment and yields reasonable absolute values depending, however, on the reaction criteria. The simulated association pathway ensembles reveal that, first, ligands are oriented in the electrostatic field of the receptor. Subsequently, a salt bridge is formed between the inhibitor's carboxyl group and neuraminidase residue Arg368, followed by adopting the native binding mode. Unexpectedly, despite oseltamivir's higher overall association rate, the rate into the intermediate salt-bridge state was found to be higher for zanamivir. The present methodology is intrinsically parallelizable and, although computationally demanding, allows systematic binding rate calculation on selected sets of potential drug molecules.

Detection of Enzyme Inhibitors in Crude Natural Extracts Using Droplet-Based Microfluidics Coupled to HPLC.

Natural product screening for new bioactive compounds can greatly benefit from low reagents consumption and high throughput capacity of droplet-based microfluidic systems. However, the creation of large droplet libraries in which each droplet carries a different compound is a challenging task. A possible solution is to use an HPLC coupled to a droplet generating microfluidic device to sequentially encapsulate the eluting compounds. In this work we demonstrate the feasibility of carrying out enzyme inhibiting assays inside nanoliter droplets with the different components of a natural crude extract after being separated by a coupled HPLC column. In the droplet formation zone, the eluted components are mixed with an enzyme and a fluorogenic substrate that permits to follow the enzymatic reaction in the presence of each chromatographic peak and identify those inhibiting the enzyme activity. Using a fractal shape channel design and automated image analysis, we were able to identify inhibitors of Clostridium perfringens neuraminidase present in a root extract of the Pelargonium sidoides plant. This work demonstrates the feasibility of bioprofiling a natural crude extract after being separated in HPLC using microfluidic droplets online and represents an advance in the miniaturization of natural products screening.

Evaluation of Abbreviated Impactor Measurements (AIM) and Efficient Data Analysis (EDA) for Dry Powder Inhalers (DPIs) Against the Full-Resolution Next Generation Impactor (NGI).

The full-resolution next generation impactor (NGI) and three abbreviated impactor systems were used to obtain the apparent aerodynamic particle size distribution (APSD) and other quality measures for marketed dry powder inhalers (DPIs) using the compendial method and efficient data analysis (EDA). APSD for the active pharmaceutical ingredient (API) in Spiriva® Handihaler®, Foradil® Aerolizer®, and Relenza® Diskhaler® was obtained using a full-resolution NGI at 39, 60, and 90 L/min, respectively. Two reduced NGI (rNGI) configurations, the filter-only configuration (rNGI-f) and the modified-cup configuration (rNGI-mc), and the fast-screening impactor (FSI) with appropriate inserts to provide a 5-µm cut size were evaluated. The fine particle dose (FPD) obtained using the FSI for Spiriva was statistically similar to that obtained using the full NGI. However, the FPD for both Foradil and Relenza obtained using the FSI was significantly different from that obtained using the full NGI. Despite this, no significant differences were observed for the fine particle fraction (FPF) obtained using the FSI relative to that obtained from the full NGI for any of the DPIs. The use of abbreviated impactor systems appears promising with good agreement observed with the full-resolution NGI, except for small differences observed for the rNGI-mc configuration. These small differences may be product- and/or flow rate-specific, and further evaluation will be required to resolve these differences.

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.

Hydration of ligands of influenza virus neuraminidase studied by the fragment molecular orbital method.

The fragment molecular orbital (FMO) method was applied to quantum chemical calculations of neuramic acid, the natural substrate of the influenza virus neuraminidase, and two of its competitive inhibitors, Oseltamivir (Tamiful(®)) and Zanamivir (Relenza(®)), to investigate their hydrated structures and energetics. Each of the three ligands was immersed in an explicit water solvent, geometry-optimized by classical MM and QM/MM methods, and subjected to FMO calculations with 2-, 3-, and 4-body corrections under several fragmentation options. The important findings were that QM/MM optimization was preferable to obtain reliable hydrated structures of the ligands, that the 3-body correction was important for quantitative evaluation of the solvation energy, and that the dehydration effect was most remarkable near the hydrophobic sections of the ligands. In addition, the hydration energy calculated by the explicit solvent was compared with the hydration free energy calculated by the implicit solvent via the Poisson-Boltzmann equation, and the two showed a fairly good correlation. These findings will serve as useful information for rapid drug design.