A thermodynamically consistent electro-chemo-mechanical theory for modeling axonal swelling

Abstract

In the present study, for the first time, a thermodynamically consistent large deformation theory is developed to model the multi physics problem of axonal swelling which is the hallmark of most of the brain diseases. To this end, the relevant axonal compartments are first explained and the corresponding model parts are introduced. Next, the problem is formulated as an open thermodynamic system and the corresponding constitutive and evolution equations are extracted utilizing the balance laws. Here, a multiplicative decomposition of the deformation gradient is used to capture the active behavior of the axonal actin cortex. Specific free energy functions are given for the model parts to complete the framework. While the developed model is general and can be extended to cover other types of axonal swellings, for the sake of brevity, two sepcific swelling modes are explored in the present study: swelling due to hypo-osmotic shocks and volume expansions because of actomyosin disruptions. To this end, three relevant sets of experimental data available in the literature are chosen verify the model: free equilibrium swelling of squid axoplasm, transient swelling of PC12 neurites due to hypo-osmotic shocks and steady state as well as transient swellings of embryonic Drosophila axons due to myosin motor blocking. The good comparison between the model predictions with realistically estimated parameters and the experimental data proves the validity of the model in predicting axonal swellings. The key findings of this study are as follows: first, it is shown that the kinetics of axonal swelling is surface limited in the presence of the axolemma barrier. Second, axonal swelling happens mainly in its radial direction due to its structural and geometrical constraints. Third, the multiphasic nature of the axoplasm can be accurately modeled as an ideal polyelectrolyte hydrogel. Finally the active contractile forces can be ignored for volume expansions due to hypo-osmotic shocks as the osmotic pressures are much bigger than the corresponding stresses due to active deformations.