TY - THES A1 - Merget, Benjamin T1 - Computational methods for assessing drug-target residence times in bacterial enoyl-ACP reductases and predicting small-molecule permeability for the \(Mycobacterium\) \(tuberculosis\) cell wall T1 - Computermethoden zur Bestimmung von Protein-Ligand Verweilzeiten in bakteriellen Enoyl-ACP Reduktasen und Vorhersage der Permeabilitätswahrscheinlichkeit kleiner Moleküle gegenüber der \(Mycobacterium\) \(tuberculosis\) Zellwand N2 - \textbf{Molecular Determinants of Drug-Target Residence Times of Bacterial Enoyl-ACP Reductases.} Whereas optimization processes of early drug discovery campaigns are often affinity-driven, the drug-target residence time $t_R$ should also be considered due to an often strong correlation with \textit{in vivo} efficacy of compounds. However, rational optimization of $t_R$ is not straightforward and generally hampered by the lack of structural information about the transition states of ligand association and dissociation. The enoyl-ACP reductase FabI of the fatty acid synthesis (FAS) type II is an important drug-target in antibiotic research. InhA is the FabI enzyme of \textit{Mycobacterium tuberculosis}, which is known to be inhibited by various compound classes. Slow-onset inhibition of InhA is assumed to be associated with the ordering of the most flexible protein region, the substrate binding loop (SBL). Diphenylethers are one class of InhA inhibitors that can promote such SBL ordering, resulting in long drug-target residence times. Although these inhibitors are energetically and kinetically well characterized, it is still unclear how the structural features of a ligand affect $t_R$. Using classical molecular dynamics (MD) simulations, recurring conformational families of InhA protein-ligand complexes were detected and structural determinants of drug-target residence time of diphenyl\-ethers with different kinetic profiles were described. This information was used to deduce guidelines for efficacy improvement of InhA inhibitors, including 5'-substitution on the diphenylether B-ring. The validity of this suggestion was then analyzed by means of MD simulations. Moreover, Steered MD (SMD) simulations were employed to analyze ligand dissociation of diphenylethers from the FabI enzyme of \textit{Staphylococcus aureus}. This approach resulted in a very accurate and quantitative linear regression model of the experimental $ln(t_R)$ of these inhibitors as a function of the calculated maximum free energy change of induced ligand extraction. This model can be used to predict the residence times of new potential inhibitors from crystal structures or valid docking poses. Since correct structural characterization of the intermediate enzyme-inhibitor state (EI) and the final state (EI*) of two-step slow-onset inhibition is crucial for rational residence time optimization, the current view of the EI and EI* states of InhA was revisited by means of crystal structure analysis, MD and SMD simulations. Overall, the analyses affirmed that the EI* state is a conformation resembling the 2X23 crystal structure (with slow-onset inhibitor \textbf{PT70}), whereas a twist of residues Ile202 and Val203 with a further opened helix $\alpha 6$ corresponds to the EI state. Furthermore, MD simulations emphasized the influence of close contacts to symmetry mates in the SBL region on SBL stability, underlined by the observation that an MD simulation of \textbf{PT155} chain A with chain B' of a symmetry mate in close proximity of the SBL region showed significantly more stable loops, than a simulation of the tetrameric assembly. Closing Part I, SMD simulations were employed which allow the delimitation of slow-onset InhA inhibitors from rapid reversible ligands. \textbf{Prediction of \textit{Mycobacterium tuberculosis} Cell Wall Permeability.} The cell wall of \textit{M. tuberculosis} hampers antimycobacterial drug design due to its unique composition, providing intrinsic antibiotic resistance against lipophilic and hydrophilic compounds. To assess the druggability space of this pathogen, a large-scale data mining endeavor was conducted, based on multivariate statistical analysis of differences in the physico-chemical composition of a normally distributed drug-like chemical space and a database of antimycobacterial--and thus very likely permeable--compounds. The approach resulted in the logistic regression model MycPermCheck, which is able to predict the permeability probability of small organic molecules based on their physico-chemical properties. Evaluation of MycPermCheck suggests a high predictive power. The model was implemented as a freely accessible online service and as a local stand-alone command-line version. Methodologies and findings from both parts of this thesis were combined to conduct a virtual screening for antimycobacterial substances. MycPermCheck was employed to screen the chemical permeability space of \textit{M. tuberculosis} from the entire ZINC12 drug-like database. After subsequent filtering steps regarding ADMET properties, InhA was chosen as an exemplary target. Docking to InhA led to a principal hit compound, which was further optimized. The quality of the interaction of selected derivatives with InhA was subsequently evaluated using MD and SMD simulations in terms of protein and ligand stability, as well as maximum free energy change of induced ligand egress. The results of the presented computational experiments suggest that compounds with an indole-3-acethydrazide scaffold might constitute a novel class of InhA inhibitors, worthwhile of further investigation. N2 - \textbf{Molekulare Determinanten von Wirkstoff-Angriffsziel Verweilzeiten bakterieller Enoyl-ACP Reduktasen.} In frühen Phasen der Wirkstoffentwicklung sind Optimierungsprozesse häufig affini\-täts\-geleitet. Darüber hinaus sollte zusätzlich die Wirkstoff-Angriffsziel Verweilzeit $t_R$ berücksichtigt werden, da diese oft eine starke Korrelation zur \textit{in vivo} Wirksamkeit der Substanzen aufweist. Rationale Optimierung von $t_R$ ist jedoch auf Grund eines Mangels an struktureller Information über den Übergangszustand der Ligandbindung und Dissoziierung nicht einfach umsetzbar. Die Enoyl-ACP Reduktase FabI der Fettsäurebio\-synthese (FAS) Typ II ist ein wichtiger Angriffspunkt in der Antibiotikaforschung. InhA ist das FabI Enzym des Organismus \textit{Mycobacterium tuberculosis} und kann durch Substanzen diverser Klassen gehemmt werden. Es wird vermutet, dass Hemmung von InhA durch langsam-bindende (``slow-onset'') Inhibitoren mit der Ordnung der flexibelsten Region des Enzyms assoziiert ist, dem Substratbindungsloop (SBL). Diphenylether sind eine InhA Inhibitorenklasse, die eine solche SBL Ordnung fördern und dadurch lange Verweilzeiten im Angriffsziel aufweisen. Obwohl diese Inhibitoren energetisch und kinetisch gut charakterisiert sind, ist noch immer unklar, wie die strukturellen Eigenschaften eines Liganden $t_R$ beeinflussen. Durch die Verwendung klassischer Molekulardynamik (MD) Simulationen wurden wiederkehrende Konformationsfamilien von InhA Protein-Ligand Komplexen entdeckt und strukturelle Determinanten der Wirkstoff-Angriffsziel Verweilzeit von Diphenylethern mit verschiedenen kinetischen Profilen beschrieben. Anhand dieser Ergebnisse wurden Richtlinien zur Wirksamkeitsoptimierung von InhA Inhibitoren abgeleitet, einschließlich einer 5'-Substitution am Diphenylether B-Ring. Die Validität dieses Vorschlags wurde mittels MD Simulationen nachfolgend analysiert. Darüber hinaus wurden ``Steered MD'' (SMD) Simulationen als MD Technik für umfangreicheres Sampling verwendet um die Liganddissoziation von Diphenylethern aus dem FabI Enzym von \textit{Staphylococcus aureus} zu untersuchen. Dieser Ansatz resultierte in einem sehr akkuraten, quantitativen linearen Regressionsmodell der experimentellen Verweilzeit $ln(t_R)$ dieser Inhibitoren als Funktion der berechneten maximalen freien Energieänderung induzierter Ligandextraktion. Dieses Modell kann genutzt werden um die Verweilzeiten neuer potentieller Inhibitoren aus Kristallstrukturen oder validen Dockingposen vorherzusagen. Die korrekte strukturelle Charakterisierung des intermediären und des finalen Zustandes (EI und EI*-Zustand) eines Enzym-Inhibitor Komplexes bei einem zweistufigen Inhibitionsmechanismus durch langsam-bindende Hemmstoffe ist essentiell für rationale Verweilzeitoptimierung. Daher wurde die gegenwärtige Ansicht des EI und EI*-Zustandes von InhA mittels Kristallstrukturanalyse, MD und SMD Simulationen erneut aufgegriffen. Insgesamt bestätigten die Analysen, dass der EI*-Zustand einer Konformation ähnlich der 2X23 Kristallstruktur (mit langsam-bindenden Inhibitor \textbf{PT70}) gleicht, während eine Drehung der Reste Ile202 und Val203 mit einer weiter geöffneten Helix $\alpha 6$ dem EI-Zustand entspricht. Des Weiteren zeigten MD Simulationen den Einfluss naher Kristallkontakte zu Symmetrie-Nachbarn in der SBL Region auf die SBL Stabilität. Dies wird durch die Beobachtung hervorgehoben, dass die Ketten A und B' eines InhA-\textbf{PT155}-Komplexes und des angrenzenden Symmetrie-Nachbars, welche in engem Kontakt in der SBL Region stehen, signifikant stabilere SBLs aufweisen, als die Ketten A und B in einer Simulation des Tetramers. Zum Abschluss von Teil I wurden SMD Simulationen angewandt, auf deren Basis es möglich war, langsam-bindende InhA Inhibitoren von schnell-reversiblen (``rapid reversible'') Liganden zu unterscheiden. \textbf{Vorhersage von \textit{Mycobacterium tuberculosis} Zellwand Permeabilität.} Die Zellwand von \textit{M.~tuberculosis} erschwert die antimycobakterielle Wirkstofffindung auf Grund ihrer einzigartigen Zusammensetzung und bietet eine intrinsische Antibiotikaresistenz gegenüber lipophilen und hydrophilen Substanzen. Um den chemischen Raum wirkstoffähnlicher Moleküle gegen diesen Erreger (``Druggability Space'') einzugrenzen, wurde eine groß angelegte Dataminingstudie durchgeführt, welche auf multivariater statistischer Analyse der Unterschiede der physikochemischen Zusammensetzung eines normalverteilten wirkstoffähnlichen chemischen Raumes und einer Datenbank von antimycobakteriellen -- und somit höchstwahrscheinlich permeablen -- Substanzen beruht. Dieser Ansatz resultierte in dem logistischen Regressionsmodell MycPermCheck, welches in der Lage ist die Permeabilitätswahrscheinlichkeit kleiner organischer Moleküle anhand ihrer physikochemischen Eigenschaften vorherzusagen. Die Evaluation von MycPermCheck deutet auf eine große Vorhersagekraft hin. Das Modell wurde als frei zugänglicher online Service und als lokale Kommandozeilenversion implementiert. Methodiken und Ergebnisse aus beiden Teilen dieser Dissertation wurden kombiniert um ein virtuelles Screening nach antimycobakteriellen Substanzen durchzuführen. Myc\-PermCheck wurde verwendet um den chemischen Permeabilitätsraum von \textit{M.~tuberculosis} anhand der gesamten ZINC12 Datenbank wirkstoffähnlicher Moleküle abzuschätzen. Nach weiteren Filterschritten mit Bezug auf ADMET Eigenschaften, wurde InhA als exemplarisches Angriffsziel ausgewählt. Docking nach InhA führte schließlich zu einer Treffersubstanz, welche in darauffolgenden Schritten weiter optimiert wurde. Die Interaktionsqualität ausgewählter Derivate mit InhA wurde daraufhin mittels MD und SMD Simulationen in Bezug auf Protein und Ligand Stabilität, sowie auch der maximalen freien Energieänderung induzierter Ligandextraktion, untersucht. Die Ergebnisse der vorgestellten computerbasierten Experimente legen nahe, dass Substanzen mit einem Indol-3-Acethydrazid Gerüst eine neuartige Klasse von InhA Inhibitoren darstellen könnten. Weiterführende Untersuchungen könnten sich somit als lohnenswert erweisen. KW - Computational chemistry KW - Arzneimitteldesign KW - Molekulardynamik KW - Permeabilität KW - Tuberkelbakterium KW - Computational drug-design KW - steered molecular dynamics KW - molecular dynamics KW - residence time KW - mycobacterium tuberculosis KW - staphylococcus aureus KW - permeability KW - InhA KW - FabI KW - Enoyl-acyl-carrier-protein-Reductase KW - Drug design KW - Computational chemistry Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-127386 ER - TY - THES A1 - Hirschbeck, Maria Wenefriede T1 - Structure-based drug design on the enoyl-ACP reductases of Yersinia pestis and Burkholderia pseudomallei T1 - Struktur-basiertes Wirkstoffdesign an Enoyl-ACP Reductase von Yersinia pestis und Burkholderia pseudomallei N2 - Spreading drug resistances among Gram-negative pathogens and the paucity of new agents on the antibacterial drug market against these tenacious bacteria create a pressing need for the development of new antibiotics. The bacterial fatty acid biosynthesis pathway FAS-II, especially the enoyl-ACP reductase catalyzing the last step of the elongation cycle, is an established drug target against tuberculosis but has not been extensively exploited for drug design against other bacterial pathogens. In this thesis the enoyl-ACP reductases of the Gram-negative biothreat organisms Burkholderia pseudomallei and Yersinia pestis were targeted in a structure-based drug design approach. The structure of the most recently identified enoyl-ACP isoenzyme FabV was characterized by X-ray crystallography and could be determined in three different states. FabV from B. pseudomallei was obtained in the apo-form of the enzyme, whereas FabV from Y. pestis was characterized in a binary complex with the cofactor NADH as well as in a ternary complex with NADH and the triclosan-based 2-pyridone inhibitors PT172 and PT173. Analysis of the FabV structure revealed the typical fold of the short chain dehydrogenase/reductase superfamily with the NADH-binding Rossmann fold and a substrate-binding pocket with a conserved active site geometry compared to the related isoenzyme FabI. Additional structural elements of FabV are located around the active site. The monomeric form of the enzyme is thereby stabilized and the substrate-binding loop is kept in a closed, helical conformation. The ternary complexes of FabV exhibited a similar inhibitor-binding mode as observed for triclosan inhibition in FabI and point to a potential substrate-binding mechanism. B. pseudomallei possesses FabI as an additional enoyl-ACP reductase isoenzyme, which was structurally characterized in the apo form and in ternary complexes with NAD+ and the diphenyl ether inhibitors triclosan, PT02, PT12 or PT404 as well as the 4-pyridone inhibitor PT155. The structural data of the ternary enoyl-ACP reductases complexes of B. pseudomallei and Y. pestis hold the promise for the possibility to develop antibacterials targeting FabV or even both isoenzymes, FabI and FabV, based on the triclosan scaffold. N2 - Die Ausbreitung von Antibiotikaresistenzen in Gram-negativen Pathogenen sowie der Mangel neuer Medikamente auf dem Arzneimittelmarkt gegen diese hartnäckigen Bakterien weist einen dringenden Bedarf an neuen Antibiotika auf. Die bakterielle Fettsäurebiosynthese (FAS-II), speziell die Enoyl-ACP-Reduktase, welche den finalen Schritt des Elongationszyklus katalysiert, ist ein etablierter Angriffspunkt in der Tuberkulosetherapie. Sie wurde jedoch bisher noch nicht für die gezielte Wirkstoffentwicklung gegen andere bakterielle Krankheitserreger genutzt. In dieser Dissertation waren die Enoyl-ACP Reduktasen aus Burkholderia pseudomallei und Yersinia pestis Gegenstand des strukturbasierten Wirkstoffdesigns. Die Struktur des zuletzt gefundenen Isoenzyms der Enoyl-ACP-Reduktase, FabV, wurde röntgenstrukturanalytisch charakterisiert und konnte in drei verschiedenen Zuständen bestimmt werden. Die Struktur des FabV Proteins aus B. pseudomallei wurde in der Apo-Form gelöst, während FabV aus Y. pestis in binären und ternären Komplexen mit NADH bzw. NADH und einem Triclosan-basierten 2-Pyridon-Inhibitor, PT172 bzw. PT173 charakterisiert wurde. FabV weist die typische Struktur eines Mitglieds der Short-Chain-Dehydrogenase/Reduktase Superfamilie mit einer NADH-bindenden Rossmann-Faltung und einer Substratbindetasche auf mit einer, im Vergleich zu dem verwandten Isoenzym FabI, konservierte Geometrie des aktiven Zentrums. Zusätzliche strukturelle Elemente sind um das aktive Zentrum gefaltet und stabilisieren damit das Enzym in seiner monomeren Form. Darüber hinaus halten sie den Substratbindeloop in einer geschlossenen helikalen Gestalt. Die ternären FabV Komplexe zeigen Übereinstimmungen mit dem bekannten Bindungsmechanismus des Inhibitors Triclosan in FabI und deuten auf einen möglichen Substratbindungsmechanismus hin. B. pseudomallei besitzt FabI als zusätzliches Isoenzym der Enoyl-ACP-Reduktasen. Dieses Isoenzym wurde in der Apo-Form und in ternären Komplexen mit NAD+ und den Diphenylether-Inhibitoren Triclosan, PT02, PT12 und PT404 sowie dem 4-Pyridon-Inhibitor PT155 strukturell charakterisiert. Die strukturellen Daten der ternären Enoyl-ACP-Reduktase Komplexe von B. pseudomallei und Y. pestis stellen die Möglichkeit in Aussicht Antibiotika zu entwickeln, welche FabV oder auch beide Isoenzyme, FabI und FabV, inhibieren. KW - Yersinia KW - Burkholderia KW - Arzneimitteldesign KW - Fettsäurestoffwechsel KW - Struktur-basiertes Wirkstoff Design KW - Fettsäurebiosynthese KW - Triclosan KW - FabI KW - FabV KW - Yersinia pestis KW - Burkholderia pseudomallei KW - FAS-II Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-70869 ER -