@phdthesis{Cappel2011, author = {Cappel, Daniel}, title = {Computersimulationen zur Untersuchung von Wassermolek{\"u}len in Protein-Ligand Komplexen am Beispiel einer Modellbindetasche}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-55003}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2011}, abstract = {Wassermolek{\"u}le spielen oft eine entscheidende Rolle bei der Bindung von Liganden an Proteine. Zum einen ist dies in ihrer Eigenschaft als Wasserstoffbr{\"u}ckendonor und -akzeptor begr{\"u}ndet, die es erm{\"o}glicht Wechselwirkung zwischen Ligand und Rezeptor zu vermitteln. Zum anderen stellen die Desolvatisierungsenthalpie und -entropie einer Bindetasche w{\"a}hrend der Ligandbindung einen entscheidenden Anteil der Bindungsaffinit{\"a}t dar. Obwohl man sich dieser Einfl{\"u}sse seit langem bewusst ist, sind aktuelle Methoden des computerbasierten Wirkstoffdesigns nur in sehr begrenztem Umfang in der Lage, die entsprechenden Effekte zu erfassen und vorherzusagen. Da experimentelle Daten {\"u}ber die Effekte von Wassermolek{\"u}len in Protein-Ligand Komplexen von Natur aus schwierig zu erhalten sind, untersucht die vorliegende Arbeit eine Modellbindetasche einer Cytochrom c Peroxidase Mutante (CCP W191G) mit Hilfe von Molecular Modeling Techniken. Diese polare und solvatisierte Kavit{\"a}t ist strukturell sehr gut charakterisiert und bindet kleine, kationische Heterozyklen zusammen mit unterschiedlichen Mengen an Wassermolek{\"u}len. F{\"u}r die Untersuchungen wurden strukturell {\"a}hnliche Liganden mit einem unterschiedlichen Wechselwirkungsmuster ausgew{\"a}hlt. Davon ausgehend wurde die M{\"o}glichkeit zweier Docking-Programme, den Grad der Wasserverdr{\"a}ngung durch den Liganden zusammen mit dem Bindungsmodus vorherzusagen, untersucht. Die dynamischen Eigenschaften der Bindetaschenwassermolek{\"u}le wurden mittels Molekulardynamiksimulationen studiert. Schließlich wurden diese rein strukturellen Betrachtungen durch eine energetische/thermodynamische Analyse komplettiert. Die Anwendung dieser unterschiedlichen Verfahren liefert einige neue Erkenntnisse {\"u}ber die untersuchte Modellbindetasche. Trotz der relativen Einfachheit der kleinen Kavit{\"a}t der CCP W191G Mutante war die vollst{\"a}ndige Charakterisierung und eine korrekte (retrospektive) Vorhersage des Wasser-Wechselwirkungsmuster der Ligand-Komplexe nicht trivial. Zusammenfassend kann man festhalten, dass insgesamt eine gute {\"U}bereinstimmung zwischen den durch Computersimulationen erhaltenen Ergebnissen und den kristallographischen Daten erzielt wurde. Unerwartete Befunde, die auf den ersten Blick mit den kristallographischen Beobachtungen nicht {\"u}bereinstimmen, k{\"o}nnen ebenso durch Limitationen in den Kristallstrukturen bedingt sein. Dar{\"u}ber hinaus gaben die Ergebnisse auch eine Hilfestellung, welches Verfahren zur Beantwortung einer Fragestellung im Rahmen von Wassermolek{\"u}len im Wirkstoffdesign geeignet sind. Schließlich wurden ebenso die Begrenzungen der jeweiligen Methoden aufgezeigt.}, subject = {Cytochromperoxidase}, language = {de} } @phdthesis{Dietl2004, author = {Dietl, Christian}, title = {Beobachtung und Steuerung molekularer Dynamik mit Femtosekunden-Laserpulsen}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-12182}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2004}, abstract = {In dieser Arbeit wurden zwei Aspekte der Femtochemie mit den Methoden der Femtosekunden--Laserspektroskopie untersucht. Dabei wurden folgende Ziele verfolgt: Einerseits sollte die j{\"u}ngst entwickelte Technik der adaptiven Pulsformung auf das Problem bindungsselektiver Photodissoziationsreaktionen angewandt werden, zum Anderen bestand die Aufgabe darin, die nichtadiabatische, photoinduzierte Dynamik am Beispiel der Photoisomerisierung von Stilben mit Hilfe der Photoelektronenspektroskopie zeitaufgel{\"o}st zu untersuchen. Die Methode der adaptiven Pulsformung wurde mit dem Ziel eingesetzt, eine bindungsselektive Photodissoziation zu verwirklichen. Dazu wurde diese Technik in Verbindung mit einem massenspektroskopischen Nachweis der Photofragmente verwendet. Die Experimente wurden an einigen Spezies der Methylhalogenide CH2XY (X,Y = Halogen) durchgef{\"u}hrt. Diese Verbindungen wurden als Modellsysteme gew{\"a}hlt, da sich gezeigt hat, dass auf Grund stark gekoppelter konkurrierender Dissoziationskan{\"a}le durch modenselektive Laseranregung keine Kontrolle erreicht werden kann. Mit dem hier durchgef{\"u}hrten Experiment an CH2ClBr wurde erfolgreich erstmals die Anwendung der adaptiven Femtosekunden-Pulsformung auf das Problem einer bindungsselektiven Photodissoziation demonstriert. Dabei konnte eine Steigerung der Dissoziation der st{\"a}rkeren Kohlenstoff-Halogen Bindung um einen Faktor zwei erreicht werden. Weiterhin konnte experimentell gezeigt werden, dass das optimierte Produktverh{\"a}ltnis nicht durch eine einfache Variation der Laserpulsdauer oder Laserpulsenergie erzielt werden kann. Es wurde ein m{\"o}glicher Mechanismus f{\"u}r die Kontrolle diskutiert, der im Gegensatz zu einem unmodulierten Laserpuls die Wellenpaketdynamik auf neutralen dissoziativen Potentialfl{\"a}chen zur Steuerung des Produktverh{\"a}ltnisses involviert. Wie sich aus einer genaueren Analyse des Fragmentspektrums ergab, wird durch den optimalen Laserpuls die Dissoziation in komplexer Weise moduliert. Dies zeigte sich z.B. auch durch eine {\"A}nderung des Isotopenverh{\"a}ltnisses in der Ausbeute des dissoziierten Br-Liganden vor und nach der Optimierung. Dieser Frage nach einer isotopenselektiven Photodissoziation wurde in einem weiteren Experiment an CH2Br2 nachgegangen. Dabei konnte jedoch nur eine geringe Variation von etwa f{\"u}nf Prozent gegen{\"u}ber dem nat{\"u}rlichen Isotopenverh{\"a}ltnis festgestellt werden. Als gr{\"o}ßtes experimentelles Problem stellte sich dabei die starke Intensit{\"a}tsabh{\"a}ngigkeit der Produktausbeuten heraus, was die Suche nach der optimalen Pulsform stark einschr{\"a}nkte. Anhand des molekularen Photodetachments CH2I2-->CH2+I2 wurde gezeigt, dass durch die Analyse der optimalen Pulsformen Informationen {\"u}ber die Dynamik dieses Prozesses gewonnen werden k{\"o}nnen. Dazu wurde zun{\"a}chst in einem Pump-Probe-Experiment die Dynamik der I2-Fragmentation nach einer Mehrphotonen-Anregung von CH2I2 mit 266nm Laserpulsen untersucht. Dieses Experiment ergab, dass das Molek{\"u}l {\"u}ber einen angeregten Zwischenzustand auf einer sehr schnellen Zeitskala {\"u}ber Dissoziationskan{\"a}len zerfallen kann. Der dominante Kanal f{\"u}hrt zu einer sequentiellen Abgabe einer der I-Liganden und resultiert in den Photoprodukten CH2I und I Im anderen Kanal, dem molekularen Photodetachment, werden die Photoprodukte I2 und CH2 gebildet. In einem Kontrollexperiment wurde dann versucht, das molekulare Photodetachment gegen{\"u}ber dem dominanten sequentiellen Kanal mit geformten 800nm Laserpulsen zu optimieren. Es wurden Optimierungen mit dem Ziel der Maximierung der Ausbeute an den Photoprodukten I2 und CH2 gegen{\"u}ber CH2I durchgef{\"u}hrt. Diese Experimente ergaben, dass f{\"u}r beide Fragmente des molekularen Photodetachments eine Steigerung des Produktverh{\"a}ltnisses um etwa einen Faktor drei m{\"o}glich ist. Dabei zeigte sich, dass eine Maximierung auf ein Produktverh{\"a}ltnis (z.B. I2/CH2I) eine Steigerung des anderen um etwa den gleichen Faktor hervorruft. Dies ist ein deutlicher Hinweis, dass beide Photoprodukte {\"u}ber denselben Dissoziationskanal gebildet werden. Ein weiterer inweis wurde aus der Analyse der optimalen Pulsformen erhalten: In beiden F{\"a}llen weisen diese eine markante Doppelpulsstruktur mit einem zeitlichen Abstand von etwa 400fs auf. Dies erinnert stark an die Situation des Pump-Probe--Experiments, wo durch die Analyse des transienten Signals ebenfalls eine optimale Verz{\"o}gerungszeit zwischen dem Pump- und Probe-Laserpuls von etwa 400fs ermittelt werden konnte, bei der die Produktverh{\"a}ltnisse gerade maximal sind. Im Vergleich zur Massenspektroskopie liefert die Photoelektronenspektroskopie in der kinetischen Energie der Photoelektronen eine zus{\"a}tzliche Messgr{\"o}ße, die direkt Informationen {\"u}ber die Kerngeometrie des Systems liefern kann. Mit dieser Technik wurde die trans-cis-Photoisomerisierung von Stilben im ersten elektronisch angeregten Zustand S1(1Bu) zeitaufgel{\"o}st untersucht. Dabei ging es speziell um die Frage nach der Existenz eines weiteren 1Bu Zustandes, der in neueren theoretischen Untersuchungen diskutiert wurde. In einem Pump-Probe-Experiment wurde dazu das im Molekularstrahl pr{\"a}parierte trans-Stilben durch einen 266nm Laserpuls angeregt und die Dynamik durch einen weiteren 266nm Laserpuls abgefragt. Im Photoelektronenspektrum konnten zwei signifikante Beitr{\"a}ge mit unterschiedlicher Dynamik gefunden werden. Das transiente Signal des ersten Beitrags weist eine Zeitkonstante von etwa 20ps auf und konnte eindeutig der Isomerisierung des S1 Zustandes zugeordnet werden. Im Gegensatz dazu zeigte das Signal des zweiten Beitrags eine Zeitkonstante von 100fs. Dieses Signal k{\"o}nnte aus der Ionisation des S2 Zustandes resultieren, welcher bislang experimentell nicht beobachtet werden konnte.}, subject = {Ultrakurzer Lichtimpuls}, language = {de} } @phdthesis{Falge2012, author = {Falge, Mirjam}, title = {Dynamik gekoppelter Elektronen-Kern-Systeme in Laserfeldern}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-72889}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2012}, abstract = {Die vorliegende Arbeit besch{\"a}ftigt sich mit der theoretischen Untersuchung zweier Themenkomplexe: der Erzeugung Hoher Harmonischer in Molek{\"u}len und dem Einfluss von gekoppelter Elektronen-Kern-Dynamik auf Ultrakurzpuls-Ionisationsprozesse und Quantenkontrolle. W{\"a}hrend bei der Untersuchung der Hohen Harmonischen die Auswirkungen der Kernbewegung auf die Spektren im Mittelpunkt des Interesses stehen, wird bei der Analyse der gekoppelter Elektronen-Kern-Dynamik das Hauptaugenmerk auf die nicht-adiabatischen Effekte gerichtet, die auftreten, wenn Kern- und Elektronenbewegung sich nicht, wie es im Rahmen der Born-Oppenheimer-N{\"a}herung in der Quantenchemie h{\"a}ufig angenommen wird, voneinander trennen lassen.}, subject = {Nichtadiabatischer Prozess}, language = {de} } @phdthesis{Gerbich2015, author = {Gerbich, Thiemo M. P.}, title = {Pikosekunden-zeitaufgel{\"o}ste Deaktivierungsprozesse in isolierten Molek{\"u}len - Fluorenon, NDCA, Me-NI und NTCDA}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-118654}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2015}, abstract = {Im Rahmen der vorliegenden Dissertation wurden die Dynamiken von strahlungslosen Deaktivierungsprozessen von vier verschiedenen Molek{\"u}len im elektronisch angeregten Zustand untersucht. Ein fundiertes Verst{\"a}ndnis der intramolekularen Energieumverteilung in isolierten pi-konjugierten Systemen ist neben dem Modellcharakter auch f{\"u}r Anwendungen in der organischen Elektronik von Interesse. Die Untersuchungen dienen zudem als optimaler Maßstab f{\"u}r theoretische Simulationen, die auf eine Nachbildung der molekularen Dynamik ausgerichtet sind. Die Inbetriebnahme des Pikosekunden-Lasersystems stellt in der Arbeitsgruppe ein großes Potential f{\"u}r die Untersuchung der Dynamik von isolierten pi-konjugierten Molek{\"u}len zur Verf{\"u}gung. Erste Experimente konnten an unterschiedlichen Heterocyclen mit interessantem zeitlichen Verhalten erfolgreich durchgef{\"u}hrt werden und lieferten bereits wichtige Erkenntnisse {\"u}ber die strahlungslose Deaktivierung auf der ps-Zeitskala. Selbst f{\"u}r große Molek{\"u}le mit geringem Dampfdruck, die nur mit hohem experimentellen Aufwand im isolierten Zustand charakterisierbar sind, konnten Relaxationszeiten der angeregten Zust{\"a}nde ermittelt werden. Der Fokus der einzelnen Studien lag in der Erforschung der isolierten Molek{\"u}le, welche durch Anwendung der Molekularstrahl-Technik mit zeitaufgel{\"o}ster REMPI-Spektroskopie anhand des ps-Systems untersucht werden sollten. Zur Kontrolle der experimentellen Ergebnisse wurden zudem Vergleichsmessungen der transienten Absorptionsspektroskopie (TA) in der Fl{\"u}ssigphase herangezogen, wodurch eine fundierte Interpretation der Dynamik m{\"o}glich wurde. Zu den wichtigen Zielen geh{\"o}rten jedoch die Vergleiche der experimentellen Ergebnisse von isolierten Molek{\"u}len mit Berechnungen der Zustandsenergien sowie Simulationen der Molek{\"u}ldynamik aus dem Theorie-Arbeitskreis von Prof. Mitric. Auf diese Weise konnten wichtige Erkenntnisse {\"u}ber die Dynamik der Deaktivierungsprozesse gewonnen werden. Die Kombination der Gasphasen-Experimente mit TA-Messungen in der Fl{\"u}ssigphase hat sich als besonders n{\"u}tzlich erwiesen, um bei mehrstufigen Deaktivierungsprozessen einen erweiterten Einblick in die Dynamik der Molek{\"u}le zu erhalten. - So konnte bei Fluorenon in Cyclohexan und Acetonitril durch Vergleich der Anregungen des S3- und S1-Zustands eine zus{\"a}tzliche Zeitkonstante von 8-16 ps beobachtet werden, welche die innere Umwandlung zum S1-Zustand dokumentiert und die Ergebnisse der Gasphasen-Messungen best{\"a}tigt. - Durch Verwendung von L{\"o}sungsmitteln unterschiedlicher Polarit{\"a}t und der damit verbundenen Verschiebung der elektronischen Zust{\"a}nde von Fluorenon konnte zudem der zweite Deaktivierungsprozess eindeutig einem ISC-Prozess mit Zeitkonstanten von 120-154 ps zugeordnet werden. In der Gasphase wurde dieser Prozess lediglich als langlebiger Offset wahrgenommen. - Unterschiedliche Anregungsenergien zeigten bei TA-Messungen von NDCA eine nahezu identische Molek{\"u}ldynamik mit ca. 200 ps, w{\"a}hrend f{\"u}r isoliertes NDCA ein starker Abfall der Lebensdauer mit zunehmender Schwingungsenergie beobachtet wurde. In der Gasphase wird somit von einer Deaktivierung {\"u}ber eine Energiebarriere ausgegangen, w{\"a}hrend in L{\"o}sung eine zu schnelle Abk{\"u}hlung durch Schwingungsrelaxation diesen Prozess verhindert. - Bei NTCDA konnten in den TA-Messungen nach Anregung des S1-Zustands eine Relaxation in die Triplett-Umgebung innerhalb von wenigen Pikosekunden beobachtet werden, was im Einklang mit der sehr schnellen Deaktivierung in der Gasphase betrachtet werden kann. Eine ausf{\"u}hrliche Vergleichsstudie von isolierten Molek{\"u}len mit computergest{\"u}tzten Rechnungen und Simulationen wurde f{\"u}r die Molek{\"u}le NDCA und Me-NI durchgef{\"u}hrt. Dabei wurde explizit auf den Einfluss von Spin-Bahn-Kopplungen und konischen Durchschneidungen eingegangen, welche zu konkurrierenden Deaktivierungsprozessen des S1-Zustands f{\"u}hren k{\"o}nnen. - Durch Simulationen der Surface-Hopping-Dynamik wurde deutlich, dass bei NDCA und Me-NI im ersten angeregten Zustand eine konische Durchschneidung (CI) zwischen dem S1- und S0-Zustand erreicht werden kann. - W{\"a}hrend die Dynamik von NDCA bei h{\"o}herer Schwingungsanregung stark durch die CI dominiert wird, spielt die direkte Relaxation in den elektronischen Grundzustand bei Me-NI offenbar keine Rolle. - In Abwesenheit der CI zeigen NDCA und Me-NI in einer mit Spin-Bahn-Kopplung erweiterten Simulation der Populationsdynamik einen signifikanten Populationstransfer in die Triplett-Umgebung (T1-T4). Eine innere Umwandlung in den Grundzustand konnte jedoch nur bei Erreichen der CI beobachtet werden. Eine weitere Verbesserung der ps-Experimente kann durch Aufbau eines Photoelektronen-Spektrometers erreicht werden, da durch diese Technik eine pr{\"a}zisere Aussage dar{\"u}ber getroffen werden kann, aus welchem elektronischen Zustand die Molek{\"u}le ionisiert wurden. Eine Unterscheidung von ISC- und IC-Prozessen k{\"o}nnte somit gew{\"a}hrleistet werden.}, subject = {REMPI}, language = {de} } @phdthesis{Kehrein2022, author = {Kehrein, Josef}, title = {Simulationsstudien zur ortsspezifischen Biokonjugation maßgeschneiderter Polymere}, doi = {10.25972/OPUS-28958}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-289589}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2022}, abstract = {Polymer-Biokonjugationen, vornehmlich mit dem Goldstandard PEG, f{\"u}hren zu einer verbesserten Pharmakokinetik, beeinflussen aber auch die konformative Stabilit{\"a}t von Proteinen. Bisherige Mutationsstudien, in denen {\"u}berwiegend (Asn)PEG4 -Konjugate der Beta-faltblattstrukturreichen, humanen Pin 1 WW-Dom{\"a}ne untersucht wurden, postulieren auf einer Proteindesolvatation beruhende Stabilisierungsmechanismen: eine St{\"a}rkung intramolekularer Salzbr{\"u}cken und NH-pi-Bindungen, sowie entropisch g{\"u}nstige Wasserverdr{\"a}ngungen um apolare Aminos{\"a}uren und Hydroxylgruppen. Ziel dieser Arbeit ist es, die Protein-Polymer-Dynamik auf molekularer Ebene zu charakterisieren, um damit rationale Ans{\"a}tze zum Design neuer Biokonjugate voranzutreiben und m{\"o}gliche PEG-Alternativen zu etablieren. Hierzu wurde eine Vielzahl an Deskriptoren mittels Molekulardynamik-Simulationen der WW-Konjugate gewonnen und mit publizierten Stabilit{\"a}tsdaten in multivariaten Regressions- und logistischen Klassifikationsmodellen korreliert. Die gewonnenen QSPR-Modelle decken im Vergleich zu einer bereits publizierten, kristallstrukturbasierten Richtlinie einen gr{\"o}ßeren und strukturell vielf{\"a}ltigeren Datensatz an Konjugaten ab und zeigen gleichzeitig, auch f{\"u}r ein Konjugat der Src SH3-Dom{\"a}ne, eine deutlich verbesserte Leistung. Die Modelldeskriptoren beschreiben sowohl eine Modulation der Solvatation als auch Protein-Polymer-Interaktionen. Metadynamik-Simulationen zeigten zudem die Polymerdynamik w{\"a}hrend einer partiellen Proteinentfaltung auf. Mithilfe weiterer Simulationen von Konjugaten des alpha-helikalen Her2-Affibodys wurde die Dynamik von PEG und verschiedener Alternativen (LPG, PEtOx, PMeOx) systematisch studiert. PEG interagierte mit positiv geladenen Lysinen und Argininen in der N{\"a}he hydrophober Aminos{\"a}uren. LPG zeigte zus{\"a}tzliche Wechselwirkungen der Hydroxylgruppen mit Aspartaten und Glutamaten. POx-Polymere interagierten mit Phenylalaninen, Tyrosinen und {\"u}ber Carbonylgruppen mit HB-Donatoren. Gr{\"o}ßere Konjugate (10 - 50 kDa PEG/LPG/PEtOx) des antiviralen Biologikums Interferon-alpha2a wurden mittels gaußbeschleunigter MDs und einer CG-Simulation analysiert. Charakteristische Wechselwirkungspartner stimmten mit den Beobachtungen zu Oligomer-Konjugaten {\"u}berein. In Einklang mit experimentellen Daten der Kooperationspartner zu den 10-kDa-Varianten deuteten zus{\"a}tzliche Constrained-Network-Analysen, welche die Proteinflexibilit{\"a}t evaluieren, auf eine thermische Destabilisierung hin. Die Bioaktivit{\"a}t der untersuchten Konjugate wurde weiterhin erfolgreich mit den Gyrationsdurchmessern der modellierten Strukturen korreliert.}, subject = {Konjugate}, language = {de} } @phdthesis{Nocker2012, author = {Nocker, Monika}, title = {Molekulardynamische Untersuchungen zur Charakterisierung von Flexibilit{\"a}t, Bindemechanismen und Bindungsaffinit{\"a}ten von Aldose Reduktase und Nukle{\"a}ren Rezeptoren}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-72110}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2012}, abstract = {Aldose Reduktase ALR2 katalysiert den ersten Schritt des Sorbitol-Stoffwechselweges. In diesem wird mit Hilfe des Kofaktors NADPH Glukose zu Sorbitol reduziert. Bei erh{\"o}htem Blutzuckerspiegel, wie dies bei Diabetes-Patienten der Fall ist, ist dieser metabolische Weg von Bedeutung. Bis zu einem Drittel der Blutglukose wird zu Sorbitol reduziert. Die Folge der Sorbitolakkumulation in den Zellen und der Verminderung der NADPH-Konzentration sind „osmotischer" sowie „oxidativer" Stress. Diese stehen in Zusammenhang mit den vielfach diskutierten Sp{\"a}tsch{\"a}den des Diabetes, wie diabetischer Katarakt, Neuro- und Nephropathie. Das Enzym ist experimentell sehr gut untersucht und eignet sich daher als Modellsystem zur Untersuchung der intrinsischen Proteinflexibilit{\"a}t und thermodynamischer Daten mit Hilfe von Computermethoden. Unter diesen Voraussetzungen steht der Gewinn eines besseren Verst{\"a}ndnisses von molekularer Erkennung und Proteinbeweglichkeit der ALR2 unter Verwendung von Molekulardynamik-Simulationen MD als prim{\"a}res Ziel im Zentrum dieser Arbeit. Dabei wurden MD-Studien zu zwei kristallographisch erhaltenen Protein-Ligand-Komplexen durchgef{\"u}hrt. Die Liganden unterscheiden sich nur geringf{\"u}gig in der L{\"a}nge einer Seitenkette, ihre Bindung f{\"u}hrt allerdings zu g{\"a}nzlich unterschiedlichen Bindemodi. Einer davon ist bislang einzigartig f{\"u}r die ALR2. Mit Hilfe von MD-Simulationen wurde versucht, eine Erkl{\"a}rung f{\"u}r diese neue Konformation der Bindetasche im Vergleich zu jener eines strukturell sehr {\"a}hnlichen Liganden zu finden. Außerdem waren {\"u}ber diese Studien Aussagen {\"u}ber besonders flexible Bereiche der ALR2-Bindetasche m{\"o}glich, die mit bereits existierenden Erkenntnissen {\"u}ber die Bindetaschenflexibilit{\"a}t verglichen werden konnten. Dar{\"u}ber hinaus gelang es, durch die Methode der gesteuerten Molekulardynamik SMD einen {\"U}bergang zwischen einer r{\"o}ntgenkristallografisch ermittelten kofaktorgebundenen Holo-Konformation und kofaktorfreien Apo-Konformation zu simulieren. Computergest{\"u}tzte Methoden erm{\"o}glichen es also, weitl{\"a}ufige Bewegungen von einer Proteinkonformation in die andere nachzuvollziehen bzw. die experimentell erhaltenen Strukturen zu best{\"a}tigen. Eine mechanistische Deutung des Kofaktorassoziations- und Kofaktordissoziationsprozesses wurde ebenfalls versucht. Daf{\"u}r war es notwendig, strukturelle Ver{\"a}nderungen im Protein zeitlich zu verfolgen und entscheidende Vorg{\"a}nge zu identifizieren. Die Methode der SMD wurde in dieser Arbeit auch auf ein weiteres, pharmakologisch interessantes System {\"u}bertragen. Dabei wurde versucht auch an zwei Vertretern der Klasse der Nukle{\"a}ren Rezeptoren NRs, dem Androgenrezeptor AR und dem Estrogenrezeptor ER, eine solche weitreichende Bewegung nachzuvollziehen. Auch bei diesen Rezeptoren sind zwei in der Position einer alpha-Helix unterschiedliche Formen bekannt. Auch hier wurden mit Hilfe der genannten Methode, relevante Ereignisse hinsichtlich der Helixmobilit{\"a}t identifiziert. Abschließend wurde auf den thermodynamischen Aspekt der Protein-Ligand-Komplexe eingegangen. Durch Berechnungen anhand der Methode der thermodynamischen Integration TI wurden relative Bindungsaffinit{\"a}ten am Modellsystem ALR2 gewonnen. Durch den Vergleich mit experimentell vorhandenen Daten konnte die Methode validiert werden. Das Verfahren der TI sollte in Zukunft eine Voraussage von Affinit{\"a}ten beliebiger, sich geringf{\"u}gig unterscheidender Inhibitoren, die aber denselben Bindemodus aufweisen, erm{\"o}glichen und damit den Prozess des Wirkstoffdesigns erleichtern. Zusammenfassend ergab sich eine gute {\"U}bereinstimmung der experimentell ermittelten Strukturen bzw. Daten mit den durch Computersimulationen erhaltenen.}, subject = {Molekulardynamik}, language = {de} } @phdthesis{Raffauf2011, author = {Raffauf, Claudia}, title = {Thermodynamische Untersuchungen und quantenchemische Berechnungen zum L{\"o}sungsverhalten von Phenylethylaminen}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-56239}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2011}, abstract = {Die Wasserl{\"o}slichkeit von Wirkstoffen ist von elementarer Bedeutung f{\"u}r deren pharmazeutische Verwendung. Im Rahmen dieser Arbeit wurde der Ansatz der thermodynamischen Betrachtung des L{\"o}sungsprozesses von einfachen Modellsubstanzen auf pharmazeutisch relevante Substanzen aus der Gruppe der Phenylethylamine {\"u}bertragen und entsprechend erweitert. Der L{\"o}sungsprozess wurde auf Grundlage des Satzes von Hess in Ersatzschritte aufgeteilt, die dann mittels kalorimetrischer Bestimmung, quantenchemischen Rechnungen und MD-Simulation bestimmt wurden, um die limitierenden Faktoren der L{\"o}slichkeit zu beschreiben und so gezielte Ans{\"a}tze f{\"u}r die L{\"o}slichkeitsverbesserung finden zu k{\"o}nnen. Der L{\"o}sungsprozess wurde zum einen durch die Schritte Sublimation und Solvatation und zum anderen durch Schmelzen und Mischen beschrieben. F{\"u}r die einzelnen Schritte sind die thermodynamischen Gr{\"o}ßen Enthalpie, Entropie und Freie Enthalpie bestimmt worden. Dies war zum Teil experimentell mithilfe von kalorimetrischen Messungen m{\"o}glich, andere Gr{\"o}ßen sind mit dem semiempirischen Solvatationsprogramm AMSOL oder mit COSMOtherm berechnet worden. Die so gewonnenen thermodynamischen Gr{\"o}ßen geben Aufschluss {\"u}ber das L{\"o}sungsverhalten der untersuchten Substanzen, die zum Teil trotz augenscheinlich hydrophilen Merkmalen schlecht l{\"o}slich in Wasser sind. Zus{\"a}tzlich wurden mit Hilfe von MD-Simulationen mit dem Programm AMBER die inter- und intramolekularen Wechselwirkungsm{\"o}glichkeiten untersucht. Die Ausbildung intermolekularer Wasserstoffbr{\"u}ckenbindungen mit Wasser wurde mit Radialen-Dichteverteilungs-Funktionen (RDF) analysiert. Die M{\"o}glichkeit der Ausbildung von intramolekularen Wasserstoffbr{\"u}ckenbindungen und deren Einfluss auf das L{\"o}sungsverhalten wurde {\"u}ber die Bindungsabst{\"a}nde und der Darstellung der polaren Oberfl{\"a}chen bestimmt.}, subject = {Thermodynamik}, language = {de} } @phdthesis{Schaefer2012, author = {Schaefer, Benjamin}, title = {Computergest{\"u}tzte Untersuchungen zur Inhibition und Dynamik der ß-Ketoacyl-ACP-Synthase I (KasA) aus Mycobacterium tuberculosis}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-74635}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2012}, abstract = {Die vorliegende Arbeit beschreibt die Durchf{\"u}hrung computergest{\"u}tzter Untersuchungen an den Wildtyp-Kristallstrukturen dieses Enzyms - einerseits zur Suche nach neuen Inhibitoren im Rahmen von virtuellen Screening- (VS-) Studien, andererseits zur Charakterisierung der strukturellen Flexibilit{\"a}t mit Hilfe von Molekulardynamik- (MD-) Simulationen. F{\"u}r ein erstes VS wurde zun{\"a}chst eine Datenbank von mehreren Millionen kommerziell erh{\"a}ltlichen Verbindungen mit Arzneistoff-{\"a}hnlichen physikochemischen Eigenschaften erstellt. Als Ausgangspunkt der Screening-Studie diente ein Teildatensatz, welcher mit Hilfe eines auf dem nativen Bindemodus des KasA-Inhibitors Thiolactomycin (TLM) beruhenden Pharmakophormodells erhalten wurde. Diese Verbindungen wurden in die Bindetasche von KasA gedockt und die Qualit{\"a}t der erhaltenen Posen in einem Re- und Konsensus-Scoring-Verfahren bewertet. Schließlich wurden 14 Substanzen k{\"a}uflich erworben und im Rahmen einer Fluoreszenz-Bindungsstudie experimentell getestet. F{\"u}r sechs Molek{\"u}le war gegen{\"u}ber KasA eine schwache, zu TLM vergleichbare Aktivit{\"a}t zu verzeichnen. Eine zweite VS-Studie befasste sich mit der Bewertung der Bindungsaffinit{\"a}t synthetisch leicht zug{\"a}nglicher Derivate von GS95, einem gegen{\"u}ber dem KasA-Wildtyp aktiven Molek{\"u}l mit einer 1-Benzyluracil-Grundstruktur. Anhand geeigneter Synthesebausteine wurde eine virtuelle Datenbank von insgesamt 16 Derivaten erstellt, welche in die Bindetasche des Enzyms gedockt wurden. F{\"u}r die vorhergesagten Bindemodi erfolgte dann eine Absch{\"a}tzung der freien Bindungsenthalpie. Nach einer Bewertung der Orientierungen auf Basis der errechneten ΔG-Werte sowie einer visuellen Analyse wurden schließlich elf Verbindungen synthetisiert und im Fluoreszenz-Experiment getestet, wobei f{\"u}r alle Uracilderivate eine Aktivit{\"a}t zu beobachten war. Die Kd-Werte fallen jedoch {\"a}hnlich hoch aus wie bei GS95. Zur Untersuchung der Strukturdynamik des KasA-Wildtyps wurden drei MD-Simulationen des homodimeren Proteins von je 15 ns L{\"a}nge durchgef{\"u}hrt. Mit Hilfe von 2D-RMSD-Berechnungen und einer hierarchischen Clusteranalyse wurden aus den drei Simulationen insgesamt zehn repr{\"a}sentative Snapshots entnommen, welche die im Rahmen der Simulationszeit von insgesamt 90 ns produzierte strukturelle Vielfalt der Bindetasche von KasA wiedergeben. Wie die Analysen zeigen, wird ein dualer Charakter hinsichtlich der Flexibilit{\"a}t der unmittelbaren Taschenreste beobachtet. Hierbei zeigt Phe404 eine besonders ausgepr{\"a}gte strukturelle Vielfalt; diese Beobachtung deckt sich mit der gatekeeper-Rolle der Aminos{\"a}ure zwischen der Malonyl-Bindetasche und dem Acyl-Bindekanal, welcher f{\"u}r die Unterbringung der wachsenden Fetts{\"a}urekette im Enzym w{\"a}hrend der Katalyse verantwortlich zeichnet. Dar{\"u}ber hinaus erkl{\"a}rt die hohe Flexibilit{\"a}t von Phe404 m{\"o}glicherweise die recht schwache Bindungsaffinit{\"a}t von TLM gegen{\"u}ber dem Wildtyp von KasA, da die gatekeeper-Aminos{\"a}ure nur in der geschlossenen Form einen stabilisierenden Effekt auf den Liganden aus{\"u}bt. Besondere Bedeutung kommt hierbei einem Wassermolek{\"u}l zu, welches als eine Art molekularer Schalter f{\"u}r die Flexibilit{\"a}t von Phe404 betrachtet werden kann und somit die Fixierung von TLM in der Bindetasche maßgeblich beeinflusst. Des Weiteren wurden innerhalb der Tasche hohe Besetzungsraten f{\"u}r je ein Wassermolek{\"u}l identifiziert. Die aus den MD-Simulationen gewonnenen Erkenntnisse wurden anschließend zur Aufstellung von Empfehlungen f{\"u}r das Design neuartiger KasA-Inhibitoren verwendet. Des Weiteren wurde die Dynamik des oben erw{\"a}hnten Acyl-Bindekanals, welcher sich aus den Aminos{\"a}uren 115-147 zusammensetzt, n{\"a}her charakterisiert. Hierbei wurden die Reste 115-119 und insbesondere Leu116 als zweiter gatekeeper identifiziert, welcher die {\"O}ffnung des Acyl-Bindekanals zur Oberfl{\"a}che des Proteins reguliert und somit eine entscheidende Funktion bei der Unterbringung der langkettigen Fetts{\"a}uresubstrate {\"u}bernimmt. Schließlich wurden zwei repr{\"a}sentative Bindetaschenkonformationen aus den MD-Simulationen hinsichtlich einer Verwendung in strukturbasierten VS-Studien n{\"a}her untersucht. Mit Hilfe von hot spot Analysen und unter Ber{\"u}cksichtigung oben genannter Empfehlungen f{\"u}r das Design neuartiger KasA-Inhibitoren wurden verschiedene Pharmakophormodelle erstellt, welche nach Durchsuchung der zu Anfang dieses Kapitels erw{\"a}hnten virtuellen Molek{\"u}ldatenbank zwischen 149 und 420 verschiedene hit-Strukturen lieferten. Folglich scheint eine Adressierung der beiden Konformationen durch arzneistoffartige Verbindungen prinzipiell m{\"o}glich. Unter den erhaltenen Verbindungen herrscht eine hohe strukturelle Vielfalt; außerdem unterscheiden sich diese im Allgemeinen deutlich von den Molek{\"u}len aus den vorangegangenen VS Studien, was das Potential der beiden Bindetaschenkonformationen zur Identifizierung von potentiellen KasA-Inhibitoren mit neuartigen Grundstrukturen zum Ausdruck bringt.}, subject = {Tuberkelbakterium}, language = {de} } @phdthesis{Suess2021, author = {S{\"u}ß, Jasmin}, title = {Theoretische Untersuchungen an molekularen Aggregaten: 2D-Spektroskopie und Exzitonendynamik}, doi = {10.25972/OPUS-24713}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-247136}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2021}, abstract = {Diese Dissertation besch{\"a}ftigt sich mit der Exzitonendynamik molekularer Aggregate, die nach Mehrphotonen-Anregung auf ultrakurzer Zeitskala stattfindet. Hierbei liegt der Fokus auf der Charakterisierung der Exziton-Exziton-Annihilierung (EEA) mithilfe von zweidimensionaler optischer Spektroskopie f{\"u}nfter Ordnung. Dazu werden zwei verschiedene Modellsysteme implementiert: Das elektronische Homodimer und das elektronische Homotrimer-Modell, wobei Letzteres eine Erweiterung des Dimer-Modells darstellt. Die Kopplung des quantenmechanischen Systems an die Umgebung wird mithilfe des Quantum-Jump-Ansatzes umgesetzt. Besonderes Interesse kommt der Analyse des Signals f{\"u}nfter Ordnung in Abh{\"a}ngigkeit der Populationszeit T zu. Anhand des Dimer-Modells als kleinstm{\"o}gliches Aggregat lassen sich bereits gute Vorhersagen auch {\"u}ber das Verhalten gr{\"o}ßerer molekularer Aggregate treffen. Der Zerfall des oszillierenden Signals f{\"u}r lange Populationszeiten korreliert mit der EEA. Dies zeigt, dass die zweidimensionale optische Spektroskopie genutzt werden kann, um den Annihilierungsprozess zu charakterisieren. Innerhalb des Modells des Dimers wird weiterhin der Einfluss der Intraband-Relaxation untersucht. Zunehmende Intraband-Relaxation verhindert den Austausch zwischen den lokalen Zust{\"a}nden, der essentiell f{\"u}r den Annihilierungsprozess ist, und die EEA wird blockiert. Das elektronische Trimer-Modell erweitert das Dimer-Modell um eine Monomereinheit. Somit befinden sich die Exzitonen im Anschluss an die Anregung nicht mehr unvermeidlich nebeneinander. Es gibt somit eine Konfiguration, bei der sich die Exzitonen zun{\"a}chst zueinander bewegen m{\"u}ssen, bevor die Startbedingung des Annihilierungsprozesses gegeben ist. Dieser zus{\"a}tzliche Schritt wird auch Exzitonendiffusion genannt. Die Ergebnisse dieser Arbeit legen nahe, dass das erwartete Verhalten nur zu sehr kurzen Zeiten im Femtosekundenbereich auftritt und somit die Zeitskala der Exzitonendiffusion im Falle des Trimers nicht sichtbar wird. Es bedarf demnach eines gr{\"o}ßeren Modellsystems, bei dem sich der Effekt der zeitverz{\"o}gert eintretenden EEA deutlich in der Zerfallsdynamik manifestieren kann.}, subject = {Molekulardynamik}, language = {de} } @phdthesis{Wehner2018, author = {Wehner, Johannes}, title = {Wellenfunktionsbasierte Analyse zweidimensionaler Spektren: Wellenpaketbewegung in Dimeren und Quantendiffusionsdynamik}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-163555}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2018}, abstract = {Diese Arbeit befasst sich mit der st{\"o}rungstheoretischen Berechnung von zweidimensionalen Photonen-Echo-Spektren f{\"u}r das elektronische und vibronische Modell eines Homo- und Hetero-Dimers sowie f{\"u}r ein vibronisches Modell eines Monomers unter dem Einfluss einer System-Bad-Wechselwirkung. Bei der Analyse der Dimerspektren steht neben der Orientierungsmittelung der Polarisation dritter Ordnung der Unterschied zwischen elektronischen und vibronischen Spektren sowie der Vergleich der Spektren von Homo- und Hetero-Dimeren im Zentrum des Interesses. Bei der Analyse der Monomer-Spektren steht die Behandlung einer dissipativen Dynamik bzw. des vibrational-coolings innerhalb eines stochastischen Ansatzes im Vordergrund. Der erste Teil dieser Arbeit konzentriert sich auf die st{\"o}rungstheoretische Berechnung der Polarisation dritter Ordnung in Dimeren. Dabei werden alle Aspekte und Ergebnisse f{\"u}r verschiedene Geometrien der {\"U}bergangsdipolmomente analysiert und diskutiert. Die Berechnungen ber{\"u}cksichtigen dabei auch die zuf{\"a}llige Anordnung der Molek{\"u}le in der Probe. Die Zusammenh{\"a}nge zwischen den 2D-Spektren und den Eigenschaften der Monomereinheiten, die Abh{\"a}ngigkeit der Intensit{\"a}ten mancher Peaks von der zeitlichen Abfolge der Pulse sowie der Einfluss der elektronischen Kopplung und verschiedener {\"U}bergangsdipolmomente erm{\"o}glichen ein grundlegendes Verst{\"a}ndnis der elektronischen Photonen-Echo-Spektren. Im elektronischen Dimer wird der Hetero-Dimer-Charakter durch verschiedene Monomeranregungsenergien sowie unterschiedliche {\"U}bergangsdipolmomente der Monomereinheiten bestimmt. Der Einfluss dieser Gr{\"o}ßen auf die Photonen-Echo-Spektren kann durch die Kombination einer detaillierten analytischen Betrachtung und numerischen Rechnungen anschaulich nachvollzogen werden. In der vibronischen Betrachtungsweise zeigt sich, dass die Spektren deutlich an Komplexit{\"a}t gewinnen. Durch die Vibrationsfreiheitsgrade vervielfachen sich die m{\"o}glichen {\"U}berg{\"a}nge im System und damit die m{\"o}glichen Peakpositionen im Spektrum. Jeder Peak spaltet in eine Vibrationssubstruktur auf, die je nach ihrer energetischen Position mit anderen {\"u}berlagern kann. Der Vergleich zwischen Homo- und Hetero-Dimer-Spektren wird durch die Wahl verschiedener Vibrationsfrequenzen und unterschiedlicher Gleichgewichtsabst{\"a}nde entlang der Vibrationskoordinaten erweitert. Die Berechnung des Orientierungsmittels erfolgt mit zwei verschiedenen Ans{\"a}tzen. Zum einen wird das Mittel durch den numerischen sampling-Ansatz berechnet. Dabei werden Azimutal- und Polarwinkel in kleinen Winkelinkrementen abgetastet und f{\"u}r jede Kombination ein 2D-Spektrum berechnet. Die Einzelspektren werden anschließend gemittelt. Diese Methode erweist sich im Dimer als sehr effektiv. Zum anderen erlaubt die analytische Auswertung der Polarisation dritter Ordnung, das gemittelte Spektrum direkt in einer einzelnen Rechnung durch winkelgemittelte Gewichtungsfaktoren zu bestimmen. Bei der Berechnung der elektronischen 2D-Spektren ist diese Methode sehr leistungsf{\"a}hig, da alle Ausdr{\"u}cke analytisch bekannt sind. F{\"u}r vibronische Systeme ist dieser Ansatz ebenfalls sehr leistungsstark, ben{\"o}tigt aber eine einmalige aufwendige Analyse vor der Berechnung. Trotz der deutlich erh{\"o}hten Anzahl an Zustandsvektoren, die propagiert werden m{\"u}ssen, ist diese Methode circa zweimal schneller als die direkte Mittelung mit der sampling-Methode. Im zweiten Teil konzentriert sich die Arbeit auf die Beschreibung eines Monomers, das sich in einer dissipativen Umgebung befindet. Dabei wird auf die L{\"o}sung einer stochastischen Schr{\"o}dingergleichung zur{\"u}ckgegriffen. Speziell wird die sogenannte quantum-state-diffusion-Methode benutzt. Dabei werden nicht nur die Erwartungswerte f{\"u}r die Energie und den Ort, sondern auch die Polarisation dritter Ordnung - eine phasensensitive Gr{\"o}ße - bestimmt. In der theoretischen Fragestellung wird dabei, ausgehend von der von-Neumann Gleichung, die Zeitentwicklung der reduzierten Dichtematrix durch die Integration einer stochastischen zeitabh{\"a}ngigen Schr{\"o}dingergleichung reproduziert. In Rechnungen koppelt die Stochastik {\"u}ber die Erwartungswerte von Ort und Impuls die verschiedenen st{\"o}rungstheoretischen Korrekturen der Wellenfunktion miteinander. Die Spektren, die aus den numerischen Simulationen erhalten werden, spiegeln das dissipative Verhalten des Systems detailliert wider. Eine Analyse der Erwartungswerte von Ort und Energie zeigt, dass sich die einzelnen elektronischen Zust{\"a}nde wie ged{\"a}mpfte harmonische Oszillatoren verhalten und jeweils einen exponentiellen Zerfall abh{\"a}ngig von der Dissipationskonstante zeigen. Dieser Teil der Arbeit erweitert vorausgehende Untersuchungen, bei denen ein vereinfachter Ansatz zu Einsatz kam, der die korrelierte Stochastik nicht ber{\"u}cksichtigte.}, subject = {Molekulardynamik}, language = {de} }