Post on 13-Feb-2016
description
ZellsimulationenNerven-Verbindung (synaptischer Spalt).Nach seiner Auschüttung bindet der Neurotransmitter Acetylcholin (hellblaue Kugeln) an die Acetylcholinrezeptoren (tassenförmige Objekte) und an Acetylcholinesterase. Doppelt besetzte Rezeptoren (gelb) leiten Strom, der eine Kaskade von Vorgängen einleitet, die zur Kontraktion des Muskelstrangs führen.
Diese und nächste Stunde werden 3 Simulationspakete behandelt, E-cell, Virtual Cell, Mcell, die 3 verschiedene Paradigmen für Zellsimulationen verkörpern
- ODE = gewöhnliche Differentialgleichungen, enthalten ∂/∂tgrösseres - PDE = partielle Differentialgleichungen, enthalten ∂/∂t und ∂/∂rDetail - explizite Simulation der Brownschen Bewegungen einzelner
Moleküle- Projekte unserer Arbeitsgruppe
Virtual Cell: Software-Umgebung für computerunterstützte Zellbiologie
Prof. Leslie Loew, University of Conneticut Health CenterNational Resource for Cell Analysis and Modeling
Virtual Cell
- wurde für Zellbiologie-Community entwickelt
- ermöglicht die Konstruktion räumlicher Modelle
- Verbindung zur quantiativen Lichtmikroskopie an lebenden Zellen
- Kann man auf der Basis des komplexen räumlichen und zeitlichen Zusammen-spiels der Zellkomponenten ein quantitatives Verständnis des gesamten zellulären Verhaltens entwickeln?
- Sind die identifizierten Komponenten notwendig und hinreichend?
- Wie sensitiv reagiert der Gesamtprozess auf Veränderungen einer Komponente?
(Zellen sind „robust“).
- Die Simulationen werden über das Internet auf einem 16-Prozessor cluster mit
Alpha-Prozessoren durchgeführt.
Design einer Virtuellen Zelle
Die `Physiologie' beinhaltet die topologische Anordnung von Kompartments und Membranen, die mit ihnen assoziierten Moleküle, und die Reaktionen zwischen den Molekülen.
Die getrennt definierte `Geometrie‘ ist eine räumliche Beschreibung der Kompartments in 0-3 Dimensionen.
Sie kann aus analytischen Ausdrücken bestehen oder aus einem experimentellen Bild abgeleitet werden, das z.B. von einem Mikroskop stammt.
Das eigentliche Modell besteht aus der Verbindung der Topologie der physiologischen Beschreibung mit einer geeigneten räumlichen Geometrie.
Virtual cell: graphische Benutzerschnittstelle (GUI)Das GUI von Virtual Cell ist als JAVA applet innerhalb eines Webbrowsers entwickelt.
Hier sieht man, wie eine Zelltopologie einer bestimmten experimentellen Geometrie zugeordnet wird.
Einzelschritte bei Erstellung von BioModellenStructure Mapping – definiert die Beziehung zwischen der Physiologie (zelluläre Strukturen) und der Geometrie des Modells. Bestimme das Verhältnis von Oberfläche zu Volumen für die Modelle der Kompartments oder für nicht aufgelöste räumliche Strukturen. Bilde die zellulären Strukturen auf geometrische Objekte ab.
Wähle zwischen unterschiedlichen Randbedingungen (Wert bzw. Ableitung am Rand = Dirichlet bzw. Neumann) für die Strukturen.
Anfangsbedingungen – Konzentrationen und Diffusionsraten können räumlich variable definiert werden. Wähle Anfangsbedingungen für Diffusion ≠ 0.
Reaction Mapping – erlaube oder verbiete Reaktionen bzw. Flüsse.
Math Viewer – prüfe die mathematische Beschreibung, die vom Programm automatisch für die Abbildung des physiologischen Modells auf ein Kompartment-Modell oder auf ein räumliches Modell erstellt wird.
Abstraktion experimenteller Geometrien
Zuordnung von Physiologie und Geometrie
Definition des Flusses durch Membranen
Zwischen verschiedenen Kompartments erzeugt das Programm automatisch Membranen.
Flüsse beziehen sich auf eine bestimmte Membran und beschreiben den Fluss eines bestimmten Moleküls durch diese Membran.
Das Molekül muss auf beiden Seiten der Membran existieren.
Flüsse werden in (mM*mm*s-1) gemes-sen und sind beliebige Funktionen der Transporter- bzw. Kanalkonzentrationen.
Definition der mathematischen Formeln
zeitabhängige Darstellung der Ergebnisse
räumliche Darstellung der Ergebnisse
BK-induzierte Calcium-Wellen
Fink, Slepchenko, Moraru, Watras, Schaff, Loew
Biophysical Journal, 79, 163 (2000)
“An Image-Based Model of Calcium Waves in Differentiated Neuroblastoma Cells”
Die intrazelluläre Calcium-Dynamik ist ein besonders geeignetes Modellsystem. Es existieren viele experimentelle Techniken um die molekularen Komponenten zu visualisieren, die die Calcium-Ausschüttung bestimmen, und um die Calcium-Konzentration selbst mit Fluoreszenz-Farbstoffen zu visualisieren.
Pathway für die Bradykinin-induzierte Ausschüttung von Calcium in differentiereten Neuroblastoma-Zellen
Sobald Bradykinin (BK) an seinen Rezeptor in der Plasmamembran bindet (BKR), setzt dies eine G-Protein Kaskade in Gang (gezeigt als αGq, βγ), die Phospholipase C (PLC) aktiviert.
PLC hydrolysiert wiederum Glycerinphosphat, das in Phosphatidylinositolbisphosphat (PIP2) gebunden ist. Das Produkt dieser Hydrolyse setzt IP3 von der Membran frei. IP3 kann frei im Zytosol diffundieren, dabei von Phosphatasen und Kinasen abgebaut werden, oder an seinen Rezeptor IP3R in der ER-Membran binden. IP3R ist ein Calcium-Kanal, der sich nach Bindung von IP3 und Ca2+ öffnet. Das so aus dem ER freigesetzte Ca2+ bindet an Ca2+ Buffer (B) im Zytosol. Schliesslich wird Ca2+ durch die SERCA Ca2+-Pumpe ins ER zurückgepumpt, bzw. durch andere Calcium-pumpen und –austauscher aus der Zelle hinaus (Caext.).
BK
IP3
Experiment: BK-induzierte Calcium-Welle
BK (500 nM) wurde extern zum Medium einer N1E-115 Nb Zelle zugegeben, die mit dem Farbstoff fura-2 gesättigt war.
Die Fluoreszenz einer 380-nm Anregung wurde mit einer CCD Kamera gemessen, in die [Ca2+] Konzentration umgewandelt, und in einer Falschfarbendarstellung gezeigt.
Relative zeitliche Änderung von [Ca2+] in zwei
interessanten Regionen im Zellsoma und in dem Neuron (im obersten Bild durch Rechtecke gezeigt):
Welle kommt etwas später im Soma an,
besitzt überall gleiche Amplitude.
Calcium-Welle (exp)
QuickTime™ and aYUV420 codec decompressor
are needed to see this picture.
Es gibt mehr ER im soma als in neuronaler Fortsetzung
(a) Electronmikroskopische Aufnahme des ERs im Soma und im Neurite von N1E-115 Neuroblastoma-Zellen.
(b) Das ER einer lebenden Zelle wird einige Minuten nach Injektion durch einen fluorezenten, lipophilen Farbstoff sichtbar gemacht.
Verteilung von InsP3R, BKR und SERCA in ER- Membran
(a) Visualisierung intrazellulärer Verteilungen mit 3D Immunofluoreszenz. Hier
wurde N1E-115 Nb-Zelle fixiert und mit einem Antikörper für SERCA2
markiert.
(b) Intrazelluläre Verteilung von BKR und ER/InsP3-R/SERCA2.
Die relative Konzentration von BKR in der Plasmamembran (äussere Schicht)
wird gemittelt für sechs verschiedene Regionen der Zelle angegeben.
Die mittlere Verteilung des ER
ist identisch mit der von InsP3R
and SERCA und wird innerhalb
der Zelle durch Graustufen
angegeben (siehe Skala rechts).
Simulationsdetails
InsP3-Produktion
InsP3-Diffusion und -Abbau
Calcium-Dynamik
Kinetik der InsP3-Kanäle
Kinetik der SERCA-Pumpe
und: Calcium-Speicherung (nicht gezeigt)
instantaner Anstieg nach BK-Reizung
Simulation von BK-induzierter Calcium-Welle
In den Spalten von links nach rechts zeigen die pseudokolorierten Bilder die Simulationsergebnisse für [Ca2+]cyt, [InsP3]cyt, und PO (die W’kt dass der
IP3-Rezeptor offen ist).
Unten sieht man den zeitlichen Verlauf in Soma und Neurite (deren Positionen sind als gelbe und grüne Punkte im Bild links oben gezeigt).
Experimentelle Werte
Calcium wave (simulation)
Ca2+ – ist exp.messbar
InsP3 – nicht
exp. messbar
Simulation in 2D dauert etwa 15 Minuten auf PC. Δt = 0.01 s.
QuickTime™ and aYUV420 codec decompressor
are needed to see this picture.Text
unterschiedliche Zellgeometrien
Experimentelle und simulierte BK-induzierte Calcium-Ausschüttung in N1E-115 Nb Zellen unterschiedlicher Morphologie.
(a) Modellierung wie zuvor, aber mit 2 Neuriten anstelle von einem.
(b) Es wurden konstante Verteilungen von BKR und ER in der Zelle angenommen.
Skalierung (a) 50 µM und (b) 20 µM.
Virtual cell: Ausblick
aktuelle Version: - ermöglicht Simulation von Reaktions-Diffusionsprozessen in beliebigen Geometrien. Anpassung notwendig für Probleme, die Änderungen der Geometrie erfordern (Zellwanderung, Zellteilung).
-behandelt nur bestimmte Sorten von stochastischen Prozessen:
Brownsche Bewegung, gerichtete Teilchenbewegung entlang von Mikrotubuli,
Reaktion einzelner Teilchen mit kontinuierlich verteilten Molekülen.
- wenn die Anzahl an wechselwirkenden Molekülen zu klein wird, braucht man statt der stochastischen Beschreibung Reaktionen zwischen diskreten Molekülen.
- Behandlung diskreter Zustände ist auch erforderlich zur Modellierung der Ströme von einzelnen Ionenkanälen und deren räumlicher Verteilung.
Virtual cell: Ausblick II
Es wird erforderlich (und möglich) sein, Simulationsmodelle zwischen verschiedenen Softwareplattformen auszutauschen und eventuell gleichzeitig auszuführen.
z.B. benutze GEPASI bzw. E-Cell um ein Modell auf der Ebene von ODEs zu optimieren und reduzieren bevor man mit Virtual Cell die räumlichen Verteilungen mittels PDEs simuliert.
Zwei Ansätze hierfür:
CellML – markup language von Physiome Sciences Inc. & Univ. of Auckland
SBML – Caltech & ERATO-Kitano Symbiotic Systems Project
M-cell: Allgemeine Monte–Carlo Simulation vonzellulären Mikrophysiologien
Thomas M. Bartol Jr. Joel R. Stiles
Computational Neurobiology Laboratory Biomedical Applications
(T. Sejnowski), Salk Institute, San Diego Pittsburgh Supercomputing Center
Ziel: quantitatives, molekulares Verständnis der Nervenfortleitung,Funktion von Nervengasen, Modulatoren, oder Autoimmunerkrankungen.
MCell: Idee + Motivation
MCell ermöglicht 3-D Monte–Carlo Simulationen für Ligandendiffusion und chemische Signalprozesse.
Biologische Strukturen wie Neuronen zeigen auf der subzellulären Ebene eine enorme Komplexität und Diversität. Die inter- und intrazelluläre Kommunikation geschieht mittels verschiedener chemischer Signalpfade.
Am Prozess der synaptischen Transmission sind z.B. Neurotransmitter und Neuromodulatoren beteiligt. Ebenfalls beteiligt sind Proteine, die die Auffüllung und Entleerung der synaptischen Vesikel mit Neurotransmitter-Molekülen beeinflussen, Rezeptorproteine, Transportproteine, sowie oxidierende und hydrolytische Enzyme.
Mit Mcell kann man alle diese Parameter in beliebig komplexen räumlichen Darstellungen der beteiligten zellulären Strukturen darstellen und variieren.
Anfangsbedingung: Eine Monte–Carlo Simulation beginnt damit, dass die Zellumgebung mit einzelnen Liganden und Liganden-bindenden Molekülen angefüllt wird.
Warum soll man Monte Carlo Algorithmen benutzen? 1. löse PDEs zwischen Voxels
Die Diffusion von Ligandenmolekülen in Lösung basiert auf Brownscher Bewegung.
Der mittlere Netto-Fluss aus einer Region des Raums in eine andere hängt von der Mobilität der Moleküle und dem räumlichen Konzentrationsunterschied der beiden Regionen ab.
Eine Methode, den räumlichen Gradienten zu berechnen, ist, den Raum in kleine, üblicherweise kubische Volumenelemente (Voxels) aufzuteilen, innerhalb derer man gute Mischung annimmt, und dann mittels eindimensionaler partieller Differentialgleichungen den mittleren Fluss durch die Verbindungsfläche zwischen angrenzenden Voxeln zu berechnen.
Sofern die Granularität der räumlichen und zeitlichen Unterteilung fein genug ist, wird eine numerische Simulation das korrekte mittlere Verhalten des Systems erzeugen.
löse PDEs innerhalb von Voxels
Man kann weitere PDEs hinzufügen um die mittlere Raten chemischer Raten Reaktionen innerhalb jedes Voxels zu beschreiben. Man erhält damit eine Simulation des räumlichen und zeitlichen Diffusionsverhaltens und der chemischen Reaktionen.
Für einfache räumliche Anordnungen kann diese Methode sehr effizient sein.
Für komplexe (d.h. realistische) Structuren werden die räumlichen Unterteilungen immer komplexer und eine grosse Anzahl an Voxeln ist erforderlich.
Auf jeden Fall liefert die Simulation keine direkten Informationen über die stochastischen Schwankungen, die auf der endlichen Anzahl an beteiligten Molekülen beruhen. Diese sind in biologischen Systemen jedoch oft von grossem Interesse.
2. völlig andere Methode: Random walk
Direkte Beschreibung der Brownschen Bewegung der einzelnen Ligandenmoleküle. Durch Verwendung von Zufallszahlen werden bei jedem Zeitschritt beliebige erlaubte Richtungen und Verschiebungen ausgewählt. Indem die Zeitschritte und Verschiebungen deutlich kleiner als die Teilchengröße gehalten werden, erreicht man eine hohe numerische Genauigkeit.
Kollisionen mit beliebigen Oberflächen werden entdeckt und gemäss von Regeln behandelt (Bindung, Transport, Reflexion etc.). Voxel sind unnötig.
Gleichsam werden Kollisionen mit möglichen Bindungsstellen entdeckt und behandelt. Für die Ausbildung von Bindungen werden Bindungswahrscheinlich-keiten festgelegt. Die momentane Entscheidung wird durch eine Zufallszahl bestimmt.
Alle möglichen Vorgänge werden auf einer Molekül-für-Molekül Basis betrachtet. Dadurch enthält die Simulation realistische stochastische Schwankungen in Abhängigkeit von der räumlichen Verteilung und der Anzahl an Molekülen.
Rechenaufwand: sehr hoch
konkrete Anwendung, Effekt eines Nervengases15,680 Simulationen für jeden Zustand.
Jede Simulation entspricht 40 Millisekunden in Realzeit und dauert je nach Parametern zwischen 10 Sekunden bis zu 6 Stunden, im Mittel etwa 1 Stunde.
gesammelte Daten: 60 Gigabyte
Die Rechnungen wurden parallel auf bis zu 1,024 Prozessoren von Blue Horizon (SDSC) durchgeführt.
Typische Vorgänge während einer MCell-Simulation
- Freisetzung von Ligandenmolekülen aus einer Struktur (z.B. einem Vesikel)
- Erzeugung oder Vernichtung von Ligandenmolekülen
(z.B. Synthese, Hydrolyse, oder Redox-Reaktionen)
- Diffusion der Liganden innerhalb des Raums zwischen beliebigen Oberflächen
(z.B. prä- und postsynaptische Membranen)
- chemische Reaktionen von Ligandenmolekülen
mit “Effektor”molekülen
(z.B. Rezeptoren oder Enzyme)
MCell: biologische Skala
Der Level an Detail von MCell Simulationen liegt zwischen denen von atomistischen Moleküldynamik-Simulationen und der Simulationen gesamter Zellen (Virtual Cell, E-cell).
Die Diffusion einzelner Ligandenmoleküle wird als Brownsche Bewegung mit einem Random–Walk–Algorithmus simuliert.
Mittlere Ratenkonstanten werden in Monte– Carlo–Wahrscheinlichkeiten für Reaktionen einzelner Moleküle pro Zeitschritt umgeformt.
Damit können die Ligandenmoleküle stochastisch reagieren sobald sie an Rezeptoren, Enzyme oder Transporter gebunden sind.
Tutorial: nichtgebundene Diffusion
In diesem einfachen Beispiel werden Liganden an zwei Punkten freigesetzt (rot und blau markiert).
Durch die Verwendung unter-schiedlicher Diffusionskonstanten wird unterschiedliches Verhalten erzeugt.
Tutorial: nichtgebundene Diffusion
file run_parameters /* run_parameters */
dt = 1.0e-6 /* Simulationszeitschritt*/
iterations = 300
viz_output_file = "viz_data“
/* ligand_parameters */
red_ligand_diffusion_constant = 6.0e-6
blue_ligand_diffusion_constant = red_ligand_diffusion_constant/2.0
number_of_ligand_molecules = 5.0e3
release_site_diameter = 0.03
red_location = [0.0, 0.0, 0.0]
blue_location = [-1.0, -1.0, 0.0]
zeitlicher Verlauf
Tutorial 2: Sampling der Liganden
Hier werden die gleichen Punkte zur Ausschüttung der Liganden verwendet wie zuvor.
Zur Analyse werden Reaktions-behälter definiert um die Liganden-moleküle zu zählen, die sich in bestimmten Gebieten befinden und diese Information in eine Datei
REACTION_DATA_OUTPUT zu schreiben.
Tutorial 2: Sampling der Liganden
Tutorial 4: einzelne Bindungsprozesse
Hier wird die Simulation von Reaktionszentren wie Effektoren eingeführt.
In der Box führen Ligandenmoleküle Brownsche Bewegung aus.
Am Boden der Box befindet sich ein Gitter von Effektormolekülen.
Die rot markierten Zylinder haben ein Ligandenmolekül gebunden.
Tutorial 6: AntiporterHier werden Antiporter-Moleküle in der Membran modelliert, die den gleichzeitigen entgegengesetzen Transport zweier Liganden L und M durch die Membran bewirken.
Blaue Antiporter haben keine Bindungsstelle besetzt, rote und grüne Antiporter haben nur einen roten bzw. grünen Liganden gebunden.
Gelbe Antiporter sind doppelt besetzt. Nun ist Transport möglich.
Iteration 1000 Iteration 5000
Iteration 0
Iteration 15000
Tutorial 9: einfache Pore
Hier wird die Grösse einer Vesikelpore während der Simulation skaliert.
Iteration 0 Iteration 250 Iteration 750
Iteration 1000
räumliche Rekonstruktion aus exp. Abbildungen einer neuromuskulären Verbindung in Maus–Sternomastoid
Simulation der in den Spalten erzeugten Ströme
Ausschnitt aus Bild auf vorheriger Seite.
Diffusion von Neurotransmittern in Synapsen
Neurotransmitter-Moleküle (Acetylcholin) diffundieren und aktivieren Rezeptoren in den Synapsen zwischen verschiedenen Zellen.Die wie Tassen aussehenden Objekte sind Acetylcholin-Rezeptoren:
blaue sind frei;rote haben 1 Ach; grüne haben 2 Ach gebunden,
sind jedoch geschlossen; gelbe haben 2 Ach gebunden,
sind offen.Die Anzahl an Ionen, die durch die gelben Rezeptoren fliesst, ergibt den elektrischen Strom an der Synapse.
Ausschnitt aus Bild auf vorheriger Seite.
zeitliche Entwicklung der Rezeptorzustände
Stiles et al., 1996, PNAS 93:5747-5752.
22 μs nachdem sich die Pore öffnet 90 μs
160 μs
Zusammenfassung
Zellsimulationen sind im Kommen!
Detaillierte, dreidimensionale Modelle sind notwendig,sobald Lokalisation z.B. an Membran stattfindet, undsobald wichtige Moleküle in kleinen Zahlen vorliegen.
Schlagwort: Systems Biology.
Messung von kinetischen Konstanten im Allgemeinen mühsam.Daher zunächst Konzentration auf Modellsysteme.
Molekulare Simulationen können sehr aufwendig sein.Ein Volumen von 100 nm3 enthält ca. 1000 Proteine.1 μm3 enthält dagegen bereits ca. 106 Proteine.
Bei Beschränkung auf Molekül-Konzentrationen kann dagegen fast in real time simuliert werden.