35.7NAMay 19
Structure-preserving local discontinuous Galerkin discretization of conformational conversion systemsPaola F. Antonietti, Mattia Corti, Sergio Gómez et al.
We investigate a two-state conformational conversion system and introduce a novel structure-preserving numerical scheme that couples a local discontinuous Galerkin space discretization with the backward Euler time-integration method. The model is first reformulated in terms of auxiliary variables involving suitable nonlinear transformations, which allow us to enforce positivity and boundedness at the numerical level. Then, we prove a discrete entropy-stability inequality, which we use to show the existence of discrete solutions, as well as to establish the convergence of the scheme by means of some discrete compactness arguments. As a by-product of the theoretical analysis, we also prove the existence of global weak solutions satisfying the system's physical bounds. Numerical results validate the theoretical results and assess the capabilities of the proposed method in practice.
24.6NAApr 20
High-fidelity and Network-based Spatio-temporal Mathematical Models of Alzheimer's Disease Progression and their Validation Against PET-SUVR Imaging DataBeatrice Caon, Mattia Corti, Francesca Bonizzoni et al.
Alzheimer's disease is the most common neurodegenerative disorder. Its pathological development is connected with the misfolding and accumulation of two toxic proteins: amyloid-beta and tau proteins. Mathematical models provide a valuable quantitative tool for monitoring disease progression. In this work, we proposed and compare a novel framework where the spatio-temporal dynamics of amyloid-beta and tau proteins is modeled based on employing either three-dimensional patient-specific geometries or through reduced network-based models defined on the brain connectome. More specifically, a high-fidelity biophysical model is proposed on three-dimensional brain geometries reconstructed from magnetic resonance imaging, whereas a network-based reduced formulation is defined on the brain connectome. For both approaches, a suitable numerical discretisation is proposed. A sensitivity analysis is presented to quantify the influence of model parameters on protein concentration patterns as well as compare the quality of the predictions. For both approaches, the results are validated against PET-SUVR clinical data using 18FAZD4694 for amyloid-beta and 18FMK6240 for tau protein. The results indicate that the three-dimensional model provides the most accurate and biologically consistent description of the disease progression, but remains computationally demanding. On the other hand, the reduced graph-based model is cheaper, but it is not always able to achieve reliable results.