Nanobiological studies on drug design using molecular mechanic method
Hooria Seyedhosseini Ghaheh1, Maryam Mousavi2, Mahmood Araghi3, Reza Rasoolzadeh1, Zahra Hosseini4
1 Young Researchers and Elites Club, Science and Research Branch, Islamic Azad University, Tehran, Iran
2 Department of Biology, Science and Research Branch, Islamic Azad University, Tehran, Iran
3 Department of Pathology, Amir-Momenin Hospital, Isfahan University of Medical Sciences, Isfahan, Iran
4 Department of Public Health, Isfahan University of Medical Sciences, Isfahan, Iran
|Date of Submission||15-Jan-2013|
|Date of Acceptance||10-Mar-2013|
|Date of Web Publication||07-Oct-2015|
Department of Biology, Science and Research Branch, Islamic Azad University, Tehran
Source of Support: None, Conflict of Interest: None
Background: Influenza H1N1 is very important worldwide and point mutations that occur in the virus gene are a threat for the World Health Organization (WHO) and druggists, since they could make this virus resistant to the existing antibiotics. Influenza epidemics cause severe respiratory illness in 30 to 50 million people and kill 250,000 to 500,000 people worldwide every year. Nowadays, drug design is not done through trial and error because of its cost and waste of time; therefore bioinformatics studies is essential for designing drugs.
Materials and Methods: This paper, infolds a study on binding site of Neuraminidase (NA) enzyme, (that is very important in drug design) in 310K temperature and different dielectrics, for the best drug design. Information of NA enzyme was extracted from Protein Data Bank (PDB) and National Center for Biotechnology Information (NCBI) websites. The new sequences of N1 were downloaded from the NCBI influenza virus sequence database. Drug binding sites were assimilated and homologized modeling using Argus lab 4.0, HyperChem 6.0 and Chem. D3 softwares. Their stability was assessed in different dielectrics and temperatures.
Result: Measurements of potential energy (Kcal/mol) of binding sites of NA in different dielectrics and 310K temperature revealed that at time step size = 0 pSec drug binding sites have maximum energy level and at time step size = 100 pSec have maximum stability and minimum energy.
Conclusions: Drug binding sites are more dependent on dielectric constants rather than on temperature and the optimum dielectric constant is 39/78.
Keywords: Binding site, influenza A (H1N1), molecular mechanic, neuraminidase enzyme
|How to cite this article:|
Ghaheh HS, Mousavi M, Araghi M, Rasoolzadeh R, Hosseini Z. Nanobiological studies on drug design using molecular mechanic method. Adv Biomed Res 2015;4:219
|How to cite this URL:|
Ghaheh HS, Mousavi M, Araghi M, Rasoolzadeh R, Hosseini Z. Nanobiological studies on drug design using molecular mechanic method. Adv Biomed Res [serial online] 2015 [cited 2020 Feb 22];4:219. Available from: http://www.advbiores.net/text.asp?2015/4/1/219/166642
| Introduction|| |
Many worldwide influenza epidemics have been recorded in the history of medicine; the first epidemic occurred in 1918 in Spain in which about 50-100 million people died:  Its reason was the H1N1 virus. The second epidemic occurred in 1957 in Asia, in which about 4 million people died and it was caused by the H2N2 virus. The third epidemic that killed about 1 million people occurred in 1968 in Hong Kong; the cause was the H3N3 virus. The most important matter about these viruses is their nature of being epidemic and pandemic. , In Iran in November of 2009 (Aban of 1388), about 2662 cases of influenza were proved using reverse transcription polymerase chain reaction (RT-PCR),  Nowadays this illness could be transmitted from human to human and it can transfer very quickly through sneezing, coughing and the secretions of the nose and mouth, and touching respiratory droplets.  This virus can stay alive for 2 to 8 hours, so it can transfer through objects, doors and public devices. It has been proven that the risk of transfer is higher in the first 4 hours after exposal.  We can recognize this illness by the use of Nasofarnex and cell culture and also RT-PCR; that would interpret good results  for special patients and groups (high risk) such as: adults, younger children, pregnant women, and diabetics. Influenza A virus belongs to the Orthomyxoviridae family and is classified as three types: A, B and C; Type A is clinically the most important one. The hosts of influenza A are human, swine, hen, horses and birds. ,,
The form of swine influenza virus is mostly spherical or ovoid [Figure 1] and sometimes long, spaghetti- like. They are composed of lipid membrane, glycoprotein molecules, HA and NA, and small amounts of trans-membrane protein, M2. ,, The genome of influenza A viruses consist of eight segments of single-stranded RNA of negative polarity, that code about 11-12 proteins. ,
Proteins that are coded by Influenza A
Protein that is coded by segment 4, Hemagglutinin (HA), and by segment 6, Neuraminidase (NA), and the M2 protein are the most important proteins that employ important roles in drug design (Gong et al., 2007, and Bauer et al., 2009). There are double lipid layer subtype viruses Hemagglutinin (HA, 1-16) and Neuraminidase (NA, 1-9) on the membrane. Humans have Hemagglutinin H1 and H3 and Neuraminidase N1 and N2 (Gong et al.).,,,,,,,
Crystallography and Neuraminidase (NA) enzyme were provided by Laver and his assistants in 1983 [Figure 2]. Neuraminidase (EC: 220.127.116.11) has an important role in the influenza virus and is the major and important goal for controlling the development of this illness. Spanish influenza that occurred in 1918 and its complex with Zanamivir (ZMR) was determined at 1.65-A and 1.45-A resolutions, respectively. The most important characteristic of this virus surface protein is its enzymic role. This protein consists of 4 subunits with asymmetric tetramer architecture which shape a β-sheet and the weight of each subunit is 200 KD.  Considering the sequence of amino acids, this enzyme has been conserved in active sites of different viral species. Babu and colleagues mentioned the conserved binding site theory in 2000. It must be noted that Ca ion enzyme has been observed near the conserved binding sites which its role is still not clear.  As has mentioned before, this enzyme has two chains of α and β. ,,,, Binding sites and active sites are conserved in different species [Figure 2]. The sequence of amino acids could be observed in [Figure 2]. Neuraminidase enzyme is responsible for decomposition or hydrolyzation of Sialic acid and the enzyme of virus enters the host's cell through this way.
|Figure 2: Crystal structure of the NA. The tetramer is composed of four identical monomers. The active site is located on top of the molecule (Xiaojin et al., 2008)|
Click here to view
Zanamivir (Releza) and Oseltamivir (Tamiflu) are sensitive Neuraminidase enzyme inhibitors. ,, Amantadin and Rimantadin are reluctant inhibitors and the M2 protein (responsible for targeting the M2 protein-hydrogen pump channel), is an inhibitor that its consumption has been forbidden by Centers for Disease Control and Prevention (CDC). ,,,
| Materials and methods|| |
In this study, data related to the classification of Neuraminidase (EC18.104.22.168), the weight of subtype 200KD, ,,,,,,,, and succession of amino acid in alpha and beta chains have been obtained from the Protein Data Bank (2010) and NCBI website. Then the whole structure of the enzyme was simulated using Argus Lab 4.0 software.
Since I PDB websites this protein is mentioned along with its inhibitor (Relenza = ZMR), during the simulation the inhibitor was simulated in connection to the active site, which in absence of the Sialic acid inhibitor the host's cell would be set in this site; on the other hand it is the inhibitor's binding site.  (Sialic acid's active site = Inhibitor's binding site)
There are 4 amino acids involved in altering or designing drugs that are: Asp151, Glu276, Arg152 and Arg371. , They are introduced as binding sites. These amino acids were stimulated with the drug (ZMR) in their own spatial situation. The obtained spectrums of NA protein proved that the amino acids had the most interactions with the drug in the 0active site. These amino acids have an important role in drug design that is the main goal of our study in WWthis research on NA binding sites.  Structural data about Neuraminidase enzyme's sequence, the quantity of amino acids, the quantity of water molecules in the proteins structure, for us,  the position of Ca ion, and drug's active sites and binding sites were obtained using Argus Lab software. There are 7009 atoms in this protein, and our study is on 4 amino acids and their inhibitors. In this study the behaviors of enzyme's binding site and active site was observed in aspect of potential energy levels in different dielectric constant and in 310 K temperature. The amino acids that are involved with the drug were studied in different organic solution such as Methanol, aqueous environment like water, and gaseous environmental like vacuum. Transferring the active site of enzyme to different phases in aspect of different dielectric constants was done using the HyperChem 6.0 software. Then the design of the organic environment was conducted through the Chem3D software and the most appropriate dimensions (X, Y, Z) for the box was calculated. The binding sites were also studied in different environments and in 310 K temperature (37DC, body temperature) at the step size of 100 psec.
| Result and Discussion|| |
After conducting the Molecular Mechanic study and gaining the potential energy by Monte Carlo method and studying the amino acids that were involved with the drug in different dielectric constant, the following results were concluded:
Binding site's behavior in vacuum
The Molecular Mechanic computations on the amino acids involved with the drug in vacuum displayed that in the primary stage of step size = 0 ps the account was in its highest level (90 kcal/mol), which decreased to energy level of 51 kcal/mol during time; in a 10 ps period of time, the range of decreasing energy could be observed clearly. The lowest level of energy was observed after 71-80 pSec [Figure 3]a.
|Figure 3: ( a ) Potential energy of binding site NA in vacuum. (b) Potential energy of binding site NA in Water box. (c) Potential energy of binding site NA in methanol|
Click here to view
Binding site's behavior in water box
The Molecular Mechanic computations on the amino acids involved with the drug in aqueous environment showed that the level of the potential energy rose from 980 Kcal/mol to 833 Kcal/mol that can be seen in [Figure 3]b.
Binding sites behavior in Methanol box
The potential energy computations on the amino acids studied in a Methanol box revealed that the level of energy during step sizes 0 ps to 100 ps was decreased from 11034 Kcal/mol to 1310 kcal/mol [Table 1] and [Figure 4]. The Monte Carlo chart also confirms that as time passes the level of energy would be decreased [Figure 3]c.
|Figure 4: (a) Comparison of potential energy levels of binding site NA in water and vacuum. (b) Comparison of potential energy levels of binding site NA in methanol and vacuum. (c) Comparison of potential energy levels of binding site NA in water and methanol|
Click here to view
|Table 1: Comparison of potential energy (Kcal/mol) at 310 K in different dielectrics |
Click here to view
Reports have shown that some drugs are synthesized in the organic environments and it is expected that some of these materials enter the body along with the drugs; therefore this is one of the reasons that the organic environment should be considered as the practical environment of the drug. Also the aqueos environment is the biologic situation of the human body in which protein's structure is hard and strong and it is the best biological solution for correct synthesizing and folding of protein structure.
On the other hand most of drugs need an aqueos environment for effectiveness. Also vacuum environment has basic importance for druggists. The dielectric constant is one of the most effective factors in determining the molecular structure and biological function. Besides statistical studies like Quantum Mechanic and Molecular Mechanic help to develop the drug design, but for good understanding of drug design we should have full and complete data about receptors, target sites and also binding sites. 
| Conclusion|| |
The study of binding sites showed that the system has the highest level of energy and the lowest stability in vacuum environment but during time the level of stability increases and the level of energy decreases, which is caused by the forces from inside the system and would make the amino acids and their aspects to change and find the best spatial conformity; in this situation binding sites have the highest stability, and it's the best condition for them. Vacuum environment has no logical effect on the stability of binding sites of Neuraminidase; therefore it has no effect on designing drug. 
The Molecular Mechanic computations (Monte Carlo) and comparison of the amino acids' stability are the main basics for designing drugs. On the other hand these 4 amino acids should find the best spatial conformity which means the highest stability level or the lowest level of energy. So this study proves that water environment is better for Neuraminidase in comparison with organic and gaseous environments and in the stable biological condition, water part's dielectric constant is fit for biomolecules. Comparing the energy levels in natural temperature and fever temperature, it was revealed that when we face the Influenza enzyme it has no significant change in aspect of EPOT or Potential Energy and we can almost predict that stability or changes of enzyme behavior depends on dielectric constant, not the changes of temperature.
| References|| |
Chen CY, Huang HJ, Tsai FJ. Drug design for Influenza A virus subtype H1N1. J Taiwan Ins Chem Eng 2010;41:8-15.
Anh Nguyen TN, Dao TT, Tung BT, Choi H, Kim E. Influenza A (H1N1) neuraminidase inhibitors from Vitis amurensis. Food Chem 2011;124:437-43.
Wang B, Dwyer DE, Blyth CC, Soedjono M, Shi H, Kesson A, et al
. Detection of the rapid emergence of the H275Y mutation associated with oseltamivir resistance in severe pandemic influenza virus A/H1N1 09 infections. Antiviral Res 2010;87:16-21.
Noorifard M, Khoshdel AR, Hosseini Shokouh SJ. A review in novel pandemic H1N1 A. J Army Univ Med Sci Iran 2009;3:228-40.
Treanor JJ. Influenza viruses including avian influenza and swine influenza. In: Mandell GL, Bennett GE, Dolin R, editors. Mandell, Douglas, and Bennett's Principles and Practice of Infectious Diseases. 7 th
ed. Philadelphia: Churchill Livingstone; 2009. p. 2265-88.
Qi-Shi Dua bd, Wang SQ, Huang RB, Chou KC. Computational 3D structures of drug-targeting proteins in the 2009-H1N1 influenza A virus Chemical Physics. Letters 2010;482:191-5.
Arias CF, Escalera-Zamudio M, Soto-Del Rio ML, Cobian-Guemes AG, Isa P, Lopez S. Molecular anatomy of 2009 influenza virus A (H1N1). Arch Med Res 2009;40:643-54.
Gong J, Xu W, Zhang J. Structure and functions of influenza virus neuraminidase. Curr Med Chem 2007;14:113-22.
Joseph NV. Development of Neuraminidase Inhibitors as Anti-Influenza Virus Drugs. Drug Dev Res 1999;46:176-96.
Xu X, Zhu X, Dwek RA, Stevens J, Wilson IA. Structural characterization of the 1918 influenza virus H1N1 neuraminidase. J Virol 2008;82:10493-501.
Bauer K, Richter M, Wutzler P, Schmidtke M. Different neuraminidase inhibitor susceptibilities of human H1N1, H1N2, and H3N2 influenza A viruses isolated in Germany from 2001 to 2005/2006. Antiviral Res 2009;82:34-41.
Masukawa KM, Kollman PA, Kuntz ID. Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations. J Med Chem 2003;46:5628-37.
Hong AK, Peg MS, von IM. Influenza virus sialidase: Effect of calcium on steady-state kinetic parameters. Biochim Biophys Acta 1991;1077:65-71.
Wang SQ, Du QS, Huang RB, Zhang DW, Chou KC. Insights from investigating the interaction of oseltamivir (Tamiflu) with neuraminidase of the 2009 H1N1 swine flu virus Biochem Biophys Res Commun 2009;386:432-6.
Okomo-Adhiambo M, Nguyen HT, Sleeman K, Sheu TG, Deyde VM, Garten RJ, et al
. Host cell selection of influenza neuraminidase variants: Implications for drug resistance monitoring in A (H1N1) viruses. Antiviral Res 2010;85:381-8.
Baranovich T, Saito R, Suzuki Y, Zaraket H, Dapat C, Caperig-Dapat I, et al
. Emergence of H274Y oseltamivir-resistant A (H1N1) influenza viruses in Japan during the 2008-2009 season. J Clin Virol 2010;47:23-8.
Roschek B Jr, Fink RC, McMichael MD, Li D, Alberte RS. Elderberry flavonoids bind to and prevent H1N1 infection in vitro
. Phytochemistry 2009;70:1255-61.
Mercader AG, Pomilio AB. QSAR study of flavonoids and biflavonoids as influenza H1N1 virus neuraminidase inhibitors. Eur J Med Chem 2010;45:1724-30.
Ives JA, Carr JA, Mendel DB, Tai CY, Lambkin R, Kelly L, et al
. The H274Y mutation in the influenza A/H1N1 neuraminidase active site following oseltamivir phosphate treatment leave virus severely compromised both in vitro
and in vivo
. Antiviral Res 2002;55:307-17.
Monajjemi M, Mahdavian L, Mollaamin F. Characterization of nanocrystalline cylicon germanium film and nanotube in adsoption gas by Monte Carlo and langevin dynamic simulation. Bull Chem Soc Ethiop Bull Chem Soc Ethiop 2008;22:227-86.
Lawrenz M, Wereszczynski J, Amaro R, Walker R, Roitberg A, McCammon JA. Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy. Proteins 2010;78:2523-32.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]