The provided code is an implementation of a modified Morris-Lecar model designed to simulate the electrophysiological behavior of neurons, with a particular focus on pyramidal neurons as noted in the Prescott et al. study. This model is an extension of the classic Morris-Lecar model, a two-dimensional model originally developed to describe the oscillatory behavior of barnacle muscle fibers. The modifications in this version are geared towards capturing the dynamics of pyramidal neurons under in vitro and in vivo conditions.
Fast Inward Sodium Current (INa):
minf(V)
.vna
) is representative of the potential where sodium would be in equilibrium, a typical characteristic of action potentials.Delayed Rectifier Potassium Current (IKdr):
w
, resembling the gating behavior of potassium channels.vk
) reflects the typical equilibrium potential for potassium ions, important for the return to resting potential.M-type Potassium Current:
zM
) is a slowly activating, voltage-dependent potassium current that contributes to subthreshold oscillations and neuronal excitability.tauzM
, betazM
, and gammazM
.Calcium-Activated Afterhyperpolarization (AHP) Current:
zAHP
) is modeled to activate during action potentials.Leak or Shunt Current (Ishunt):
gshunt
), the shunt current maintains the resting membrane potential.Noise:
This model aims to capture the dynamic properties of pyramidal neurons, such as spike generation, potential oscillations, and different firing regimes. By doing so, it provides a mechanistic understanding of how these neurons might behave under various physiological conditions, aligning with the known electrophysiological behaviors of pyramidal cells in both isolated and intact systems.