Discovery of Human HMG-Coa Reductase Inhibitors Using Structure-Based Pharmacophore Modeling Combined with Molecular Dynamics Simulation Methodologies

Authors: Minky Son, Chanin Park, Ayoung Baek, Shalini John, Keun Woo Lee

Abstract:

3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) catalyzes the conversion of HMG-CoA to mevalonate using NADPH and the enzyme is involved in rate-controlling step of mevalonate. Inhibition of HMGR is considered as effective way to lower cholesterol levels so it is drug target to treat hypercholesterolemia, major risk factor of cardiovascular disease. To discover novel HMGR inhibitor, we performed structure-based pharmacophore modeling combined with molecular dynamics (MD) simulation. Four HMGR inhibitors were used for MD simulation and representative structure of each simulation were selected by clustering analysis. Four structure-based pharmacophore models were generated using the representative structure. The generated models were validated used in virtual screening to find novel scaffolds for inhibiting HMGR. The screened compounds were filtered by applying drug-like properties and used in molecular docking. Finally, four hit compounds were obtained and these complexes were refined using energy minimization. These compounds might be potential leads to design novel HMGR inhibitor.

Keywords: Anti-hypercholesterolemia drug, HMGR inhibitor, Molecular dynamics simulation, Structure-based pharmacophore modeling.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1333334

References:

[1] E.S. Istvan, J. Deisenhofer, "The structure of the catalytic portion of human HMG-CoA reductase", Biochimica et Biophysica Acta (BBA)-Molecular and Cell Biology of Lipids. vol.1529, pp. 9-18, 2000.
[2] L. Tabemero, D.A. Bochar, V.W. Rodwell, C.V. Stauffacher, "Substrate-induced closure of the flap domain in the ternary complex structures provides insights into the mechanism of catalysis by 3-hydroxy-3-methylglutaryl—CoA reductase", Proceedings of the National Academy of Sciences, vol. 96, pp. 7167, 1999.
[3] E.S. Istvan, J. Deisenhofer, "Structural mechanism for statin inhibition of HMG-CoA reductase", Science, vol. 292, pp. 1160-1164, 2001.
[4] E. Boumendil, P. Tubert-Bitter, "Depression-induced absenteeism in relation to antihyperlipidemic treatment: a study using GAZEL cohort data", Epidemiology, pp. 322-325, 1995.
[5] T. Phan, J. McLeod, J. Pollard, 0. Peiris, A. Rohan, J. Halpern, "Peripheral neuropathy associated with simvastatin", Journal of Neurology, Neurosurgery & Psychiatry, vol. 58, pp.625-628, 1995.
[6] E. Proksch, "Antilipemic drug-induced skin manifestations", Hautarzt, vol. 46, pp. 76, 1995.
[7] S. John, S. Thangapandian, K.W. Lee, "Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies", J. Biomol. Struct. Dyn., vol. 29, pp. 921-936, 2012.
[8] J.A. Pfefferkom, C. Choi, Y. Song, B.K. Trivedi, S.D. Larsen, V. Askew, L. Dillon, J.C. Hanselman, Z. Lin, G. Lu, "Design and synthesis of novel, conformationally restricted HMG-CoA reductase inhibitors", Bioorg. Med. Chem. Lett., vol. 17, pp. 4531-4537, 2000.
[9] J.A. Pfefferkorn, C. Choi, S.D. Larsen, B. Auerbach, R. Hutchings, W. Park, V. Askew, L. Dillon, J.C. Hanselman, Z. Lin, "Substituted Pyrazoles as Hepatoselective HMG-CoA Reductase Inhibitors: Discovery of (3 R, 5R) — 7 -
[ 2 - ( 4 - F luoro — phenyl ) — 4 — isopropyl —5 - ( 4 — methyl-benzylcarbamoyl)-2 H-pyrazol-3-yl]-3, 5-dihydroxyheptanoic Acid (PF-3052334) as a Candidate for the Treatment of Hypercholesterolemia", J. Med. Chem., vol. 51, pp. 31-45, 2007.
[10] H.J.C. Berendsen, D. van der Spoel, R. van Drunen, "GROMACS: A message-passing parallel molecular dynamics implementation", Comput. Phys. Commun., vol. 91, pp. 43-56, 1995.
[11] D. Van Der Spoel, E. Lindahl, B. Hess, G. Groenhof, A.E. Mark, H.J.C. Berendsen, "GROMACS: fast, flexible, and free", Journal of computational chemistry., vol. 26, pp. 1701-1718, 2005.
[12] V. Zoete, M.A. Cuendet, A. Grosdidier, 0. Michielin, "SwissParam: A fast force field generation tool for small organic molecules", Journal of computational chemistry., 2011.
[13] H. Berendsen, J. Postma, W. Van Gunsteren, J. Hermans, "Interaction models for water in relation to protein hydration" Intermolecular forces., vol. 331, pp. 331, 1981.
[14] W.L. Jorgensen, J. Chandrasekhar, J.D. Madura, R.W. Impey, M.L. Klein, "Comparison of simple potential functions for simulating liquid water", J. Chem. Phys., vol. 79, pp. 926, 1983.
[15] J.P. Ryckaert, G. Ciccotti, H.J.C. Berendsen, "Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of< i> n-alkanes", Journal of Computational Physics., vol. 23, pp. 327-341, 1977.
[16] B. Hess, H. Bekker, H.J.C. Berendsen, J.G.E.M. Fraaije, "LINCS: a linear constraint solver for molecular simulations", Journal of computational chemistry, vol. 18, pp. 1463-1472, 1997.
[17] Miyamoto, P.A. Kollman, "SETTLE: an analytical version of the SHAKE and RATTLE algorithm for rigid water models", Journal of computational chemistry., vol. 13, pp. 952-962, 1992.
[18] Darden, D. York, L. Pedersen, "Particle mesh Ewald: An N. log (N) method for Ewald sums in large systems", J. Chem. Phys., vol. 98, pp. 10089, 1993.
[19] Essmann, L. Perera, M.L. Berkowitz, T. Darden, H. Lee, L.G. Pedersen, "A smooth particle mesh Ewald method", J. Chem. Phys., vol. 103, pp. 8577, 1995.
[20] M. Parrinello, A. Rahman, "Polymorphic transitions in single crystals: A new molecular dynamics method", J. Appl. Phys., vol. 52, pp. 7182-7190, 1981.
[21] H.J.C. Berendsen, J.P.M. Postma, W.F. Van Gunsteren, A. DiNola, J. Haak, "Molecular dynamics with coupling to an external bath", J. Chem. Phys., vol. 81, pp. 3684, 1984.
[22] C.A. Lipinski, F. Lombardo, B.W. Dominy, P.J. Feeney, "Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settingsl", Adv. Drug Deliv. Rev., vol. 46, pp. 3-26, 2001.
[23] W.P. Walters, M.A. Murcko, "Prediction of 'drug-likeness"', Adv. Drug Deliv. Rev., vol. 54 , pp. 255-271, 2002.
[24] G. Jones, P. Willett, R.C. Glen, A.R. Leach, R. Taylor, " Development and validation of a genetic algorithm for flexible dockingl", J. MoL Biol. 267, vol. 267, pp. 727-748, 1997.
[25] M.L. Verdonk, J.C. Cole, M.J. Hartshorn, C.W. Murray, R.D. Taylor, "Improved protein—ligand docking using GOLD", Proteins: Structure, Function, and Bioinformatics., vol. 52, pp. 609-623, 2003.
[26] G.M. Morris, D.S. Goodsell, R.S. Halliday, R. Huey, W.E. Hart, R.K. Belew, A.J. Olson, "Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function", Journal of computational chemistry., vol. 19, pp. 1639-1662, 1998.
[27] G.M. Morris, R. Huey, W. Lindstrom, M.F. Sanner, R.K. Belew, D.S. Goodsell, A.J. Olson, "AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility", Journal of computational chemistry., vol. 30, pp. 2785-2791, 2009..
[28] A.B. Wagner, "SciFinder Scholar 2006: an empirical analysis of research topic query processing", Journal of chemical information and modeling., vol. 46, pp. 767-774, 2006.
[29] Y. Wang, E. Bolton, S. Dracheva, K. Karapetyan, B.A. Shoemaker, T.O. Suzek, J. Wang, J. Xiao, J. Zhang, S.H. Bryant, "An overview of the PubChem BioAssay resource", Nucleic Acids Res., vol. 38, pp. D255-D266, 2010.