how to build a synapse with mcell/cellblender: an ......2. molecular weight of water 㲓 18g/mole....
TRANSCRIPT
![Page 1: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/1.jpg)
How to Build a Synapsewith
MCell/CellBlender:An Environment for Spatially RealisticSimulation of Cellular Microphysiology
Thomas M. Bartol
Computational Neurobiology LaboratoryThe Salk Institute
MMBioS Cell Modeling WorkshopJune1-3, 2016
Tuesday, June 7, 16
![Page 2: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/2.jpg)
Scientific Computing EnablesThe Scientific Discovery Cycle
“What I cannot create, I do not understand” -- Richard Feynman
Philosophical Approach: Understanding through simulation
Tuesday, June 7, 16
![Page 3: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/3.jpg)
How do synapses work?How does synaptic structure influence synaptic function?
Motivation
•Topology of synaptic and perisynaptic space (neuropil !!!)•Organization of pre- and postsynaptic cytoplasmic organelles•Distributions of receptors, ion channels, enzymes, transporters...•Biochemical reaction networks and their dynamics•Cells are not well mixed & numbers of molecules can be small• Simulations must include stochastic spatiotemporal dynamics
Tuesday, June 7, 16
![Page 4: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/4.jpg)
Cells are Not Well Mixed Bags of Cytoplasm
Biochemistry should be explored using 3D reaction/diffusion simulation to put reaction networks in their natural context.
Tuesday, June 7, 16
![Page 5: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/5.jpg)
MCell: Monte Carlo Simulator of Cellular Microphysiology
•Models realistic 3D reaction/diffusion systems•Rigorously validated and highly optimized stochastic Monte Carlo
methods•Tracks Brownian dynamics diffusion and reaction of individual
particles in 3D volumes and on 2D surfaces embedded in 3D•Arbitrarily complex 3D geometry -- triangle surface meshes•Arbitrarily complex reaction networks -- Markov processes and
Network-free rule-based specification modeled as discrete event-driven point processes
•High-level, flexible Model Description Language•New model building, visualization and analysis environment --
CellBlender
Tuesday, June 7, 16
![Page 6: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/6.jpg)
Realistic Model of Synaptic Signaling in NeuropilRealistic synaptic morphology from ssTEM reconstruction of hippocampal neuropil(collaboration with Kristen Harris, UT Austin)
Tuesday, June 7, 16
![Page 7: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/7.jpg)
Tuesday, June 7, 16
![Page 8: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/8.jpg)
Tuesday, June 7, 16
![Page 9: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/9.jpg)
Tuesday, June 7, 16
![Page 10: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/10.jpg)
Tuesday, June 7, 16
![Page 11: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/11.jpg)
Tuesday, June 7, 16
![Page 12: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/12.jpg)
• Molecular: molecular modeling and simulations, with focus on identifying functional substates and interaction mechanisms for proteins and their complexes (ProDy, GTKDynamo, WESTPA)
• Cellular: cell modeling, with an emphasis on developing tools to handle spatial and molecular complexity inherent to neuronal signal transmission (MCell, CellBlender, BioNetGen)
• Tissue: Image processing and analysis, with an emphasis on analysis of cell and tissue organization in support of modeling (AlignTK, CellOrganizer)
Our 5 year mission: Technology Research and Development to link (at least) 3 Scales:
Tuesday, June 7, 16
![Page 13: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/13.jpg)
Challenges in Spatially Realistic Cell Modeling
•Creation of computational quality 3D models is too hard•Managing the complexity of reaction pathways is too hard•Managing simulations, data, analysis, and visualization in high-
dimensional parameter space is too hard
PS: Did I forget to mention, it’s too HARD
Tuesday, June 7, 16
![Page 14: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/14.jpg)
SBML+SBML-spatial
SBM
L
Experiments &
Literature
BioNetGen
BNG
Spa
tial
Para
met
ers
Spatially dependent response
BNG
Spa
tial
Para
met
ers
Step 1: Systems modeling
Step 2: Generative spatial modeling
Step 3: Spatially resolved simulationsGeometry-specific
Initializations (SBML)
Devin SullivanTuesday, June 7, 16
![Page 15: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/15.jpg)
CellBlender Addon for Blender•A complete IDE for realistic reaction-diffusion modeling•MCell acts as external Physics Engine for data-driven visualization •We envision/hope that others in community will plugin as well
Tuesday, June 7, 16
![Page 16: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/16.jpg)
CellBlender Capabilities
All of Blender’s features plus:
•Model Specification•model parameters•molecules• reactions•output options• visualization options• run-time options
• Import geometry from many sources:•CellOrganizer•VolRoverN
• Import whole models from:• SBML Spatial•CellOrganizer as SBML Spatial•BioNetGen as reaction network (as rules, coming soon)
Tuesday, June 7, 16
![Page 17: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/17.jpg)
Let’s Build the Model:Model Constraints
Realistic synaptic morphology from ssTEM reconstruction of hippocampal neuropil(collaboration with Kristen Harris, UT Austin)
Tuesday, June 7, 16
![Page 18: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/18.jpg)
Model ConstraintsExperimental estimates of reaction kinetics and subcellular distributions for important molecules
Tuesday, June 7, 16
![Page 19: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/19.jpg)
Model Constraints
Experimental observations of evoked vesicular release:
Goda & Stevens, PNAS,1994
Dobrunz, Huang & Stevens, PNAS,1997
single stimulus paired-pulse stimulus
Tuesday, June 7, 16
![Page 20: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/20.jpg)
Model Constraint: kinetics of calcium sensor in SNARE complex
Tuesday, June 7, 16
![Page 21: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/21.jpg)
Diffusion: Brownian Motion
Thermal Velocity:
For Water:
Tuesday, June 7, 16
![Page 22: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/22.jpg)
But, how far do the water molecules go between collisions?
1. Density of water ≅ 1g/cc or 1000g/liter.
2. Molecular weight of water ≅ 18g/mole.
3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter4. 55.56 x 6.022 x 1023 molecules/mole = 3.3458 * 1025 or 3.3458 * 1022 molecules/cc5. Average volume of each molecule is: 1/ 3.3458 * 1022 ≅ 3.0 * 10-23 cc/molecule
6. Assuming that each molecule corresponds to a spherical space:
Diffusion: Brownian Motion
in 0.3 ps !!!
Tuesday, June 7, 16
![Page 23: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/23.jpg)
Diffusion: Fick’s Second Law
Applies when concentration in volume is changing in time (i.e. Jin ≠ Jout)
Tuesday, June 7, 16
![Page 24: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/24.jpg)
Free Diffusion: Concentration in Space and Time
Tuesday, June 7, 16
![Page 25: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/25.jpg)
Diffusion: Random Directions
Tuesday, June 7, 16
![Page 26: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/26.jpg)
Grid-Free Random Walk Diffusion
Tuesday, June 7, 16
![Page 27: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/27.jpg)
Surface Mesh Representation
Tuesday, June 7, 16
![Page 28: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/28.jpg)
Diffusion with Boundaries
Tuesday, June 7, 16
![Page 29: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/29.jpg)
Unimolecular Reactions
Distribution of first-order decay lifetimes is:
Tuesday, June 7, 16
![Page 30: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/30.jpg)
Bimolecular Reactions: Rate of Encounter
Tuesday, June 7, 16
![Page 31: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/31.jpg)
Bimolecular Reactions: Volume/Surface
Probability of reaction between diffusing volume molecule and a surface molecule:
From rate of encounter: From Mass Action:
Tuesday, June 7, 16
![Page 32: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/32.jpg)
Bimolecular Reactions: Collision Detection
Tuesday, June 7, 16
![Page 33: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/33.jpg)
Segmentation, Mesh Generation, and Mesh Improvement
• RECONSTRUCT• VolRoverN• GAMer
Tuesday, June 7, 16
![Page 34: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/34.jpg)
Mesh Improvement
Tuesday, June 7, 16
![Page 35: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/35.jpg)
Tuesday, June 7, 16
![Page 36: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/36.jpg)
Future Directions
•Advanced simulation control•parameter sweep•parameter fitting/optimization/estimation
• Integration/coupling of other Physics Engines•Cytomechanics•Electrophysiology/electrodiffusion•Hybrid reaction-diffusion: ODE/PDE/SSA
• Space-filling molecules -- macromolecular crowding, self-assembly of filaments, scaffolds etc...
•Dynamic cell geometry
•Multi-scale simulations spanning molecules to neural circuits -- Blue Brain Project +++
Tuesday, June 7, 16
![Page 37: How to Build a Synapse with MCell/CellBlender: An ......2. Molecular weight of water 㲓 18g/mole. 3. From 1 and 2, the molarity of water = 1000/18 = 55.56 moles/liter 4. 55.56 x 6.022](https://reader031.vdocument.in/reader031/viewer/2022011913/5fb286ae58b49c4eca1c0836/html5/thumbnails/37.jpg)
Acknowledgements
SalkTerry Sejnowski
Bob Kuczewski
Jed Wing
Ben RegnerJay CogganCourtney LopreoreDan KellerJustin KinneyKevin FranksEugene KimElaine Zhang UCSD
Mark EllismanMaryann MartoneJeff BushStephan Lamont
PSCMarkus Dittrich
Jacob Czech
Art Wetzel
CaltechMary Kennedy
Melanie StephanShirley Pepke
UT AustinKristen HarrisChandra Bajaj
NIH & NSF
UPittJames FaederJose-Juan TapiaDipak Barua
CMUBob MurphyDevin Sullivan
Tuesday, June 7, 16