TY - THES A1 - Schmidt, Thomas Christian T1 - Theoretical Investigations on the Interactions of Small Compounds with their Molecular Environments T1 - Theoretische Untersuchungen der Wechselwirkungen Kleiner Moleküle mit deren Molekularen Umgebungen N2 - Im ersten Teil dieser Arbeit wird eine Kombination theoretischer Methoden für die strukturbasierte Entwicklung neuer Wirkstoffe präsentiert. Ausgehend von der Kristallstruktur eines kovalenten Komplexes einer Modellverbindung mit dem Zielprotein wurde mit Hilfe von quantenmechanischen und QM/MM Rechnungen die genaue Geometrie des vorausgehenden nicht-kovalenten Komplexes betimmt. Letztere ist der bestimmende Faktor für die Reaktivität des Inhibitors gegenüber der katalytisch aktiven Aminosäure und damit für die Ausbildung einer kovalenten Bindung. Aus diesem Grund wurde diese Geometrie auch für die Optimierung der Substitutionsmusters des Ihnibitors verwendet, um dessen Affinität zum Zielenzyme zu verbessern ohne dass dieser seine Fähigkeit kovalent an das aktive Zentrum zu binden verliert. Die Optimierung des Substitutionsmuster wurde doch Methode des Molekularen Dockings unterstützt, das diese optimal dazu geeignet sind, Bindungsaffinitäten vorherzusagen, die durch eine Modifikation der chemischen Struktur entstehen. Eine Auswahl der besten Strukturen wurde anschließend verwendet, um zu überprüfen, ob die veränderten Moleküle noch genügen Reaktivität gegenüber dem Zielprotein aufweisen. Moleküldynamik Simulationen der neuen Verbindungen haben jedoch gezeigt, dass die veränderten Verbindungen nur so and das Protein binden, dass die Bilung eine kovalenten Bindung zum Enzym nicht mehr möglich ist. Daher wurden in einem weiteren Schritt die Modellverbindungen weiter modifiziert. Neben Änderungen im Substitutionsmuster wurde auch die chemische Struktur im Kern verändert. Die Bindungsaffinitäten wurde wieder mittels Docking überprüft. Für die besten Bindungsposen wurden wieder Simulationen zur Moleküldynamik durchgeführt, wobei diesmal die Ausbildung einer kovalenten Bindung zum Enzyme möglich erscheint. In einer abschließenden Serie von QM/MM Rechnungen unter Berücksichtigung verschiedener Protonierungszustände des Inhibitors und des Proteins konnten Reaktionspfade und zugehörige Reaktionsenergien bestimmt werden. Die Ergebnisse lassen darauf schließen, dass eines der neu entwickelten Moleküle sowohl eine stark verbesserte Bindungsaffinität wie auch die Möglichkeit der kovalenten Bindung an Enzyme aufweist. Der zweite Teil der Arbeit konzentriert sich auf die Umgebungseinflüsse auf die Elektronenverteilung eines Inhibitormodells. Als Grundlage dient ein vinylsulfon-basiertes Moekül, für das eine experimentell bestimmte Kristallstruktur sowie ein theoretisch berechneter Protein Komplex verfügbar sind. Ein Referendatensatz für diese Systeme wurde erstellt, indem der Konformationsraum des Inhibitors nach möglichen Minimumsstrukturen abgesucht wurde, welche später mit den Geometrien des Moleküls im Kristall und im Protein verglichen werden konnten. The Geometrie in der Kristallumgebung konnte direkt aus den experimentellen Daten übernommen werden. Rechnungen zum nicht-kovalenten Protein Komplex hingegen haben gezeigt, dass für das Modellsystem mehrere Geometrien des Inhibiors sowie zwei Protonierungszustände für die katalytisch aktiven Aminosäuren möglich sind. Für die Analyse wurden daher alle möglichen Proteinkomplexe mit der Kristallstruktur verglichen. Ebenso wurden Vergleiche mit der Geometrie des isolierten Moleküls im Vakuum sowie der Geometrie in wässriger Lösung angestellt. Für die Geometrie des Moleküls an sich ergab sich eine gute Übereinstimmung für alle Modellsysteme, für die Wechselwirkungen mit der Umgebung jedoch nicht. Die Ausbildung von Dimeren in der Kristallumgebung hat einen stark stablisierenden Effekt und ist einer der Gründe, warum dieser Kristall so gut wie keine Fehlordungen aufweist. In den Proteinkomplexen hingegen ergibt sich eine Abstoßung zwischen dem Inhibitor und einer der katalytisch aktiven Aminosäuren. Als Ursache für diese Abstoßung konnte die Einführung der Methylaminfunktion ausgemacht werden. Vermutlicherweise führt diese strukturelle Änderung auch dazu, dass der Modellinhibitor nicht in der Lage ist, so wie die Leitstruktur K11777 an das aktive Zentrum des Enzyms zu binden. N2 - In the first part of this work, a combination of theoretical methods for the rational design of covalent inhibitor is presented. Starting from the crystal structure of the covalent complex of a lead compound, quantum mechanical and QM/MM calculations were used to derive the exact geometry of the preceeding non-covalent enzyme inhibitor complex. The geometry of the latter mainly determines the reactivity of the inhibitor against its target enzyme concerning the formation of the covalent bond towards an active site residue. Therefore, this geometry was used as starting point for the optimization of the substitution pattern of the inhibitor such as to increase its binding affinity without loosing its ability to covalently bind to the target protein. The optimization of the chemical structure was supported by using docking procedures, which are best suited to estimate binding affinities that arise from the introduced changes. A screening of the novel substitution patterns resulted in a first generation of model compounds which were further tested for their reactivity against the target. Dynamic simulations on the novel compounds revealed that the orientation that compounds adopt within the active site are such that a covalent interaction with the enzyme is no longer possible. Hence, the chemical structure was further modified, including not only changes in the substituents but also within the core of the molecule. Docking experiments have been conducted to assure sufficiently high binding affinities and to obtain the most favored binding poses. Those have then again been used for dynamic simulations which resulted in structures, for which the bond formation process appeared feasible. A final series of QM/MM calculations considering various protonation states was computed to estimate the reaction energies for the covalent attachment of the inhibitor to the enzyme. The theoretical results indicate a reasonable high inhibition potency of the novel compounds. The second part concentrates on the environmental influences on the electron density of an inhibitor molecule. Therefore, a vinylsulfone-based model compound was selected for which an experimental crystal structure for the pure compound as well as a theoretically determined enzyme-inhibitor complex have been available. To provide reference data for the larger systems, the conformational space of the isolated molecule was screened for favorable geometries which were later compared to those within the crystal and protein surrounding. The geometry of the crystal structure could readily be taken from the experimental data whereas calculations on the protein complex revealed four potential non-covalent complexes exhibiting different arrangements of the molecule within the active site of the protein as well as two possible protonation states of the catalytic dyad. Hence, all four protein complexes have been compared to the crystal structure of the molecule as well as against the more favorable geometries of the isolated molecule being determined within vacuum or aqueous surrounding. Whereas the molecule itself was found to adopt comparable geometries within all investigated environments, the interactions pattern between the crystal surrounding and the protein differed largely from each other. The favorable formation of dimers within the crystal has a strong stabilizing effect and explains the extraordinarily good quality of the crystal. Within the protein however, repulsive forces have been found between the protein and the inhibitor. The origin of the repulsion could be traced back to effect of on of the substituents to the vinyl scaffold. The difference in the chemical structure in comparison to a well known inhibitor might also explain the experimentally found loss of activity for the model compound in comparison to K11777. KW - Theoretische Chemie KW - theoretical chemistry KW - electron density KW - inhibition KW - Elektronendichte KW - Inhibitor Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-127860 ER - TY - THES A1 - Heilos, Anna T1 - Mechanistic Insights into the Inhibition of Cathepsin B and Rhodesain with Low-Molecular Inhibitors T1 - Mechanistische Untersuchungen zur Inhibition von Cathepsin B und Rhodesain mit niedermolekularen Inhibitoren N2 - Cysteine proteases play a crucial role in medical chemistry concerning various fields reaching from more common ailments like cancer and hepatitis to less noted tropical diseases, namely the so-called African Sleeping Sickness (Human Arfican Trypanosomiasis). Detailed knowledge about the catalytic function of these systems is highly desirable for drug research in the respective areas. In this work, the inhibition mechanisms of the two cysteine proteases cathepsin B and rhodesain with respectively one low-molecular inhibitor class were investigated in detail, using computational methods. In order to sufficiently describe macromolecular systems, molecular mechanics based methods (MM) and quantum mechanical based method (QM), as well as hybrid methods (QM/MM) combining those two approaches, were applied. For Cathespin B, carbamate-based molecules were investigated as potential inhibitors for the cysteine protease. The results indicate, that water-bridged proton-transfer reactions play a crucial role for the inhibition. The energetically most favoured pathway (according to the calculations) includes an elimination reaction following an E1cB mechanism with a subsequent carbamylation of the active site amino acid cysteine. Nitroalkene derivatives were investigated as inhibitors for rhodesain. The investigation of structurally similar inhibitors showed, that even small steric differences can crucially influence the inhibition potential of the components. Furthermore, the impact of a fluorination of the nitroalkene inhibitors on the inhibition mechanism was investigated. According to experimental data measured from the working group of professor Schirmeister in Mainz, fluorinated nitroalkenes show – in contrast to the unfluorinated compounds – a time dependent inhibition efficiency. The calculations of the systems indicate, that the fluorination impacts the non-covalent interactions of the inhibitors with the enzymatic environment of the enzyme which results in a different inhibition behaviour. N2 - Cysteinproteasen spielen eine wichtige Rolle in der medizinischen Chemie. Nicht nur im Bereich bekannterer Krankheiten wie Krebs oder Hepatitis, sondern auch bezüglich weniger verbreiteter, tropischer Krankheiten wie der sogenannten afrikanischen Schlafkrankheit (Afrikanische Trypanosomiasis) haben diese Enzyme eine große Bedeutung. Im Bereich der Wirkstofffindung ist ein detailliertes Wissen über die katalytische Funktion der an einer Krankheit beteiligten Enzyme unabdingbar .In der vorliegenden Arbeit wurden die Inhibitionsmechanismen der beiden Cysteinproteasen Cathepsin B und Rhodesain in Verbindung mit zwei niedermolekularen Inhibitorklassen anhand theoretischer Berechnungen untersucht. Um die makromolekularen Systeme ausreichend genau beschreiben zu können, wurden neben molekularmechanischen (MM) und quantenmechanischen (QM) Ansätzen auch Hybridmethoden verwendet, welche beide Ansätze (QM/MM) verbinden. Für Cathepsin B wurden Carbamat-basierte Moleküle als potenzielle Inhibitoren der Cysteinprotease untersucht. Die Ergebnisse weisen darauf hin, dass wasser-verbrückte Protonentransferreaktionen eine entscheidende Rolle für die Inhibition spielen. Der laut den Rechnungen energetisch günstigste Mechanismus beinhaltet eine Eliminierungsreaktion nach einem E1cB Mechanismus gefolgt von der Carbamylierung der Aminosäure Cystein in der aktiven Tasche des Enzyms. Nitroalken-Derivate wurden als potenzielle Rhodesain Inhibitoren untersucht. Der Vergleich strukturell ähnlicher Verbindungen weist darauf hin, dass schon kleine sterische Veränderungen einen großen Einfluss auf das Inhibitionspotenzial der Nitroalkene haben können. Außerdem wurde der Einfluss einer Fluorierung der Inhibitoren anhand von Berechnungen untersucht. Messungen der Arbeitsgruppe von Prof. Schirmeister in Mainz zu fluorierten und unfluorierten Nitroalkenen zeigen, dass die fluorierten Verbindungen ein zeitabhängiges Inhibitionspotenzial in Rhodesain aufweisen. Die Berechnungen der Systeme deuten darauf hin, dass die Fluorierung die nicht-kovalenten Wechselwirkungen der Inhibitoren mit der enzymatischen Umgebung des Systems beeinflussen, was zu einem unterschiedlichen Inhibitionsverhalten führt. KW - Cysteinproteasen KW - Inhibitor KW - Mechanismus KW - Berechnung KW - Inhibition Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-178228 ER - TY - THES A1 - Paasche, Alexander T1 - Mechanistic Insights into SARS Coronavirus Main Protease by Computational Chemistry Methods T1 - Mechanistische Einblicke in die SARS Coronavirus Hauptprotease mit computerchemischen Methoden N2 - The SARS virus is the etiological agent of the severe acute respiratory syndrome, a deadly disease that caused more than 700 causalities in 2003. One of its viral proteins, the SARS coronavirus main protease, is considered as a potential drug target and represents an important model system for other coronaviruses. Despite extensive knowledge about this enzyme, it still lacks an effective anti-viral drug. Furthermore, it possesses some unusual features related to its active-site region. This work gives atomistic insights into the SARS coronavirus main protease and tries to reveal mechanistic aspects that control catalysis and inhibition. Thereby, it applies state-of-the-art computational methods to develop models for this enzyme that are capable to reproduce and interpreting the experimental observations. The theoretical investigations are elaborated over four main fields that assess the accuracy of the used methods, and employ them to understand the function of the active-site region, the inhibition mechanism, and the ligand binding. The testing of different quantum chemical methods reveals that their performance depends partly on the employed model. This can be a gas phase description, a continuum solvent model, or a hybrid QM/MM approach. The latter represents the preferred method for the atomistic modeling of biochemical reactions. A benchmarking uncovers some serious problems for semi-empirical methods when applied in proton transfer reactions. To understand substrate cleavage and inhibition of SARS coronavirus main protease, proton transfer reactions between the Cys/His catalytic dyad are calculated. Results show that the switching between neutral and zwitterionic state plays a central role for both mechanisms. It is demonstrated that this electrostatic trigger is remarkably influenced by substrate binding. Whereas the occupation of the active-site by the substrate leads to a fostered zwitterion formation, the inhibitor binding does not mimic this effect for the employed example. The underlying reason is related to the coverage of the active-site by the ligand, which gives new implications for rational improvements of inhibitors. More detailed insights into reversible and irreversible inhibition are derived from in silico screenings for the class of Michael acceptors that follow a conjugated addition reaction. From the comparison of several substitution patterns it becomes obvious that different inhibitor warheads follow different mechanisms. Nevertheless, the initial formation of a zwitterionic catalytic dyad is found as a common precondition for all inhibition reactions. Finally, non-covalent inhibitor binding is investigated for the case of SARS coranavirus main protease in complex with the inhibitor TS174. A novel workflow is developed that includes an interplay between theory and experiment in terms of molecular dynamic simulation, tabu search, and X-ray structure refinement. The results show that inhibitor binding is possible for multiple poses and stereoisomers of TS174. N2 - Das Schwere Akute Respiratorische Syndrom (SARS) wird durch eine Infektion mit dem SARS Virus ausgelöst, dessen weltweite Verbreitung 2003 zu über 700 Todesfällen führte. Die SARS Coronavirus Hauptprotease stellt ein mögliches Wirkstoffziel zur Behandlung dar und hat Modellcharakter für andere Coronaviren. Trotz intensiver Forschung sind bis heute keine effektiven Wirkstoffe gegen SARS verfügbar. Die vorliegende Arbeit gibt Einblicke in die mechanistischen Aspekte der Enzymkatalyse und Inhibierung der SARS Coronavirus Hauptprotease. Hierzu werden moderne computerchemische Methoden angewandt, die mittels atomistischer Modelle experimentelle Ergebnisse qualitativ reproduzieren und interpretieren können. Im Zuge der durchgeführten theoretischen Arbeiten wird zunächst eine Fehlereinschätzung der Methoden durchgeführt und diese nachfolgend auf Fragestellungen zur aktiven Tasche, dem Inhibierungsmechanismus und der Ligandenbindung angewandt. Die Einschätzung der quantenchemischen Methoden zeigt, dass deren Genauigkeit teilweise von der Umgebungsbeschreibung abhängt, welche als Gasphasen, Kontinuum, oder QM/MM Modell dargestellt werden kann. Letzteres gilt als Methode der Wahl für die atomistische Modellierung biochemischer Reaktionen. Die Vergleiche zeigen für semi-empirische Methoden gravierende Probleme bei der Beschreibung von Proton-Transfer Reaktionen auf. Diese wurden für die katalytische Cys/His Dyade betrachtet, um Einblicke in Substratspaltung und Inhibierung zu erhalten. Dem Wechsel zwischen neutralem und zwitterionischem Zustand konnte hierbei eine zentrale Bedeutung für beide Prozesse zugeordnet werden. Es zeigt sich, dass dieser „electrostatic trigger“ von der Substratbindung, nicht aber von der Inhibitorbindung beeinflusst wird. Folglich beschleunigt ausschließlich die Substratbindung die Zwitterionbildung, was im Zusammenhang mit der Abschirmung der aktiven Tasche durch den Liganden steht. Dies gibt Ansatzpunkte für die Verbesserung von Inhibitoren. Aus in silico screenings werden genauere Einblicke in die reversible und irreversible Inhibierung durch Michael-Akzeptor Verbindungen gewonnen. Es wird gezeigt, dass unterschiedlichen Substitutionsmustern unterschiedliche Reaktionsmechanismen in der konjugierten Additionsreaktion zugrunde liegen. Die vorangehende Bildung eines Cys-/His+ Zwitterions ist allerdings für alle Inhibierungsmechanismen eine notwendige Voraussetzung. Letztendlich wurde die nicht-kovalente Bindung eines Inhibitors am Beispiel des TS174-SARS Coronavirus Hauptprotease Komplexes untersucht. Im Zusammenspiel von Theorie und Experiment wurde ein Prozess, bestehend aus Molekulardynamik Simulation, Tabu Search und Röntgenstruktur Verfeinerung ausgearbeitet, der eine Interpretation der Bindungssituation von TS174 ermöglicht. Im Ergebnis zeigt sich, dass der Inhibitor gleichzeitig in mehreren Orientierungen, als auch in beiden stereoisomeren Formen im Komplex vorliegt. KW - SARS KW - Inhibitor KW - Enzym KW - Computational chemistry KW - Coronaviren KW - SARS KW - Protease KW - Mechanismus KW - Inhibitor KW - Computerchemie KW - SARS KW - protease KW - mechanism KW - inhibitor KW - computational chemistry Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-79029 ER -