Interaction of Andrographolide with Neuraminidase of H5N1 Virus and Nuclear Factor-kappaB


Indo. J. Phar. Scie. Tech. Vol. 2, No. 1, 1-7 (2015).
Jutti Levita, Monika Novryanti, Masripah Masripah, Driyanti Rahayu, Shelvy E Suherman

Show more

Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Jawa Barat, Indonesia
Korespondensi = (Jutti Levita)


Download citation | PDF (283 kb)


Kata kunci : Andrografolida, neuraminidase, NF-kappaB, studi pemodelan molekul, virus H5N1


Figure 1 Structure of Andrographolide


One of Andrographis paniculata (Burm.F) Nees’s active compounds is andrographolide. It is a major constituent that is extracted from leaves, and it is empirically used as medicine plant.2

Andrographolide’s activities have been reported by some researches, which can be used as antiinflammatory.3 ̶ 5 Not only as antiinflamatory, andrographolide also have antihyperglicemia activity and antioxidant activity.6 Andrographolide can inhibit the proliferation of T-cell and release cytokine as a response to allogenic stimulation.7 In dose of andrographolide was 100 mg/kg bodyweight of rats showed that it’s has antioxidant activity and can decrease lipid peroxidation.6

Nuclear factor-kappa B (NF-kappaB) beginning with its discovery by Baltimore (1986) and continuing to the present, has attracted widespread interest based on its unsual regulation, the kind of stimuli that activate it, biological responses that it controls, the diverse genes and the striking evolutionary conservation of structure and function among family members, and its apparent involvement in a kind of human diseases. Importantly, consistent with the last point, NF-kappaB has been shown to be the target of several antiinflammatory and anticancer drugs.8

NF-kappaB is inducible transcription factor of Rel family and it is sequestered in the cytoplasm by the IkappaB family of proteins. NF-kappaB can exists in several dimeric forms, but the predominant one is p65/p50 heterodimer. Activation of NF-kappaB by wide range of stimuli including viral products and oxidative stress, leads to phosphorilation and proteasome dependent degradation of IkappaB, lead to the release of free NF-kappaB. This free NF-kappaB then binds to its target sites in the DNA to initiate the transcription.9 These kappaB sites are also present in the Long Terminal Repeat (LTR) of HIV-1, and hence NF-kappaB (p50 subunit) binding to LTR-DNA is critical poin in viral replication.10 Targeting direct p50-DNA binding, in this regard, is a recent approach of antiHIV gene expression inhibitors design, which do not have problem of resistance unlike in other antiHIV strategies.11

A palindromic kappaB site is bound the 2.3 Ǻ crystal structure of the transcription factor NF-kappaB p50 homodimer (Figure 2a) reveals that the Rel homology region folds into two distinct domains, similar to those in immunoglobulin super family. The p50 dimer envelopes an undistorted B-DNA helix, and making specific contacts along the 10 base pair kappaB recognition site mainly pass through loops connecting secondary structure of elements in both domains. The carboxy terminal domains form a dimerization of interface between beta-sheets using residues that are strongly conserved in the Rel family.12

Influenza pandemics, a global outbreak of the disease due to viruses with novel antigenic subtypes, have caused death on human population. Pandemics were caused by the hybrid viruses that contained avian and human viral genes. Main contributors of human influenza pandemics are avian influenza viruses.13

Influenza viruses, which belong to the Orthomyxoviridae family, are classified as A, B, and C based on antigenic differences in their nucleoprotein (NP) and matrix (M1) protein. All avian influenza viruses are classified as type A. Further subtyping is based on the antigenicity of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Nowadays, 16 HA and 9 NA subtypes have been identified among influenza A viruses.13


Figure 2 3D Structure of NF-kappaB p50 Homodimer12 (a), H5N1 Neuraminidase16 (b)

Two glycoprotein spikes, HA and NA, and the M2 protein are embedded in the lipid bilayers derived from the host plasma membrane. HA elicits virus neutralizing antibodies. NA is a sialidase that prevents virion aggregation by removing cell and virion surface sialic acid.14

Hence NA is antiviral potential target due to its key role in the infection step of the virus to other cells.15 By inhibiting NA, the virus will not be able to move outward from the cell and finally it will be die.16

The interaction of andrographolide with NF-kappaB and H5N1 neuraminidase (Figure 2b) can be studied by molecular modeling technique, example molecular docking. The focus of molecular docking is to get the optimal conformation, both for the protein and the ligand. In the case of molecular docking, particular arrangement of a ligand and a protein can be defined by a set of values describing the translation, orientation, and conformation of the ligand with respect to the protein. These are the ligand’s state variables and, in the genetic algorithm, each state variable corresponds to a gene. The ligand’s state corresponds to the genotype, the other way the phenotype corresponds to its atomic coordinates. In molecular docking, the fitness is the total interaction energy of the ligand with the protein, and is evaluated using the energy function.17



Instruments that used for this study were personal computer with Mobile Dual Core Intel Core 2 Duo T7250, 1978 MHz (10×198) (, harddisk 250 GB, graphic ATI Mobility Radeon HD 3400 Series (256 MB) and system memory 2044 MB (DDR2-667 DDR2 SDRAM). Softwares used in this study were MOE 2007.09 (Chemical Computing Group,, and HyperChem v7.03 (

The research start with macromolecule preparation. The 3D structures of NF-kappaB crystallized using X-ray diffraction by Gosh G, et al.8 The resolution is 2.30 Å, downloaded from online PDB (Protein Data Bank) code 1NFK, chain A was selected and getting seperated from the macromolecule by using MOE sequence editor. Hydrogens were added to the macromolecule PDB crystal structures followed by calculating the partial charges by using Amber 99. Site finder tool was held to every crystal structure molecule to predict the most likely sites for the ligand-macromolecule interactions. The amino acid sequence in the macromolecules was analized using MOE sequence editor.

3D structures of neuraminidase H5N1 cocrystallized with zanamivir (PDB access code: 2HTQ, resolution 2.2 Å crystallized by Russell RJ).16 Neuraminidase H5N1 cocrystallized with oseltamivir (the PDB access code: 2HT8 resolution 2.6 Å) were downloaded from online PDB.12

The 2D structure of andrographolide was copied from Fujita.18 The 2D structure built using MOE-2007.09.02 builder, 2D structure of aurine tricarboxylic acid and gallic acid copied from Vineet Pande.11

Figure 3 The Electrostatic Map of Andrographolide (a), Pharmacophores of Andrographolide [Showed as Circles] (b)

The 2D structure then built by using MOE-2007.09.02 builder. Energy minimization of each molecule was carried out by using MMFF94x a molecular mechanics force field. Stochastic conformational search with molecular mechanics method was held to find the lowest energy conformer.

The interaction of andrographolide with the macromolecules was analyzed and also compared with the inhibitors. Docking was performed by London dG scoring function and Alpha Triangle placement method.


The andrographolide pharmacophores and also the result electrostatic map of andrographolide showed in Figure 3. The visualizations docking of andrographolide with DBR of p50 NF kappaB and neuraminidase showed in Figure 4. The scoring calculations from andrographolide with p50 NF-kappaB and neuraminidase showed in Table 1.

Figure 4 Visualization of the Docking of Andrographolide with DBR of p50 NF kappaB (a) and Neuraminidase (b)


3D structures of the compounds were built by using MOE-2007.09.02 builder, which constructed the molecules from preexisting fragment libraries. Fragment libraries can be utilized like an electronic 3D tool kit, which is easy to handle. Because of the preoptimized standard geometries of all entries in the fragment pool, the resulting 3D structures already have an acceptable geometry. Mostly only torsion angles have to be cleared to avoid atom overlapping or close van der Waals contacts.

Molecular structures generated using preexisting fragment libraries described previously should always be geometry optimized to find the individual energy minimum state. This is normally done by applying a molecular mechanics method. The molecular mechanics considers the atomic composition of a molecule to be a collection of masses interacting with each other via harmonic forces. As a result of this simplification, molecular mechanics is a relatively fast computational method practicable for small molecules as well as for larger molecules and even systems of oligomolecular. In the calculation, the total energy is minimized with respect to atomic coordinates, where:

Etotal     = Estr + Ebend + Etors + Evdw + Eelec +…

Etotal is total energy of the molecule, Estr is bond stretching energy term, Ebend is angle bending energy term, Etors is torsional energy term, Evdw is van der Waals energy term, and Eelec is electrostatic energy term.

In minimization procedure, molecule structure will be relaxed. Therefore, all of the above corresponding energy terms are altered in a force field optimization.

Molecules are not rigid. At room temperature, the motional energy is large enough to let all atoms in a molecule move permanently. A stochastic conformational search using molecular mechanics method of MOE-2007.09.02 software was held to find the individual energy minimum state of andrographolide and known inhibitors of NF-kappaB conformational search produces eleven conformers. The molecule is flexible due to its many single bonds.

The andrographolide structure shows that the molecule is likely form hydrogen bonds considering its carbonyl, lactone and hydroxyls (see Figure 1). Andographolide is a lipophylic compound which has log P value 2.24 (calculated by using Hyperchem Release7).

The docking was performed by London dG scoring function and Alpha Triangle placement method. Andrographolide can interacts with DBR p50 NF-kappaB via four hydrogen bondings with Asp336, Arg51, Asp206, and Asn244 with its scoring value is -9.99 kcal/mol and the scoring value describes its docking energy.

The same interactions are also showed by aurine tricarboxylic acid and gallic acid, known inhibitors of p50 NF-kappaB. This docking indicates that andrographolide is potential to be developed as inhibitor of

According to both interactions with p50 NF-kappaB and neuraminidase, the pharmacophores of andrographolide are the lactone ring and the hydroxyls attached to C-14 and C-19 as showed in Figure 3b.

Table 1 Scoring Calculations Andrographolide with p50 NF-kappaB and Neuraminidase


From this research can be concluded that andrographolide is able to be proposed as inhibitors of NF-kappaB and H5N1 neuraminidase. This compound interacts with H5N1 via hydrogen bonds and van der Waals interactions. The lactone ring and hydroxyl group are the pharmacophore of andrographolide.



The authors would like to thank DIKTI for defrayed this research via Competitive Grant Project (Penelitian Hibah Bersaing) 2009.


  1. Bimolendu D, Chinmay C, Deepak K, Rupashree S, Rajneeta R, Padma D, and Mitali C. Synthesis, cytotoxicity, and structure–activity relationship (SAR) studies of andrographolide analogues as anti-cancer agent. Bioorganic & Medicinal Chemistry Letters. 2010;20(23):6947–6950.
  2. Subramanian R, Zaini MA, and Amirin In vitro α-glucosidase and α-amylase enzyme inhibitory effects of Andrographis paniculata extract and andrographolide. Acta Biochimica Polonica. 2008;55(2):391–398.
  3. Levita J, Cahya NA, As’ari N, Abdul M, dan      Slamet                       Radioiodinasi andrografolid     dan     biodistribusinya dalam mencit untuk perunut inflamasi. Indonesian             Journal  of           Pharmacy. 2010;21(4):258–265.
  4. Cheng YH, Ming JH, George H, Yi HW, Chi WH, Shiuan WC, Thanasekaran J, Pei TC, Yi HC, and Joen    Andrographolide  enhances nuclear factor- b subunit p65 Ser536 dephosphorylation through activation of protein phosphatase 2A in vascular smooth muscle cells. Journal of Biological    Chemistry.    2011;286(8):5942–5955.
  5. Low M, Cheang SK , Gerald M, Suresh G, and Nikolaus JS. An in vitro study of anti-inflammatory activity of standardised Andrographis paniculata extracts and pure andrographolide. BMC Complementary and Alternative Medicine. 2015;15:18–27.
  6. Ojha SK, Mukesh N, Santosh K, and Arya DS. Antioxidant activity of Andrographis paniculata in ischemic myocardium of rats. Global Journal of Pharmacology. 2009;3(3):154–157.
  7. Michelsen KS, Michelle HW, Brian K, Lisa ST, Deepti D, and Stephan RT. HMPL-004 (Andrographis paniculata extract) prevents development of murine colitis by inhibiting T cell proliferation and TH1/TH17 responses. Inflammatory Bowel Diseases. 2013;19 (1):151–164.
  8. Lee CH, Jeon YT, Kim SH, Song YS. NF-kappaB as a potential molecular target for cancer therapy. BioFactors. 2007;29(1):19–35.
  9. Gilmore TD. Introduction to NF-κB: players, pathways, perspectives. Oncogene. 2006;25:6680–6684.
  10. Williams SA, Lin FC, Kwon H, Carmen MR, Eric V and Warner CG. NF-κB p50 promotes HIV latency through HDAC recruitment and repression of transcriptional initiation. The EMBO Journal. 2006;25(1):139– 149.
  11. Sharma RK, Chopra S, Sharma SD, Pande V, Ramos MJ, Meguro K, Inuoue J, Otsuka M. Biological evaluation, chelation, and molecular modeling studies of novel metal- chelating inhibitors NF-kappaB-DNA binding: structure activity relationship. Journal Med. Chem. 2006;49: 3595-3601
  12. Ghosh G, van Duyne G, Ghosh S, and Sigler. 1995 ( code 1NFK).
  13. Garcia JM, Jimmy CL, Thomas H, Ka TC, Hui LY, Joseph SM, Mark VI, and John MN. Investigation of the binding and cleavage characteristics of N1 neuraminidases from avian, seasonal, and pandemic influenza viruses using saturation transfer difference nuclear magnetic resonance. Influenza and Other Respiratory 2014;8(2):235–242.
  14. Nayak DP, Rilwan AB, Hiroshi Y, Hong ZZ, and Subrata B. Influenza virus morphogenesis and budding. Virus Res. 2009;143(2):147–161.
  15. Jefferson T, Mark J, Peter D, and Chris DM. Neuraminidase inhibitors for preventing and treating influenza in healthy adults: systematic review and meta-analysis. BMJ. 2009;339:5106- 5114.
  16. Russell RJ, Haire LF, Stevens DJ, Collins PJ, Lin YP, Blackburn GM, Hay AJ, Gamblin SJ, Skehel, John J. The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design. Nature. 2006;443:45-49.
  17. Amaro RE, Riccardo B, and J. Andrew An improved relaxed complex scheme for receptor flexibility in computer-aided drug design. Journal Compuert Aided Molecular. 2008;22: 693–705.
  18. Fujita T, Fujitani R, Takeda Y, Takaishi Y, Yamada T, Kido M, and Miura On the Diterpenoids of Andrographis paniculata: X-Ray Crystallographic Analysis of Andrographolide and Structure Determination of New Minor Diterpenoids. Chem. Pham. Bull. 1984;32(6):2117-2125.