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)
Institute
- Biochemie und Chemie (36)
- Biochemie, Chemie und Pharmazie (26)
- Biowissenschaften (2)
- Extern (1)
- Georg-Speyer-Haus (1)
- Pharmazie (1)
- Physik (1)
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.
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.
Die vorliegende Dissertation mit dem Titel “Structural dynamics of eukaryotic H/ACA RNPs from Saccharomyces cerevisiae & Structural dynamics of the Guanidine-II riboswitch from Escherichia coli” besteht aus zwei Projekten. Das erste Projekt befasst sich mit den eukaryotischen H/ACA Ribonukleoproteinen (RNP) aus der Hefe. Diese können sequenzspezifisch in der RNA ein Uridin Nukleotid in das Rotationsisomer Pseudouridin (Ψ) umwandeln. Die H/ACA RNPs bestehen aus einer Leit-RNA und vier Proteinen, der katalytisch aktiven Pseudouridylase Cbf5, Nhp2, Gar1 und Nop10. Die Leit-RNA besteht in Eukaryoten konserviert aus zwei Haarnadelstrukturen, die von einem H-Box oder ACA-Box Sequenzmotiv gefolgt sind. In jeder dieser Haarnadeln befindet sich ein ungepaarter Bereich, die sogenannte Pseudouridylierungstasche, wo durch komplementäre Basenpaarung die Ziel-RNA gebunden wird. Fehlerhafte H/ACA RNPs können beim Menschen zu schweren Krankheiten wie verschiedenen Krebsarten oder dem Knochenmarksversagen Dyskeratosis congenita führen, aber sie bieten auch Möglichkeiten zum Einsatz als Therapiemethode. In dieser Arbeit wurde hauptsächlich der zweiteilige Aufbau der H/ACA RNPs untersucht.
Dafür wurden zunächst die einzelnen Komponenten hergestellt werden. Cbf5, Nop10 und Gar1 wurden zusammen heterolog in E. coli exprimiert und gereinigt. Außerdem wurden mehrere Deletionsvarianten von Gar1 hergestellt. Zusätzlich wurde die Leit-RNA unmarkiert über T7 Transkription synthetisiert, sowie sechs verschiedene FRET-Konstrukte mit verschiedenen Markierungschemas der Fluorophore Cy3 und Cy5 über DNA-geschiente Ligation. Anschließend wurde über Größenausschlusschromatographie und radioaktiven Aktivitätsassays geprüft, dass sich die aktiven H/ACA RNPs in vitro aus den einzelnen Komponenten rekonstituieren lassen.
In smFRET Experimenten wurden einzelne Haarnadelstrukturen mit dem zweiteiligen Komplexen verglichen. Dabei konnte gezeigt werden, dass die H3 Haarnadel durch die Anwesenheit von H5 dynamischer und heterogener wurde, während H5 überwiegend unbeeinflusst war. Außerdem konnte die dreidimensionale Orientierung der Haarnadelstrukturen in verschiedenen Assemblierungsschritten mittels smFRET untersucht werden. Hier deutete sich an, dass in Abwesenheit von Proteinen beide Haarnadeln eher entgegengesetzt stehen als in einer parallelen Konformation. Cbf5 scheint den Linker zwischen den Beiden auszustrecken bzw. zu orientieren und die Haarnadelstrukturen etwas gegeneinander zu neigen. Ein Zusammenspiel von Nhp2 und Gar1 war nötig um die oberen Bereiche der Haarnadeln zusammenzuziehen. Es konnte auch ein Modell für den vollen H/ACA RNP vorgeschlagen werden. Im kompletten Komplex könnte das Zusammenziehen der Haarnadelstrukturen durch Nhp2 und Gar1 mit dem Effekt von Cbf5 konkurrieren und könnte hauptsächlich den oberen Bereich von H3 betreffen. Zum Schluss wurde das Zusammenspiel von Gar1 und Nhp2 auf eine Abhängigkeit von den RGG Domänen von Gar1 hin untersucht. Hier besteht möglicherweise eine Hierarchie, die eine Kooperativität von den N- und C-terminalen Domänen benötigt.
Das zweite Projekt befasst sich mit dem Guanidin-II Riboschalter aus E. coli. Der Riboschalter kann das toxische Molekül Guanidinium (Gdm+) spezifisch in seiner Aptamerdomäne binden und dadurch die Genexpression von Proteinen zur Detoxifizierung von Gdm+ aktivieren. Der Riboschalter besteht aus zwei Haarnadelstrukturen, mit einer Schleife, die aus der Sequenz ACGR besteht, wobei R ein Purin ist. In einem vorgeschlagenen Modell soll die Ribosomenbindestelle (Shine-Dalgarno Sequenz) in Abwesenheit von Ligand mit dem Linker komplementär Basenpaaren und so die Translation verhindern. Mit Ligand würde sich dann eine Schleifen-Schleifen Interaktion mit den beiden CG Basen ausbilden, wodurch die Anti-Shine-Dalgarno Sequenz nicht mehr zugänglich wäre. Bisherige Studien arbeiteten zumeist nur mit der Aptamerdomäne, den einzelnen Haarnadeln oder noch kleineren Elementen. In dieser Arbeit wurden die Strukturdynamiken von verschiedenen Längen, auch mit der Expressionsplatform, untersucht. Außerdem wurden verschiedene Mutationen analysiert und die Effekte auf den Riboschalter in seiner natürlichen Umgebung in E. coli.
Zunächst mussten insgesamt 24 FRET-Konstrukte hergestellt werden, die sich in Länge, Markierungsschema und Mutationen unterschieden. Hierfür wurde DNA-geschiente Ligation verwendet. Dank der verschiedenen Fluorophorpositionen konnte ein konformationelles Modell für die Aptamerdomäne vorgeschlagen werden. In diesem Modell könnte in Abwesenheit von Ionen das Aptamer offen vorliegen. Durch Mg2+ würde sich bereits eine lockere Schleifen-Schleifen Interaktion ausbilden. Zusätzlich deuten die Ergebnisse auf eine neue Konformation hin, der stabilisierten Schleifen-Schleifen Interaktion, bei der der Linker zusätzlich mit den Haarnadelstrukturen interagiert, beispielswese mit den Purinen an der vierten Schleifenposition...
N6-methyladenosine (m6A) is the most abundant and well understood modification in eukaryotic mRNA and was first identified in polyadenylated parts of the mRNA.The distinct distribution of m6A in the transcriptome with special enrichment in long internal exons, 39UTRs and around stop codons was uncovered by early biochemical work and later on antibody based sequencing techniques. The so called m6A writer, reader and eraser machinery is responsible for the dynamic and with that regulatory nature of the m6A modification. As m6A writer, the human N6-methyltransferase complex (MTC) cotranscriptionally methylates the central adenine within a RRACH (preferably GGACU) sequence context to form m6A in the nascent RNA chain.9–15 The catalytic core of the complex is formed by the two proteins METTL3 and METTL14, with the active site located in the methyltransferase domain (MTD) of METTL3.16–18 The DPPW motif near the methyl donor S-adenosylmethionine (SAM) binding site in this MTD was postulated to bind the target adenine during catalysis. Moreover, a positively charged groove in the METTL3-METTL14 interface, the C-terminal RGG domain in METTL14 and the zinc finger motifs in METTL3 were identified as important domains for RNA binding. However, to date there are no full-length or substrate-RNA-bound structures of the catalytic METTL3-METTL14 complex.
In addition, a set of accessory proteins assembles to the METTL3-METTL14 heterodimer to form the full MTC, mediated by WTAP that firmly binds to the N-terminal leader helix in METTL3.20 WTAP was shown to locate the whole complex to the nuclear speckles and can modulate m6A deposition to specific sites in the RNA. Moreover, WTAP acts as binding platform for other accessory proteins including VIRMA, RBM15, ZC3H13 and HAKAI that are mostly identified to mediate position specific methylation. For example, RBM15 was shown to mediates region-selective methylation in a WTAP dependent manner, directing specificity towards U-rich sequences.
The observed specificity of the methyltransferase complex to methylate only site specific DRACH sequenced is still poorly understood. Some possible modulators like the role of the accessory proteins are under investigation, however, the structural context of the RNA methylation sites or a structural preference of the complex have been mainly neglected so far. Moreover, the structural dynamics of this methylation process still remain elusive. This thesis contributes to the afore-mentioned aspects by analysis of the methylation process regarding RNA structure sensitivity with enzymatic activity assays and its dynamic nature by implementing a smFRET approach.
We hypothesized the target RNA secondary structure to be an additional important modulator of methylation efficiency, based on the RNA binding elements of the complex (positively charged binding groove, zinc finger domain, RGG domain) and the supposed target adenine binding in the active site. Here, we postulated the possibility for a flipped-out adenine to be of special relevance, which is closely related to the local stability of the target adenine containing structure. Moreover, efficient binding of the protein complex to the RNA should require the ability to anchor the RNA on both sides of the target sequence.
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.
Diseases such as cardiac arrhythmias, CPVT and other issues of the human heart still remain largely unexplored. To contribute to this field of research, it is necessary to create tools to control the spatial and temporal release and reuptake of Ca2+ from the sarcoplasmic/endoplasmic reticulum (SR/ER). Ca2+ release and uptake by the ryanodine receptor (RyR) and Sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA), respectively, are essential for the function of excitable cells. In this process, the rapid Ca2+ release from the SR/ER and the associated contraction in muscle cells is modulated by RyR. However, diseases due to calcium leakage, such as cardiac arrhythmias, seizures and contractile dysfunction, are also caused by RyR. The resting Ca2+ concentration in the cytosol, which is important for the cell, is kept in balance by Ca2+ release and reuptake into the SR/ER. This reuptake is controlled quite considerably by SERCA. SERCA is important for development and muscle function in both nematodes such as C. elegans and mammals, though there is also a great need for tools that can help study precise function.
To advance towards the goal of developing tools for optogenetic stimulation of intracellular Ca2+ release from the SR/ER, the model organism C. elegans was chosen. Its advantages are the fully sequenced genome and the neural network connectome. In addition, the ease of maintenance, self-fertilisation, transparency and rapid generation cycles, as well as the fact that it is a eutelic animal, are advantages for the application of the optogenetic approach.
So far, tools for light-induced Ca2+ release (LICR) have already been developed, involving the creation of ChR2 versions with higher Ca2+ conductivity based on the "CatCh" variant and further improving their conductivity through several established mutations. In addition, the pharynx of C. elegans was modified to produce an optogenetically stimulated muscle pump that resembles mammalian cardiac muscle cells. In this work, both optoUNC-68 (optically excitable RyR) and SERCA/LOV2 were generated in different variants by CRISPR/Cas9 and plasmid-based genome editing to achieve light-driven manipulation of calcium homeostasis in C. elegans. Here, LICR was triggered by LOV2 domains in an opto-mechanical manipulation of RyR as well as SERCA. This approach was made possible by recently published high-resolution cryoEM structural images. In addition, alternative approaches using Ca2+ conductance-optimised channelrhodopsin variants were tested in C. elegans body wall muscle cells.
By inserting ChR-XXM into C. elegans and subsequent fluorescence microscopy of the co-introduced GFP, an expression in body wall muscle cells could be detected. Furthermore, in contraction assays, ChR-XXM was demonstrated to induce contractions of the animals of up to 16% compared to the original body length in both medium (0.8mW/mm²) and high (1.4mW/mm²) stimulation at 470nm. ChR-XXM was thus identified as an excellent candidate for the development of an optogenetic tool, as it exhibits significantly increased Ca2+ conductivity compared to other ChR2 variants.
The use of CRISPR/Cas9 to insert AsLOV2 domains (L404-L546) into different insertion sites of RyR allowed the generation of a transgenic strain of C. elegans that could be stimulated to elongate during 0.3mW/mm² photostimulation. This demonstrated that RyR can be manipulated by photostimulation, spatiotemporally through conformational changes in the LOV2 domain and the resulting disruption of the pore region.
The CRISPR/Cas9 method was also used to insert LOV2 domains into SERCA. Here it could be demonstrated that a conformational change of the LOV2 domains induced by photostimulation leads to a stop or impairment of Ca2+ ion translocation by SERCA from the cytosol into the SR/ER. In contrast to LOV2 in RyR, this resulted in a contraction of C. elegans body length.
The data presented here indicate that the intracellular Ca2+ cycle involving the SR/ER and cytosol can be successfully manipulated by the introduction of optogenetic tools. It turned out that the manipulation/impairment of individual components of this system, such as RyR or SERCA, is usually insufficient to achieve a clear response. Therefore, simultaneous manipulation of the two main actors RyR and SERCA is arguably the best way to take another step towards creating optogenetic tools for light-stimulated manipulation of Ca2+ release and reuptake from the SR/ER.