# American Institute of Mathematical Sciences

July  2017, 37(7): 4131-4158. doi: 10.3934/dcds.2017176

## Mathematical analysis of an in vivo model of mitochondrial swelling

 1 Institute for Computational Biology, Helmholtz Zentrum München, Ingolstäder Landstr. 1, 85764 Neuherberg, Germany 2 Department of Applied Phsyics, School of Science and Engineering, Waseda University 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-855, Japan 3 Department of Mathematics and Statistics, University of Guelph, Guelph ON, N1G2W1, Canada

* Corresponding author: Messoud Efendiev

Received  October 2016 Revised  February 2017 Published  April 2017

Fund Project: M.O. is partly supported by the Grant-in-Aid for Scientific Research #15K13451, the Ministry of Education, Culture, Sports, Science, and Technology, Japan; H.J.E. is partly supported by the Natrural Science and Engineering Researc Council of Canada through a Discovery Grant.

We analyze the effect of Robin boundary conditions in a mathematical model for a mitochondria swelling in a living organism. This is a coupled PDE/ODE model for the dependent variables calcium ion contration and three fractions of mitochondria that are distinguished by their state of swelling activity. The model assumes that the boundary is a permeable 'membrane', through which calcium ions can both enter or leave the cell. Under biologically relevant assumptions on the data, we prove the well-posedness of solutions of the model and study the asymptotic behavior of its solutions. We augment the analysis of the model with computer simulations that illustrate the theoretically obtained results.

Citation: Messoud Efendiev, Mitsuharu Ôtani, Hermann J. Eberl. Mathematical analysis of an in vivo model of mitochondrial swelling. Discrete & Continuous Dynamical Systems - A, 2017, 37 (7) : 4131-4158. doi: 10.3934/dcds.2017176
##### References:
 [1] S. Brenner, Poincaré-Friedrichs inequalities for piecewise $H^1$ functions, SIAM J. Numer. Anal., 41 (2003), 306-324. doi: 10.1137/S0036142902401311. Google Scholar [2] H. Brézis, Opérateurs Maximaux Monotones et Semigroupes de Contractions dans un Espace de Hilbert, North Holland, Amsterdam, The Netherlands, 1973.Google Scholar [3] H. Brézis, Monotonicity methods in Hilbert spaces and some applications to nonlinear partial differential equations, Contributions to Nonlinear Functional Analysis (ed. E. H. Zarantonello), Academic Press, (1971), 101-179. Google Scholar [4] M. A. Efendiev, M. Ôtani and H. J. Eberl, A coupled PDE/ODE model of mitochondrial swelling: Large-time behavior of homogeneous Dirichlet problem, Journal of Coupled Systems and Multiscale Dynamics, 3 (2015), 1-13. doi: 10.1166/jcsmd.2015.1070. Google Scholar [5] S. Eisenhofer, A coupled system of ordinary and partial differential equations modeling the swelling of mitochondria, PhD Dissertation, TU Munich, 2013.Google Scholar [6] S. Eisenhofer, M. A. Efendiev, M. Ôtani, S. Schulz and H. Zischka, On a ODE-PDE coupling model of the mitochondrial swelling process, Discrete and Continuous Dynamical Syst. Ser. B, 20 (2015), 1031-1057. doi: 10.3934/dcdsb.2015.20.1031. Google Scholar [7] S. Eisenhofer, F. Toókos, B. A. Hense, S. Schulz, F. Filbir and H. Zischka, A mathematical model of mitochondrial swelling BMC Research Notes, 3 (2010), p67. doi: 10.1186/1756-0500-3-67. Google Scholar [8] G. Kroemer, L. Galluzzi and C. Brenner, Mitochondrial membrane permeabilization in cell death, Physiological Reviews, 87 (2007), 99-163. Google Scholar [9] M. Ôtani, Nonmonotone perturbations for nonlinear parabolic equations associated with subdifferential operators, Cauchy problems, J. Differential Equations, 46 (1982), 268-299. doi: 10.1016/0022-0396(82)90119-X. Google Scholar [10] V. Petronilli, C. Cola, S. Massari, R. Colonna and P. Bernardi, Physiological effectors modify voltage sensing by the cyclosporin A-sensitive permeability transition pore of mitochondria, Journal of Biological Chemistry, 268 (1993), 21939-21945. Google Scholar [11] R. Rizzuto and T. Pozzan, Microdomains of intracellular $\textrm{Ca}^{2+}$: Molecular determinants and functional consequences, Physiological Reviews, 86 (2006), 369-408. Google Scholar [12] R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physics Springer-Verlag, New York, 1997. doi: 10.1007/978-1-4612-0645-3. Google Scholar [13] H. Triebel, Interpolation Theory, Function Spaces, Differential Operators, J. A. Barth, 1995. Google Scholar [14] H. Zischka, N. Larochette, F. Hoffmann, D. Hamöoller, N. Jägemann, J. Lichtmannegger, L. Jennen, J. Müller-Höcker, F. Roggel, M. Göttlicher, A. M. Vollmar and G. Kroemer, Electrophoretic analysis of the mitochondrial outer membrane rupture induced by permeability transition, Analytical Chemistry, 80 (2008), 5051-5058. doi: 10.1021/ac800173r. Google Scholar

show all references

##### References:
 [1] S. Brenner, Poincaré-Friedrichs inequalities for piecewise $H^1$ functions, SIAM J. Numer. Anal., 41 (2003), 306-324. doi: 10.1137/S0036142902401311. Google Scholar [2] H. Brézis, Opérateurs Maximaux Monotones et Semigroupes de Contractions dans un Espace de Hilbert, North Holland, Amsterdam, The Netherlands, 1973.Google Scholar [3] H. Brézis, Monotonicity methods in Hilbert spaces and some applications to nonlinear partial differential equations, Contributions to Nonlinear Functional Analysis (ed. E. H. Zarantonello), Academic Press, (1971), 101-179. Google Scholar [4] M. A. Efendiev, M. Ôtani and H. J. Eberl, A coupled PDE/ODE model of mitochondrial swelling: Large-time behavior of homogeneous Dirichlet problem, Journal of Coupled Systems and Multiscale Dynamics, 3 (2015), 1-13. doi: 10.1166/jcsmd.2015.1070. Google Scholar [5] S. Eisenhofer, A coupled system of ordinary and partial differential equations modeling the swelling of mitochondria, PhD Dissertation, TU Munich, 2013.Google Scholar [6] S. Eisenhofer, M. A. Efendiev, M. Ôtani, S. Schulz and H. Zischka, On a ODE-PDE coupling model of the mitochondrial swelling process, Discrete and Continuous Dynamical Syst. Ser. B, 20 (2015), 1031-1057. doi: 10.3934/dcdsb.2015.20.1031. Google Scholar [7] S. Eisenhofer, F. Toókos, B. A. Hense, S. Schulz, F. Filbir and H. Zischka, A mathematical model of mitochondrial swelling BMC Research Notes, 3 (2010), p67. doi: 10.1186/1756-0500-3-67. Google Scholar [8] G. Kroemer, L. Galluzzi and C. Brenner, Mitochondrial membrane permeabilization in cell death, Physiological Reviews, 87 (2007), 99-163. Google Scholar [9] M. Ôtani, Nonmonotone perturbations for nonlinear parabolic equations associated with subdifferential operators, Cauchy problems, J. Differential Equations, 46 (1982), 268-299. doi: 10.1016/0022-0396(82)90119-X. Google Scholar [10] V. Petronilli, C. Cola, S. Massari, R. Colonna and P. Bernardi, Physiological effectors modify voltage sensing by the cyclosporin A-sensitive permeability transition pore of mitochondria, Journal of Biological Chemistry, 268 (1993), 21939-21945. Google Scholar [11] R. Rizzuto and T. Pozzan, Microdomains of intracellular $\textrm{Ca}^{2+}$: Molecular determinants and functional consequences, Physiological Reviews, 86 (2006), 369-408. Google Scholar [12] R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physics Springer-Verlag, New York, 1997. doi: 10.1007/978-1-4612-0645-3. Google Scholar [13] H. Triebel, Interpolation Theory, Function Spaces, Differential Operators, J. A. Barth, 1995. Google Scholar [14] H. Zischka, N. Larochette, F. Hoffmann, D. Hamöoller, N. Jägemann, J. Lichtmannegger, L. Jennen, J. Müller-Höcker, F. Roggel, M. Göttlicher, A. M. Vollmar and G. Kroemer, Electrophoretic analysis of the mitochondrial outer membrane rupture induced by permeability transition, Analytical Chemistry, 80 (2008), 5051-5058. doi: 10.1021/ac800173r. Google Scholar
Model simulation with $\alpha=10<C^-$: Shown are $u, N_1, N_2, N_3$ for selected times.
Model simulation with $\alpha=10<C^-$: Shown is $N_1$ for selected times.
Simulation to illustrate partial swelling in Theorem 5.2, using initial data (ref{T2init:eq}): shown is the minimum value of $N_2$ as a function of time for different base calcium ion concentrations $u_{base}$ (top left), along with the steady state distributions for $N_1$ (top right), $N_2$ (bottom left), and $N_3$ (bottom right) in the case $u_{base}=100$.
Mitochondria populations $N_1$ and $N_2$ as a function of time in three points of the domain on a line through the center point: A (close to the boundary), B (half way between boundary and center), C (in the center), for six different values of the external calcium ion concentration $\alpha$.
Default parameter values, cf also [5]
 parameter symbol value remark lower (initiation) swelling threshold $C^-$ 20 (varied) upper (maximum) swelling threshold $C^+$ 200 maximum transition rate for $N_1\rightarrow N_2$ $f^\ast$ 1 maximum transition rate for $N_2\rightarrow N_3$ $g^\ast$ 1 diffusion coefficient $d_1$ 0.2 (varied) feedback parameter $d_2$ 30
 parameter symbol value remark lower (initiation) swelling threshold $C^-$ 20 (varied) upper (maximum) swelling threshold $C^+$ 200 maximum transition rate for $N_1\rightarrow N_2$ $f^\ast$ 1 maximum transition rate for $N_2\rightarrow N_3$ $g^\ast$ 1 diffusion coefficient $d_1$ 0.2 (varied) feedback parameter $d_2$ 30
 [1] Sabine Eisenhofer, Messoud A. Efendiev, Mitsuharu Ôtani, Sabine Schulz, Hans Zischka. On an ODE-PDE coupling model of the mitochondrial swelling process. Discrete & Continuous Dynamical Systems - B, 2015, 20 (4) : 1031-1057. doi: 10.3934/dcdsb.2015.20.1031 [2] Mustapha Cheggag, Angelo Favini, Rabah Labbas, Stéphane Maingot, Ahmed Medeghri. Complete abstract differential equations of elliptic type with general Robin boundary conditions, in UMD spaces. Discrete & Continuous Dynamical Systems - S, 2011, 4 (3) : 523-538. doi: 10.3934/dcdss.2011.4.523 [3] Antonio Suárez. A logistic equation with degenerate diffusion and Robin boundary conditions. Communications on Pure & Applied Analysis, 2008, 7 (5) : 1255-1267. doi: 10.3934/cpaa.2008.7.1255 [4] Adam Bobrowski, Katarzyna Morawska. From a PDE model to an ODE model of dynamics of synaptic depression. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2313-2327. doi: 10.3934/dcdsb.2012.17.2313 [5] Guowei Dai, Ruyun Ma, Haiyan Wang, Feng Wang, Kuai Xu. Partial differential equations with Robin boundary condition in online social networks. Discrete & Continuous Dynamical Systems - B, 2015, 20 (6) : 1609-1624. doi: 10.3934/dcdsb.2015.20.1609 [6] Fatiha Alabau-Boussouira. On the influence of the coupling on the dynamics of single-observed cascade systems of PDE's. Mathematical Control & Related Fields, 2015, 5 (1) : 1-30. doi: 10.3934/mcrf.2015.5.1 [7] Peter Brune, Björn Schmalfuss. Inertial manifolds for stochastic pde with dynamical boundary conditions. Communications on Pure & Applied Analysis, 2011, 10 (3) : 831-846. doi: 10.3934/cpaa.2011.10.831 [8] Satoshi Kosugi, Yoshihisa Morita, Shoji Yotsutani. A complete bifurcation diagram of the Ginzburg-Landau equation with periodic boundary conditions. Communications on Pure & Applied Analysis, 2005, 4 (3) : 665-682. doi: 10.3934/cpaa.2005.4.665 [9] Zhongkai Guo. Invariant foliations for stochastic partial differential equations with dynamic boundary conditions. Discrete & Continuous Dynamical Systems - A, 2015, 35 (11) : 5203-5219. doi: 10.3934/dcds.2015.35.5203 [10] Minoo Kamrani. Numerical solution of partial differential equations with stochastic Neumann boundary conditions. Discrete & Continuous Dynamical Systems - B, 2019, 24 (10) : 5337-5354. doi: 10.3934/dcdsb.2019061 [11] Kota Kumazaki, Adrian Muntean. Local weak solvability of a moving boundary problem describing swelling along a halfline. Networks & Heterogeneous Media, 2019, 14 (3) : 445-469. doi: 10.3934/nhm.2019018 [12] Raffaela Capitanelli. Robin boundary condition on scale irregular fractals. Communications on Pure & Applied Analysis, 2010, 9 (5) : 1221-1234. doi: 10.3934/cpaa.2010.9.1221 [13] Shin-Ichiro Ei, Toshio Ishimoto. Effect of boundary conditions on the dynamics of a pulse solution for reaction-diffusion systems. Networks & Heterogeneous Media, 2013, 8 (1) : 191-209. doi: 10.3934/nhm.2013.8.191 [14] Davide Guidetti. Partial reconstruction of the source term in a linear parabolic initial problem with Dirichlet boundary conditions. Discrete & Continuous Dynamical Systems - A, 2013, 33 (11&12) : 5107-5141. doi: 10.3934/dcds.2013.33.5107 [15] Rinaldo M. Colombo, Francesca Marcellini. Coupling conditions for the $3\times 3$ Euler system. Networks & Heterogeneous Media, 2010, 5 (4) : 675-690. doi: 10.3934/nhm.2010.5.675 [16] Gideon Simpson, Michael I. Weinstein, Philip Rosenau. On a Hamiltonian PDE arising in magma dynamics. Discrete & Continuous Dynamical Systems - B, 2008, 10 (4) : 903-924. doi: 10.3934/dcdsb.2008.10.903 [17] Anna-Lena Horlemann-Trautmann, Alessandro Neri. A complete classification of partial MDS (maximally recoverable) codes with one global parity. Advances in Mathematics of Communications, 2020, 14 (1) : 69-88. doi: 10.3934/amc.2020006 [18] Kingshook Biswas. Complete conjugacy invariants of nonlinearizable holomorphic dynamics. Discrete & Continuous Dynamical Systems - A, 2010, 26 (3) : 847-856. doi: 10.3934/dcds.2010.26.847 [19] Youshan Tao, J. Ignacio Tello. Nonlinear stability of a heterogeneous state in a PDE-ODE model for acid-mediated tumor invasion. Mathematical Biosciences & Engineering, 2016, 13 (1) : 193-207. doi: 10.3934/mbe.2016.13.193 [20] Ciro D'Apice, Peter I. Kogut, Rosanna Manzo. On relaxation of state constrained optimal control problem for a PDE-ODE model of supply chains. Networks & Heterogeneous Media, 2014, 9 (3) : 501-518. doi: 10.3934/nhm.2014.9.501

2018 Impact Factor: 1.143