A dynamical systems model for the total fission rate in Drp1-dependent mitochondrial fission

Mitochondrial hyperfission in response to cellular insult is associated with reduced energy production and programmed cell death. Thus, there is a critical need to understand the molecular mechanisms coordinating and regulating the complex process of mitochondrial fission. We develop a nonlinear dynamical systems model of dynamin related protein one (Drp1)-dependent mitochondrial fission and use it to identify parameters which can regulate the total fission rate (TFR) as a function of time. The TFR defined from a nondimensionalization of the model undergoes a Hopf bifurcation with bifurcation parameter [Formula: see text] where [Formula: see text] is the total concentration of mitochondrial fission factor (Mff) and k+ and k- are the association and dissociation rate constants between oligomers on the outer mitochondrial membrane. The variable μ can be thought of as the maximum build rate over the disassembling rate of oligomers. Though the nondimensionalization of the system results in four dimensionless parameters, we found the TFR and the cumulative total fission (TF) depend strongly on only one, μ. Interestingly, the cumulative TF does not monotonically increase as μ increases. Instead it increases with μ to a certain point and then begins to decrease as μ continues to increase. This non-monotone dependence on μ suggests interventions targeting k+, k-, or [Formula: see text] may have a non-intuitive impact on the fission mechanism. Thus understanding the impact of regulatory parameters, such as μ, may assist future therapeutic target selection.