Istanbul Technical University Graduate School of Science
Total Page:16
File Type:pdf, Size:1020Kb
ISTANBUL TECHNICAL UNIVERSITY GRADUATE SCHOOL OF SCIENCE ENGINEERING AND TECHNOLOGY IN SILICO DESIGN OF HERG NON-BLOCKER COMPOUNDS WITH RETAINED PHARMACOLOGICAL ACTIVITY USING MULTI-SCALE MOLECULAR MODELING APPLICATIONS Ph.D. THESIS Gülru KAYIK Chemistry Department Chemistry Programme DECEMBER 2017 ISTANBUL TECHNICAL UNIVERSITY GRADUATE SCHOOL OF SCIENCE ENGINEERING AND TECHNOLOGY IN SILICO DESIGN OF HERG NON-BLOCKER COMPOUNDS WITH RETAINED PHARMACOLOGICAL ACTIVITY USING MULTI-SCALE MOLECULAR MODELING APPLICATIONS Ph.D. THESIS Gülru KAYIK (509112008) Chemistry Department Chemistry Programme Thesis Advisor: Prof. Dr. Nurcan TÜZÜN Thesis Co-Advisor: Assoc. Prof. Dr. Serdar DURDAĞI DECEMBER 2017 İSTANBUL TEKNİK ÜNİVERSİTESİ FEN BİLİMLERİ ENSTİTÜSÜ HERG BLOKER OLMAYAN FARMAKOLOJİK AKTİVİTESİ KORUNMUŞ BİLEŞİKLERİN ÇOK BOYUTLU MOLEKÜLER MODELLEME UYGULAMALARI İLE İN SİLİKO TASARIMI DOKTORA TEZİ Gülru KAYIK (509112008) Kimya Anabilim Dalı Kimya Programı Tez Danışmanı: Prof. Dr. Nurcan TÜZÜN Eş Danışman: Doç. Dr. Serdar DURDAĞI ARALIK 2017 Gülru KAYIK, a Ph.D. student of İTU Graduate School of Science Engineering and Technology student ID 509112008, successfully defended the thesis entitled “IN SILICO DESIGN OF HERG NON-BLOCKER COMPOUNDS WITH RETAINED PHARMACOLOGICAL ACTIVITY USING MULTI-SCALE MOLECULAR MODELING APPLICATIONS”, which she prepared after fulfilling the requirements specified in the associated legislations, before the jury whose signatures are below. Thesis Advisor : Prof. Dr. Nurcan TÜZÜN .............................. Istanbul Technical University Co-advisor : Assoc. Prof. Dr. Serdar DURDAĞI .............................. Bahçeşehir University Jury Members : Prof. Dr. Mine YURTSEVER ............................. Istanbul Technical University Prof. Dr. Kemal YELEKÇİ .............................. Kadir Has University Assis. Prof. Dr. Bülent BALTA .............................. Istanbul Technical University Assoc. Prof. Dr. Fethiye Aylin SUNGUR .............................. Istanbul Technical University Prof. Dr. Safiye ERDEM .............................. Marmara University Date of Submission : 26 October 2017 Date of Defense : 01 December 2017 v vi FOREWORD First of all, I would like to thank my Ph.D. thesis advisor Prof. Dr. Nurcan Tüzün and co-advisor Assoc. Prof. Dr. Serdar Durdağı for their kind concern, recommendations and supports during the course of my Ph.D. studies. I would like to present my acknowledgements to Istanbul Technical University Research Fund BAP (Project numbers: 38208 and 30492) and the National Center for High Performance Computing of Turkey (UHEM) under Grant 10982010 for supporting this thesis and providing the related computer resources. The numerical calculations reported in this thesis were partially performed at TUBITAK ULAKBIM, High Performance and Grid Computing Center (TRUBA resources). I would also like to thank The Scientific and Technological Research Council of Turkey (TUBITAK) for granting me the 2214-A Research Grant and providing financial support during the course of my Ph.D. thesis. December 2017 Gülru KAYIK (Chemical Engineer) vii viii TABLE OF CONTENTS Page FOREWORD ............................................................................................................ vii TABLE OF CONTENTS .......................................................................................... ix ABBREVIATIONS ................................................................................................. xiii SYMBOLS ................................................................................................................ xv LIST OF TABLES ................................................................................................. xvii LIST OF FIGURES ................................................................................................ xix SUMMARY ............................................................................................................ xxv ÖZET......... ............................................................................................................ xxvii 1. INTRODUCTION .................................................................................................. 1 2. IN SILICO DESIGN OF NOVEL HERG-NEUTRAL SILDENAFIL LIKE PDE5 INHIBITORS ............................................................................... 3 2.1 Introduction ........................................................................................................ 3 2.2 Methods .............................................................................................................. 5 2.2.1 Molecular docking simulations ................................................................... 6 2.2.2 Fragment-based de novo drug design & virtual screening .......................... 7 2.2.3 MD simulations and post-processing MD analyses .................................... 8 2.2.4 MD simulations of the target receptor: PDE5 in its apo state and bound with its inhibitors ................................................................................................. 8 2.2.5 MD simulations of hERG K+ ion channel: Apo state and bound with PDE5 inhibitors .................................................................................................... 9 2.2.6 Molecular Mechanics/Generalized Born surface area (MM/GBSA) calculations ......................................................................................................... 11 2.3 Results and Discussion ..................................................................................... 11 2.3.1 Analysis of the key interactions of sildenafil with the target receptor (PDE5) ............................................................................................................... 13 2.3.2 Comparison of used docking tools in terms of predicting the binding positions of “Sildenafil” in the central cavities of hERG1 channel ................... 15 2.3.2.1 GOLD ................................................................................................. 16 2.3.2.2 AutoDock ........................................................................................... 21 2.3.2.3 MOE ................................................................................................... 24 2.3.3 In silico Alanine mutagenesis study .......................................................... 25 2.3.4 General statements on the binding energy predictions derived from GOLD, AutoDock and MOE ............................................................................. 27 2.3.5 Binding interactions of Vardenafil and Tadalafil with the hERG K+ ion channel ............................................................................................................... 27 2.3.6 Virtual screening results ............................................................................ 29 2.3.7 MM/GBSA analyses ................................................................................. 30 2.4 Conclusions ...................................................................................................... 30 3. INVESTIGATION OF PDE5/PDE6 AND PDE5/PDE11 SELECTIVE POTENT TADALAFIL-LIKE PDE5 INHIBITORS USING COMBINATION OF MOLECULAR MODELING APPROACHES, MOLECULAR FINGERPRINT-BASED VIRTUAL SCREENING ix PROTOCOLS AND STRUCTURE-BASED PHARMACOPHORE DEVELOPMENT ................................................................................. 37 3.1 Introduction ...................................................................................................... 37 3.2 Methods ............................................................................................................ 41 3.2.1 Ligand and protein preparations ................................................................ 41 3.2.2 Virtual library screening............................................................................ 41 3.2.3 Flexible molecular docking simulations ................................................... 41 3.2.4 Molecular Dynamics simulations .............................................................. 42 3.2.5 Molecular Mechanics Generalized Born Solvation (MM/GBSA) Calculations ........................................................................................................ 42 3.3 Results and discussion ...................................................................................... 43 3.3.1 Validation of the docking methodology .................................................... 43 3.3.2 Constructing the homology models of the catalytic domains of PDE6 (amino acid residues: 482-816) and PDE11 (amino acid residues: 587-910) .... 43 3.3.3 Binding affinity and binding pattern analysis of the hit compounds and tadalafil with PDE5, PDE6 and PDE11 ............................................................. 47 3.3.4 MD simulations of apo and holo states of PDE5, PDE6 and PDE11 bound with the selected hit compounds (ZINC02120502 and ZINC16031243) and tadalafil ............................................................................................................... 59 3.3.5 MM-GBSA calculations ............................................................................ 62 3.3.6 hERG K+ ion channel activity of the compounds ..................................... 63 3.3.7 E-Pharmacophore studies .......................................................................... 65 3.4 Conclusions ...................................................................................................... 68 4. STRUCTURAL