Simulating growth of TDP-43 cytosolic inclusion bodies in neurons

Andrey V Kuznetsov
DOI: https://doi.org/10.1101/2023.11.28.569118
2024-08-28
Abstract:This paper introduces a mathematical model for the growth of transactive response DNA binding protein of 43 kDa (TDP-43) inclusion bodies in neuron soma. The model's equations enable numerical determination of the concentrations of TDP-43 dimers, monomers, and aggregates. Assuming all aggregates integrate into the inclusion bodies, the model predicts the size of TDP-43 inclusion bodies. An approximate solution to the model equations is derived for the scenario in which protein degradation machinery is dysfunctional, resulting in infinite half-lives for TDP-43 dimers, monomers, and aggregates. This solution, valid for large times, predicts that the inclusion body's radius increases proportionally to the cube root of time. To the best of the author's knowledge, this study presents the first attempt to model the relationship between the size of TDP-43 inclusion bodies and time. The sensitivity analysis of the approximate solution indicates that the concentrations of TDP-43 monomers and aggregates, as well as inclusion body radii, are independent of the kinetic constants. While the numerical solution of the full mathematical model continues to work with finite half-lives, the approximate solution becomes invalid for scenarios with physiologically relevant (finite) half-lives of TDP-43 dimers, monomers, and aggregates. In contrast to the situation with infinite half-lives, the numerical solution shows that for various values of kinetic constants, the curves representing concentrations of monomers and aggregates, as well as the curves depicting inclusion body radii, converge to distinct constant values.
Biophysics
What problem does this paper attempt to address?