XPP model

This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6.

# XPP scripts for supplementary Figure S5A of
# Maex R, Grinevich VP, Grinevich V, Budygin E, Bencherif M, Gutkin B (2014)
# Understanding the role a7 nicotinic receptors play in dopamine efflux
# in nucleus accumbens. ACS Chemical Neuroscience 5, 1032-1040.
#
# Here I want to implement the facilitation of the receptor
# at very low (partial) agonist concentration, as described 
# by Smulders et al. (2005) and Prickaerts et al. (2012).

# The graph plots the fraction of open a7 channels across time.

# Naming conventions 
#
# V_      = membrane voltage of (neuron or neuron population)
# R_      = release of (transmitter)
# C_      = concentration of (transmitter)
# I_      = membrane current of (channel or receptor)
# P_      = presynatic membrane current of (receptor)
# tau_    = time-constant of
# tmin_ = minimum value of time-constant of
# tmax_ = maximum value of time-constant of
# inf_    = steady-state value of
# act_    = level of activation of (receptor)
# des_    = level of desensitisation of (receptor)
# w_      = weight of
# EC50_   = halfmax concentration of steady-state response
# Ktau_   = halfmax concentration of time-constant
# stim_   = stimulation by

# glu  = glutamate or glu-ergic neuron (population)
# dop  = dopamine or dopamine-ergic neuron (population)
# gab  = gaba or gaba-ergic neuron (population)
# ach  = acetylcholine
# nic  = nicotine
# a7   = alpha7-type nicACh receptor
# a4b2 = alpha4-beta2-type nicACh receptor
# bas  = basal
#


# Units 
# time seconds
# concentration micromolar


# Transfer functions (input-output) 
# sigmoid
f(x) = 1/(1 + exp(-x)) 
# halfwave rectification
hwr(x) = heav(x) * x    
# clipped at 0 and 1
clip(x) = (1 + (x-1) * (heav (1-x))) * heav(x) 
# Hill equation
# Hill(x,K,n) = 1 / (1 + (K/x)^n)
Hill(x,K,n) = x^n / (x^n + K^n) 
# competitive Hill 
compHill(x,y,z,wx,wy,wz,Kx,Ky,Kz,n) = \
  (wx * x^n/(x^n + Kx^n * (1 + (y/Ky)^n + (z/Kz)^n))) + \
  (wy * y^n/(y^n + Ky^n * (1 + (x/Kx)^n + (z/Kz)^n))) + \
  (wz * z^n/(z^n + Kz^n * (1 + (x/Kx)^n + (y/Ky)^n)))
# steady-state Hill function with cross-terms representing the effect of
# the simultaneous binding of two different compounds, for instance ACh and EVP 
Hillcross(x,y,z,wx,wy,wz,Kx,Ky,Kz,n,fxy,fyz,fxz) = \
  ((wx * (x/Kx)^n) + \
   (wy * (y/Ky)^n) + \
   (wz * (z/Kz)^n) + \
   (fxy * (x/Kx) * (y/Ky)) + \
   (fxz * (x/Kx) * (z/Kz)) + \
   (fyz * (y/Kz) * (z/Kz))) / (1 + (z/Kz)^n + (x/Kx)^n + (y/Ky)^n + \
                                   2 * (z/Kz) + 2 * (x/Kx) + 2 * (y/Ky) + \
                                   2 * (x/Kx) * (y/Ky) + \
                                   2 * (x/Kx) * (z/Kz) + \
                                   2 * (y/Ky) * (z/Kz) )
# rectangular pulse
rect_pulse(t,from,to,amp) = amp * (heav (t-from) - heav (t-to))



# THE MODEL 

# We now have two neuron populations sharing inputs from
# glu neurons (mostly a7-driven) and through a4b2 receptors.
# The parameter coefficients r and s specify the balance 
# between these inputs.

   par s=0.8,r=0.8


# the dynamics of the dopaminergic neuron population V_dop

   V_dop' = (-V_dop - I_gab + hwr (I_basDop + s*I_glu + r*I_b2)) / tau_Vdop

   par tau_Vdop=0.02,I_basDop=0.1
   aux aux_hwr=hwr (I_basDop + s*I_glu + r*I_b2)


# the dynamics of the gaba-ergic neuron population V_gab

   V_gab' = (-V_gab + hwr ((1-s)*I_glu + (1-r)*I_b2)) / tau_Vgab

   par tau_Vgab=0.02


# the (stationary) glutamatergic input I_glu to dopamine (and GABA)neurons

   I_glu = w_glu * clip (V_glu + P_glu)

   par w_glu=1,V_glu=0.1


# the (varying) gaba-ergic input I_gab to dopamine neurons

   I_gab = w_gab * clip (V_gab + I_basGab)
   aux auxI_gab=I_gab
 
   par w_gab=1.5
   par I_basGab=-0.05
   aux auxI_gab=I_gab


# the presynaptic facilitation of glu input to dopamine (and GABA) neurons

   P_glu = P_a7


# the dynamics of presynaptic a7 receptors (same single variable for 
# glu inputs to dopamine cells in VTA and medium spiny neurons in striatum)

   P_a7 = act_a7 * (1 - des_a7)

   aux auxP_a7=P_a7
   aux auxact_a7=act_a7
   aux auxsen_a7=1-des_a7

#   /* activation of alpha7 */

   act_a7' = (- act_a7 + inf_actA7) / tau_actA7

   inf_actA7 = Hillcross(C_ach + C_chol,        C_nic,        C_agA7, \
                                      1.0,   w_actNicA7,     w_actAgA7, \
                             EC50_A7ach,   EC50_A7nic,     EC50_A7ag, Hill_actA7, \
                                      0,            0,    crossactA7) 
   par crossactA7=0.3

   aux   ainf_actA7=inf_actA7

   par EC50_A7=80
   par EC50_A7ach=68
   par EC50_A7nic=13
   par Hill_actA7=1.73
   par EC50_A7ag=0.03 
   par w_actNicA7=0.8
   par w_actAgA7=0
   par tau_actA7=0.005


#   /* desensitisation of alpha7 */

   des_a7' = (- des_a7 + inf_desA7) / tau_desA7

   inf_desA7 = compHill(    0,      C_nic,      C_agA7, \
                            1,          1,           0, \
                            1, IC50_A7nic,   IC50_A7ag, Hill_desA7)

   par  crosssenA7=0

   aux ainf_desA7=inf_desA7

   par IC50_A7ach=1.3
   par IC50_A7nic=1.3
   par IC50_A7ag=0.002
   par Hill_desA7=2

   tau_desA7 = tmin_desA7 + \ 
               tmax_desA7 * (1 - inf_desA7)
   aux aux_tdA7=tau_desA7

   par tmin_desA7=0.05
   par tmax_desA7=120


# the dynamics of somatic a4b2 receptors on the soma/dendrite of dopamine neurons

   I_b2 = act_b2 * (1 - des_b2)

   aux auxI_b2=I_b2
   aux auxact_b2=act_b2
   aux auxsen_b2=1-des_b2

#   /* activation of alpha4beta2 */

   act_b2' = (- act_b2 + inf_actB2) / tau_actB2

   inf_actB2 = compHill(C_ach,      C_nic,        C_agB2, \
                            1, w_actNicB2,   w_actAgB2, \
                   EC50_B2ach, EC50_B2nic,   EC50_B2ag, Hill_actB2)


   par EC50_actB2=30
   par EC50_B2ach=30
   par EC50_B2ag=30
   par EC50_B2nic=0.23
   par Hill_actB2=1.05
   par w_actNicB2=2
   par w_actAgB2=1
   par tau_actB2=0.005

#   /* desensitisation */

   des_b2' = (- des_b2 + inf_desB2) / tau_desB2

   inf_desB2 = Hill (C_nic + C_agB2, EC50_desB2, Hill_desB2)

   tau_desB2 = tmin_desB2 + \ 
               tmax_desB2 * Hill (Ktau_desB2, C_nic + C_agB2, Htau_desB2)
   aux aux_tdB2=tau_desB2

   par EC50_desB2=0.061
   par Hill_desB2=0.5
   par tmin_desB2=0.5
   par tmax_desB2=600
   par Ktau_desB2=0.11
   par Htau_desB2=3


# the dynamics of dopamine release AND RE-UPTAKE

   R_dop' = (C_basDop * (1 + (V_dop-ss_Vdop) / ss_Vdop) - R_dop) / tau_Rdop - R_dop * Vmax_Rdop / (R_dop + EC50_Rdop)

   par tau_Rdop=0.2,P_a7bas=0.01
# the assumed basal DA production, denotated C_basDop, in the absence of re-uptake 
   par C_basDop=0.1
# the following values are from Chen and Budygin 2007
   par Vmax_Rdop=1.3
   par EC50_Rdop=0.2


# the physiological (ach) and pharmacological (nic) stimuli

#   /* acetylcoline HERE I USE CHOLINE BECAUSE THIS SELECTIVELY ACTIVATES THE A7s */
 
    C_ach' = (-C_ach + stim_ach) / 1

    stim_chol = bas_chol + rect_pulse (t, 1, 10, dose_chol)  + rect_pulse (t, 31, 35, dose_chol)  \
                    + rect_pulse (t, 61, 65, dose_chol)  \
         	    + rect_pulse (t, 91, 95, dose_chol)  \
		    + rect_pulse (t, 121, 125, dose_chol) \
		    + rect_pulse (t, 151, 155, dose_chol) \
		    + rect_pulse (t, 181, 185, dose_chol) \
		    + rect_pulse (t, 211, 215, dose_chol) \
		    + rect_pulse (t, 241, 245, dose_chol) \
		    + rect_pulse (t, 271, 275, dose_chol) \
		    + rect_pulse (t, 301, 305, dose_chol) \
                    + rect_pulse (t, 331, 335, dose_chol) \
		    + rect_pulse (t, 361, 365, dose_chol) \
                    + rect_pulse (t, 391, 395, dose_chol) \
		    + rect_pulse (t, 421, 425, dose_chol) \
                    + rect_pulse (t, 481, 485, dose_chol) \
                    + rect_pulse (t, 541, 545, dose_chol) \
                    + rect_pulse (t, 571, 575, dose_chol) \
                    + rect_pulse (t, 601, 605, dose_chol) \
                    + rect_pulse (t, 631, 635, dose_chol) 
 
    C_chol' = (-C_chol + stim_chol) / 1

    par bas_ach=0
    par bas_chol=0
    par dose_ach=0
    par stim_ach=0

# FOR THE UPPER PLOT PULSES OF 1 microM ARE APPLIED
    par dose_chol=1
# FOR THE LOWER PLOT PULSES OF 60 microM ARE APPLIED
#    par dose_chol=60


#   /* nicotine */
    stim_nic = bas_nic + rect_pulse (t, 66, 67, dose_Nic)

    par dose_Nic=0

# we now simulate an alpha function as a second-order ode
# the alpha-function is area normalized
# to calculate its peak, take amplitude and divide by 10*e

    C_nic_nic' = (-C_nic_nic + stim_nic)/tau1_nic
    C_nic' = (-C_nic + C_nic_nic)/tau2_nic

    par tau1_nic=10,tau2_nic=10,bas_nic=0 


#   /* alpha7 agonist */
    stim_a7 = rect_pulse (t, 170, 510, dose_Aga7)

    par dose_Aga7=0.003

    C_agA7A7' = (-C_agA7A7 + stim_a7) / tau1_agA7
    C_agA7' = (-C_agA7 + C_agA7A7) / tau2_agA7

    par tau1_agA7=10,tau2_agA7=10 


#   /* alpha4 beta2 agonist */
    stim_b2 = rect_pulse (t, 60, 61, 0)
    C_agB2' = (-C_agB2 + stim_b2) / tau_agB2

    par tau_agB2=10


# /* calculating the steady-states (in reverse order) */
    ss_Ib2 = compHill(   bas_ach,          0,           0, \
                            1, w_actNicB2,   w_actAgB2, \
                   EC50_B2ach, EC50_B2nic,   EC50_B2ag, Hill_actB2)
    ss_Pa7 = compHill(   bas_ach + bas_chol,          0,           0, \
                            1, w_actNicA7,   w_actAgA7, \
                   EC50_A7ach, EC50_A7nic,   EC50_A7ag, Hill_actA7)
    ss_Pglu = ss_Pa7
    ss_Iglu = w_glu * clip (V_glu + ss_Pglu)
    ss_Vgab = hwr ((1-s)*ss_Iglu + (1-r)*ss_Ib2)
    ss_Igab = w_gab * clip (ss_Vgab + I_basGab)
    ss_Vdop = - ss_Igab + hwr (I_basDop + s*ss_Iglu + r*ss_Ib2)
# aux auxss_Vdop=ss_Vdop


# the initial conditions
act_a7(0)=0
des_a7(0)=0
act_b2(0)=0
des_b2(0)=0
V_dop(0)=0
V_gab(0)=0
C_nic(0)=bas_nic
C_nic_nic(0)=bas_nic


# the numerical parameters
@ MAXSTOR=500000
@ BOUNDS=1e+8
@ DT=0.005
@ TOTAL=700
@ xlo=0
@ xhi=700
@ ylo=0
@ yhi=0.001
@ METHOD=rk
@ NJP=200
@ YP=auxP_a7

done

This file has been produced by sbmlutils.

Terms of use

Redistribution and use of any part of this model, with or without modification, are permitted provided that the following conditions are met:

Redistributions of this SBML file must retain the above copyright notice, this list of conditions and the following disclaimer.
Redistributions in a different form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.

This model is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

Model :

id
name
time
substance
extent
volume
area
length

L3V1

type	FunctionDefinitions [11]	name	math
FunctionDefinition	max max <functionDefinition id="max" name="minimum"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <bvar> <ci> y </ci> </bvar> <piecewise> <piece> <ci> x </ci> <apply> <gt/> <ci> x </ci> <ci> y </ci> </apply> </piece> <otherwise> <ci> y </ci> </otherwise> </piecewise> </lambda> </math> </functionDefinition>	minimum	$x y x x y y$
FunctionDefinition	min min <functionDefinition id="min" name="maximum"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <bvar> <ci> y </ci> </bvar> <piecewise> <piece> <ci> x </ci> <apply> <lt/> <ci> x </ci> <ci> y </ci> </apply> </piece> <otherwise> <ci> y </ci> </otherwise> </piecewise> </lambda> </math> </functionDefinition>	maximum	$x y x x y y$
FunctionDefinition	heav heav <functionDefinition id="heav" name="heavyside"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <piecewise> <piece> <cn type="integer"> 0 </cn> <apply> <lt/> <ci> x </ci> <cn type="integer"> 0 </cn> </apply> </piece> <piece> <cn> 0.5 </cn> <apply> <eq/> <ci> x </ci> <cn type="integer"> 0 </cn> </apply> </piece> <piece> <cn type="integer"> 1 </cn> <apply> <gt/> <ci> x </ci> <cn type="integer"> 0 </cn> </apply> </piece> <otherwise> <cn type="integer"> 0 </cn> </otherwise> </piecewise> </lambda> </math> </functionDefinition>	heavyside	$x 0 x 0 0.5 x 0 1 x 0 0$
FunctionDefinition	mod mod <functionDefinition id="mod" name="modulo"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <bvar> <ci> y </ci> </bvar> <piecewise> <piece> <apply> <minus/> <ci> x </ci> <apply> <times/> <ci> y </ci> <apply> <ceiling/> <apply> <divide/> <ci> x </ci> <ci> y </ci> </apply> </apply> </apply> </apply> <apply> <xor/> <apply> <lt/> <ci> x </ci> <cn type="integer"> 0 </cn> </apply> <apply> <lt/> <ci> y </ci> <cn type="integer"> 0 </cn> </apply> </apply> </piece> <otherwise> <apply> <minus/> <ci> x </ci> <apply> <times/> <ci> y </ci> <apply> <floor/> <apply> <divide/> <ci> x </ci> <ci> y </ci> </apply> </apply> </apply> </apply> </otherwise> </piecewise> </lambda> </math> </functionDefinition>	modulo	$x y x y x y x 0 y 0 x y x y$
FunctionDefinition	f f <functionDefinition id="f"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <apply> <divide/> <cn type="integer"> 1 </cn> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <exp/> <apply> <minus/> <ci> x </ci> </apply> </apply> </apply> </apply> </lambda> </math> </functionDefinition>		$x 1 1 x$
FunctionDefinition	hwr hwr <functionDefinition id="hwr"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <apply> <times/> <apply> <ci> heav </ci> <ci> x </ci> </apply> <ci> x </ci> </apply> </lambda> </math> </functionDefinition>		$x heav x x$
FunctionDefinition	clip clip <functionDefinition id="clip"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <apply> <times/> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <times/> <apply> <minus/> <ci> x </ci> <cn type="integer"> 1 </cn> </apply> <apply> <ci> heav </ci> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> x </ci> </apply> </apply> </apply> </apply> <apply> <ci> heav </ci> <ci> x </ci> </apply> </apply> </lambda> </math> </functionDefinition>		$x 1 x 1 heav 1 x heav x$
FunctionDefinition	hill hill <functionDefinition id="hill"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <bvar> <ci> k </ci> </bvar> <bvar> <ci> n </ci> </bvar> <apply> <divide/> <apply> <power/> <ci> x </ci> <ci> n </ci> </apply> <apply> <plus/> <apply> <power/> <ci> x </ci> <ci> n </ci> </apply> <apply> <power/> <ci> k </ci> <ci> n </ci> </apply> </apply> </apply> </lambda> </math> </functionDefinition>		$x k n x n x n k n$
FunctionDefinition	comphill comphill <functionDefinition id="comphill"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <bvar> <ci> y </ci> </bvar> <bvar> <ci> z </ci> </bvar> <bvar> <ci> wx </ci> </bvar> <bvar> <ci> wy </ci> </bvar> <bvar> <ci> wz </ci> </bvar> <bvar> <ci> kx </ci> </bvar> <bvar> <ci> ky </ci> </bvar> <bvar> <ci> kz </ci> </bvar> <bvar> <ci> n </ci> </bvar> <apply> <plus/> <apply> <divide/> <apply> <times/> <ci> wx </ci> <apply> <power/> <ci> x </ci> <ci> n </ci> </apply> </apply> <apply> <plus/> <apply> <power/> <ci> x </ci> <ci> n </ci> </apply> <apply> <times/> <apply> <power/> <ci> kx </ci> <ci> n </ci> </apply> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <power/> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> <ci> n </ci> </apply> <apply> <power/> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> <ci> n </ci> </apply> </apply> </apply> </apply> </apply> <apply> <divide/> <apply> <times/> <ci> wy </ci> <apply> <power/> <ci> y </ci> <ci> n </ci> </apply> </apply> <apply> <plus/> <apply> <power/> <ci> y </ci> <ci> n </ci> </apply> <apply> <times/> <apply> <power/> <ci> ky </ci> <ci> n </ci> </apply> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <power/> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <ci> n </ci> </apply> <apply> <power/> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> <ci> n </ci> </apply> </apply> </apply> </apply> </apply> <apply> <divide/> <apply> <times/> <ci> wz </ci> <apply> <power/> <ci> z </ci> <ci> n </ci> </apply> </apply> <apply> <plus/> <apply> <power/> <ci> z </ci> <ci> n </ci> </apply> <apply> <times/> <apply> <power/> <ci> kz </ci> <ci> n </ci> </apply> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <power/> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <ci> n </ci> </apply> <apply> <power/> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> <ci> n </ci> </apply> </apply> </apply> </apply> </apply> </apply> </lambda> </math> </functionDefinition>		$x y z wx wy wz kx ky kz n wx x n x n kx n 1 y ky n z kz n wy y n y n ky n 1 x kx n z kz n wz z n z n kz n 1 x kx n y ky n$
FunctionDefinition	hillcross hillcross <functionDefinition id="hillcross"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> x </ci> </bvar> <bvar> <ci> y </ci> </bvar> <bvar> <ci> z </ci> </bvar> <bvar> <ci> wx </ci> </bvar> <bvar> <ci> wy </ci> </bvar> <bvar> <ci> wz </ci> </bvar> <bvar> <ci> kx </ci> </bvar> <bvar> <ci> ky </ci> </bvar> <bvar> <ci> kz </ci> </bvar> <bvar> <ci> n </ci> </bvar> <bvar> <ci> fxy </ci> </bvar> <bvar> <ci> fyz </ci> </bvar> <bvar> <ci> fxz </ci> </bvar> <apply> <divide/> <apply> <plus/> <apply> <times/> <ci> wx </ci> <apply> <power/> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <ci> n </ci> </apply> </apply> <apply> <times/> <ci> wy </ci> <apply> <power/> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> <ci> n </ci> </apply> </apply> <apply> <times/> <ci> wz </ci> <apply> <power/> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> <ci> n </ci> </apply> </apply> <apply> <times/> <ci> fxy </ci> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> </apply> <apply> <times/> <ci> fxz </ci> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> </apply> <apply> <times/> <ci> fyz </ci> <apply> <divide/> <ci> y </ci> <ci> kz </ci> </apply> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> </apply> </apply> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <power/> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> <ci> n </ci> </apply> <apply> <power/> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <ci> n </ci> </apply> <apply> <power/> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> <ci> n </ci> </apply> <apply> <times/> <cn type="integer"> 2 </cn> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> </apply> <apply> <times/> <cn type="integer"> 2 </cn> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> </apply> <apply> <times/> <cn type="integer"> 2 </cn> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> </apply> <apply> <times/> <cn type="integer"> 2 </cn> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> </apply> <apply> <times/> <cn type="integer"> 2 </cn> <apply> <divide/> <ci> x </ci> <ci> kx </ci> </apply> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> </apply> <apply> <times/> <cn type="integer"> 2 </cn> <apply> <divide/> <ci> y </ci> <ci> ky </ci> </apply> <apply> <divide/> <ci> z </ci> <ci> kz </ci> </apply> </apply> </apply> </apply> </lambda> </math> </functionDefinition>		$x y z wx wy wz kx ky kz n fxy fyz fxz wx x kx n wy y ky n wz z kz n fxy x kx y ky fxz x kx z kz fyz y kz z kz 1 z kz n x kx n y ky n 2 z kz 2 x kx 2 y ky 2 x kx y ky 2 x kx z kz 2 y ky z kz$
FunctionDefinition	rect_pulse rect_pulse <functionDefinition id="rect_pulse"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <lambda> <bvar> <ci> t </ci> </bvar> <bvar> <ci> from </ci> </bvar> <bvar> <ci> to </ci> </bvar> <bvar> <ci> amp </ci> </bvar> <apply> <times/> <ci> amp </ci> <apply> <minus/> <apply> <ci> heav </ci> <apply> <minus/> <ci> t </ci> <ci> from </ci> </apply> </apply> <apply> <ci> heav </ci> <apply> <minus/> <ci> t </ci> <ci> to </ci> </apply> </apply> </apply> </apply> </lambda> </math> </functionDefinition>		$t from to amp amp heav t from heav t to$

type	Parameters [106]	name	constant	value	unit	derived unit
Parameter	s s <parameter id="s" name="s = 0.8" value="0.8" constant="false"/>	s = 0.8	F	0.8		None
Parameter	r r <parameter id="r" name="r = 0.8" value="0.8" constant="false"/>	r = 0.8	F	0.8		None
Parameter	tau_vdop tau_vdop <parameter id="tau_vdop" name="tau_vdop = 0.02" value="0.02" constant="false"/>	tau_vdop = 0.02	F	0.02		None
Parameter	i_basdop i_basdop <parameter id="i_basdop" name="i_basdop = 0.1" value="0.1" constant="false"/>	i_basdop = 0.1	F	0.1		None
Parameter	tau_vgab tau_vgab <parameter id="tau_vgab" name="tau_vgab = 0.02" value="0.02" constant="false"/>	tau_vgab = 0.02	F	0.02		None
Parameter	w_glu w_glu <parameter id="w_glu" name="w_glu = 1" value="1" constant="false"/>	w_glu = 1	F	1.0		None
Parameter	v_glu v_glu <parameter id="v_glu" name="v_glu = 0.1" value="0.1" constant="false"/>	v_glu = 0.1	F	0.1		None
Parameter	w_gab w_gab <parameter id="w_gab" name="w_gab = 1.5" value="1.5" constant="false"/>	w_gab = 1.5	F	1.5		None
Parameter	i_basgab i_basgab <parameter id="i_basgab" name="i_basgab = -0.05" value="-0.05" constant="false"/>	i_basgab = -0.05	F	-0.05		None
Parameter	crossacta7 crossacta7 <parameter id="crossacta7" name="crossacta7 = 0.3" value="0.3" constant="false"/>	crossacta7 = 0.3	F	0.3		None
Parameter	ec50_a7 ec50_a7 <parameter id="ec50_a7" name="ec50_a7 = 80" value="80" constant="false"/>	ec50_a7 = 80	F	80.0		None
Parameter	ec50_a7ach ec50_a7ach <parameter id="ec50_a7ach" name="ec50_a7ach = 68" value="68" constant="false"/>	ec50_a7ach = 68	F	68.0		None
Parameter	ec50_a7nic ec50_a7nic <parameter id="ec50_a7nic" name="ec50_a7nic = 13" value="13" constant="false"/>	ec50_a7nic = 13	F	13.0		None
Parameter	hill_acta7 hill_acta7 <parameter id="hill_acta7" name="hill_acta7 = 1.73" value="1.73" constant="false"/>	hill_acta7 = 1.73	F	1.73		None
Parameter	ec50_a7ag ec50_a7ag <parameter id="ec50_a7ag" name="ec50_a7ag = 0.03" value="0.03" constant="false"/>	ec50_a7ag = 0.03	F	0.03		None
Parameter	w_actnica7 w_actnica7 <parameter id="w_actnica7" name="w_actnica7 = 0.8" value="0.8" constant="false"/>	w_actnica7 = 0.8	F	0.8		None
Parameter	w_actaga7 w_actaga7 <parameter id="w_actaga7" name="w_actaga7 = 0" value="0" constant="false"/>	w_actaga7 = 0	F	0.0		None
Parameter	tau_acta7 tau_acta7 <parameter id="tau_acta7" name="tau_acta7 = 0.005" value="0.005" constant="false"/>	tau_acta7 = 0.005	F	0.005		None
Parameter	crosssena7 crosssena7 <parameter id="crosssena7" name="crosssena7 = 0" value="0" constant="false"/>	crosssena7 = 0	F	0.0		None
Parameter	ic50_a7ach ic50_a7ach <parameter id="ic50_a7ach" name="ic50_a7ach = 1.3" value="1.3" constant="false"/>	ic50_a7ach = 1.3	F	1.3		None
Parameter	ic50_a7nic ic50_a7nic <parameter id="ic50_a7nic" name="ic50_a7nic = 1.3" value="1.3" constant="false"/>	ic50_a7nic = 1.3	F	1.3		None
Parameter	ic50_a7ag ic50_a7ag <parameter id="ic50_a7ag" name="ic50_a7ag = 0.002" value="0.002" constant="false"/>	ic50_a7ag = 0.002	F	0.002		None
Parameter	hill_desa7 hill_desa7 <parameter id="hill_desa7" name="hill_desa7 = 2" value="2" constant="false"/>	hill_desa7 = 2	F	2.0		None
Parameter	tmin_desa7 tmin_desa7 <parameter id="tmin_desa7" name="tmin_desa7 = 0.05" value="0.05" constant="false"/>	tmin_desa7 = 0.05	F	0.05		None
Parameter	tmax_desa7 tmax_desa7 <parameter id="tmax_desa7" name="tmax_desa7 = 120" value="120" constant="false"/>	tmax_desa7 = 120	F	120.0		None
Parameter	ec50_actb2 ec50_actb2 <parameter id="ec50_actb2" name="ec50_actb2 = 30" value="30" constant="false"/>	ec50_actb2 = 30	F	30.0		None
Parameter	ec50_b2ach ec50_b2ach <parameter id="ec50_b2ach" name="ec50_b2ach = 30" value="30" constant="false"/>	ec50_b2ach = 30	F	30.0		None
Parameter	ec50_b2ag ec50_b2ag <parameter id="ec50_b2ag" name="ec50_b2ag = 30" value="30" constant="false"/>	ec50_b2ag = 30	F	30.0		None
Parameter	ec50_b2nic ec50_b2nic <parameter id="ec50_b2nic" name="ec50_b2nic = 0.23" value="0.23" constant="false"/>	ec50_b2nic = 0.23	F	0.23		None
Parameter	hill_actb2 hill_actb2 <parameter id="hill_actb2" name="hill_actb2 = 1.05" value="1.05" constant="false"/>	hill_actb2 = 1.05	F	1.05		None
Parameter	w_actnicb2 w_actnicb2 <parameter id="w_actnicb2" name="w_actnicb2 = 2" value="2" constant="false"/>	w_actnicb2 = 2	F	2.0		None
Parameter	w_actagb2 w_actagb2 <parameter id="w_actagb2" name="w_actagb2 = 1" value="1" constant="false"/>	w_actagb2 = 1	F	1.0		None
Parameter	tau_actb2 tau_actb2 <parameter id="tau_actb2" name="tau_actb2 = 0.005" value="0.005" constant="false"/>	tau_actb2 = 0.005	F	0.005		None
Parameter	ec50_desb2 ec50_desb2 <parameter id="ec50_desb2" name="ec50_desb2 = 0.061" value="0.061" constant="false"/>	ec50_desb2 = 0.061	F	0.061		None
Parameter	hill_desb2 hill_desb2 <parameter id="hill_desb2" name="hill_desb2 = 0.5" value="0.5" constant="false"/>	hill_desb2 = 0.5	F	0.5		None
Parameter	tmin_desb2 tmin_desb2 <parameter id="tmin_desb2" name="tmin_desb2 = 0.5" value="0.5" constant="false"/>	tmin_desb2 = 0.5	F	0.5		None
Parameter	tmax_desb2 tmax_desb2 <parameter id="tmax_desb2" name="tmax_desb2 = 600" value="600" constant="false"/>	tmax_desb2 = 600	F	600.0		None
Parameter	ktau_desb2 ktau_desb2 <parameter id="ktau_desb2" name="ktau_desb2 = 0.11" value="0.11" constant="false"/>	ktau_desb2 = 0.11	F	0.11		None
Parameter	htau_desb2 htau_desb2 <parameter id="htau_desb2" name="htau_desb2 = 3" value="3" constant="false"/>	htau_desb2 = 3	F	3.0		None
Parameter	tau_rdop tau_rdop <parameter id="tau_rdop" name="tau_rdop = 0.2" value="0.2" constant="false"/>	tau_rdop = 0.2	F	0.2		None
Parameter	p_a7bas p_a7bas <parameter id="p_a7bas" name="p_a7bas = 0.01" value="0.01" constant="false"/>	p_a7bas = 0.01	F	0.01		None
Parameter	c_basdop c_basdop <parameter id="c_basdop" name="c_basdop = 0.1" value="0.1" constant="false"/>	c_basdop = 0.1	F	0.1		None
Parameter	vmax_rdop vmax_rdop <parameter id="vmax_rdop" name="vmax_rdop = 1.3" value="1.3" constant="false"/>	vmax_rdop = 1.3	F	1.3		None
Parameter	ec50_rdop ec50_rdop <parameter id="ec50_rdop" name="ec50_rdop = 0.2" value="0.2" constant="false"/>	ec50_rdop = 0.2	F	0.2		None
Parameter	bas_ach bas_ach <parameter id="bas_ach" name="bas_ach = 0" value="0" constant="false"/>	bas_ach = 0	F	0.0		None
Parameter	bas_chol bas_chol <parameter id="bas_chol" name="bas_chol = 0" value="0" constant="false"/>	bas_chol = 0	F	0.0		None
Parameter	dose_ach dose_ach <parameter id="dose_ach" name="dose_ach = 0" value="0" constant="false"/>	dose_ach = 0	F	0.0		None
Parameter	stim_ach stim_ach <parameter id="stim_ach" name="stim_ach = 0" value="0" constant="false"/>	stim_ach = 0	F	0.0		None
Parameter	dose_chol dose_chol <parameter id="dose_chol" name="dose_chol = 1" value="1" constant="false"/>	dose_chol = 1	F	1.0		None
Parameter	dose_nic dose_nic <parameter id="dose_nic" name="dose_nic = 0" value="0" constant="false"/>	dose_nic = 0	F	0.0		None
Parameter	tau1_nic tau1_nic <parameter id="tau1_nic" name="tau1_nic = 10" value="10" constant="false"/>	tau1_nic = 10	F	10.0		None
Parameter	tau2_nic tau2_nic <parameter id="tau2_nic" name="tau2_nic = 10" value="10" constant="false"/>	tau2_nic = 10	F	10.0		None
Parameter	bas_nic bas_nic <parameter id="bas_nic" name="bas_nic = 0" value="0" constant="false"/>	bas_nic = 0	F	0.0		None
Parameter	dose_aga7 dose_aga7 <parameter id="dose_aga7" name="dose_aga7 = 0.003" value="0.003" constant="false"/>	dose_aga7 = 0.003	F	0.003		None
Parameter	tau1_aga7 tau1_aga7 <parameter id="tau1_aga7" name="tau1_aga7 = 10" value="10" constant="false"/>	tau1_aga7 = 10	F	10.0		None
Parameter	tau2_aga7 tau2_aga7 <parameter id="tau2_aga7" name="tau2_aga7 = 10" value="10" constant="false"/>	tau2_aga7 = 10	F	10.0		None
Parameter	tau_agb2 tau_agb2 <parameter id="tau_agb2" name="tau_agb2 = 10" value="10" constant="false"/>	tau_agb2 = 10	F	10.0		None
Parameter	act_a7 act_a7 <parameter id="act_a7" name="act_a7 = 0" value="0" constant="false"/>	act_a7 = 0	F	0.0		None
Parameter	des_a7 des_a7 <parameter id="des_a7" name="des_a7 = 0" value="0" constant="false"/>	des_a7 = 0	F	0.0		None
Parameter	act_b2 act_b2 <parameter id="act_b2" name="act_b2 = 0" value="0" constant="false"/>	act_b2 = 0	F	0.0		None
Parameter	des_b2 des_b2 <parameter id="des_b2" name="des_b2 = 0" value="0" constant="false"/>	des_b2 = 0	F	0.0		None
Parameter	v_dop v_dop <parameter id="v_dop" name="v_dop = 0" value="0" constant="false"/>	v_dop = 0	F	0.0		None
Parameter	v_gab v_gab <parameter id="v_gab" name="v_gab = 0" value="0" constant="false"/>	v_gab = 0	F	0.0		None
Parameter	c_nic c_nic <parameter id="c_nic" name="c_nic" value="0" constant="false"/>	c_nic	F	0.0		None
Parameter	c_nic_nic c_nic_nic <parameter id="c_nic_nic" name="c_nic_nic" value="0" constant="false"/>	c_nic_nic	F	0.0		None
Parameter	r_dop r_dop <parameter id="r_dop" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	c_ach c_ach <parameter id="c_ach" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	c_chol c_chol <parameter id="c_chol" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	c_aga7a7 c_aga7a7 <parameter id="c_aga7a7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	c_aga7 c_aga7 <parameter id="c_aga7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	c_agb2 c_agb2 <parameter id="c_agb2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	aux_hwr aux_hwr <parameter id="aux_hwr" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	i_glu i_glu <parameter id="i_glu" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	i_gab i_gab <parameter id="i_gab" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxi_gab auxi_gab <parameter id="auxi_gab" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	p_glu p_glu <parameter id="p_glu" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	p_a7 p_a7 <parameter id="p_a7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxp_a7 auxp_a7 <parameter id="auxp_a7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxact_a7 auxact_a7 <parameter id="auxact_a7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxsen_a7 auxsen_a7 <parameter id="auxsen_a7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	inf_acta7 inf_acta7 <parameter id="inf_acta7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ainf_acta7 ainf_acta7 <parameter id="ainf_acta7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	inf_desa7 inf_desa7 <parameter id="inf_desa7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ainf_desa7 ainf_desa7 <parameter id="ainf_desa7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	tau_desa7 tau_desa7 <parameter id="tau_desa7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	aux_tda7 aux_tda7 <parameter id="aux_tda7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	i_b2 i_b2 <parameter id="i_b2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxi_b2 auxi_b2 <parameter id="auxi_b2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxact_b2 auxact_b2 <parameter id="auxact_b2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	auxsen_b2 auxsen_b2 <parameter id="auxsen_b2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	inf_actb2 inf_actb2 <parameter id="inf_actb2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	inf_desb2 inf_desb2 <parameter id="inf_desb2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	tau_desb2 tau_desb2 <parameter id="tau_desb2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	aux_tdb2 aux_tdb2 <parameter id="aux_tdb2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	stim_chol stim_chol <parameter id="stim_chol" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	stim_nic stim_nic <parameter id="stim_nic" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	stim_a7 stim_a7 <parameter id="stim_a7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	stim_b2 stim_b2 <parameter id="stim_b2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_ib2 ss_ib2 <parameter id="ss_ib2" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_pa7 ss_pa7 <parameter id="ss_pa7" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_pglu ss_pglu <parameter id="ss_pglu" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_iglu ss_iglu <parameter id="ss_iglu" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_vgab ss_vgab <parameter id="ss_vgab" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_igab ss_igab <parameter id="ss_igab" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	ss_vdop ss_vdop <parameter id="ss_vdop" value="0" units="dimensionless" constant="false"/>		F	0.0	dimensionless	None
Parameter	t t <parameter id="t" name="model time" value="0" units="dimensionless" constant="false"/>	model time	F	0.0	dimensionless	None

type	InitialAssignments [2]	name		assignment	derived units	sbo	cvterm
InitialAssignment	c_nic		=	$bas_nic$	None
InitialAssignment	c_nic_nic		=	$bas_nic$	None

type	Rules [49]		assignment	derived units
Rule	d v_dop/dt v_dop <rateRule variable="v_dop"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <apply> <minus/> <ci> v_dop </ci> </apply> <ci> i_gab </ci> </apply> <apply> <ci> hwr </ci> <apply> <plus/> <ci> i_basdop </ci> <apply> <times/> <ci> s </ci> <ci> i_glu </ci> </apply> <apply> <times/> <ci> r </ci> <ci> i_b2 </ci> </apply> </apply> </apply> </apply> <ci> tau_vdop </ci> </apply> </math> </rateRule>	=	$v_dop i_gab hwr i_basdop s i_glu r i_b2 tau_vdop$	None
Rule	d v_gab/dt v_gab <rateRule variable="v_gab"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> v_gab </ci> </apply> <apply> <ci> hwr </ci> <apply> <plus/> <apply> <times/> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> s </ci> </apply> <ci> i_glu </ci> </apply> <apply> <times/> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> r </ci> </apply> <ci> i_b2 </ci> </apply> </apply> </apply> </apply> <ci> tau_vgab </ci> </apply> </math> </rateRule>	=	$v_gab hwr 1 s i_glu 1 r i_b2 tau_vgab$	None
Rule	d act_a7/dt act_a7 <rateRule variable="act_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> act_a7 </ci> </apply> <ci> inf_acta7 </ci> </apply> <ci> tau_acta7 </ci> </apply> </math> </rateRule>	=	$act_a7 inf_acta7 tau_acta7$	None
Rule	d des_a7/dt des_a7 <rateRule variable="des_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> des_a7 </ci> </apply> <ci> inf_desa7 </ci> </apply> <ci> tau_desa7 </ci> </apply> </math> </rateRule>	=	$des_a7 inf_desa7 tau_desa7$	None
Rule	d act_b2/dt act_b2 <rateRule variable="act_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> act_b2 </ci> </apply> <ci> inf_actb2 </ci> </apply> <ci> tau_actb2 </ci> </apply> </math> </rateRule>	=	$act_b2 inf_actb2 tau_actb2$	None
Rule	d des_b2/dt des_b2 <rateRule variable="des_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> des_b2 </ci> </apply> <ci> inf_desb2 </ci> </apply> <ci> tau_desb2 </ci> </apply> </math> </rateRule>	=	$des_b2 inf_desb2 tau_desb2$	None
Rule	d r_dop/dt r_dop <rateRule variable="r_dop"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <minus/> <apply> <divide/> <apply> <minus/> <apply> <times/> <ci> c_basdop </ci> <apply> <plus/> <cn type="integer"> 1 </cn> <apply> <divide/> <apply> <minus/> <ci> v_dop </ci> <ci> ss_vdop </ci> </apply> <ci> ss_vdop </ci> </apply> </apply> </apply> <ci> r_dop </ci> </apply> <ci> tau_rdop </ci> </apply> <apply> <divide/> <apply> <times/> <ci> r_dop </ci> <ci> vmax_rdop </ci> </apply> <apply> <plus/> <ci> r_dop </ci> <ci> ec50_rdop </ci> </apply> </apply> </apply> </math> </rateRule>	=	$c_basdop 1 v_dop ss_vdop ss_vdop r_dop tau_rdop r_dop vmax_rdop r_dop ec50_rdop$	None
Rule	d c_ach/dt c_ach <rateRule variable="c_ach"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_ach </ci> </apply> <ci> stim_ach </ci> </apply> <cn type="integer"> 1 </cn> </apply> </math> </rateRule>	=	$c_ach stim_ach 1$	None
Rule	d c_chol/dt c_chol <rateRule variable="c_chol"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_chol </ci> </apply> <ci> stim_chol </ci> </apply> <cn type="integer"> 1 </cn> </apply> </math> </rateRule>	=	$c_chol stim_chol 1$	None
Rule	d c_nic_nic/dt c_nic_nic <rateRule variable="c_nic_nic"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_nic_nic </ci> </apply> <ci> stim_nic </ci> </apply> <ci> tau1_nic </ci> </apply> </math> </rateRule>	=	$c_nic_nic stim_nic tau1_nic$	None
Rule	d c_nic/dt c_nic <rateRule variable="c_nic"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_nic </ci> </apply> <ci> c_nic_nic </ci> </apply> <ci> tau2_nic </ci> </apply> </math> </rateRule>	=	$c_nic c_nic_nic tau2_nic$	None
Rule	d c_aga7a7/dt c_aga7a7 <rateRule variable="c_aga7a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_aga7a7 </ci> </apply> <ci> stim_a7 </ci> </apply> <ci> tau1_aga7 </ci> </apply> </math> </rateRule>	=	$c_aga7a7 stim_a7 tau1_aga7$	None
Rule	d c_aga7/dt c_aga7 <rateRule variable="c_aga7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_aga7 </ci> </apply> <ci> c_aga7a7 </ci> </apply> <ci> tau2_aga7 </ci> </apply> </math> </rateRule>	=	$c_aga7 c_aga7a7 tau2_aga7$	None
Rule	d c_agb2/dt c_agb2 <rateRule variable="c_agb2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <divide/> <apply> <plus/> <apply> <minus/> <ci> c_agb2 </ci> </apply> <ci> stim_b2 </ci> </apply> <ci> tau_agb2 </ci> </apply> </math> </rateRule>	=	$c_agb2 stim_b2 tau_agb2$	None
Rule	aux_hwr aux_hwr <assignmentRule variable="aux_hwr"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> hwr </ci> <apply> <plus/> <ci> i_basdop </ci> <apply> <times/> <ci> s </ci> <ci> i_glu </ci> </apply> <apply> <times/> <ci> r </ci> <ci> i_b2 </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$hwr i_basdop s i_glu r i_b2$	None
Rule	i_glu i_glu <assignmentRule variable="i_glu"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <times/> <ci> w_glu </ci> <apply> <ci> clip </ci> <apply> <plus/> <ci> v_glu </ci> <ci> p_glu </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$w_glu clip v_glu p_glu$	None
Rule	i_gab i_gab <assignmentRule variable="i_gab"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <times/> <ci> w_gab </ci> <apply> <ci> clip </ci> <apply> <plus/> <ci> v_gab </ci> <ci> i_basgab </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$w_gab clip v_gab i_basgab$	None
Rule	auxi_gab auxi_gab <assignmentRule variable="auxi_gab"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> i_gab </ci> </math> </assignmentRule>	=	$i_gab$	None
Rule	p_glu p_glu <assignmentRule variable="p_glu"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> p_a7 </ci> </math> </assignmentRule>	=	$p_a7$	None
Rule	p_a7 p_a7 <assignmentRule variable="p_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <times/> <ci> act_a7 </ci> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> des_a7 </ci> </apply> </apply> </math> </assignmentRule>	=	$act_a7 1 des_a7$	None
Rule	auxp_a7 auxp_a7 <assignmentRule variable="auxp_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> p_a7 </ci> </math> </assignmentRule>	=	$p_a7$	None
Rule	auxact_a7 auxact_a7 <assignmentRule variable="auxact_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> act_a7 </ci> </math> </assignmentRule>	=	$act_a7$	None
Rule	auxsen_a7 auxsen_a7 <assignmentRule variable="auxsen_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> des_a7 </ci> </apply> </math> </assignmentRule>	=	$1 des_a7$	None
Rule	inf_acta7 inf_acta7 <assignmentRule variable="inf_acta7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> hillcross </ci> <apply> <plus/> <ci> c_ach </ci> <ci> c_chol </ci> </apply> <ci> c_nic </ci> <ci> c_aga7 </ci> <cn> 1 </cn> <ci> w_actnica7 </ci> <ci> w_actaga7 </ci> <ci> ec50_a7ach </ci> <ci> ec50_a7nic </ci> <ci> ec50_a7ag </ci> <ci> hill_acta7 </ci> <cn type="integer"> 0 </cn> <cn type="integer"> 0 </cn> <ci> crossacta7 </ci> </apply> </math> </assignmentRule>	=	$hillcross c_ach c_chol c_nic c_aga7 1 w_actnica7 w_actaga7 ec50_a7ach ec50_a7nic ec50_a7ag hill_acta7 00 crossacta7$	None
Rule	ainf_acta7 ainf_acta7 <assignmentRule variable="ainf_acta7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> inf_acta7 </ci> </math> </assignmentRule>	=	$inf_acta7$	None
Rule	inf_desa7 inf_desa7 <assignmentRule variable="inf_desa7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> comphill </ci> <cn type="integer"> 0 </cn> <ci> c_nic </ci> <ci> c_aga7 </ci> <cn type="integer"> 1 </cn> <cn type="integer"> 1 </cn> <cn type="integer"> 0 </cn> <cn type="integer"> 1 </cn> <ci> ic50_a7nic </ci> <ci> ic50_a7ag </ci> <ci> hill_desa7 </ci> </apply> </math> </assignmentRule>	=	$comphill 0 c_nic c_aga7 1101 ic50_a7nic ic50_a7ag hill_desa7$	None
Rule	ainf_desa7 ainf_desa7 <assignmentRule variable="ainf_desa7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> inf_desa7 </ci> </math> </assignmentRule>	=	$inf_desa7$	None
Rule	tau_desa7 tau_desa7 <assignmentRule variable="tau_desa7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <plus/> <ci> tmin_desa7 </ci> <apply> <times/> <ci> tmax_desa7 </ci> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> inf_desa7 </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$tmin_desa7 tmax_desa7 1 inf_desa7$	None
Rule	aux_tda7 aux_tda7 <assignmentRule variable="aux_tda7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> tau_desa7 </ci> </math> </assignmentRule>	=	$tau_desa7$	None
Rule	i_b2 i_b2 <assignmentRule variable="i_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <times/> <ci> act_b2 </ci> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> des_b2 </ci> </apply> </apply> </math> </assignmentRule>	=	$act_b2 1 des_b2$	None
Rule	auxi_b2 auxi_b2 <assignmentRule variable="auxi_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> i_b2 </ci> </math> </assignmentRule>	=	$i_b2$	None
Rule	auxact_b2 auxact_b2 <assignmentRule variable="auxact_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> act_b2 </ci> </math> </assignmentRule>	=	$act_b2$	None
Rule	auxsen_b2 auxsen_b2 <assignmentRule variable="auxsen_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> des_b2 </ci> </apply> </math> </assignmentRule>	=	$1 des_b2$	None
Rule	inf_actb2 inf_actb2 <assignmentRule variable="inf_actb2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> comphill </ci> <ci> c_ach </ci> <ci> c_nic </ci> <ci> c_agb2 </ci> <cn type="integer"> 1 </cn> <ci> w_actnicb2 </ci> <ci> w_actagb2 </ci> <ci> ec50_b2ach </ci> <ci> ec50_b2nic </ci> <ci> ec50_b2ag </ci> <ci> hill_actb2 </ci> </apply> </math> </assignmentRule>	=	$comphill c_ach c_nic c_agb2 1 w_actnicb2 w_actagb2 ec50_b2ach ec50_b2nic ec50_b2ag hill_actb2$	None
Rule	inf_desb2 inf_desb2 <assignmentRule variable="inf_desb2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> hill </ci> <apply> <plus/> <ci> c_nic </ci> <ci> c_agb2 </ci> </apply> <ci> ec50_desb2 </ci> <ci> hill_desb2 </ci> </apply> </math> </assignmentRule>	=	$hill c_nic c_agb2 ec50_desb2 hill_desb2$	None
Rule	tau_desb2 tau_desb2 <assignmentRule variable="tau_desb2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <plus/> <ci> tmin_desb2 </ci> <apply> <times/> <ci> tmax_desb2 </ci> <apply> <ci> hill </ci> <ci> ktau_desb2 </ci> <apply> <plus/> <ci> c_nic </ci> <ci> c_agb2 </ci> </apply> <ci> htau_desb2 </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$tmin_desb2 tmax_desb2 hill ktau_desb2 c_nic c_agb2 htau_desb2$	None
Rule	aux_tdb2 aux_tdb2 <assignmentRule variable="aux_tdb2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> tau_desb2 </ci> </math> </assignmentRule>	=	$tau_desb2$	None
Rule	stim_chol stim_chol <assignmentRule variable="stim_chol"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <plus/> <ci> bas_chol </ci> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 1 </cn> <cn type="integer"> 10 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 31 </cn> <cn type="integer"> 35 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 61 </cn> <cn type="integer"> 65 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 91 </cn> <cn type="integer"> 95 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 121 </cn> <cn type="integer"> 125 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 151 </cn> <cn type="integer"> 155 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 181 </cn> <cn type="integer"> 185 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 211 </cn> <cn type="integer"> 215 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 241 </cn> <cn type="integer"> 245 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 271 </cn> <cn type="integer"> 275 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 301 </cn> <cn type="integer"> 305 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 331 </cn> <cn type="integer"> 335 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 361 </cn> <cn type="integer"> 365 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 391 </cn> <cn type="integer"> 395 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 421 </cn> <cn type="integer"> 425 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 481 </cn> <cn type="integer"> 485 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 541 </cn> <cn type="integer"> 545 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 571 </cn> <cn type="integer"> 575 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 601 </cn> <cn type="integer"> 605 </cn> <ci> dose_chol </ci> </apply> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 631 </cn> <cn type="integer"> 635 </cn> <ci> dose_chol </ci> </apply> </apply> </math> </assignmentRule>	=	$bas_chol rect_pulse t 110 dose_chol rect_pulse t 3135 dose_chol rect_pulse t 6165 dose_chol rect_pulse t 9195 dose_chol rect_pulse t 121125 dose_chol rect_pulse t 151155 dose_chol rect_pulse t 181185 dose_chol rect_pulse t 211215 dose_chol rect_pulse t 241245 dose_chol rect_pulse t 271275 dose_chol rect_pulse t 301305 dose_chol rect_pulse t 331335 dose_chol rect_pulse t 361365 dose_chol rect_pulse t 391395 dose_chol rect_pulse t 421425 dose_chol rect_pulse t 481485 dose_chol rect_pulse t 541545 dose_chol rect_pulse t 571575 dose_chol rect_pulse t 601605 dose_chol rect_pulse t 631635 dose_chol$	None
Rule	stim_nic stim_nic <assignmentRule variable="stim_nic"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <plus/> <ci> bas_nic </ci> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 66 </cn> <cn type="integer"> 67 </cn> <ci> dose_nic </ci> </apply> </apply> </math> </assignmentRule>	=	$bas_nic rect_pulse t 6667 dose_nic$	None
Rule	stim_a7 stim_a7 <assignmentRule variable="stim_a7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 170 </cn> <cn type="integer"> 510 </cn> <ci> dose_aga7 </ci> </apply> </math> </assignmentRule>	=	$rect_pulse t 170510 dose_aga7$	None
Rule	stim_b2 stim_b2 <assignmentRule variable="stim_b2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> rect_pulse </ci> <ci> t </ci> <cn type="integer"> 60 </cn> <cn type="integer"> 61 </cn> <cn type="integer"> 0 </cn> </apply> </math> </assignmentRule>	=	$rect_pulse t 60610$	None
Rule	ss_ib2 ss_ib2 <assignmentRule variable="ss_ib2"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> comphill </ci> <ci> bas_ach </ci> <cn type="integer"> 0 </cn> <cn type="integer"> 0 </cn> <cn type="integer"> 1 </cn> <ci> w_actnicb2 </ci> <ci> w_actagb2 </ci> <ci> ec50_b2ach </ci> <ci> ec50_b2nic </ci> <ci> ec50_b2ag </ci> <ci> hill_actb2 </ci> </apply> </math> </assignmentRule>	=	$comphill bas_ach 001 w_actnicb2 w_actagb2 ec50_b2ach ec50_b2nic ec50_b2ag hill_actb2$	None
Rule	ss_pa7 ss_pa7 <assignmentRule variable="ss_pa7"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> comphill </ci> <apply> <plus/> <ci> bas_ach </ci> <ci> bas_chol </ci> </apply> <cn type="integer"> 0 </cn> <cn type="integer"> 0 </cn> <cn type="integer"> 1 </cn> <ci> w_actnica7 </ci> <ci> w_actaga7 </ci> <ci> ec50_a7ach </ci> <ci> ec50_a7nic </ci> <ci> ec50_a7ag </ci> <ci> hill_acta7 </ci> </apply> </math> </assignmentRule>	=	$comphill bas_ach bas_chol 001 w_actnica7 w_actaga7 ec50_a7ach ec50_a7nic ec50_a7ag hill_acta7$	None
Rule	ss_pglu ss_pglu <assignmentRule variable="ss_pglu"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <ci> ss_pa7 </ci> </math> </assignmentRule>	=	$ss_pa7$	None
Rule	ss_iglu ss_iglu <assignmentRule variable="ss_iglu"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <times/> <ci> w_glu </ci> <apply> <ci> clip </ci> <apply> <plus/> <ci> v_glu </ci> <ci> ss_pglu </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$w_glu clip v_glu ss_pglu$	None
Rule	ss_vgab ss_vgab <assignmentRule variable="ss_vgab"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <ci> hwr </ci> <apply> <plus/> <apply> <times/> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> s </ci> </apply> <ci> ss_iglu </ci> </apply> <apply> <times/> <apply> <minus/> <cn type="integer"> 1 </cn> <ci> r </ci> </apply> <ci> ss_ib2 </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$hwr 1 s ss_iglu 1 r ss_ib2$	None
Rule	ss_igab ss_igab <assignmentRule variable="ss_igab"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <times/> <ci> w_gab </ci> <apply> <ci> clip </ci> <apply> <plus/> <ci> ss_vgab </ci> <ci> i_basgab </ci> </apply> </apply> </apply> </math> </assignmentRule>	=	$w_gab clip ss_vgab i_basgab$	None
Rule	ss_vdop ss_vdop <assignmentRule variable="ss_vdop"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <apply> <plus/> <apply> <minus/> <ci> ss_igab </ci> </apply> <apply> <ci> hwr </ci> <apply> <plus/> <ci> i_basdop </ci> <apply> <times/> <ci> s </ci> <ci> ss_iglu </ci> </apply> <apply> <times/> <ci> r </ci> <ci> ss_ib2 </ci> </apply> </apply> </apply> </apply> </math> </assignmentRule>	=	$ss_igab hwr i_basdop s ss_iglu r ss_ib2$	None
Rule	t t <assignmentRule variable="t"> <math xmlns="http://www.w3.org/1998/Math/MathML"> <csymbol encoding="text" definitionURL="http://www.sbml.org/sbml/symbols/time"> time </csymbol> </math> </assignmentRule>	=	$time$	None