
MBW:The MorrisLecar Model for Excitable SystemsFrom MathBioThis is a summary of the paper Morris, C. and Lecar, H. (1981) Voltage oscillations in the barnacle giant muscle fiber. Biophys. J. 35: 193  213.
ContentsOverviewThe biological system studied was the Barnacle muscle through its membrane potential and the concentrations which had an effect on it. This model used a simplified version of the HodgkinHuxley model. A reduced system of equations was created that produce the same results for most cases. In order to approximate the muscle reactions appropriately, the models are based off of one you would use to follow an equivalent circuit design. For more on the HodgkinHuxley model, please visit MBW:HodgkinHuxley. IntroductionIn 1981, Morris and Lecar proposed their eponymous equation as an reduced excitation model in the vein of the HogdkinHuxley equations, which constitutes an exact description of the conductancebased model for the giant squid neuron. In particular, they begin with an "already reduced" system of equations taken from the voltage clamp studies of Keynes, et. all [1], which proposes that the relevant state variables for the barnacle muscle are, along with the membrane potential, the proportion of voltage dependent Ca+ and K+ ion channels open at a given time (which is to say, the Na+ conduction plays no significant role in dynamics):
From here, they make various assumptions about independence (or uncoupling of their dynamics) of the Ca++ and K+ channels, as well as assumptions about their relative recovery rates in order to further reduce the dimension of the model while preserving the essential character of the system's long term behavior. While this reduced model exhibits a rich variety of behaviors, not all of the possible behavior observed in excitable systems, such as bursting and chaos, can described it. Parameter list and A Qualitative Discussion of the Nonreduced modelParameter List= applied current () = leak, Ca++, and K+ currents, respectively () = maximum or instantaneous conductance values for leak, Ca++, and K+ pathways,respectively () = conductance constant for nonlinear () = membrane potential = equilibrium potential corresponding to leak, Ca++, and K+ conductances, respectively (mV) = fraction of open Ca++ channels = fraction of open K+ channels = fraction of open Ca++ and K+ channels, at steady state = rate constant for opening of Ca++ and K+ channels = maximum rate constants for Ca++ and K+ channel opening = potential at which = reciprocal of slope of voltage dependence of Mz. () = potential at which = reciprocal of slope of voltage dependence of NCO () = membrane capacitance () where
. Discussion of nonreduced modelIn order to see how equations (1) might be an intuitive model for voltage dependence ion channel, fix . Consider the 2nd (1); if , then we have that , so that the solution is decreasing. On the other hand, if , then we have that , so that the solution is increasing. Thus, is a stable equilibrium as expected. A similar sequence of ideas applies to , and we can conclude that in absence of voltage fluctuation, the number Ca++ and K+ ion channel tends to a constant determined by its Boltzmann distribution () for the potential.
Now consider nonconstant voltage; in light, of above we can think of and "following" the moving voltage. However, since the voltage itself depends on the number of open Ca++ and K+ channels, this ODE describes a feed back loop whereby the voltage gradient increases as if the net number of open ion pathways is decreasing, while an increasing voltage forces more ion channels to open and produces more current which depletes the voltage once V becomes greater than or , creating a kind of pushpull, ebb and flow mechanism. It is easy to see that this pushpull mechanism is partially dependent (in addition to the the applied current) on the values of the Calcium and Potassium conduction, and . This along with the fact that the Ca++ K+ have opposite signs in the equation for voltage, suggests: given a large enough ratio of relative to (and with an adequate amount of applied current) sustained or dampened oscillations might occur. Indeed this case, as we shall see in Existence of Limit Cycles section. The ExperimentMorris and Lecar employed a technique known as current clamping and space clamping, which essentially eliminated spatial variations the traveling action potential and fixed the current across the membrane to set parameter which the experimenter could control [2]. Then in order to reduce the model, they submerged the barnacle muscle in various solutions in order to eliminate either Ca+ and Na+ conductance, or in the case of the MorrisLecar equation, simultaneously optimize them. By doing this, they reduced the dynamics of the system to two state variables (depending on which solution the barnacle is submersed in), which considerably simplifies possible behavior the system can exhibit due to special topological properties of the phase plane that is not enjoy by higher dimensional Euclidean space (namely that a simple closed curve divides the plane into two distinct connected components. See [1] [2]. Pure K+ conductanceWhen the barnacle is submersed in a Ca+ free solution, Ca+ conductance is eliminates and eq (1) reduces to:
This system has precisely one fixedpoint which is always an attractor. Here is a typical orbit for this system: Pure Ca+ condunctanceLikewise, when the barnacle is submersed in the solution that minimizes the K+ conductance, we have
The situation for pure Ca+ conductance is more complicated. Here, in addition to an attracting fixedpoint, the fixed point can be a saddle point and a second stable fixed point appears. See the section Bifurcations for details. MorrisLecar equationThe main idea behind the derivation of the MorrisLecar equation is that relative to the K+ ion channel, the recovery of Ca++ state variable is instaneous, i.e. . The we obtain
Existence of Limit CyclesMorris and Lecar prove the existence of a limit cycle for certain values of the applied current and the ion conductances and , in the cannonical way. That is, they construct a trapping region which contains no fixed point, and then apply the PoincareBendixson theorem. Their argument is reproduced below: First, we note that, intuitively, this physical system must have finite energy, and therefore is bounded in some sense. To see this, we note the state variable N, can only takes on value between 0 and 1, and thus remains in this interval. Likewise, given a fixed applied current, the analogous circuit cannot exceed the total energy within it, and is bounded by
which is derived from the fact that when or , and thus must be bounded by these two limiting potentials. Thus we have shown that any orbit with the rectangle must remain in . Second, we linearize system (4) by finding its Jacobian matrix of partial derivatives and evaluate it at the fixed point (which is contained in the rectangle T, as shown by the nullclines ). From here, we obtain the characteristic polynomial, the solutions of which are the eigenvalues of the linearized ODE: . In order to have a stable limit cycle, both roots must be positive if they are real or have a positive real part if they are complex. Thus we could apply the HartmanGrobman theorem, and conclude that the MorrisLecar equation is conjugate (call it ) to the unstable linearization. Thus, say if (4) is conjugate in a neighborhood of the equilibrium, then any initial condition (for the linearized system) on the closed curve can never enter the region . Therefore, forms the inner boundary of our trapping region for equation 4, and hence created a trapping region which contains no fixed points. Thus we can conclude by PoincareBendixson that orbits with in this region must approach a periodic orbit. In order to guarantee, that the real parts of the eigenvalues of the linearized system are positive, we must have that: + . Substituting equation (4) into the above we get ,
which gives the sufficient conditions, and defines a parameter space, under which limit cycles occur. BifurcationsFor fixed , as the applied voltage increased past certain thresholds, the MorrisLecar equation can under go a series of bifurcation which exhibits the entire range of behavior discussed above:
Also, for a fixed applied current, we can construct a bifurcation diagram in terms of the ion conductance and : Here we can see a transition from the stable node behavior and the bistable saddle Notice that as increases beyond the bifurcation point (for ). two stable equilibrium appear, each basin of which having a common boundary of the separatrix formed by the stable manifold of the saddle point equilibrium. ConclusionThe MorrisLecar model for excitable systems is a simplified form of the HodkinsHuxley equations that captures much of the behavior of the higherorder system. By reducing the system to two state variables, analysis become more tenable and the theory of bifurcations in the plane can applied in order to rigorously show the existence of mono and bistable orbits, Limit cycles, and dampened oscillations for certain parameters. See also, The Behavior of Traveling Waves in excitable media Citations[1] Keynes, R. D., E. Rojas, R. E. Taylor, and J. Vergara. 1973. Calcium and potassium systems of a giant barnacle muscle fibre under membrane potential control. J. Physiol. (Lond.) 229:409455. [2] N. Britton, Essential Mathematical Biology. (2003) 