GABA_BSynapse.h

Go to the documentation of this file.

#ifndef _GABA_BSYNAPSE_H_
#define _GABA_BSYNAPSE_H_

#include "GABA_BResponse.h"
#include "GABA_BResponseODE.h"

#include "GenericStaticSpikingSynapse.h"
#include "GenericConductanceBasedSpikingSynapse.h"
// NEEDS TO BE TESTED


/* DATA from Alain Destexhe

Cmax_GABAb = 1      // (mM) max transmitter concentration
Cdur_GABAb = 1      // (ms) transmitter duration (rising phase)
K1_GABAb = 0.09     // (/ms mM) forward binding rate to receptor
K2_GABAb = 0.0012   // (/ms) backward (unbinding) rate of receptor
K3_GABAb = 0.18     // (/ms) rate of G-protein production
K4_GABAb = 0.034    // (/ms) rate of G-protein decay
KD_GABAb = 100      // dissociation constant of K+ channel
n_GABAb  = 4        // nb of binding sites of G-protein on K+
Erev_GABAb = -95    // (mV) reversal potential (E_K)
Prethresh_GABAb = 0 // (mV) voltage level nec for release
Deadtime_GABAb = 1  // (ms) mimimum time between release events
gmax  = 0.0001      // (umho) maximum conductance 
*/


// with ODE Solver

class GABA_BSynapseODE : public GenericStaticSpikingSynapse<
                                GenericConductanceBasedSpikingSynapse<
                                GABA_BResponseODE > >
{
    SIMOBJECT(GABA_BSynapseODE, AdvancePhase::SpikeDriven)

public:
    GABA_BSynapseODE(
        const float W         = 1e-9,
        const float delay     = 2e-4,
        const float Erev      = -90e-3,
        const float k1        = 0.09*10e4,   // [M^-1 s^-1]
        const float k2        = 1.2,         // [s^-1]
        const float k3        = 180.0,       // [s^-1]
        const float k4        = 34.0,        // [s^-1]
        const float kd        = 100.0,
        const float dead_time = 1e-3,
        const float Cdur      = 1.0e-3,      // [s]
        const float Cmax      = 0.5e-3,      // [M]
        const int   n         = 4)
    {
        this->W = W;
        this->delay = delay;
        this->Erev = Erev;

        this->g_max = W;
        this->k1=k1;
        this->k2=k2;
        this->k3=k3;
        this->k4=k4;
        this->kd=kd;
        this->dead_time=dead_time;
        this->Cdur=Cdur;
        this->Cmax=Cmax;
        this->n=n;
    };

};

// euler


class GABA_BSynapse : public GenericStaticSpikingSynapse<
                             GenericConductanceBasedSpikingSynapse<
                             GABA_BResponse > >
{
    SIMOBJECT(GABA_BSynapse, AdvancePhase::SpikeDriven)

public:
    GABA_BSynapse(
        const float W         = 1e-9,
        const float delay     = 2e-4,
        const float Erev      = -90e-3,
        const float k1        = 0.09*10e6,   // [M^-1 s^-1]
        const float k2        = 1.2,         // [s^-1]
        const float k3        = 180.0,       // [s^-1]
        const float k4        = 34.0,        // [s^-1]
        const float kd        = 100.0,
        const float dead_time = 1e-3,
        const float Cdur      = 0.5e-3,      // [s]
        const float Cmax      = 1.0e-3,      // [M]
        const int   n         = 4)
    {
        this->W = W;
        this->delay = delay;
        this->Erev = Erev;

        this->g_max = W;
        this->k1=k1;
        this->k2=k2;
        this->k3=k3;
        this->k4=k4;
        this->kd=kd;
        this->dead_time=dead_time;
        this->Cmax=Cmax;
        this->Cdur=Cdur;
        this->n=n;
    };

};




// Nest style


class GABA_BSynapseNest : public GenericStaticSpikingSynapse<
                                 GenericConductanceBasedSpikingSynapse<
                                 GABA_BResponseNest > >
{
    SIMOBJECT(GABA_BSynapseNest, AdvancePhase::SpikeDriven)

public:
    GABA_BSynapseNest(
        const float W     = 1e-9,
        const float delay = 2e-4,
        const float Erev  = -90e-3,

        const float k1    = 0.09*10e6,   // [M^-1 s^-1]
        const float k2    = 1.2,         // [s^-1]
        const float k3    = 180.0,       // [s^-1]
        const float k4    = 34.0,        // [s^-1]
        const float kd    = 100.0,       // [mM]
        const float pulse = 1.105*10e-3, // [s]
        const float T     = 1e-3,        // [mM]
        const int   n     = 4)
    {
        this->W = W;
        this->delay = delay;
        this->Erev = Erev;

        this->g_max = W;
        this->k1=k1;
        this->k2=k2;
        this->k3=k3;
        this->k4=k4;
        this->kd=kd;
        this->T=T;
        this->pulse=pulse;
        this->n=n;
    };

};


#endif //_GABA_BSYNAPSE_H_

---