Refine
Has Fulltext
- yes (4)
Is part of the Bibliography
- yes (4)
Document Type
- Doctoral Thesis (4)
Keywords
- QM/MM (4) (remove)
Institute
Two thematic complexes were addressed within this work. One part is related to improvements and new implementations into the CAST program package. Thereby the main focus laid on the delivery of a tool which can be used to characterize complex reactions and their mechanisms. But also within the new force field (FF) method (SAPT-FF) within the CAST program, several improvements were made. The second topic is related to the description of dye molecules and their spectral properties. The main focus within these studies was set on the influence of the environment on these properties. In the first topic improvements of the local acting NEB (nudged elastic band) methods were included and the number of available methods was extended. The initial pathway generation was improved by implementing the IDPP (image dependent pair potential) method and a new method was implemented for describing temperature dependent pathways. Additionally, improvements have been made to the optimization routines (global NEB). As a second part the Pathopt (PO) method was considerably improved. In the beginning of the work the original PO idea was used. In this approach one starts with a global optimization on one n-1 dimensional hyperplane which divides the reaction into two sub-areas for obtaining guesses of TSs (transition states). These found TS guesses were used to optimize to the ”true” TS. Starting from the optimized ones a relaxation to the next connected minima is done. This idea has been automatically implemented and extended to several number of hyperplanes. In this manner a group of pathsegments is obtained which needs to be connected, but within this work it was realized that such a procedure might be not very efficient. Therefore, a new strategy was implemented which is founded on the same constrained global optimization scheme (MCM) for which the user defines the number of hyperplanes generated. The number of such generated hyperplanes should be large enough
134
to describe the space between the concerning reactants in a sufficient way. The found minima are directly used to built up the reaction pathway. For this purpose a RMSD (root mean square deviation) criterion is used to walk along ways of minimal change from one to another hyperplane. To prove the implementations various test calculations were carried out and extensions included to prove the capabilities of the new strategy. Related to these tests a new strategy for applying the move steps in MCM (Monte Carlo with minimization) was realized which is also related to the question of the coordinates representation. We were able to show that the hopping steps in MCM can be improved by applying Cartesian steps in combination of random dihedral moves with respect to the constraint. In this way it was possible to show that a large variety of systems can be treated. An additional chapter shows the improvements of the SAPT-FF implementation and related test cases. It was possible to treat benzene dimer and cluster systems of different sizes consistently also in accordance with high level ab initio based approaches. Furthermore, we showed that the SAPT-FF with the right parameters outperforms the standard AMOEBA implementation which is the basis of the SAPT-FF implementation. In the last three chapters deal with the description of perlyene-based dyes. In the first smaller chapter ground state chemistry description of macro cycles of PBI (perylene bisimide) derivatives were investigated. Therefore, AFM (atomic force microscopy) based pictures were explained within our study. The methods to explain aggregation behavior in dependency of the ring size were MD simulations and configuration studies. The last two chapters deal with opto-electronic or photo-physical properties of PBI and PTCDA (perylene-3,4,9,10-tetracarboxylic dianhydride). In detail, we investigated the role of the environment and the aggregate or crystal surrounding by applying different models. In that way implicit and explicit solvation models, the size of aggregates and vibration motions were used. In the case of PBI the recent work is found on preliminary studies related to my bachelor thesis and extends it. It was shown that the direct influence of a polarizable surrounding, as well as explicit inclusion of solvent molecules on the overall description of the excitations and nature of the excited states is weaker as one might expect. However the inclusion of intra-molecular degrees of freedom showed a stronger influence on the state characteristics and can induce a change of the order of states within the dimer picture. For the PTCDA molecule the main focus was set on the description of the absorption spectrum of crystalline thin films. Related to this older works exist which already gave a description and assignment of the absorption band, but are based on different approaches compared to the one used in this work. We used the supermolecule ansatz, whereas the environment and different aggregate sizes were investigated. Within the dimer based approach we were able to show that using continuum solvation (IEFPCM/COSMO) based description for the environment the relative order of states remains unchanged. Similar to the PBI calculations the influence of the vibrational motions /distortions is larger. The simulation of the crystal environment by using QM/MM (quantum mechanics/molecular mechanics) approaches delivered that an asymmetric charge distribution might induce a localization of the excitation and a stronger mixing of states. For obtaining further insights we go beyond the dimer picture and aggregates of different sizes were used, whereas the simulations up to the octadecamer mono- and even dual-layer stack were carried out. Within these calculations it was shown that the H-coupling is dominating over a weaker J-coupling between different stacks. Additionally the calculations based on DFT (density functional theory) and semi-empirics showed that the lowest state in terms of energy are mostly of Frenkel type, whereas the higher lying states are CT ones which mix with embedded Frenkel type states. The first band of the absorption spectrum was explained by inclusion of vibrational motions within the stacks which induce an intensity gain of the first excited state. This intensity was not explainable by using the undistorted stacks. Also relaxations at the crystal surface might play a role, but are experimentally not explainable.
Nowadays, computational-aided investigations become an essential part in the chemical, biochemical or pharmaceutical research. With increasing computing power, the calculation of larger biological systems becomes feasible. In this work molecular mechanical (MM) and quantum mechanical approaches (QM) and the combination of both (QM/MM) have been applied to study several questions which arose from different working groups. Thus, this work comprises eight different subjects which deals with chemical reactions or proton transfer in enzymes, conformational changes of ligands or proteins and verification of experimental data.
This work firstly deals with reaction mechanisms of aromatic inhibitors of cysteine proteases which can be found in many organisms. These enzymes are responsible for various cancer or diseases as for example Human African Trypanosomiasis (HAT) or the Chagas disease. Aromatic SNAr-type electrophiles might offer a new possibility to covalently modify these proteases. Quantum mechanical calculations have been performed to gain insights into the energetics and possible mechanisms.
The next chapter also deals with Trypanosomiasis but the focus was set on a different enzyme. The particularity of Trypanosomiasis is the thiol metabolism which can also be modified by covalent inhibitors. In this context, the wild type and point mutations of the enzyme tryparedoxin have been investigated via molecular dynamic (MD) simulations to examine the influence of specific amino acids in regard to the inhibitor. Experimental data showed that a dimerization of the enzyme occurs if the inhibitor is present. Simulations revealed that the stability of the dimer decreases in absence of the inhibitor and thus confirms these experiments.
Further investigations concerning cysteine proteases such as cruzain and rhodesain have been conducted with respect to experimental kinetic data of covalent vinylsulfone inhibitors. Several approaches such as QM or QM/MM calculations and docking, MD or MMPBSA/MMGBSA simulations have been applied to reproduce these data. The utilization of force field approaches resulted in a qualitatively accurate prediction.
The kinase AKT is involved in a range of diseases and plays an important role in the formation of cancer. Novel covalent-allosteric inhibitors have been developed and crystallized in complex with AKT. It was shown that depending on the inhibitor a different cysteine residue is modified. To investigate these differences in covalent modification computational simulations have been applied.
Enoyl-(acyl carrier) (ENR) proteins are essential in the last step of the fatty acid biosynthesis II (FAS) and represent a good target for inhibition. The diphenylether inhibitor SKTS1 which was originally designed to target the ENR’s of Staphylococcus aureus was also crystallized in InhA, the ENR of Mycobacterium tuberculosis (TB). Crystal structures indicate a change of the inhibitor's tautomeric form. This subject was investigated via MD simulations. Results of these simulations confirmed the tautomerization of the inhibitor.
This work also deals with the development of a covalent inhibitor originating from a non-covalent ligand. The target FadA5 is an essential enzyme for the degradation of steroids in TB and is responsible for chronic tuberculosis. This enzyme was crystallized in complex with a non-covalent ligand which served as starting point for this study. Computations on QM or QM/MM level and docking and MD simulations have been applied to evaluate potential candidates.
The next chapter focuses on the modification of the product spectrum of Bacillus megaterium levansucrase, a polymerase which catalyzes the biosynthesis of fructans. The covalent modification of the wild type or mutants of the enzyme lead to an accumulation of oligosaccharides but also to polymers with higher polymerization degree. To understand these changes in product spectra MD simulations have been performed.
Finally, the proton transfer in catalytic cysteine histidine dyads was investigated. The focus was set on the influence of the relaxation of the protein environment to the reaction. Calculations of the enzymes FadA5 and rhodesain revealed that the preferred protonation state of the dyade depends on the protein environment and has an impact on the reaction barrier. Furthermore, the adaptation of the environment to a fixed protonation state was analyzed via MD simulations.
Die Entwicklung maßgeschneiderter Proteinliganden ist ein integraler Bestandteil unterschiedlicher wissenschaftlicher Disziplinen, wie z.B. Wirkstoffentwicklung. Die vorliegende Arbeit befasst sich mit der reversiblen Inhibition in Form von kovalent gebundenen Enzym-Inhibitor-Komplexen der humanen Golgi-alpha-Mannosidase II (GM II) und der Cysteinprotease Rhodesain. Beide Enzyme sind erfolgversprechende Targets in der Bekämpfung von zwei sehr unterschiedlichen Erkrankungen. Einerseits die Golgi-alpha-Mannosidase bei der Behandlung der Tumorprogression und andererseits die Cysteinprotease Rhodesain bei der Behandlung der Afrikanischen Schlafkrankheit. Die Arbeit an den zwei Enzymen unterteilt sich in zwei Teilprojekte.
Die Entwicklung von maßgeschneiderten kovalent-reversiblen Inhibitoren für die genannten Enzyme wurde im Rahmen eines in-house entwickeltes Protokolls zwecks des rationalen Designs kovalenter Inhibitoren, durchgeführt. Dieses Protokoll basiert auf einer sich gegenseitig unterstützenden Zusammenarbeit zwischen Theorie und Experiment. Die vorliegende Arbeit befasst sich mit den theoretischen Untersuchungen mit Hilfe der quantenmechanischen (QM) als auch mit Hilfe der kombinierten quantenmechanisch/molekülmechanischen (QM/MM) Methoden zu den genannten Enzymen. In einem ersten Schritt des Protokolls geht es um die Anwendung von Screeningverfahren. In einem Screening werden Leitstrukturen, zunächst in Lösung (Schritt I), für eine weitere Untersuchung im Enzym (Schritt II) evaluiert. So können die Inhibitoren, für die experimentelle Mess- oder theoretische Dockingdaten vorhanden sind, als eine Leitstruktur betrachtet werden. Durch das Screening unter Verwendung der quantenmechanischen (QM-Modell) Methode kann eine Reihe von Inhibitoren nach einem sich konsistent veränderndem Muster erstellt werden und auf Bindungsparameter hin untersucht werden (Schritt I). Diese Parameter sind Reaktionsenergien und Höhen der Reaktionsbarriere einer Inhibitionsreaktion. Reaktionsenergien werden in dieser Betrachtung quantenmechanisch innerhalb der Born-Oppenheimer-(BO)-Näherung und im Rahmen des Konzeptes der Potentialhyperflächen (PES) als relative Energien zwischen den optimierten Geometrien der Produkte und der Edukte auf einer Potentialhyperfläche für die Inhibitionsreaktion ermittelt. Die Höhen der Reaktionsbarrieren werden durch die relativen Energien zwischen den Geometrien der Edukte und der Zwischenstufen oder Übergangszustände abgeschätzt. Unter Inhibitionsreaktion wird eine chemische Reaktion verstanden, bei der eine kovalente Bindung zwischen dem Inhibitormolekül und den Aminosäuren in der aktiven Tasche eines Enzyms ausgebildet wird. Für den Schritt I werden die Aminosäuren der aktiven Tasche durch kleine Moleküle, wie Essigsäure und Methanthiol, angenähert. Die kovalent-reversiblen Inhibitoren sollten in dieser Betrachtung nur leicht exotherme Reaktionen mit den relativen Energien im Bereich -5 bis -10 kcal/mol aufweisen. Der experimentelle Teil liefert währenddessen die Synthese der neuen Inhibitoren und die Nachweise zur kovalenten Bindung mit Hilfe massenspektrometrischer Messungen (Schritt I). Die passenden Kandidaten aus dem ersten Schritt des Protokolls, d.h. Inhibitoren mit gewünschten Bindungsparametern, werden durch die QM/MM-Berechnungen im Enzym (Schritt II) und durch die experimentellen Messungen an den Enzym-Inhibitor-Komplexen in Assays (Schritt II) analysiert. Die Untersuchungen für die Stufe II des Protokolls umfassen die Berechnungen der Reaktionsprofile und Minimumenergiereaktionspfade für die chemischen Reaktionen von Inhibitoren im Zielenzym. Ein Pfad minimaler potentieller Energie, der zwei Minima (hier Edukt und Produkt) verbindet, stellt ein Reaktionsprofil für eine chemische Reaktion dar. In der vorliegenden Arbeit wird dies auch als Minimumenergiepfad (MEP) bezeichnet. Der Letztere lässt sich durch die Nudged Elastic Band (NEB)-Methode und mittels Potentialhyperflächen darstellen. Die Reversibilität der Inhibitoren wurde anhand der berechneten chemischen Reaktionen in Form von erstellten Reaktionsprofilen analysiert und diskutiert. Durch Protein-Ligand Docking (Schritt III) wird ein Screening von variierbaren Erkennungseinheiten der neuen Inhibitoren durchgeführt. Die Ergebnisse der Untersuchungen aus dem dritten Schritt liefern Hinweise zur Weiterentwicklung der ausgewählten Inhibitoren. Die letzte Stufe des in-house Protokolls besteht in der erneuten Untersuchung der optimierten Inhibitoren mit Hilfe von Theorie und Experiment (Schritt IV). Die theoretische Untersuchung anhand von QM/MM-Berechnungen überprüft, ob die Inhibitionsreaktion der reaktiven Kopfgruppe nach der Änderung der Erkennungseinheit des Inhibitors weiterhin effektiv und nach dem gleichen Mechanismus mit der aktiven Seite des Enzyms ablaufen kann. Die experimentelle Untersuchung liefert, ähnlich wie im Schritt II, die messbaren Ergebnisse der Inhibition in Hinblick auf die Bindungseigenschaften und die Entstehung der Nebenprodukte.
Die Untersuchungen am System Mannosidase GM II wurden in Zusammenarbeit mit den Arbeitskreisen von Prof.Dr. J. Seibel und Prof.Dr. T. Schirmeister durchgeführt. Die Leitstruktur zur Entwicklung des kovalent-reversiblen Inhibitors stellt der cyclische O,O-Acetal-Inhibitor (bestimmt anhand von Dockingexperimenten an der beta-L-Anhydrogulose durch Arbeitskreis Prof.Dr. J. Seibel) dar.
Die Ergebnisse der theoretischen Studie liefern für den ersten Schritt im Rahmen des Protokolls den geeigneten Kandidaten aus einer Menge von insgesamt 22 modellierten Inhibitoren für die reversible Inhibition der Mannosidase GM II durch die Ausbildung einer kovalenten Bindung. Hierzu zählen zunächst die thermodynamischen Modellberechnungen der Inhibitionsreaktion, welche die Reaktionsenergien für alle Kandidaten des Screenings liefern. Die Inhibitionsreaktion wird in diesem Schritt als Additionsreaktion von Essigsäure an den Inhibitor-Kandidaten modelliert. Für die Leitstruktur resultiert eine thermoneutrale Beschreibung der Reaktion mit Essigsäure und dient im Weiteren als Referenz. Der Inhibitor Nr.7 (der cyclische N,O-Acetal-Inhibitor) zeigt mit -7,7 kcal/mol eine leicht exotherme Reaktion und somit eine bessere Triebkraft der untersuchten Reaktion im Vergleich zur Referenz. Die beiden Inhibitoren wurden dann für Stufe 2 des Protokolls untersucht, in der eine Analyse der Reaktionsprofile im Enzym mit Hilfe der QM/MM-Methodik durchgeführt wurde. Die Ergebnisse des zweiten Teils der Studie zeigen, dass der cyclische N,O-Acetal-Inhibitor eine deutlich bessere Affinität zur aktiven Seite der GM II im Vergleich zu seiner Leitstruktur aufweisen sollte. Dies zeigt sich auch in der deutlich höheren Triebkraft der Inhibitionsreaktion von ca. -13 kcal/mol. Dieser Energiebeitrag ist klein genug, um eine Reversibilität der Inhibitionsreaktion gewährleisten zu können. Das bedeutet auch, dass der N,O-Acetal-Inhibitor im Vergleich zur Referenzstruktur eine deutlich stärkere Inhibition bedingen sollte. Berücksichtigt man dann noch, dass die Reaktion laut Berechnungen nur leicht exotherm sein sollte, erhält man die Möglichkeit einer reversibel stattfinden kovalenten Hemmung. Zusammenfassend liefert dieser Teil der Arbeit, der mit Hilfe der QM- und QM/MM-Berechnungen durchgeführt wurde, ein reaktives molekulares Gerüst mit den gewünschten Eigenschaften.
Durch die theoretischen Untersuchungen (MD-Simulationen am Enzym-Inhibitor-Komplex) konnte außerdem eine zur Komplexbildung geeignete Konformation der Leitstruktur sowie des neuen Inhibitors gefunden werden. Die reversibel agierenden Acetal-Inhibitoren befinden sich in der aktiven Tasche in einer energetisch höher liegenden Twist-Boot-Konformation und begünstigen mit zwei entstehenden Bindungen zum Zn2+-Ion die oktaedrische Koordination im Enzym. Als Teil dieser Arbeit wurden NEB-Berechnungen zur Bestimmung von Minimumenergiepfadendurchgeführt. Dies lieferte erweiterte Einblicke in der Berechnung von Reaktionsmechanismen jeweils auch in Kombination von 2- bzw. 3-dimensionalen Scans. Auch in der Beschreibung von Protonenübertragungsreaktionen nach Grotthus, die einem Umklappen der kovalenten Bindungen entsprechen, erhält man hier Geometrien für die Teilschritte und somit eine detaillierte Beschreibung des Vorgangs.
Der Mechanismus der Inhibition von GM II durch die Leitstruktur beinhaltet einen Wasser-katalysierten (oder auch Wasser-vermittelten) Ringöffnungsschritt in der Tasche des Enzyms. Die Testrechnungen zum Protontransfer haben gezeigt, dass der Protontransfer über ein oder mehrere Wassermoleküle unter Verwendung von Standard-PES-Berechnungen nicht spontan stattfindet. Die Berechnung des MEP durch das Erstellen einer 3-dimensionalen Potentialhyperfläche kann nur dann sinnvolle Ergebnisse liefern, wenn der Protontransfer vom Aspartat Asp341 zum Inhibitor über zwei Wassermoleküle explizit berücksichtigt wird. In diesem Fall ist die Berechnung der PES kein Standard und erfordert eine zusätzliche Variation der Bindungsabstände O-H der beteiligten Moleküle des Protontransfers. Die Details für die zusätzliche Variation der Bindungsabstände O-H bei der Berechnung der 3-dimensionalen PES haben die NEB-Berechnungen geliefert. Der NEB-Formalismus hat sich in der Beschreibung dieser komplexen Reaktionskoordinaten als besser geeignet erwiesen und wurde in dieser Arbeit aus diesem Grund hauptsächlich verwendet. Die Berechnung des Protonentransfers während einer Hemmungsreaktion durch zwei Wassermoleküle mit der NEB-Methode hat den MEP ermittelt, welcher zunächst nicht auf der Grundlage eines 3-dimensionalen Scans ermittelt werden konnte. Solche QM/MM-Rechnungen wurden im Rahmen des in-house Protokolls zum ersten Mal durchgeführt. Dieser Protontransfer ist mit dem Grotthus-Mechanismus konform und kann plausibel anhand einer Klapp-Mechanismus-Betrachtung nachvollzogen werden.
Mit Hilfe der NEB-Methode ist es möglich MEPs effektiv und relativ schnell zu ermitteln. Es werden sowohl die Geometrien entlang des Pfades wie auch die einzelnen relativen Energien erhalten. Zur Überprüfung der gefundenen Übergangszustände wurden die einzelnen Strukturen mit Hilfe der Normalmodenanalyse weiter untersucht und konnten verifiziert werden.
Die MEP-Berechnungen für den Inhibitor Nr.1 ermöglichen die Etablierung eines Protokolls zur Berechnung eines Reaktionspfades über mehrere Moleküle, welches anschließend zur Berechnung des MEP für den Inhibitor Nr.7 angewendet wird. Das Protokoll beinhaltet in seiner einfachen Form die Ermittlung der Two-End-Komponenten einer chemischen Reaktion - Geometrien von Reaktant und Produkt. Betrifft dies eine Reaktion, die über mehrere Moleküle, z.B. Wassermoleküle oder deren Netzwerk, stattfindet, wird die Aufgabe komplexer. In diesem Fall ist eine Berechnung mit Hilfe des NEB-Moduls wesentlich produktiver als die Charakterisierung mit Hilfe der 3-dimensionalen PES. Der Vorteil liegt in der kollektiven Beschreibung der Reaktionskoordinaten, sodass die entscheidenden Reaktionskoordinaten und Variablen für die Durchführung von Scans nicht einzeln bestimmt werden müssen. Dennoch kann es hier bei komplexen Reaktionskoordinaten auch zu Konvergenzproblemen bzw. zu langwierigen Optimierungszyklen kommen.
Als weiteres Resultat liefern die durchgeführten MEP-Berechnungen Einblicke in die katalytischen Eigenschaften der Wassermoleküle für den Protonübertragungsmechanismus nach Grotthus. Die Daten zeigen, dass die Barriere am niedrigsten wird, wenn zwei Wassermoleküle beim Protontransfer beteiligt sind. Wenn nur ein oder gar kein Wassermolekül die Ringöffnung katalysiert, steigt die Barriere auf 12 und 17 kcal/mol.
Die Untersuchung in diesem Teil der Arbeit lässt zudem Einblicke in die nukleophile Substitution der Vollacetale in der Enzym-Tasche der GM II erlangen. Die Rechnungen deuten darauf hin, dass die Vollacetal-Inhibitoren durch Wassermoleküle in der Tasche aktiviert werden. Die ausgebildeten Wasserstoffbrückenbindungen begünstigen die Geometrie des Enzym-Inhibitor-Komplexes. Dies befördert die Ringöffnungreaktion gleichzeitig mit dem nucleophilen Angriff des Aspartatrestes an dem C1-Atom des Inhibitors. Im Falle des gemischten Acetal-Inhibitors hingegen wird die Treibkraft bereits durch die Einführung des Stickstoffatoms deutlich erhöht. Durch die richtig angeordneten Grotthus-Wassermoleküle ist in diesem Fall die Barriere der Protonübertragung durch das Aspartat-Aspartat-System der GM II (Asp341/Asp240) sekundär. Betrachtet man die Schwingungsbewegung entlang der imaginären Moden der Übergangszustände, sind diese in beiden E-I-Komplexen ähnlich. Hierbei wird eine synergistische Bewegung der Bindungsabstände OD2-C1-O6 (Inhibitor Nr.1) bzw. OD2-C1-N (Inhibitor Nr.7) beobachtet.
Die Entwicklung der kovalent-reversiblen Inhibitoren für das Enzym Rhodesain wurde in Zusammenarbeit mit dem Arbeitskreis Prof.Dr. T. Schirmeister durchgeführt. Als Leitstruktur zur Entwicklung des neuen kovalent-reversiblen Vinylsulfon-Inhibitors 4-Pyridyl-Phenylalanyl-Homophenylalanyl-alpha-Fluor-Phenylvinylsulfon dient in diesem Projekt der kovalent-irreversibel bindende Inhibitor K777, für den kristallographische Daten bekannt sind. Im Rahmen des Protokolls wurde eine Reihe von Inhibitoren untersucht, in denen ein alpha-H-Atom der Vinylsulfon-Einheit (im Weiteren VS für Vinylsulfon) durch verschiedene Gruppen X substituiert wurde. Für den zunächst vorgeschlagenen Cyano-Substituent (CN) in einem VS-Inhibitor ergab sich bei Berechnungen in einem polaren Lösungsmittel eine relativ niedrige Reaktionsenergie, d.h. es wurde eine reversible Reaktion vorhergesagt. Dies wurde experimentell bestätigt. Die theoretischen und experimentellen Ergebnisse von Schritt II widersprechen sich aber. Während experimentell eine schwache reversible Hemmung gefunden wurde, sagten die Berechnungen keine Hemmung voraus. Tatsächlich zeigte sich im Nachhinein, dass die experimentell gefundene Hemmung nicht-kompetitiv ist, d.h. nicht in der aktiven Tasche stattfindet. Im Rahmen des Protokolls werden dagegen nur die kompetitiv interagierenden Inhibitoren ausgewertet. An dieser Stelle lassen sich die anhand theoretischen Methoden erhaltenen Daten über die Reversibilität der Hemmung (Reaktion im aktiven Zentrum) mit den Ergebnissen aus den experimentellen Untersuchungen (Reaktion außerhalb des aktiven Zentrums) nicht vergleichen. Durch den Wechsel von CN zu Halogenen wurde schließlich eine neue Reihe von Inhibitoren auf VS-Basis entwickelt. Die Berechnungen von Reaktionsenergien in Lösung und im Enzym haben für diese Inhibitoren eine reversible Hemmung vorhergesagt. Allerdings findet man eine einfache Additionsreaktion an der Doppelbindung der Vinylsulfon-Gruppe. Für X=CN wurde von einer SN2-Reaktion ausgegangen. Für X=Br fand man, dass sich nach der Addition ein HBr-Molekül abspaltet, sodass die Hemmung insgesamt irreversibel ist. Da die Substitutionsreaktion ein irreversibler Prozess ist und die Freisetzung von Bromwasserstoff durch die experimentellen Untersuchungen bestätigt werden konnte, scheint Fluor ein geeigneter Substituent zu sein (X=F). Hier konnte man auch experimentell eine deutlich bessere Hemmung messen. Es wurden daher die Berechnungen im Enzym für Systeme mit den Inhibitoren K777-X mit X=F und X=H (K777-F- und K777-H-Inhibitor) durchgeführt und analysiert. In der vorliegenden Arbeit wurde versucht, die Reversibilität des K777-F-Inhibitors gegen die Irreversibilität des K777-H-Inhibitors durch die quantenmechanischen Berechnungen im Rahmen des Protokolls darzulegen.
Die QM/MM-Berechnungen unterteilen sich in zwei Bereiche. Zunächst wurde das Reaktionsprofil (auch Reaktionspfad) der Additionsreaktion des K777-X-Inhibitors an die aktive Tasche von Rhodesain ausgehend von der vorhandenen Kristallstruktur (PDB-Datei) berechnet. Im Folgenden wird dieses Teilergebnis als XP-Berechnung (im Weiteren XP für X-ray-Pfad) bezeichnet. Alle vier PES (X=H, F, Br und Cl) weisen prinzipiell die gleiche Form auf. Es ergeben sich aber Unterschiede in den berechneten Reaktionsenergien der Additionsreaktion (-20, -16, -10 und -11 kcal/mol). Die Reaktionsenergien der Substituenten Brom und Chlor entsprechen dem Bereich für reversible Reaktionen (ca. -10 kcal/mol), wobei Fluor mit -16 kcal/mol einen Grenzfall darstellt. Die Konturen der beiden PES (X=H vs. X=F) sind allerdings sehr ähnlich: In beiden Fällen findet sich für das anionische Intermediat kein Minimum. In der Potentialhyperfläche für X=F steigt die Barriere der Rückreaktion zwischen dem Intermediat und dem nicht-kovalenten Komplex auf etwa 5 kcal/mol an, die Rückreaktion ist im Vergleich zu dem X=H mit ca. 1,5 kcal/mol leicht exotherm. Das veränderte Verhältnis zwischen der Höhe der Reaktionsbarriere und dem Betrag der Reaktionsenergie (der Übergang von der endothermen zur exothermen Reaktion) auf diesem Abschnitt der PES könnte dazu beitragen, dass die Gesamtreaktion insgesamt reversibel ablaufen kann. Die Reversibilität des Inhibitors mit dem Substituenten Fluor lässt sich auf diesem Schritt der Untersuchung durch die Absenkung der Reaktionsenergie der Additionsreaktion auf etwa 16 kcal/mol erklären, denn die irreversible Reaktionen wurden bisher mit deutlich höheren Reaktionsenergien assoziiert. Die erhaltenen nicht-kovalenten Enzym-Inhibitor-Komplexe der XP-Berechnungen wurden in einem zweiten Teilergebnis weiter verwendet, indem der Reaktionspfad der Additionsreaktion des K777-X-Inhibitors vom nicht-kovalenten Enzym-Inhibitor-Komplex zum kovalenten Enzym-Inhibitor-Komplex hin berechnet wurde. Im Folgenden wird dieses Teilergebnis als NP-Berechnung (NP für Nicht-kovalente-Pfad) bezeichnet.
Die Berechnung der Reaktionsprofile der Additionsreaktion des VS-Inhibitors für X=H und X=F am alpha-Kohlenstoffatom der VS-Kopfgruppe lieferte konsistente Ergebnisse in Bezug auf die Reaktionsenergien. Ähnlich den XP-Berechnungen, wurde ebenfalls die Tendenz der Absenkung der Reaktionsenergie von -7 kcal/mol (X=H) auf -4,3 (X=F) und -0,9 kcal/mol (X=Br) beobachtet. Die Thermodynamik der Additionsreaktion wurde durch einen Wechsel des Substituenten X von H nach F in der VS-Kopfgruppe des K777-X Inhibitors beeinflusst, indem die niedrigere Energiedifferenz zwischen den Edukten und den Produkten erzielt werden konnte. Für beide Teile der Arbeit (XP- und NP-Berechnungen) implizierte dies einen Wechsel von einem irreversiblen zum einem reversiblen Verlauf in der Beschreibung der Reaktionsprofile. Die Ergebnisse des zweiten Teils der Arbeit (NP-Berechnungen) liefern nicht nur die konsistent geringeren Reaktionsenergien (Thermodynamik) sondern auch die höheren Reaktionsbarrieren der Additionsreaktion im Vergleich zu den Ergebnissen der XP-Berechnungen. Die Änderungen der Reaktionsbarrieren im NP-Ansatz weisen zusätzliche Diskrepanzen auf, wenn diese jeweils mittels der PES-Scan- und CI-NEB-Dimer-Methodik berechnet werden. Während die Barriere des irreversiblen Inhibitors K11777 mit dem NEB-Ansatz ca. 11 kcal/mol beträgt und durch den PES-Scan nur um 4 kcal/mol höher liegt, ergibt sich eine umgekehrte Situation beim Übergang zu Fluor als Substituent: Durch die NEB-Berechnung liegt die Barriere bei ca. 18 kcal/mol und durch den PES-Scan ergibt sich eine um 4 kcal/mol niedrig liegende Barriere. Um die Ergebnisse der NP-Berechnungen zu überprüfen, wurden diese QM/MM-Rechnungen wiederholt durchgeführt. In den beiden neu durchgeführten Berechnungen für die Inhibitoren K777-X mit X=H und X=F wurden nur sehr kleine Abweichungen gefunden, die kleiner als die Fehler der Berechnung sind. Die Startstrukturen für die Berechnung des MEP stammten aus der erneut durchgeführten MD-Simulation an der Geometrie des nicht-kovalenten Enzym-Inhibitor-Komplexes, welche die XP-Berechnung resultierte. Die Gesamtdauer der MD-Simulation wurde zu einem Wert von 9 ns gewählt, welche insgesamt 900 Startstrukturen entlang der Simulation lieferte. Die Berechnung ergab die Reaktionsenergie von -8,4 kcal/mol (-7,0 kcal/mol als erstes Ergebnis) und die relative Energie des Int-Komplexes von 13,2 kcal/mol. Somit beträgt die Barriere der Rückreaktion zur Freisetzung des Inhibitors K11777 (X=H) in Form von einem nicht-kovalenten Enzym-Inhibitor-Komplex einen Wert von 21,6 kcal/mol. In analoger Vorgehensweise wurde die Evaluierung der NP-Berechnung für den Inhibitor K777-X mit X=F durchgeführt. Die Reaktionsenergien in den beiden Berechnungen unterscheiden sich in einem marginalen Abstieg zu den Werten von -2,9 kcal/mol und -3,3 kcal/mol (-4,3 kcal/mol als erstes Ergebnis). Beide Berechnungen liefern zudem die relativ kleinen Anstiege der Reaktionsbarriere zu den Werten von 19,8 kcal/mol und 20,9 kcal/mol.
Für die Inhibitoren K777-X mit X=H und X=F entsprechen die gefundenen Barrieren einer verzerrten Konformation des nicht-kovalenten Enzym-Inhibitor-Komplexes, die als eine bioaktive Konformation bezeichnet werden kann. Der anionische Übergangszustand Int*, der oft in der Literatur als ein anionisches Intermediat der Additionsreaktion bezeichnet wird, wurde nur für den Inhibitor mit dem Substituenten Brom (K777-X mit X=Br) identifiziert. Da der Übergangszustand (Int* mit der relativen Energie von 11,1 kcal/mol) nur 1,5 kcal/mol über der bioaktiven Konformation (Int mit der relativen Energie von 12,6 kcal/mol) liegt und die NEB-Reaktionspfade alleine die Barrieren überschätzen, besitzen die anionischen Übergangszustände der Additionsreaktion der Inhibitoren mit X=H und X=F eine geschätzte relative Energie mit vergleichbaren Abweichungen von ca. 2 kcal/mol zu den identifizierten Int-Geometrien. Die durchgeführten Berechnungen zeigen, dass die Substituenten X=H und X=F im Vergleich zum Brom die anionischen Geometrien der nicht-kovalenten Enzym-Inhibitor-Komplexe jedoch mangelhaft bis ausreichend stabilisieren können. Zusätzlich liegt die geschätzte Energiedifferenz zwischen den Geometrien Int* und Int unter der möglichen Fehlergrenze der Berechnungen (ca. 3-4 kcal/mol). Aus diesem Grund misslang die Optimierung in Richtung der metastabilen anionischen Geometrien Int* mit Hilfe der CI-NEB-Dimer-Methodik im Fall der VS-Inhibitoren K777-X mit X=H und X=F.
Der direkte Vergleich der geometrischen Parameter der nicht-kovalenten Enzym-Inhibitor-Komplexe für den Inhibitor K777-F aus den XP-Berechnungen mit solchen aus den NP-Berechnungen lässt darauf schließen, dass die Geometrien der Enzym-Inhibitor-Komplexe der XP-Berechnung nur die lokalen Minima mit der verzerrten Geometrie des Inhibitors auf der PES darstellen und die Gesamtinformation über die Barrieren der Reaktion durch die Ergebnisse aus der NP-Berechnung ergänzt werden sollten.
Zusammenfassend sagen die Berechnungen für die reaktiven Kopfgruppen der Substanzklasse der halogenierten Vinylsulfone K777-X (X=Br, Cl und F) im Vergleich zur Leitstruktur des Vinylsulfon-Inhibitors K11777 deutlich geringere exotherme Additionsreaktionen im aktiven Zentrum von Rhodesain voraus. Darüberhinaus konnte anhand der QM/MM-Berechnungen ein experimentell gemessenen verlangsamten Verlauf der reversiblen Inhibition im Falle von X=F (Inhibitor K777-X) durch die relativ erhöhte Reaktionsbarriere im Vergleich zur Leitstruktur erklärt werden. Dieser Inhibitor dient zunächst als ein erfolgreich selektiertes reaktives Gerüst des neuen Inhibitors K777-X-S3 mit X=F und S3=4-Pyridyl (K777-F-Pyr), welcher mit Hilfe des Docking-Experiments (Schritt III durch die Arbeitsgruppe Prof.Dr. T. Schirmeister) deutlich verbessert werden konnte. Die Affinität des durch Docking verbesserten VS-Inhibitors mit Fluor als Substituent durch die eingeführte Seitenkette S3=4-Pyridyl (4-Pyridyl-Phenylalanyl-Homophenylalanyl-(Phenyl)-alpha-F-Vinylsulfon) stieg im Rhodesain von 190 nM zu 32 nM (Schritt IV, experimenteller Teil). Gleichzeitig konnte durch die QM/MM-Berechnungen in Schritt IV gezeigt werden, dass die Reaktion der reaktiven Kopfgruppe im neuen Inhibitor immer noch eine kovalent-reversible Hemmung von Rhodesain darstellt, auch wenn die Erkennungseinheit geändert wurde. Hierfür kann man die beiden Reaktionsprofile der NP-Berechnungen vergleichen. Die beiden fluorierten VS-Inhibitoren weisen eine Ähnlichkeit bezüglich der Barrierenhöhe und der Reaktionsenergie auf. Der fluorierte Vinylsulfon-Inhibitor K777-F wurde somit als ein neuer kovalent-reversibler Vinylsulfon-Inhibitor der Cysteinprotease Rhodesain erfolgreich eingefügt.
Within this work, an additive and a subtractive QM/MM interface were implemented into CAST. The interactions between QM and MM system are described via electrostatic embedding. Link atoms are used to saturate dangling bonds originating from the separation of QM and MM system. Available energy evaluation methods to be combined include force fields (OPLSAA and AMBER), semi-empirical programs (Mopac and DFTB+), and quantum-chemical methods (from Gaussian, Orca, and Psi4). Both the additive and the subtractive interface can deal with periodic boundary conditions. The subtractive scheme was extended to enable QM/QM, three-layer, and multi-center calculations. Another feature only available within the subtractive interface is the microiteration procedure for local optimizations.
The novel QM/MM methods were applied to the investigation of the reaction path for the complex formation between rhodesain and K11777. Benchmark calculations show a very good agreement with results from Gaussian-ONIOM. When comparing the relative energies obtained with different options to a computation where the whole system was treated with the “QM method” DFTB3, the electrostatic embedding scheme with option “delM3” gives the best results. “delM3” means that atoms with up to three bonds distance to the QM region are ignored when creating the external charges. This is done in order to avoid a double counting of Coulomb interactions between QM and MM system. The embedding scheme for the inner system in a three-layer calculation, however, does not have a significant influence on the energies. The same is true for the choice of the coupling scheme: Whether the additive or the subtractive QM/MM interface is applied does not alter the results significantly. The choice of the QM region, though, proved to be an important factor. As can be seen from the comparison of two QM systems of different size, bigger is not always better here. Instead, one has to make sure not to separate important (polar) interactions by the QM/MM border.
After this benchmark study with singlepoint calculations, the various possibilities of CAST were used to approximate the solution of a remaining problem: The predicted reaction energy for the formation of the rhodesain-K11777 complex differs significantly depending on the starting point of the reaction path.
The reason for this is assumed to be an inadequate adjustment of the environment during the scans, which leads to a better stabilization of the starting structure in comparison to the final structure. The first approach to improve this adjustment was performing the relaxed scan with a bigger QM region instead of the minimal QM system used before. While the paths starting from the covalent complex do not change significantly, those starting from the non-covalent complex become more exothermic, leading to a higher similarity of the two paths. Nevertheless, the difference of the reaction energy is still around 15 kcal/mol, which is far from a perfect agreement. For this reason, Umbrella Samplings were run. Here, the adjustment of the environment is not done by local optimizations like in the scans, but by MD simulations. This has the advantage that the system can cross barriers and reach different local minima. The relative free energies obtained by Umbrella Samplings with suitable QM regions are nearly identical, independently of the starting point of the calculation. Thus, \(\Delta A\) evaluated by these computations can be assumed to reproduce the real energy change best. An MD simulation that was started from the transition state in order to mimic a “real-time” reaction indicates a very fast adjustment of the environment during the formation of the complex. This confirms that Umbrella Sampling is probably better suitable to describe the reaction path than a scan, where the environment can never move strong enough to leave the current local minimum.