Antiviral activity of stachyflin on influenza A viruses of different hemagglutinin subtypes
Virology Journal volume 10, Article number: 118 (2013)
The hemagglutinin (HA) of influenza viruses is a possible target for antiviral drugs because of its key roles in the initiation of infection. Although it was found that a natural compound, Stachyflin, inhibited the growth of H1 and H2 but not H3 influenza viruses in MDCK cells, inhibitory activity of the compound has not been assessed against H4-H16 influenza viruses and the precise mechanism of inhibition has not been clarified.
Inhibitory activity of Stachyflin against H4-H16 influenza viruses, as well as H1-H3 viruses was examined in MDCK cells. To identify factors responsible for the susceptibility of the viruses to this compound, Stachyflin-resistant viruses were selected in MDCK cells and used for computer docking simulation.
It was found that in addition to antiviral activity of Stachyflin against influenza viruses of H1 and H2 subtypes, it inhibited replication of viruses of H5 and H6 subtypes, as well as A(H1N1)pdm09 virus in MDCK cells. Stachyflin also inhibited the virus growth in the lungs of mice infected with A/WSN/1933 (H1N1) and A/chicken/Ibaraki/1/2005 (H5N2). Substitution of amino acid residues was found on the HA2 subunit of Stachyflin-resistant viruses. Docking simulation indicated that D37, K51, T107, and K121 are responsible for construction of the cavity for the binding of the compound. In addition, 3-dimensional structure of the cavity of the HA of Stachyflin-susceptible virus strains was different from that of insusceptible virus strains.
Antiviral activity of Stachyflin was found against A(H1N1)pdm09, H5, and H6 viruses, and identified a potential binding pocket for Stachyflin on the HA. The present results should provide us with useful information for the development of HA inhibitors with more effective and broader spectrum.
Influenza A virus is widely distributed in birds and mammals, including humans, and causes seasonal and pandemic influenza. For the prevention and therapy of influenza, anti-influenza drugs play an important role in addition to vaccination. Currently available anti-influenza virus drugs are M2 inhibitors (Amantadine and Rimantadine) and neuraminidase (NA) inhibitors (Oseltamivir, Zanamivir, Peramivir, and Laninamivir) [1, 2]. The M2 transfers protons into the infecting virus in the endosome at low pH, and the M1 is dissociated from the genome-transcriptase complex . M2 inhibitors block ion channel activity and inhibit the influx of protons, thereby exhibiting antiviral activity against influenza A viruses . At the end of the virus life cycle, the NA catalyses the cleavage of terminal sialic acid from glycoconjugates on the host cell surface to release progeny virions . Of these anti-influenza virus drugs, the NA inhibitors, which interfere with the release of the virus, are used clinically since they are broadly effective ; however, concern has been raised because of the isolation of NA inhibitor-resistant viruses from clinical samples . Therefore, it is expected to develop drugs targeting other virus proteins than the NA and M2.
Hemagglutinin (HA) is a surface glycoprotein of influenza A virus, and is a possible target of antiviral drugs because of its key roles in the initiation of infection. Each monomer of the trimeric HA is composed of 2 subunits, HA1 and HA2. The HA1 has a receptor binding domain, and the HA2 mediates the fusion of the virus envelope with the cellular membrane . Several studies have identified compounds which inhibit viral infection by blocking the binding of the HA to sialic acid receptor on the host cell surface (cyanovirin-N and trisphenol-sialyllactose) or fusion step (TBHQ, BMY-27709, CL-385319, and N-carboxamide) [8–12]; however, for many of these inhibitors, the antiviral spectrum is limited to the HA of certain subtypes, so that they have not been used clinically. To develop more effective HA inhibitors, further investigations of these HA inhibitors and the analysis of the attachment and fusion steps of influenza virus infection in the host cells are needed.
It was found that a sesquiterpene derivative, Stachyflin, inhibited replication of H1 and H2 influenza A viruses in vitro[13, 14] and in vivo[15, 16]. Although Stachyflin is postulated to inhibit the fusion step, its precise mechanism has not been clarified. In the present study, it is revealed that Stachyflin inhibit the growth of H1, H2, H5, and H6 influenza viruses by binding the site of the HA2 and preventing the HA from fusion of the virus envelope with cellular membrane.
Antiviral activity of stachyflin in vitro and in vivo
The antiviral spectrum of Stachyflin was determined by measuring its inhibitory effect on the replication of 31 influenza virus strains of H1-H16 subtypes in Madin-Darby canine kidney (MDCK) cells. The antiviral effects were evaluated in various concentrations of Stachyflin up to 6.50 μM by virus-induced cytopathic effects (CPE). Stachyflin inhibited the replication of H1 including A(H1N1)pdm09 virus, H2, H5, and H6 subtype influenza virus strains, but not that of the other subtype strains. Susceptibility of the viruses to Stachyflin varied with the strains (Table 1). In all of the viruses tested, AWSN/1933 (H1N1) (WSN) showed the highest susceptibility to Stachyflin.
The antiviral activity of Stachyflin in mice challenged with A/Kumamoto/5/1967 (H2N2) was evaluated in previous study [15, 16]. To confirm whether the antiviral activity in vitro can be applied to in vivo, the virus titers in the lungs of mice infected with WSN or A/chicken/Ibaraki/1/2005 (H5N2) (Ibaraki) 72 h post-inoculation were evaluated. WSN showed efficient replication and was lethal to mice, while infection with Ibaraki, a low pathogenic H5 avian influenza virus, was not lethal. Although mice treated with Stachyflin did not reduce the weight loss induced by virus challenge, only both groups of mice treated with Stachyflin at 100 mg/kg/day showed significantly lower mean virus titers in the lungs than each control group (Figure 1). In mice infected with WSN, the mean virus titer in the lungs of mice treated with Stachyflin at 100 mg/kg/day was 103.5 50% tissue culture infectious dose per gram (TCID50/g) whereas that of control mice was 106.2 TCID50/g. In mice infected with Ibaraki, the mean virus titer in the lungs of mice treated with Stachyflin at 100 mg/kg/day was 103.4 TCID50/g, whereas that of control mice was 105.4 TCID50/g. Thus, Stachyflin showed inhibitory activity on virus growth in vitro and in vivo. In these experiments, CPE induced by Stachyflin was not observed even in the highest concentration, 6.50 μM in MDCK cells and weight changes of mice were not observed even in administration of 100 mg/kg/day for 3 days.
Selection of stachyflin-resistant virus clones
To determine the amino acids which contribute to the susceptibility of the viruses to Stachyflin, Stachyflin-resistant virus clones were selected from WSN, A/Puerto Rico/8/1934 (H1N1) (PR8), Ibaraki, and A/duck/Taiwan/4801/1990 (H6N5) (Taiwan). Six Stachyflin-resistant virus clones were selected from WSN (WSN R1-R6) and 2 clones from PR8 (PR8 R1-R2) by single passage in MDCK cells in the presence of 0.52 or 1.30 μM of Stachyflin. The frequency of the Stachyflin-resistant virus clones was around 10-3.0-10-4.0. Meanwhile from Ibaraki and Taiwan, Stachyflin-resistant virus clones (Ibaraki R1 and Taiwan R1) were selected by 3 passages in the presence of 1.30 μM of Stachyflin. The clones were plaque-purified and propagated in MDCK cells in the presence of Stachyflin. The replication of these Stachyflin-resistant virus clones was not inhibited even with 6.50 μM Stachyflin (Table 2).
Nucleotide sequences of the HA genes of the wild type and Stachyflin-resistant virus clones were determined. All of the mutants had a single amino acid substitution in the HA2 (Table 2). The number of amino acid residue of HA2 starts from GLF motif and is common in these strains. Of these amino acid substitutions, K51R was common in the HAs of Stachyflin-resistant virus clones of WSN, Ibaraki, and Taiwan. These amino acid substitutions were mapped on the structure of the HA monomer (Figure 2). Although each of the amino acid substitutions was observed in the stem region of the HA2, it was impossible to form a binding site for Stachyflin by the distance. To confirm that each single mutation was responsible for Stachyflin resistance, rgWSN mutants, which have one amino acid substitution of the mutants, were generated by reverse genetics and site-directed mutagenesis and were characterized. The replication of rgWSN mutants was not inhibited by 6.50 μM Stachyflin, indicating that all the amino acid substitutions were responsible for Stachyflin resistance (Table 3).
Optimal pH for hemolysis of stachyflin-resistant virus clones
Influenza virus mediates the hemolysis of chicken red blood cells (cRBC), which has been thought to represent the fusion of the virus envelope with cellular membrane . Using a hemolysis assay, the effect of Stachyflin on the fusion of WSN wild type and Stachyflin-resistant virus clones was assessed. Stachyflin inhibited the hemolysis of cRBC induced by the wild type virus but not that by the mutants. In addition, optimal pH for fusion, at which 50% hemolysis occurred, shifted from 6.0 for the wild type virus to as follows: WSN R1: 6.3, R2: 5.7, R4: 6.2, and R5: 5.8 (Table 2).
Relationship of amino acid substitution and the structure of potential binding pocket for stachyflin in the HA
To predict the possible docking model for Stachyflin in the HA trimer of WSN, PR8, Ibaraki, and Taiwan, computer docking simulations of Stachyflin with these HAs were performed. On the surface of these HA trimers, first, the binding pocket for Stachyflin was supposed to be located in a large cavity on the HA2 region, because most amino acid substitutions found on the HA of Stachyflin-resistant virus clones were in this cavity. In the cavity, we found a potential binding pocket for Stachyflin, which was formed by helix A and helix D of the HA2 subunit (Figure 3A, B). This binding pocket contained the residues, D37, K51, and T107, which were substituted in the HAs of Stachyflin-resistant virus clones (Figure 3B). In addition, a residue identified as a Stachyflin-resistant mutation previously , K121, was also contained in the region of the binding pocket (Figure 3B). It was also found that K51 and T107 made a hydrogen bond between helix A and helix D, which may stabilize the structure of the binding pocket.
Using computer docking modeling, it was investigated that how Stachyflin makes bonds with the amino acid around the binding pocket. In the present study, 2 possible docking models of Stachyflin and the HA were proposed (Figure 3B). In one model represented by orange-colored Stachyflin, Stachyflin bound to the site in the vicinity of T107 in the binding pocket, which is similar to that in a previous report  (Figure 3B). In the other model represented by yellow-colored Stachyflin, Stachyflin bound directly to D37 and K121 (Figure 3B). Both models were different from that in a previous study which postulated that Stachyflin forms a hydrogen bond with both K51 (helix A) and K121 (helix D) .
Anti-influenza virus drugs are important for the prevention and treatment of seasonal and pandemic influenza. The HA inhibitor is a candidate drug which inhibits virus attachment to or penetration into the host cells. Most fusion inhibitors hitherto reported had H1 and H2 or H3 subtype-specific antiviral activity and have not been examined their activities against H4-H16 viruses [10, 11, 19], except for CL-385319 . Stachyflin was also reported as H1 and H2 subtype-specific fusion inhibitor and its 50% inhibitory concentration (IC50) was 0.2-0.6 μM . The results of the present study revealed that Stachyflin had subtype-specific antiviral activities against not only H1 and H2 viruses, but also H5, including highly pathogenic avian influenza viruses, and H6 viruses, as well as A(H1N1)pdm09 virus in MDCK cells and its IC50 was 0.05-4.7 μM (Table 1). It is revealed that cytotoxicity of Stachyflin to the MDCK monolayers did not appear up to a concentration of 75 μM , however, for its insolubility , antiviral activity in vitro was assessed up to a concentration of 6.5 μM. In the present study, it was found that WSN strain was the most susceptible to Stachyflin, and then Ibaraki, an avian H5N2 virus isolated from chicken. Then, antiviral activities of the compound were evaluated against these viruses in a mouse model. It was previously revealed that the activity of Stachyflin was limited as about 40% of viruses were recovered from lungs of mice injected intraperitoneally with 2 mg/mouse/day (about 100 mg/kg/day) of Stachyflin compared to non-injected mice after the challenge of A/Kumamoto/5/1967 (H2N2) and 400 mg/kg/day of Stachyflin by intraperitoneal injection was not toxic to mice [15, 16]. Stachyflin showed antiviral activity to reduce 102.0–3.0 virus titer in lungs of mice against H1 and H5 viruses at a dose of 100 mg/kg/day. NA inhibitors, which are used clinically, showed efficient antiviral activity in mice at a dose of 20 mg/kg/day  whereas Stachyflin showed the same effect at a dose of 100 mg/kg/day, which is considered an overdose; therefore, in addition to the poor pharmacokinetic of Stachyflin [15, 16] and limited spectrum, it may be difficult to apply Stachyflin in clinical use in the present form. However, Stachyflin may be clinically used in combination with some NA inhibitor such as Oseltamivir.
Antiviral activity of Stachyflin was related with the structure of the HA. The structure of H1, H2, and H5 HAs, which are susceptible to Stachyflin, closely resemble each other  and these HAs including H6 were identified as group 1 (H1, H2, H5, H6, H8, H9, H11, H12, H13, and H16) by phylogenetic groupings of HA. The viruses used in this study have a similar sequence and structure of the binding pocket for the compound on the HA; for example, the structure of the binding pocket of the H1 HA is similar to that of the H5 HA compared to that of the H3 HA (Figure 4A, B). There are 6 different amino acids between the H1 and H5 HA around this region (Figure 4A). In particular, inside the binding pocket, only one amino acid at position 43 in the HA2 is different: WSN: asparagine, Ibaraki: lysine, which is assumed to cause the difference in the susceptibility to Stachyflin due to the difference in the size and charge of their side chains. For example, lysine has a larger side chain than asparagine and may make it more difficult for Stachyflin to enter into the binding pocket; therefore, the susceptibility to Stachyflin of Ibaraki was lower than that of WSN (Table 1). On the other hand, the HAs of the virus strains insusceptible to Stachyflin have different amino acid sequences in the binding pocket from that of the susceptible ones. For instance, 14 amino acids were different between the H1 and H3 HAs in the vicinity of the binding pocket, which cause a structural difference between these HAs (Figure 4B).
In the previous reports, 2 docking models of the HA and Stachyflin were suggested [14, 18]. In one of the models, Stachyflin was postulated to make hydrogen bonds with K51 and K121 in the HA2 directly. The structure of the H3 HA, which is not susceptible to Stachyflin, was used to make the docking model; therefore, the model differed from our model, which uses the HA of susceptible viruses to Stachyflin. The other docking model was similar to one of those proposed in the present study, which is indicated as orange-colored Stachyflin (Figure 3B). However, since the software program for docking simulation was different from that used in the present study, several differences were found between these models. Indeed, the model which differed from both the previous docking models was also shown in the present study. In this study, we were unable to judge which of these models is more feasible. To further clarify these discrepancies, it is necessary to perform an X-ray crystallographic analysis of Stachyflin complex with the HAs to define the binding site for Stachyflin.
In most studies of HA inhibitors, mutations found in the HAs of inhibitor-resistant virus are explained as the cause to reduce their binding affinity with the compound and stabilization of the HA [10, 19]. In the present study, it was indicated that D37, K51, D85, I91, L98, and T107 were involved in binding affinity with Stachyflin of the HA by the selection of Stachyflin-resistant virus clones. On the basis of the computer modeling in this study, K51 and T107 are postulated to make a hydrogen bond, which may stabilize the structure of the binding pocket for Stachyflin; therefore, mutations of these amino acids should lead to loss of this hydrogen bond, which may destabilize and distort the binding pocket and decrease the binding affinity of the HA to Stachyflin. Indeed, we selected Stachyflin-resistant virus clones which have the amino acid substitutions, K51R and T107I. In addition, we also identified amino acid substitutions, D37N and K121E, which are located close to each other on the HA. Interestingly, some of the crystal structure of the HA shows the possibility that D37 and K121 make a hydrogen bond via a water molecule. Then, similar to K51 and T107, the binding between D37 and K121 stabilizes the structure of the binding pocket and mutations of these amino acids lead to distortion of the binding pocket,decreasing the binding affinity of the HA to Stachyflin. Indeed, it was indicated by shift of fusion pH of the mutants that the amino acid substitutions, D37N, K51R, and T107I, may change the stability of the HA  (Table 2). The other possibility is that D37 and K121 are predicted to bind directly to Stachyflin based on the computer docking model (Figure 3B).
Although mutations of D85H, I91F and L98V were responsible for Stachyflin resistance, their locations on the HA were far from the region of the binding pocket for Stachyflin, leading us to investigate the effect of these mutations. Three-dimensional structure analysis showed that D85 and K83 of another HA2 subunit made a salt bridge, which stabilizes the structure of the HA strongly (Figure 3C) . Amino acid substitution of D85H may abolish the interaction of the salt bridge and cause structural changes and destabilization of the HA, including the binding pocket. In addition, the substitution of aspartic acid to histidine causes structural and pH-dependent instability of the HA because of its electrostatic force . Histidine collects protons around its side chain under low pH conditions, which might cause electrostatic repulsion between the 2 subunits . Additionally, the Stachyflin-resistant virus clone with an amino acid substitution of D85H showed a weak hemolysis of cRBC, which may indicate its structural destabilization; therefore, it is assumed that the structural change of the HA occur by D85H, then Stachyflin is unable to bind to the HA of the D85H mutant.
The mutations of I91F and L98V are distant from the predicted binding pocket for Stachyflin whereas the residue of L98 is reported to be involved in the formation of the binding pocket for TBHQ, which is located upper that for Stachyflin, and have a hydrophobic interaction with TBHQ ; however, our docking simulation showed that Stachyflin was not likely to enter into the cavity for TBHQ (data not shown), suggesting that L98 may not interact directly with Stachyflin. Likewise, I91 may not interact with Stachyflin directly since there is no cavity around I91; therefore, amino acid substitutions of L98V and I91F are postulated to conformational change the structure of the HA, leading to a change in the structure or stability of the binding pocket for Stachyflin.
The frequencies of the Stachyflin-resistant virus clones selected from WSN and PR8 were 10-3.0-10-4.0. RNA viruses lack the ability to detect and repair mistakes during replication and, as a result, the mutation rate can be as high as 1 mutation per each 103-105 bases copied per replication cycle . Based on these data, selection of these resistant variants was not ‘rapid’ but same as other drugs. Resistant variants from Ibaraki and Taiwan were selected in the presence of Stachyflin by 3 passage, indicating that the frequency of the virus clones were much lower than those of WSN and PR8. Structural analysis of the HA indicated that the lower stability of the HA by amino acid substitution for the Stachyflin-resistant virus lead to inefficient virus replication. These results indicate that the resistant virus clones exist in virus population and were isolated only under limited conditions.
In the present study, we found antiviral activity of Stachyflin against A(H1N1)pdm09, H5, including highly pathogenic avian influenza viruses, and H6 viruses, and identified a potential binding pocket for Stachyflin, which differs from that previously proposed [12, 14]. We hereby propose that molecular structures in the potential binding site for Stachyflin depend on the HA of different subtypes, affecting the susceptibility to this compound. Additionally, the present results suggest that further precise analysis of fusion inhibitors reveals their unidentified activities and more suitable docking poses with the HA, contributing to the further development of effective and broad-spectrum fusion inhibitors.
Materials and methods
Stachyflin obtained from Stachybotrys sp. has already been purified and characterized at Discovery Research Laboratories, Shionogi & Co., Ltd . Stachyflin powder was dissolved in dimethyl sulfoxide (DMSO) (Nacalai Tesque, Kyoto, Japan) and was further diluted in each test medium.
Cells and viruses
MDCK cells  were grown in Minimum Essential Medium (MEM) (Nissui Pharmaceutical, Tokyo, Japan) supplemented with 0.3 mg/ml L-glutamine, 5% fetal bovine serum (SAFC Biosciences, Street Lenexa, KS, U.S.A.), 100 U/ml penicillin G, 0.1 mg/ml streptomycin, and 8 μg/ml gentamicin. Human embryonic kidney (293T) cells were maintained in Dulbecco's Modified Eagle's Medium (DMEM) (Life Technologies Japan, Tokyo, Japan) supplemented with 0.3 mg/ml L-glutamine, 10% FBS (Cambrex, Grand Island, NY, U.S.A.), and antibiotics. Both cell lines were maintained at 37°C in a 5% CO2 atmosphere. The influenza virus strains used in the present study are listed in Table 1. All viruses were propagated in the allantoic cavities of 10-day-old embryonated chicken eggs at 35°C for 30–48 h. Before the infectious allantoic fluids were harvested, the eggs were chilled at 4°C overnight, and the harvested allantoic fluids were stored at −80°C.
In vitro antiviral assay
Anti-influenza virus activity of Stachyflin was evaluated by its inhibition of virus-induced CPE in MDCK cells. The virus was inoculated onto confluent monolayers of MDCK cells at the titer of 100 TCID50/ml and virus was adsorbed to the cells at 4°C for 1 h. Unbound viruses were removed by washing the cells with PBS. MEM containing 1.0% DMSO and various concentrations of Stachyflin, from 0.004 to 6.50 μM, were added to the cells and incubated at 35°C. After 72 h, antiviral activity was evaluated by virus-induced CPE and expressed as 50% effective concentration of the compound (EC50).
In vivo antiviral assay
Under anesthesia, 30 μl of WSN or Ibaraki containing 10 50% mouse infectious dose (MID50) was intranasally inoculated into 4-week-old female BALB/c mice (Japan SLC, Shizuoka, Japan). For the anesthesia, a mixture of tiletamine hydrochloride (20 mg/kg) (United States Pharmacopeia, Rockville, MD, U.S.A.), zolazepam hydrochloride (20 mg/kg) (United States Pharmacopeia), and xylazine (20 mg/kg) (Bayer HealthCare, Leverkusen, Germany) was injected intraperitoneally into mice. After the virus inoculation, 100 μl of the solutions containing various concentrations of Stachyflin in polyethylene glycol 400 (Nacalai Tesque) were intraperitoneally administered to each group every 12 h for 72 h. Control mice were injected with only polyethylene glycol 400 after the challenge. At 72 h post-inoculation, mice were euthanized and the lungs were collected for virus recovery. The supernatants of 10% lung homogenates were inoculated onto confluent monolayers of MDCK cells and the virus titers were calculated using the method of Reed and Muench and expressed as TCID50/g of tissue samples . All animal experiments were carried out in self-contained isolator units (Tokiwa Kagaku, Tokyo, Japan) at the BSL-2 or BSL-3 facility of the Graduate School of Veterinary Medicine, Hokkaido University, Japan. The institutional animal care and use committee of the Graduate School of Veterinary Medicine authorized this animal experiment (approval numbers: 10–1052) and all experiments were performed according to the guidelines of this committee.
Selection and characterization of stachyflin-resistant virus clones
WSN, Ibaraki, PR8, and Taiwan were diluted 10-fold series and inoculated on MDCK cells in the presence of various concentrations (0.26, 0.52, 1.30, and 6.50 μM) of Stachyflin. After 72 hours incubation at 35°C, the supernatant of the highest dilution series in the wells in which CPE was observed was collected. EC50 of the viruses were determined as described above and their susceptibilities to the compound were evaluated compared with the parent virus. If the recovered virus did not show increase of EC50, the virus was passaged by same method. For cloning of the viruses, each passaged virus was inoculated on MDCK cells and the cells were then overlaid with MEM containing 1% bacto-agar (Becton, Dickinson, and Company, Franklin Lakes, NJ, U.S.A.) in the presence of 6.50 μM Stachyflin. After 48 h of incubation at 35°C, the cells were stained with 0.014% neutral red (Kanto Chemical, Tokyo, Japan) and the plaques were collected. Individual clones were incubated on MDCK cells in the presence of 6.50 μM Stachyflin at 35°C. After 72 h incubation, each supernatant was collected and stored at −80°C.
Sequence analysis of virus genes
Viral RNA was extracted from the allantoic fluid of embryonated chicken eggs or the supernatant of MDCK cells by TRIzol LS Reagent (Invitrogen, Carlsbad, CA, U.S.A.) and reverse-transcribed with the Uni12 primer  and M-MLV Reverse Transcriptase (Invitrogen). The full-length cDNA of the 8 gene segments was amplified by polymerase chain reaction (PCR) with gene-specific primer sets reported previously  or designed in the present study. The sequences of primers designed in the present study are as follows: PB2-826F: GTTAGGAGAGCAACAGTATCAG, PB2-922R: CAGCTTGCTCTTCTGTTGG, PB1-1240F: GGAATGATGATGGGCATGTT, PB1-1472R: CATCAGACGATTGGAGACCG, PA-723F: CATTGAGGGCAAGCTTTCTC, PA-1110R: CATGTTCTCACCTAATGCCC. Direct sequencing of all 8 gene segments was performed using an auto sequencer, 3500 Genetic Analyzer (Applied Biosystems, Foster City, CA, U.S.A.). To identify amino acid substitutions which should contribute to the susceptibility of viruses, the nucleotide sequences of Stachyflin-resistant virus clones were proofread, and the deduced amino acid sequences were compared with the wild-type virus using GENETYX-WIN version 10 (Genetyx, Tokyo, Japan).
WSN and their mutants were generated by reverse genetics (rg) according to the procedure reports [30, 31], which were named rgWSN, rgR1, rgR2, rgR3, and rgR4, respectively (Table 3). Briefly, viral RNA was extracted and amplified by RT-PCR. The PCR product of each gene segment was cloned into pHW2000 plasmid . Eight genome sets of plasmid were transfected to MDCK and 293T cells and incubated at 37°C for 30 h and then 35°C. After 48 h, rgWSN was collected. All of the collected viruses were propagated in MDCK cells at 35°C and collected after 48 h.
Stachyflin-resistant virus clones with the amino acid substitutions were generated by site-directed mutagenesis as described previously . Briefly, the residue of amino acid substitutions in the HA2 were introduced into the HA genes of WSN using a Quik-Change II site-directed mutagenesis kit (Agilent, Santa Clara, CA, U.S.A.) according to the manufacturer’s instructions. The mutant viruses, rgR1, rgR2, rgR3, and rgR4, were rescued by reverse genetics as described above, and the entire genomes of the 8 gene segments were sequenced to confirm the existence of the introduced mutations and the absence of undesired mutations.
Hemolysis assay was performed as described previously . Briefly, WSN and Stachyflin-resistant virus clones were centrifuged at 25,000 rpm for 1.5 h and the pellets were resuspended in PBS (pH 7.2). Virus concentrates were added to 1 ml of 1% cRBC in saline buffered with 0.1 M citric acid-sodium citrate at a final concentration of 200 HA unit and incubated on ice for 1 h. After the incubation at 37°C for 1 h with mixing every 10 min, the cells were sedimented by centrifugation and the supernatants were measured for hemoglobin at 540 nm.
Protein and ligand structures
Three dimensional models of the H1 HA (WSN) and H5 HA (Ibaraki) molecules were constructed based on the HA crystal structures of PR8 and A/Vietnam/1194/2004 (H5N1), respectively (PDB codes: 1RU7 and 2IBX). After 100 models of the HA trimer were generated using MODELLER 9v6 , a model was chosen by a combination of the MODELLER objective function value and the discrete optimized protein energy (DOPE) statistical potential score . The HA model was evaluated using PROCHECK  and VERIFY3D . The structure of Stachyflin (CID: 493326) was downloaded from the PubChem database.
Molecular docking simulations of the HA and Stachyflin were performed using the Molegro Virtual Docker (MVD) with the default parameter settings .
Hay AJ, Wolstenholme AJ, Skehel JJ, Smith MH: The molecular basis of the specific anti-influenza action of Amantadine. EMBO J 1985, 4: 3021-3024.
Moscona A: Medical management of influenza infection. Annu Rev Med 2008, 59: 397-413. 10.1146/annurev.med.59.061506.213121
Pinto LH, Holsinger LJ, Lamb RA: Influenza virus M2 protein has ion channel activity. Cell 1992, 69: 517-528. 10.1016/0092-8674(92)90452-I
Palese P, Tobita K, Ueda M, Compans RW: Characterization of temperature sensitive influenza virus mutants defective in neuraminidase. Virology 1974, 61: 397-410. 10.1016/0042-6822(74)90276-1
Gubareva LV, Kaiser L, Hayden FG: Influenza virus neuraminidase inhibitors. Lancet 2000, 355: 827-835. 10.1016/S0140-6736(99)11433-8
Kiso M, Mitamura K, Sakai-Tagawa Y, Shiraishi K, Kawakami C, Kimura K, Frederick GH, Sugaya N, Kawaoka Y: Resistant influenza a viruses in children treated with oseltamivir: descriptive study. Lancet 2004, 364: 759-765. 10.1016/S0140-6736(04)16934-1
Shental-Bechor D, Danieli T, Henis Y, Ben-Tal N: Long-range effects on the binding of the influenza HA to receptors are mediated by changes in the stability of a metastable HA conformation. Biochim Biophys Acta 2002, 1565: 81-89. 10.1016/S0005-2736(02)00511-4
Feng F, Miura N, Isoda N, Sakoda Y, Okamatsu M, Kida H, Nishimura S: Novel trivalent anti-influenza reagent. Bioorg Med Chem Lett 2010, 20: 3772-3776. 10.1016/j.bmcl.2010.04.060
O'Keefe BR, Smee DF, Turpin JA, Saucedo CJ, Gustafson KR, Mori T, Blakeslee D, Buckheit R, Boyd MR: Potent anti-influenza activity of cyanovirin-N and interactions with viral hemagglutinin. Antimicrob Agents Chemother 2003, 47: 2518-2525. 10.1128/AAC.47.8.2518-2525.2003
Luo G, Torri A, Harte WE, Danetz S, Cianci C, Tiley L, Day S, Mullaney D, Yu KL, Ouellet C, Dextraze P, Meanwell N, Colonno R, Krystal M: Molecular mechanism underlying the action of a novel fusion inhibitor of influenza a virus. J Virol 1997, 71: 4062-4070.
Vanderlinden E, Göktas F, Cesur Z, Froeyen M, Reed ML, Russell CJ, Cesur N, Naesens L: Novel inhibitors of influenza virus fusion: structure-activity relationship and interaction with the viral hemagglutinin. J Virol 2010, 84: 4277-4288. 10.1128/JVI.02325-09
Plotch SJ, O'Hara B, Morin J, Palant O, LaRocque J, Bloom JD, Lang SA Jr, DiGrandi MJ, Bradley M, Nilakantan R, Gluzman Y: Inhibition of influenza a virus replication by compounds interfering with the fusogenic function of the viral hemagglutinin. J Virol 1999, 73: 140-151.
Yoshimoto J, Kakui M, Iwasaki H, Fujiwara T, Sugimoto H, Hattori N: Identification of a novel HA conformational change inhibitor of human influenza virus. Arch Virol 1999, 144: 865-878. 10.1007/s007050050552
Yoshimoto J, Kakui M, Iwasaki H, Sugimoto H, Fujiwara T, Hattori N: Identification of amino acids of influenza virus HA responsible for resistance to a fusion inhibitor, stachyflin. Microbiol Immunol 2000, 44: 677-685.
Yagi S, Ono J, Yoshimoto J, Sugita K, Hattori N, Fujioka T, Fujiwara T, Sugimoto H, Hirano K, Hashimoto N: Development of anti-influenza virus drugs I: improvement of oral absorption and in vivo anti-influenza activity of stachyflin and its derivatives. Pharm Res 1999, 16: 1041-1046. 10.1023/A:1018983715982
Yoshimoto J, Yagi S, Ono J, Sugita K, Hattori N, Fujita T, Fujiwara T, Sugimoto H, Hashimoto N: Development of anti-influenza drugs: II. Improvement of oral and intranasal absorption and the anti-influenza activity of stachyflin derivatives. J Pharm Pharmacol 2000, 52: 1247-1255. 10.1211/0022357001777225
Kida H, Webster RG, Yanagawa R: Inhibition of virus-induced hemolysis with monoclonal antibodies to different antigenic areas on the hemagglutinin molecule of a/seal/Massachusetts/1/80 (H7N7) influenza virus. Arch Virol 1983, 76: 91-99. 10.1007/BF01311693
Yanagita H, Yamamoto N, Fuji H, Liu X, Ogata M, Yokota M, Takaku H, Hasegawa H, Odagiri T, Tashiro M, Hoshino T: Mechanism of drug resistance of hemagglutinin of influenza virus and potent scaffolds inhibiting its function. ACS Chem Biol 2012, 7: 552-562. 10.1021/cb200332k
Russell RJ, Kerry PS, Stevens DJ, Steinhauer DA, Martin SR, Gamblin SJ, Skehel JJ: Structure of influenza hemagglutinin in complex with an inhibitor of membrane fusion. Proc Natl Acad Sci USA 2008, 105: 17736-17741. 10.1073/pnas.0807142105
Liu S, Li R, Zhang R, Chan CC, Xi B, Zhu Z, Yang J, Poon VK, Zhou J, Chen M: CL-385319 inhibits H5N1 avian influenza a virus infection by blocking viral entry. Eur J Pharmacol 2011, 660: 460-467. 10.1016/j.ejphar.2011.04.013
Govorkova EA, Ilyushina NA, McClaren JL, Naipospos TS, Douangngeun B, Webster RG: Susceptibility of highly pathogenic H5N1 influenza viruses to the neuraminidase inhibitor oseltamivir differs in vitro and in a mouse model. Antimicrob Agents Chemother 2009, 53: 3088-3096. 10.1128/AAC.01667-08
Rachakonda PS, Veit M, Korte T, Ludwig K, Böttcher C, Huang Q, Schmidt MF, Herrmann A: The relevance of salt bridges for the stability of the influenza virus hemagglutinin. FASEB J 2007, 21: 995-1002. 10.1096/fj.06-7052hyp
Mueller DS, Kampmann T, Yennamalli R, Young PR, Kobe B, Mark AE: Histidine protonation and the activation of viral fusion proteins. Biochem Soc Trans 2008, 36: 43-45. 10.1042/BST0360043
Drake JW, Holland JJ: Mutation rates among RNA viruses. Proc Natl Acad Sci USA 1999, 96: 13910-13913. 10.1073/pnas.96.24.13910
Minagawa K, Kouzuki S, Tani H, Ishii K, Tanimoto T, Terui Y, Kamigauchi T: Novel stachyflin derivatives from stachybotrys sp. RF-7260. Fermentation, isolation, structure elucidation and biological activities. J Antibiot (Tokyo) 2002, 55: 239-248. 10.7164/antibiotics.55.239
Tobita K: Permanent canine kidney (MDCK) cells for isolation and plaque assay of influenza B viruses. Med Microbiol Immunol 1975, 162: 23-27. 10.1007/BF02123574
Reed LJ, Muench H: A simple method of estimating fifty per cent endpoints. Am J Hyg 1938, 27: 493-497.
Skehel JJ, Hay AJ: Nucleotide sequences at the 5’ termini of influenza virus RNAs and their transcripts. Nucleic Acids Res 1978, 5: 1207-1219. 10.1093/nar/5.4.1207
Hoffmann E, Neumann G, Hobom G, Webster RG, Kawaoka Y: "Ambisense" approach for the generation of influenza a virus: vRNA and mRNA synthesis from one template. Virology 2000, 267: 310-317. 10.1006/viro.1999.0140
Hoffmann E, Neumann G, Kawaoka Y, Hobom G, Webster RG: A DNA transfection system for generation of influenza a virus from eight plasmids. Proc Natl Acad Sci USA 2000, 97: 6108-6113. 10.1073/pnas.100133697
Manzoor R, Sakoda Y, Nomura N, Tsuda Y, Ozaki H, Okamatsu M, Kida H: PB2 Protein of a highly pathogenic avian influenza virus strain a/chicken/Yamaguchi/7/2004 (H5N1) determines its replication potential in pigs. J Virol 2009, 83: 1572-1578. 10.1128/JVI.01879-08
Tsuda Y, Isoda N, Sakoda Y, Kida H: Factors responsible for plaque formation of a/duck/Siberia/272/1998 (H13N6) influenza virus on MDCK cells. Virus Res 2009, 140: 194-198. 10.1016/j.virusres.2008.12.005
Eswar N, John B, Mirkovic N, Fiser A, Ilyin VA, Pieper U, Stuart AC, Marti-Renom MA, Madhusudhan MS, Yerkovich B, Sali A: Tools for comparative protein structure modeling and analysis. Nucleic Acids Res 2003, 31: 3375-3380. 10.1093/nar/gkg543
Shen MY, Sali A: Statistical potential for assessment and prediction of protein structures. Protein Sci 2006, 15: 2507-2524. 10.1110/ps.062416606
Laskowsk RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst 1993, 26: 283-291. 10.1107/S0021889892009944
Eisenberg D, Lüthy R, Bowie JU: VERIFY3D: assessment of protein models with three-dimensional profiles. Methods Enzymol 1997, 277: 396-404.
Thomsen R, Christensen MH: MolDock: a new technique for high-accuracy molecular docking. J Med Chem 2006, 49: 3315-3321. 10.1021/jm051197e
We thank Dr. J. Yoshimoto and Dr. M. Kobayashi, Shionogi & Co., Ltd., for valuable advice on this study and Dr. K. Minagawa, Shionogi & Co., Ltd., for providing Stachyflin. We also thank our colleagues who have made major contributions to the development of this study. This work was funded by Shionogi & Co., Ltd., Japan. We are grateful for the partial support of the Japan Initiative for Global Research Network on Infectious Diseases (J-GRID) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan and the Global Center of Excellence (GCOE) Program of the Graduate School of Veterinary Medicine, Hokkaido University.
The authors declare that they have no competing interests.
YM drafted the manuscript and carried out in vitro and in vitro experiments in this study. MI performed the computer analyses in this study. NY and MO generated rgWSN and their mutants. TN and RY prepared the compound. MO, YS, KI, and HK participated in the coordination of the study. All authors read and approved the final manuscript.
Yurie Motohashi, Manabu Igarashi contributed equally to this work.
About this article
Cite this article
Motohashi, Y., Igarashi, M., Okamatsu, M. et al. Antiviral activity of stachyflin on influenza A viruses of different hemagglutinin subtypes. Virol J 10, 118 (2013). https://doi.org/10.1186/1743-422X-10-118
- Anti-influenza drug
- HA inhibitor
- Docking model