Refine
Has Fulltext
- yes (17)
Is part of the Bibliography
- yes (17)
Year of publication
Document Type
- Doctoral Thesis (17) (remove)
Keywords
- Arzneimitteldesign (17) (remove)
Die Interaktion des onkogenen Transkriptionsfaktors MYCN mit der Ser/Thr Kinase Aurora-A verhindert
dessen Abbau über das Ubiquitin Proteasomsystem indem die Rekrutierung des SCF FbxW7 Komplexes
verhindert wird. Die Kinase nimmt mit der Bindung an MYCN eine aktive Konformation ein und erhält
somit die Fähigkeit zur Kinaseaktivität ohne die sonst notwendige Phosphorylierung von Thr288 oder
die Anwesenheit eines Aktivators wie TPX2. Da hohe MYCN Konzentrationen Tumore wie
Neuroblastome antreiben, ist die Störung der Komplexbildung mit Aurora-A eine valide Strategie zur
Entwicklung von Chemotherapeutika. Einige Inhibitoren von Aurora-A wie Alisertib (MLN8237) sind in
der Lage, eine Konformationsänderung in der Kinase zu verursachen, die mit der Bindung von MYCN
inkompatibel ist und auf diese Weise den Abbau des Transkriptionsfaktors induziert. Da Aurora-A
wichtige Funktionen in der Mitose übernimmt, könnte eine direkte Adressierung des Komplexes anstelle
einer systemischen Inhibition der Kinase vielversprechender sein.
Ziel des Projektes war die Identifizierung von Molekülen, die selektiv an das Interface des
Aurora-A – MYCN Komplexes binden und weiter optimiert werden können, um einen gezielten Abbau
des Transkriptionsfaktors über einen PROTAC Ansatz zu ermöglichen. Virtuelle Screenings und
molekulardynamische Simulationen wurden durchgeführt, um kommerziell erhältliche Verbindungen zu
identifizieren, welche mit einer Bindetasche des Komplexes interagieren, die nur zustande kommt, wenn
beide Proteine miteinander interagieren. Aus einem ersten Set von zehn potentiellen Liganden wurde
für vier eine selektive Interaktion mit dem Protein – Protein Komplex gegenüber Aurora-A oder MYCN
alleine in STD-NMR Experimenten bestätigt. Zwei der Hits besaßen ein identisches Grundgerüst und
wurden als Ausganspunkt für die Optimierung zu potenteren Liganden genutzt. Das Gerüst wurde
fragmentweise vergrößert und in Richtung besserer in-silico Ergebnisse und Funktionalisierung zur
Anbringung von E3-Ligase-Liganden optimiert. Neun dieser Liganden der zweiten Generation wurden
synthetisiert.
Um quantitative Bindungsdaten zu erhalten, wurde ein kovalent verknüpftes Aurora-A – MYCN
Konstrukt entworfen. Die strukturelle und funktionale Integrität wurde in STD-NMR und BLI
Experimenten mit bekannten Aurora-A Inhibitoren bestätigt, sowie in NMR-basierten ATPase Assays.
Zusätzlich konnte die Kristallstruktur des Konstrukts gelöst und damit die Validität des Designs bestätigt
werden. Quantitative Messungen der synthetisierten Moleküle identifizierten HD19S als Hit mit einer
zehnfach höheren Affinität für das Aurora-A – MYCN Konstrukt im Vergleich zu der Kinase allein.
Zusätzlich wurden in-silico Untersuchungen zu PROTACs der Aurora-A Kinase durchgeführt.
Interaktionen zwischen Aurora-A, der E3-Ligase Cereblon und den Liganden wurden modelliert und für
die Erklärung unterschiedlicher Aktivitäten der eingesetzten PROTACs verwendet. Zudem zeigte das
aktivste PROTAC eine hohe Selektivität für Aurora-A gegenüber Aurora-B, obwohl die verwendete
Erkennungseinheit (Alisertib) an beide Aurora-Proteine bindet. Dieser Umstand konnte durch
energetische Analysen von molekulardynamischen Simulationen der ternären Komplexe erklärt werden.
Optimierungsmöglichkeiten für eine effizientere Degradation von Aurora-A durch die PROTACs wurden
basierend auf modifizierten Erkennungseinheiten und verbesserten Linkern untersucht.
Mycobacterium tuberculosis is the causative agent of tuberculosis and responsible for more than eight million new infections and about two million deaths each year. Novel chemotherapeutics are urgently needed to treat the emerging threat of multi drug resistant and extensively drug resistant strains. Cell wall biosynthesis is a widely used target for chemotherapeutic intervention in bacterial infections. In mycobacteria, the cell wall is comprised of mycolic acids, very long chain fatty acids that provide protection and allow the bacteria to persist in the human macrophage. The type II fatty acid biosynthesis pathway in Mycobacterium tuberculosis synthesizes fatty acids with a length of up to 56 carbon atoms that are the precursors of the critical mycobacterial cell wall components mycolic acids. KasA, the mycobacterial ß-ketoacyl synthase and InhA, the mycobacterial enoyl reductase, are essential enzymes in the fatty acid biosynthesis pathway and validated drug targets. In this work, KasA was expressed in Mycobacterium smegmatis, purified and co-crystallized in complex with the natural thiolactone antibiotic thiolactomycin (TLM). High-resolution crystal structures of KasA and the C171Q KasA variant, which mimics the acyl enzyme intermediate of the enzyme, were solved in absence and presence of bound TLM. The crystal structures reveal how the inhibitor is coordinated by the enzyme and thus specifically pinpoint towards possible modifications to increase the affinity of the compound and develop potent new drugs against tuberculosis. Comparisons between the TLM bound crystal structures explain the preferential binding of TLM to the acylated form of KasA. Furthermore, long polyethylene glycol molecules are bound to KasA that mimic a fatty acid substrate of approximately 40 carbon atoms length. These structures thus provide the first insights into the molecular mechanism of substrate recognition and reveal how a wax-like substance can be accommodated in a cytosolic environment. InhA was purified and co-crystallized in complex with the slow, tight binding inhibitor 2-(o-tolyloxy)-5-hexylphenol (PT70). Two crystal structures of the ternary InhA-NAD+-PT70 were solved and reveal how the inhibitor is bound to the substrate binding pocket. Both structures display an ordered substrate binding loop and corroborate the hypothesis that slow onset inhibition is coupled to loop ordering. Upon loop ordering, the active site entrance is more restricted and the inhibitor is kept inside more tightly. These studies provide additional information on the mechanistic imperatives for slow onset inhibition of enoyl ACP reductases.
Malaria and HIV are among the most important global health problems of our time and together are responsible for approximately 3 million deaths annually. These two diseases overlap in many regions of the world including sub-Saharan Africa, Southeast Asia and South America, leading to a higher risk of co-infection. In this study, we generated and characterized hybrid molecules to target P. falciparum and HIV simultaneously for a potential HIV/malaria combination therapy. Hybrid molecules were synthesized by covalent fusion between azidothymidine (AZT) and dihydroartemisinin (DHA), tetraoxane or chloroquine (CQ); and a small library was generated and tested for antiviral and antimalarial activity. Our data suggest that dihyate is the most potent molecule in vitro, with antiplasmodial activity comparable to that of DHA (IC50 = 26 nM, SI > 3000), a moderate activity against HIV (IC50 = 2.9 µM; SI > 35) and safe to HeLa cells at concentrations used in the assay (CC50 > 100 µM). Pharmacokinetic studies further revealed that dihyate is metabolically unstable and is cleaved following an O-dealkylation once in contact with cytochrome P450 enzymes. The later further explains the uneffectiveness of dihyate against the CQ-sensitive P. berghei N strain in mice when administered by oral route at 20 mg/kg. Here, we report on a first approach to develop antimalarial/anti-HIV hybrid molecules and future optimization efforts will aim at producing second generation hybrid molecules to improve activity against HIV as well as compound bioavailability. With the emergence of resistant parasites against all the counterpart drugs of artemisinin derivatives used in artemisinin based combination therapies (ACTs), the introduction of antibiotics in the treatment of malaria has renewed interest on the identification of antibiotics with potent antimalarial properties. In this study we also investigated the antiplasmodial potential of thiostrepton and derivatives, synthesized using combinations of tail truncation, oxidation, and addition of lipophilic thiols to the terminal dehydroamino acid. We showed that derivatives SS231 and SS234 exhibit a better antiplasmodial activity (IC50 = 1 µM SI > 59 and SI > 77 respectively) than thiostrepton (IC50 = 8.95 µM, SI = 1.7). The antiplasmodial activity of these derivatives was observed at concentrations which are not hemolytic and non-toxic to human cell lines. Thiostrepton and derivatives appeared to exhibit transmission blocking properties when administered at their IC50 or IC90 concentrations and our data also showed that they attenuate proteasome activity of Plasmodium, which resulted in an accumulation of ubiquitinated proteins after incubation with their IC80 concentrations. Our results indicate that the parasite’s proteasome could be an attractive target for therapeutic intervention. In this regard, thiostrepton derivatives are promising candidates by dually acting on two independent targets, the proteasome and the apicoplast, with the capacity to eliminate both intraerythrocytic asexual and transmission stages of the parasite. To further support our findings, we evaluated the activity of a new class of antimalarial and proteasome inhibitors namely peptidyl sulfonyl fluorides on gametocyte maturation and analogues AJ34 and AJ38 were able to completely suppress gametocytogenesis at IC50 concentrations (0.23 µM and 0.17 µM respectively) suggesting a strong transmission blocking potential. The proteasome, a major proteolytic complex, responsible for the degradation and re-cycling of non-functional proteins has been studied only indirectly in P. falciparum. In addition, an apparent proteasome-like protein with similarity to bacterial ClpQ/hslV threonine-peptidases was predicted in the parasite. Antibodies were generated against the proteasome subunits alpha type 5 (α5-SU), beta type 5 (β5-SU) and pfhslV in mice and we showed that the proteasome is expressed in both sexual and asexual blood stages of P. falciparum, where they localize in the nucleus and in the cytoplasm. However, expression of PfhslV was only observed in trophozoites and shizonts. The trafficking of the studied proteasome subunits was further investigated by generating parasites expressing GFP tagged proteins. The expression of α5-SU-GFP in transgenic parasite appeared to localize abundantly in the cytoplasm of all blood stages, and no additional information was obtained from this parasite line. In conclusion, our data highlight two new tools towards combination therapy. Hybrid molecules represent promising tools for the cure of co-infected individuals, while very potent antibiotics with a wide scope of activities could be useful in ACTs by eliminating resistant parasites and limiting transmission of both, resistances and disease.
A continuous arms race between the development of novel antibiotics and the evolution of corresponding resistance mechanisms in bacteria has been observed, since antibiotic agents like arsphenamines (e.g. Salvarsan, developed by Paul Ehrlich [1]), sulphonamides (e.g. Prontosil, Gerhard Domagk [2]) and penicillin (Alexander Fleming [3]) were first applied to effectively cure bacterial infections in the early 20th century. The rapid emergence of resistances in contrast to the currently lagging discovery of antibiotics displays a severe threat to human health. Some serious infectious diseases, such as tuberculosis or melioidosis, which were either thought to be an issue only in Third-World countries in case of tuberculosis, or regionally restricted with respect to melioidosis, are now on the rise to expand to other areas. In contrast, methicillin-resistant Staphylococcus aureus (MRSA) is already present in clinical setups all over the world and causes severe infections in immunocompromised patients. Thus, there is an urgent need for new and effective antimicrobial agents, which impair vital functions of the pathogen’s metabolism.
One central metabolic pathway is represented by the bacterial fatty-acid synthesis pathway (FAS II), which is essential for the synthesis of long and branched-chain fatty acids, as well as mycolic acids. These substances play a major role as modulating components of the properties of the most important protective barrier – the cell envelope. The integrity of the bacterial cell wall and the associated membrane(s) is crucial for cell growth and for protection against physical strain, intrusion of antibiotic agents and regulation of uptake of ions and other small molecules. Thus, this central pathway represents a promising target for antibiotic action against pathogens to combat infectious diseases. The last and rate-limiting step is catalysed by the trans-2-enoyl-ACP reductase (ENR) FabI or InhA (in mycobacteria), which has been demonstrated to be a valuable target for drug design and can be addressed, amongst others, by diphenyl ether (DPE) compounds, derived from triclosan (TCL) – the first one of this class which was discovered to bind to ENR enzymes [4, 5].
Based on this scaffold, inhibitors containing different combinations of substituents at crucial positions, as well as a novel type of substituent at position five were investigated regarding their binding behaviour towards the Burkholderia pseudomallei and Mycobacterium tuberculosis ENR enzymes bpFabI and InhA, respectively, by structural, kinetic and in-vivo experiments. Generally, substitution patterns modulate the association and dissociation velocities of the different ENR inhibitors in the context of the two-step slow-onset binding mechanism, which is observed for both enzymes. These alterations in the rapidity of complex formation and decomposition have a crucial impact on the residence time of a compound and hence, on the pharmacokinetic properties of potential drug candidates. For example, the substituents at the 2’-position of the DPE scaffold influence the ground- and transition state stability during the binding process to bpFabI, whereas 4’-substituents primarily alter the transition state [6]. The novel triazole group attached to the 5-position of the scaffold, targeting the hydrophobic part of the substrate-binding pocket in InhA, significantly enhances the energy barrier of the transition state of inhibitor binding [7] and decelerates the association- as well as the dissociation processes. Combinations with different substituents at the 2’-position can enhance or diminish this effect, e.g. by ground-state stabilisation, which will result in an increased residence time of the respective inhibitor on InhA.
Further structural investigations carried out in this work, confirm the proposed binding mode of a customised saFabI inhibitor [8], carrying a pyridone moiety on the DPE scaffold to expand interactions with the protein environment. Structural and preliminary kinetic data confirm the binding of the same inhibitor to InhA in a related fashion. Comparisons with structures of the ENR inhibitor AFN-1252 [9] bound to ENR enzymes from other organisms, addressing a similar region as the pyridone-moiety of the DPE inhibitor, suggest that also the DPE inhibitor bears the potential to display binding to homologues of saFabI and InhA and may be optimised accordingly.
Both of the newly investigated substituents, the pyridone moiety at the 4’-position as well as the 5-triazole substituent, provide a good starting point to modify the DPE scaffold also towards improved kinetic properties against ENR enzymes other than the herein studied and combining both groups on the DPE scaffold may have beneficial effects. The understanding of the underlying binding mechanism is a crucial factor to promote the dedicated design of inhibitors with superior pharmacokinetic characteristics.
A second target for a structure-based drug-design approach is the interaction surface between ENR enzymes and the acyl-carrier protein (ACP), which delivers the growing acyl chain to each distinct enzyme of the dissociated FAS-II system and presumably recognises its respective interaction partner via electrostatic contacts. The interface between saACP and saFabI was investigated using different approaches including crosslinking experiments and the design of fusion constructs connecting the ACP and the FabI subunits via a flexible linker region of varying lengths and compositions. The crosslinking studies confirmed a set of residues to be part of the contact interface of a previously proposed complex model [10] and displayed high crosslinking efficiency of saACP to saFabI when mutated to cysteine residues. However, crystals of the complex obtained from either the single components, or of the fusion constructs usually displayed weak diffraction, which supports the assumption that complex formation is highly transient. To obtain ordered crystals for structural characterisation of the complex it is necessary to trap the complex in a fixed state, e.g. by a high-affinity substrate attached to ACP [11], which abolishes rapid complex dissociation. For this purpose, acyl-coupled long-residence time inhibitors might be a valuable tool to elucidate the detailed architecture of the ACP-FabI interface. This may provide a novel basis for the development of inhibitors that specifically target the FAS-II biosynthesis pathway.
With 9.6 million new cases and 1.5 million deaths in 2014, tuberculosis (TB) is alongside with AIDS the most deadly infection. Foremost, the increased prevalence of resistant strains of M. tuberculosis among the TB-infected population represents a serious thread. Hence, in the last decades, novel drug targets have been investigated worldwide. So far a relatively unexplored target is the cell wall enzyme β-ketoacyl-ACP-synthase “KasA”, which plays a crucial role in maintaining the membrane impermeability and hence the cell ability to resist to the immune response and drug therapy. KasA is a key enzyme in the fatty acid synthase “FAS-II” elongation cycle, responsible for the extension of the growing acyl chain within the biosynthesis of precursors for the most hydrophobic constituents of the cell wall – mycolic acids. Design of the novel KasA inhibitors, performed in the research group of Prof. Sotriffer by C. Topf and B. Schaefer, was based on the recently published crystal structure of KasA in complex with its known inhibitor thiolactomycin (TLM). Considering the essential ligand-enzyme interactions, a pharmacophore model was built and applied in the virtual screening of a modified ZINC database. Selected hits with the best in silico affinity data have been reported by Topf and Schaefer.
In this work, two of the obtained hits were synthesized and their structure was systematically varied. First, a virtual screening hit, chromone-2-carboxamide derivative GS-71, was modified in the amide part. Since the most of the products possessed a very low solubility in the aqueous buffer medium used in biological assays, polar groups (nitro, succinamidyl and trimethyl-amino substituent in position 6 of the chromone ring or hydroxyl group on the benzene ring in the amide part have been inserted to the molecule. Further variations yielded diaryl ketones, diaryl ketone bearing a succinamidyl substituent, carboxamide bearing a methylpiperazinyl-4-oxobutanamido group and methyl-malonyl ester amides. Basically, the essential structural features necessary for the ligand-enzyme interactions have been maintained. The latter virtual screening hit, a pyrimidinone derivative VS-8 was synthesized and the structure was modified by substitution in positions 2, 4, 5 and 6 of the pyrimidine ring. Due to autofluorescence, detected in most of the products, this model structure was not further varied.
Simultaneously, experiments on solubilization of the first chromone-2-carboxamides with cyclodextrins, cyclic oligosacharides known to form water-soluble inclusion complexes, were performed. Although the assessed solubility of the chromone 3b/DIMEB (1:3) mixture exceeded 14-fold the intrinsic one, the achieved 100 µM solubility was still not sufficient to be used as a stock solution in the binding assay. The experiments with cyclodextrin in combination with DMSO were ineffective. Owing to high material costs necessary for the appropriate cyclodextrin amounts, the aim focused on structural modification of the hydrophobic products.
Precise structural data have been obtained from the solved crystal structures of three chromone derivatives: the screening hit GS-71 (3b), its trimethylammonium salt (18) and 6-nitro-substituted N-benzyl-N-methyl-chromone-2-carboxamide (9i). The first two compounds are nearly planar with an anti-/trans-rotamer configuration. In the latter structure, the carboxamide bridge is bent out of the chromone plane, showing an anti-rotamer, too. Considering the relatively low partition coefficient of compound 3b (cLogP = 2.32), the compound planarity and correlating tight molecular packing might be the factors significantly affecting its poor solubility.
Regarding the biological results of the chromone-based compounds, similar structure-activity correlations could be drawn from the binding assay and the whole cell activity testing on M. tuberculosis. In both cases, the introduction of a nitro group to position 6 of the chromone ring and the presence of a flexible substituent in the amide part showed a positive effect. In the binding study, the nitro group at position 4 on the N-benzyl residue was of advantage, too. The highest enzyme affinity was observed for N-(4-nitrobenzyl)-chromone-2-carboxamide 4c (KD = 34 µM), 6-nitro substituted N-benzyl-chromone-2-carboxamide 9g (KD = 40 µM) and 6‑nitro-substituted N-(4-nitrobenzyl)-chromone-2-carboxamide 9j (KD = 31 µM), which could not be attributed to the fluorescence quenching potential of the nitro group. The assay interference potential of chromones, due to a covalent binding on the enzyme sulfhydryl groups, was found to be negligible at the assay conditions. Moderate in vivo activity was detected for 6‑nitro-substituted N-benzyl-chromone-2-carboxamide 9g and its N-benzyl-N-methyl-, N‑furylmethyl-, N-cyclohexyl- and N-cyclohexylmethyl derivatives 9i, 9d, 9e, 9f, for which MIC values 20 – 40 µM were assessed. Cytotoxicity was increased in the N‑cyclohexylmethyl derivative only. None of the pyrimidine-based compounds showed activity in vivo. The affinity of the model structure, VS-8, surpassed with KD = 97 µM the assessed affinity of TLM (KD = 142 µM).
Since for the model chromone compound GS-71 no reliable KasA binding data could be obtained, a newly synthesized chromone derivative 9i was docked into the KasA binding site, in order to derive correlation between the in silico and in vitro assessed affinity. For the 6‑nitro-derivative 9i a moderate in vivo activity on M. tuberculosis was obtained. The in silico predicted pKi values for TLM and 9i were higher than the corresponding in vitro results, maintaining though a similar tendency, i.e., the both affinity values for compound 9i (pKi predicted = 6.64, pKD experimental = 4.02) surpassed those obtained for TLM (pKi predicted = 5.27, pKD experimental = 3.84). Nevertheless, the experimental pKD values are considered preliminary results.
The binding assay method has been improved in order to acquire more accurate data. Owing to the method development, limited enzyme batches and solubility issues, only selected compounds could be evaluated. The best hits, together with the compounds active on the whole cells of M. tuberculosis, will be submitted to the kinetic enzyme assay, in order to confirm the TLM-like binding mechanism. Regarding the in vivo testing results, no correlations could be drawn between the predicted membrane permeability values and the experimental data, as for the most active compounds 9e and 9f, a very low permeability was anticipated (0.4 and 0.7 %, respectively). Further biological tests would be required to investigate the action- or transport mode.
Bakterielle und parasitäre MIP-Proteine stellen wichtige Virulenzfaktoren dar, deren Inhibition das Überleben der Erreger sowie deren Penetration in menschliche Zellen stark einschränken kann. In dieser Arbeit standen die MIP-Proteine von Burkholderia pseudomallei (Auslöser der Melioidose) und Legionella pneumophila (Legionärskrankheit) im Fokus. Außerdem wurde das MIP-Protein von Trypanosoma cruzi (Chagas-Krankheit) untersucht. Die strukturverwandten humanen FKB-Proteine FKBP12 und FKBP52 sind relevante „off-targets“, wie Experimente mit Knockout-Mäusen gezeigt haben.
Ziel dieser Arbeit war die Verbesserung von bekannten MIP-Inhibitoren im Hinblick auf ihre Affinität und Selektivität für MIP-Proteine gegenüber den beiden genannten FKB-Proteinen bei gleichzeitig verbesserter Löslichkeit, mit Hilfe von in silico Methoden. Ausgangspunkt waren hierbei zwei von Dr. Christina Juli und Dr. Florian Seufert entwickelte Leitstrukturen, welche ein Pipecolinsäuregrundgerüst aufweisen. Diese Referenzliganden beinhalten einen 3,4,5-Trimethoxyphenylring (TMPR, vgl. Ref_t) bzw. einen Pyridinylring (Ref_p).
Beim Vergleich von insgesamt 32 MIP- und FKB-Proteinen konnten in zwei Loop-Bereichen, welche 50er bzw. 80er Loop genannt werden, relevante Unterschiede in der Aminosäuresequenz identifiziert werden. Die Nummerierung bezieht sich stets auf FKBP12. Diese Unterschiede ließen sich zum Design von vergleichsweise selektiv an MIP-Proteine bindenden Molekülen nutzen.
Der 50er Loop ist in nahezu allen MIP-Proteinen (jedoch nicht in BpsMIP) im Vergleich zu den FKB-Proteinen um zwei Aminosäuren verkürzt. Dadurch befindet sich das Proteinrückgrat von LpnMIP (Gln49) und TcrMIP (Arg49) näher am Zentrum der Bindetasche (definiert als Ile56, welches durch die Pipecolinsäureesterfunktion der Liganden adressiert wird). MD-Simulationen der beiden Apoproteine belegten, dass die geringere Distanz nicht durch Artefakte beim Modellieren der Strukturen bedingt ist. Aufbauend auf dieser Erkenntnis wurde gezeigt, dass der Pyridinylring von Ref_p eine Wasserstoffbrücke zu Gln49 ausbildet. Experimentell wurde dieser Befund durch eine entsprechende chemische Verschiebung der Aminosäure im NMR-Experiment von Dr. Kristian Schweimer bestätigt. Durch Überbrückung des Pipecolinsäurerings (Ligand 6bp) konnte die Wasserstoffbrücke in MD-Simulationen weiter stabilisiert werden. Durch Rechnungen zur Abschätzung der freien Bindungsenthalpien (mittels LIE und MM/GBSA) wurde eine erhöhte Affinität von 6bp im Vergleich zu Ref_p in LpnMIP ermittelt.
Im Laufe der Arbeit wurde anhand von pIC50-Werten, welche von Dr. Mathias Weiwad bestimmt wurden, erkannt, dass Liganden mit Pyridinylring oftmals eine bessere Affinität in LpnMIP aufweisen als die entsprechenden Liganden mit TMPR. Durch MD Simulationen wurde nachgewiesen, dass der TMPR in LpnMIP nur schwer an der in den anderen Proteinen bevorzugten Position binden kann. Grund hierfür ist die Mutation einer Aminosäure (zu Pro57) in diesem Bereich von LpnMIP: Diese verfügt über eine wenig flexible Seiten-kette, an welche sich der TMPR auf Grund seiner Rigidität nicht anpassen kann, was die Interaktion zwischen Protein und Ligand stört. Der Pyridinylring von Ref_p ist hiervon nicht betroffen, da er bevorzugt an einer anderen Stelle (Gln49, s. o.) bindet.
Der 80er Loop weist in vielen MIP-Proteinen deutlich hydrophobere Aminosäuren auf als in FKB-Proteinen. Von besonderem Interesse ist die Position 90, da hier in BpsMIP und LpnMIP sterisch weniger anspruchsvolle Aminosäuren (Val, Pro) vorliegen als in den bei-den FKB-Proteinen (Ile, Lys). Dieser Unterschied wurde mit kleinen hydrophoben Substituenten am Phenylring der Liganden adressiert. Bereits im Docking zeigten sich die positiven Effekte der para-Substitution durch Halogenatome oder eine Methylgruppe. Die von Dr. Mathias Weiwad und Dr. Mirella Vivoli ermittelten pIC50- bzw. pKi-Werte bestätigten diesen Trend. Zugleich nahm die Affinität zu FKBP12 deutlich ab. Bei der Untersuchung der Referenzliganden sowie deren Chlor- und Bromderivate in MD-Simulationen zeigte sich, dass der Phenylring der Liganden in den MIP-Proteinen bevorzugt in Richtung des 80er Loops orientiert ist; in den FKB-Proteinen liegt er hingegen um etwa 110° gedreht vor und kann somit schlechter mit der Bindetasche interagieren. Besonders ausgeprägt ist dieser Effekt in FKBP12. Basierend auf diesen Ergebnissen wurde der Phenylring durch einen 4-Bromo-1H-imidazol-2-ylsubstituenten ersetzt (Ligand 8ap). Dieser ist in der Lage, in der erwarteten Orientierung im Bereich des 80er Loops von BpsMIP zu binden und gleichzeitig eine stabile Wasserstoffbrücke zu Asp37 auszubilden. Hieraus resultiert für den Liganden eine deutlich höhere Affinität in LIE- und MM/GBSA-Rechnungen; in FKBP12 blieb sie auf Grund der dort instabilen Interaktion unverändert.
Die berechneten Energien können unmittelbar für einen relativen Vergleich verschiedener Liganden in einer Bindetasche verwendet werden. Für die Vorhersage von pKi- bzw. pIC50-Werten in den verschiedenen Proteinen ist eine Kalibrierung gegen die gemessenen Affinitäten erforderlich. Dies wurde für BpsMIP durchgeführt, indem eine lineare Korrelation zwischen den pKi- bzw. pIC50-Werten und den mit MM/GBSA ermittelten Energien aufgestellt wurde. Für LIE wurde auf publizierte Werte von Lamb et al. zurückgegriffen. Die berechneten Affinitäten stimmen für die bereits getesteten Inhibitoren gut mit den experimentellen pKi- und pIC50-Werten überein. Anhand der Modelle werden für 8ap Werte vorhergesagt, die besser als die experimentellen Affinitäten bekannter Liganden sind.
Idealerweise können auch aus den Scores, die durch Docking erhalten werden, bereits Rückschlüsse auf die Affinitäten der Liganden gezogen werden. Für die untersuchten Proteine war dies, auf Grund des engen Bereichs der experimentell ermittelten pKi- und pIC50-Werte, nicht mit hinreichender Richtigkeit möglich. Um die Scores dennoch für die Beurteilung neuer Liganden verwenden zu können, wurden logistische Regressionsmodelle erstellt. Anhand dieser kann abgeschätzt werden, ob ein Molekül in BpsMIP submikromolare Affinität aufweist. Die Richtigkeit dieser Vorhersagemodelle konnte durch die Berücksichtigung dreier weiterer Deskriptoren (Konfiguration am Stereozentrum der Pipecolinsäure, Molekulargewicht und logD-Wert) deutlich verbessert werden, wobei die AUC der entsprechenden ROC-Kurven Werte bis zu 0.9 erreichte. Diese Modelle können für die Postprozessierung eines Dockings angewendet werden, um die vielversprechendsten Kandidaten zu identifizieren und anschließend in rechnerisch anspruchsvolleren MD-Simulationen genauer zu untersuchen.
Mit dieser Arbeit wurde zur Weiterentwicklung der Leitstrukturen Ref_t und Ref_p beigetragen. Viele der getesteten Derivate wiesen deutlich verbesserte Löslichkeit bei gleichbleibender Affinität auf. Ferner wurden erstmalig detailliert die Unterschiede in den Bindetaschen zwischen 32 MIP- und FKB-Proteinen evaluiert. Hiervon wurden fünf in MD-Simulationen als Apoprotein und im Komplex mit verschiedenen Inhibitoren verglichen. Anhand dieser Simulationen wurde nachgewiesen, dass jeweils eine Aminosäure in BpsMIP und LpnMIP im Vergleich zum wichtigsten „off-target“ FKBP12 selektiv durch eine Wasserstoffbrücke adressiert werden kann. Durch LIE- und MM/GBSA-Rechnungen konnte gezeigt werden, dass in diesen hochkonservierten Bindetaschen eine bedeutende Modulation der Affinität zugunsten von BpsMIP möglich ist.
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.
Kovalente Inhibition stellt einen effektiven Weg dar, die Verweildauer des Liganden innerhalb einer Bindetasche zu erhöhen. In dieser Arbeit wurden theoretische Methoden angewendet, um die Reaktivität und den nichtkovalenten Zustand vor der Reaktion zu modellieren. Im Rahmen einer Fallstudie zu Cathepsin K wurden nichtkovalente Modelle von kovalenten Inhibitoren generiert. Für verschiedene Komplexe aus Cathepsin K und einem kovalent gebundenem Liganden wurde der Zustand vor der Reaktion modelliert und dessen Stabilität im Rahmen einer klassischen MD-Simulation überprüft. Die Stabilität des Warheads in der Bindetasche hing hauptsächlich vom gewählten Protonierungszustand der katalytischen Aminosäuren ab. Für eine Reihe von Inhibitoren der ChlaDUB1 wurde ein Protokoll aus quantenmechanischen Rechnungen genutzt, um die Reaktivität verschiedener Warheads abzuschätzen. Die erhaltenen Aktivierungsenergien korrelierten mit experimentell bestimmten Raten zur Inaktivierung des Enzyms. Im Rahmen eines Wirkstoffdesign-Projektes zur Deubiquitinase USP28 wurden von unpublizierten Kristallstrukturen ausgehend erste Docking-Experimente durchgeführt. Es konnte gezeigt werden, dass ein literaturbekannter Inhibitor von USP28 mit einem Warhead so modifiziert werden kann, dass die reaktive Einheit in direkter Nachbarschaft zu einem Cystein positioniert wird. Für diese Warheads wurden ebenfalls quantenmechanische Rechnungen zur Bestimmung der Aktivierungsenergie durchgeführt. Um besser nachvollziehen zu können, warum bei einem Photoswitch-Inhibitor der Butyrylcholin-Esterase der cis-Zustand des Moleküls besser inhibiert als der trans-Zustand, wurde eine Docking-Studie des Zustandes vor der Reaktion durchgeführt. Es konnte ein qualitatives Modell aufgestellt werden, das zeigt, dass der trans-Zustand aufgrund seiner längeren Form mit wichtigen Aminosäuren am Eingang der Bindungstasche kollidiert.
The high infection rates and recent emergence of extremely drug resistant forms of
Mycobacterium tuberculosis pose a significant challenge for global health. The NADH-
dependent enoyl-ACP-reductase InhA of the type II mycobacterial fatty acid biosynthesis
pathway is a well-validated target for inhibiting mycobacterial growth. InhA has been
shown to be inhibited by a variety of compound series. Prominent classes of InhA
inhibitors from literature include diaryl ethers, pyrrolidine carboxamides and arylamides
which can be subjected to further development. Despite the progress in this area, very
few compounds are in clinical development phase. The present work involves a detailed
computational investigation of the binding modes and structure-based optimisation of
pyrrolidine carboxamides as InhA inhibitors.
With substituents of widely varying bulkiness, the pyrrolidine carboxamide dataset
presented a challenge for prediction of binding mode as well as affinity. Using advanced
docking protocols and in-house developed pose selection procedures, the binding modes
of 44 compounds were predicted. The poses from docking were used in short molecular
dynamics (MD) simulations to ascertain the dominant binding conformations for the
bulkier members of the series. Subsequently, an activity-based classification strategy
could be developed to circumvent the affinity prediction problems observed with this
dataset. The prominent motions of the bound ligand and the active site residues were
then ascertained using Essential Dynamics (ED). The information from ED and literature
was subsequently used to design a total of 20 compounds that were subjected to extensive
in-silico evaluations. Finally, the molecular determinants of rapid-reversible binding of
pyrrolidine carboxamides were investigated using long MD simulations.
Die AAA+ ATPase p97 ist ein essenzielles Protein, das an einer Vielzahl zellulärer Prozesse beteiligt ist und eine Schlüsselrolle in der Protein-Homöostase spielt. Die funktionale Diversität von p97 beruht auf der Interaktion zahlreicher unterschiedlicher Kofaktoren, die vorwiegend an die N-Domäne von p97 binden. Aufgrund seiner Bedeutung in der Regulierung diverser physiologischer und pathologischer Prozesse stellt p97 eine interessante Zielstruktur für die Entwicklung neuer Wirkstoffe dar, die insbesondere in der Krebstherapie von Bedeutung sein könnte. Bekannte p97-Inhibitoren greifen vor allem die ATPase-Funktion des Proteins an. Ein neuer pharmakologischer Ansatz stellt die Inhibierung der Kofaktorbindung an die N-Domäne dar. Ein solcher Protein-Protein-Interaktionsinhibitor wäre nicht nur von therapeutischem Interesse, sondern hätte auch einen besonderen Nutzen für die Entschlüsselung molekularer und zellulärer Funktionen von p97-Kofaktoren. In dieser Arbeit wurde ein fragmentbasierter Ansatz für die Identifizierung von chemischen Startstrukturen für die Entwicklung eines Protein-Protein- Interaktionsinhibitors verfolgt. Als Zielstruktur wurde die SHP-Bindestelle in der N-Domäne gewählt. Die Identifizierung von Liganden erfolgte sowohl durch computergestützte Methoden (insbesondere virtuelles Screening und Molekulardynamik-Simulationen) als auch experimentell durch biophysikalische Techniken (wie Biolayer-Interferometrie, Röntgenstrukturanalyse und ligandbasierte NMR-Techniken). Die Grundlage des computerbasierten Designs stellte eine Analyse der bekannten Kristallstrukturen der p97-Komplexe mit den SHP-Motiven der Kofaktoren UFD1 und Derlin-1 dar. Darüber hinaus dienten Molekulardynamik-Simulationen der Analyse der Wassereigenschaften innerhalb der SHP-Bindestelle. Darauf aufbauend wurden verschiedene Pharmakophormodelle entwickelt, die die Grundlage des im Anschluss durchgeführten virtuellen Screenings und Dockings bildeten. Anhand der Ergebnisse von Molekulardynamik-Simulationen wurden zehn Verbindungen für die experimentelle Validierung ausgewählt. Hiervon konnten zwei Fragmente in STD-NMR- und Biolayer-Interferometrie-Experimenten als Liganden bestätigt werden. In einem parallel durchgeführten biophysikalischen Fragmentscreening mittels Biolayer-Interferometrie wurden unter mehr als 650 Verbindungen 22 identifiziert, die an die N-Domäne binden. 15 dieser Fragmente wurden durch einen orthogonalen STD-NMR-Assay bestätigt. Fünf dieser Verbindungen zeigten Affinitäten mit KD-Werten kleiner 500μMund günstigen Ligandeffizienzen. Des Weiteren konnte die Bindungskinetik und Affinität des in der Literatur als p97-Inhibitor berichteten Naturstoffes Xanthohumol bestimmt und eine Bindung an die N-Domäne bestätigt werden. Zur Identifizierung möglicher Bindestellen dieser fünf Fragmente wurden mixed-solvent Molekulardynamik-Simulationen durchgeführt. Diese ergaben, dass alle Verbindungen die SHP-Bindestelle in der N-Domäne adressieren. Die Regionen fielen mit hot spots der Kofaktorwechselwirkungen zusammen und stellen somit mögliche Ankerpunkte für die Weiterentwicklung dar. Für zwei Fragmente konnten die postulierten Bindestellen mittels Röntgenstrukturanalyse bzw. STD-NMR-Messungen an p97-Alanin-Mutanten bestätigt werden. Die erhaltene Röntgenstruktur ist die erste p97-Struktur, die ein gebundenes Fragment an der N-Domäne zeigt.
Entwicklung, Validierung und Anwendung einer interpretierbaren und alignment-freien 4D-QSAR Methodik
(2006)
Die vorliegende Arbeit beschreibt die Entwicklung, Validierung und erfolgreiche Anwendung der interpretierbaren 4D-QSAR Methodik xMaP. Die neue Methode benötigt weder die Auswahl des vermuteten bioaktiven Konformers noch eine Überlagerung der Moleküle im Raum, sie ist also alignment-frei. xMaP ist invariant gegenüber Rotation, Translation und kodiert die Flexibilität der Moleküle. Dadurch wird der Einfluss durch den Benutzer praktisch ausgeschaltet.
Dockingbasierte Ansätze zählen zu den wichtigsten Komponenten im virtuellen Screening. Sie dienen der Vorhersage der Ligandposition und -konformation in der Bindetasche sowie der Abschätzung der Bindungsaffinität zum Protein. Bis heute stellt die korrekte Identifizierung proteingebundener Ligandkonformationen ein noch nicht vollständig gelöstes Problem für Scoringfunktionen dar. Der erste Teil der vorliegenden Arbeit ist daher der Entwicklung computergestützter Methoden zur Bewertung von Docking-Lösungen gewidmet.
Der Fokus eines ersten Teilprojektes lag auf der Berücksichtigung der Absättigung vergrabener Wasserstoffbrückenakzeptoren (HBA) und -donoren (HBD) bei der Bewertung von Docking-Lösungen. Nicht-abgesättigte vergrabene HBA und HBD stellen einen der Bindungsaffinität abträglichen Beitrag dar, der bis dato aufgrund fehlender Struktur- bzw. Affinitätsdaten in Scoringfunktionen vernachlässigt wird. Im Rahmen der vorliegenden Arbeit wurde auf der Basis einer detaillierten Untersuchung zur Häufigkeit vergrabener nicht-abgesättigter HBA und HBD in hochaufgelösten Protein-Ligand-Komplexen des Hartshorn-Datensatzes eine empirische Filterfunktion („vnaHB“-Filterfunktion) entwickelt, die unerwünschte Ligandbindeposen erkennt und von der Bewertung mittels Scoringfunktionen ausschließt. Der praktische Nutzen der empirischen Filterfunktion wurde für die Scoringfunktionen SFCscore und DSX anhand vorgenerierter Docking-Lösungen des Cheng-Datensatzes untersucht. Die Häufigkeitsuntersuchung zeigt, dass eine Absättigung vergrabener polarer Gruppen in Protein-Ligand-Komplexen für eine hochaffine Protein-Ligand-Bindung notwendig ist, da vergrabene nicht-abgesättigte HBA und HBD nur selten auftreten. Eine vollständige Absättigung durch entsprechende Proteinpartner wird für ca. 48 % der untersuchten Komplexe beobachtet, ca. 92 % weisen weniger als drei hauptsächlich schwache, nicht-abgesättigte HBA bzw. HBD (z. B. Etherfunktionen) auf. Unter Einbeziehung von Wassermolekülen in die Häufigkeitsanalyse sind sogar für ca. 61 % aller Komplexe alle wasserstoffbrückenbindenden Gruppen abgesättigt. Im Gegensatz zu DSX werden für SFCscore nach Anwendung der empirischen Filterfunktion erhöhte Erfolgsraten für das Auffinden einer kristallnahen Pose (≤ 2.0 Å Abweichung) unter den am besten bewerteten Docking-Posen erzielt. Für die beste SFCscore-Funktion (SFCscore::229m) werden Steigerungen dieses als „Docking Power“ bezeichneten Kriteriums für die Top-3-Posen (Erfolgsrate für die Identifizierung einer kristallnahen 2.0 Å Pose unter den besten drei Docking-Lösungen) von 63.1 % auf 64.2 % beobachtet.
In einem weiteren Teilprojekt wurden repulsive Protein-Ligand-Kontakte infolge sterischer Überlappungen der Bindungspartner bei der Bewertung von Docking-Lösungen berücksichtigt. Die adäquate Einbeziehung solcher repulsiver Kontakte im Scoring ist für die Identifizierung proteingebundener Ligandkonformationen entscheidend, jedoch aufgrund fehlender Affinitäts- bzw. Strukturdaten problematisch. Im Rahmen der vorliegenden Arbeit wurde auf der Basis des Lennard-Jones-Potentiales des AMBER-Kraftfeldes zunächst ein neuer Deskriptor zur Beschreibung repulsiver Kontakte („Clash“-Deskriptor) entwickelt und zur Untersuchung der Häufigkeit ungünstiger Protein-Ligand-Kontakte in hochaufgelösten Protein-Ligand-Komplexen des Hartshorn-Datensatzes herangezogen. Eine aus der Häufigkeitsverteilung abgeleitete empirische Filterfunktion („Clash“-Filterfunktion) wurde anschließend der Bewertung von Docking-Lösungen des Cheng-Datensatzes mittels der Scoringfunktionen SFCscore und DSX vorgeschaltet, um unerwünschte Ligandbindeposen auszuschließen. Die Häufigkeitsuntersuchung zeigt, dass vorwiegend schwache repulsive Kontakte in Protein-Ligand-Komplexen auftreten. So werden in 75 % der Komplexe des Hartshorn-Datensatzes abstoßende Potentiale unter 0.462 kcal/mol beobachtet. Zwar betragen die ungünstigen Beiträge pro Komplex für 50 % aller Strukturen ca. 0.8 kcal/mol bis 2.5 kcal/mol, jedoch können diese auf Ungenauigkeiten der Kristallstrukturen zurückzuführen sein bzw. durch günstige Protein-Ligand-Wechselwirkungen kompensiert werden. Die Anwendung der „Clash“-Filterfunktion zeigt signifikante Verbesserungen der Docking Power für SFCscore. Für die beste SFCscore-Funktion (SFCscore::frag) werden Steigerungen der Erfolgsraten für das Auffinden einer kristallnahen Pose unter den drei am besten bewerteten Docking-Lösungen von 61.4 % auf 86.9 % erzielt, was an die Docking Power der bis dato besten Scoringfunktionen aus der Literatur (z. B. DSX, GlideScore::SP) heranreicht (Docking Power (DSX): 92.6 %; Docking Power (GlideScore::SP): 86.9 %). Die „Clash“-Filterfunktion allein ist auch der Kombination der „Clash“- und der „vnaHB“-Filterfunktion überlegen.
Ein weiterer Schwerpunkt der vorliegenden Arbeit wurde auf die Einbeziehung von Decoy-Daten (Struktur- und Affinitätsdaten schwach affiner und inaktiver Liganden) im Zuge der Entwicklung computergestützter Methoden zur Bewertung von Docking-Lösungen gelegt. Dadurch soll eine adäquate Berücksichtigung ungünstiger Beiträge zur Bindungsaffinität ermöglicht werden, die für die Richtigkeit und Zuverlässigkeit ermittelter Vorhersagen essentiell ist. In der vorliegenden Arbeit wurden binäre Klassifizierungsmodelle zur Bewertung von Docking-Lösungen entwickelt, die die Einbeziehung von Decoy-Daten ohne die Verfügbarkeit von Affinitätsdaten erlauben. Der Random-Forest-Algorithmus (RF), SFCscore-Deskriptoren, der neu entwickelte „Clash“-Deskriptor, und die Decoy-Datensätze von Cheng und Huang (Trainingsdaten) bilden die Grundlage des leistungsfähigsten Klassifizierungsmodells. Der praktische Nutzen des „besten“ RF-Modells wurde nach Kombination mit der Scoringfunktion DSX anhand der Docking Power für das Auffinden einer kristallnahen Pose auf Rang 1 am unabhängigen Cheng-/Huang- (Komplexe, die nicht in den Trainingsdaten enthalten sind) und CSAR-2012-Testdatensatz untersucht. Gegenüber einer alleinigen Anwendung von DSX werden an beiden Testdatensätzen weitere Verbesserungen der Docking Power erzielt (Cheng-/Huang-Testdatensatz: DSX 84.24 %, RF 87.27 %; CSAR-2012-Testdatensatz: DSX 87.93 %, RF 91.38 %). Das „beste“ Modell zeichnet sich durch die zuverlässige Vorhersage richtig-positiver Docking-Lösungen für einige wenige Komplexe aus, für die DSX keine kristallnahe Ligandkonformation identifizieren kann. Ein visueller Vergleich der jeweils am besten bewerteten RF- und DSX-Pose für diese Komplexe zeigt Vorteile des RF-Modells hinsichtlich der Erkennung für die Protein-Ligand-Bindung essentieller Wechselwirkungen. Die Untersuchung der Bedeutung einzelner SFCscore-Deskriptoren für die Klassifizierung von Docking-Lösungen sowie die Analyse der Misserfolge nach Anwendung des Modells geben wertvolle Hinweise zur weiteren Optimierung der bestehenden Methode. Hinsichtlich der zu bewertenden Eigenschaften ausgeglichenere Trainingsdaten, Weiterentwicklungen bestehender SFCscore-Deskriptoren sowie die Implementierung neuer Deskriptoren zur Beschreibung bis dato nicht-berücksichtigter Beiträge zur Bindungsaffinität stellen Ansatzpunkte zur Verbesserung dar.
Der zweite Teil der vorliegenden Arbeit umfasst die Anwendung dockingbasierter Methoden im Rahmen der Entwicklung neuer Inhibitoren des „Macrophage Infectivity Potentiator“-(MIP)-Proteins von Legionella pneumophila und Burkholderia pseudomallei.
Das MIP-Protein von Legionella pneumophila stellt einen wichtigen Virulenzfaktor und daher ein attraktives Zielprotein für die Therapie der Legionellose dar. Im Rahmen der vorliegenden Arbeit erfolgten systematische Optimierungen des Pipecolinsäure-Sulfonamides 1, des bis dato besten niedermolekularen MIP-Inhibitors (IC50 (1): 9 ± 0.7 µM). Nach Hot-Spot-Analysen der Bindetasche wurden Docking-Studien zur Auswahl aussichtsreicher Kandidaten für die Synthese und Testung auf MIP-Inhibition durchgeführt. Die Ergebnisse der Hot-Spot-Analysen zeigen günstige Wechselwirkungsbereiche für Donorgruppen und hydrophobe Substituenten in meta-Position sowie Akzeptorgruppen in para-Position des Benzylringes von 1 auf. Die Einführung einer Nitrofunktion in para-Position des Benzylringes von 1 (2h) resultiert in einer erhöhten MIP-Inhibition (IC50 (2h): 5 ± 1.5 µM), was wahrscheinlich auf die Ausbildung einer zusätzlichen Wasserstoffbrücke zu Gly116 zurückzuführen ist. Selektivitätsverbesserungen gegenüber dem strukturverwandten humanen FKBP12-Protein werden insbesondere für das para-Aminoderivat von 1 (2n) erzielt (Selektivitätsindex (1): 45, Selektivitätsindex (2n): 4.2; mit Selektivitätsindex = IC50 (MIP)/IC50 (FKBP12)). Der Ersatz des hydrophoben Trimethoxyphenylrestes von 1 durch einen Pyridinring (2s) führt zu einer verbesserten Löslichkeit bei vergleichbarer MIP-Inhibition.
Das MIP-Protein von Burkholderia pseudomallei spielt eine wichtige Rolle in der Pathogenese der Melioidose und stellt daher ein attraktives Zielprotein für die Entwicklung neuer Arzneistoffe dar. In der vorliegenden Arbeit erfolgten Optimierungen des bis dato besten niedermolekularen MIP-Inhibitors 1. Ausgehend von einem Strukturvergleich von Burkholderia pseudomallei MIP mit Legionella pneumophila MIP und einer Hot-Spot-Analyse der Burkholderia pseudomallei MIP-Bindetasche wurden Docking-Studien zur Auswahl aussichtsreicher Kandidaten für die Synthese und Testung auf MIP-Inhibition durchgeführt. Der Strukturvergleich zeigt eine hohe Homologie beider Bindetaschen. Größere konformelle Änderungen werden lediglich für den von Ala94, Gly95, Val97 und Ile98 geformten Bindetaschenbereich beobachtet, was unterschiedliche Optimierungsstrategien für 1 erforderlich macht. Günstige Wechselwirkungsbereiche der Burkholderia pseudomallei MIP-Bindetasche finden sich einerseits für Donorgruppen oder hydrophobe Substituenten in para-Position des Benzylringes (Region A) von 1, andererseits für Akzeptor- bzw. Donorgruppen in para- bzw. meta-/para-Position des Trimethoxyphenylringes (Region B). Anhand von Docking-Studien konnten sowohl für Variationen in Region A als auch in Region B aussichtsreiche Kandidaten identifiziert werden. Initiale MIP-Inhibitionsmessungen der bis dato synthetisierten Derivate deuten auf erhöhte Hemmungen im Vergleich zu 1 hin. Der Ersatz des hydrophoben Trimethoxyphenylrestes von 1 durch einen Pyridinring führt auch hier zu vergleichbarer MIP-Inhibition bei verbesserter Löslichkeit. Derzeit sind weitere Synthesen und Testungen aussichtsreicher Liganden durch die Kooperationspartner geplant. Die Ergebnisse der Inhibitionsmessungen sollen deren Nutzen als MIP-Inhibitoren aufzeigen und wertvolle Informationen für weitere Zyklen des strukturbasierten Wirkstoffdesigns liefern.
Diese Dissertation beschreibt Methoden zur Lösung wichtiger anwendungsorientierter Aspekte des struktur- und ligandbasierten in-silico-Wirkstoffdesigns. Dabei liegt der Fokus auf der Entwicklung chemometrischer Verfahren und der Überprüfung ihrer Leistungsfähigkeit. Die vorgeschlagenen Algorithmen werden mit entsprechenden etablierten Techniken verglichen. Die folgenden Abschnitte fassen die Vorgehensweisen und Resultate in den einzelnen Projektbereichen zusammen. Identifizierung von Outliern. Die Untersuchung eines QSAR-Datensatzes mit dem Ziel der Outlier-Identifizierung wird in der Praxis häufig vernachlässigt. Dabei ist es offensichtlich, daß kein QSAR-Modell auf jede nur denkbare chemische Verbindung anwendbar sein kann. Vielmehr handelt es sich um empirische mathematische Modelle, die nur innerhalb jenes Datenraums Gültigkeit besitzen, der von den Trainingsobjekten aufgespannt wird. Daher ist jedes Modell auf gewisse Grenzen beschränkt, außerhalb derer eine verläßliche Vorhersage unmöglich ist. Die in dieser Arbeit entwickelte Methode ODD dient der Ermittlung dieser Grenzen und damit der Identifizierung von Outliern, also Objekten außerhalb des Anwendungsbereichs des Modells. Ziel der Entwicklung war ein nur auf den unabhängigen Variablen (X-Daten) basierendes Verfahren, das auch auf hochdimensionaleDatensätze anwendbar ist undweitestgehend auf den Eingriff des Benutzers (etwa die Definition von Grenzwerten) verzichtet. Ebenfalls wünschenswert war die Fähigkeit zur Identifikation von Inliern. Eine ausreichend hohe Geschwindigkeit sollte die Einsetzbarkeit im virtuellen Screening gewährleisten. Die Methode mußte der Überprüfung standhalten, den Vorhersagefehler eines Modells bei Vorhandensein extremer Outlier zu reduzieren, gleichzeitig aber unkritische Datensätze unbeeinflußt zu lassen. ODD basiert auf der Beurteilung der euklidischen Distanz eines Testobjekts zu seinem am nächsten benachbarten Trainingsobjekt. Der Schwellenwert für die Betrachtung eines Objekts als Outlier wird dabei aus der Verteilung der Nächster-Nachbar-Distanzen der Trainingsobjekte berechnet. Durch dieses intrinsische Maß ergibt sich die gewünschte Dimensionsunabhängigkeit und vor allem die automatische Anpassung des Grenzwerts an die Charakteristik des Kalibrierdatensatzes ohne Eingriff des Benutzers. Die Validierung zeigt, daß ODD extreme Outlier zuverlässig erkennt und sich gleichzeitig durch eine im Vergleich zu anderen gebräuchlichen Verfahren geringere Anzahl falsch positiver Identifizierungen auszeichnet. Ensemble-Techniken. In einer vergleichenden Studie wurde die Leistungsfähigkeit verschiedener Ensemble-Techniken hinsichtlich ihres Einflusses auf den Vorhersagefehler untersucht. Dazu wurden umfangreiche Simulationen anhand mehrerer realer QSAR-Datensätze durchgeführt. Die Verwendung von Ensembles (d. h. einer Sammlung vielerModelle, diemit geringfügigmanipulierten Varianten des Trainingsdatensatzes kalibriert wurden) wirkt sich im allgemeinen positiv auf den Vorhersagefehler (RMSEP) aus. Diese Reduzierung des RMSEP wurde hier ermittelt und für verschiedenen Ansätze zur Ensemble-Generierung verglichen. Insgesamt betrachtet erwiesen sich die Methoden der konvexen Pseudodaten und des Baggings als die effektivsten Verfahren zur Ensemble-Generierung, da sie den Vorhersagefehler am deutlichsten verbesserten. Die konvexen Pseudodaten wurden erstmalig zur Erzeugung von Ensembles in der QSAR-Analyse eingesetzt; sie werden als neuer Standard zur Reduzierung des RMSEP bei QSAR-Problemen vorgeschlagen, die Regressionsmodelle auf Basis von latenten Variablen verwenden. Darüber hinaus bieten die Studien eine Abschätzung dermit Hilfe von Ensembles zu erzielenden Reduktion des Vorhersagefehlers bei typischen QSAR-Datensätzen. Virtuelles Screening. Beim virtuellen Screening handelt es sich um eine Technik zum Durchsuchen großer (virtueller)Molekülbibliotheken—oftmehrere Millionen Verbindungen — nach den aussichtsreichsten Wirkstoffkandidaten. Dies kann sowohl durch strukturbasierte als auch mit Hilfe ligandbasierter Verfahren geschehen. Es wurden umfangreiche Simulationen anhand sechs verschiedener Targets und einer Bibliothek von mehr als 90 000 Molekülen durchgeführt, um das Potential strukturbasierter (Docking mit FLEXX) und ligandbasierter (Ähnlichkeitssuchemitmehreren Referenzen) Verfahren zu vergleichen. Darüber hinauswurde durch Berechnung von Interaktionsfingerprints eineMöglichkeit geschaffen, die Information der beiden sonst getrennten Herangehensweisen zu kombinieren. Um den Einfluß des Klassifizierungsalgorithmus zu untersuchen, wurden verschiedene statistische Methoden zur Datenauswertung herangezogen. Als Bewertungskriterium für die Leistungsfähigkeit eines Verfahrens diente jeweils die Anzahl der wiedergefundenen aktiven Moleküle in der simulierten Screeningdatenbank. Die Resultate führen zu dem Schluß, daß ligandbasierte Verfahren, die einfacher einzusetzen sind aber mehr a-priori -Information benötigen, dem strukturbasierten virtuellen Screening hinsichtlich der Datenbankanreicherung überlegen sind. Weiterhin konnte gezeigt werden, wie nutzbringend die Zusammenführung von strukturbasierter Information und solcher über das Interaktionsmuster bekanntermaßen aktiver Verbindungen für die Erhöhung der Wiederfindungsrate ist. Bei der Datenanalyse stellte sich heraus, daß im Mittel bestimmte statistische Methoden (minimale euklidische Distanz ED/Min bzw. Tanimoto-Ähnlichkeit der Integer-Fingerprints Int/Min) zu bevorzugen sind. Kovalentes Docking von Cathepsin-Inhibitoren. Die Cysteinproteasen Cathepsin B und L sind interessante pharmakologische Targets. Geeignete Inhibitoren stammen u. a. aus der Strukturklasse der Aziridine. Ein nukleophiler Angriff des Cysteinrests des Enzyms auf den elektrophilen Aziridinring führt hier zur Ausbildung einer kovalenten Ligand-Rezeptor-Bindung. Praktisch alle erhältlichen Dockingprogramme konzentrieren sich jedoch auf nicht-kovalente Ligand-Rezeptor-Interaktionen und lassen kein uneingeschränktes kovalentes Docking zu. Daher wurde für FLEXX ein Dockingprotokoll entworfen, das den entscheidenden nicht-kovalenten Zustand vor Ausbildung der kovalenten Bindung simulieren kann. Auf dieseWeise konnte untersucht werden, ob sich die Reaktionszentren von Ligand und Enzym ausreichend nahe für die Ausbildung einer kovalenten Bindung kommen. Der vorgestellte Ansatz läßt sich leicht auf andere kovalente Ligand-Rezeptor- Systeme übertragen und bietet somit eine breite Anwendbarkeit. Weiterhin wurde die Parametrisierung der in FLEXX vorgesehenen Interaktionsgeometrien an die strukturellen Eigenheiten der zu dockenden Aziridide angepaßt. Diese weisen nämlich formal eine Amidbindung auf, deren geometrische und elektronische Eigenschaften jedoch deutlich von den Werten eines typischen Amids abweichen. Die Ergebnisse der Dockingstudien liefern wertvolle Einblicke für das Verständnis der Selektivität der untersuchten Liganden bezüglich Cathepsin B beziehungsweise L. Umgekehrt erbringt die gute Übereinstimmung der FLEXX-Resultate mit den experimentell bestimmten Inhibitionskonstanten den Nachweis für die Validität des verwendeten Dockingprotokolls.
Im Mittelpunkt dieser Arbeit stand die Entwicklung neuer Wirkstoffe gegen Tuberkulose, einer schwerwiegenden bakteriellen Infektionskrankheit, die am häufigsten die Lunge befällt. Die Entwicklung neuer Arzneistoffe gegen diese Erkrankung ist immens wichtig, da nach Angaben der WHO weltweit jährlich über 1 Million Menschen an den Folgen der Tuberkulose sterben, derzeit kein effizienter Impfstoff zur Verfügung steht und sich die Therapiemöglichkeiten auf wenige Arzneistoffe beschränken. Zudem steigt weltweit das Auftreten von arzneistoff- und totalresistenten Tuberkuloseformen. Tuberkulose wird vorwiegend durch das Mycobacterium tuberculosis erregt. Eine Besonderheit des M. tuberculosis stellt die mykobakterielle Zellwand dar, da diese durch einen hohen Anteil an Fettsäuren besonders wachsartig und dick ist. Die mykobakterielle Fettsäuresynthese unterscheidet sich signifikant von der Synthese eukaryotischer Fettsäuren. Daher besteht die Möglichkeit, Inhibitoren der mykobakteriellen Fettsäuresynthese als effektive und selektive neue Antituberkulotika zu entwickeln. Zielstruktur dieser Arbeit ist KasA (β-Ketoacylsynthase), ein Enzym der mykobakteriellen Fettsäuresynthese II, das die Kondensation zwischen der wachsenden Fettsäurekette und Malonyl-ACP katalysiert. Ein literaturbekannter KasA-Inhibitor ist Thiolactomycin, ein Thiolacton-Derivat mit einer schwachen inhibitorischen Aktivität (IC50: 242 µM; Kd 226 µM), für den eine KasA-Komplexstruktur verfügbar ist. Ziel der Arbeit war es, mittels computergestützten Wirkstoffdesigns neue Leitstrukturen für KasA-Inhibitoren zu entwickeln und davon abgeleitet Substanzbibliotheken kleiner Moleküle zu synthetisieren. Zur Bestimmung der In-vitro-Aktivitäten sollte KasA exprimiert und ein Assay etabliert werden. Theoretische und experimentelle Affinitäten sollten anschließend analysiert und bewertet werden. Zur Identifizierung neuer potenzieller KasA-Inhibitoren wurde mit Hilfe des Thiolactomycin-Bindemodus ein Pharmakophor-Modell erstellt. In diesem wurden die essentielle Wasserstoffbrücke zwischen den Histidinen und dem Carbonyl-Sauerstoff des Thiolactonrings, zwei hydrophobe Bereiche und ein verbindendes Strukturelement definiert und das Volumen des Pharmakophors begrenzt. Das Screening der Datenbank erfolgte mit GOLD4.0 und GOLDscore. Zur Identifizierung der 16 aussichtsreichsten Verbindungen wurden Rescorings mit ChemScore und sfc_score290m durchgeführt, sowie verschiedene physikochemische Deskriptoren und der errechnete Bindungsmodus einbezogen. Ausgewählte Verbindungen des Screenings wurden synthetisiert. Weitere Variationen wurden durch Einführung von Substituenten und Bromierung und Nitrierung der Grundgerüste erhalten. Zur biologischen Testung dieser Verbindungen konnte KasA in M. smegmatis exprimiert werden. Die Reinigung des Proteins erfolgte mittels Affinitäts- und Größenausschlusschromatographie. Affinitätswerte an KasA konnten mit einem Fluoreszenzassays bestimmt werden, da in jedem KasA-Monomer vier Tryptophane zur intrinsischen Fluoreszenz beitragen. Die Bindung eines Inhibitors in die TLM-Bindetasche führte zum Quenching der Fluoreszenz von KasA und konnte unter Berücksichtigung von Verdünnungs- und inneren Filtereffekten zur Berechnung der Dissoziationskonstante Kd herangezogen werden. Die In-vitro-Untersuchungen der Inhibitoren von KasA zeigten im Vergleich zu TLM eine Verbesserung der Affinität bis zu einem Faktor von 11, die beste Verbindung war das Nitroisatin-Derivat 2l (22.1 µM). Einen Hinweis auf Hemmung des Wachstums von Mykobakterien war für die Verbindungen 2e (5-Nitro-1-phenethyl-2,3-indolindion) und 3a (5,7-Dibrom-1-(4-chlorbenzyl)indolin-2,3-dion) ersichtlich. Die übrigen Verbindungen zeigten keine Aktivität, was dadurch bedingt sein kann, dass sie Substanzen nicht lipophil genug sind (clogP-Werte zwischen 1 und 3), um die mykobakterielle Zellwand zu durchdringen. Analog dem Docking im Rahmen des virtuellen Screenings wurde ein Docking mit GOLD4.0 und GOLDscore für die Substanzbibliothek durchgeführt. Verglichen mit den In-vitro-Affinitäten konnte eine gute Übereinstimmung in der Differenzierung der Substanzklassen gefunden werden. Da kleine Moleküle mit großer biologischer Aktivität zu bevorzugen sind, wurde die „ligand efficiency“, die inhibitorische Potenz unabhängig vom Molekulargewicht, für die Verbindungen berechnet. Für die Substanzbibliothek wurde eine gute Korrelation von „ligand efficiency“ und GOLDscore pro Schweratom erzielt (R2=0.65), beste Substanzgruppen waren monoalkylierte Uracil- und Isatin-Derivate. Der beste Wert wurde für das Isatin-Derivat 1a erzielt. Mit den erarbeiteten theoretischen und experimentellen Ergebnissen und den etablierten Methoden bietet diese Arbeit eine wichtige Grundlage, um erste „hits“ von KasA-Inhibitoren zu neuen Leitstrukturen für Wirkstoffe gegen Mycobakterium tuberculosis zu entwickeln.
Die vorliegende Arbeit beschreibt die Durchführung computergestützter Untersuchungen an den Wildtyp-Kristallstrukturen dieses Enzyms – einerseits zur Suche nach neuen Inhibitoren im Rahmen von virtuellen Screening- (VS-) Studien, andererseits zur Charakterisierung der strukturellen Flexibilität mit Hilfe von Molekulardynamik- (MD-) Simulationen. Für ein erstes VS wurde zunächst eine Datenbank von mehreren Millionen kommerziell erhältlichen Verbindungen mit Arzneistoff-ähnlichen physikochemischen Eigenschaften erstellt. Als Ausgangspunkt der Screening-Studie diente ein Teildatensatz, welcher mit Hilfe eines auf dem nativen Bindemodus des KasA-Inhibitors Thiolactomycin (TLM) beruhenden Pharmakophormodells erhalten wurde. Diese Verbindungen wurden in die Bindetasche von KasA gedockt und die Qualität der erhaltenen Posen in einem Re- und Konsensus-Scoring-Verfahren bewertet. Schließlich wurden 14 Substanzen käuflich erworben und im Rahmen einer Fluoreszenz-Bindungsstudie experimentell getestet. Für sechs Moleküle war gegenüber KasA eine schwache, zu TLM vergleichbare Aktivität zu verzeichnen. Eine zweite VS-Studie befasste sich mit der Bewertung der Bindungsaffinität synthetisch leicht zugänglicher Derivate von GS95, einem gegenüber dem KasA-Wildtyp aktiven Molekül mit einer 1-Benzyluracil-Grundstruktur. Anhand geeigneter Synthesebausteine wurde eine virtuelle Datenbank von insgesamt 16 Derivaten erstellt, welche in die Bindetasche des Enzyms gedockt wurden. Für die vorhergesagten Bindemodi erfolgte dann eine Abschätzung der freien Bindungsenthalpie. Nach einer Bewertung der Orientierungen auf Basis der errechneten ΔG-Werte sowie einer visuellen Analyse wurden schließlich elf Verbindungen synthetisiert und im Fluoreszenz-Experiment getestet, wobei für alle Uracilderivate eine Aktivität zu beobachten war. Die Kd-Werte fallen jedoch ähnlich hoch aus wie bei GS95. Zur Untersuchung der Strukturdynamik des KasA-Wildtyps wurden drei MD-Simulationen des homodimeren Proteins von je 15 ns Länge durchgeführt. Mit Hilfe von 2D-RMSD-Berechnungen und einer hierarchischen Clusteranalyse wurden aus den drei Simulationen insgesamt zehn repräsentative Snapshots entnommen, welche die im Rahmen der Simulationszeit von insgesamt 90 ns produzierte strukturelle Vielfalt der Bindetasche von KasA wiedergeben. Wie die Analysen zeigen, wird ein dualer Charakter hinsichtlich der Flexibilität der unmittelbaren Taschenreste beobachtet. Hierbei zeigt Phe404 eine besonders ausgeprägte strukturelle Vielfalt; diese Beobachtung deckt sich mit der gatekeeper-Rolle der Aminosäure zwischen der Malonyl-Bindetasche und dem Acyl-Bindekanal, welcher für die Unterbringung der wachsenden Fettsäurekette im Enzym während der Katalyse verantwortlich zeichnet. Darüber hinaus erklärt die hohe Flexibilität von Phe404 möglicherweise die recht schwache Bindungsaffinität von TLM gegenüber dem Wildtyp von KasA, da die gatekeeper–Aminosäure nur in der geschlossenen Form einen stabilisierenden Effekt auf den Liganden ausübt. Besondere Bedeutung kommt hierbei einem Wassermolekül zu, welches als eine Art molekularer Schalter für die Flexibilität von Phe404 betrachtet werden kann und somit die Fixierung von TLM in der Bindetasche maßgeblich beeinflusst. Des Weiteren wurden innerhalb der Tasche hohe Besetzungsraten für je ein Wassermolekül identifiziert. Die aus den MD-Simulationen gewonnenen Erkenntnisse wurden anschließend zur Aufstellung von Empfehlungen für das Design neuartiger KasA-Inhibitoren verwendet. Des Weiteren wurde die Dynamik des oben erwähnten Acyl-Bindekanals, welcher sich aus den Aminosäuren 115-147 zusammensetzt, näher charakterisiert. Hierbei wurden die Reste 115-119 und insbesondere Leu116 als zweiter gatekeeper identifiziert, welcher die Öffnung des Acyl-Bindekanals zur Oberfläche des Proteins reguliert und somit eine entscheidende Funktion bei der Unterbringung der langkettigen Fettsäuresubstrate übernimmt. Schließlich wurden zwei repräsentative Bindetaschenkonformationen aus den MD-Simulationen hinsichtlich einer Verwendung in strukturbasierten VS-Studien näher untersucht. Mit Hilfe von hot spot Analysen und unter Berücksichtigung oben genannter Empfehlungen für das Design neuartiger KasA-Inhibitoren wurden verschiedene Pharmakophormodelle erstellt, welche nach Durchsuchung der zu Anfang dieses Kapitels erwähnten virtuellen Moleküldatenbank zwischen 149 und 420 verschiedene hit-Strukturen lieferten. Folglich scheint eine Adressierung der beiden Konformationen durch arzneistoffartige Verbindungen prinzipiell möglich. Unter den erhaltenen Verbindungen herrscht eine hohe strukturelle Vielfalt; außerdem unterscheiden sich diese im Allgemeinen deutlich von den Molekülen aus den vorangegangenen VS Studien, was das Potential der beiden Bindetaschenkonformationen zur Identifizierung von potentiellen KasA-Inhibitoren mit neuartigen Grundstrukturen zum Ausdruck bringt.
In this thesis, computational structure-based design approaches were employed to target the HIV-1 integrase and the macrophage infectivity potentiator (MIP) of Legionella pneumophila. The thesis yields valuable information about the mechanism of action of a known class of integrase inhibitors and a novel approach towards enzyme inhibition, which still is mainly unaddressed in current integrase research. For the MIP enzyme, two small-molecule MIP inhibitors were discovered. The computational studies of HIV-1 integrase have provided valuable information for IN inhibitor design. Docking experiments supported the hypothesis that the well-known diketo acid inhibitors enter the IN active site not as free ligands, but rather as metal complexes. These results help to reveal the mechanism of action of this important class of IN inhibitors.To give an impulse for the development of a novel class of inhibitors, a new strategy towards IN inhibition was introduced: An alternative binding site, the dimerization interface of an IN catalytic core domain monomer, was explored for inhibitor design. The lack of structural data of the free monomer was overcome by extensive MD studies. Snapshots derived from the MD simulation were used as protein input structures in a docking study with the inhibitory peptide YFLLKL to reveal its potential binding mode. The docking procedure showed that the peptidic ligand binds to a dimerization interface conformation which shows a Y-shaped binding site.. The next step was to address this protein conformation with small, non-peptidic molecules. The first strategy towards finding small-molecule interface binders was to create a pharmacophore model with hydrophobic features and shape constraints, aiming to find molecules with a good complementarity to the Y-shaped dimerization interface. Virtual screening yielded a total of 10 compounds, which all displayed good shape complementarity and favorable hydrophobic interactions. Unfortunately, none of the compounds showed a reproducible inhibitory activity in biological assays. Some doubts remain about the validity of the assay results: The use of BSA was critical, since it is not unlikely that BSA “intercepted” the hydrophobic candidate compounds. The first strategy towards finding small-molecule dimerization inhibitors was reconsidered: In the second approach, the satisfaction of hydrogen bonding residues at the dimerization interface, was of major interest. Two pharmacophore models were employed, which retrieved several hundred hit molecules. However, docking of these molecules showed that still many hydrogen bonding groups of the protein remained unaddressed by the ligands. Eventually, after visual inspection, only eight molecules were selected as candidate compounds for further testing (results pending). This small “yield” underlines the difficulties in finding interface binders: The IN dimerization interface is a peculiar target with frequently alternating basic, acidic, and hydrophobic residues. It is not a well-ordered binding site with continuous hydrophobic areas and distinct hydrogen bond donors / acceptors. Other protein-protein interfaces show such well-ordered binding sites. Accordingly, the peculiarity of the IN dimerization interface, in addition to the delicate task of disrupting protein-protein interactions at all, makes the development of IN dimerization inhibitors very challenging. For MIP, the studies revealed two experimentally validated MIP inhibitors, which significantly reduce MIP enzymatic activity. To our knowledge, no small-molecule MIP inhibitor has been reported in the literature so far. A detailed analysis of the available structural data of MIP and a comparison to the human PPIase counterpart, FKBP12, pointed out a conformational diversity among the MIP structures and a crucial difference between the two PPIases, which could be traced to mainly one residue (Tyr109). The detailed comparison of FKBP12 and MIP complex structures made it possible to give an explanation, why a ketoacyl-substituted pipecoline derivative most probably does not bind to MIP, but a sulfone-substituted pipecoline derivative does bind to MIP. Knowledge of Legionella MIP inhibitors could be transferred also to other organisms (e.g. trypanosoms), where homologous MIP proteins are also pathological factors.
\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.