notebook.community

Edit and run



In [1]:

    
import moose
import pylab
import rdesigneur as rd

This function sets up a simple oscillatory chemical system within
the script. The reaction system is::

    s ---a---> a  // s goes to a, catalyzed by a.
    s ---a---> b  // s goes to b, catalyzed by a.
    a ---b---> s  // a goes to s, catalyzed by b.
    b -------> s  // b is degraded irreversibly to s.

in sum, **a** has a positive feedback onto itself and also forms **b**.
**b** has a negative feedback onto **a**.
Finally, the diffusion constant for **a** is 1/10 that of **b**.



In [2]:

    
rdes = rd.rdesigneur(
    turnOffElec = True,
    diffusionLength = 1e-3, # Default diffusion length is 2 microns to make it essentially single compartment
    chemProto = [['makeChemOscillator()', 'osc']],
    chemDistrib = [['osc', 'soma', 'install', '1' ]],
    plotList = [['soma', '1', 'dend/a', 'conc', 'a Conc'],
        ['soma', '1', 'dend/b', 'conc', 'b Conc']]
)



In [3]:

    
rdes.buildModel()
b = moose.element( '/model/chem/dend/b' )
b.concInit *= 5









    



Rdesigneur: Elec model has 1 compartments and 0 spines on 1 compartments.
Chem part of model has  1 dendrite voxels X 3 pools,



In [4]:

    
moose.reinit()
moose.start( 200 )

rdes.display()