The code provided represents a computational model of an NMDA receptor-mediated synapse, using a two-state kinetic scheme to describe the synaptic conductance. These biological aspects are modeled in detail:
NMDA Receptors:
mgblock
function, which models the relief of the magnesium block with depolarization.Kinetic Scheme:
A
and G
to describe the transition of the synaptic state: rise (taur
) and decay (taud
) times of synaptic currents.Glutamate as a Neurotransmitter:
gluti
) is the key neurotransmitter driving NMDA receptor activation, represented by its binding dynamics and impact on synaptic conductance.Ionic Current Calculations:
ina
, ik
, ica
) flowing through the receptor channel due to ionic gradients of Na⁺, K⁺, and Ca²⁺, using the Goldman-Hodgkin-Katz (GHK) current equation with the ghk
function.cai
, cao
, nai
, nao
, ki
, ko
) and utilizing valence information to determine the flow of ions given the membrane potential.Voltage Dependence and Temperature:
mgblock
function considers the dependence on voltage (v
) and extracellular magnesium concentration (mg
), a characteristic of NMDA receptor function.R
and considerations for body temperature (celsius
).Pmax
, the model can experiment with different scenarios of synaptic strength and explore conditions that lead to long-term potentiation (LTP) or depression (LTD).This model provides a detailed simulation of NMDA receptor behavior and its contribution to synaptic transmission, encapsulating the receptive dynamics to biochemical and physiological variables critical for neuronal communication.