Influenza type A is considered as a severe public health concern. The mechanism of drugs applied for the control of this virus depends on two surface glycoproteins with antigenic properties, namely hemagglutinin (HA) and neuraminidase (NA). HA aids the virus to penetrate cells in the early stage of infection and NA is an enzyme with the ability to break glycoside bonds, which enables virion spread through the host cell membrane. Since NA contains a relatively preserved active site, it has been an important target in drug design. Oseltamivir is a common drug used for the treatment of influenza infections, for which cases of resistance have recently been reported, giving rise to health concerns. Flavonoids are natural polyphenolic compounds with potential blocking effects in the neuraminidase active site. Based on their antiviral effect, the flavonoids quercetin, catechin, naringenin, luteolin, hispidulin, vitexin, chrysin and kaempferol were selected in the present study and compared alongside oseltamivir on molecular docking, binding energy and active site structure, in order to provide insight on the potential of these compounds as targeted drugs for the control and treatment of influenza type A. The molecular characterization of flavonoids with binding affinity was performed using AutoDock Vina software. The results indicated that these compounds may effectively block the NA active site. Therefore, these natural compounds derived from fruits have the potential for development into drugs for controlling influenza, which may aid overcome the clinical challenge of the H1N1 strain epidemic.
influenza type Aneuraminidaseflavonoidsmolecular dockingIntroduction
Influenza type A virus is historically among the most aggressive virus' to have been identified, having caused several global epidemics. Between 1918 and 1919 this virus caused a mortality rate accounting for 5% of the world population (1–3). Unfortunately, the potential hazard and spread of the virus remain to be considered as a threat to human health (4). Surface glycoproteins with antigenic properties including hemagglutinin (HA) and neuraminidase (NA) account for the emergence of novel subtypes of type A influenza and of novel forms of influenza. To date, 18 subtypes of HA (H1-H18) and 11 subtypes of NA (N1-N11) have been identified (5,6). Hemagglutinin binds to sialic acid receptor for viral entry into the host cell (7). NA is an enzyme with the ability to break glycoside bonds, which may block cell-surface sialic acid receptor binding and allow virion spread across the host cell membrane. In the absence of this mechanism, virions accumulate on the cell membrane (8). In this manner NA may serve a facilitating role in the early stage of infection into endothelial cells (9).
Due to the importance of NA in virus distribution and its relatively preserved active site, it is an important target of drug design for control of influenza infections (10). Among previously developed drugs, oseltamivir, which targets NA, has been identified as the most effective, with an inhibitory effect on the N1, N2 and N9 subtypes (11,12). However, several cases of resistance reported for this drug raises concern, since lack of disease control can be fatal with the potential for epidemic spread (3,13–15).
Flavonoids are naturally occurring, polyphenolic compounds of relatively low molecular weight that exhibit benzopyrone structure, produced as secondary metabolites in plants to stimulate growth during the state of tension (16,17). Due to the reported antioxidant, antibacterial and antiviral effects, presence in normal daily diet and minimal side-effects of flavonoids, they are considered useful resources for drug design (8,18–22). Initial reports on the use of flavonoids as pharmaceutical compounds emerged in the 1940s, when the potential of these compounds was evaluated in the prevention and treatment of coronary heart disease and for their anti-cancer, anti-inflammatory, antibacterial and antiviral effects (23–29). Due to the diversity of flavonoids, the use of computational designs and molecular information on flavonoid compounds enables the selection of the optimum compounds in investigations aimed at identifying alternative drugs with greater impact.
Compared with laboratory techniques, molecular data-based methods present the possibility to identify compounds with potent inhibitory effects, predict their structure and design effective drugs in a rapid and high-throughput approach. Compounds predicted through such computational methods require approval via classic experimental procedures (30–32). In fact, the study of molecular data offers the potential to select the most suitable compounds for treatment. Such a benefit is not possible without the need for a glossary of information, which would otherwise be time-consuming. In this regard, docking as the first medical computational design tool provides a method based on the prediction of binder energy interactions (using platforms including MetaPocket, GRAMM-X and HEX Protein Docking). The optimum orientation and interaction of the molecules with the lowest binding energy can be thus determined with the three-dimensional information index of the protein-protein binding (33–35).
Due to the abundance of research on the antiviral effects of flavonoids, in the present study, a molecular analysis of the binding affinity of various flavonoids was conducted by the docking method to determine the potential of these compounds as drug candidates for the control and treatment of influenza type A. The results indicated that the binding energy of the majority of flavonoids studied was high at the active site of neuraminidase N1 compared with oseltamivir, which may provide a basis for the replacement of oseltamivir with such flavonoids-based compounds.
Materials and methodsMolecular analysis protocol
Protein-ligand reactions, binding mode and binding affinity of neuraminidase N1 with oseltamivir and the selected flavonoids were determined via the docking method through the following steps:
Proteins
The Protein Data Bank (PDB) file of the region of neuraminidase N1 which interacts with oseltamivir (PDB code:3ti6) was obtained. Water molecules were removed and non-polar hydrogen molecules were added to the protein structure in Python molecule viewer setting of PyMOL v1.9 (http://www.pymol.org), which obtained the active site of the protein with a size of 24 × 34 × 22 at 1Å grid spacing.
Ligands
CID files and molecular ligand data were obtained from PubChem on the NCBI database (https://pubchem.ncbi.nlm.nih.gov/search/search.cgi) (Table I). Ligands were converted to pdbqt format in Open Babel (https://openbabel.org/) GUI 2.3.2a setting.
Protein-ligand docking
Docking was performed in AutoDock Vina (http://vina.scripps.edu/), which predicts interactions between small molecules and proteins.
Determination of protein active site
The amino acids relevant to binding in the active site of the protein were determined using MetaPocket 2.0 (projects.biotec.tu-dresden.de/metapocket).
Analysis
Results were viewed and analyzed with PyMOL v1.9.
ResultsCompound binding
The results indicated that all tested flavonoid compounds (Fig. 1) exhibited relatively high potential for binding to the active site of neuraminidase N1 (Table II). Notably their binding energy was higher in comparison to that of oseltamivir. It is noteworthy that the docking site of these compounds and oseltamivir were similar
Compound docking
Results of compound docking in PyMol viewer revealed that all the tested flavonoid compounds have affinity for the active site of neuraminidase N1, similar to oseltamivir (Fig. 2).
Based on previous literature, several flavonoid compounds were initially selected (16). Compounds including caffeic acid and caffeine were also included in the computational investigation, but the obtained results were insignificant compared with the reported compounds, and only the docking results of the compounds with highest binding energy were reported (Table II). The obtained docking results confirmed the capability of the selected compounds in terms of effective molecular interaction with the investigated protein, when compared with oseltamivir; selected compounds were determined to have stronger affinities in binding with the active site of neuraminidase.
Protein active site
Amino acid residues involved in the interaction between the protein region and studied ligands were determined using MetaPocket 2.0 software (Fig. 3 and Table III).
Since the active site of neuraminidase N1, which interacts with oseltamivir, has a lower affinity than with the flavonoids, the selection of flavonoids in this study appeared successful. The selection of more optimal molecules exhibit a binding position with greater affinity may provide novel pharmaceutical approaches. In addition, the replacement of oseltamivir with flavonoids such as those currently tested when regarding the interaction with neuraminidase N1 may be a topic of interest, because of the importance of antimicrobials in the treatment of influenza type A, and flavonoids have been found to be effective against Gram-positive bacteria (8,16).
Discussion
In recent years, among subtypes of influenza type A, H1N1 strain has been considered a challenge in terms of disease control and is a concern of national and international health organizations (13,29). Several cases of mutations in the active site of neuraminidase have been reported, and may lead to resistance to oseltamivir, as a drug commonly administered for influenza (14,15). These cases highlight the need for broader research aimed at determining alternative therapeutic methods.
The aim of the present study was to provide a basis for in vitro experiments from the perspective of protein-ligand molecular computing, for which AutoDock Vina software was used. Despite the rapidity of calculations, the quality of results obtained from AutoDock Vina software is comparable to those from AutoDock software (31,32). Another aim was to shed light on the public notion regarding the potential medicinal effects of flavonoid consumption on influenza type A (36,37).
Flavonoids are currently identified not as medicine, but as necessary elements of daily diet that aid functioning of the immune system (19,20). Due to their antioxidant activity, flavonoid compounds may regulate systemic anti-inflammatory reactions (16). Thus, applying them as a classical medicine does not appear a priority. Therefore, the present study aimed to highlight the overall effects of flavonoids as a dietary factor on the prevention, control and treatment of influenza (38,39).
With the development of nanobiotechnological drug design, targeted delivery of flavonoid compounds to target tissues now appears a possibility (40). Via this method, the systemic side-effects of flavonoids during drug intake may be prevented and the therapeutic dose of flavonoids may be delivered directly to the targeted tissue, leading to improved efficiency of these drugs in lower doses.
Determining precisely the interactions in an active analogous position with antibiotics and candidate flavonoids may aid to better determine the alternatives. Since the variety of antibiotic sources useful for treatment is limited, and in recent decades there has been a lack of studies introducing novel antibiotics, the currently available antibiotics are becoming increasingly less effective due to the high rates of mutation in viruses and bacteria (41–44). The performance of comprehensive studies to identify novel sources with the potential of replacing or at least complementing current antibiotics is considered a priority in the pharmaceutical industry. There are numerous antioxidant sources with antibiotic potential against viruses in nature, including fruits, vegetables and honey (45). So far, their established antioxidant effects have often been limited to applications in traditional medicine (46–49). Among these potentially antioxidant compounds are flavonoids; the study of these compounds is thus warranted to select antibiotic candidates (50).
The establishment of novel drugs against influenza type A, for application when mutations in the binding site of the drug occur, appears imperative. Since the binding energies and active site interactions of compounds studied herein were comparable to oseltamivir, further evaluation is recommended with the ultimate aim of introducing a novel generation of drugs in this field.
Acknowledgements
The authors thank Dr Ruhollah Dorostkar at Baqiyatallah University of Medical Sciences, Tehran, Iran, for providing suggestions on the article and servers used.
Funding
No funding was received.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Authors' contributions
SMS conceived and designed the study. NG analyzed the data. SR wrote the article. MSH performed the laboratory and server procedures. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
ReferencesAlmondDIs the 1918 Influenza pandemic over? Long-term effects of in utero Influenza exposure in the post-1940 US populationJ Polit Econ114672712200610.1086/507154TaubenbergerJKMorensDM1918 Influenza: The mother of all pandemicsEmerg Infect Dis121522200610.3201/eid1209.05-097916494711NeumannGNodaTKawaokaYEmergence and pandemic potential of swine-origin H1N1 influenza virusNature459931939200910.1038/nature0815719525932NeumannGKawaokaYTransmission of influenza A virusesVirology 479-480234246201510.1016/j.virol.2015.03.009MallajosyulaVVCitronMFerraraFLuXCallahanCHeideckerGJSarmaSPFlynnJATempertonNJLiangXInfluenza hemagglutinin stem-fragment immunogen elicits broadly neutralizing antibodies and confers heterologous protectionProc Natl Acad Sci USA111E2514E2523201410.1073/pnas.140276611124927560TrottOOlsonAJAutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreadingJ Comput Chem31455461201019499576WeisWBrownJHCusackSPaulsonJCSkehelJJWileyDCStructure of the influenza virus haemagglutinin complexed with its receptor, sialic acidNature333426431198810.1038/333426a03374584SibiGDhananjayaKBhimanagoudaRPMolecular docking and receptor based analysis of H1N1 antiviral drugsBioresearch Bull21091132011WangYTChanCHSuZYChenCLHomology modeling, docking, and molecular dynamics reveal HR1039 as a potent inhibitor of 2009 A(H1N1) influenza neuraminidaseBiophys Chem1477480201010.1016/j.bpc.2009.12.00220045243SeniyaCKhanGJMisraRVyasVKaushikSIn-silico modelling and identification of a possible inhibitor of H1N1 virusAsian Pac J Trop Dis4S467S476201410.1016/S2222-1808(14)60492-8GovorkovaEAIlyushinaNAMcClarenJLNaiposposTSDouangngeunBWebsterRGSusceptibility of highly pathogenic H5N1 influenza viruses to the neuraminidase inhibitor oseltamivir differs in vitro and in a mouse modelAntimicrob Agents Chemother5330883096200910.1128/AAC.01667-0819349520ParkJWJoWHComputational design of novel, high-affinity neuraminidase inhibitors for H5N1 avian influenza virusEur J Med Chem45536541201010.1016/j.ejmech.2009.10.04019914748DuQSWangSQHuangRBChouKCComputational 3D structures of drug-targeting proteins in the 2009-H1N1 influenza A virusChem Phys Lett485191195201010.1016/j.cplett.2009.12.037GuptaCLAkhtarSBajpaibPKandpalKNDesaiGSTiwariAKComputational modeling and validation studies of 3-D structure of neuraminidase protein of H1N1 influenza A virus and subsequent in silico elucidation of piceid analogues as its potent inhibitorsEXCLI J12215225201326417228BloomJDGongLIBaltimoreDPermissive secondary mutations enable the evolution of influenza oseltamivir resistanceScience32812721275201010.1126/science.118781620522774KumarSPandeyAKChemistry and biological activities of flavonoids: an overviewScientificWorldJournal2013162750201310.1155/2013/16275024470791DixonRADeyPMLambCJPhytoalexins: Enzymology and molecular biologyAdv Enzymol Relat Areas Mol Biol55113619836353887CushnieTPLambAJAntimicrobial activity of flavonoidsInt J Antimicrob Agents26343356200510.1016/j.ijantimicag.2005.09.00216323269OrhanDDOzçelikBÖzgenSErgunFAntibacterial, antifungal, and antiviral activities of some flavonoidsMicrobiol Res165496504201010.1016/j.micres.2009.09.00219840899CushnieTPLambAJRecent advances in understanding the antibacterial properties of flavonoidsInt J Antimicrob Agents3899107201110.1016/j.ijantimicag.2011.02.01421514796HavsteenBHThe biochemistry and medical significance of the flavonoidsPharmacol Ther9667202200210.1016/S0163-7258(02)00298-X12453566MiddletonEJrKandaswamiCTheoharidesTCThe effects of plant flavonoids on mammalian cells: Implications for inflammation, heart disease, and cancerPharmacol Rev52673751200011121513BandarukYMukaiRTeraoJCellular uptake of quercetin and luteolin and their effects on monoamine oxidase-A in human neuroblastoma SH-SY5Y cellsToxicol Rep1639649201410.1016/j.toxrep.2014.08.01628962277LiBQFuTDongyanYMikovitsJARuscettiFWWangJMFlavonoid baicalin inhibits HIV-1 infection at the level of viral entryBiochem Biophys Res Commun276534538200010.1006/bbrc.2000.348511027509AlamMASubhanNRahmanMMUddinSJRezaHMSarkerSDEffect of citrus flavonoids, naringin and naringenin, on metabolic syndrome and their mechanisms of actionAdv Nutr5404417201410.3945/an.113.00560325022990GanugapatiJMukkavalliSSahithiADocking studies of green tea flavonoids as insulin mimeticsInt J Comput Appl3048522011JeongHJRyuYBParkS-JKimJHKwonH-JKimJHParkKHRhoMCLeeWSNeuraminidase inhibitory activities of flavonols isolated from Rhodiola rosea roots and their in vitro anti-influenza viral activitiesBioorg Med Chem1768166823200910.1016/j.bmc.2009.08.03619729316GanugapatiJBaldwaALalaniSMolecular docking studies of banana flower flavonoids as insulin receptor tyrosine kinase activators as a cure for diabetes mellitusBioinformation8216220201210.6026/9732063000821622493522LuSJChongFCCombining molecular docking and molecular dynamics to predict the binding modes of flavonoid derivatives with the neuraminidase of the 2009 H1N1 influenza A virusInt J Mol Sci1344964507201210.3390/ijms1304449622605992HuangSYZouXAdvances and challenges in protein-ligand dockingInt J Mol Sci1130163034201010.3390/ijms1108301621152288SeeligerDde GrootBLLigand docking and binding site analysis with PyMOL and Autodock/VinaJ Comput Aided Mol Des24417422201010.1007/s10822-010-9352-620401516MorrisGMHueyROlsonAJUsing AutoDock for ligand-receptor dockingCurr Protoc Bioinformatics248.14.18.14.402008RamharackPSolimanMESBioinformatics-based tools in drug discovery: the cartography from single gene to integrative biological networksDrug Discov Today2316581665201810.1016/j.drudis.2018.05.04129864527MunsamyGSolimanMESHomology modeling in drug discovery-an update on the last decadeLett Drug Des Discov1410991111201710.2174/1570180814666170110122027TsujikawaHSatoKWeiCSaadGSumikoshiKNakamuraSTeradaTShimizuKDevelopment of a protein-ligand-binding site prediction method based on interaction energy and sequence conservationJ Struct Funct Genomics173949201610.1007/s10969-016-9204-227400687LiuALWangHDLeeSMWangYTDuGHStructure-activity relationship of flavonoids as influenza virus neuraminidase inhibitors and their in vitro anti-viral activitiesBioorg Med Chem1671417147200810.1016/j.bmc.2008.06.04918640042LuriaSEDelbrückMMutations of bacteria from virus sensitivity to virus resistanceGenetics28491511194317247100RussoMSpagnuoloCTedescoIBilottoSRussoGLThe flavonoid quercetin in disease prevention and therapy: Facts and fanciesBiochem Pharmacol83615201210.1016/j.bcp.2011.08.01021856292KuoSMDietary flavonoid and cancer prevention: Evidence and potential mechanismCrit Rev Oncog84769199710.1615/CritRevOncog.v8.i1.309516086ChenYWuQSongLHeTLiYLiLSuWLiuLQianZGongCPolymeric micelles encapsulating fisetin improve the therapeutic effect in colon cancerACS Appl Mater Interfaces7534542201510.1021/am506689325495760DrakeJWThe distribution of rates of spontaneous mutation over viruses, prokaryotes, and eukaryotesAnn N Y Acad Sci870(1 MOLECULAR STR)100107199910.1111/j.1749-6632.1999.tb08870.x10415476ElenaSFLenskiREEvolution experiments with microorganisms: The dynamics and genetic bases of adaptationNat Rev Genet4457469200310.1038/nrg108812776215DuffySShackeltonLAHolmesECRates of evolutionary change in viruses: Patterns and determinantsNat Rev Genet9267276200810.1038/nrg232318319742YogeshMSaurabhSFlavonoid: a pharmacologically significant scaffoldJ Pharm Pharm Sci64885042017Eteraf-OskoueiTNajafiMTraditional and modern uses of natural honey in human diseases: A reviewIran J Basic Med Sci16731742201323997898Viuda-MartosMRuiz-NavajasYFernández-LópezJPérez-ÁlvarezJAFunctional properties of honey, propolis, and royal jellyJ Food Sci73117124200810.1111/j.1750-3841.2008.00966.xGuptaPDBirdiTJDevelopment of botanicals to combat antibiotic resistanceJ Ayurveda Integr Med8266275201710.1016/j.jaim.2017.05.00428869082StankovićNMihajilov-KrstevTZlatkovićBStankov-JovanovićVMitićVJovićJČomićLKocićBBernsteinNAntibacterial and antioxidant activity of traditional medicinal plants from the Balkan PeninsulaNJAS Wagening J Life Sci782128201610.1016/j.njas.2015.12.006LeeJHChoSPaikHDChoiCWNamKTHwangSGKimSKInvestigation on antibacterial and antioxidant activities, phenolic and flavonoid contents of some thai edible plants as an alternative for antibioticsAsian-Australas J Anim Sci2714611468201410.5713/ajas.2013.1362925178298GonçalvesEMoreiraCGrossoPAndradeBValentãoPRomanoAPhenolic profile, antioxidant activity and enzyme inhibitory activities of extracts from aromatic plants used in Mediterranean dietJ Food Sci Technol54219227201710.1007/s13197-016-2453-z28242919
Chemical formulas of the studied flavonoids and oseltamivir.
Docking results obtained with AutoDock Vina using PyMOL indicated that similar to oseltamivir, the selected flavonoids have affinity to bind with the activity site of neuraminidase (based on the results in Table II).
Active site of protein with oseltamivir determined using MetaPocket 2.0.
Molecular data of studied flavonoids and oseltamivir.