/* * Role for G protein G-beta-gamma isoform specificity in synaptic * signal processing: a computational study * * Model Status * * This CellML model represents the post-synaptic cell. The model * runs in both OpenCell and COR, but the membrane only depolarises * to -5 mV once as oposed to repeatedly depolarising to ~40 mV. * The lack of repeated spikes is due to there being no simulus * input. The post-synaptic model has to be coupled to the pre-synaptic * for this to occur. The units of the CellML model have been checked * and there are no inconsistencies (just equivalences). * * Model Structure * * ABSTRACT: Computational modeling is used to investigate the * functional impact of G protein-mediated presynaptic autoinhibition * on synaptic filtering properties. It is demonstrated that this * form of autoinhibition, which is relieved by depolarization, * acts as a high-pass filter. This contrasts with vesicle depletion, * which acts as a low-pass filter. Model parameters are adjusted * to reproduce kinetic slowing data from different Gbetagamma * dimeric isoforms, which produce different degrees of slowing. * With these sets of parameter values, we demonstrate that the * range of frequencies filtered out by the autoinhibition varies * greatly depending on the Gbetagamma isoform activated by the * autoreceptors. It is shown that G protein autoinhibition can * enhance the spatial contrast between a spatially distributed * high-frequency signal and surrounding low-frequency noise, providing * an alternate mechanism to lateral inhibition. It is also shown * that autoinhibition can increase the fidelity of coincidence * detection by increasing the signal-to-noise ratio in the postsynaptic * cell. The filter cut, the input frequency below which signals * are filtered, depends on several biophysical parameters in addition * to those related to Gbetagamma binding and unbinding. By varying * one such parameter, the rate at which transmitter unbinds from * autoreceptors, we show that the filter cut can be adjusted up * or down for several of the Gbetagamma isoforms. This allows * for great synapse-to-synapse variability in the distinction * between signal and noise. * * The original paper reference is cited below: * * Role for G Protein G-Beta-Gamma Isoform Specificity in Synaptic * Signal Processing: A Computational Study, Richard Bertram, Michelle * I. Arnot, and Gerald W. Zamponi, 2002,Journal of Neurophysiology * , 87, 2612-2623. PubMed ID: 11976397 * * reaction diagram * * [[Image file: bertram_2002.png]] * * Schematic diagram of the presynaptic model. */ import nsrunit; unit conversion on; unit s=1 second^1; unit ms=.001 second^1; unit uM=1E-3 meter^(-3)*mole^1; unit uA=1E-6 ampere^1; unit pS=1E-12 kilogram^(-1)*meter^(-2)*second^3*ampere^2; unit mV=.001 kilogram^1*meter^2*second^(-3)*ampere^(-1); unit cm2=1E-4 meter^2; unit per_cm2=1E4 meter^(-2); unit per_ms=1E3 second^(-1); unit uM_per_ms=1 meter^(-3)*second^(-1)*mole^1; unit per_uM_per_ms=1E6 meter^3*second^(-1)*mole^(-1); unit per_mV_per_ms=1E6 kilogram^(-1)*meter^(-2)*second^2*ampere^1; unit uM2_per_ms=.001 meter^(-6)*second^(-1)*mole^2; unit mV_per_uM=1 kilogram^1*meter^5*second^(-3)*ampere^(-1)*mole^(-1); unit uF_per_cm2=.01 kilogram^(-1)*meter^(-4)*second^4*ampere^2; unit uA_per_cm2=.01 meter^(-2)*ampere^1; unit pS_per_cm2=1E-8 kilogram^(-1)*meter^(-4)*second^3*ampere^2; unit mS_per_cm2=10 kilogram^(-1)*meter^(-4)*second^3*ampere^2; // unit nanometre predefined unit uM2_nanometre=1E-15 meter^(-5)*mole^2; unit uM_per_cm2=10 meter^(-5)*mole^1; unit millijoule_per_mole_kelvin=.001 kilogram^1*meter^2*second^(-2)*kelvin^(-1)*mole^(-1); unit coulomb_per_mole=1 second^1*ampere^1*mole^(-1); unit uM_per_coulomb=1E-3 meter^(-3)*second^(-1)*ampere^(-1)*mole^1; unit micrometre2_per_second=1E-12 meter^2*second^(-1); unit per_micrometre3=1E18 meter^(-3); math main { realDomain time ms; time.min=0; extern time.max; extern time.delta; real R_gas_const millijoule_per_mole_kelvin; R_gas_const=8314.41; real Temp kelvin; Temp=310; real F coulomb_per_mole; F=96485; real R(time) dimensionless; when(time=time.min) R=0; real kr_plus per_uM_per_ms; kr_plus=0.15; real kr_minus per_ms; kr_minus=2.5; real Ca(time) uM; real T(time) uM; real T_bar uM; T_bar=4000.0; real Ca_ex uM; Ca_ex=2000.0; real Ca_open(time) uM; real Dc micrometre2_per_second; Dc=220; real r nanometre; r=10; real sigma(time) uM_per_ms; real i_V(time) uA; real g_Ca pS; g_Ca=1.2; real P mV_per_uM; P=0.006; real V_post(time) mV; when(time=time.min) V_post=-65; real O(time) dimensionless; real alpha(time) per_ms; real alpha_(time) per_ms; real beta(time) per_ms; real beta_(time) per_ms; real kG_plus(time) per_ms; real b(time) dimensionless; when(time=time.min) b=0; real kb_plus per_uM_per_ms; kb_plus=2000.0; real kb_minus per_ms; kb_minus=1.0; real kG_minus per_ms; kG_minus=0.00025; real kG2_minus per_ms; kG2_minus=0.016; real kG3_minus per_ms; kG3_minus=1.024; real C1(time) dimensionless; when(time=time.min) C1=1; real C2(time) dimensionless; when(time=time.min) C2=0; real C_G1(time) dimensionless; when(time=time.min) C_G1=0; real C3(time) dimensionless; when(time=time.min) C3=0; real C_G2(time) dimensionless; when(time=time.min) C_G2=0; real C4(time) dimensionless; when(time=time.min) C4=0; real C_G3(time) dimensionless; when(time=time.min) C_G3=0; real C_G(time) dimensionless; real Cm uF_per_cm2; Cm=1.0; real i_syn(time) uA_per_cm2; real i_Na_post(time) uA_per_cm2; real i_K_post(time) uA_per_cm2; real i_leak_post(time) uA_per_cm2; real g_syn mS_per_cm2; g_syn=0.2; real V_syn mV; V_syn=0; real x_infinity(time) dimensionless; real alpha_x(time) dimensionless; real beta_x(time) dimensionless; real n_post(time) dimensionless; when(time=time.min) n_post=0; real alpha_n(time) per_ms; real beta_n(time) per_ms; // // // R:time=(kr_plus*Ca*(1-R)-kr_minus*R); T=(T_bar*R); // Ca_open=(sigma/((2 per_micrometre3)*Dc*r*3.141592653589793)); sigma=((-1)*(5.182 uM_per_coulomb)*i_V); i_V=(g_Ca*P*2*F*V_post/(R_gas_const*Temp)*Ca_ex/(1-exp(2*F*V_post/(R_gas_const*Temp)))); Ca=(O*Ca_open+(.1 uM)); // alpha=((.45 per_ms)*exp(V_post/(22 mV))); alpha_=(alpha/8); beta=((.015 per_ms)*exp((-1)*V_post/(14 mV))); beta_=(beta*8); b:time=(kb_plus*T*(1-b)-kb_minus*b); kG_plus=((3 per_ms)*b/(680+320*b)); // C1:time=(beta*C2+kG_minus*C_G1-C1*(4*alpha+kG_plus)); // C2:time=(4*alpha*C1+2*beta*C3+kG2_minus*C_G2-C2*(beta+3*alpha+kG_plus)); // C3:time=(3*alpha*C2+3*beta*C4+kG3_minus*C_G3-C3*(2*beta+2*alpha+kG_plus)); // C4:time=(2*alpha*C3+4*beta*O-C4*(3*beta+alpha)); // O=(1-C1-C2-C3-C4-C_G1-C_G2-C_G3); C_G=(C_G1+C_G2+C_G3); // C_G1:time=(beta_*C_G2+kG_plus*C1-C_G1*(4*alpha_+kG_minus)); // C_G2:time=(4*alpha_*C_G1+2*beta_*C_G3+kG_plus*C2-C_G2*(beta_+3*alpha_+kG2_minus)); // C_G3:time=(3*alpha_*C_G2+kG_plus*C3-C_G3*(2*beta_+kG3_minus)); // V_post:time=((-1)*(i_Na_post+i_K_post+i_leak_post+i_syn)/Cm); // i_syn=(g_syn*b*(V_post-V_syn)); // i_leak_post=((.3 pS_per_cm2)*(V_post+(54 mV))); // i_Na_post=((120 pS_per_cm2)*x_infinity^3*(1-n_post)*(V_post-(120 mV))); x_infinity=(alpha_x/(alpha_x+beta_x)); alpha_x=(.2*(V_post+(40 mV))/((1 mV)-(1 mV)*exp((-1)*(V_post+(40 mV))/(10 mV)))); beta_x=(8*exp((1 mV)/((-1)*(V_post+(65 mV)/18)))); // i_K_post=((36 pS_per_cm2)*n_post^4*(V_post+(77 mV))); // n_post:time=(alpha_n*(1-n_post)-beta_n*n_post); alpha_n=((.02 per_mV_per_ms)*(V_post+(55 mV))/(1-exp((-1)*(V_post+(55 mV))/(10 mV)))); beta_n=((.25 per_ms)*exp((-1)*(V_post+(65 mV))/(80 mV))); }