- Open Access
Agonistic effect of polyunsaturated fatty acids (PUFAs) and its metabolites on brain-derived neurotrophic factor (BDNF) through molecular docking simulation
Lipids in Health and Diseasevolume 11, Article number: 109 (2012)
Brain-derived neurotrophic factor (BDNF) is a potent neurotrophic factor that is implicated in the regulation of food intake and body weight. Polyunsaturated fatty acids (PUFAs) localised in cell membranes have been shown to alter the levels of BDNF in the brain, suggesting that PUFAs and BDNF could have physical interaction with each other. To decipher the molecular mechanism through which PUFAs modulates BDNF’s activity, molecular docking was performed for BDNF with PUFAs and its metabolites, with 4-Methyl Catechol as a control.
Inferring from molecular docking studies, lipoxin A4 (LXA4), and a known anti-inflammatory bioactive metabolite derived from PUFAs, with a binding energy of −3.98 Kcal/mol and dissociation constant of 1.2mM showed highest binding affinity for BDNF in comparison to other PUFAs and metabolites considered in the study. Further, the residues Lys 18, Thr 20, Ala 21, Val 22, Phe 46, Glu 48, Lys 50, Lys 58, Thr 75, Gln 77, Arg 97 and Ile 98 form hot point motif, which on interaction enhances BDNF’s function.
These results suggest that PUFAs and their metabolites especially, LXA4, modulate insulin resistance by establishing a physical interaction with BDNF. Similar interaction(s) was noted between BDNF and resolvins and protectins but were of lesser intensity compared to LXA4.
Neurotrophins protect and promote the growth and development of nerves . They have potent effects on neuronal differentiation, survival, neurite outgrowth, synaptic formation, and plasticity [2, 3]. Of all, BDNF and its receptor TrkB are the most abundantly expressed in hippocampus and cerebral cortex region of the brain . BDNF plays a crucial role in controlling body weight and energy homeostasis . Increased levels of BDNF reduce food consumption and maintain energy balance . Both central and peripheral administration of BDNF decreased food intake, increased energy expenditure and ameliorated hyperinsulinaemia and hyperglycaemia in diabetic db/db mice [7–10].
TrkB activation by BDNF is essential for appetite regulation and energy homeostasis . BDNF mediated TrkB signalling is an important downstream target for MC4R-mediated signalling which participates in the regulation of energy balance and feeding behaviour . Mutations in the TrkB receptor resulted in hyperphagia and morbid obesity in humans and rodents. Conversely, peripheral or central stimulation of TrkB by its natural ligands BDNF or NT4 reduced body weight and food intake in mice, supporting the idea that TrkB is a key anorexigenic signal downstream of the melanocortin-4 receptor (MC4R) system .
Polyunsaturated fatty acids such as linoleic acid (LA), arachidonic acid (AA), α-linolenic acid (ALA), eicosapentaneoic acid (EPA) and docosahexaenoic acid are present in almost all cell membranes and regulate their (cell membrane) fluidity, receptor number and affinity of receptor to their respective hormones, peptides and growth factors [12–15]. In this context, it is interesting to note that dietary supplementation with omega-3 PUFAs reportedly normalizes BDNF levels which are reduced following brain injury . Omega-3 enriched dietary supplement provides protection against reduced plasticity and impaired learning ability associated with brain injury in rats. High-fat diets, particularly those rich in saturated fats, adversely affect insulin action and alter homeostasis model assessment (HOMA) , which can be prevented by n-3 PUFAs rich in eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) by reducing insulin resistance .
Furthermore, PUFAs are needed for synapse formation, neurite outgrowth and also have neuroprotective actions [19, 20]. Nerve growth cones are highly enriched with AA-releasing phospholipases, which have been implicated in neurite outgrowth. Cell membrane expansion occurs through the fusion of transport organelles with plasma membrane, and syntaxin 3, a plasma membrane protein that is important in the growth of neurites, is a direct target for AA, DHA and other PUFAs. These (AA, DHA, and other PUFAs) polyunsaturated fatty acids but not saturated and monounsaturated fatty acids activate syntaxin 3.
SNAP25 (synaptosomal-associated protein of 25 kDa), a syntaxin partner implicated in neurite outgrowth, interacts with syntaxin 3 only in the presence of AA that allowed the formation of the binary syntaxin 3-SNAP 25 complex. AA stimulated syntaxin 3 to form the ternary SNARE complex (soluble N-ethylmaleimide-sensitive factor attachment protein receptor), which is needed for the fusion of plasmalemmal precursor vesicles into the cell surface membrane that leads to membrane fusion . Thus, AA and DHA change the α-helical syntaxin structure to expose SNARE motif for immediate SNAP 25 engagement that facilitates neurite outgrowth. This could be one important mechanism by which AA, EPA and DHA are able to enhance neurogenesis and bring about their beneficial actions. In this context, it is interesting to note that BDNF and TrkB are expressed in the retinal ganglion cell layer and are able to influences the morphological differentiation of the cells and act as cell survival factors [22, 23].
Since both BDNF and n-3 PUFAs levels are known to be altered in diabetes and several other neurological conditions, our current study aims to investigate the effect of PUFAs and their metabolites on BDNF binding towards the receptor, tropomyosin-related kinase-B (TrkB) and its function through molecular docking simulation. This will provide an insight on the bioactive conformation of the protein through simulation studies and the possible high affinity binding orientations of the ligands, thus paving way for understanding the modulating mechanism of PUFAs and its metabolites on BDNF’s affinity.
Molecular Dynamics simulation and Optimization
BDNF is a 247 amino acid protein, secreted in the extracellular space. The 3D (three dimensional) structure of the heterodimeric form of BDNF with Neurotrophin 3 [PDBID: 1BND] has been elucidated using X-ray crystallography at 2.30Å resolution . To optimize and understand the conformational changes of the monomeric form of BDNF in the presence of explicit solvent, molecular dynamics simulation was performed using GROMACS (Groningen Machine for Chemical Simulations) 4.3.1 package . The model system was solvated with three-point transferable intermolecular potential (TIP3P) molecules  in a cubic water box with the periodic boundary conditions set. To neutralize the system, eight CL − ions were added and subjected to energy minimization through Optimized Potentials for Liquid Simulations (OPLS) force field using the steepest descent integrator. After reaching the convergence limit, the BDNF in the system was subjected to position-restrains by keeping number of particles (N), pressure (P) and temperature (T) constant (NPT). Using Particle Mesh Ewald (PME) electrostatics method under NPT conditions, final MD simulation of 5,000,000 was performed for 10,000 ps (10 ns). Finally, the optimized structure was utilized for further docking studies.
Small molecule optimization
3D Structural co-ordinates of 4-MC [CID 9958], LXA4 [CID 5280914], AA [CID 444899], EPA [CID 446284], DHA [CID 445580], RVE1 [CID 10473088] and NPD1 [CID 16042541] were retrieved from Pubchem database  and were geometry optimized using PRODRG server with full charges .
Molecular Docking simulation of the ligands with the monomeric form of BDNF was carried out using AutoDock 4.2 in combination with Lamarckian genetic algorithm (LGA) to find the optimal conformation of the ligand. The optimized 3D structure of BDNF was prepared by removing water molecules and CL- ions; polar hydrogens were added. Further, Kollman united atom charges and Gasteiger charges were added to the receptor and ligand, respectively. The atomic solvation parameters were assigned using the ADDSOL utility of Auto Dock 4.2. Flexibility of the ligands was assigned based on its torsional degrees of freedom through Autotors option, with the protein torsional angles were kept fixed throughout the process of docking simulation [29–31].
Grid box covering the complete surface of BDNF was constructed and used for molecular docking with a grid spacing of 0.514 Å and the number of points in each grid dimension set to 126 × 126 × 126. Further, Grid maps were generated for each atom within the ligands using Autogrid. Using LGA, molecular docking was carried out with default parameters, except for the number of GA runs which was set to 100, Maximum number of generations was set to 500 and maximum number of evaluations set to 250 000 .
Further, cluster analysis was performed by sorting the binding modes based on the binding energy within the specified RMSD threshold, which was set to 2Å. The optimal binding orientation of the ligand is the one with the lowest binding energies picked up from the cluster, which has maximum number of similar conformations. However, agonistic mode was screened based on the binding energy, dissociation constant and through interaction studies . Hot point interactions mediating BDNF - ligand complex were analyzed and visualized using PDBsum  and PyMOL (http://www.pymol.org), respectively.
Molecular Dynamics simulation of BDNF
MD simulation of the monomeric form of BDNF was analyzed to gain a better insight into the stability and dynamics nature of the molecule. Initial analysis was performed to ascertain the quality of MD simulation performed. The quality checks involves convergence of thermodynamic parameters: temperature, potential and kinetic energy. Analyzing the graphs, convergence in all these parameters were observed, which is suggestive of system’s equilibrium.
Further, Systems equilibrium in terms of structure was inspected by comparing the Root Mean Square Deviation (RMSD) of the C-alpha with respect to the starting structure. Convergence of RMSD was inferred from the graph towards the process of simulation (Figure 1a). Root Mean Square Fluctuations (RMSF) was examined to infer the flexibility in the residues during the course of simulation (Figure 1b). These structural level analyses were indicative of the system reaching equilibrium. And, the final optimized monomeric form of BDNF was taken for further docking analysis.
Molecular Docking studies with a known reference enhancer, 4-methyl catechol
Atsumi Nitta et al. proposed the stimulating effects of 4-MC on the synthesis of BDNF in neurons and astrocyte cells . However, the agonistic binding mode of 4-MC on BDNF has not been studied. Therefore, molecular docking studies have been performed for BDNF with 4-MC. Among 69 conformers in Cluster1, the one with a lowest binding energies of −4.36Kcal/mol and dissociation constant of 635.36 μM was considered for the further analysis (Table 1). Inferring from interaction analysis, Tyr47, Glu48 and Asp99 showed hydrogen bond interactions; wherein Val 13, Lys 18, Ala 21, Val 22, Phe 46, Arg 97 and Ile 98 showed non-bonded interactions with 4-MC (Figure 2) (Table 1) and it could be a probable hot point residues which on interaction stimulates BDNF activity.
Molecular Docking studies with PUFAs
The binding conformation of PUFAs such as DHA, EPA and AA has been inferred through molecular docking studies. Among the PUFAs considered in the study, EPA, an omega-3 fatty acid with a binding energy of −3.67 Kcal/mol and dissociation constant of 2.06mM showed highest binding affinity for BDNF by forming a hydrogen bond interaction with Lys 18 and Lys 50; non – bonded interactions with Thr 20, Ala 21, Val 22, Phe 46, Glu 48, Arg 97, Ile 98 and Asp 99. However, DHA and AA also showed relatively higher binding affinity than the metabolites with a binding energy of −2.9 Kcal/mol and −2.89 Kcal/mol, respectively (Table 1) (Figure 2).
Docking interaction analysis of BDNF with metabolites REV1, NPD1 and LXA4
Molecular docking studies were performed to infer the mechanism through which metabolites of PUFAs mediates BDNFs activity. LXA4, an endogenously synthesized nonclassic eicosanoid, showed highest binding affinity for BDNF in comparison to the PUFAs and metabolites considered in the study with a binding energy of −3.98 Kcal/mol and dissociation constant of 1.2 mM. And, LXA4 established a network of hydrogen bond interactions with Cys 51 and Lys 58; hydrophobic interaction with Thr 20, Ala 21, Phe 46, Lys 50, Cys 73, Arg 74, Thr 75, Gln 77,Arg 97, Ile 98 and Asp 99 of BDNF. Wherein, NPD1 and REV1 showed comparatively lower binding affinity for BDNF with the binding energy of −2.51 Kcal/mol and −2.1 Kcal/mol, respectively (Table 1) (Figure 2).
Inspecting at the interacting residues, consensus in the residues mediating interaction with PUFAs and its metabolites was observed. Further, same residues were found to establish interaction with 4-MC and, these residues involve Lys 18, Thr 20, Ala 21, Val 22, Phe 46, Glu 48, Lys 50, Lys 58, Thr 75, Gln 77, Arg 97 and Ile 98. Thus, suggestive of these residues forming a probable hot point motif mediating the BDNFs activity (Figure 3).
Molecular docking studies with the optimized model of BDNF elicited the potential binding modes of agonists. The active binding mode of the catechol compound, 4-MC has been inferred through molecular docking studies. Thus confirming the 4-MC’s enhancing effect on BDNF mRNA expression and secretion. However, the inability of 4-MC to cross the Blood Brain Barrier limits it therapeutic efficacy, which in turn requires novel agonist for BDNF neurotrophic activity.
In case of endogenous agonists considered in the study, LXA4 showed highest binding affinity for BDNF by forming network of bonded and non- bonded interactions with Thr 20, Ala 21, Phe 46, Lys 50, Cys 51, Lys 58, Cys 73, Arg 74, Thr 75, Gln 77, Arg 97, Ile 98 and Asp 99 with a significant binding energy and disassociation constant.
EPA showed highest binding affinity with a binding energy of −3.67Kcal/mol and disassociation of 2.06mM among other PUFAs studied. Molecular interaction analysis showed polar and non-polar contacts with the binding cavity, which is shared by LXA4 involving the residues: Lys 18, Thr 20, Ala 21, Val 22, Phe 46, Glu 48, Lys 50, Arg 97, Ile 98 and Asp 99. Wherein, NPD1 showed lowest binding affinity for BDNF in comparison to other metabolites considered in the study.
Further, the present study reveals the common binding mode of PUFAs, metabolites and 4-MC to BDNF, inferred through molecular docking studies (Figure 3). Hence, the region of BDNF harbouring this common binding mode may be the probable site where small molecular enhancers bind, and thereby promoting neurotrophic activity.
Moreover, this site also does not interfere with the TrkB binding pockets in BDNF, which spans at the n-terminal region . Hence, it could be hypothesized that the binding of PUFAs to BDNF shall exert allosteric effect on the TrkB binding pocket, thereby modulating the activity cascade. However, in vitro studies need to be performed to ascertain the in silico predictions.
The results of the present study are consistent with the hypothesis that PUFAs and their metabolites such as lipoxins, resolvins and protectins interact with BDNF and possibly, other neurotrophic factors and thus bring about their beneficial actions both in diabetes mellitus and neurological conditions such as depression [37–39]. Alternatively, the beneficial actions of BDNF in these clinical conditions could also be attributed to its interaction with PUFAs. Thus, a better understanding of the close interaction(s) between BDNF and PUFAs and their metabolites may pave way for the development of newer therapeutic strategies in diabetes mellitus, depression and other clinical conditions in which they are believed to play a significant role.
Brain Derived Neurotrophic Factor
Poly Unsaturated Fatty Acids
Tropomyosin-Related Kinase B
homeostasis model assessment
Groningen Machine for Chemical Simulations
Three-Point Transferable Intermolecular Potential
Optimized Potentials for Liquid Simulations
Particle Mesh Ewald
Lamarckian Genetic Algorithm
Root Mean Square Deviation.
Lewin GR, Barde YA: Physiology of the neurotrophins. Annu Rev Neurosci. 1996, 19: 289-317. 10.1146/annurev.ne.19.030196.001445
Hennigan A, O'Callaghan RM, Kelly AM: Neurotrophins and their receptors: roles in plasticity, neurodegeneration and neuroprotection. Biochem Soc Trans. 2007, 35: 424-427. 10.1042/BST0350424
Rosenberg SS, Ng BK, Chan JR: The quest for remyelination: a new role for neurotrophins and their receptors. Brain Pathol. 2006, 16: 288-294. 10.1111/j.1750-3639.2006.00035.x
Huang EJ, Reichardt LF: Neurotrophins: roles in neuronal development and function. Annu Rev Neurosci. 2001, 24: 677-736. 10.1146/annurev.neuro.24.1.677
Wisse BE, Schwartz MW: The skinny on neurotrophins. Nat Neurosci. 2003, 6: 655-656. 10.1038/nn0703-655
Xu B, Goulding EH, Zang K, Cepoi D, Cone RD, Jones KR, Tecott LH, Reichardt LH: Brain-derived neurotrophic factor regulates energy balance downstream of melanocortin-4 receptor. Nat Neurosci. 2003, 6: 736-742. 10.1038/nn1073
Ono M, Ichihara J, Nonomura T, Itakura Y, Taiji M, Nakayama C, Noguchi H: Brain-derived neurotrophic factor reduces blood glucose level in obese diabetic mice but not in normal mice. Biochem Biophys Res Commun. 1997, 238: 633-637. 10.1006/bbrc.1997.7220
Nakagawa T, Tsuchida A, Itakura Y, Nonomura T, Ono M, Hirota F, Inoue T, Nakayama C, Taiji M, Noguchi H: Brain-derived neurotrophic factor regulates glucose metabolism by modulating energy balance in diabetic mice. Diabetes. 2000, 49: 436-444. 10.2337/diabetes.49.3.436
Nonomura T, Tsuchida A, Ono-Kishino M, Nakagawa T, Taiji M, Noguchi H: Brain-derived neurotrophic factor regulates energy expenditure through the central nervous system in obese diabetic mice. Int J Exp Diabetes Res. 2001, 2: 201-209. 10.1155/EDR.2001.201
Tsuchida A, Nonomura T, Ono-Kishino M, Nakagawa T, Taiji M, Noguchi H: Acute effects of brain-derived neurotrophic factor on energy expenditure in obese diabetic mice. Int J Obes Relat Metab Disord. 2001, 25: 1286-1293. 10.1038/sj.ijo.0801678
Lin JC, Tsao D, Barras P, Bastarrachea RA, Boyd B, Chou J, Rosete R, Long H, Forgie A, Abdiche Y, Dilley J, Stratton J, Garcia C, Sloane DL, Comuzzie AG, Rosenthal A: Appetite enhancement and weight gain by peripheral administration of TrkB agonists in non-human primates. PLoS One. 2008, 3 (Suppl 4): 1900-1907.
Das UN: Essential fatty acids: Biochemistry, physiology, and pathology. Biotechnology J. 2006, 1: 420-439. 10.1002/biot.200600012.
Das UN: Essential fatty acids- a review. Curr Pharm Biotechnol. 2006, 7: 467-482. 10.2174/138920106779116856
Das UN: Metabolic Syndrome Pathophysiology: The Role of Essential Fatty Acids Fatty Acids and their metabolites. 2010, Ames, IA, USA: Wiley-Blackwell Publishers.
Das UN: Molecular Basis of Health and Disease. 2011, New York: Springer.
Wu A, Ying Z, Gomez-Pinilla F: Dietary omega-3 fatty acids normalize BDNF levels, reduce oxidative damage, and counteract learning disability after traumatic brain injury in rats. J Neurotrauma. 2004, 21: 1457-1467. 10.1089/neu.2004.21.1457
Taouis M, Dagou C, Ster C, Durand G, Pinault M, Delarue J: n-3 polyunsaturated fatty acids prevent the defect of insulin receptor signaling in muscle. Am J Physiol Endocrinol Metab. 2002, 282: E664-671.
Stettler R, Ith M, Acheson KJ, Decombaz J, Boesch C, Tappy L, Binnert C: Interaction between dietary lipids and physical inactivity on insulin sensitivity and on intramyocellular lipids in healthy men. Diabetes Care. 2005, 28: 1404-1409. 10.2337/diacare.28.6.1404
Zhang W, Li P, Hu X, Zhang F, Chen J, Gao Y: Omega-3 polyunsaturated fatty acids in the brain: metabolism and neuroprotection. Front Biosci. 2011, 17: 2653-2670.
Wurtman RJ: Synapse formation and cognitive brain development: effect of docosahexaenoic acid and other dietary constituents. Metabolism. 2008, 57 (Suppl 2): S6-S10.
Darios F, Davletov B: Omega-3 and omega-6 fatty acids stimulate cell membrane expansion by acting on syntaxin 3. Nature. 2006, 440: 813-817. 10.1038/nature04598
Perez MT, Caminos E: Expression of brain-derived neurotrophic factor and of its functional receptor in neonatal and adult rat retina. Neurosci Lett. 1995, 183: 96-99. 10.1016/0304-3940(94)11123-Z
Zheng ZL, Ling QZ, Fernette FE: Critical Role of TrkB and Brain-Derived Neurotrophic Factor in the Differentiation and Survival of Retinal Pigment Epithelium. J Neurosci. 1997, 17: 8749-8755.
Robinson RC, Radziejewski C, Stuart DI, Jones EY: Structure of the brain-derived neurotrophic factor/neurotrophin 3 heterodimer. Biochemistry. 1995, 34: 4139-4146. 10.1021/bi00013a001
Van Der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJ: GROMACS: fast, flexible, and free. J Comput Chem. 2005, 26: 1701-18. 10.1002/jcc.20291
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983, 79: 926-935. 10.1063/1.445869.
Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH: PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009, 37: W623-633. 10.1093/nar/gkp456
Schüttelkopf AW, Van Aalten DM: PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallogr D: Biol Crystallogr. 2004, 60: 1355-63. 10.1107/S0907444904011679.
Goodsell DS, Morris GM, Olson AJ: Automated docking of flexible ligands: applications of AutoDock. J Mol Recognit. 1996, 9: 1-5. 10.1002/(SICI)1099-1352(199601)9:1<1::AID-JMR241>3.0.CO;2-6
Corbeil CR, Englebienne P, Moitessier N: Docking ligands into flexible and solvated macromolecules 1. Development and validation of FITTED 1.0. J Chem Inf Model. 2007, 47: 435-449. 10.1021/ci6002637
Rarey M, Kramer B, Lengauer T, Klebe G: A fast flexible docking method using an incremental construction algorithm. J Mol Biol. 1996, 261: 470-489. 10.1006/jmbi.1996.0477
Chandra N, Bhagava R, Sharma E, Sreekanthreddy P, Somasundaram K: Virtual screening, identification and experimental testing of novel inhibitors of PBEF1/Visfatin/NMPRTase for glioma therapy. J Clin Bioinforma. 2011, 1: 5- 10.1186/2043-9113-1-5
Taly A, Colas C, Malliavin T, Blondel A, Nilges M, Corringer PJ, Joseph D: Discrimination of agonists versus antagonists of nicotinic ligands based on docking onto AChBP structures. J Mol Graph Model. 2011, 30: 100-109.
Wallace AC, Laskowski RA, Thornton JM: LIGPLOT a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995, 8: 127-134. 10.1093/protein/8.2.127
Nitta A, Ito M, Fukumitsu H, Ohmiya H, Sometani A, Nomoto H, Furukawa Y, Furukawa S: 4-methylcatechol increases brain-derived neurotrophic factor content and mRNA expression in cultured brain cells and in rat brain in vivo. J Pharmacol Exp Ther. 1999, 291: 1276-1283.
Cazorla M, Prémont J, Mann A, Girard N, Kellendonk C, Rognan D: Identification of a low-molecular weight TrkB antagonist with anxiolytic and antidepressant activity in mice. J Clin Invest. 2011, 121: 1846-1857. 10.1172/JCI43992
Das UN: Obesity: Genes, Brain, Gut and Environment. Nutrition. 2010, 26: 459-473. 10.1016/j.nut.2009.09.020
Pae CU, Chiesa A, Porcelli S, Han C, Patkar AA, Lee SJ, Park MH, Serretti A, De Ronchi D: Influence of BDNF variants on diagnosis and response to treatment in patients with major depression, bipolar disorderand schizophrenia. Neuropsychobiology. 2011, 65 (Suppl 3): 1-11.
Yary T, Aazami S: The association between polyunsaturated fatty acids and depression among Iranian postgraduate students in Malaysia. Lipids Health Dis. 2011, 10: 151- 10.1186/1476-511X-10-151
This study was funded by the Department of Biotechnology, Government of India (DBT No. BT/PR11627/MED/30/157/2010). Dr. U N Das was in receipt of Ramalingaswami Fellowship of the Department of Biotechnology, New Delhi, India during the tenure of this study. This study was funded, in part, by a grant from the Department of Science and Technology to UND (No. IR/SO/LU/03/2008/1) under Intensification of Research in High Priority Areas (IRPHA).
The authors declare that they have no competing interests.
UV and SBR performed the bioinformatics and molecular docking studies. KK and JM performed the literature survey, helped in the bioinformatics studies. UND originated the idea and proposed the hypothesis that an interaction occurs between PUFAs and their metabolites and BDNF, initiated the bioinformatics studies and supervised the study. AN supervised and co-ordinated the study. All authors drafted the manuscript, read and approved the final version.