Refine
Has Fulltext
- yes (6)
Is part of the Bibliography
- yes (6)
Document Type
- Doctoral Thesis (6)
Keywords
- Chemische Bindung (6) (remove)
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.
Die vorliegende Arbeit beschäftigt sich mit der Anwendung des Sulfonyl-stabilisierten Methandiids 20-Li2 als Ligand in Übergangsmetallkomplexen. Dabei konnte 20-Li2 mit ver-schiedenen Übergangsmetallhalogeniden in Salzmetathesereaktionen umgesetzt werden. Insgesamt wiesen die synthetisierten Methandiid-basierten Komplexe flexible Bindungsverhältnisse bezüglich der MC-Bindung und unterschiedliche Koordinationsmodi der Sulfonyl-Gruppe auf, die die Stabilität und Reaktivität der Komplexe signifikant beeinflussten. In Abhängigkeit von der chemischen Natur des Metallfragmentes und der Co Liganden konnten Carbenkomplexe mit einer ylidischen M-C-Wechselwirkung (A) und solche mit einer echten M=C-Doppelbindung (B) zugänglich gemacht werden. Dabei gelang die Etablierung einer Vielzahl an neuen Komplexen sowohl mit frühen (Zirkonium) als auch späten (Palladium, Ruthenium, Iridium) Übergangsmetallen. Die synthetisierten Verbindungen zeigten dabei unterschiedliche strukturelle und elektronische Eigenschaften, was zu deren Unterteilung in die zwei Komplexklassen A und B führte. So konnte bei der Umsetzung von Methandiid 20-Li2 mit Zirkonocendichlorid die selektive Bildung des Zirkonocenkomplexes 50 beobachtet werden, bei dem NMR spektroskopische (z.B. Hochfeldverschiebung des 13C NMR-Signals des Carben-Kohlenstoffatoms) und röntgenstrukturanalytische (z.B. Pyramidalisierung des „Carben“-Kohlenstoffatoms) Untersuchungen erste Hinweise darauf lieferten, dass sich 50 nicht als Carbenkomplex mit einer Zr=C-Doppelbindung beschreiben lässt. Dies konnte durch quantenchemische Rechnungen bestätigt werden, wobei die „Natural Bond Orbital“-Analyse (NBO-Analyse) eine deutliche negative Ladung am zentralen Kohlenstoffatom (qc = 1.42) und somit dessen nukleophilen Charakter aufdeckte. Zusätzlich lieferten die Rechnungen eine deutlich positive Ladung am Zirkoniumatom (qZr = 1.35), weshalb die Zr-C-Interaktion in 50 am besten mit einer ylidischen Wechselwirkung beschrieben wird. Ähnliche Resultate konnten auch bei den aus den Umsetzungen von 20-Li2 mit [(PPh3)2PdCl2] bzw. [(PPh3)3RuCl2] erhaltenen Komplexen 51a bzw. 52-Int beobachtet werden. Wie für Verbindung 50 ergab die NBO-Analyse von 51a bzw. 52-Int zwar eine Bindung zwischen Metall- und Kohlenstoffatom, interessanterweise aber keine Wechselwirkung. Aufbauend auf der elektronischen Struktur von 51a bzw. 52-Int zeichnen sich die beiden Komplexe durch eine hohe Instabilität und Reaktivität aus. Dabei bildete 51a in Lösung diverse Zersetzungsprodukte, während der Ruthenium-Carbenkomplex 52-Int selektiv die Phenylgruppe des Sulfonyl-Substituentens in ortho Position unter Ausbildung der cyclometallierten Spezies 52 intramolekular deprotonierte. Das Cyclometallierungsprodukt 52 konnte in einer Ausbeute von 62% isoliert und vollständig charakterisiert werden.
Die schwache -Interaktion zwischen Metall- und Kohlenstoffatom konnte im Falle der Palladium- und Rutheniumkomplexe auf den Elektronenreichtum der späten Übergangsmetalle zurückgeführt werden, welcher durch die guten Donor- und schlechten Akzeptoreigenschaften der Phosphan-Liganden zusätzlich verstärkt wurde. Durch Austausch der Triphenylphosphan-Liganden in der Rutheniumdichlorid-Vorstufe gegen das Aren p-Cymol konnte die elektronische Natur am Metallfragment derartig beeinflusst werden, dass ein selektiver Zugang zu Ruthenium-Carbenkomplex 53 gelang. Verbindung 53 konnte in einer guten Ausbeute von 86% in Form eines dunkelvioletten Feststoffes isoliert und vollständig charakterisiert werden. Dass es sich bei 53 tatsächlich um einen Carbenkomplex mit einer M=C-Doppelbindung handelt, konnte mithilfe der Molekülstruktur im Festkörper, den NMR-spektroskopischen Daten und der berechneten elektronischen Struktur bestätigt werden. So wies 53 eine kurze Ru=C-Bindung und eine planare Koordinationsumgebung des zentralen Kohlenstoffatoms [Winkelsumme: 358.9(1) Å] auf. Zusätzlich sprachen die im Vergleich zu Methandiid 20-Li2 verlängerten P-C- und C-S-Abstände für geschwächte elektrostatische Wechselwirkungen im Ligand-Rückgrat und somit für einen effizienten Elektronentransfer vom Methandiid zum Metall. Die NBO-Analyse ergab sowohl eine - als auch -Wechselwirkung der M-C-Bindung mit einer nur leichten Polarisierung zum Kohlenstoffatom. Ähnliche Beobachtungen (kurzer Ir-C-Abstand, Planarität am Kohlenstoffatom, reduzierte elektrostatische Wechselwirkungen im Ligand-Rückgrat, NBO-Analyse) wurden ebenfalls für den Iridium-Carbenkomplex 53 gemacht.Die negativere Ladung am Carben-Kohlenstoffatom wies hierbei allerdings auf einen leicht ylidischeren Charakter der MC-Bindung als im Ruthenium-Analogon 53 hin.
Aufbauend auf der elektronischen Natur der M=C-Bindung ergaben sich unterschiedliche Reaktivitäten der Carbenkomplexe. Während der Zirkonocenkomplex 50 gegenüber Aldehyden, Ketonen und Disulfiden entweder keine Reaktivität oder Zersetzung zum zweifach protonierten Liganden zeigte, erfolgte ausgehend von Ruthenium-Carbenkomplex 52-Int die intramolekulare CH-Aktivierung zu 52. Im Gegensatz dazu konnte der Ruthenium-Carbenkomplex 53 in einer Vielzahl von EH-Bindungsaktivierungen eingesetzt werden.
Dabei konnten zahlreiche E-H-Bindungen bei Raumtemperatur aktiviert und das nicht-unschuldige Verhalten des Methandiid-Liganden unter Beweis gestellt werden. So konnten die O-H- und N-H-Bindungen in einer Serie von Alkoholen und Aminen, die P-H-Bindung in sekundären Phosphanoxiden und die hydridischen SiH- und BH-Bindungen in Silanen und Boranen durch 53 gespalten werden. Durch röntgenstrukturanalytische Aufklärung der Molekülstrukturen im Festkörper konnte gezeigt werden, dass die Bindungsaktivierung im Allgemeinen unter 1,2-Addition der Substrate auf die Ru=C-Doppelbindung unter Bildung der entsprechenden cis-Additionsprodukte erfolgte. Die Aufhebung der Metall-Kohlenstoffdoppel- zu einer -einfachbindung machte sich in einer Verlängerung der Ru=C-Bindung von 1.965(2) Å in 53 auf etwa 2.2 Å bemerkbar. Zudem konnte in allen Molekülstrukturen der Aktivierungsprodukte eine Pyramidalisierung des ehemals planaren Carben-Kohlenstoffatoms detektiert werden.
Bezüglich der Regioselektivität verliefen die Umsetzungen mit Substraten, in denen das Wasserstoffatom einen protischen (O-H, N-H-Bindungen) bzw. mäßig protischen/hydridischen (P-H-Bindungen) Charakter aufweist, erwartungsgemäß unter Protonierung des nukleophilen Carben-Kohlenstoffatoms. Interessanterweise führten die O-H- und N-H-Aktivierungsreaktionen z.T. zur Ausbildung eines Gleichgewichts zwischen Carbenkomplex und Additionsprodukt. Dabei konnte ein derartiger Gleichgewichtsprozess in der Chemie Methandiid-basierter Carbenkomplexe bisher nicht beobachtet werden, was die außerordentliche Stabilität des Rutheniumkomplexes 53 unterstreicht. Diese Reversibilität wurde bspw. anhand der Umsetzung von Komplex 53 mit p Methoxyphenol mittels VT-NMR-Studien untersucht. Dabei konnte gezeigt werden, dass sich das Gleichgewicht beim Abkühlen auf 80 °C gemäß entropischer Effekte fast vollständig auf die Seite des Additionsproduktes verschieben lässt, während beim Erwärmen auf Raumtemperatur das Gleichgewicht auf der Seite des Carbenkomplexes liegt. Ähnliche Gleichgewichtsprozesse konnten bei der N-H-Aktivierung beobachtet werden. Bei der Aktivierung von Ammoniak konnte das Additionsprodukt 60 nicht isoliert werden, da auch hier ein stark temperaturabhängiges Gleichgewicht vorlag, wobei erst ab 90 °C das Gleichgewicht vollständig auf der Seite des Aktivierungskomplexes 60 lag. Daher konnte 60 nicht isoliert und eindeutig identifiziert werden. In folgenden Arbeiten sollte die Isolierung von 60 im Festkörper angestrebt und somit dessen Existenz nachgewiesen werden. Zudem könnten auch hier Übertragungsreaktionen des aktivierten Ammoniaks auf ungesättigte Substrate durchgeführt werden.
Überraschenderweise zeigte die Si-H-Bindungsaktivierung von unterschiedlich substituierten aliphatischen und aromatischen Silanen ein analoges Reaktionsmuster und führte zur selektiven Bildung der entsprechenden Silylkomplexe 66a-66f anstelle der aufgrund der Polaritäten zu erwartenden Hydrido-Spezies. Mittels DFT-Rechnungen konnte gezeigt werden, dass der Reaktionsmechanismus der SiH-Aktivierung nicht über eine konzertierte 1,2-Addition, sondern über einen zweistufigen Prozess verläuft. Dabei ermöglichen die flexiblen Koordinationseigenschaften des Liganden in 53 eine oxidative Addition der Si-H-Bindung an das Ruthenium-Zentrum, auf die ein Hydrid-Transfer zum Methandiid-Kohlenstoffatom folgt. Neben einfachen Bindungsaktivierungen wurde das Potential der synthetisierten Silylkomplexe 66a-66c in Hydrosilylierungsreaktionen untersucht. In diesem Zusammenhang wurde die Hydrosilylierung von Norbornen angestrebt. Während bei Raumtemperatur keine Reaktion stattfand, konnte nach Erhöhung der Temperatur auf 80 °C das gewünschte Hydrosilylierungsprodukt 68 zwar mittels GC-MS-Analytik nachgewiesen werden, jedoch entstand bei den gewählten Reaktionsbedingungen das ROMP-Produkt 69 als Hauptprodukt. In weiterführenden Arbeiten müssen noch Optimierungsversuche der Reaktionsbedingungen zu einem selektiveren Umsatz zum Hydrosilylierungsprodukt durchgeführt werden.
Interessante Ergebnisse lieferten zudem die Umsetzungen des Ruthenium-Carbenkomplexes 53 mit Boranen und verschiedenen Boran-Lewis-Basen-Addukten. Dabei führte die Reaktion von 53 mit Catecholboran zur Bildung des Hydridokomplexes 73, dessen Molekülstruktur im Festkörper bestimmt werden konnte. Jedoch konnte der Komplex aufgrund seiner Instabilität in Lösung bisher nicht vollständig NMR-spektroskopisch und mittels Elementaranalytik charakterisiert werden. Im Gegensatz zur Si-H-Aktivierung findet hier die Addition entsprechend der Polarität der B-H-Bindung statt. Erstaunlicherweise führte die BH Bindungsaktivierung in Pinakolboran jedoch nicht zu einer zu 73 analogen Hydrid-Spezies. Der NMR-spektroskopische Verlauf der Umsetzung deutete zunächst auf die Bildung des BH-Additionsproduktes unter Protonierung des PCS-Rückgrats hin, welches sich in eine andere, bisher nicht identifizierbare Spezies umwandelte. Wiederum zu einem anderen Ergebnis führte die Umsetzung von 53 mit BH3∙SMe2. Durch Insertion eines Borans in die Thiophosphoryl-Einheit unter Aktivierung der B-H-Bindung wurde hierbei Komplex 76 gebildet, der als zentrales Strukturmotiv einen P–B–S–Ru–C-Fünfring aufwies.
Neben der Spaltung polarer E-H-Bindungen gelang außerdem die Aktivierung der unpolaren Bindung in Diwasserstoff unter Bildung des Hydridokomplexes 77. Mittels Röntgenstrukturanalyse konnte auch hier eine cis-Addition von H2 auf die RuC-Doppel-bindung bestätigt und das Signal des hydridischen Wasserstoffatom eindeutig im 1H NMR-Spektrum der Verbindung bei H = 6.62 ppm detektiert werden. Interessanterweise konnte Verbindung 77 ebenfalls durch Dehydrierung von iso Propanol bzw. Ameisensäure (HCOOH) unter Abspaltung von Aceton bzw. CO2 synthetisiert werden.
Aufbauend auf der beobachteten Dehydrierung von iPrOH unter Bildung des Hydridokomplexes 77 wurde der Frage nach einer möglichen Anwendung des Carbenkomplexes 53 in der katalytischen Transferhydrierung von Ketonen zu Alkoholen nachgegangen. Obgleich die Aktivierung von H2 bzw. die Dehydrierung von iPrOH keine Reversibilität aufwies, sollte ein Katalysezyklus basierend auf einem Wechselspiel zwischen Carben- 53 und Hydridokomplex 77 mit iPrOH als Wasserstoffquelle realisierbar sein. Diesbezüglich lieferten erste Reduktionsversuche von Acetophenon zu 1 Phenylethanol mit 53 und KOtBu als Hilfsbase allerdings schlechte Alkohol-Ausbeuten im Vergleich zu literaturbekannten, übergangsmetallkatalysierten Transferhydrierungen. Ein Katalyseansatz mit 0.50 mol-% 53 und 19 mol-% KOtBu ergab nach 24 h bei 75 °C eine Alkohol-Ausbeute von gerade einmal 55%. Zudem konnte eine starke Abhängigkeit der Umsätze von der eingesetzten Basenmenge beobachtet werden, was auf eine konkurrierende, Basen-induzierte Reduktion hindeutete. Eine Optimierung der Katalysebedingungen gelang durch Zugabe von Triphenylphosphan. Mithilfe des Additivs konnte innerhalb von 12 h bei 75 °C mit 0.50 mol-% 53, 6.20 mol-% KOtBu und 6.20 mol-% PPh3 ein nahezu quantitativer Umsatz (94%) von Acetophenon zu 1-Phenylethanol beobachtet werden. Sogar eine Verringerung der Basen- und Phosphanmenge auf 1.60 und 1.10 mol-% reichte aus, um Ausbeuten von 90% zu erreichen (Abb. 4.5., rechts). Dabei konnte Rutheniumkomplex 53 als erster Methandiid-basierter Carbenkomplex mit katalytischem Potential in Transferhydrierungen etabliert werden. Außerdem beschränkte sich die katalytische Aktivität von 53/PPh3 nicht nur auf die Reduktion von Acetophenon, sondern konnte auch erfolgreich auf weitere aromatische und aliphatische Ketone übertragen werden.
Mittels NMR-spektroskopischer Untersuchungen des Katalyseverlaufs gelang ein Nachweis der katalytisch aktiven Spezies im Katalysezyklus. So konnte bei 75 °C zunächst die erwartungsgemäße Entstehung des Hydridokomplexes 77 beobachtet werden. Dieser setzte sich anschließend mit PPh3 zum cyclometallierten Phosphankomplex 52 um. Aufbauend auf diesen Beobachtungen wurde ebenfalls Komplex 52 hinsichtlich seines katalytischen Potentials in der Reduktion von Acetophenon untersucht, wobei noch bessere Umsätze als mit dem Katalysator 53/PPh3 beobachtet wurden. Hierbei konnte bereits nach 3 h mit 0.50 mol-% 52 und 1.60 mol-% KOtBu eine Ausbeute von 95% erzielt werden. Zudem führten Ansätze mit 52 auch ohne Zugabe einer Base zu Umsätzen von ca. 40%. Eine Übertragung der Katalysebedingungen auf die Reduktion weiterer Keton-Derivate lieferte ebenfalls gute Ergebnisse und ergab Alkohol-Ausbeuten zwischen 72% und 96%.
Die für Ruthenium-Carbenkomplex 53 gefundene Reaktivität und das nicht-unschuldige Verhalten des Methandiid-Liganden konnten außerdem auch für Iridium-Carbenkomplex 55 beobachtet werden. So konnten analoge NH, PH- und SiH-Additionsprodukte selektiv synthetisiert und in guten Ausbeuten (etwa 60-90%) analysenrein erhalten werden.
In Analogie zu Rutheniumkomplex 53 führte die Aktivierung von Substraten mit unterschiedlichen E-H-Bindungen entsprechend der Ladungsverteilung im Ir+C--Fragment zur Protonierung der PCS-Brücke in 55. Dabei wiesen auch hier die Additionsprodukte im Allgemeinen eine cis-Anordnung der vorherigen E-H-Einheit auf. Einzige Ausnahme stellte das mit p-Nitroanilin gebildete NH-Aktivierungsprodukt 61b dar. Hierbei konnte mittels Röntgenstrukturanalyse eine trans-Anordnung der Amido-Einheit und des PCHS-Brückenprotons detektiert werden, die durch Ausbildung einer Wasserstoffbrückenbindung zwischen der Amido-NH- Einheit und dem Sauerstoffatom des Sulfonyl-Substituentens begünstigt wird. Zudem konnte für die Bildung von 61b ein bei Raumtemperatur reversibler Reaktionsprozess unter Rückbildung des Carbenkomplexes 55 und Abspaltung von p-Nitroanilin beobachtet werden. In künftigen Experimenten sollte untersucht werden, ob aufgrund der Reversibilität katalytische Hydroaminierungen mit 61b realisierbar sind.
Trotz des hydridischen Charakters des Si-H-Wasserstoffatoms in Silanen wurden auch mit Carbenkomplex 55 ausschließlich die SiH-Bindungsaktivierungskomplexe 71a-71c gebildet. Zudem konnte bei der Aktivierung von Triphenylsilan zwar das Additionsprodukt 71a mittels NMR-spektroskopischer Untersuchungen in der Reaktionslösung nachgewiesen werden, jedoch setzte sich dieses bereits bei Raumtemperatur zum cyclometallierten Komplex 72 um. Interessanterweise resultierten die Aktivierung von H2 und die Dehydrierung von iPrOH ebenfalls in 72. Mittels NMR-spektroskopischer Untersuchungen des Reaktionsverlaufes konnte hierbei gezeigt werden, dass die Cyclometallierung ausgehend von dem in situ gebildeten Iridium-Hydridokomplex 79 stattfindet. Deuterierungsexperimente mit iPrOH-d8 belegten außerdem, dass die Protonierung der PCS Brücke durch iPrOH und nicht durch direkte ortho-C-H-Aktivierung der Sulfonyl-Phenyl-Gruppe erfolgt. Die Isolierung des Iridium-Hydridokomplexes 79 war aufgrund seiner schnellen Umsetzung zu 72 daher nicht möglich.
Die Nukleophilie des Carben-Kohlenstoffatoms und die ausgezeichnete M=C-Wechselwirkung in Ruthenium-Carbenkomplex 53 ermöglichten neben EH-Bindungsaktivierungen außerdem [2+2]-Cycloadditionsreaktionen mit Iso- und Thioisocyanaten. In diesem Zusammenhang konnten mit tert-Butyl- und Phenylisocyanat die Cycloadditionsprodukte 80a und 80b synthetisiert, in guten Ausbeuten isoliert (79% bzw. 80%) und vollständig charakterisiert werden. Die mittels Röntgenstrukturanalyse durchgeführte Aufklärung der Molekülstruktur von 80a im Festkörper bestätigte die Ausbildung eines C-Ru-N-C-Vierringes als zentrales Strukturmotiv, was mit literaturbekannten Umsetzungen dieser Art übereinstimmt. Mit tert-Butyl- und Phenylthioisocyanat hingegen wurden die Iminkomplexe 81a/b unter Addition der Ruthenium-Kohlenstoff-Doppelbindung an das CS-Fragment im Thioisocyanat erhalten. Dabei konnte die Ausbildung eines C-Ru-S-C-Vierringes als zentrales Strukturmotiv beobachtet werden.
Insgesamt folgte die Selektivität der gebildeten [2+2]-Cycloadditionsprodukte 80a/b und 81a/b den Prinzipien des HSAB-Konzeptes, wonach jeweils das weichere Atom des Heteroallens an das Ruthenium-Zentrum bindet. Obgleich die Reaktivität Methandiid-basierter Carbenkomplexe mit verschiedenen Heteroallenen bereits in der Literatur beschrieben wurde, stellte die Umsetzung mit Thioisocyanaten zu 81a/b ein bisher unbekanntes Reaktionsverhalten dieser Verbindungsklasse dar.
Neben der Anwendung des Methandiids 20-Li2 als Ligand für die Synthese neuer Übergangsmetallkomplexe erwies sich das Dianion außerdem als geeignet für die Darstellung des Li/Cl-Carbenoids 83. Dabei konnte 83 zum einen durch Oxidation von 20-Li2 mit Hexachlorethan (C2Cl6) und zum anderen durch Metallierung des chlorierten Liganden 82 synthetisiert und in guten Ausbeuten (67-82%) als farbloser, kristalliner Feststoff isoliert werden. Verbindung 83 erwies sich dabei als ein seltenes, bei Raumtemperatur stabiles Li/Cl-Carbenoid. Aufgrund der Stabilität im Festkörper als auch in Lösung bei Raumtemperatur konnte 83 zudem NMR-spektroskopisch und mittels Elementaranalytik vollständig charakterisiert werden. Ebenfalls gelang die Aufklärung der Molekülstruktur von 83 im Festkörper. Diese zeigte keinen direkten Kontakt zwischen dem Carbenoid-Kohlenstoff- und Lithiumatom und lieferte damit neben der elektronischen Stabilisierung eine Erklärung für die beobachtete Stabilität von 83. Dabei beteiligt sich das Ligandsystem durch Koordination der Sulfonyl-Gruppen an das Lithiumatom erheblich an der Stabilisierung, sodass eine Lithiumchlorid-Eliminierung erschwert wird. Außerdem zeigte die Molekülstruktur keine Verlängerung der C-Cl-Bindung, wie es für unstabilisierte Carbenoide in der Literatur beschrieben wird. Diese Tatsache und die im 13C-NMR-Spektrum beobachtete Abschirmung des Carbenoid-Kohlenstoffatoms im Vergleich zur chlorierten Vorstufe 82 lieferten erste Anzeichen für einen geringen carbenoiden Charakter von 83. Außerdem bestätigten quantenchemische Rechnungen keine signifikante Polarisierung der CCl-Bindung.
Die durch die Stabilisierung resultierende Verringerung des carbenoiden Charakters und somit der Ambiphilie spiegelte sich auch in der Reaktivität von 83 wider. So konnte Verbindung 83 nicht als Cyclopropanierungsreagenz verwendet werden, wie es zumeist für klassische Carbenoide der Fall ist. Gegenüber Elektrophilen wie Methyliodid oder Chlordiphenylphosphan reagierte 83 in Analogie zu Organolithiumbasen zu den Verbindungen 84a und 84b.
Jedoch konnte 83 als Carbenvorstufe zur Synthese des Palladium-Carbenkomplexes 51a unter LiCl-Eliminierung eingesetzt werden, was den leicht vorhandenen carbenoiden Charakter von 83 wiedergibt. Zudem wurde 83 hinsichtlich seines Aktivierungspotentials von EE-Bindungen untersucht. Während die Aktivierung der BH-Bindung in Boranen und die BB-Bindung in Diboranen nicht gelang, konnte die SS-Bindung in 2,2‘-Dipyridyl- und 4,4‘-Dipyridyldisulfid gespalten und Verbindung 90 analysenrein erhalten werden (Schema 4.8.). Studien zur Aufklärung dieses Reaktionsverhaltens stehen jedoch noch aus.
Bezüglich der Aktivierung von P-H-Bindungen in unterschiedlich substituierten aromatischen Phosphanen konnte für 83 eine zu einem Silyl-stabilisierten Carbenoid analoge Reaktivität gefunden werden. Hierbei erfolgte keine Addition der P-H-Bindung an das carbenoide Kohlenstoffatom, sondern die selektive Dehydrokupplung der Phosphane zu Diphosphanen unter LiCl-Eliminierung. Diese überraschende und bis dato für Carbenoide unbekannte Reaktivität erfolgte unter milden Reaktionsbedingungen (Raumtemperatur) und ohne Einsatz von Übergangsmetallkatalysatoren.
Insgesamt konnte für Verbindung 83 ein vielfältiges Reaktionsverhalten gefunden werden. Neben dessen Eignung als Carbenvorstufe bei der Synthese von Übergangsmetall-Carbenkomplexen, konnte die Spezies in der Aktivierung von SS- und PH-Bindungen eingesetzt werden. In zukünftigen Reaktivitätsstudien sollte das beobachtete Potential auf weitere Substrate übertragen werden.
Die vorliegende Arbeit befasst sich mit Untersuchungen von Element-Element- Bindungsaktivierungsreaktionen des dinuklearen Nickel(0)-NHC-Komplex [Ni2(iPr2Im)4(COD)] A mit verschiedenen reaktionsträgen Substraten, die ihrerseits wichtige Ausgangsstoffe für katalytische Anwendungen sind. Die Arbeit gliedert sich dabei in vier verschiedene Teile.
Summary The nature of the chemical bond is a topic under constant debate. What is known about individual molecular properties and functional groups is often taught and rationalized by explaining Lewis structures, which, in turn, make extensive use of the valence concept. The valence concept distinguishes between electrons, which do not participate in chemical interactions (core electrons) and those, which do (single, double, triple bonds, lone-pair electrons, etc.). Additionally, individual electrons are assigned to atomic centers. The valence concept is of paramount success: It allows the successful planning of chemical syntheses and analyses, it explains the behavior of individual functional groups, and, moreover, it provides the “language” to think of and talk about molecular structure and chemical interactions. The resounding success of the valence concept may be misleading to forget its approximative character. On the other hand, quantum mechanics provide in principle a quantitative description of all chemical phenomena, but there is no discrimination between electrons in quantum mechanics. From the quantum mechanical point of view there are only indistinguishable electrons in the field of the nuclei, i.e., it is impossible to assign a given electron to a particular center or to ascribe a particular purpose to individual electrons. The concept of indistinguishability of micro particles is founded on the Heisenberg uncertainty relation, which states, that wavepackets diverge in the 6N dimensional phase space, such that individual trajectories can not be identified. Hence it is a deep-rooted and approved physical concept. As an introduction to the present work density partitioning schemes were discussed, which divide the total molecular density into chemically meaningful areas. These partitioning schemes are intimately related to either the concepts of bound atoms in a molecule (as in the Atoms In Molecules theory (AIM) according to Bader or as in the Hirshfeld partitioning scheme) or to the concept of chemical structure in the sense of Lewis structures, which divide the total molecular density into core and valence density, where the valence density is split up again into bonding and non-bonding electron densities. Examples are early and recent loge theories, the topological analysis by means of the Electron Localization Function (ELF), and the Natural Bond Orbital (NBO) approach. Of these partitioning schemes, the theories according to Bader (AIM), to Becke and Edgecomb (ELF) and according to Weinhold (NBO and Natural Resonance Theory, NRT), respectively, were reviewed in detail critically. Points of criticism were explicated for each of the mentioned theories. Since theoretically derived electron densities are to be compared to experimentally derived densities, a brief introduction into the theory of X-ray di®raction experiments was given and the multipole formalism was introduced. The procedure of density refinement was briefly discussed. Various suggestions for improvements were developed: One strategy would be the employment of model parameters, which are to a maximum degree mutually orthogonal, with the object of minimizing correlations among the model parameters, e.g., to introduce nodal planes into the radial functions of the multipole model. A further suggestion involves the guidance of the iterative refinement procedure by an extremum principle, which states, that when di®erent solutions to the least squares minimization problem are available with about the same statistical measures of quality and with about the same residual density, then the solution is to prefer, which yields a minimum density at the bond critical point (BCP) and a maximum polarity in terms of the ratio of distances between the BCP and the nuclei. This suggestion is based on the well known fact, that the bond polarity (in terms of the ratio of distances between the BCP and the respective nuclei) is underestimated in the experiment. Another suggestion for including physical constraints is the explicit consideration of the virial theorem, e.g., by evaluating the integration of the Laplacian over the entire atomic basins and comparing this value to zero and to the value obtained from the integration of the electron gradient field over the atomic surface. The next suggestion was to explicitly use the electrostatic theorem of Feynman (often also denoted as Hellmann-Feynman theorem), which states, that the forces onto the nuclei can be calculated from the purely classical electrostatic forces of the electron distribution and the nuclei distribution. For a stationary system, these forces must add to zero. This also provides an internal quality criterion of the density model. This can be performed in an iterative way during the refinement procedure or as a test of the final result. The use of the electrostatic theorem is expected to reduce significantly correlations among static density parameters and parameters describing vibrations, since it is a valuable tool to discriminate between physically reasonable and artificial static electron densities. All of these mentioned suggestions can be applied as internal quality criteria. The last suggestion is based on the idea to initiate the experimental refinement with a set of model parameters, which is, as much as possible close to the final solution. This can be achieved by performing periodic boundary conditions calculations, from which theoretically created files are obtained, which contain the Miller indices (h, k, l) and the respective intensity I. This file is used for a model parameter estimation (refinement), which excludes vibrations. The resulting parameters can be used for the experimental refinement, where, in a first step, the density parameters are fixed to determine the parameters describing vibrations. For a fine tuning, again the electrostatic theorem and the other above mentioned suggestions could be applied. Theoretical predictions should not be biased by the method of computation. Therefore the dependence of the density analyzing tools on the level of calculation (method of calculation/basis set) and on the substituents in complex chemical bonding situations were evaluated in the second part of the present work. A number of compounds containing formal single and double sulfur nitrogen bonds was investigated. For these compounds, experimental data were also available. The calculated data were compared internally and with the experimental results. The internal comparison was drawn with regard to questions of convergency as well as with regard to questions of consistency: The resulting molecular properties from NBO/NRT analyses were found to be very stable, when the geometries were optimized at the respective level of theory. This stability is valid for variations in the methods of calculation as well as for variations in the basis set. Only the individual resonance weights of the contributing Natural Lewis Structures differed considerably depending on the level of calculation and depending on the substituents. However, the deviations were in both cases to a large extent within a limit which preserves the descending order of the leading resonance structure weights. The resulting bond orders, i.e., the total, covalent and ionic bond order from NRT calculations, were not affected by the shift in the resonance weights. The analysis of the bond topological parameters resulted in a discrimination between insensitive parameters and sensitive parameters. The stable parameters do neither depend strongly on the method of calculation nor on the basis set. Only minor variation occurs in the numerical values of these parameters, when the level of calculation is changed or even when other functional groups (H, Me, or tBu) are employed, as long as the methods of calculation do not drop considerably below a standard level. The bond descriptors of the sulfur nitrogen bonds were found to be also stable with respect to the functional groups R = H, R = Me, and R = tBu. Stable parameters are the bond distance, the density at the bond critical point (BCP) and the ratio of distances between the BCP and the nuclei A and B, which varies clearly when considering the formal bond type. For very small basis sets like the 3-21G basis set, this characteristic stability collapses. The sensitive parameters are based on the second derivatives of the density with respect to the coordinates. This is in accordance with the well known fact, that the total second derivative of the density with respect to the coordinates is a strongly oscillating function with positive as well as negative values. A profound deviation has to be anticipated as a consequence of strong oscillations. lambda3, which describes the local charge depletion in the direction of the interaction line, is the most varying parameter. A detailed analysis revealed that the position of the BCP in the rampant edge of the Laplacian distribution is responsible for the sensitivity of the numerical value of lambda3 in formal double bonds. Since the slope of the Laplacian assumes very high values in its rampant edge, a tiny displacement of the BCP leads already to a considerable change in lambda3. This instability is not a failure of the underlying theory, but it yields de facto to a considerable dependence of sensitive bond topological properties on the method of calculation and on the applied basis sets. Since the total second derivative is important to judge on the nature of the bond in the AIM theory (closed shell interactions versus shared interactions), the changes in lambda3 can lead to differing chemical interpretations. The comparison of theoretically derived bond topological properties of various sulfur nitrogen bonds provides the possibility to measure the self consistency of this data set. All data sets clearly exhibit a linear correlation between the bond distances and the density at the BCP on one hand and between the bond distances and the Laplacian values at the BCP on the other hand. These correlations were almost independent of the basis set size. In this context, the linear regression has to be regarded exclusively as a descriptive statistics tool. There is no correlation anticipated a priori. The formal bond type was found to be readily deducible from the theoretically obtained bond topological descriptors of the model systems. In this sense, the bond topological properties are self consistent despite of the numerical sensitivity of the derivatives, as exemplified above. Often, calculations are performed with the experimentally derived equilibrium geometries and not with optimized ones. Applying this approach, the computationally costly geometry optimizations are saved. Following this approach the bond topological properties were calculated using very flexible basis sets and employing the fixed experimental geometry (which, of course, includes the application of tBu groups). Regression coe±cients similar to those from optimized geometries were obtained for correlations between bond distances and the densities at the BCP as well as for the correlation between bond distances and the Laplacian at the BCP, i.e. the approach is valid. However, the data points scattered less and the coe±cient of correlation was clearly increased when geometry optimizations were performed beforehand. The comparison between data obtained from theory and experiment revealed fundamental discrepancies: In the data set of bond topological parameters from the experiment, the behavior of only 2 out of 3 insensitive parameters was comparable to the behavior of the theoretically obtained values, i.e. theoretical and experimental bond distances as well as theoretical and experimental densities at the BCP correlate. From the theoretically obtained data it was easy to deduce the formal bond type from the position of the BCP, since it changed in a systematic manner. The respective experimentally obtained values were almost constant and did not change systematically. For the SN bonds containing compounds, the total second derivative assumes exclusively negative values in the experiment. Due to the different internal behavior, experimentally and theoretically sensitive bond topological values could not be compared directly. The qualitative agreement in the Laplacian distribution, however, was excellent. In the third and last part of this work, the application to chemical systems follows. Formal hypervalent molecules, i.e. molecules where some atoms are considered to hold more than 8 electrons in their valence shell, were investigated. These were compounds containing sulfur nitrogen bonds (H(NtBu)2SMe, H2C{S(NtBu)2(NHtBu)}2, S(NtBu)2 and S(NtBu)3) and a highly coordinated silicon compound. The set of sulfur nitrogen compounds also contained a textbook example for valence expansion, the sulfur triimide. For these molecules, experimental reference values were available from high resolution X-ray experiments. The experimental results were in the case of the sulfur triimide not unique. Furthermore, from the experimental bond topological data no definite conclusion about the formal bonding type could be drawn. The situation of sulfur nitrogen bonds in the above mentioned set of molecules was analyzed in terms of a geometry discussion and by means of a topological analysis. The methyl-substituted isolated molecules served as model compounds. For the interpretation of the bonding situation additional NBO/NRT calculations were preformed for the sulfur nitrogen compounds and an ELF calculation and analysis was performed for the silicon compound. The ELF analysis included not only the presentation and discussion of the ELF-isosurfaces (eta = 0.85), but also the investigation of populations of disynaptic valence basins and the percentage contributions to these populations of the individual atoms when the disynaptic valence basins are split into atomic contributions according to Bader’s partitioning scheme. The question of chemical interest was whether hypervalency is present in the set of molecules or not. In the first case the octet rule would be violated, in the second case Pauling’s verdict would be violated. While the concept of hypervalency is well established in chemistry, the violation of Pauling’s verdict is not. The quantitative numbers of the sensitive bond topological values from theory and experiment were not comparable, since no systematic relationship between the experimentally and theoretically determined sensitive bond descriptors was found. However, the insensitive parameters are in good agreement and the qualitative Laplacian distribution is, with few exceptions, in excellent agreement. The formal bonding type was deduced from experimental and theoretical topological data by considering the number and shape of valence shell charge concentrations in proximity to the sulfur and nitrogen centers. The results from NBO/NRT calculations confirmed the findings. All employed density analyzing tools AIM, ELF and NBO/NRT coincided in describing the bonding situation in the formally hypervalent molecules as highly polar. A comparison and analysis of experimentally and theoretically derived electron densities led consistently to the result, that regarding this set of molecules, hypervalency has to be excluded unequivocally.
The main aim of this work was the classification of highly polar E–N (E = Al, Si, P) and Li–E’ (E’ = C, N, O) bonds in terms of ionic (closed-shell) or covalent (shared) interactions. To answer this question the experimentally determined electron density was analyzed using Bader’s theory of ‘Atoms in Molecules’ (AIM). This allows a quantitative evaluation of properties derived from the electron density, such as the Laplacian, the ellipticitiy and the ratio of the highest charge concentration perpendicular to the bond path, to the largest charge depletion along the bonding vector. Most of these properties were monitored along the entire bonding region and not limited to the BCP as in former studies. The analyses are completed by the calculation of the electronic energy densities Hl at the BCPs and the integration of atomic basins also defined within the AIM theory. The electrostatic potential (ESP) was computed from the multipole parameters to reveal preferred reactive sites of the structures under investigation. Apart from that, the multipole formalism was applied to problematic crystal structures in order to open this method for twinned samples or those including disordered groups in the molecule.