Abstract
We present a simplified reaction network in a single well-mixed volume that captures the general features of CaMKII dynamics observed during both synaptic input and spine depolarization. Our model can also account for the greater-than-control CaMKII activation observed with added EGTA during depolarization. Calcium input currents are modeled after experimental observations, and existing models of calmodulin and CaMKII autophosphorylation are used. After calibration against CaMKII activation data in the absence of chelators, CaMKII activation dynamics due to synaptic input via n-methyl-d-aspartate receptors are qualitatively accounted for in the presence of the chelators EGTA and BAPTA without additional adjustments to the model. To account for CaMKII activation dynamics during spine depolarization with added EGTA or BAPTA, the model invokes the modulation of CaV2.3 (R-type) voltage-dependent calcium channel (VDCC) currents observed in the presence of EGTA or BAPTA. To our knowledge, this is a novel explanation for the increased CaMKII activation seen in dendritic spines with added EGTA, and suggests that differential modulation of VDCCs by EGTA and BAPTA offers an alternative or complementary explanation for other experimental results in which addition of EGTA or BAPTA produces different effects. Our results also show that a simplified reaction network in a single, well-mixed compartment is sufficient to account for the general features of observed CaMKII dynamics.