Theoretical model of the interactions between Ca2+, calmodulin and myosin light chain kinase

Active Ca²⁺/calmodulin (CaM)-dependent myosin light chain kinase (MLCK) plays an important role in the process of MLC phosphorylation and consecutive smooth muscle contraction. Here, we propose a mathematical model of a detailed kinetic scheme describing interactions among Ca²⁺, CaM and MLCK and taking into account eight different aggregates. The main model result is the prediction of the Ca²⁺ dependent active form of MLCK, which is in the model taken as proportional to the concentration of Ca₄CaM · MLCK complex. Wegscheider’s condition is additionally applied as a constraint enabling the prediction of some parameter values that have not yet been obtained by experiments.