DiffEqs

From Sean_Carver
Revision as of 22:54, 10 February 2009 by Carver (talk | contribs) (Add Leak Channels)
Jump to: navigation, search

Plain Bilayer Membrane

Remember the capacitor equation:

  C \frac{dV(t)}{dt} = I(t)

 V(0) = 0

Here C is a constant (a parameter): the capacitance of the membrane.

Next V(t) is the membrane potential. Typically V(t) is the quantity of interest -- and it is the quantity that is usually measure.

Next I(t) is the capacitive current into the cell. In general the capacitive current involves the ionic currents into the cell as well as the injected current, but for a plain lipid bilayer, it is just the injected current.

Finally V(0) is the initial membrane potential of the cell, the constant of integration required for a definite solution to the equation.

We can solve this equation

 V(t) = V(0) + \frac{1}{C} \int_0^t I(s) ds

Add Na+/K+ Pump

It seems clear from the Neurons in Action Tutorial that they modeled the addition of the pump as a change in the initial condition of the cell.

  C \frac{dV(t)}{dt} = I(t)

 V(0) = -75

 V(t) = V(0) + \frac{1}{C} \int_0^t I(s) ds

Add Leak Channels

 C \frac{dV}{dt} = I(t) - g_L (V - E_L)

Leak channels add a resistor in parallel to the capacitor.

 g_L is the conductance of the channels (which is the reciprocal of the resistance of the channels).

 E_L is the reversal potential of the leak current.

 V - E_L is the driving potential.

 g_L (V - E_L) is the current across the resistor (Ohm's Law).

The resistive current is minus its contribution to the capacitive current (Kirchov's Law).

Add HH Channels

 C \frac{dV}{dt} = I(t) - g_L (V(t) - E_L) - g_{Na} m^3(t) h(t) (V(t) - E_{Na}) - g_K n^4(t) (V(t) - E_K)