UNIVERSITY of CALIFORNIA, IRVINE Development And
Total Page:16
File Type:pdf, Size:1020Kb
UNIVERSITY OF CALIFORNIA, IRVINE Development and application of molecular modeling methods for characterization of drug targets in Mycobacterium tuberculosis DISSERTATION submitted in partial satisfaction of the requirements for the degree of DOCTOR OF PHILOSOPHY in Pharmaceutical Sciences by Kalistyn H. Burley Dissertation Committee: Professor David L. Mobley, Chair Professor Celia W. Goulding Professor Thomas L. Poulos 2020 Chapter 2 © 2019 American Chemical Society Chapter 3 © 2019 American Chemical Society All other materials © 2020 Kalistyn H. Burley DEDICATION To my children, Aubrey and Baron – being your mother is, and will always be, my greatest accomplishment. To Scott, an amazing partner in life and in parenthood, who has never doubted me, and has also never let me settle. To my mom, dad, and brothers Kristopher and Jacob, who gave me the confidence to “jump off cliffs and build my wings on the way down.” - Ray Bradbury, 1986 And to my mentors and role models, particularly the women before me who have paved the road for my successes, I am forever indebted and committed to paying it forward. ii TABLE OF CONTENTS LIST OF FIGURES iv LIST OF TABLES vi ACKNOWLEDGMENTS vii VITA ix ABSTRACT OF THE DISSERTATION xii CHAPTER 1: An introduction to tuberculosis and 1 computational methods for structure-based drug design References 16 CHAPTER 2: Enhancing side chain sampling using 19 non-equilibrium candidate Monte Carlo References 59 CHAPTER 3: Structure of a Mycobacterium tuberculosis heme- 63 degrading protein, MhuD, variant in complex with its product References 92 CHAPTER 4: Insights from the structure and conformational 97 dynamics of Mycobacterium tuberculosis malic enzyme, MEZ References 134 CHAPTER 5: Summary and Conclusions - Complementing structural 139 information using computational methods for drug-design References 146 APPENDIX A: Chapter 2 Supporting Information 148 APPENDIX B: Chapter 3 Supporting Information 154 APPENDIX C: Chapter 4 Supporting Information 157 iii LIST OF FIGURES Page Figure 1.1 Percentage of new TB cases that are MDR or rifampicin 2 resistant (RR-TB) Figure 1.2 Mtb cell wall architecture 5 Figure 1.3 Mtb iron uptake 7 Figure 1.4 Computational methods for drug development 12 Figure 1.5 Approximate timescales of molecular motions 14 Figure 2.1 Cycling of atom interactions during execution of NCMC 26 side chain move Figure 2.2 Move biasing based on known side chain rotamer states 30 Figure 2.3 Workflow illustration of BLUES side chain proposals 32 Figure 2.4 Rotameric states of valine side chain 35 Figure 2.5 Umbrella sampling of valine-alanine peptide in explicit 38 solvent Figure 2.6 Valine-Alanine rotamer transition data where the 40 torsional force constant, k = 3.8 kcal/mol Figure 2.7 Valine-Alanine rotamer transition data where 42 k = 10 kcal/mol. Figure 2.8 Val111 χ1 rotamer data for BLUES simulations of 47 p-xylene bound T4 lysozyme L99A in explicit solvent Figure 2.9 Val111 χ1 transition rates in T4 lysozyme L99A bound to 47 p-xylene for BLUES and MD Figure 2.10 Backbone RMSDs for T4 lysozyme L99A simulations 50 Figure 2.11 Overlay of backbone RMSDs of T4 lysozyme L99A for 51 BLUES and MD (first 100ns) 6 Figure 2.12 Val111 χ1 state transitions per 10 FEs from simulations 53 from 3 distinct starting conformations of p-xylene bound T4 lysozyme L99A Figure 3.1 Structures of tetrapyrroles and WT-MhuD-monoheme 65 Figure 3.2 Heme and αBV affinity of MhuD and the R26S variant 73 Figure 3.3 Structure of the MhuD-R26S-αBV complex 76 Figure 3.4 Structural comparison of MhuD-heme-CN and 81 MhuD-R26S-αBV Figure 3.5 Helical stability in MD simulations of MhuD 82 Figure 3.6 Orientation of His75 in MD simulations of MhuD 83 Figure 3.7 Position of Arg79 side chain during MD simulations of 85 MhuD Figure 4.1 Malic enzyme structural classes and oligomeric states 100 Figure 4.2 Structural domains of MEZ 103 Figure 4.3 Biophysical and biochemical analyses of MEZ 109 Figure 4.4 Different modes of NAD(P)+ binding after docking with Mn2+ 112 and malate iv Figure 4.5 Divalent cation stabilizes malate in binding site 114 Figure 4.6 Average RMSF per MEZ residue 115 Figure 4.7 Conformational dynamics of MEZ active site 116 Figure 4.8 Putative NAD(P)+ binding modes in MEZ after MD 121 simulations v LIST OF TABLES Page Table 2.1 Data summary for simulations of T4 Lysozyme L99A with 53 p-xylene from alternate starting conformations. Table 3.1 Heme and αBV binding affinities (Kd) for MhuD 72 Table 3.2 Statistics for X-ray diffraction data collection and atomic 75 refinement for the MhuD-R26S-αBV complex. Table 4.1. Data collection and refinement statistics for MEZ crystal 106 structure. vi ACKNOWLEDGMENTS In my time as a graduate student, I have been fortunate to have had the guidance and mentorship of two outstanding research advisors: Dr. Celia Goulding and Dr. David Mobley. They are each highly revered in their respective fields and are deeply committed to their responsibility to nurture and train future scientists. Thank you, Celia, for applying steady pressure to move me forward in my research, persistently advocating on my behalf, and for helping me navigate the highs, lows, and idiosyncrasies of life as a graduate student. Thank you to David, for always making yourself available, encouraging and providing opportunities to engage in the wider computational chemistry community, looking out for my general welfare and for trusting in my ability to work independently. Thank you both for unequivocally supporting and respecting my decision to start and grow my family while in graduate school. To that end, I would also like to acknowledge and thank my husband Scott - I could not have done this without you. Thank you for being an amazing father, partner, and coach. Finally, a special thank you to Mary and Henry Bugielski - for providing outstanding care and a warm and nurturing environment for our children, each and every weekday. Chapter 2 is a minimally modified reprint of the published work: Burley KH, Gill SC, Lim NM, Mobley DL. Enhancing Side chain sampling using non-equilibrium Monte Carlo. J Chem Theory Comput. 2019 Jan 24;15(3):1848-1862. I would like to acknowledge financial support from the National Science Foundation (CHE 1352608) and the National Institutes of Health (1R01GM108889-01 and T32GM108561), and the Vertex Fellowship as well as computing support from the UCI GreenPlanet cluster, supported in part by NSF Grant CHE- 0840513 and the infrastructure support from the Triton Shared Computing Cluster (TSCC) at the San Diego Supercomputing Center (SDSC). I particularly thank my co-authors, Sam Gill for developing and pioneering the BLUES method and Nathan Lim for helping transform BLUES into a developer-friendly and user-friendly software package. I would also like to recognize Sukanya Sasmal (UCI) for help with docking, Gaetano Calabró (UCI) for basic MD training when joining the lab, and Victoria Lim (UCI) for help with umbrella sampling. Chapter 3 is a minimally modified reprint of the published work: Chao A, Burley KH, Sieminski PJ, de Miranda R, Chen X, Mobley DL, Goulding CW. Structure of Mycobacterium tuberculosis heme-degrading protein, MhuD, variant in complex with its product. Biochemistry. 2019 Oct 22;58(46):4610-20. This work and its authors were financially supported by the National Institutes of Health (NIH) (P01-AI095208, T32GM108561), the National Science Foundation (DGE-1321846). The authors are grateful to Tom Poulos for discussions and critical reading of the manuscript, as well as Dr. Fengyun Ni, Baylor College of Medicine for advice on the more difficult aspects of structure refinement. The authors also acknowledge the Advanced Light Source at Berkeley National Laboratories (ALS) and the Stanford Synchrotron Radiation Lightsource (SSRL) for their invaluable help in data collection, and The Triton Shared Computing Cluster (TSCC) at the San Diego Supercomputing Center (SSDC) for their computing support. Regarding my co- authors, I especially thank Alex Chao for solving the structure, Xiaorui Chen for her vii contributions in data processing and refinement, as well as Paul Sieminski and Rodger de Miranda for their biochemical characterization of the MhuD variant with its product. Chapter 4 is part of a manuscript in preparation. I would like to acknowledge Bonnie Cuthbert for her tenacity and expertise in refining the crystal structure of MEZ, her analysis of the structure, as well as her moral support and friendship, Ervin Irimpan for his assistance in preparing protein and setting up crystal trays, Ilona Foik for the well-kept records of her prior crystallization attempts and purification methods, and Robert Quechol for his work in performing critical biochemical assays. I also appreciate financial support from the National Institutes of Health (NIH) (P01-AI095208, T32GM108561) and the National Science Foundation (DGE-1321846), as well as computing support from The Triton Shared Computing Cluster (TSCC) at the San Diego Supercomputing Center (SSDC) and UCI GreenPlanet cluster, supported in part by NSF Grant CHE-0840513. viii VITA Kalistyn Hope Burley EDUCATION Ph.D. in Pharmaceutical Sciences June 2020 University of California, Irvine Irvine, CA Co-Advisors: Dr. David Mobley, Dr. Celia Goulding A.B. in Biomedical Engineering June 2009 Dartmouth College Hanover, NH PROFESSIONAL EXPERIENCE Graduate Student Researcher September 2015 – April 2020 Mobley Lab, UC Irvine, Department of Pharmaceutical Sciences Irvine, CA Graduate Student Researcher January 2016 – April 2020 Goulding Lab,