Switch branches/tags
Nothing to show
Find file
Fetching contributors…
Cannot retrieve contributors at this time
155 lines (139 sloc) 12.9 KB
<!DOCTYPE html>
<html xmlns="">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8"/>
<title>GFRD background</title>
<link rel="stylesheet" href="/css/reset.css" type="text/css"/>
<link rel="stylesheet" href="/css/style.css" type="text/css"/>
<link rel="stylesheet" href="/css/syntax.css" type="text/css"/>
<link rel="stylesheet" href="/css/menu.css" type="text/css"/>
<link rel="stylesheet" href="/css/footer.css" type="text/css"/>
<script src="/html5media/html5media.min.js"></script>
<!--[if lt IE 9]>
<script src="/ie7/IE9.js"></script>
<div class="background">
<div id="globalheader">
<!--googleoff: all-->
<ul id="globalnav">
<li id="gn-home"><a href="/">Home</a></li>
<li id="gn-algorithm"><a href="/algorithm">Algorithm</a></li>
<li id="gn-background"><a href="/background">Background</a></li>
<li id="gn-downloads"><a href="/downloads">Downloads</a></li>
<li id="gn-documentation"><a href="/documentation">Documentation</a></li>
<li id="gn-support"><a href="/support">Support</a></li>
<li class="gn-empty"><a></a></li>
<!--googleon: all-->
<div id="lastbutton"> </div>
<h1>Background information</h1>
<div id="container">
<h2>The importance of modeling networks in time and space at the particle level</h2>
<p>In many networks of interacting components, the spatial distribution of the
components and the stochastic character of their interactions are of critical
importance for the macroscopic behavior of the system. Examples are to be
found in biological evolution and population dynamics, but arguably the most
studied and best characterized examples are biochemical networks. Biochemical
networks consist of biomolecules, such as proteins and DNA, which chemically
and physically interact with one another. These biochemical networks are the
analog computers of life. They allow living cells to detect, amplify and
integrate signals, as well as transmit signals from one place to another.
Indeed, biochemical networks can perform a large number of computational tasks
analogous to electronic circuits. However, their design principles are
markedly different. In a biochemical network, the components often move in an
erratic fashion, namely by diffusion, and they interact with each other in a
stochastic manner – chemical and physical interactions are probabilistic in
nature. To understand the design principles that allow biochemical networks to
process information reliably, we have to take into account the spatial
distribution of the components and the stochastic character of their
interactions<sup><span id="r1"><a href="#c1">1</a></span></sup>. Moreover, we have recently shown that
spatio-temporal correlations at molecular length scales can drastically change
the macroscopic behavior of systems that are <em>uniform</em> at cellular length
scales<sup><span id="r2"><a href="#c2">2</a></span></sup>. Finally, many biological systems exhibit
macroscopic concentration gradients. Well-known examples are morphogen
gradients, which provide positional information to differentiating cells
during embryonic development. To understand how these gradients are read out
precisely<sup><span id="r3"><a href="#c3">3</a></span></sup>, the system has to be described at the particle
level in time and space<sup><span id="r4"><a href="#c4">4</a></span></sup>.</p>
<h2>Overview of simulation techniques</h2>
<p>In principle, computer simulations are ideally suited for studying
reaction-diffusion systems. The conventional approach is to write down the
macroscopic reaction-diffusion equations and solve the corresponding
differential equations either analytically or numerically. This approach
resolves the system in time and space, but ignores the discrete nature of the
reactants and the stochastic character of their interactions. Schemes based on
the zero-dimensional chemical master equation, such as the Gillespie
algorithm<sup><span id="r5"><a href="#c5">5</a></span></sup> or the Gibson-Bruck scheme<sup><span id="r6"><a href="#c6">6</a></span></sup>,
incorporate the latter, but assume that at each instant the particles are
uniformly distributed in space.</p>
<p>In the past years, a number of techniques have been developed that take into
account both the spatial distribution of the components and the stochastic
nature of their interactions. One class of techniques is based upon the
reaction-diffusion master equation<sup><span id="r7"><a href="#c7">7</a></span>,<span id="r8"><a href="#c8">8</a></span></sup>;
examples are SmartCell<sup><span id="r9"><a href="#c9">9</a></span></sup>, MesoRD<sup><span id="r10"><a href="#c10">10</a></span></sup> and
URMD<sup><span id="r11"><a href="#c11">11</a></span></sup>. The main idea of these techniques is to divide the
reaction volume into a number of subvolumes; particles can react within the
subvolumes, but also diffuse from one subvolume to the next. Importantly, it
is assumed that within each subvolume the particles are well mixed. These
techniques thus rely on the existence of a length and time scale on which the
system is spatially uniform. When the concentrations are fairly high, as in
spatio-temporal calcium oscillations, such a scale possibly exists, but at
lower concentrations this assumption is likely to fail.</p>
<p>Another class of techniques simulates the network in time and space at the
particle level. One natural approach is to use Brownian Dynamics; examples are
Smoldyn<sup><span id="r12"><a href="#c12">12</a></span></sup>, MCell<sup><span id="r13"><a href="#c13">13</a></span></sup>, Reaction Brownian
Dynamics<sup><span id="r14"><a href="#c14">14</a></span></sup> and GridCell<sup><span id="r15"><a href="#c15">15</a></span></sup>. In essence,
the particles undergo a random walk, and when two reaction partners happen to
meet each other, they can react with a probability consistent with the rate
constant. However, under biologically relevant conditions, Brownian Dynamics
is not very efficient, since much CPU time is wasted on propagating the
particles towards each other. Moreover, it is not exact, since a finite time
step is used.</p>
<p>Another approach to simulate a network at the particle level in time and space
is to use an event-driven algorithm. Green’s Function Reaction Dynamics
(GFRD)<sup><span id="r16"><a href="#c16">16</a></span>,<span id="r17"><a href="#c17">17</a></span>,<a href="#c2">2</a></sup> and First Passage
Kinetic Monte Carlo (FPKMC)<sup><span id="r18"><a href="#c18">18</a></span>,<span id="r19"><a href="#c19">19</a></span></sup> are
examples of such an approach. The idea is to exploit the solution of the
Smoluchowski equation — the Green’s function — to set up an event-driven
algorithm that concatenates the propagation of the particles in space with the
reactions between them. These algorithms thus alleviate the need to propagate
the particles toward each other to let them react, as in conventional Brownian
Dynamics: even when the reactants are far apart from one another, the
algorithm can immediately jump to the next reaction event. This event-driven
nature makes these schemes highly powerful, especially when the concentrations
of the components are low, as in most signal transduction pathways and gene
regulation networks.</p>
<div class="references"><h3>References</h3>
<p><a href="#r1"><span id="c1">1</span></a> van Zon JS, Morelli M, Tanase-Nicola S, ten Wolde PR (2006) Diffusion of transcription factors can drastically enhance the noise in gene expression. <em>Biophys J</em>, 91:4350 — 4367. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r2"><span id="c2">2</span></a> Takahashi K, Tanase-Nicola S, Ten Wolde PR (2010), Spatio-temporal correlations can drastically change the response of a MAPK pathway. <em>Proc Natl. Acad Sci USA</em>, 107:2473 — 2478. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r3"><span id="c3">3</span></a> Gregor T, Tank DW, Wieschaus EF, Bialek W (2007) Probing the limits to positional information. <em>Cell</em>, 130: 153 — 164. (<a href="">doi</a>)<br/>
<a href="#r4"><span id="c4">4</span></a> Erdmann T, Howard M, Ten Wolde PR (2009) Role of spatial averaging in the precision of gene expression patterns. <em>Phys Rev Lett</em>, 103: 258101. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r5"><span id="c5">5</span></a> Gillespie DT (1977) Exact stochastic simulation of coupled chemical reactions. <em>J Phys Chem</em>, 81:2340 — 2361. (<a href="">doi</a>)<br/>
<a href="#r6"><span id="c6">6</span></a> Gibson MA, Bruck J (2000) Efficient exact stochastic simulation of chemical systems with many species and many channels. <em>J Phys Chem A</em>, 104: 1876 — 1888. (<a href="">doi</a>)<br/>
<a href="#r7"><span id="c7">7</span></a> Baras F, Mansour M (1996) Reaction-diffusion master equation: A comparison with microscopic simulations. <em>Phys Rev E</em>, 54:6139 — 6148. (<a href="">doi</a>)<br/>
<a href="#r8"><span id="c8">8</span></a> Stundzia AB, Lumsden CJ (1996) Stochastic simulation of coupled reaction-diffusion processes. <em>J Comp Phys</em>, 127: 196 — 207. (<a href="">doi</a>)<br/>
<a href="#r9"><span id="c9">9</span></a> Ander M et al. (2004) SmartCell, a framework to simulate cellular processes that combines stochastic approximation with diffusion and localization: analysis of simple networks. <em>Sys Biol</em>, 1:129 — 138. (<a href="">doi</a>)<br/>
<a href="#r10"><span id="c10">10</span></a> Hattne J, Fange D, Elf J (2005) Stochastic reaction-diffusion simulation with MesoRD. <em>Bioinformatics</em>, 21: 2923 — 2924. (<a href="">doi</a>)<br/>
<a href="#r11"><span id="c11">11</span></a> Engblom S, Ferm L, Hellander A, Lötstedt P (2009) Simulation of stochastic reaction-diffusion processes on unstructured meshes. <em>SIAM J Sci Comput</em>, 31:1774 — 1797. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r12"><span id="c12">12</span></a> Andrews S, Bray D (2004) Stochastic simulation of chemical reactions with spatial resolution and single molecule detail. <em>Phys Biol</em>, 1:137 — 151. (<a href="">doi</a>)<br/>
<a href="#r13"><span id="c13">13</span></a> Coggan JS, Bartol TM, Esquenazi E, Stiles JR, Lamont S, Martone ME, Berg DK, Ellisman MH, Sejnowski TJ (2005) Evidence for ectopic neurotransmission at a neauronal synapse. <em>Science</em>, 309:446 &mdash; 451. (<a href="">doi</a>)<br/>
<a href="#r14"><span id="c14">14</span></a> Morelli MJ, Ten Wolde PR (2008) Reaction Brownian Dynamics and the effect of spatial fluctuations on the gain of a push-pull network. <em>J Chem Phys</em>, 129: 054112. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r15"><span id="c15">15</span></a> Boulianne L, Al Assaad S, Dumontier M, Gross WJ (2008) GridCell: a stochastic particle-based biological system simulator. <em>BMC Sys Biol</em>, 2:66. (<a href="">doi</a>)<br/>
<a href="#r16"><span id="c16">16</span></a> van Zon JS, tenWolde PR (2005) Simulating biochemical networks at the particle level and in time and space: Green’s function reaction dynamics. <em>Phys Rev Lett</em>, 94:128103. (<a href="">doi</a>)<br/>
<a href="#r17"><span id="c17">17</span></a> van Zon JS, ten Wolde PR (2005) Green’s-function reaction dynamics: A particle-based approach for simulating biochemical networks in time and space. <em>J Chem Phys</em>, 123:234910. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r18"><span id="c18">18</span></a> Opplestrup T, Bulatov VV, Gilmer GH, Kalos MH, Sadigh B (2006) First-passage Monte Carlo algorithm: diffusion without all the hops. <em>Phys Rev Lett</em>, 97:230602. (<a href="">doi</a>, <a href="">arXiv</a>)<br/>
<a href="#r19"><span id="c19">19</span></a> Opplestrup T, Bulatov VV, Donev A, Kalos MH, Gilmer GH, Sadigh B (2009) First passage kinetic Monte Carlo. <em>Phys Rev E</em>, 80: 066701. (<a href="">doi</a>, <a href="">arXiv</a>)</p></div>
<div id="globalfooter">
<p>Copyright &copy; 2010 FOM institute AMOLF</p>
<ul class="piped">
<li><a href="/about">About</a></li>
<li><a href="/credits">Credits</a></li>