Article Information

• PDF (715 kb)
• Full Text with Large Figures

PubMed

• Articles by Panita Decha
• Articles by Supot Hannongbua

Related Articles

• Strategies for developing vaccines against H5N1 influenza A ...

  Strategies for developing vaccines against H5N1 influenza A viruses Trends in Molecular Medicine, Volume 12, Issue 11, 1 November 2006, Pages 506-514 Taisuke Horimoto and Yoshihiro Kawaoka Abstract Recent outbreaks of highly pathogenic avian influenza A virus (H5N1 subtype) infections in poultry and humans (through direct contact with infected birds) have raised concerns that a new influenza pandemic might occur in the near future. Effective vaccines against H5N1 virus are, therefore, urgently needed. Reverse-genetics-based inactivated vaccines have been prepared according to World Health Organization (WHO) recommendations and are now undergoing clinical evaluation in several countries. Here, we review the current strategies for the development of H5N1 influenza vaccines, and future directions for vaccine development. Abstract | Full Text | PDF (716 kb)

• The Genesis of a Pandemic Influenza Virus

  The Genesis of a Pandemic Influenza Virus Cell, Volume 123, Issue 3, 4 November 2005, Pages 368-371 Charles J. Russell and Robert G. Webster Summary Pandemic influenza viruses pose a significant threat to public health worldwide. In a recent paper, now report remarkable similarities between the polymerase genes of the influenza virus that caused the 1918 Spanish influenza pandemic and those of avian influenza viruses. Meanwhile, reporting in show that the reconstructed 1918 Spanish influenza virus kills mice faster than any other influenza virus so far tested. Summary | Full Text | PDF (103 kb)

• The continued pandemic threat posed by avian influenza virus...

  The continued pandemic threat posed by avian influenza viruses in Hong Kong Trends in Microbiology, Volume 10, Issue 7, 1 July 2002, Pages 340-344 Masato Hatta and Yoshihiro Kawaoka Abstract After killing 18 people in 1997, the H5N1 avian influenza virus was thought to have been eliminated, but its descendents have repeatedly been detected in live poultry markets in Hong Kong. Continued surveillance is therefore crucial to prevent further pandemics. Abstract | Full Text | PDF (159 kb)

• …more

Copyright © 2008 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 95, Issue 1, 128-134, 1 July 2008

doi:10.1529/biophysj.107.127456

Biophysical Theory and Modeling

Source of High Pathogenicity of an Avian Influenza Virus H5N1: Why H5 Is Better Cleaved by Furin

Panita Decha^*, Thanyada Rungrotmongkol^*, Pathumwadee Intharathep^*, Maturos Malaisree^*, Ornjira Aruksakunwong^†, Chittima Laohpongspaisan^*, Vudhichai Parasuk^*, Pornthep Sompornpisut^*, Somsak Pianwanit^*, Sirirat Kokpol^* and Supot Hannongbua^*, ,  

^* Department of Chemistry, Faculty of Science, Chulalongkorn University, Patumwan, Bangkok 10330, Thailand
^† Department of Chemistry, Faculty of Science, Rangsit University, Pathumtani 12000, Thailand

Address reprint requests to Supot Hannongbua.

Proteolytic activation of the hemagglutinin (HA) is essential for viral infectivity and for spread of the avian influenza virus through the host's body 1,2,3. This process is determined by a cleavage reaction at the HA cleavage site, a conserved arginine, by host proteases 4,5,6. Insertion of the –RRRKK– residues into the low pathogenic avian influenza (LPAI) cleavage site is known to potentially activate infectivity of viruses, i.e., the LPAI viruses, which then become high pathogenic avian influenza (HPAI) viruses, allowing highly virulent strains to be cleaved by furin (FR), an ubiquitously expressed protease. The goal of this study is to understand why FR cleaves the inserted hemagglutinin strains better than the noninserted strains.

The HA is synthesized initially as an inactive precursor (HA0) that is then activated through a controlled proteolytic cleavage by protease into HA1 and HA2 subunits. HA1 mediates virus binding to host-cell receptors whereas HA2 promotes the release of the viral RNA complexed with polymerase through membrane fusion 1,2,3,4,5,6. Without proteolysis, the fusion peptide cannot occur and therefore the virus is essentially noninfectious.

The HA cleavage site relies on the presence of basic amino acids and relates directly to influenza virus pathogenicity, i.e., the LPAI viruses possess a single arginine at the cleavage site, while the HPAI viruses containing the polybasic insertion upstream of the cleavage site of the H5 subtype have potential to cause devastating pandemics in the future 3,7. In addition, in vitro cleavage of a series of peptide substrates showed that the insertion of the –RRRKK– residues at cleavage site led to an increase in cleavage by FR protease 2,8,9.

FR, a subtilisin-like serine endoprotease, seems to be a highly specific enzyme, cleaving pro-protein precursors at specific consensus sequence –RXK/RR–, usually to produce biologically active products 10,11,12,13,14,15. Based on the x-ray structure of mouse FR complexed with the dec-RVKR-cmk inhibitor, the catalytic ability of FR is considered to have originated from the catalytic triad residues (Ser³⁶⁸, His¹⁹⁴, and Asp¹⁵³) and the oxyanion hole (formed by the carboxamide nitrogens of Asn²⁹⁵ and backbone nitrogen of Ser³⁶⁸ to carbonyl oxygen of the inhibitor/substrate's centered arginine) at the active site 12. The proposed cleavage mechanism is shown in Fig. 116,17.

Display large version of this figure

Figure 1

(A) Proposed cleavage mechanism of HA by furin and definitions of d₁–d₆. (B) Loop of HPAI H5 (ball and stick model) in the electrostatic surface of furin.

This work seeks the source of high pathogenicity in avian influenza A virus subtype H5 in comparison with LPAI subtypes H5 and H3. Three molecular dynamics simulations were carried out for the three complexes, HPH5–FR, LPH5–FR, and LPH3–FR. The investigation was focused to intra- and intermolecular interactions and geometries of the substrate-furin complex potentially involved to the cleavage mechanism.

The crystal structure of FR with bound dec-RVKR-cmk inhibitor 12 used as receptor model for MD simulations was obtained from Protein Data Bank (PDB), code: 1P8J. This study covers the loop of HA substrate features of subtypes H5 (both HPAI and LPAI) and H3 (known as LPAI). The sequence of the HPAI subtype H5 isolated during the 2006 influenza outbreaks in Thailand was taken from Genbank LOCUS ID ABK13784 (A/chicken/Thailand/PC-168/2006(H5N1)) 18. The sequence of the LPAI subtype H5 was obtained from Protein Data Bank (PDB entry code: 1JSM) 19.

The initial model for the HPAI_H5 loop (RERRRKKRGL) was built up using the sequence alignments and the atomic coordinates of the x-ray structure (residues 322–331: NVPEKQTQGL) of the HA0 of H3 (PDB entry code: 1HA0) 20 and dec-RVKR-cmk inhibitor of FR 12 as a template, by using the homology module in the Insight II program (Accelrys, San Diego, CA) 21. The Needleman-Wunsch algorithm 22 was used for pairwise alignment to identify the correspondence region. The coordinates of the structurally conserved regions (SCRs) of the template —in which amino acid sequences in the template and the model are identical—were copied to be those of the model. In addition to the similar SCRs, only the backbone coordinates were replicated and the side chain atoms were added in a manner that preserves the residue type of the model protein.

For LPAI subtype H5, the initial structure of the cleavage loop (NVPQRETRGL) 19 was constructed using the backbone atoms of the HPH5 loop built previously. The side chains were generated by the LEaP module of AMBER 7 23. With this module, the H3 loop above was mutated at residue 329 from glutamine (Q) to arginine (R).

The cleavage loop of HA complexed with FR was generated by the following procedures (Fig. 1): i), all HA heavy atoms of S1–S4 were superimposed with the crystal structure of the dec-RVKR-cmk inhibitor whereas the HA backbone atoms of S5–S8 and S1′–S2′ were superimposed with the HA0 loop of H3. The coordinates of S1–S4 residues were changed according to the –RVKR– coordinates of inhibitor leading to a newly formed conformation of the HA loop; ii), this loop was then placed manually into the FR active site by superimposition between the –R_KR– residues of the HA loop and the dec-RVKR-cmk with a creation of the furin-substrate complex; iii), the complex was minimized by keeping the S1–S4 residues and the FR fixed; and iv), finally, three steps of the restrained MD simulations at 298K were carried out for relaxing the modeled systems with the restrain factors of 10, 5, and 2.5kcalmol⁻¹Å² for 200ps, 200ps, and 200ps, respectively. Therefore, the conformation of each cleavage loop was adapted (from the initial model) to fit better with the FR cavity. The last snapshot obtained from restrained MD procedure was used as the starting structure of the substrate-enzyme complex for the next MD simulations with all atoms allowed to move freely for 2ns.

Three MD simulations for the HA cleavage loop complexed with furin, HPH5–FR, LPH5–FR, and LPH3–FR, were carried out where their initial structures were generated as described above. The simulated systems were neutralized by 6, 11, and 11 Na⁺ ions and solvated by TIP3P water molecules leading to total atoms of 54,500, 56,663, and 56,667 for the HPH5–FR, LPH5–FR, and LPH3–FR complexes, respectively. The dimensions of the obtained simulation boxes for the three systems are 86Å×90Å×88Å. The periodic boundary condition with the NPT ensemble was used. Energy minimization and MD simulations were carried out using the SANDER module of AMBER 7 23 with the Cornell force field 24. A Berendsen coupling time of 0.2ps was used to maintain the temperature and standard pressure of the system 25. The SHAKE algorithm 26 was applied to constrain all bonds involving hydrogens. The simulation time step of 2 fs was used. All MD simulations were run with a 12Å residue-based cutoff for nonbonded interactions and the particle mesh Ewald method was used for an adequate treatment of long-range electrostatic interactions 27. The MD simulations were carried out for 2.0ns where the production phase starts from 0.75ns to 2.0ns. The convergences of energies, temperature, pressure, and global root mean-square displacement (RMSD) were used to verify the stability of the systems. The MD trajectories were collected every 0.2ps. Analysis of all MD trajectories, i.e., RMSD, distances, hydrogen bonds, etc., were carried out using the Ptraj, CARNAL, and MMGB/SA modules of AMBER 7.

To monitor the stability of the three simulated systems, HPH5–FR, LPH5–FR, and LPH3–FR, RMSD of the structures obtained during the 2ns MD simulations relative to the initial structure of all heavy atoms in the substrate-enzyme complex and the 10 residues of HA substrate were evaluated and plotted in Fig. 2. It can be seen from the plot of RMSD plots versus simulation time that all three systems were found to reach equilibrium at 0.75ns.

Display large version of this figure

Figure 2

RMSDs relative to the initial structure for all heavy atoms of the substrate-furin complexes (black) and the HA loops (gray) for the three systems studied: HPH5–FR, LPH5–FR, and LPH3–FR.

The main results refer to the schematic representation shown in Fig. 1 where the conformation structures of the complexes are described in terms of d₁–d₆ (Fig. 3) in which the reaction coordinates involved in the acylation in three regions, the catalytic triad of FR, nucleophilic attack, and oxyanion hole, were described by the d₁–d₃, d₄, and d₅–d₆, respectively (Fig. 1). Accordingly, percentage H-bond occupations were evaluated using the CARNAL module in AMBER7 based on the following criteria: i), proton donor-acceptor distance ≤3.5Å; and ii), donor-H-acceptor bond angle ≥120°. The results were summarized in Table 1.

Display large version of this figure

Figure 3

Distributions of the d₁–d₆ distances defined in Fig. 1 for the three simulated systems, sampling from 0.75 to 2.0ns in MD simulations.

At the catalytic triad (Asp¹⁵³, His¹⁹⁴, and Ser³⁶⁸) region of the HPH5–FR, the simulated distances of d₁=d₂=d₃=2.85Å (Fig. 3) are comparable to those found experimentally for the inhibitor-furin complex of 2.53Å, 3.10Å, and 2.95Å, respectively 12. This indicates that the three catalytic residues, Asp¹⁵³, His¹⁹⁴, and Ser³⁶⁸, are in the configuration suitable for initiating the nucleophilic reaction. The situation is different for the bound LPAI, H5, and H3, hemagglutinin loops, where only d₃ of both systems and d₁ of LPH5–FR fall within the range whereby the reaction can take place.

The nucleophilic attack is determined by the distance between the O³-hydroxyl oxygen of Ser³⁶⁸ and the carbonyl carbon of the S1-Arg, d₄ in Fig. 1. The d₄ bond-making distances were found to be at 3.05Å, 3.45Å, and 3.75Å for the HPH5–FR, LPH5–FR, and LPH3–FR systems, respectively (Fig. 3). In addition to a short d₄, the detected sharp and narrow peak was found only in the HPH5–FR complex, signifying the rigidity of the complex that thus serves as a more appropriate configuration for the nucleophilic attack. In contrast, for the other two complexes the d₄ peak shows a broad distribution and takes place at significantly longer distances.

At the oxyanion hole region of the FR active site (d₅ and d₆ in Fig. 1), the hole is generally formed by the backbone nitrogen of Ser³⁶⁸ and the carboxamide nitrogen of Asn²⁹⁵, which specifically binds with the carbonyl group of the S1 reacting residue 12. In Fig. 3, the d₅ distance of 2.95Å for the HPH5–FR (solid line) is slightly shorter than that of 3.15Å for the LPH5–FR (dashed line) with the presence of a strong hydrogen bond to Ser³⁶⁸, 100% and 86% occupations (d₅ in Table 1), respectively. The percentage occupation of 15% for the LPH3–FR confirms a very weak interaction and highly flexible of the complex. This conclusion was also confirmed by a broad distribution of d₅ distance varying from 3.00Å to 5.85Å (Fig. 3). In terms of d₆ distance, the distribution plot for the HPH5–FR complex shows two apparent peaks centered at 2.85Å and 5.45Å where the sharper first peak indicated a preferential interaction. For LPAI systems, the d₆ distance is detected at 3.05Å in H3 whereas this distance is significantly longer in H5. This means that interaction between S1 residue and Asn²⁹⁵ was almost lost in the LPH5–FR system. The observed intermolecular distances agree very well with the percentage pair of hydrogen bond through the d₅ and d₆ in which those values (d₅=100% and d₆=34%) for the HPH5–FR are higher than those (d₅=86% and d₆=6%) for the LPH5–FR. On the contrary, bonding via these two distances almost disappear for the LPH3–FR complex. Note that the experimental d₅ and d₆ distances for the inhibitor-furin complex are 3.34Å and 2.78Å, respectively 12.

To assess to molecular interaction at the cleavage loop of hemagglutinin, percentage, and number of hydrogen bond between each HA residue and the active site residues of FR, were evaluated and plotted in Fig. 4. The criteria used are the same as those mentioned earlier. Here, the experimentally detected hydrogen bonds for the inhibitor-furin complex 12 were also given for comparison, the residue with a box around the amino acid label. Interest was focused on S1-Arg of HA and FR's flat groove—a cavity formed by residues Pro²⁵⁶, Asp²⁵⁸, Ala²⁹², and Asp³⁰⁶. The latter was found to play an important role in facilitating the cleavage reaction by forming hydrogen bonds to the guanidinium side chain (Fig. 1) of P1-Arg of the inhibitor 12 (Si denotes a residue of the substrate whereas Pi represents one of the inhibitors, where i=1, 2, 3, …) 28. Strong hydrogen bonds between the four flat groove residues and S1-Arg were highly detected for the HPH5–FR complex (S1 in Fig. 4). The cleavage reaction in the HPH5–FR system was also promoted by the other hydrogen bonds from S1-Arg to Ser³⁶⁸ (two bonds, one of which is d₅ in Fig. 1) and Ser²⁵³ (one bond) of FR with almost 100% occupation. These observations are strongly supported by the experimentally detected hydrogen bonds also shown in Fig. 1. This is in contrast with LPH5–FR in which no hydrogen bond with the flat groove was detected and only a single bond with the catalytic Ser³⁸⁶ (d₅ in Fig. 1) was formed. The LPH3–FR loop shows much less interaction with the FR residues, forming no hydrogen with either the flat groove or the catalytic Ser³⁶⁸ residue.

Display large version of this figure

Figure 4

Percent occupations of hydrogen bonds between furin and the ten HA residues (defined in Fig. 1) where the residues with a box around the label represent experimentally detected bonds for the inhibitor-furin complex.

Some comments could be made concerning the two hydrogen bonds formed between the reacting S1-Arg and the catalytic Ser³⁸⁶ of HPH5–FR (Fig. 4). Besides the bond defined by d₅, another bond is S³⁶⁸_O³…NH_S1-Arg (71% occupation). This hydrogen bond that has never been found in the oxyanion hole region of any other HA subtype or in the inhibitor-furin complex 12, could be a reason for the shortening and rigidity (narrow and sharp peaks shown in Fig. 3) of the intermolecular d₄ distance, which could directly facilitate the nucleophilic reaction in HPH5–FR.

Considering the role of the –RRRKK– insertion (Fig. 1), more hydrogen bonds and a higher percentage occupation between the S2-S6 residues of HA and the surrounding residues of FR (Fig. 4) were found for HPH5–FR in comparison with the two LPAI systems. This means that the –RRRKK– insertion can directly help to hold the substrate in place.

Taking into account all the hydrogen bond data given above, the hemagglutinin loop of HPH5–FR is observed to bind much more tightly into the catalytic site of FR than the LPH5–FR and LPH3–FR systems.

To seek the fundamental basis of substrate-furin interactions, the interaction energy between individual FR residue and the three HA loops and vice versa were calculated using the decomposition energy module of AMBER 7. Plots of the decomposition energies (DC) of the HA residues, as well as the selected FR residues that are located within 5Å from the 10 residues of HA loop, were shown in Fig. 5. The overall DCs of FR residues were observed in the following order: HPH5–FR ≫ LPH5–FR ∼ LPH3–FR (Figure 5B). For the HA loops, the DCs for the residues of HPH5–FR are significantly lower than for the corresponding residues of LPH5–FR and LPH3–FR (Figure 5A). This is especially true for the S1, S4, and S6 residues. The calculated results are consistent with experimental results for the inhibitor-furin complex, which stated that P1, P4, and P6 residues of the inhibitor were observed to interact strongly with FR 11,29. In addition, these two P4 and P6 residues are required to dramatically increase the cleavage reaction 30. The DCs for the FR and HA residues agree very well with the hydrogen bond data shown (Fig. 4) and discussed above.

Display large version of this figure

Figure 5

Decomposition energies (DC) of (A) 10 HA residues (defined in Fig. 1), and (B) selected furin residues within 5Å of those 10HA residues.

Hydrophilicity that is determined by the ligand solution in the cavity region of enzyme-substrate complex is known to play an essential role in the catalytic mechanism of biological systems. To seek this information, the plot of the radial distribution functions (RDFs, g_xy(r))—the probability of finding the particle of type y in a spherical radii, r, around the particle of type x—was evaluated. Here, the RDFs from the acceptor atoms of FR pocket (O¹, O²:Asp¹⁵³; N¹, N²:His¹⁹⁴; N⁴:Asn²⁹⁵, and N³, O³:Ser³⁶⁸) and the reactive HA residue (O^γ:S1-Arg) (see label in Fig. 1), to oxygen atom of water was calculated. The results for the three simulated systems are shown in Fig. 6 together with the corresponding running integration numbers up to the first minimum of the RDF plot.

Display large version of this figure

Figure 6

Radial distribution functions, g(r), from acceptor atoms of furin catalytic residues and HA loops (see Fig. 1 for definition) to oxygen atoms of water molecules for the three simulated systems, HPH5–FR, LPH5–FR, and LPH3–FR complexes.

The RDFs for all oxygen atoms, O¹, O², O³, and O^γ, of the LPH3–FR complex show sharp and narrow first peaks at ∼2.9Å with the integration number (coordination number, CN) of 1.5, 1, 1 and 2 water molecules, respectively (Fig. 6). The CN denotes the number of water molecules that locate in the first hydration shell around the central atom. This is in contrast for the LPH5–FR complex where only O² can be accessed by water molecule whereas no water was detected around the four oxygen atoms of HPH5–FR.

The RDFs centered on N¹ and N⁴ atoms of the three systems are almost similar, the first broad peaks at ∼3.8Å and sharp peaks at ∼3Å, respectively. Note that, N³ atom of only LPH3–FR was solvated by water molecule. Interest is focused to N² atom of His¹⁹⁴ where the catalytic water was observed and proposed to play role in the reaction mechanism of serine protease 14. In good agreement with our simulated results, one water (CN=1) was detected to preferentially coordinate (first peak centered at ∼2.9Å) to N² atom of only for the HPH5–FR complex. This is not the case for the other two low-pathogenic systems, LPH3–FR and LPH5–FR, where less water, 0.5 and 0.2 were, respectively, was found (Fig. 6, dashed and dotted lines). In addition, distance to the LPH3-FR water molecule of ∼3.45Å is significantly longer than that of the HPH5–FR complex.

Taking into account the solvation data shown above, the number of overall water molecules around the selected central atoms are in the following order: LPH3–FR>LPH5–FR>HPH5–FR. The obtained data indicates that more water molecules can take place in the catalytic pocket of FR in the complex whose substrate does not fit well with the FR active site, giving a subsequent result of more cavities or less binding between substrate and enzyme. This suggestion agrees very well with the conformational structures of the complexes (Fig. 3) and the hydrogen bond data (Fig. 4) shown and discussed previously.

The technique of molecular dynamics simulations has been applied to look for detailed information focused on the intra- and intermolecular interactions and geometries of the three substrate-furin complexes, HPH5–FR, LPH5–FR, and LPH3–FR. The results obtained from this study indicate that the –RRRKK– insertion in the HPH5–FR (two arginines at S4 and S6 positions in particular) helps to directly hold the HA cleavage loop in place by forming many strong hydrogen bonds between the residues of HA and FR. This data is well supported by the corresponding strong-decomposition interaction energies and leads consequently to the accumulation of less water molecules accessible to the FR cavity. An active conformation of the HPH5–FR complex suitable for the acylation reaction by FR was formed. In conclusion, the results from this study provide a clear answer to the question of why HPAI hemagglutinin H5 is better cleaved by FR than the other two HA subtypes.