Refine
Year of publication
Document Type
- Doctoral Thesis (67)
Has Fulltext
- yes (67)
Is part of the Bibliography
- no (67)
Keywords
- NMR-Spektroskopie (7)
- RNA (6)
- NMR (4)
- NMR spectroscopy (4)
- Proteinfaltung (3)
- Amyloid (2)
- DNA (2)
- Dynamik (2)
- G-Quadruplex (2)
- Kinetik (2)
- Ligand (2)
- Magnetische Kernresonanz (2)
- Molekülstruktur (2)
- NMR-spectroscopy (2)
- Peptide (2)
- RNS (2)
- dynamics (2)
- molecular structure (2)
- protein folding (2)
- Acute myeloid leukemia (AML) (1)
- Alzheimer's disease (1)
- Alzheimer-Krankheit (1)
- Aminosäurensequenz (1)
- Amyloidkern (1)
- Aptamer (1)
- Biophysik (1)
- Biopysikalische Chemie (1)
- CIDNP (1)
- Calcium (1)
- Chemie und Pharmazie (1)
- DEAD-Box Helicase (1)
- Dreidimensionale NMR-Spektroskopie (1)
- Enantioselektive Katalyse (1)
- Faltungsdynamik (1)
- Fettsäuresynthase (1)
- Fibrillen (1)
- Fibrils (1)
- Force Field Development (1)
- Frankfurt <Main> / Universität / Fachbereich Biochemie (1)
- H/ACA-RNP (1)
- Harnstoff (1)
- Helicase (1)
- Inhibition (1)
- Inhibitor-binding (1)
- J coupling constants (1)
- Konformation (1)
- Kopplungskonstante (1)
- Kopplungskonstanten (1)
- Lactalbumin (1)
- Ligand <Biochemie> (1)
- Lysozym (1)
- MD Simulation (1)
- MD simulations (1)
- MD-Simulation (1)
- Magnesium (1)
- Mehrdimensionale NMR-Spektroskopie (1)
- Metabolism (1)
- Molekulardynamik (1)
- NMR Spectroscopy (1)
- Nukleinsäuren (1)
- Oligonukleotid (1)
- Organische Synthese (1)
- Organokatalyse (1)
- Photochemie (1)
- Photolabile Schutzgruppen (1)
- PrP27-30 (1)
- Prion Protein (1)
- Prionprotein (1)
- Protein (1)
- Protein-Tyrosin-Phosphatase (1)
- Proteintyrosinphosphatase (1)
- Pseudouridine (1)
- Pseudouridylation (1)
- Puromycin (1)
- RDCs (1)
- RNA-folding (1)
- RbfA (1)
- Regulation der Genexpression (1)
- RhlB (1)
- Ribonucleoproteins (1)
- Ribosom (1)
- Ribosome (1)
- Riboswitch, RNA, translation, NMR, smFRET (1)
- Ryanodin (1)
- SERCA (1)
- Spin-Spin-Kopplungskonstante (1)
- Spinlabel (1)
- TSE (1)
- Tau-Protein (1)
- Wirkstoff-Rezeptor-Bindung (1)
- Zellzyklus (1)
- Zyklisierung (1)
- aggregopathy (1)
- cMYC (1)
- cell cycle (1)
- elegans (1)
- hTERT (1)
- hen egg white lysozyme (1)
- ligand interactions (1)
- multidimensional nmr spectroscopy (1)
- nicht-native Proteine (1)
- non-native proteins (1)
- peptides (1)
- photoschalbare (1)
- rRNA (1)
- residual dipolar couplings (1)
- skalare Kopplungskonstanten (1)
- smFRET (1)
- tau protein (1)
- unfolded proteins (1)
Institute
- Biochemie und Chemie (36)
- Biochemie, Chemie und Pharmazie (26)
- Biowissenschaften (2)
- Extern (1)
- Georg-Speyer-Haus (1)
- Pharmazie (1)
- Physik (1)
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.
In this thesis the three dimensional solution strucutre of the RbfA protein from Thermotoga maritima was solved using multidimensional heteronuclear NMR spectroscopy. The RbfA protein binds to the helix I region of the 16S rRNA. To gain insights into the binding mode of RbfA to its target, a second RbfA construct from Helicobacter pylori was used. Comparison of the RbfA proteins with the published structure of RbfA from Escherichia coli, led to studies concerning the differences between proteins from thermophile and mesophile systems. In the second part of this thesis the native binding motive of the RbfA protein was identified. The RbfA protein binds to an alternate helix fold within the pre-sequence of the immature 16S rRNA.
One of the most important tasks in chemistry and especially in structural biology has always been the elucidation of three-dimensional molecular structures - either of small molecules or large biopolymers. Among the (bio)physical methods to acquire structural data at atomic resolution electron paramagnetic resonance (EPR) spectroscopy is the most valuable technique for obtaining structural information about many different kinds of paramagnetic species. In biological systems, either paramagnetic metal ions/clusters, transient paramagnetic intermediates in electron transfer processes or artificially attached stable spin labels can be found. The usual approach to interpret EPR spectra is to perform simulations based on the so-called spin Hamiltonian (SH). This means that the well-defined numerical parameters (tensors) in the SH representing different types of interaction are obtained by fitting the experimental data. The SH parameters include electronic g-values, hyperfine coupling (HFC) and quadrupole coupling (&C) constants, zero-field splittings and constants to describe exchange and dipolar interactions between electron spin systems. However, since the SH only contains spin degrees of freedom, a direct translation of the SH EPR parameters into structural information is not straightforward. Therefore, methods to predict such SH interaction parameters starting from molecular structures are required. In this thesis it was investigated whether quantum chemical calculations of EPR parameters based on density functional theory (DFT) methods may be employed to overcome these problems thus enabling a correlation of experimental EPR data with molecular structure. It was the central goal of this work to point out the potential of a fruitful interplay between quantum chemistry and experiment and to study how both can benefit from each other. For this purpose DFT methods were applied to a variety of organic radical or transition metal systems to calculate different EPR parameters. Using the 'broken symmetry' formalism it was possible to compute the exchange coupling constant for a nitroxide biradical and furthermore decompose the exchange mechanism in different through-bond and through-space interactions. Spin density distributions, 14N and 1H HFC constants as well as dipole moments and polarizabilities were computed for a number of aromatic nitroxides to examine their properties and select promising candidates which may serve as DNA-intercalating spin labels. Systematic investigations of the influence of hydrogen bond geometry on the 14N QC parameters for imidazole-water and methylimidazole-benzosemiquinone complexes lead to the conclusion that especially the imidazole amino nitrogen &C parameters are very sensitive probes of the bond geometry, in particular of the hydrogen bond length. The results of this study may be applied to biological systems, e.g. to gain structural information about quinone binding sites. Moreover, quantum chemical methods were applied to elucidate the structure of a nitrogen-centered radical intermediate in the inhibition process of ribonucleotide reductase (RNR). It was possible to find a molecular structure in accordance with all experimentally available data, thus revealing the longsought structure of the No radical and providing evidence for the trapping of a 3'-ketonucleotide in the reduction process catalyzed by RNR. To test the capability of modern DFT methods to predict g- and molybdenum HFC tensors for MoV complexes, validation studies were carried out. Comparison of computed EPR parameters of a number of MoV compounds with corresponding experimental values showed that g- and HFC tensors could be predicted in good accuracy, although some systematic errors of the computational methods have to be considered for such heavy 4d1 transition meta1 systems. Furthermore, DFT calculations on a Mn2+ binding site model of the hammerhead ribozyme allowed to conclude that the structure of the binding site as studied by EPR spectroscopy in frozen solution is very likely to be identical to the site found occupied by Mn2+ in crystals. Finally, computational methods were employed to aid in the structural characterization of the Mn2+ binding site in Ras (rat sarcoma protein) by providing accurate starting parameters for spectral simulations and furthermore helping to interpret the experimental data. In conclusion, it was demonstrated in this thesis that the combination of sophisticated experimental and quantum chemical methods represents a powerful approach in the field of EPR spectroscopy and that it may be essential to employ EPR parameter computations to extract the full information content from EPR spectra. Therefore, great potential lies in future applications of DFT methods to the large number of systems where detailed and reliable experimental data is available but where an unequivocal correlation of these data with structural information is still lacking.
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.
This cumulative thesis discusses the development of optimized force field parameters for Magnesium and resulting improved simulations of Magnesium-RNA interactions, including the in silico exploration of binding sites. This thesis is based on four publications as well as unpublished data. A fifth publication that was written during the time of the Ph.D. is discussed in the Appendix. This publication analyzes monovalent ion-specific effects at mica surfaces.
Nucleic acids in general and RNA in particular are fundamental to life itself. Especially in the folding and function of RNA, metal cations are crucial to screen the negatively charged nucleic acid backbones to allow for complex functional structures. They stabilize the tertiary structure of RNA and even drive its folding. Furthermore, similarly to proteins, RNAs can catalyze multiple reactions, rather than consisting of the 20 amino acids of a protein, RNA constitues of only four different building blocks. Metal cations play an important role here as additional cofactors. One essential ion is Magnesium (Mg2+), commonly referred to as the most important cofactor for nucleic acids. Mg2+ carries two positive charges. Its comparably small size and high charge result in a high charge density that has strong polarizing effects on its surroundings. Furthermore, Mg2+ forms a sharply defined first hydration shell with an integer number of coordinating water molecules. As a result, an exclusion zone exists around the ion within which no water molecules are observed. Moreover, Mg2+ displays a high solvation free energy and a low exchange rate of waters from its first hydration shell. Finally, it contains a strong preference towards oxygens . Together, this makes Mg2+ a particularly well suited interaction partner for the charged non-bridging phosphate oxygens on nucleic acid backbones and explains its crucial biological role.
The immense number of physiological and technological functions and applications indicates the significant scientific attention Mg2+ received. In experimental studies, however, severe difficulties arise for multiple reasons: Mg2+ is spectroscopically silent and cannot be detected directly by resonance techniques like NMR or EPR. Indirect observation is possible, either by detecting changes in the overall RNA structure with and without bound Mg2+, or by replacing the Mg2+ ion with another spectroscopically visible ion. In the latter, however, it cannot be guaranteed that the altered ion does not also alter the interaction site or even the whole structure. Another detection method is X-ray crystallography, but here challenges arise from Mg2+ being almost indistinguish- able from other ions as well as from water if not for very high resolutions and precise stereochemical considerations.
Alternatively, molecular dynamics (MD) simulations can be performed, with the power of adding atomistic insight to the interplay of metal cations and nucleic acids. MD simulations, however, are only as accurate as their underlying interaction models and the development of accurate models for the description of Mg2+ faces challenges especially in describing three properties:
(i) Polarizability. Commonly used simple models like the 12-6 type Lennard-Jones model typically fail to reproduce simultaneously thermodynamic and structural properties of a single ion in water. Alternative strategies include the use of a 12-6-4 type Lennard-Jones potential as proposed by Li and Merz, where the additional r−4 term explicitly accounts for polarization effects. The resulting Lennard-Jones potential is thereby more attractive and more long-ranged than for typical models of the 12-6 type.
(ii) Kinetics. Most Mg2+ models either fully ignore considerations about the timescales on which water exchanges from the first hydration shell of the ion or use inappropriate methodology to calculate the underlying kinetics. A realistic characterization of the involved timescales is imperative to be able to describe a seemingly simple process like the transition from inner-to-outer sphere binding and vice versa. This transition governs most biochemical reactions involving Mg2+ and therefore subsequent processes can only by as fast as the transition itself. However, already the previous step – the exchange of a water from the first hydration shell of the ion – is described my current Mg2+ models up to four orders of magnitude too slowly, which makes the observation of such events on the timescale of a typical simulation difficult or even impossible. Alln ́er et al. [48] as well as Lemkul and MacKerell explicitly considered the exchange rate into their parameter optimization procedure. To compute the rate, both studies applied Transition State Theory along a single reaction coordinate – the distance towards one of the exchanging waters. However, it could be shown that the water exchange from the first hydration shell requires at least the consideration of both exchanging water molecules in order to be able to realistically record the underlying rate using Transition State Theory. Furthermore, the model of Alln ́er et al. significantly underestimates the free energy of solvation of the ion.
(iii) Interactions between Mg2+ and nucleic acids. Typically, ionic force field parame- terization concentrates on the optimization of solution properties. The trans- ferability of these solution optimized parameters towards interactions with biomolecules, however, often fails.
Das Tau-Protein bildet im Verlauf von zahlreichen Demenzen, mit Morbus Alzheimer als die bekannteste unter ihnen, Aggregate, die sogenannten „neurofibrillären Geflechte“, die aus Fibrillen, den „paired helical filaments“ (PHFs) bestehen. Außerdem ist das Tau-Protein ein essentielles „microtubule-associated protein“, welches für die Aufrechterhaltung des neuronalen Zellmetabolismus notwendig ist. Dies veranlasste uns dazu, das Tau-Protein als Monomer, in der Mikrotubuli-gebundenen Form und als Fibrille zu charakterisieren. Die Technik, die wir hierzu verwendeten war die NMR-Spektroskopie, die als einzige strukturaufklärende Technik mit atomarer Auflösung dazu in der Lage ist, intrinsisch ungeordnete Proteine, wie das Tau-Protein, zu charakterisieren. Zwar war die Signalzuordnung der NMR-Spektren eine große Herausforderung, dennoch war es möglich praktisch alle Rückgratresonanzen sogar für die längste Tau-Isoform htau40 mit 441 Resten erfolgreich eindeutig zu identifizieren. Mit Hilfe der Zuordnung war es möglich das Tau-Protein bezüglich residualer Strukturelemente, Rückgratdynamik und Bindungsverhalten zu untersuchen. Wir konnten zeigen, dass in der C-terminalen Hälfte des Tau-Proteins, in welcher eine charakteristische Domäne vorliegt, die durch vier imperfekte „repeat“-Regionen (Länge ist jeweils ca. 31 Reste) gekennzeichnet ist, partieller beta-Strukturcharakter vorliegt. Ebenfalls weist diese Region eine verhältnismäßig hohe Rigidität auf. Aus diesem Grund betrachten wird diesen sequentiellen Bereich als den Aggregationskeim, was auch durch die Beobachtung verstärkt wird, das genau diese Zone den rigiden Teil der Fibrillen bildet. Diese aggregationsanfällige Region bindet gleichzeit stark an Mikrotubuli, wodurch ihre Pathogenität im gesunden biologischen System blockiert sein sollte. Mutationen oder Instabilität in den Mikrotubuli können jedoch dazu führen, dass immer höhere Mengen an Tau in freier gelöster Form vorliegen, sich Dimere ausbilden, welche dann weiter aggregieren und schließlich PHFs bilden, die eine starke cross-beta-Struktur aufweisen. Die übrigen Bereiche, wie die N-terminale Hälfe oder die äußersten 50 C-terminalen Reste weisen hingegen einen partiellen alpha-helikalen Charakter auf und eine höhere Peptidrückgratflexibilität. Deshalb kann man diese Elemente als aggregationsinhibierend betrachten. Das genauere Zusammenspiel dieser Elemente muss noch aufbauend auf der vorliegenden Dissertation verstärkt im Detail untersucht werden.
RNA research is very important since RNA molecules are involved in various gene regulatory mechanisms as well as pathways of cell physiology and disease development.1 RNAs have evolved from being considered as carriers of genetic information from DNA to proteins, with the three major types of RNA involved in protein synthesis, including messenger RNA (mRNA), transfer RNA (tRNA), and ribosomal RNA (rRNA).2 In addition to the RNAs involved in protein synthesis numerous regulatory non-coding RNAs (ncRNAs) have been discovered in the transcriptome. The regulatory ncRNAs are classified into small ncRNAs (sncRNAs) with transcripts less than 200 nucleotides (nt) and long non-coding RNAs (lncRNAs) with more than 200 nt.3
LncRNAs represent the most diverse and versatile class of ncRNAs that can regulate cellular functions of chromatin modification, transcription, and post-transcription through multiple mechanisms.4 They are involved in the formation of RNA:protein, RNA:RNA and RNA:DNA complexes as part of their gene regulatory mechanism.4,5 The RNA:DNA interactions can be divided into RNA:DNA heteroduplex formation, also called R-loops, and RNA:DNA:DNA triplex formation. In triplex formation, RNA binds to the major groove of double-stranded DNA through Hoogsteen or reverse Hoogsteen hydrogen bonding, resulting in parallel or anti-parallel triplexes, respectively. In vitro studies have confirmed the formation of RNA:DNA:DNA triplexes.6 However, the extent to which these interactions occur in cells and their effects on cellular function are still not understood, which is why these structures are so exciting to study (Chapter I RNA:DNA:DNA Triplexes).
This cumulative thesis investigates several functional and regulatory important RNAs. The first project involves the improved biochemical and biophysical characterization of RNA:DNA:DNA triplex formation between lncRNAs of interest and their target genes. Triplex formation was confirmed by a series of experiments including electromobility shift assays (EMSA), thermal melting assays, circular dichroism (CD), and liquid state nuclear magnetic resonance (NMR) spectroscopy. The following is a summary of the main findings of these publications.
In research article 5.1, the oxygen-sensitive HIF1α-AS1 was identified as a functionally important triplex-forming lncRNA in human endothelial cells using a combination of bioinformatics techniques, RNA/DNA pulldown, and biophysical experiments. Through RNA:DNA:DNA triplex formation, endogenous HIF1α-AS1 decreases the expression of several genes, including EPH receptor A2 (EPHA2) and adrenomedullin (ADM), by acting as an adaptor for the repressive human silencing hub (HUSH) complex, which has been studied by our collaborators in the groups of Leisegang and Brandes.
2) Triplex formation between HIF1α-AS1 and the target genes EPHA2 and ADM was investigated in biochemical and biophysical studies. The EMSA results indicated that HIF1α-AS1 forms a low mobility RNA:DNA:DNA triplex complex with the EPHA2 DNA target sequence. The CD spectrum of the triplex showed distinct features compared to the EPHA2 DNA duplex and the RNA:DNA heteroduplex. Melting curve analysis revealed a biphasic melting transition for triplexes, with a first melting point corresponding to the dissociation of the RNA strand with melting of the Hoogsteen hydrogen bonds. The second, higher melting temperature corresponds to the melting of stronger Watson-Crick base pairing. Stabilized triplexes were formed using an intramolecular EPHA2 DNA duplex hairpin construct in which both DNA strands were attached to a 5 nucleotide (nt) thymidine linker. This approach allowed improved triplex formation with lower RNA equivalents and higher melting temperatures. By NMR spectroscopy, the triplex characteristic signals were observed in the 1H NMR spectrum, the imino signals in a spectral region between 9 and 12 ppm resulting from the Hoogsteen base pairing. To elucidate the structural and sequence specific Hoogsteen base pairs 2D 1H,1H-NOESY measurements of the EPHA2 DNA duplex and the HIF1α-AS1:EPHA2 triplex were performed. The 1H,1H-NOESY spectrum of the HIF1α-AS1:EPHA2 triplex with a 10-fold excess of RNA was semi-quantitatively analyzed for changes in the DNA duplex spectrum. We discovered, strong and moderate attenuation of cross peak intensities in the imino region of the NOESY spectrum. This attenuation was proposed to result from weakening of Watson-Crick base pairing by Hoogsteen hydrogen bonding induced by RNA binding. The Hoogsteen interactions can be mapped based on the analysis of the cross peak attenuation in the NOESY spectra, which we used to generate a structural model of the RNA:DNA:DNA triplex. These biophysical results support the physiological function of HIF1α as a triplex-forming lncRNA that recruits the HUSH-epigenetic silencing complex to specific target genes such as EPHA2 and ADM, thereby silencing their gene expression through RNA:DNA:DNA triplex formation.