Design of multi-binding-site inhibitors, ligand efficiency, and consensus screening of avian influenza H5N1 wild-type neuraminidase and of the oseltamivir-resistant H274Y variant

The binding sites of wild-type avian influenza A H5N1 neuraminidase, as well as those of the Tamiflu (oseltamivir)-resistant H274Y variant, were explored computationally to design inhibitors that target simultaneously several adjacent binding sites of the open conformation of the virus protein. The compounds with the best computed free energies of binding, in agreement by two docking methods, consensus scoring, and ligand efficiency values, suggest that mimicking a polysaccharide, beta-lactam, and other structures, including known drugs, could be routes for multibinding site inhibitor design. This new virtual screening method based on consensus scoring and ligand efficiency indices is introduced, which allows the combination of pharmacodynamic and pharmacokinetic properties into unique measures.     

Virtual Screening and Molecular Dynamics Simulation Study of Influenza Polymerase PB2 Inhibitors

Influenza A virus is the main cause of worldwide epidemics and annual influenza outbreaks in humans. In this study, a virtual screen was performed to identify compounds that interact with the PB2 cap-binding domain (CBD) of influenza A polymerase. A virtual screening workflow based on Glide docking was used to screen an internal database containing 8417 molecules, and then the output compounds were selected based on solubility, absorbance, and structural fingerprints. Of the 16 compounds selected for biological evaluation, six compounds were identified that rescued cells from H1N1 virus-mediated death at non-cytotoxic concentrations, with EC50 values ranging from 2.5–55.43 μM, and that could bind to the PB2 CBD of H1N1, with Kd values ranging from 0.081–1.53 μM. Molecular dynamics (MD) simulations of the docking complexes of our active compounds revealed that each compound had its own binding characteristics that differed from those of VX-787. Our active compounds have novel structures and unique binding modes with PB2 proteins, and are suitable to serve as lead compounds for the development of PB2 inhibitors. An analysis of the MD simulation also helped us to identify the dominant amino acid residues that play a key role in binding the ligand to PB2, suggesting that we should focus on increasing and enhancing the interaction between inhibitors and these major amino acids during lead compound optimization to obtain more active PB2 inhibitors.     

Binding interaction analysis of the active site and its inhibitors for neuraminidase (N1 subtype) of human influenza virus by the integration of molecular docking, FMO calculation and 3D-QSAR CoMFA modeling

Recently, the worldwide spread of A/H5N1 avian influenza with high virulence has highlighted the potential threat of human influenza pandemic. Tamiflu and Relenza are currently the only two anti-influenza drugs targeting the neuraminidase (NA) enzyme of human influenza virus. Reports of the emergence of drug resistance further make the development of new potent anti-influenza inhibitors a priority. The X-ray crystallographic study of A/H5N1 avian influenza NA subtypes (Russell, R. J. Nature 2006, 443, 45-49) has demonstrated that there exist two genetically distinct groups, group-1 (N1, N4, N5 and N8) and group-2 (N2, N3, N6, N7 and N9), whose conformations are substantially different. The detailed comparison of their active sites has established, heretofore, the most accurate and solid molecular basis of structure and mechanism for the development of new anti-influenza drugs. In the present study, a three-dimensional structure of N1 subtype of human influenza type A virus (N1hA) has been generated by homology modeling using the X-ray crystallographic structure of N1 subtype of avian influenza virus (N1aA) as the template. Binding interaction analysis between the active site and its inhibitors has been performed by combining ab initio fragment molecular orbital (FMO) calculations and three-dimensional quantitative structure-activity relationship with comparative molecular field analysis (3D-QSAR CoMFA) modeling. Integrated with docking-based 3D-QSAR CoMFA modeling, molecular surface property (electrostatic and steric) mapping and FMO pair interaction analysis, a set of new receptor-ligand binding models and bioaffinity predictive models for rational design and virtual screening of more potent inhibitors of N1hA are established. In addition, the flexibility of the loop-150 of N1hA and N1aA has been examined by a series of molecular dynamics simulations.     

单管高灵敏度等温扩增技术快速检测甲型H1N1流感病毒

建立了单管逆转录环介导等温扩增法(RT-LAMP)快速检测甲型H1N1流感病毒的方法.针对甲型H1N1流感病毒的M基因和HA基因的保守区,设计了两组特异性引物,分别用于筛选甲型流感病毒及鉴定甲型H1N1流感病毒.对反应体系中的关键因素进行优化,反应结果可直接通过浊度或者SYBR Green Ⅰ荧光进行判定.本方法最低可...     

Identification of Novel Influenza Polymerase PB2 Inhibitors Using a Cascade Docking Virtual Screening Approach

To discover novel inhibitors that target the influenza polymerase basic protein 2 (PB2) cap-binding domain (CBD), commercial ChemBridge compound libraries containing 384,796 compounds were screened using a cascade docking of LibDock–LigandFit–GOLD, and 60 compounds were selected for testing with cytopathic effect (CPE) inhibition assays and surface plasmon resonance (SPR) assay. Ten compounds were identified to rescue cells from H1N1 virus-mediated death at non-cytotoxic concentrations with EC₅₀ values ranging from 0.30 to 67.65 μM and could bind to the PB2 CBD of H1N1 with Kd values ranging from 0.21 to 6.77 μM. Among these, four compounds (11D4, 12C5, 21A5, and 21B1) showed inhibition of a broad spectrum of influenza virus strains, including oseltamivir-resistant ones, the PR/8-R292K mutant (H1N1, recombinant oseltamivir-resistant strain), the PR/8-I38T mutant (H1N1, recombinant baloxavir-resistant strain), and the influenza B/Lee/40 virus strain. These compounds have novel chemical scaffolds and relatively small molecular weights and are suitable for optimization as lead compounds. Based on sequence and structure comparisons of PB2 CBDs of various influenza virus subtypes, we propose that the Phe323/Gln325, Asn429/Ser431, and Arg355/Gly357 mutations, particularly the Arg355/Gly357 mutation, have a marked impact on the selectivities of PB2 CBD-targeted inhibitors of influenza A and influenza B.     

Differential receptor binding affinities of influenza hemagglutinins on glycan arrays

A library of 27 sialosides, including seventeen 2,3-linked and ten 2,6-linked glycans, has been prepared to construct a glycan array and used to profile the binding specificity of different influenza hemagglutinins (HA) subtypes, especially from the 2009 swine-originated H1N1 and seasonal influenza viruses. It was found that the HAs from the 2009 H1N1 and the seasonal Brisbane strain share similar binding profiles yet different binding affinities toward various α2,6 sialosides. Analysis of the binding profiles of different HA subtypes indicate that a minimum set of 5 oligosaccharides can be used to differentiate influenza H1, H3, H5, H7, and H9 subtypes. In addition, the glycan array was used to profile the binding pattern of different influenza viruses. It was found that most binding patterns of viruses and HA proteins are similar and that glycosylation at Asn27 is essential for receptor binding.     

Pharmacophore Mode lingand Virtual Screening to Design the Potential Influenza Virus Endonuclease Inhibitors

Influenza virus endonuclease is an attractive target for antiviral therapy in the treatment of influenza infection. The purpos e of this study is to design a novel antiviral agent with improved biological activities against the influenza virus endonuclease. In this study, chemical feature‐based 3D pharmacophore models were developed from 41 known influenza virus endonuclease inhibitors. The best quantitative pharmacohore model (Hypo 1), which consists of two hydrogen‐bond acceptors and two hydrophobic features, yields the highest correlation coefficient (R = 0.886). Hypo 1 was further validated by the cross validation method and the test set compounds. The application of this model for predicting the activities of 11 known influenza virus endonuclease inhibitors in the test set shows great success. The correlation coefficient of 0.942 and a cross validation of 95;% confidence level prove that this model is reliable in identifying structurally diverse compounds for influenza virus endonuclease inhibition. The most active compound (compound 1) from the training set was docked into the active site of the influenza virus endonuclease as an additional verification that the pharmacophore model is accurate. The docked conformation showed important hydrogen bond interactions between the compound and two amino acids, Lys 134 and Lys 137. After validation, this model was used to screen the NCI chemical database to identify new influenza virus endonuclease inhibitors. Our study shows that the to pranking compound out of the 10 newly identified compounds using fit value ranking has an estimated activity of 0.049 μM. These newly identified lead compounds can be further experimentally validated using in vitro techniques.     

Molecular dynamics simulation directed rational design of inhibitors targeting drug-resistant mutants of influenza A virus M2

Influenza A virus M2 (A/M2) forms a homotetrameric proton selective channel in the viral membrane. It has been the drug target of antiviral drugs such as amantadine and rimantadine. However, most of the current virulent influenza A viruses carry drug-resistant mutations alongside the drug binding site, such as S31N, V27A, and L26F, etc., each of which might be dominant in a given flu season. Among these mutations, the V27A mutation was prevalent among transmissible viruses under drug selection pressure. Until now, V27A has not been successfully targeted by small molecule inhibitors, despite years of extensive medicinal chemistry research efforts and high throughput screening. Guided by molecular dynamics (MD) simulation of drug binding and the influence of drug binding on the dynamics of A/M2 from earlier experimental studies, we designed a series of potent spirane amine inhibitors targeting not only WT, but also both A/M2-27A and L26F mutants with IC(50)s similar to that seen for amantadine's inhibition of the WT channel. The potencies of these inhibitors were further demonstrated in experimental binding and plaque reduction assays. These results demonstrate the power of MD simulations to probe the mechanism of drug binding as well as the ability to guide design of inhibitors of targets that had previously appeared to be undruggable.     

Top leads for swine influenza A/H1N1 virus revealed by steered molecular dynamics approach

Since March 2009, the rapid spread of infection during the recent A/H1N1 swine flu pandemic has raised concerns of a far more dangerous outcome should this virus become resistant to current drug therapies. Currently oseltamivir (tamiflu) is intensively used for the treatment of influenza and is reported effective for 2009 A/H1N1 virus. However, as this virus is evolving fast, some drug-resistant strains are emerging. Therefore, it is critical to seek alternative treatments and identify roots of the drug resistance. In this paper, we use the steered molecular dynamics (SMD) approach to estimate the binding affinity of ligands to the glycoprotein neuraminidase. Our idea is based on the hypothesis that the larger is the force needed to unbind a ligand from a receptor the higher its binding affinity. Using all-atom models with Gromos force field 43a1 and explicit water, we have studied the binding ability of 32 ligands to glycoprotein neuraminidase from swine flu virus A/H1N1. The electrostatic interaction is shown to play a more important role in binding affinity than the van der Waals one. We have found that four ligands 141562, 5069, 46080, and 117079 from the NSC set are the most promising candidates to cope with this virus, while peramivir, oseltamivir, and zanamivir are ranked 8, 11, and 20. The observation that these four ligands are better than existing commercial drugs has been also confirmed by our results on the binding free energies obtained by the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. Our prediction may be useful for the therapeutic application.     

Free energy calculations on the two drug binding sites in the M2 proton channel

Two alternative binding sites of adamantane-type drugs in the influenza A M2 channel have been suggested, one with the drug binding inside the channel pore and the other with four drug molecule S-binding to the C-terminal surface of the transmembrane domain. Recent computational and experimental studies have suggested that the pore binding site is more energetically favorable but the external surface binding site may also exist. Nonetheless, which drug binding site leads to channel inhibition in vivo and how drug-resistant mutations affect these sites are not completely understood. We applied molecular dynamics simulations and potential of mean force calculations to examine the structures and the free energies associated with these putative drug binding sites in an M2-lipid bilayer system. We found that, at biological pH (~7.4), the pore binding site is more thermodynamically favorable than the surface binding site by ~7 kcal/mol and, hence, would lead to more stable drug binding and channel inhibition. This result is in excellent agreement with several recent studies. More importantly, a novel finding of ours is that binding to the channel pore requires overcoming a much higher energy barrier of ~10 kcal/mol than binding to the C-terminal channel surface, indicating that the latter site is more kinetically favorable. Our study is the first computational work that provides both kinetic and thermodynamic energy information on these drug binding sites. Our results provide a theoretical framework to interpret and reconcile existing and often conflicting results regarding these two binding sites, thus helping to expand our understanding of M2-drug binding, and may help guide the design and screening of novel drugs to combat the virus.     

Photoinduced electron transfer and fluorophore motion as a probe of the conformational dynamics of membrane proteins: application to the influenza a M2 proton channel

The structure and function of the influenza A M2 proton channel have been the subject of intensive investigations in recent years because of their critical role in the life cycle of the influenza virus. Using a truncated version of the M2 proton channel (i.e., M2TM) as a model, here we show that fluctuations in the fluorescence intensity of a dye reporter that arise from both fluorescence quenching via the mechanism of photoinduced electron transfer (PET) by an adjacent tryptophan (Trp) residue and local motions of the dye molecule can be used to probe the conformational dynamics of membrane proteins. Specifically, we find that the dynamics of the conformational transition between the N-terminal open and C-terminal open states of the M2TM channel occur on a timescale of about 500 μs and that the binding of either amantadine or rimantadine does not inhibit the pH-induced structural equilibrium of the channel. These results are consistent with the direct occluding mechanism of inhibition which suggests that the antiviral drugs act by sterically occluding the channel pore.     

Redesign of protein domains using one-bead-one-compound combinatorial chemistry

A novel combinatorial strategy for the redesign of proteins based on the strength and specificity of intra- and interprotein interactions is described. The strategy has been used to redesign the hydrophobic core of the B domain of protein A. Using one-bead-one-compound combinatorial chemistry, 300 analogues of the C-terminal alpha-helix of the B domain, H3x, have been synthesized using a biocompatible resin and the HMFS linker, allowing the screening to occur while the peptides were bound to the resin. The screening was based on the ability of the H3x analogues to interact with the N-terminal helices of the B domain, H1-H2, and retain the native B domain activity, that is binding to IgG. Eight different analogues containing some nonconservative mutations were obtained from the library, the two most frequent of which, H3P1 and H3P2, were studied in detail. CD analysis revealed that the active analogues interact with H1-H2. To validate the redesign strategy the covalent modified domains H1-H2-H3P1 and H1-H2-H3P2 were synthesized and characterized. CD and NMR analysis revealed that they had a unique, stable, and well-defined three-dimensional structure similar to that for the wild-type B domain. This combinatorial strategy allows us to select for redesigned proteins with the desired activity or the desired physicochemical properties provided the right screening test is used. Furthermore, it is rich in potential for the chemical modification of proteins overcoming the drawbacks associated with the total synthesis of large protein domains.     

Catalytic and biochemical features of a monoclonal antibody heavy chain, JN1-2, raised against a synthetic peptide with a hemagglutinin molecule of influenza virus

It has long been an important issue to produce a catalytic antibody that possesses the ability to lose the infectivity of a bacteria or virus. The monoclonal antibody JN1-2 was generated using a synthetic peptide (TGLRNGITNKVNSVIEKAA) conjugated with human IgG. The peptide sequence includes the conserved region of the hemagglutinin molecule (HA(1) and HA(2) domains), which locates on the envelope of the influenza virus and plays an important role in influenza A virus infection. The monoclonal antibody specifically reacted with the HA2 domain, not only of H2 but also of an H1 strain of the H1N1 subtype (H1 strain). The heavy chain (JN1-2-H) isolated from the parent antibody showed catalytic activity cleaving the above antigenic peptide with very high turnover (kcat = 26 min(-1)), and it could slowly degrade the recombinant HA(2) domain by the catalytic function. Interestingly, the heavy chain exhibited the ability to reduce the infectivity of type A H1N1 but not type B, indicating specificity to type A. This characteristic monoclonal catalytic antibody heavy chain could suppress the infection of the influenza virus in vitro assays.     

Structure-based virtual screening of influenza virus RNA polymerase inhibitors from natural compounds: Molecular dynamics simulation and MM-GBSA calculation

The resistances of matrix protein 2 (M2) protein inhibitors and neuraminidase inhibitors for influenza virus have attracted much attention and there is an urgent need for new drug. The antiviral drugs that selectively act on RNA polymerase are less prone to resistance and possess fewer side effects on the patient. Therefore, there is increased interest in screening compounds that can inhibit influenza virus RNA polymerase. Three natural compounds were found by using molecular docking-based virtual screening, which could bind tightly within the polymerase acidic protein-polymerase basic protein 1 (PA-PB1) subunit of influenza virus polymerase. Firstly, their drug likeness properties were evaluated, which showed that the hepatotoxicity values of all the three compounds indicating they had less or no hepatotoxicity, and did not have the plasma protein biding (PPB) ability, the three compounds needed to be modified in some aspects, like bulky molecular size. The stability of the complexes of PA-hits was validated through molecular dynamics (MD) simulation, revealing compound 2 could form more stable complex with PA subunit. The torsional conformations of each rotatable bond of the ligands in PA subunit were also monitored, to investigate variation in the ligand properties during the simulation, compound 3 had fewer rotatable bonds, indicating that the molecule had stronger rigidity. The bar charts of protein–ligand contacts and contacts over the course of trajectory showed that four key residues (Glu623, Lys643, Asn703 and Trp706) of PA subunit that participated in hydrogen-bond, water bridge and hydrophobic interactions with the hit compounds. Finally, the binding free energy and contributed energies were calculated by using MM-GBSA method. Out of the three compounds, compound 1 showed the lowest total binding free energy. Among all the interactions, the contribution of the covalent binding and the van der Waals energy were more than other items, compound 1 formed more stable hydrogen bonds with the residues of PA subunit binding pocket. This study smoothed the path for the development of novel lead compounds with improved binding properties, high drug likeness, and low toxicity to humans for the treatment of influenza, which provided a good basis for further research on novel and effective influenza virus PA-PB1 interaction inhibitors.     

Synthesis and Antiviral Evaluation of Nucleoside Analogues Bearing One Pyrimidine Moiety and Two D-Ribofuranosyl Residues

A series of 1,2,3-triazolyl nucleoside analogues in which 1,2,3-triazol-4-yl-β-d-ribofuranosyl fragments are attached via polymethylene linkers to both nitrogen atoms of the heterocycle moiety (uracil, 6-methyluracil, thymine, quinazoline-2,4-dione, alloxazine) or to the C-5 and N-3 atoms of the 6-methyluracil moiety was synthesized. All compounds synthesized were evaluated for antiviral activity against influenza virus A/PR/8/34/(H1N1) and coxsackievirus B3. Antiviral assays revealed three compounds, 2i, 5i, 11c, which showed moderate activity against influenza virus A H1N1 with IC₅₀ values of 57.5 µM, 24.3 µM, and 29.2 µM, respectively. In the first two nucleoside analogues, 1,2,3-triazol-4-yl-β-d-ribofuranosyl fragments are attached via butylene linkers to N-1 and N-3 atoms of the heterocycle moiety (6-methyluracil and alloxazine, respectively). In nucleoside analogue 11c, two 1,2,3-triazol-4-yl-2′,3′,5′-tri-O-acetyl-β-d-ribofuranose fragments are attached via propylene linkers to the C-5 and N-3 atoms of the 6-methyluracil moiety. Almost all synthesized 1,2,3-triazolyl nucleoside analogues showed no antiviral activity against the coxsackie B3 virus. Two exceptions are 1,2,3-triazolyl nucleoside analogs 2f and 5f, in which 1,2,3-triazol-4-yl-2′,3′,5′-tri-O-acetyl-β-d-ribofuranose fragments are attached to the C-5 and N-3 atoms of the heterocycle moiety (6-methyluracil and alloxazine respectively). These compounds exhibited high antiviral potency against the coxsackie B3 virus with IC₅₀ values of 12.4 and 11.3 µM, respectively, although both were inactive against influenza virus A H1N1. According to theoretical calculations, the antiviral activity of the 1,2,3-triazolyl nucleoside analogues 2i, 5i, and 11c against the H1N1 (A/PR/8/34) influenza virus can be explained by their influence on the functioning of the polymerase acidic protein (PA) of RNA-dependent RNA polymerase (RdRp). As to the antiviral activity of nucleoside analogs 2f and 5f against coxsackievirus B3, it can be explained by their interaction with the coat proteins VP1 and VP2.     

Inhibitory effect and possible mechanism of action of patchouli alcohol against influenza A (H2N2) virus

In the present study, the anti-influenza A (H2N2) virus activity of patchouli alcohol was studied in vitro, in vivo and in silico. The CC?? of patchouli alcohol was above 20 μM. Patchouli alcohol could inhibit influenza virus with an IC?? of 4.03 ± 0.23 μM. MTT assay showed that the inhibition by patchouli alcohol appears strongly after penetration of the virus into the cell. In the influenza mouse model, patchouli alcohol showed obvious protection against the viral infection at a dose of 5 mg/kg/day. Flexible docking and molecular dynamic simulations indicated that patchouli alcohol was bound to the neuraminidase protein of influenza virus, with an interaction energy of -40.38 kcal mol?1. The invariant key active-site residues Asp151, Arg152, Glu119, Glu276 and Tyr406 played important roles during the binding process. Based on spatial and energetic criteria, patchouli alcohol interfered with the NA functions. Results presented here suggest that patchouli alcohol possesses anti-influenza A (H2N2) virus properties, and therefore is a potential source of anti-influenza agents for the pharmaceutical industry.     

A practical synthesis of zanamivir phosphonate congeners with potent anti-influenza activity

Two phosphonate compounds 1a (4-amino-1-phosphono-DANA) and 1b (phosphono-zanamivir) are synthesized and shown more potent than zanamivir against the neuraminidases of avian and human influenza viruses, including the oseltamivir-resistant strains. For the first time, the practical synthesis of these phosphonate compounds is realized by conversion of sialic acid to peracetylated phosphono-DANA diethyl ester (5) as a key intermediate in three steps by a novel approach. In comparison with zanamivir, the high affinity of 1a and 1b can be partly attributable to the strong electrostatic interactions of their phosphonate groups with the three arginine residues (Arg118, Arg292, and Arg371) in the active site of neuraminidases. These phosphonates are nontoxic to the human 293T cells; they protect cells from influenza virus infection with EC(50) values in low-nanomolar range, including the wild-type WSN (H1N1), the 2009 pandemic (H1N1), the oseltamivir-resistant H274Y (H1N1), RG14 (H5N1), and Udorn (H3N2) influenza strains.     

Novel algorithm for phylogenetic analysis of proteins: application to analysis of the evolution of H5N1 influenza viruses

The highly pathogenic avian influenza virus (HPAIV) A subtype H5N1 is causing threat to human health over the years. Phylogenetic analysis is an important tool for analyzing the evolution of influenza. A novel phylogenetic algorithm based on a new protein distance measure derived from the informational spectrum method (ISM) has been presented. The new phylogenetic approach allows assessment of functional evolution of protein sequences. The new ISM-based phylogenetic approach has been found to overcome some drawbacks of other phylogenetic approaches, particularly concerning sensitivity to a single mutation, deletion and the position of the mutation. The ISM-based approach applied to hemagglutinin subunit 1 protein (HA1) of HPAIV A subtype H5N1 viruses in Egypt between 2006 and 2011, revealed clear clustering in two groups, with one growing group of H5N1 viruses after 2009 with increased number of human infections with H5N1. Four group-specific mutations are identified which are important for increased human tropism and the pandemic potential.     

Pharmacophore modelling,virtual screening and molecular docking studies on PLD1 inhibitors

Lipid metabolism plays a significant role in influenza virus replication and subsequent infection. The regulatory mechanism governing lipid metabolism and viral replication is not properly understood to date, but both Phospholipase D (PLD1 and PLD2) activities are stimulated in viral infection. In vitro studies indicate that chemical inhibition of PLD1 delays viral entry and reduction of viral loads. The current study reports a three-dimensional pharmacophore model based on 35 known PLD1 inhibitors. A sub-set of 25 compounds was selected as the training set and the remaining 10 compounds were kept in the test set. One hundred and twelve pharmacophore models were generated; a six-featured pharmacophore model (AADDHR.57) with survival score (2.69) produced a statistically significant three-dimensional quantitative structure–activity relationship model with r² = 0.97 (internal training set), r² = 0.71 (internal test set) and Q² = 0.64. The predictive power of the pharmacophore model was validated with an external test set (r² = 0.73) and a systematic virtual screening work-flow was employed showing an enrichment factor of 23.68 at the top 2% of the dataset (active and decoys). Finally, the model was used for screening of the filtered PubChem database to fetch molecules which can be proposed as potential PLD1 inhibitors for blocking influenza infection.