N-recognins UBR1 and UBR2 as central ER stress sensors in mammals.

In eukaryotes, a primary protein quality control (PQC) process involves the destruction of conformationally misfolded proteins through the ubiquitin-proteasome system. Because approximately one-third of eukaryotic proteomes fold and assemble within the endoplasmic reticulum (ER) before being sent to their destinations, the ER plays a crucial role in PQC. The specific functions and biochemical roles of several E3 ubiquitin ligases involved in ER-associated degradation in mammals, on the other hand, are mainly unknown. We identified 2 E3 ligases, ubiquitin protein ligase E3 component N-recognin 1 (UBR1) and ubiquitin protein ligase E3 component N-recognin 2 (UBR2), which are the key N-recognins in the N-degron pathway and participate in the ER stress response in mammalian cells by modulating their stability. Cells lacking UBR1 and UBR2 are hypersensitive to ER stress-induced apoptosis. Under normal circumstances, these proteins are polyubiquitinated through Lys48-specific linkages and are then degraded by the 26S proteasome. In contrast, when cells are subjected to ER stress, UBR1 and UBR2 exhibit greater stability, potentially as a cellular adaptive response to stressful conditions. Although the precise mechanisms underlying these findings require further investigation, our findings show that cytoplasmic UBR1 and UBR2 have anti-ER stress activities and contribute to global PQC in mammals. These data also reveal an additional level of complexity within the mammalian ER-associated degradation system, implicating potential involvement of the N-degron pathway.

A Multicenter Study of Clinicopathology and Immunohistochemical Distinction between Adult and Pediatric Large B-Cell Lymphoma.

Introduction: Pediatric DLBCL is considered a homogenous group and has superior outcomes compared to adults. This study investigated the clinical pathology and immunohistochemical distinction between adult and pediatric large B-cell lymphoma. Methods: A cross-sectional study of 314 NHLs with the morphology of diffuse pattern, large B-cell, and CD20 expression was investigated. Results: Of 314 cases, there were 6 cases of pleomorphic MCL (all in adults), 19 cases of Burkitt lymphoma (all in children), and 289 cases of DLBCL. Pediatric DLBCL had many striking differences: More frequency in extra-nodal sites; a higher proportion of centroblastic morphology; a predominance of GCB-type; a high proliferation rate; an infrequency of Bcl2 protein expression, and a lack of double-expresser lymphoma. Conclusions: Our study demonstrated the significant biological differences between adult and pediatric DLBCL.

iANP-EC: Identifying Anticancer Natural Products Using Ensemble Learning Incorporated with Evolutionary Computation.

Cancer is one of the most deadly diseases that annually kills millions of people worldwide. The investigation on anticancer medicines has never ceased to seek better and more adaptive agents with fewer side effects. Besides chemically synthetic anticancer compounds, natural products are scientifically proved as a highly potential alternative source for anticancer drug discovery. Along with experimental approaches being used to find anticancer drug candidates, computational approaches have been developed to virtually screen for potential anticancer compounds. In this study, we construct an ensemble computational framework, called iANP-EC, using machine learning approaches incorporated with evolutionary computation. Four learning algorithms (k-NN, SVM, RF, and XGB) and four molecular representation schemes are used to build a set of classifiers, among which the top-four best-performing classifiers are selected to form an ensemble classifier. Particle swarm optimization (PSO) is used to optimise the weights used to combined the four top classifiers. The models are developed by a set of curated 997 compounds which are collected from the NPACT and CancerHSP databases. The results show that iANP-EC is a stable, robust, and effective framework that achieves an AUC-ROC value of 0.9193 and an AUC-PR value of 0.8366. The comparative analysis of molecular substructures between natural anticarcinogens and nonanticarcinogens partially unveils several key substructures that drive anticancerous activities. We also deploy the proposed ensemble model as an online web server with a user-friendly interface to support the research community in identifying natural products with anticancer activities.

Prevalence and risk factors for human papillomavirus infection among female sex workers in Hanoi and Ho Chi Minh City, Viet Nam: a cross-sectional study.

Objective: Female sex workers (FSWs) are at high risk of human papillomavirus (HPV) infections and cervical cancer due to their high number of sexual partners. The objectives of this study were to determine the prevalence of HPV and identify risk factors for high-risk HPV infection among FSWs in Hanoi and Ho Chi Minh City (HCMC), Viet Nam. Methods: A cross-sectional study was conducted in Hanoi and HCMC between December 2017 and May 2018. We surveyed and screened 699 FSWs aged 318 years for HPV infection and abnormal cytology. A multivariable modified Cox regression model was used to determine risk factors for high-risk HPV infection. Results: The overall prevalence of any HPV, high-risk HPV and HPV-16/18 infection in the 699 FSWs was 26.3%, 17.6% and 4.0%, respectively, and were similar in both cities. Multiple infections were identified in 127 participants (69.0%). HPV-52 was the most prevalent (7%), followed by HPV-58 (6%). Abnormal cytology was detected in 91 participants (13.0%). FSWs who are divorced (adjusted prevalence ratio [aPR]: 1.96, 95% confidence interval [CI]: 1.01-3.81), widowed (aPR: 3.26, 95% CI: 1.49-7.12) or living alone (aPR: 1.85, 95% CI: 1.01-3.39) were associated with a higher prevalence of high-risk HPV infection. Discussion: Almost one in five FSWs in Viet Nam are infected with high-risk HPV. This highlights the importance of prevention strategies such as HPV vaccination and screening in this high-risk group.

Vietnam Association of Gastroenterology (VNAGE) consensus on the management of Helicobacter pylori infection.

Helicobacter pylori (H. pylori) infection is prevalent and has a rapidly increasing antibiotic resistance rate in Vietnam. Reinfection is quite common, and gastric carcinoma remains one of the most common malignancies, which is not uncommon to develop after successful eradication. The purpose of this consensus is to provide updated recommendations on the management of H. pylori infection in the country. The consensus panel consisted of 32 experts from 14 major universities and institutions in Vietnam who were invited to review the evidence and develop the statements using the Delphi method. The process followed the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) system. The consensus level was defined as ≥80% for agreement on the proposed statements. Due to the limited availability of high-quality local evidence, this consensus was also based on high-quality evidence from international studies, especially those conducted in other populations in the Asia-Pacific region. The panel finally reached a consensus on 27 statements after two voting rounds, which consisted of four sections (1) indications for testing and selection of diagnostic tests (2), treatment regimens, (3) post-treatment confirmation of H. pylori status, and (4) reinfection prevention methods and follow-up after eradication. Important issues that require further evidence include studies on third-line regimens, strategies to prevent H. pylori reinfection, and post-eradication follow-up for precancerous gastric lesions. We hope this consensus will help guide the current clinical practice in Vietnam and promote multicenter studies in the country and international collaborations.

Formation of Non-Nucleoplasmic Proteasome Foci during the Late Stage of Hyperosmotic Stress.

Cellular stress induces the formation of membraneless protein condensates in both the nucleus and cytoplasm. The nucleocytoplasmic transport of proteins mainly occurs through nuclear pore complexes (NPCs), whose efficiency is affected by various stress conditions. Here, we report that hyperosmotic stress compartmentalizes nuclear 26S proteasomes into dense nuclear foci, independent of signaling cascades. Most of the proteasome foci were detected between the condensed chromatin mass and inner nuclear membrane. The proteasome-positive puncta were not colocalized with other types of nuclear bodies and were reversibly dispersed when cells were returned to the isotonic medium. The structural integrity of 26S proteasomes in the nucleus was slightly affected under the hyperosmotic condition. We also found that these insulator-body-like proteasome foci were possibly formed through disrupted nucleus-to-cytosol transport, which was mediated by the sequestration of NPC components into osmostress-responding stress granules. These data suggest that phase separation in both the nucleus and cytosol may be a major cell survival mechanism during hyperosmotic stress conditions.

Development and Validation of a Scoring System to Predict Severe Acute Lower Gastrointestinal Bleeding in Vietnamese.

BACKGROUND/AIMS: The prevalence of acute lower gastrointestinal bleeding (ALGIB) has progressively increased worldwide but there are few studies in Asian populations. This study aimed to develop and validate a scoring system to predict severe ALGIB in Vietnamese. METHODS: Risk factors for severe ALGIB were identified by multiple logistic regression analysis using data from a retrospective cohort of 357 patients admitted to a tertiary hospital. These factors were weighted to develop the severe acute lower gastrointestinal bleeding (SALGIB) score to predict severe ALGIB. The performance of SALGIB was validated in a prospective cohort of 324 patients admitted to 6 other hospitals using area under the receiver operating characteristics curve (AUC) analysis. RESULTS: There were four factors at admission independently associated with severe ALGIB in the derivation cohort: heart rate ≥ 100/min, systolic blood pressure < 100 mmHg, hematocrit < 35%, and platelets ≤ 150 × 103/µL. The SALGIB score determined severe ALGIB with AUC values of 0.91 and 0.86 in the derivation and validation cohorts, respectively. A SALGIB score < 2 associated with low risk of severe ALGIB in both cohorts (3.7% and 1.2%; respectively). CONCLUSIONS: The SALGIB score has good performance in discriminating risk of severe ALGIB in Vietnamese.

Bioactive Compounds in Anti-Diabetic Plants: From Herbal Medicine to Modern Drug Discovery.

Natural products, including organisms (plants, animals, or microorganisms) have been shown to possess health benefits for animals and humans. According to the estimation of the World Health Organization, in developing countries, 80% of the population has still depended on traditional medicines or folk medicines which are mostly prepared from the plant for prevention or treatment diseases. Traditional medicine from plant extracts has proved to be more affordable, clinically effective and relatively less adverse effects than modern drugs. Literature shows that the attention on the application of phytochemical constituents of medicinal plants in the pharmaceutical industry has increased significantly. Plant-derived secondary metabolites are small molecules or macromolecules biosynthesized in plants including steroids, alkaloids, phenolic, lignans, carbohydrates and glycosides, etc. that possess a diversity of biological properties beneficial to humans, such as their antiallergic, anticancer, antimicrobial, anti-inflammatory, antidiabetic and antioxidant activities Diabetes mellitus is a chronic disease result of metabolic disorders in pancreas ß-cells that have hyperglycemia. Hyperglycemia can be caused by a deficiency of insulin production by pancreatic (Type 1 diabetes mellitus) or insufficiency of insulin production in the face of insulin resistance (Type 2 diabetes mellitus). The current medications of diabetes mellitus focus on controlling and lowering blood glucose levels in the vessel to a normal level. However, most modern drugs have many side effects causing some serious medical problems during a period of treating. Therefore, traditional medicines have been used for a long time and play an important role as alternative medicines. Moreover, during the past few years, some of the new bioactive drugs isolated from plants showed antidiabetic activity with more efficacy than oral hypoglycemic agents used in clinical therapy. Traditional medicine performed a good clinical practice and is showing a bright future in the therapy of diabetes mellitus. World Health Organization has pointed out this prevention of diabetes and its complications is not only a major challenge for the future, but essential if health for all is to be attained. Therefore, this paper briefly reviews active compounds, and pharmacological effects of some popular plants which have been widely used in diabetic treatment. Morphological data from V-herb database of each species was also included for plant identification.

Computational Determination of Potential Inhibitors of SARS-CoV-2 Main Protease.

The novel coronavirus (SARS-CoV-2) has infected several million people and caused thousands of deaths worldwide since December 2019. As the disease is spreading rapidly all over the world, it is urgent to find effective drugs to treat the virus. The main protease (Mpro) of SARS-CoV-2 is one of the potential drug targets. Therefore, in this context, we used rigorous computational methods, including molecular docking, fast pulling of ligand (FPL), and free energy perturbation (FEP), to investigate potential inhibitors of SARS-CoV-2 Mpro. We first tested our approach with three reported inhibitors of SARS-CoV-2 Mpro, and our computational results are in good agreement with the respective experimental data. Subsequently, we applied our approach on a database of ∼4600 natural compounds, as well as 8 available HIV-1 protease (PR) inhibitors and an aza-peptide epoxide. Molecular docking resulted in a short list of 35 natural compounds, which was subsequently refined using the FPL scheme. FPL simulations resulted in five potential inhibitors, including three natural compounds and two available HIV-1 PR inhibitors. Finally, FEP, the most accurate and precise method, was used to determine the absolute binding free energy of these five compounds. FEP results indicate that two natural compounds, cannabisin A and isoacteoside, and an HIV-1 PR inhibitor, darunavir, exhibit a large binding free energy to SARS-CoV-2 Mpro, which is larger than that of 13b, the most reliable SARS-CoV-2 Mpro inhibitor recently reported. The binding free energy largely arises from van der Waals interaction. We also found that Glu166 forms H-bonds to all of the inhibitors. Replacing Glu166 by an alanine residue leads to ∼2.0 kcal/mol decreases in the affinity of darunavir to SARS-CoV-2 Mpro. Our results could contribute to the development of potential drugs inhibiting SARS-CoV-2.

Structural and functional characterization of a novel cold-active S-formylglutathione hydrolase (SfSFGH) homolog from Shewanella frigidimarina, a psychrophilic bacterium.

BACKGROUND: S-Formylglutathione is hydrolyzed to glutathione and formate by an S-formylglutathione hydrolase (SFGH) (3.1.2.12). This thiol esterase belongs to the esterase family and is also known as esterase D. SFGHs contain highly conserved active residues of Ser-Asp-His as a catalytic triad at the active site. Characterization and investigation of SFGH from Antarctic organisms at the molecular level is needed for industrial use through protein engineering. RESULTS: A novel cold-active S-formylglutathione hydrolase (SfSFGH) from Shewanella frigidimarina, composed of 279 amino acids with a molecular mass of ~ 31.0 kDa, was characterized. Sequence analysis of SfSFGH revealed a conserved pentapeptide of G-X-S-X-G found in various lipolytic enzymes along with a putative catalytic triad of Ser148-Asp224-His257. Activity analysis showed that SfSFGH was active towards short-chain esters, such as p-nitrophenyl acetate, butyrate, hexanoate, and octanoate. The optimum pH for enzymatic activity was slightly alkaline (pH 8.0). To investigate the active site configuration of SfSFGH, we determined the crystal structure of SfSFGH at 2.32 Å resolution. Structural analysis shows that a Trp182 residue is located at the active site entrance, allowing it to act as a gatekeeper residue to control substrate binding to SfSFGH. Moreover, SfSFGH displayed more than 50% of its initial activity in the presence of various chemicals, including 30% EtOH, 1% Triton X-100, 1% SDS, and 5 M urea. CONCLUSIONS: Mutation of Trp182 to Ala allowed SfSFGH to accommodate a longer chain of substrates. It is thought that the W182A mutation increases the substrate-binding pocket and decreases the steric effect for larger substrates in SfSFGH. Consequently, the W182A mutant has a broader substrate specificity compared to wild-type SfSFGH. Taken together, this study provides useful structure-function data of a SFGH family member and may inform protein engineering strategies for industrial applications of SfSFGH.

A dual synergistic of curcumin and gelatin on thermal-responsive hydrogel based on Chitosan-P123 in wound healing application.

This study aimed to fabricate the potential therapeutic scaffold to efficiently and safely fastening skin wound healing. A biocompatible grafting polymer-based thermal sensitive hybrid hydrogel (Chitosan-P123, CP) containing gelatin and curcumin was designed to be suitable stiffness for tissue regeneration. A detailed in the rheological study found that the encapsulated agents induced the change in the stiffness of the hydrogel from the hard to the soft. Especial, the thermally induced phase transition of CP hydrogel was governed by the participant of gelatin rather than curcumin. For example, at 25 wt% gelatin, CP hydrogel exhibited a unique gel-sol-gel transition following the function of temperature. Moreover, in vitro investigation revealed that the hybrid hydrogel provides the capacity of especially induced curcumin release with a sustainable rate as well as the excellent biocompatibility scaffold. Altogether with in vivo study, the hybrid hydrogel highlighted the advance of the dual synergistic of curcumin and gelatin in development of smart scaffold system, which promoted the efficacy in the regeneration of the structure and the barrier's function of damaged skin such as wound or skin cancer.

Modulating α-synuclein fibril formation using DNA tetrahedron nanostructures.

The small presynaptic protein α-synuclein (α-syn) is involved in the etiology of Parkinson's disease owing to its abnormal misfolding. To date, little information is known on the role of DNA nanostructures in the formation of α-syn amyloid fibrils. Here, the effects of DNA tetrahedrons on the formation of α-syn amyloid fibrils were investigated using various biochemical and biophysical methods such as thioflavin T fluorescence assay, atomic force microscopy, light scattering, transmission electron microscopy, and cell-based cytotoxicity assay. It has been shown that DNA tetrahedrons decreased the level of oligomers and increased the level of amyloid fibrils, which corresponded to decreased cellular toxicity. The ability of DNA tetrahedron to facilitate the formation of α-syn amyloid fibrils demonstrated that structured nucleic acids such as DNA tetrahedrons could modulate the process of amyloid fibril formation. Our study suggests that DNA tetrahedrons could be used as an important facilitator toward amyloid fibril formation of α-synuclein, which may be of significance in finding therapeutic approaches to Parkinson's disease and related synucleinopathies.

Discovery of benzo[e]pyridoindolones as kinase inhibitors that disrupt mitosis exit while erasing AMPK-Thr172 phosphorylation on the spindle.

Aurora kinases play an essential role in mitotic progression and are attractive targets in cancer therapy. The first generation of benzo[e]pyridoindole exhibited powerful aurora kinase inhibition but their low solubility limited further development. Grafting a pyperidine-ethoxy group gives rise to a hydrosoluble inhibitor: compound C5M.C5M could efficiently inhibit the proliferation of cells from different origins. C5M prevented cell cycling, induced a strong mitotic arrest then, cells became polyploid and finally died. C5M did not impair the spindle checkpoint, the separation of the sister chromatids and the transfer of aurora B on the mid-zone. C5M prevented histone H3 phosphorylation at mitotic entry and erased AMPK-Thr172 phosphorylation in late mitosis. With this unique profile of inhibition, C5M could be useful for understanding the role of phospho-Thr172-AMPK in abscission and the relationship between the chromosomal complex and the energy sensing machinery.C5M is a multikinase inhibitor with interesting preclinical characteristics: high hydro-solubility and a good stability in plasma. A single dose prevents the expansion of multicellular spheroids. C5M can safely be injected to mice and reduces significantly the development of xenograft. The next step will be to define the protocol of treatment and the cancer therapeutic field of this new anti-proliferative drug.

Suggestion of suitable animal models for in vivo studies of protein tyrosine phosphatase 1b (PTP1B) inhibitors using computational approaches.

PTP1B is a prototypic enzyme of the superfamily protein tyrosine phosphatases (PTPs) which are critical regulators of tyrosine phosphorylation-dependent signaling events. It is a highly plausible candidate for designing therapeutic inhibitors of obesity and type 2 diabetes (T2D). In this study, a detailed comparative analysis to reveal the evolutionary relationship of human PTP1B among related vertebrates has been addressed. The phylogenetic trees were constructed with maximum likelihood algorithm by PhyML package on the basis of multiple sequence alignment (MSA) by ClustalΩ and T-coffee. Mutational variability of the sequences corresponding to the 3D structure (pdb: 2vev) was analyzed with Consurf software. The comparative analysis by inhibitor docking to different models was made to confirm the suitability of models. As a result, the PTP1B or PTP non-receptor type 1 homologies show high conservativity where about 70% positions on primary structures are conserved. Within PTP domain (3-277), the most variable positions are 12, 13, 19 and 24 which is a part of the second aryl binding site. Moreover, there are important evolutional mutations that can change the conformation of the proteins, for instance, hydrophilic N139 changed to hydrophobic Gly (mPTP1B); E132 to proline in the hydrophobic core structure or Y46 to cystein in pTyr recognition loop. These variations/differences should be taken into account for rational inhibitor design and in choosing suitable animal models for drug testing and evaluation. Moreover, our study suggests critically potential models which are Heterocephalus glaber, Tupaia chinensis, Sus scrofa, and Rattus norvegicus in addition to the best one Macaca fascicularis. Among these models, the H.glaber and R.norvegicus are preferable over M.musculus thanks to their similarity in binding affinity and binding modes to investigated PTP1B inhibitors.

Basal aurora kinase B activity is sufficient for histone H3 phosphorylation in prophase.

Histone H3 phosphorylation is the hallmark of mitosis deposited by aurora kinase B. Benzo[e]pyridoindoles are a family of potent, broad, ATP-competitive aurora kinase inhibitors. However, benzo[e]pyridoindole C4 only inhibits histone H3 phosphorylation in prophase but not in metaphase. Under the C4 treatment, the cells enter into mitosis with dephosphorylated histone H3, assemble chromosomes normally and progress to metaphase, and then to anaphase. C4 also induces lagging chromosome in anaphase but we demonstrated that these chromosome compaction defects are not related to the absence of H3 phosphorylation in prophase. As a result of C4 action, mitosis lasts longer and the cell cycle is slowed down. We reproduced the mitotic defects with reduced concentrations of potent pan aurora kinase as well as with a specific aurora B ATP-competitive inhibitor; we therefore propose that histone H3 phosphorylation and anaphase chromosome compaction involve the basal activity of aurora kinase B. Our data suggest that aurora kinase B is progressively activated at mitosis entry and at anaphase onset. The full activation of aurora kinase B by its partners, in prometaphase, induces a shift in the catalytic domain of aurora B that modifies its affinity for ATP. These waves of activation/deactivation of aurora B correspond to different conformations of the chromosomal complex revealed by FRAP. The presence of lagging chromosomes may have deleterious consequences on the daughter cells and, unfortunately, the situation may be encountered in patients receiving treatment with aurora kinase inhibitors.

Hydrosoluble benzo[e]pyridoindolones as potent inhibitors of aurora kinases.

Aurora kinases play an essential role in mitotic progression and are potentially druggable targets in cancer therapy. We identified benzo[e]pyridoindoles (BePI) as powerful aurora kinase inhibitors. Their efficiency was demonstrated both in enzymatic inhibition studies and in cell culture assays. New BePI molecules were synthesized, and a structure-activity relationship study was conducted with the aim of improving the activity and solubility of the lead compound. Tetracyclic BePI derivatives are characterized by a particular curved shape, and the presence of an oxo group on the pyridine ring was found to be required for aurora kinase B inhibition. New hydrosoluble benzo[e]pyridoindolones were subsequently designed, and their efficacy was tested by a combination of cell-cycle analysis and time-lapse experiments in live cells. The most active BePI derivative, 13 b, inhibited the cell cycle, drove cells to polyploidy, and eventually induced apoptosis. It exhibited high antiproliferative activity in HeLa cells with an IC(50) value of 63 nM. Relative to compounds tested in clinical trials, this antiproliferative potency places 13 b among the top 10 aurora kinase inhibitors. Our results justify further in vivo evaluation in preclinical animal models of cancer.

Study on phylogenetic relationships, variability, and correlated mutations in M2 proteins of influenza virus A.

M2 channel, an influenza virus transmembrane protein, serves as an important target for antiviral drug design. There are still discordances concerning the role of some residues involved in proton transfer as well as the mechanism of inhibition by commercial drugs. The viral M2 proteins show high conservativity; about 3/4 of the positions are occupied by one residue in over 95%. Nine M2 proteins from the H3N2 strain and possibly two proteins from H2N2 strains make a phylogenic cluster closely related to 2RLF. The variability range is limited to 4 residues/position with one exception. The 2RLF protein stands out by the presence of 2 serines at the positions 19 and 50, which are in most other M2 proteins occupied by cysteines. The study of correlated mutations shows that there are several positions with significant mutational correlation that have not been described so far as functionally important. That there are 5 more residues potentially involved in the M2 mechanism of action. The original software used in this work (Consensus Constructor, SSSSg, Corm, Talana) is freely accessible as stand-alone offline applications upon request to the authors. The other software used in this work is freely available online for noncommercial purposes at public services on bioinformatics such as ExPASy or NCBI. The study on mutational variability, evolutionary relationship, and correlated mutation presented in this paper is a potential way to explain more completely the role of significant factors in proton channel action and to clarify the inhibition mechanism by specific drugs.

Top-hits for H1N1pdm Identified by Virtual Screening Using Ensemble-based Docking.

A list of 27 promising antiviral drugs is proposed for use against the H1N1pdm strain. Since the binding site of the H1N1pdm neuraminidase is similar to that of the bird flu H5N1, an effective means to quickly identify top candidates for use against H1N1pdm is to use known bird-flu drugs and the 27 compounds from the NCI diversity set which bind best to H5N1 neuraminidase. These compounds serve as viable candidates for docking against the H1N1pdm neuraminidase, using ensembles extracted from molecular dynamics simulations of the H1N1pdm system. The ranking order of these top candidates was found to be different from the previously published results for H5N1. The results indicated that the Oseltamivir (Tamiflu) and Peramivir drugs have higher ranking than Zanamivir (Relenza). However, six drug candidates were found to bind more effectively to H1N1pdm neuraminidase than Tamiflu. Detailed hydrogen bond network analysis for these six candidates is also provided.