Refine
Has Fulltext
- yes (12)
Is part of the Bibliography
- yes (12)
Year of publication
Document Type
- Doctoral Thesis (12)
Keywords
- Enzym (12) (remove)
Institute
- Graduate School of Life Sciences (3)
- Institut für Physikalische und Theoretische Chemie (3)
- Theodor-Boveri-Institut für Biowissenschaften (2)
- Institut für Hygiene und Mikrobiologie (1)
- Institut für Organische Chemie (1)
- Institut für Pharmazie und Lebensmittelchemie (1)
- Klinik und Poliklinik für Psychiatrie, Psychosomatik und Psychotherapie (1)
- Lehrstuhl für Biochemie (1)
- Rudolf-Virchow-Zentrum (1)
Charakterisierung von NadR : das essentielle Enzym der NAD-Synthese bei Haemophilus influenzae
(2005)
I Zusammenfassung Haemophilus influenzae, ein Gram-negatives, Bakterium der Familie Pasteurellaceae, kann beim Menschen eine Vielzahl an Erkrankungen auslösen: Die bekapselte Stämme, v. a. mit Typ b Kapsel können Cellulitis, septische Arthritis, Epiglottitis und Meningitis verursachen. Die nicht-bekapselte Stämme können Otitis media, Sinusitis, Pneumonie und in selteneren Fällen Bakterämie verursachen. Ein besonderes Merkmal des Metabolismus von H. influenzae ist dessen Unfähigkeit Nikotinamid-Adenin-Dinukleotid (NAD+) de novo zu synthetisieren. Daher sind die Enzyme bzw. Transporter, die an NAD+ Aufnahme und Resynthese beteiligt sind, als putative antimikrobielle Ziele von Interesse. In unserer Arbeitsgruppe konnte gezeigt werden, dass NAD+ zu Nikotinamidribosyl degradiert werden muss, bevor es in die Zelle aufgenommen werden kann. Auch Proteine, die an der Degradation des exogenen NAD+ zu Nikotinamidribosyl und dessen anschließender Aufnahme in die Zelle verantwortlich sind, konnten identifiziert und charakterisiert werden. Wie Nikotinamidribosyl im Cytoplasma wiederum zu NAD+ synthetisiert wird, ist auch erst kürzlich geklärt worden: für NadR konnte sowohl eine Ribosyl-Nukleotid-Kinase (RNK) Aktivität als auch eine Nikotinamid-Mononukleotid-Adenylyltransferase (NMNAT) Aktivität in vitro gezeigt werden. Die Kristallstruktur von hiNadR im Komplex mit NAD+ wurde auch aufgeklärt. In dieser Arbeit sollte NadR, insbesondere dessen RNK Domäne, in vivo und in vitro näher charakterisiert werden. Um zu untersuchen, ob beide Domänen in vivo essentiell sind, wurden Deletionsmutanten erzeugt, bei welchen die komplette bzw. der C-terminale Teil der RNK Domäne fehlten. Diese Deletionen konnten im nadV+ Hintergrund erzeugt werden. Die Deletionen konnten in H. influenzae nur zusammen mit dem nadV-Gen transferiert werden oder alternativ nur in die Zellen, die mit pNadRKan Plasmid transformiert wurden. Dies verdeutlicht, dass nicht nur die NMNAT Domäne sondern auch die RNK Domäne bzw. sogar nur wenige C-terminal fehlende Aminosäuren des NadR Proteins essentiell für die Lebensfähigkeit von H. influenzae sind. Gleichzeitig zeigen diese Experimente, dass die RNK-Domäne in Anwesenheit von NadV redundant ist. Ein weiterer Phänotyp der RNK-Deletionsmutante zeigte sich beim Nikotinamidribosyl-Transport. Im Gegensatz zum Wt, welcher ca. 60-80% des 14C-Nikotinamidribosyls aufnahm, konnte für die RNK-Deletionsmutante nur 2-5% Aufnahme gemessen werden. Dies konnte durch das pNadRKan Plasmid komplementiert werden. Weiterhin wurde festgestellt, dass spontan Aminopyridin-resistente H. influenzae Zellen Mutationen im nadR Gen haben, insbesondere im Walker A-Motif (P-Loop) der RNK Domäne. Zusätzlich konnte in dieser Arbeit gezeigt werden, dass NadR aus Aminopyridin und ATP Aminopyridin-Adenin-Dinukleotid synthetisieren kann. Somit konnte gezeigt werden, dass die wachstumshemmende Wirkung eigentlich durch das aus Aminopyridin synthetisierte Aminopyridin-Adenin-Dinukleotid entsteht, welches NAD+ in Redox-Reaktionen verdrängt, wodurch es letztendlich zum Stillstand des Metabolimus kommt. Durch Einführen von gezielten AS-Substitutionen im Walker A und B Motif und in der LID-Domäne von NadR, konnten einige Aminosäuren identifiziert werden, welche essentiell für die Aktivität der RNK Domäne sind. Alle Aminosäuren-Substitutionen führten zum Verlust der RNK Aktivität, die NMNAT Aktivität jedoch war nicht beeinträchtigt. Desweiteren wurden diese NadR Punktmutanten in vivo untersucht. Für alle konnte eine signifikante Defizienz in der Nikotinamidribosyl-Aufnahme beobachtet werden, die gemessene Aufnahme lag im Bereich der RNK-Deletionsmutante. Dadurch konnte eine direkte Korrelation zwischen der RNK Aktivität und der Nikotinamidribosyl-Aufnahme gezeigt werden. In weiteren in vitro Experimenten konnte für NadR eine Feedback-Inhibition durch das NAD+ gezeigt werden, wobei NAD+ in erster Linie die RNK Domäne von NadR inhibiert. Eine graduelle Erhöhung der NAD+ Konzentration führte in den in vitro Assays zu einer graduellen Abnahme der RNK. Bei der NMNAT Aktivität jedoch zeigte sich keine signifikante Inhibition in Anwesenheit von NAD+. Begleitende in vivo Experimente, zeigten eine 2/3 Reduktion der Nikotinamidribosyl-Aufnahme bei den Zellen, die mit NAD+ inkubiert wurden, d. h. höhere intrazelluläre NAD+ Konzentration hatten. Für die genauere Analyse der Feedback-Inhibition durch NAD+ wurden weitere Punktmutanen hergestellt. Bei zwei der Punktmutanten wurde eine Beeinträchtigung der NadR-Aktivität beobachtet, daher wurden diese Punktmutanten von weiteren Analysen im Bezug auf NAD+-Feedback Inhibition ausgeschlossen. Eine Mutante (NadRW256F) jedoch, zeigte ähnliche Aktivität wie das Wt-NadR. In Anwesenheit von NAD+ wurde die RNK Aktivität dieser Punktmutante, im Gegensatz zum Wt-Protein, kaum gehemmt. Dadurch konnte W256 als eine der Aminosäuren identifiziert werden, die an der Vermittlung der NAD+-bedingten Inhibition der RNK-Domäne beteiligt ist.
KasA is a key enzyme which plays an essential part in the biosynthetic pathway of mycolic acids, the building block of cell wall in Mycobacterium tuberculosis. Its importance was demonstrated by the finding that the depletion of KasA leads to the cell lysis of Mycobacterium tuberculosis. Since Mycobacterium tuberculosis is a pathogen of tuberculosis, the second leading cause of death from an infectious disease worldwide, KasA has drawn attention as one of the attractive drug targets against tuberculosis. Due to the emergence of extensively drug-resistant strains which make most of the known antibiotics for treating tuberculosis ineffective, it became an urgent issue to develop new drugs against tuberculosis. In chapter 3.1, the protonation state of the catalytic residues in the resting state was mainly addressed. The FEP computation and MD simulations were employed for this investigation, and the results showed that the zwitterionic state is most probable. To underpin this conclusion with more solid data, The PESs for the proton transfer between the neutral and zwitterionic state were computed in the context of QM/MM. However, due to the strong dependency of the QM/MM optimization on the initial structure, it was not possible to obtain consistent results from these computations. To circumvent this problem, QM/MM based umbrella sampling was carried out with a semi-empirical method (RM1), and the resulting PMF surface indicated that the zwitterionic state is more stable than the neutral state. In chapter 3.2, the protonation state of significant residues in the acyl-enzyme state was investigated. Unlike other catalytic residues, the protonation state of His311 is ambiguous in the acyl-enzyme state, and different decarboxylation mechanisms can be derived depending on the protonation state of His311 in the acyl-enzyme state. Therefore, FEP computations were carried out to find most probable protonation state of His311 in terms of free energy, and the results showed that the pKa value at Nδ is considerably lowered by the enzyme environment while that of Nε is not. Additionally, the PMF profiles for the proton transfer between Lys340 and Glu354 were computed using QM/MM based umbrellas sampling method, and the results showed that the property of the Lys340/Glu354 pair is neutral rather than ionic when His311 is protonated at Nε. Moreover, a relatively larger ionic character of the Lys340/Glu354 pair when His311 is doubly protonated provides a valuable insight into how the Lys340/Glu354 pair plays a role in shifting the protonated state from Nδ to Nε in His311 after the acyl-transfer step. Overall, the results demonstrated that His311 is neutral and protonated at Nε, and the Lys340/Glu354 pair is also neutral in the acyl-enzyme state. Those computational results lead to the conclusion that the decarboxylation reaction is facilitated by an oxyanion hole which is comprised of two catalytic histidines. In chapter 3.3, the protonation state of catalytic residues in the resting state was revisited because a recent benchmark study showed that the employed semi-empirical method (RM1) in chapter 3.1 tends to overestimate the stabilization of the zwitterionic state. Furthermore, the Lys340/Glu354 pair was considered as purely ionic in chapter 3.1, while it actually has a mixed neutral and ionic character as demonstrated in chapter 3.2. The new investigations employed a larger QM region including the Lys340/Glu354 pair with the BLYP/6-31G** approach, which was proven to be accurate enough for the present purpose by benchmark computations. The new results from the QM/MM MD and FEP computations indicated the catalytic residues to be neutral most probably in the resting state, and this in turn brought up the question how KasA can be activated to initiate the catalytic reaction. On the basis of the results from the MD simulations and FEP computations for the His311Ala mutant in chapter 3.1, we hypothesized that the open conformation of Phe404 would trigger the activation of the catalytic residues by the formation of a strong hydrogen bond. The QM/MM MD simulation proved that the activation of the catalytic residues can indeed be accomplished by the open conformation of Phe404 we suggested, and the corresponding force field based PMF profile also indicated that this conformational change is energetically feasible. The distribution of hydrophilic and hydrophobic residues in the malonyl binding pocket in conjunction with our computational results further provided a valuable insight into the detailed process how the catalytic residues is activated upon the substrate entering.
The SARS virus is the etiological agent of the severe acute respiratory syndrome, a deadly disease that caused more than 700 causalities in 2003. One of its viral proteins, the SARS coronavirus main protease, is considered as a potential drug target and represents an important model system for other coronaviruses. Despite extensive knowledge about this enzyme, it still lacks an effective anti-viral drug. Furthermore, it possesses some unusual features related to its active-site region. This work gives atomistic insights into the SARS coronavirus main protease and tries to reveal mechanistic aspects that control catalysis and inhibition. Thereby, it applies state-of-the-art computational methods to develop models for this enzyme that are capable to reproduce and interpreting the experimental observations. The theoretical investigations are elaborated over four main fields that assess the accuracy of the used methods, and employ them to understand the function of the active-site region, the inhibition mechanism, and the ligand binding. The testing of different quantum chemical methods reveals that their performance depends partly on the employed model. This can be a gas phase description, a continuum solvent model, or a hybrid QM/MM approach. The latter represents the preferred method for the atomistic modeling of biochemical reactions. A benchmarking uncovers some serious problems for semi-empirical methods when applied in proton transfer reactions. To understand substrate cleavage and inhibition of SARS coronavirus main protease, proton transfer reactions between the Cys/His catalytic dyad are calculated. Results show that the switching between neutral and zwitterionic state plays a central role for both mechanisms. It is demonstrated that this electrostatic trigger is remarkably influenced by substrate binding. Whereas the occupation of the active-site by the substrate leads to a fostered zwitterion formation, the inhibitor binding does not mimic this effect for the employed example. The underlying reason is related to the coverage of the active-site by the ligand, which gives new implications for rational improvements of inhibitors. More detailed insights into reversible and irreversible inhibition are derived from in silico screenings for the class of Michael acceptors that follow a conjugated addition reaction. From the comparison of several substitution patterns it becomes obvious that different inhibitor warheads follow different mechanisms. Nevertheless, the initial formation of a zwitterionic catalytic dyad is found as a common precondition for all inhibition reactions. Finally, non-covalent inhibitor binding is investigated for the case of SARS coranavirus main protease in complex with the inhibitor TS174. A novel workflow is developed that includes an interplay between theory and experiment in terms of molecular dynamic simulation, tabu search, and X-ray structure refinement. The results show that inhibitor binding is possible for multiple poses and stereoisomers of TS174.
Deubiquitinases are regulators of the ubiquitin proteasome system that counteract the ubiquitination cascade by removing ubiquitin from substrates and cleaving ubiquitin chains. Due to their involvment in various important pathways, they are associated with several diseases and may thus present promising drug targets. The two related ubiquitin specific proteases USP25 and USP28 share a highly conserved amino acid sequence but perform distinct biological functions. USP28 plays roles in cell cycle regulation and was also linked to several types of cancer. It adopts oncogenic functions by rescuing the oncoproteins MYC and JUN from proteasomal degradation, which is induced by the E3-ligase SCF (FBW7). Opposingly, USP28 also regulates the stability of the tumor suppressor FBW7 itself. USP25 contributes to a balanced innate immune system by stabilizing TRAF3 and TRAF6 and lately was found to promote Wnt-signaling by deubiquitinating TNKS.
Due to the high level of identity of both proteases, a recent attempt to inhibit USP28 led to cross reactivity against USP25. In our study, we characterized both USP25 and USP28 structurally and functionally using x-ray crystallography, biochemical as well as biophysical approaches to determine similarities and differences that can be exploited for the development of specific inhibitors.
The crystal structure of the USP28 catalytic domain revealed a cherry-couple like dimer that mediates self-association by an inserted helical subdomain, the USP25/28 catalytic domain inserted domain (UCID). In USP25, the UCID leads to formation of a tetramer composed of two interlinked USP28-like dimers. Structural and functional analysis revealed that the dimeric USP28 is active, whereas the tetrameric USP25 is auto inhibited. Disruption of the tetramer by a cancer-associated mutation or a deletion-variant activates USP25 through dimer formation in in vitro assays and leads to an increased stability of TNKS in cell studies. Furthermore, in vitro data showed that neither ubiquitin nor substrate binding led to the activation of the USP25 tetramer construct. With the structure of the C-terminal domain of USP25, we determined the last unknown region in the enzyme as a separately folded domain that mediates substrate interactions.
Combined the structures of the USP25 and USP28 catalytic domains and the functional characterization of both enzymes provide novel insights into the regulation of USPs by oligomerization. Furthermore, we identified individual features of each protease that might be explored for the development of specific small molecule inhibitors.
Ungeachtet der enormen Entwicklung in Krebsdiagnostik und -Therapie in den letzten Jahren, sind vollständige Heilungsaussichten weiterhin gering und die aktuellen Behandlungsmethoden oftmals mit schwerwiegenden Nebeneffekten verbunden. Aufgrund dessen sind alternative Behandlungsmethoden unbedingt erforderlich und führten zu einer zunehmenden Bedeutung des Vaccinia-Virus als onkolytisches Virus in der Krebstherapie. In der vorliegenden Arbeit wurden zwei mögliche Therapieansätze zur Verstärkung der onkolytischen Effekte in humanen Tumormodellen untersucht. Die Kombination einer gene-directed enzyme prodrug Therapie (GDEPT) mit dem onkolytischen Vaccinia-Virus GLV 1h68 sollte zur Selektivitätssteigerung eines ß-Galaktosidase-aktivierbaren, cytotoxisch-aktiven Drugs führen. Darüber hinaus diente das für MCP-1 codierende Vaccinia-Virus GLV-1h80, zielend auf eine Cytokin-vermittelten Immuntherapie, als Vektor zur spezifischen Beeinflussung des intratumoralen Chemokin-Netzwerks. Im Zuge der GDEPT wurde in dieser Arbeit ein, durch enzymatische Deglykosylierug aktivierbares Prodrug, basierend auf dem cytotoxischem Antibiotikum Duocarmycin SA verwendet. Durch eine Infektion mit GLV-1h68 und einer resultierenden Expression des aktivierenden Enzyms ß-Galaktosidase, sollte eine Umwandlung des Prodrugs in ein cytotoxisches Drug erfolgen. In vitro Infektionsstudien zeigten ein nahezu identisches Replikationsverhalten des Vaccinia-Virus GLV-1h68 und des als Kontrollvirus dienenden rVACV GLV-1h43 in humanen GI-101A-Brustkrebszellen. Die Expression der beiden Reporter-Gene Ruc-GFP sowie ß-Galaktosidase konnten auf Protein-Ebene und mittels RT-PCR nach Infektion mit GLV-1h68 nachgewiesen werden. GLV-1h43-Infektion von GI-101A-Zellen führte zu GFP-Expression, jedoch nicht zur Expression des Enzyms ß Galaktosidase. Untersuchung der Enzym-Aktivität in Zelllysaten und Zellkultur-Überständen zeigten nach Infektion mit GLV 1h68 steigende Menge zellulär assoziierter und freier ß-Galaktosidase. Des Weiteren wurde durch Koinkubation von GI-101A-Zellen mit Virus-freien, ß Galaktosidase-haltigen Zelllysaten bzw. –überständen und Prodrug eine Aktivierung des Prodrugs durch das Virus codierte Enzym nachgewiesen. Diese Koinkubation führte zur Abtötung der Zellen. Nach Inkubation mit Proben mock- oder GLV 1h43-infizierter Zellen konnte keiner Veränderung der Proliferationsrate von GI-101A-Zellen gefunden werden. Kombinierte Behandlung von GI 101A-Zellen mit Viren des Stammes GLV 1h68 und Prodrug führte zu starken Synergieeffekten bei der Abtötung der Zellen und wies einen Bystander Effekt der Kombinationstherapie nach. Dieser konnte in 4 weiteren humanen und 2 Hunde-Brustkrebszellen bestätigt werden. Der erzielte Bystander-Effekt zeigt, dass es nach Virus-induzierter ß-Galaktosidase-Expression in GLV 1h68-infizierten Zellen zu einer enzymatischen Spaltung des Prodrugs in das cytotoxische seco-Analogon des Antibiotikums Duocarmycin SA kommt. Durch die Membrangängigkeit des Drugs konnte auch in angrenzenden uninfizierten Zellen eine Wirkung erzielt werden. Anhand von Expressionsanalysen an Apoptose-assoziierten Proteinen, wie PARP und Caspasen, wurde eine Wirkung des Prodrugs über den intrinsischen Apoptose-Signalweg nachgewiesen. In athymischen Nude-Mäusen durchgeführte Replikationsanalysen und X-Gal-Färbungen GLV 1h68 infizierter Tumore nach Prodrug-Behandlung zeigten, dass GLV-1h68 ungeachtet der simultanen Behandlung mit Prodrug im Tumorgewebe repliziert und es nicht zur Anreicherung lacZ-negativer Virusmutanten kommt. Es konnten, durch Prodrug-Behandlung und einer simultanen Expression aktiver ß Galaktosidase, starke synergistische Effekte und eine signifikante Steigerung der Tumorregression erzielt werden. Da die Kombinationstherapie zu keinerlei Unterschieden in Gewicht und Gesundheitszustand behandelter Versuchstiere führte, konnte eine systemische Toxizität außerhalb des Tumorgewebes ausgeschlossen werden. Verschiedene Zelllinien weisen Unterschiede in ihrer Sensitivität gegenüber der onkolytischen Aktivität von Vaccinia-Virus GLV-1h68 auf. Während einige Zelllinien trotz Virus-Behandlung unverändertes Proliferationsverhalten zeigen (non- oder poor-responder), führt diese Behandlung in anderen Zelllinien zu einer vollständigen Tumorregression (responder). In Anbetracht dieser Unterschiede wurden in dieser Arbeit die Effekte einer induzierten Expression des murinen Chemokins MCP-1 in GI-101A-Tumoren (responder) und HT29-CBG-Tumoren (poor-responder) untersucht. MCP-1 zeichnet sich durch seine chemotaktischen Eigenschaften gegenüber mononukleärer Zellen aus und führt zu pleiotropen Tumor-Effekten. Replikationsstudien am Virus GLV-1h80 und des als Kontrollvirus dienenden rVACV GLV-1h68 zeigten, dass aus der Expression des Fremd-Gens mcp-1 sowohl in vitro als auch in vivo keinerlei negativen Effekte auf das Replikationsverhalten in humanen GI-101A- und HT29-CBG-Zellen resultieren. Durch Real-time Monitoring der GFP-Expression im Tumorgewebe lebender Tiere konnte zunächst eine mit dem Infektionsverlauf zunehmende Signalstärke beobachtet werden, welche dann 42 dpi an Intensität verlor. Toxizität und schädliche Nebeneffekte durch Infektion mit den beiden rVACV konnten anhand der viralen Titer in den Organen der Maus ausgeschlossen werden. Die Titer wiesen auf eine ausschließlich auf das Tumorgewebe begrenzte Replikation der Viren nach Injektion in Tumor-tragende Tiere hin. Die Expression des Chemokins MCP-1 wurde sowohl auf transkriptioneller als auch auf translationeller Ebene in GLV-1h80-inifzierten Zellen und im Tumorgewebe GLV 1h80-injizierter Mäuse nachgewiesen. Nach Infektion mit GLV-1h80 konnte eine mit dem Infektionsverlauf zunehmende MCP-1-Expression gezeigt werden. Dabei wurde zudem deutlich, dass nicht nur eine GLV-1h80-Infektion in vivo zu einer Zunahme der intratumoralen MCP-1-Expression führte, sondern eine Vaccinia-Virus-Infektion allein einen Anstieg des Chemokins zu bewirken vermag. Eine Quantifizierung durch ELISA machte Konzentrationsunterschiede von MCP-1 zwischen den Tumormodellen GI-101A und HT29-CBG deutlich. Sowohl in vitro als auch in vivo führte ein GLV-1h80-Infektion zu deutlich niedrigeren Konzentrationen im HT29-CBG-Kolon-Adenokarzinommodell. Ein Nachweis murinen MCP-1 in Blutseren Tumor-tragender Tiere zeigte eine für therapeutische Effekte erwünschte systemische Freisetzung des intratumoral durch die Infektion mit GLV-1h80 gebildeten Chemokins MCP-1. Durch immunhistologische Untersuchungen GLV-1h80-infizierter Zellen und Tumoren konnte diese, mit dem Infektionsverlauf zunehmende MCP-1-Expression bestätigt werden. Die funktionelle Aktivität des rekombinanten Proteins wurde anhand TNF-α-spezifischer ELISA-Analysen überprüft. Dabei zeigte sich eine erhöhte Expression dieses proinflammatorischen Cytokins in GI-101A-Tumoren nach Infektion mit GLV-1h80. Dagegen konnte keine Steigerung der Expression im HT29-CBG-Tumorgewebe nachgewiesen werden. Ein Nachweis des durch proinflammatorische Immunzellen exprimierten Oberlflächenproteins CD14 zeigte ebenfalls einen Anstieg nach Infektion mit GLV-1h80. Auch diese veränderte Expression blieb im poor-Responder-Modell HT29-CBG aus. Die steigende intratumorale Expression der beiden Proteine in GI-101A-Tumoren nach GLV 1h80-Infektion lässt auf eine Zunahme pro-inflammatorischer Immunzellen, basierend auf einer Virus-induzierten MCP-1-Expression schließen. Ein Monitoring der Tumorprogression nach Implantation von GI 101A-Zellen und Injektion der rVACV GLV-1h80 und GLV-1h68 bzw. einer PBS-Injektion führte nach einer anfänglichen Zunahme des Tumorwachstums schließlich bei beiden Viren zu einer Tumorregression. Jedoch konnte durch die GLV-1h80-vermittelte MCP-1-Expression eine Verstärkung der onkolytischen Effekte erzielt werden, welche sich durch eine signifikante Abnahme des Tumorvolumens zeigte. Im HT29-CBG-Modell führten die therapeutischen Effekte durch rVACV GLV-1h80 zwar zu keiner Regression des Tumors, jedoch zeigte sich auch in diesem humanen Tumormodell eine Verstärkung der onkolytischen Effekte nach GLV-1h80-Infektion im Vergleich zu einer GLV 1h68-Behandlung. Durch die GLV-1h80-induzierte Expression des Chemokins MCP-1 konnte somit eine Hemmung des Tumorwachstums auch im poor-Responder-Modell HT29-CBG erzielt werden. Sowohl die Verwendung eines ß-Galaktosidase-aktivierbaren Prodrugs im Zuge einer GDEPT, als auch die Beeinflussung des intratumoralen Chemokin-Netzwerks durch Expression des Chemokins MCP-1 führten in dieser Arbeit zu positiven Synergismus-Effekten in der onkolytischen Virustherapie. Durch künftige Konstruktion eines rVACV, welches sowohl die Expression des Chemokins MCP-1, als auch des prodrug-aktivierenden Enzyms ß-Galaktosidase im Tumorgewebe induziert, könnte in Kombination mit einer Prodrug-Behandlung eine zusätzliche Verstärkung der Effekte erzielt und möglicherweise eine erfolgreiche Virustherapie in bisher schwach ansprechenden poor- bzw. non-Responder-Modellen ermöglicht werden.
Im Katabolismus methylverzweigter Fettsäuren spielt die alpha-Methylacyl-CoA-Racemase eine wichtige Rolle, indem sie die (R)- und (S)-Isomere von alpha-methylverzweigten Fettsäuren als Coenzym A Thioester racemisiert. Methylverzweigte Fettsäuren entstehen beim Abbau von Isoprenoiden und werden darüber hinaus auch von vielen Organismen, wie z.B. Mycobakterien, synthetisiert. Die Hauptaufgabe der Racemase ist aber vermutlich in der Biosynthese von Gallensäuren zu sehen. Das Ziel der vorliegenden Arbeit war es, die alpha-Methylacyl-CoA-Racemase aus humanem Gewebe zu reinigen und zu charakterisieren sowie ihre physiologische Rolle im Katabolismus verzweigtkettiger Fettsäuren und der Gallensäurebiosynthese zu untersuchen. Die alpha-Methylacyl-CoA-Racemase wurde aus humanem Gewebe zur Homogenität gereinigt, umfassend biochemisch charakterisiert und zur genauen molekularbiologischen Analyse in E.coli kloniert. Die Aktivität der Racemase wurde anhand der [³H]H2O-Freisetzung aus [alpha-³H]-a-Methylacyl-CoAs bestimmt. Die humane Racemase ist in der aktiven Form ein monomeres Protein und besteht aus 382 Aminosäuren. Als Substrate akzeptiert das Enzym ein breites Spektrum von alpha-Methylacyl-CoAs. Neben den Coenzym A-Thioestern alpha-methylverzweigter Fettsäuren, wie Pristansäure, werden auch CoA-Ester von Steroidderivaten, z.B. des Gallensäureintermediats Trihydroxycoprostansäure, und aromatischen Phenylpropionsäuren, wie dem Analgetikum Ibuprofen, umgesetzt. Freie Fettsäuren, geradkettige oder beta-methylverzweigte Acyl-CoAs werden nicht racemisiert. Die alpha-Methylacyl-CoA-Racemase ist im Menschen zu ca. 80 Prozent auf die Peroxisomen und ca. 20 Prozent auf die Mitochondrien verteilt, wobei entsprechende peroxisomale (PTS 1) und mitochondriale (MTS) Transportsignale die Lokalisation bestimmen. Die vollständige cDNA-Sequenz der humanen a-Methylacyl-CoA-Racemase hat eine Gesamtlänge von 2039 Basenpaaren mit einem offenen Leseraster von 89 - 1237 bp. Das Startcodon ATG ist in eine klassische Kozak-Sequenz zum Translationsstart eingebettet. Die Protein endet am C-Terminus mit dem Sequenzmotiv –KASL, das dem peroxisomalen Transportsignal (PTS I) einiger Säugetierkatalasen entspricht. Aufgrund alternativer Polyadenylierung sind in allen untersuchten menschlichen Geweben Transkripte von 1,6 kb bzw. 2,0 kb zu finden. Es liegt keine gewebsabhängige Polyadenylierung vor, die Racemase wird aber gewebsspezifisch exprimiert (besonders stark in Leber und Niere). Das humane Racemasegen liegt auf dem kurzen Arm des Chromosoms 5 nahe am Centromer (5p1.3), im Intervall von D5S651 (46,6 cM) und D5S634 (59.9 cM).
Several lines of evidence implicate a dysregulation of tryptophan hydroxylase (TPH)-dependent serotonin (5-HT) synthesis in emotions and stress and point to their potential relevance to the etiology and pathogenesis of various neuropsychiatric disorders. However, the differential expression pattern of the two isoforms TPH1 and TPH2 which encode two forms of the rate-limiting enzyme of 5-HT synthesis is controversial. Here, a comprehensive spatio-temporal analysis clarifies TPH1 and TPH2 expression during pre- and postnatal development of the mouse brain and in adult human brain as well as in peripheral organs including the pineal gland. Four different methods (real time PCR, in situ hybridization, immunohistochemistry and Western blot analysis) were performed to systematically control for tissue-, species- and isoform-specific expression on both the pre- and posttranslational level. TPH2 expression was consistently detected in the raphe nuclei, as well as in fibres in the deep pineal gland and in the gastrointestinal tract. Although TPH1 expression was found in these peripheral tissues, no significant TPH1 expression was detected in the brain, neither during murine development, nor in mouse and human adult brain. Also under conditions like stress and clearing the tissue from blood cells, no changes in expression levels were detectable. Furthermore, the reuptake of 5-HT into the presynaptic neuron by the serotonin transporter (SERT) is the major mechanism terminating the neurotransmitter signal. Thus, mice with a deletion in the Sert gene (Sert KO mice) provide an adequate model for human affective disorders to study lifelong modified 5-HT homeostasis in interaction with stressful life events. To further explore the role of TPH isoforms, Tph1 and Tph2 expression was studied in the raphe nuclei of Sert deficient mice under normal conditions as well as following exposure to acute immobilization stress. Interestingly, no statistically significant changes in expression were detected. Moreover, in comparison to Tph2, no relevant Tph1 expression was detected in the brain independent from genotype, gender and treatment confirming expression in data from native animals. Raphe neurons of a brain-specific Tph2 conditional knockout (cKO) model were completely devoid of Tph2-positive neurons and consequently 5-HT in the brain, with no compensatory activation of Tph1 expression. In addition, a time-specific Tph2 inducible (i) KO mouse provides a brain-specific knockdown model during adult life, resulting in a highly reduced number of Tph2-positive cells and 5-HT in the brain. Intriguingly, expression studies detected no obvious alteration in expression of 5-HT system-associated genes in these brain-specific Tph2 knockout and knockdown models. The findings on the one hand confirm the specificity of Tph2 in brain 5-HT synthesis across the lifespan and on the other hand indicate that neither developmental nor adult Tph2-dependent 5-HT synthesis is required for normal formation of the serotonergic system, although Tph1 does not compensate for the lack of 5-HT in the brain of Tph2 KO models. A further aim of this thesis was to investigate the expression of the neuropeptide oxytocin, which is primarily produced in the hypothalamus and released for instance in response to stimulation of 5-HT and selective serotonin reuptake inhibitors (SSRIs). Oxytocin acts as a neuromodulator within the central nervous system (CNS) and is critically involved in mediating pain modulation, anxiolytic-like effects and decrease of stress response, thereby reducing the risk for emotional disorders. In this study, the expression levels of oxytocin in different brain regions of interest (cortex, hippocampus, amygdala, hypothalamus and raphe nuclei) from female and male wildtype (WT) and Sert KO mice with or without exposure to acute immobilization stress were investigated. Results showed significantly higher expression levels of oxytocin in brain regions which are involved in the regulation of emotional stimuli (amygdala and hippocampus) of stressed male WT mice, whereas male Sert KO as well as female WT and Sert KO mice lack these stress-induced changes. These findings are in accordance with the hypothesis of oxytocin being necessary for protection against stress, depressive mood and anxiety but suggest gender-dependent differences. The lack of altered oxytocin expression in Sert KO mice also indicates a modulation of the oxytocin response by the serotonergic system and provides novel research perspectives with respect to altered response of Sert KO mice to stress and anxiety inducing stimuli.
Ubiquitination is an important post-translational modification that maintains cellular homeostasis by regulating various biological processes. Deubiquitinases (DUBs) are enzymes that reverse the ubiquitination process by catalyzing the removal of ubiquitin from a substrate. Abnormal expression or function of DUBs is often associated with the onset and progression of various diseases, including cancer. Ubiquitin specific proteases (USPs), which constitute the largest family of DUBs in humans, have become the center of interest as potential targets in cancer therapy as many of them display increased activity or are overexpressed in a range of malignant tumors or the tumor microenvironment.
Two related members of the USP family, USP28 and USP25, share high sequence identities but play diverse biological roles. USP28 regulates cell proliferation, oncogenesis, DNA damage repair and apoptosis, whereas USP25 is involved in the anti-viral response, innate immunity and ER-associated degradation in addition to carcinogenesis. USP28 and USP25 also exhibit different oligomeric states – while USP28 is a constitutively active dimer, USP25 assumes an auto-inhibited tetrameric structure. The catalytic domains of both USP28 and USP25 comprise the canonical, globular USP-domain but contain an additional, extended insertion site called USP25/28 catalytic domain inserted domain (UCID) that mediates oligomerization of the proteins. Disruption of the USP25 tetramer leads to the formation of an activated dimeric protein. However, it is still not clear what triggers its activation.
Due to their role in maintaining and stabilizing numerous oncoproteins, USP28 and USP25 have emerged as interesting candidates for anti-cancer therapy. Recent advances in small-molecular inhibitor development have led to the discovery of relatively potent inhibitors of USP28 and USP25. This thesis focuses on the structural elucidation of USP28 and the biochemical characterization of USP28/USP25, both in complex with representatives of three out of the eight compound classes reported as USP28/USP25-specific inhibitors. The crystal structures of USP28 in complex with the AZ compounds, Vismodegib and FT206 reveal that all three inhibitor classes bind into the same allosteric pocket distant from the catalytic center, located between the palm and the thumb subdomains (the S1-site). Intriguingly, this binding pocket is identical to the UCID-tip binding interface in the USP25 tetramer, rendering the protein in a locked, inactive conformation. Formation of the binding pocket in USP28 requires a shift in the helix α5, which induces conformational changes and local distortion of the binding channel that typically accommodates the C-terminal tail of Ubiquitin, thus preventing catalysis and abrogating USP28 activity. The key residues of the USP28-inhibitor binding pocket are highly conserved in USP25. Mutagenesis studies of these residues accompanied by biochemical and biophysical assays confirm the proposed mechanism of inhibition and similar binding to USP25.
This work provides valuable insights into the inhibition mechanism of the small molecule compounds specifically for the DUBs USP28 and USP25. The USP28-inhibitor complex structures offer a framework to develop more specific and potent inhibitors.
The present work illustrates the structural and biochemical characterization of two diverse proteins, BadI and MenD from Rhodopseudomonas palustris and Staphylococcus aureus, respectively.
BadI or 2-ketocyclohexanecarboxyl-CoA is one of the key enzymes involved in the anaerobic degradation of aromatic compounds. The degradation of aromatic compounds is a vital process for the maintenance of the biogeochemical carbon cycle and bioremediation of xenobiotic compounds, which if present at higher concentrations can cause potential hazards to humans. Due to the relatively inert nature of aromatic compounds, enzymes catalyzing their degradation are of special interest for industrial applications. BadI is one of the key enzymes involved in the anaerobic degradation of aromatic compounds into an aliphatic moiety.
The major focus of this study was to provide mechanistic insights into the reaction catalyzed by BadI. BadI belongs to the crotonase superfamily and shares high sequence homology with the family members of MenB or dihydroxynaphthoate synthase. BadI is known to catalyze the cleavage of the cyclic ring of 2-ketocyclohexane carboxyl-CoA by hydrolyzing the C-C bond leading to the formation of the aliphatic compound pimelyl CoA. On the other hand MenB catalyzes the condensation reaction of o-succinylbenzoyl-CoA to dihydroxylnaphthoyl-CoA. A comprehensive amino acid sequence analysis between BadI and MenB showed that the active site residues of MenB from Mycobacterium tuberculosis (mtMenB) are conserved in BadI from Rhodopseudomonas palustris. MenB is involved in the menaquinone biosynthesis pathway and is a potential drug target against Mycobacterium tuberculosis as it has no known human homologs. Due to the high homology between MenB and BadI and the inability to obtain MenB-inhibitor complex structures we extended our interest to BadI to explore a potential substitute model for mtMenB as a drug target.
In addition, BadI possesses some unique mechanistic characteristics. As mentioned before, it hydrolyzes the substrate via a retro Dieckmann’s reaction contrasting its closest homolog MenB that catalyzes a ring closing reaction through a Dieckmann’s reaction. Nevertheless the active site residues in both enzymes seem to be highly conserved. We therefore decided to pursue the structural characterization of BadI to shed light on the similarities and differences between BadI and MenB and thereby provide some insights how they accomplish the contrasting reactions described above.
We determined the first structures of BadI, in its apo and a substrate mimic bound form. The crystal structures revealed that the overall fold of BadI is similar to other crotonase superfamily members. However, there is no indication of domain swapping in BadI as observed for MenB. The absence of domain swapping is quite remarkable because the domain swapped C-terminal helical domain in MenB provides a tyrosine that is imperative for catalysis and is also conserved in the BadI sequence. Comparison of the active sites revealed that the C-terminus of BadI folds onto its core in such a way that the conserved tyrosine is located in the same position as in MenB and can form interactions with the ligand molecule. The structure of BadI also confirms the role of a serine and an aspartate in ligand interaction, thus validating that the conserved active site triad participates in the enzymatic reaction. The structures also reveal a noteworthy movement of the active site aspartate that adopts two major conformations. Structural studies further illuminated close proximity of the active site serine to a water and chlorine molecule and to the carbon atom at which the carbonyl group of the true substrate would reside. Biochemical characterization of BadI using enzyme kinetics validated that the suggested active site residues are involved in substrate interaction. However, the role of these residues is very distinct, with the serine assuming a major role. Thus, the present work ascertain the participation of putative active site residues and demonstrates that the active site residues of BadI adopt very distinctive roles compared to their closest homolog MenB.
The MenD protein also referred to as SEPHCHC (2-succinyl-5-enolpyruvyl-6- hydroxy-3-cyclohexene-1-carboxylic acid) synthase is one of the enzymes involved in menaquinone biosynthesis in Staphylococcous aureus. Though S. aureus is usually considered as a commensal it can act as a remarkable pathogen when it crosses the epithelium, causing a wide spectrum of disorders ranging from skin infection to life threatening diseases. Small colony variants (SCVs), a slow growing, small sized subpopulation of the bacteria has been associated with persistent, recurrent and antibiotic resistant infections. These variants show autotrophy for thiamine, menaquinone or hemin. Menaquinone is an essential component in the electron transport pathway in gram-positive organisms. Therefore, enzymes partaking in this pathway are attractive drug targets against pathogens such as Mycobacterium tuberculosis and Bacillus subtilis. MenD, an enzyme catalyzing the first irreversible step in the menaquinone biosynthetic pathway has been implicated in the SCV phenotype of S. aureus. In the present work we explored biochemical and structural properties of this important enzyme.
Our structural analysis revealed that despite its low sequence identity of 28%, the overall fold of staphylococcal MenD (saMenD) is similar to Escherichia coli MenD (ecMenD) albeit with some significant disparities. Major structural differences can be observed near the active site region of the protein and are profound in the C-terminal helix and a loop near the active site. The loop contains critical residues for cofactor binding and is well ordered only in the ecMenD-ThDP structure, while in the apo and substrate bound structures of ecMenD the loop is primarily disordered. In our saMenD structure the loop is for the first time completely ordered in the apo form and displays a novel conformation of the cofactor-binding loop. The loop adopts an unusual open conformation and the conserved residues, which are responsible for cofactor binding are located too far away to form a productive complex with the cofactor in this conformation. Additionally, biochemical studies in conjugation with the structural data aided in the identification of the substrate-binding pocket and delineated residues contributing to its binding and catalysis. Thus the present work successfully divulged the unique biochemical and structural characteristics of saMenD.
With 9.6 million new cases and 1.5 million deaths in 2014, tuberculosis (TB) is alongside with AIDS the most deadly infection. Foremost, the increased prevalence of resistant strains of M. tuberculosis among the TB-infected population represents a serious thread. Hence, in the last decades, novel drug targets have been investigated worldwide. So far a relatively unexplored target is the cell wall enzyme β-ketoacyl-ACP-synthase “KasA”, which plays a crucial role in maintaining the membrane impermeability and hence the cell ability to resist to the immune response and drug therapy. KasA is a key enzyme in the fatty acid synthase “FAS-II” elongation cycle, responsible for the extension of the growing acyl chain within the biosynthesis of precursors for the most hydrophobic constituents of the cell wall – mycolic acids. Design of the novel KasA inhibitors, performed in the research group of Prof. Sotriffer by C. Topf and B. Schaefer, was based on the recently published crystal structure of KasA in complex with its known inhibitor thiolactomycin (TLM). Considering the essential ligand-enzyme interactions, a pharmacophore model was built and applied in the virtual screening of a modified ZINC database. Selected hits with the best in silico affinity data have been reported by Topf and Schaefer.
In this work, two of the obtained hits were synthesized and their structure was systematically varied. First, a virtual screening hit, chromone-2-carboxamide derivative GS-71, was modified in the amide part. Since the most of the products possessed a very low solubility in the aqueous buffer medium used in biological assays, polar groups (nitro, succinamidyl and trimethyl-amino substituent in position 6 of the chromone ring or hydroxyl group on the benzene ring in the amide part have been inserted to the molecule. Further variations yielded diaryl ketones, diaryl ketone bearing a succinamidyl substituent, carboxamide bearing a methylpiperazinyl-4-oxobutanamido group and methyl-malonyl ester amides. Basically, the essential structural features necessary for the ligand-enzyme interactions have been maintained. The latter virtual screening hit, a pyrimidinone derivative VS-8 was synthesized and the structure was modified by substitution in positions 2, 4, 5 and 6 of the pyrimidine ring. Due to autofluorescence, detected in most of the products, this model structure was not further varied.
Simultaneously, experiments on solubilization of the first chromone-2-carboxamides with cyclodextrins, cyclic oligosacharides known to form water-soluble inclusion complexes, were performed. Although the assessed solubility of the chromone 3b/DIMEB (1:3) mixture exceeded 14-fold the intrinsic one, the achieved 100 µM solubility was still not sufficient to be used as a stock solution in the binding assay. The experiments with cyclodextrin in combination with DMSO were ineffective. Owing to high material costs necessary for the appropriate cyclodextrin amounts, the aim focused on structural modification of the hydrophobic products.
Precise structural data have been obtained from the solved crystal structures of three chromone derivatives: the screening hit GS-71 (3b), its trimethylammonium salt (18) and 6-nitro-substituted N-benzyl-N-methyl-chromone-2-carboxamide (9i). The first two compounds are nearly planar with an anti-/trans-rotamer configuration. In the latter structure, the carboxamide bridge is bent out of the chromone plane, showing an anti-rotamer, too. Considering the relatively low partition coefficient of compound 3b (cLogP = 2.32), the compound planarity and correlating tight molecular packing might be the factors significantly affecting its poor solubility.
Regarding the biological results of the chromone-based compounds, similar structure-activity correlations could be drawn from the binding assay and the whole cell activity testing on M. tuberculosis. In both cases, the introduction of a nitro group to position 6 of the chromone ring and the presence of a flexible substituent in the amide part showed a positive effect. In the binding study, the nitro group at position 4 on the N-benzyl residue was of advantage, too. The highest enzyme affinity was observed for N-(4-nitrobenzyl)-chromone-2-carboxamide 4c (KD = 34 µM), 6-nitro substituted N-benzyl-chromone-2-carboxamide 9g (KD = 40 µM) and 6‑nitro-substituted N-(4-nitrobenzyl)-chromone-2-carboxamide 9j (KD = 31 µM), which could not be attributed to the fluorescence quenching potential of the nitro group. The assay interference potential of chromones, due to a covalent binding on the enzyme sulfhydryl groups, was found to be negligible at the assay conditions. Moderate in vivo activity was detected for 6‑nitro-substituted N-benzyl-chromone-2-carboxamide 9g and its N-benzyl-N-methyl-, N‑furylmethyl-, N-cyclohexyl- and N-cyclohexylmethyl derivatives 9i, 9d, 9e, 9f, for which MIC values 20 – 40 µM were assessed. Cytotoxicity was increased in the N‑cyclohexylmethyl derivative only. None of the pyrimidine-based compounds showed activity in vivo. The affinity of the model structure, VS-8, surpassed with KD = 97 µM the assessed affinity of TLM (KD = 142 µM).
Since for the model chromone compound GS-71 no reliable KasA binding data could be obtained, a newly synthesized chromone derivative 9i was docked into the KasA binding site, in order to derive correlation between the in silico and in vitro assessed affinity. For the 6‑nitro-derivative 9i a moderate in vivo activity on M. tuberculosis was obtained. The in silico predicted pKi values for TLM and 9i were higher than the corresponding in vitro results, maintaining though a similar tendency, i.e., the both affinity values for compound 9i (pKi predicted = 6.64, pKD experimental = 4.02) surpassed those obtained for TLM (pKi predicted = 5.27, pKD experimental = 3.84). Nevertheless, the experimental pKD values are considered preliminary results.
The binding assay method has been improved in order to acquire more accurate data. Owing to the method development, limited enzyme batches and solubility issues, only selected compounds could be evaluated. The best hits, together with the compounds active on the whole cells of M. tuberculosis, will be submitted to the kinetic enzyme assay, in order to confirm the TLM-like binding mechanism. Regarding the in vivo testing results, no correlations could be drawn between the predicted membrane permeability values and the experimental data, as for the most active compounds 9e and 9f, a very low permeability was anticipated (0.4 and 0.7 %, respectively). Further biological tests would be required to investigate the action- or transport mode.
Macromolecular complexes, also termed molecular machines, facilitate a large spectrum of biological reactions and tasks crucial to the survival of cells. These complexes are composed of either protein only, or proteins bound to nucleic acids (DNA or RNA). Prominent examples for each class are the proteosome, the nucleosome and the ribosome. How such units are assembled within the context of a living cell is a central question in molecular biology. Earlier studies had indicated that even very large complexes such as ribosomes could be reconstituted from purified constituents in vitro. The structural information required for the formation of macromolecular complexes, hence, lies within the subunits itself and, thus, allow for self- assembly. However, increasing evidence suggests that in vivo many macromolecular complexes do not form spontaneously but require assisting factors (“assembly chaperones”) for their maturation. In this thesis the assembly of RNA-protein (RNP) complexes has been studied by a combination of biochemical and structural approaches. A resourceful model system to study this process is the biogenesis pathway of the uridine-rich small nuclear ribonucleoproteins (U snRNPs) of the spliceosome. This molecular machine catalyzes pre-mRNA splicing, i.e. the removal of non-coding introns and the joining of coding exons to functional mRNA. The composition and architecture of U snRNPs is well defined, also, the nucleo-cytoplasmic transport events enabling the formation of these particles in vivo have been analyzed in some detail. Furthermore, recent studies suggest that the formation of U snRNPs in vivo is mediated by an elaborate assembly machinery consisting of protein arginine methyltransferase (PRMT5)- and survival motor neuron (SMN)-complexes. The elucidation of the reaction mechanism of cellular U snRNP assembly would serve as a paradigm for our understanding of how RNA-protein complexes are formed in the cellular environment. The following key findings were obtained as part of this study: 1) Efforts were made to establish a full inventory of the subunits of the SMN-complex. This was achieved by the biochemical definition and characterization of an atypical component of this complex, the unrip protein. This protein is associated with the SMN-complex exclusively in the cytoplasm and influences its subcellular localization. 2) With a full inventory of the components in hand, the architecture of the SMN-complex was defined on the basis of an interaction map of all subunits. This study elucidated that the proteins SMN, Gemin7 and Gemin8 form a backbone, onto which the remaining subunits adhere in a modular manner. 3) The two studies mentioned above formed the basis to elucidate the reaction mechanism of cellular U snRNP assembly. Initially, an early phase in the SMN-assisted formation of U snRNPs was analyzed. Two subunits of the U7 snRNP (LSm10 and 11) were found to interact with the PRMT5-complex, without being methylated. This report suggests that the stimulatory role of the PRMT5-complex is independent of its methylation activity. 4) Key reaction intermediates in U snRNP assembly were found and characterized by a combination of biochemistry and structural studies. Initially, a precursor to U snRNPs with a sedimentation coefficient of 6S is formed by the pICln subunit of the PRMT5-complex and Sm proteins. This intermediate was shown to constitute a kinetic trap in the U snRNP assembly reaction. Progression towards the assembled U snRNP depends on the activity of the SMN-complex, which acts as a catalyst. The formation of U snRNPs is shown to be structurally similar to the way clamps are deposited onto DNA to tether poorly processive polymerases. 5) The human SMN-complex is composed of several subunits. However, it is unknown whether all subunits of this entity are essential for U snRNP assembly. A combination of bioinformatics and biochemistry was applied to tackle this question. By mining databases containing whole-genome assemblies, the SMN-Gemin2 heterodimer is recognized as the most ancestral form of the SMN-complex. Biochemical purification of the Drosophila melanogaster SMN-complex reveals that this complex is composed of the same two subunits. Furthermore, evidence is provided that the SMN-Gemin2 heterodimer is necessary and sufficient to promote faithful U snRNP assembly. Future studies will adress further details in the reaction mechanism of cellular U snRNP assembly. The results obtained in this thesis suggest that the SMN and Gemin2 subunits are sufficient to promote U snRNP formation. What then is the function of the remaining subunits of the SMN-complex? The reconstitution schemes established in this thesis will be instrumental to address this question. Furthermore, additional mechanistic insights into the U snRNP assembly reaction will require the elucidation of structures of the assembly machinery trapped at various states. The prerequisite for these structural studies, the capability to generate homogenous complexes in sufficient amounts, has been accomplished in this thesis.
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.