Refine
Document Type
- Doctoral Thesis (8)
Has Fulltext
- yes (8)
Is part of the Bibliography
- no (8)
Keywords
- Molekulardynamik (2)
- Pump-Probe-Technik (2)
- Computersimulation (1)
- Dissipation (1)
- Doorway-Window-Formalismus (1)
- Doorway-window-formalism (1)
- Freie Energie (1)
- Hairpins (1)
- Hauptkomponentenanalyse (1)
- J coupling constants (1)
Institute
This work presents a contribution to the literature on methods in search of lowdimensional models that yield insight into the equilibrium and kinetic behavior of peptides and small proteins. A deep understanding of various methods for projecting the sampled configurations of molecular dynamics simulations to obtain a low-dimensional free energy landscape is acquired. Furthermore low-dimensional dynamic models for the conformational dynamics of biomolecules in reduced dimensionality are presented. As exemplary systems, mainly short alanine chains are studied. Due to their size they allow for performing long simulations. They are simple, yet nontrivial systems, as due to their flexibility they are rapidly interconverting conformers. Understanding these polypeptide chains in great detail is of considerable interest for getting insight in the process of protein folding. For example, K. Dill et al. conclude in their review [28] about the protein folding problem that "the once intractable Levinthal puzzle now seems to have a very simple answer: a protein can fold quickly and solve its large global optimization puzzle simply through piecewise solutions of smaller component puzzles".
Seit langem gibt es ein großes Interesse daran, den zeitlichen Ablauf von chemischen Reaktionen zu untersuchen. Ein wichtiges Instrument hierfür ist die Kurzzeitspektroskopie, die in jüngerer Vergangenheit parallel zum rasanten technischen Fortschritt mit immer kürzeren Laserpulsen bis hinein in den Femtosekunden-Bereich einen neuen Boom erlebt hat. Durch die extrem hohe Zeitauflösung ist es heutzutage möglich geworden, die Reaktionsdynamik eines Moleküls quasi in Echtzeit zu verfolgen, weil man durch die Femtosekunden-Pulse die Bewegung der Atomkerne wie mit einer schnellen Kamera beobachten kann. Mit wachsender Anzahl der Experimente gibt es gleichzeitig einen steigenden Bedarf, die Vielzahl der gemessenen und in vielen Aspekten noch unverstandenen Signale auf mikroskopischer Ebene theoretisch zu beschreiben und zu erklären. Besonders interessant dabei ist, dass dies u.a. für komplexe Systeme möglich wird, weil auf den extrem kurzen Zeitskalen auch in sehr großen Molekülen nur einige wenige Freiheitsgrade eine Rolle spielen. Thema dieser Arbeit war die überwiegend klassische Modellierung von nicht-adiabatischer Kurzzeitdynamik in molekularen Quantensystemen, deren Beschreibung über die Born-Oppenheimer-Näherung hinausgeht. Dabei wurden im Rahmen der vorliegenden Arbeit neue Methoden und Rechenverfahren entwickelt, die einerseits eine verbesserte. physikalisch intuitive Anschauung vermitteln und andererseits aufgrund geeigneter Näherungen deutlich Rechenzeit einsparen und dadurch den Zugang zu Modellen größerer Moleküle ermöglichen können. Zunächst wurde die Simulation von zeitaufgelösten Pump-Probe-Spektren mit Hilfe der Franck-Condon-Approximation durchgeführt, bei der die nukleare. Dynamik während des Laserpulses vernachlässigt wird. Diese bekannte Näherung konnte im Rahmen der Arbeit erstmals auf nichtadiabatisch gekoppelte Potentialflächen angewandt werden. Anschließend wurde für das interessierende Quantensystem mit Hilfe des Mapping-Formalismus ein klassisches Analogon eingeführt und dessen Dynamik eingehend studiert. Dies ist deshalb bemerkenswert, weil dadurch erstmals eine Anwendung von klassischen Methoden aus dem Bereich der nichtlinearen Dvnamik auf nichtadiabatische Quantensvsteme möglich wurde. Schließlich konnten klassische periodische Bahnen des Systems identifiziert werden. Dabei handelte es sich im hier betrachteten Fall nichtadiabatisch gekoppelter Potentiale um eine völlig neue Art von vibronischen periodischen Orbits, die sowohl aus einem nuklearen als auch einem elektronischen Freiheitsgrad bestehen. Dadurch konnte in einem nächsten Schritt mit Hilfe einiger weniger Orbits die Kurzzeitdynamik des Systems modelliert werden. Den Schlusspunkt dieser Untersuchungen bildete die klassische Simulation von zeitaufgelösten Pump-Probe-Spektren. Nachdem diese neue Methode mit periodischen Orbits an einem einfachen, stark idealisierten Modell erfolgreich getestet werden konnte, stellte sich jedoch schnell heraus, dass es in mehrdimensionalen Modellen mit einer Torsionsmode zur Beschreibung von Photoisomerisierungs-Prozessen nicht mehr so leicht möglich ist, solche einfachen periodischen Bahnen zu finden, weil die Zeitskalen von nuklearer und elektronischer Bewegung in diesen Systemen sehr unterschiedlich sind. Es ist deshalb in diesem Fall sinnvoll, das Konzept periodischer Orbits dahingehend zu erweitern, die klassische Dynamik in unterschiedliche allgemeine Bewegungstypen einzuteilen, die dann zur Interpretation von quantenmechanischen Ergebnissen herangezogen werden können. Dies stellt eine neue Möglichkeit dar, die komplizierte Wellenpaketdynamik auf physikalisch intuitive Weise zu veranschaulichen. Zusammenfassend kann man sagen, dass diese Arbeit Beiträge für das Verständnis nichtadiabatischer Quantendynamik liefert: Zum einen werden Fragestellungen von prinzipiellem Interesse diskutiert, wie z.B. der klassische Limes von Quantensystemen insbesondere bei chaotischem Verhalten, zum anderen werden durch die Anwendung von Analyseverfahren der klassischen nichtlinearen Dynamik auf vibronisch gekoppelte Molekülmodelle neue Gebiete erschlossen, die ein verbessertes theoretisches Verständnis experimenteller Ergebnisse im Bereich der Kurzzeitspektroskopie ermöglichen.
Die vorliegende Arbeit behandelt die Entwicklung und Überprüfung von Modellen zur Berechnung von Schwingungspektren von Peptiden und Proteinen. Solche Modelle verbinden die Konformationsstruktur eines Moleküls mit seinen Schwingungseigenschaften und sind demzufolge wichtig für die Interpretation der Schwingungspektren. Die im Rahmen dieser Arbeit durchgeführte theoretische Erforschung dieses Gebietes beschränkt sich auf die Betrachtung der Amide-I-Moden, welche aufgrund ihrer physikalischen Eigenschaften sich zur Untersuchung der Peptidkonformationen eignen. Die Arbeit kann prinzipiell in zwei Teile separiert werden. In dem ersten Teil werden Fragen betrachtet, die mit der Entwicklung des Schwingungshamiltonian verbunden sind. Im zweiten Teil wurden die erhaltenen Hamiltonian für die Berechnung der Schwingungspektren verwendet. Bei der Berechnung der Schwingungspektren wurden verschiedene spektroskopische Näherungen verwendet und erforscht. Die Entwicklung des Schwingungshamiltonian beinhaltet zwei Aufgaben. Die ab initio Parametrisierung des Schwingungshamiltonian von Dipeptiden, sowie die Analyse der Entwicklungsmethoden für Schwingungshamiltonian von Polypeptiden. Die Entwicklungsmethoden stützen sich auf ab initio berecheten Schwingungseigenschaften von Dipeptiden und/oder elektrostatische Modelle. Die ab initio Parametrisierung basiert auf einer Geometrieoptimierung und anschließender Berechnung von Normalmoden. Hierbei wurde die Abhängigkeit der Ergebnisse vom theoretischen Niveau und dem verwendeten Basissatz untersucht. Die Transformation der errechneten Normalmoden lieferte die Schwingungseigenschaften der lokale Amide-I-Mode. Die Lokalisierung der Normalmode folgt diversen Kriterien. Sie ist von der Wahl der Lokalmoden und somit implizit auch von der Art der Geometrieoptimierung abhängig. Mit dieser Arbeit konnte die Abhängigkeit der Ergebnisse von der Parameterwahl weitgehend aufgeklärt und eine für das Amide-I-System geeignet Parametrisierung gefunden werden. Im nächsten Arbeitsschritt wurde die Abhängigkeit der Amide-I-Schwingungseigenschaften von den Peptidseitenketten und terminalen Gruppen untersucht. Desweiteren wurden Methoden zur Formulierung der Hamiltonian für Polypeptide konzeptionell entwickelt. Diese Untersuchung ist außerordentlich wichtig, da direkte quantenmechanische Berechnungen von Polypeptiden zu zeitaufwendig sind. Solche Methoden beruhen auf dem sogenannten “Building-Block”-Ansatz und verschiedenen elektrostatischen Modellen. In dieser Arbeit wurden sowohl die einzelnen Methoden als auch ihre Kombination für die Entwicklung des Hamiltonians verwendet. Zur Abschätzung der Genauigkeit der verwendeten Methoden wurden Vergleichsrechnungen durchgeführt. Im zweiten Teil dieser Arbeit wurden die erhaltenen Schwingungshamiltonian zur Berechnung von Schwingungsspektren diverser gelöster Peptide angewandt. In diesem Zusammenhang konnte die Genauigkeit unterschiedlicher spektroskopischer Approximationen überprüft werden. Auf Grundlage der erhaltenen Ergebnisse können wir sagen, dass eine angemessene Beschreibung der konformationellen Verteilung und eine korrekte Berechnung des dynamischen Absorptionsspektrum gewährleistet ist. Was noch fehlt, ist ein hinreichend genaues quantenchemisches Modell für die Schwingungsfrequenzen eines gelösten Peptids. Diese Aufgabe stellt zur Zeit ein aktives Forschungsgebiet dar. Zuletzt wurde das Schwingungsspektrum eines sogenanten “Photoschaltbaren”-Peptids simuliert. Mit Hilfe des dafür aufgestellten Hamiltonians ist man in der Lage spektroskopische Beobachtungen auf Konformationsänderungen direkt zu übertragen.
Für das Verständnis der Proteinfaltung ist es von Interesse, die phi,psi-Torsionswinkelverteilung und deren Abhängigkeiten innerhalb einer Polypeptidkette zu kennen. Mit der in dieser Arbeit verwendeten Kombination aus MD-Simulation und NMR-Spektroskopie wird die Abhängigkeit der Konformationsverteilung kurzer alaninbasierter Modellpeptide mit einer Genauigkeit von 5 % bestimmt. Die Berechnung der thermischen Populationen der einzelnen Konformationen beruht auf einer Minimierung der Differenz aus experimentellen und berechneten skalaren Kopplungskonstanten. Trialanin populiert überwiegend den Bereich der Polyprolin Typ II Helix (~ 90 %) und daneben den beta-Faltblattbereich mit ca. 10%, jedoch nicht den alphaR-helicalen Bereich. Diese Konformationsverteilung ändert sich nicht signifikant mit zunehmender Kettenlänge in der Peptidreihe Ala3 bis Ala7. Das in der Seitenkette verzweigte Trivalin populiert dagegen alle drei Konformationsbereiche signifikant. Aufgrund der Periodizität der Torsionswinkel populiert Triglycin einen zusammenhängenden Bereich, der sich an den vier Ecken des Ramachandran-Diagramms befindet. Zudem befindet es sich in einem langsamen konformationellen Gleichgewicht zwischen der cis- und trans-Konformation der Peptidbindung. Die Temperaturabhängigkeit der Konformationsverteilung wird am Beispiel von Trialanin untersucht. Die 3J(HN,Ha) Kopplungskonstanten nehmen linear mit der Temperatur zu. Dies ist auf eine Zunahme des beta-Faltblattanteils zurückzuführen und kann theoretisch beschrieben werden. Die Konformationsverteilung der Trialaninsequenz innerhalb einer heteropolymeren Aminosäuresequenz ist von der Kettenlänge der an dem N- und C-Terminus angefügten heteropolymeren Aminosäuresequenz abhängig. Dies wird an zwei Peptiden, abgeleitet von der Sequenz des Proteins Lysozym aus Hühnereiweiß, gezeigt. Das kürzere Peptid hat an beiden Enden jeweils drei Aminosäurereste angefügt, das längere jeweils acht Aminosäurereste. Die Konformationsverteilung der Trialanisequenz des kürzeren Peptids entspricht nahezu der in der Peptidreihe Ala3 bis Ala7. Die Verteilung des längeren Peptids ist dagegen deutlich verschieden (~ 35% alphaR-helicaler Anteil). Die 1HN und 15N chemischen Verschiebungen der Trialaninsequenz des längeren Peptids sind mit denen des entfalteten Lysozym-Proteins identisch und demzufolge aller wahrscheinlichkeit nach auch die Konformationsverteilung. Kurze homopolymere Peptide eignen sich deshalb nicht als Modell für Aminosäuresequenzen in längeren heteropolymeren Peptiden.
Molecular dynamics (MD) simulation serves as an important and widely used computational tool to study molecular systems at an atomic resolution. No experimental technique is capable of generating a complete description of the dynamical structure of the biomolecules in their native solution environment. MD simulations allow us to study the dynamics and structure of the system and, moreover, helps in the interpretation of experimental observations. MD simulation was first introduced and applied by Alder and Wainwright in 1957 \cite{Alder57}. However, the first MD simulation of a macromolecule of biological interest was published 28 years ago \cite{McCammon77}. The simulation was concerned with the bovine pancreatic trypsin inhibitor (BPTI) protein, which has served as the hydrogen molecule'' of protein dynamics because of its small size, high stability, and relatively accurate X-ray structure available in 1977 \cite{Deisenhofer75}. This method is now widely used to tackle larger and more complex biological systems \cite{Groot01,Roux02} and has been facilitated by the development of fast and efficient methods for treating the long-range electrostatic interactions \cite{Essmann95}, the availability of faster parallel computers, and the continuous development of empirical molecular mechanical force fields \cite{Langley98,Cheatham99,Foloppe00}. It took several years until the first MD simulations of nucleic acid systems were performed \cite{Levitt83,Tidor83,Prabhakaran83,Nilsson86}. These investigations, which were also performed in vacuo, clearly demonstrated the importance of proper handling of electrostatics in a highly charged nucleic acid system, and different approaches, such as reduction of the phosphate charges and addition of hydrated counterions, have been applied to remedy this shortcoming and to maintain stable DNA structures. A few years later, the first MD simulation of a DNA molecule, including explicit water molecules and counterions was published \cite{Seibel85}. Various MD simulations on fully solvated RNA molecules with explicit inclusion of mobile ions indicated the importance of proper treatment of the environment of highly charged nucleic acids \cite{Lee95,Zichi95,Auffinger97,Auffinger99}. Given the central roles of RNA in the life of cells, it is important to understand the mechanism by which RNA forms three dimensional structures endowed with properties such as catalysis, ligand binding, and recognition of proteins. Furthermore, the increasing awareness of the essential role of RNA in controlling viral replication and in bacterial protein synthesis emphazises the potential of ribonucleicacids as targets for developing new antibacterial and new antiviral drugs. Driven by fruitful collaborations in the Sonderforschungsbereich RNA-Ligand interactions" the model RNA systems in this study include various RNA tetraloops and HIV-1 TAR RNA. For the latter system, the binding sites of heteroaromatic compounds have been studied employing automated docking calculations \cite{Goodsell90}. The results show that it is possible to use this tool to dock small rigid ligands to an RNA molecule, while large and flexible molecules are clearly problematic. The main part of this work is focused on MD simulations of RNA tetraloops.
Die vorliegende Arbeit beinhaltet die Untersuchung des Schwingungsenergietransfers im N-Methylazetamid (NMA) Molekül und von Wärmetransportprozessen im lichtschaltbaren 310-Aib-Oktapeptid. Dazu unternahmen wir Nichtgleichgewichts-MD-Simulationen und analysierten den zeitlichen Energiefluss in den jeweiligen Biomolekülen. Die Simulationsergebnisse wurden dann im Rahmen der in Kapitel 2 vorgestellten Theorien zum Energietransfer diskutiert. ...
Im Mittelpunkt der Arbeit steht die theoretische Beschreibung von ultraschnellen nichtadiabatischen cis-trans-Photoisomerisierunen in kondensierter Phase. Zu diesem Zweck wurde auf einen etablierten Modell-Hamilton-Operator zur Darstellung des isomerisierenden Systems zurückgegriffen und die Wechselwirkung mit der Umgebung im Rahmen der Redfield-Theorie behandelt. Eine gängige Näherung im Rahmen der Redfield-Theorie ist die Säkularnäherung, welche eine deutliche Reduktion des numerischen Aufwands bewirkt. Allerdings liefert die Säkularnäherung keine korrekte Beschreibung der Dynamik für Systeme, welche Regelmäßigkeiten im Eigenwertspektrum aufweisen, was für die in dieser Arbeit benutzten Isomerisierungsmodelle mit harmonischen Schwingungsmoden zutrifft. Andererseits verbietet sich für diese Systeme aus numerischer Sicht der Redfield-Algorithmus mit vollem Relaxationstensor. Daher wurde in dieser Arbeit ein nichtsäkularer Algorithmus entwickelt, der durch die Berücksichtigung der wichtigsten nichtsäkularen Terme eine adäquate Beschreibung der Dynamik im Rahmen der Redfield-Theorie liefert und gleichzeitig zu einer durchschnittlichen Reduktion der Rechenzeit auf ein Zehntel gegenüber der vollen Redfield-Rechnung führt. Im Rahmen der Redfield-Theorie wurden dann Dekohärenz- und dissipative Effekte für ein zweidimensionales Isomerisierungsmodell untersucht, wobei die unterschiedliche nichtadiabtische Dynamik an einer konischen Durchschneidung versus einer vermiedenen Kreuzung im Mittelpunkt des Interesses stand. Fazit dieser Studie ist, dass die konische Durchschneidung ein schnelles Abklingen anfänglicher Kohärenzen und eine effiziente Energieabgabe an die Umgebung bewirkt, was zu einer um eine Größenordnung schnelleren Isomerisierung gegenüber der vermiedenen Kreuzung führt. Daraus kann gefolgert werden, dass der photochemische Trichter tatsächlich der bevorzugte Reaktionsweg bei ultraschnellen internen Konversionsprozessen ist. Ein weiteres Anliegen dieser Arbeit war die Simulation von zeit- und frequenzaufgelösten Pump-Probe-Spektren für Photoreaktionen in dissipativer Umgebung. Hierzu wurde der Doorway-Window-Formalismus herangezogen, bei dem die Wechselwirkung der Pump- und Probepulse mit dem System im Doorway- bzw. Window-Operator enthalten ist. Für diese wurden unter der Annahme gaussförmiger Laserpulse durch analytische Integration der zweizeitigen Antwortfunktionen explizite Ausdrücke erhalten, die in den Redfield-Algorithmus integriert wurden. Somit existiert nun eine Methode zur Berechnung von Pump-Probe-Spektren, deren Skalierungsverhalten durch den Redfield-Algorithmus bestimmt wird. Diese Methode wurde dann angewandt für eine umfangreiche Modellstudie zu PumpProbe-Spektren von Isomerisierungsreaktionen in dissipativer Umgebung. Dabei wurden potentielle Probleme bei der Interpretation von transienten Spektren durch die Überlagerung und teilweisen Auslöschung der spektralen Beiträge zum Gesamtsignal aufgezeigt und diskutiert. In einer weiteren Studie wurde die Methode benutzt, um am Beispiel eines Morse-Oszillators zeitaufgelöste IR-Experimente zu simulieren, wie sie zur Gewinnung von modenselektiven Informationen über eine Reaktion durchgeführt werden.
Das genetische Material der Zellen besteht aus Molekülketten der Desoxyribonukleinsäure (DNA), die ein Träger der Erbinformation ist. In normalen Körperzellen wird die Erbinformation der DNA in eine andere Molekülkette, die sogenannte Ribonukleinsäure (RNA), übersetzt. Die RNA reguliert die Bildung von neuem Protein in der Zelle. Dass die RNA nicht bloß ein „Stempel“ ist, der die Informationen der DNA weitervermittelt, darin sind sich die Experten heute einig. RNA-Moleküle können Informationen speichern, katalytische Aktivitäten entfalten, sich perfekt tarnen, und sie regulieren auch als Produkt ihre eigene Synthese. Manche Viren enthalten ebenfalls RNA (oder DNA) und können so den Produktionsapparat der Zelle täuschen. Erkenntnisse über die Wechselwirkung dieser RNA mit natürlichen und synthetischen Liganden können zur Suche nach potentiellen Wirkstoffen beitragen. Nukleinsäuren sind lineare Biopolymere von grundlegenden Untereinheiten, die Nukleotide genannt werden und aus Adenin (A), Cytosin (C), Guanin (G), Urazil (U), und Thymin (T) zusammengesetzt sind. Sie sind jedoch in der Lage sich zu falten und so eine Doppel-Helixstruktur auszubilden. Diese besteht größtenteils aus den bekannten "Watson-Crick-Basenpaaren" (G-C und A-U oder A-T), die zur Stabilität der Struktur beitragen, sowie aus den weniger stabilen G-U-Paaren. Durch die Wechselwirkung zwischen verschiedenen Sekundärstrukturelementen entstehen Tertiärstrukturelemente, deren Struktur und Dynamik oft nur schwer experimentell zu bestimmen sind. Fortschritte in der RNA-Strukturanalyse wurden durch Röntgenkristallographie und Kernresonanzspektroskopie (NMR) möglich. Durch die Röntgenkristallographie wurden viele RNA-Eigenschaften festgestellt. Allerdings besteht keine Kristallstruktur für alle mögliche Einzelnfaser-RNA-Haarnadeln, weil diese immer dazu neigen, in eine linearen doppelte Faserform zu kristallisieren, die geringe biologische Bedeutung hat. Außerdem wurde mit Hilfe der NMR-Spektroskopie das dynamische Verhalten von RNA, z.B. Entfaltungsprozesse bei ansteigender Temperatur, beobachtet. Jedoch erlauben diese experimentellen Daten oft keine direkte mikroskopische Beschreibung der molekularen Prozesse. Molekulardynamik (MD)-Simulationen von biologischen Systemen ermöglichen es hingegen, diese Prozesse in atomischem Detail zu untersuchen. Die MD-Simulation beschreibt ein molekulares System auf atomarer Ebene mit Hilfe der klassischen Mechanik. Kräfte werden von empirischen Potentialen abgeleitet. Sie liefern zeitabhängige Trajektorien, die sich aus den Newton'schen Bewegungsgleichungen ergeben. Durch verbesserte Computerleistung, bessere Kraftfelder, und neu entwickelte genauere Methoden stimmen heutzutage MD-Simulationen von RNA mit experimentellen Daten immer besser überein. In meiner Doktorarbeit wurden MD-Simulationen durchgeführt um die Dynamik, die Struktur und insbesondere die Stabilität von RNA-Hairpins theoretisch zu beschreiben, um so ein erweitertes Verständnis für die dynamischen Vorgänge zu erhalten. Auch der SFB 579 der Universität Frankfurt beschäftigt sich mit RNA-Systemen. Erforscht wird unter anderem der D-Loop des Coxsackievirus B3 (CVB3), der Virenmyocarditis verursacht. Die Interpretation dieser experimentellen Daten wird durch MD-Simulation möglich. In dieser Arbeit wurden das GROMACS Software-Paket und das AMBER Kraftfeld verwendet, um das strukturelle, dynamische und thermische Verhalten der RNA-Hairpins mit Hilfe von MD-Simulationen auf atomarer Ebene zu untersuchen. Betrachtet wurden die 14-mer RNA-Hairpins, uCACGg und cUUCGg. Die verfügbaren NMR-Strukturen zeigen, dass das uCACGg-Tetraloop auffallend ähnlich in der gesamten Geometrie und den Wasserstoffbindungen zu der experimentellen Struktur des cUUCGg-Tetraloop ist, obwohl die schließende Basenpaarsequenz der beiden Tetraloops unterschiedlich sind. Trotz beachtlicher struktureller Ähnlichkeit unterscheiden sich allerdings die uCACGg und cUUCGg Tetraloops in Funktionalität und Thermostabilität. Zunächst orientiert sich unser erstes Bemühen an der Frage nach einem guten Modell für RNA-Hairpins und Simulationsbedingungen, um die zu untersuchenden RNA-Hairpins in Wasser möglichst realitätsnah zu simulieren. Erstens werden drei Versionen des biomolekularen AMBER-Kraftfelds geprüft, indem man die 60 ns Simulationen des 14-mer uCACGg-Hairpins durchführt. Die simulierten strukturellen Eigenschaften und Atomfluktuationen zeigen hohe Ähnlichkeiten in den drei Kraftfeldern. Darüber hinaus stimmen die von MD-Simulationen berechneten Atomkernabstände mit den experimentellen NMR-Daten gut überein. Die gute Übereinstimmung zwischen den Simulationen und den strukturellen NMR Daten belegt die Fähigkeit des AMBER-Kraftfelds zur Beschreibung der strukturellen Eigenschaft von kleinen RNA-Hairpins. Anschließend werden die Einflüsse der Methoden, welche die langreichweitigen, elektrostatischen Wechselwirkungen beschreiben, auf die strukturellen Eigenschaften untersucht. Insbesondere werden die Ergebnisse der Reaktionfeld-Methode mit denen der Particle Mesh Ewald (PME)-Methode verglichen. Es zeigt sich, dass die PME-Methode die elektrostatischen Wechselwirkungen am besten beschreibt, auch wenn die Simulationen der beiden Methoden Ähnlichkeit in der Struktur-Stabilität und der Atomfluktuation bei niedriger Natriumkonzentration aufweisen. Drittens wird der Kationseffekt auf die RNA-Stabilität untersucht. Betrachtet wurden zwei unterschiedliche Kationen (ein- und zweiwertig) und verschiedene Konzentrationen. Die Simulationen weisen darauf hin, dass sich die Metallionen in der Affinität zum RNA-Hairpin unterscheiden, wenn Na+ und/oder Mg2+ als Gegenionen verwendet werden. Weiterhin wird gezeigt, dass sich die bevorzugten Positionen der Na+-Ionen in der großen Furche (major groove) des RNA-Hairpins befinden. Insbesondere die Anlagerungsort der Na+-Ionen liegen in der Nähe des schließenden Basenpaar U5-G10. Im Vergleich zu Na+-Ionen lagern sich Mg2+-Ionen sowohl an die RNA-Basen U3, A4-U11, und die Phosphat-Gruppe, als auch an das schließenden Basenpaar U5-G10 an. Bestätigt werden die Modelle und Simulationsbedingungen durch den Vergleich von Parametern, die sowohl experimentell als auch durch Simulationen ermittelt werden können. Ferner erlauben MD-Simulationen Einblick in das System, indem sie detallierte Konformations- und andere Verteilungen liefern. In der vorliegenden Arbeit wurden die Einflüsse der Loopsequenz und des schließenden Basenpaares auf die Verteilung der Konformationen, der internen Bewegungen, und auf die Thermostabilität von zwei RNA-Hairpins mit Hilfe dieser Modelle untersucht. Zunächst wurden die strukturellen Eigenschaften bei Raumtemperatur ausgewertet. Die starken strukturellen Ähnlichkeiten und die gute Übereinstimmung mit NMR-Daten bestätigen die Hypothese, dass die zwei Tetraloops zur gleichen “erweiterten” RNA-Familie gehören. Diese zwei Hairpins haben ähnliche Lösemittelzugängliche Oberflächen (solvent accessible surface), wobei deren Lösemittel zugänglichen funktionellen Gruppen unterschiedlich sind. Weiterhin weist das uCACGg-Hairpin eine stärkere Tendenz auf Wasserstoffe abzugeben als das cUUCGg-Hairpin, was in den unterschiedlichen Bindungsaffinitäten zwischen diesen Hairpins und der viralen Protease begründet liegt. Darüber hinaus wurde der Faltungs- und Entfaltungsprozess mit Hilfe der Replica-Exchange-Molekulardynamik-Simulationen untersucht. Diese Untersuchung zielt auf das bessere Verständnis der unterschiedlichen Thermostabilität der Hairpins, indem sie die möglichen Zwischenprodukte im atomaren Detail liefern. Sowohl experimentell als auch von den MD-Simulationen ergibt sich eine Differenz in den Schmelztemperaturen der beiden Hairpins von ungefähr 20 K. Allerdings sind die von MD beobachteten Schmelztemperaturen 20 % höher als die von Experiment zu ansehende Wert. Die Ergebnisse machen deutlich, dass die Schmelztemperaturdifferenz nicht auf die Unterschiede in der Sequenz, in der Struktur, oder in der Dynamik der Loops zurückführen sind, sondern auf die Unterschiede der Basenpaaren in den Stämmen. Weiterhin wird gezeigt, dass sich das uCACGg-Hairpin einerseits kooperativ entfaltet, und die Entfaltung des cCACGg-Hairpins anderseits weniger kooperativ stattfindet. Um die schnelle interne Dynamik der uCACGg- und cUUCGg-Hairpins zu untersuchen, erlauben die Simulationen von 50 ns eine akurate Beschreibung der schnellen internen Bewegung der RNA-Hairpin, obwohl der den Hairpins zugängliche Konformationsraum nicht vollständig abgedeckt wird. Die NMR-Relaxationsparameter, die mit Hilfe der MD-Simulationen zurückgerechnet wurden, bestätigen das Modell und die Simulationsbedingungen der MD-Simulationen. Im Hinblick auf die Übereinstimmung kann man den besten Ansatz zur Berechnung der NMR-Ordnungsparameter bestimmen. In dieser Arbeit wurden drei verschiedene Ansätze angewandt, nämlich das Fitting von 100 ps auf modellfreiem Ansatz nach Lipari-Szabo, equilibrium average, und das Gaussian Axial Fluctuation (GAF)-Modell. Die zwei letzteren können nur qualitativ mit den experimentellen Daten übereinstimmen. Die NMR-Ordnungsparameter können mit Hilfe des Modells von Lipari-Szabo richtig ermittelt werden, wenn sich die interne Bewegung in kleineren Zeitskalen als zur Gesamtbewegung vollzieht. Vorausetzung für die Berechnung dieses Modells ist aber, dass das Fitting der internen Korrelationsfunktionen nur auf den ersten Teil von 100 ps der Korrelationsfunktionen eingesetzt wird. Die berechneten Ordnungsparameter deuten auf ein unterschiedliches Verhalten der beiden Hairpins besonders im Loop-Bereich hin. Die konformationelle Umordnung, die beim UUCG-Loop beobachtet wurde, tritt beim CACG-Loop nicht ein. Zusammenfassend lässt sich sagen, dass es durch den Einsatz von MD Simulationen ermöglicht wird, die strukturellen und dynamischen Eigenschaften der RNA-Systeme auf atomarer Ebene zu untersuchen. Als Schlussfolgerung zeigt diese Doktorarbeit, dass sich die Studie der konformationell Dynamik der RNA-Systeme durch die Kombination aus MD-Simulation und NMR-Spektroskopie sowie der Leistungsfähigkeit der MD-Simulationen, die die interne Bewegungen deutlich beschreiben können, untersuchen lässt.