Cartilage Dysfunction in ALS Patients as Side Effect of Motion Loss: 3D Mechano-Electrochemical Computational Model

Amyotrophic lateral sclerosis (ALS) is a debilitating motor neuron disease characterized by progressive weakness, muscle atrophy, and fasciculation. This fact results in a continuous degeneration and dysfunction of articular soft tissues. Specifically, cartilage is an avascular and nonneural connective tissue that allows smooth motion in diarthrodial joints. Due to the avascular nature of cartilage tissue, cells nutrition and by-product exchange are intermittently occurring during joint motions. Reduced mobility results in a change of proteoglycan density, osmotic pressure, and permeability of the tissue. This work aims to demonstrate the abnormal cartilage deformation in progressive immobilized articular cartilage for ALS patients. For this aim a novel 3D mechano-electrochemical model based on the triphasic theory for charged hydrated soft tissues is developed. ALS patient parameters such as tissue porosity, osmotic coefficient, and fixed anions were incorporated. Considering different mobility reduction of each phase of the disease, results predicted the degree of tissue degeneration and the reduction of its capacity for deformation. The present model can be a useful tool to predict the evolution of joints in ALS patients and the necessity of including specific cartilage protectors, drugs, or maintenance physical activities as part of the symptomatic treatment in amyotrophic lateral sclerosis

A Mathematical Modelling Study of Chemotactic Dynamics in Cell Cultures: The Impact of Spatio-temporal Heterogeneity

As motivated by studies of cellular motility driven by spatiotemporal chemotactic gradients in microdevices, we develop a framework for constructing approximate analytical solutions for the location, speed and cellular densities for cell chemotaxis waves in heterogeneous fields of chemoattractant from the underlying partial differential equation models. In particular, such chemotactic waves are not in general translationally invariant travelling waves, but possess a spatial variation that evolves in time, and may even oscillate back and forth in time, according to the details of the chemotactic gradients. The analytical framework exploits the observation that unbiased cellular diffusive flux is typically small compared to chemotactic fluxes and is first developed and validated for a range of exemplar scenarios. The framework is subsequently applied to more complex models considering the chemoattractant dynamics under more general settings, potentially including those of relevance for representing pathophysiology scenarios in microdevice studies. In particular, even though solutions cannot be constructed in all cases, a wide variety of scenarios can be considered analytically, firstly providing global insight into the important mechanisms and features of cell motility in complex spatiotemporal fields of chemoattractant. Such analytical solutions also provide a means of rapid evaluation of model predictions, with the prospect of application in computationally demanding investigations relating theoretical models and experimental observation, such as Bayesian parameter estimation

Cell dynamics in microfluidic devices under heterogeneous chemotaxis and growth conditions: a mathematical study

As motivated by studies of cellular motility driven by spatiotemporal chemotactic gradients in microdevices, we develop a framework for constructing approximate analytical solutions for the location, speed and cellular densities for cell chemotaxis waves in heterogeneous fields of chemoattractant from the underlying partial differential equation models. In particular, such chemotactic waves are not in general translationally invariant travelling waves, but possess a spatial variation that evolves in time, and may even may oscillate back and forth in time, according to the details of the chemotactic gradients. The analytical framework exploits the observation that unbiased cellular diffusive flux is typically small compared to chemotactic fluxes and is first developed and validated for a range of exemplar scenarios. The framework is subsequently applied to more complex models considering the full dynamics of the chemoattractant and how this may be driven and controlled within a microdevice by considering a range of boundary conditions. In particular, even though solutions cannot be constructed in all cases, a wide variety of scenarios can be considered analytically, firstly providing global insight into the important mechanisms and features of cell motility in complex spatiotemporal fields of chemoattractant. Such analytical solutions also provide a means of rapid evaluation of model predictions, with the prospect of application in computationally demanding investigations relating theoretical models and experimental observation, such as Bayesian parameter estimation.Comment: 32 pages, 11 figure

A new reliability-based data-driven approach for noisy experimental data with physical constraints

Data Science has burst into simulation-based engineering sciences with an impressive impulse. However, data are never uncertainty-free and a suitable approach is needed to face data measurement errors and their intrinsic randomness in problems with well-established physical constraints. As in previous works, this problem is here faced by hybridizing a standard mathematical modeling approach with a new data-driven solver accounting for the phenomenological part of the problem, with the aim of finding a solution point, satisfying some constraints, that minimizes a distance to a given data-set. However, unlike such works that are established in a deterministic framework, we use the Mahalanobis distance in order to incorporate statistical second order uncertainty of data in computations, i.e. variance and correlation. We develop the underlying stochastic theoretical framework and establish the fundamental mathematical and statistical properties. The performance of the resulting reliability-based data-driven procedure is evaluated in a simple but illustrative unidimensional problem as well as in a more realistic solution of a 3D structural problem with a material with intrinsically random constitutive behavior as concrete. The results show, in comparison with other data-driven solvers, better convergence, higher accuracy, clearer interpretation, and major flexibility besides the relevance of allowing uncertainty management with low computational demand

Non-adiabatic small polaron hopping in the n=3 Ruddlesden-Popper compound Ca4Mn3O10

Magnetotransport properties of the compound Ca4Mn3O10 are interpreted in terms of activated hopping of small magnetic polarons in the non-adiabatic regime. Polarons are most likely formed around Mn3+ sites created by oxygen substoichiometry. The application of an external field reduces the size of the magnetic contribution to the hopping barrier and thus produces an increase in the conductivity .We argue that the change in the effective activation energy around TN is due to the crossover to VRH conduction as antiferromagnetic order sets in.Comment: 29 pages, 7 figure

Understanding glioblastoma invasion using physically-guided neural networks with internal variables.

Microfluidic capacities for both recreating and monitoring cell cultures have opened the door to the use of Data Science and Machine Learning tools for understanding and simulating tumor evolution under controlled conditions. In this work, we show how these techniques could be applied to study Glioblastoma, the deadliest and most frequent primary brain tumor. In particular, we study Glioblastoma invasion using the recent concept of Physically-Guided Neural Networks with Internal Variables (PGNNIV), able to combine data obtained from microfluidic devices and some physical knowledge governing the tumor evolution. The physics is introduced in the network structure by means of a nonlinear advection-diffusion-reaction partial differential equation that models the Glioblastoma evolution. On the other hand, multilayer perceptrons combined with a nodal deconvolution technique are used for learning the go or grow metabolic behavior which characterises the Glioblastoma invasion. The PGNNIV is here trained using synthetic data obtained from in silico tests created under different oxygenation conditions, using a previously validated model. The unravelling capacity of PGNNIV enables discovering complex metabolic processes in a non-parametric way, thus giving explanatory capacity to the networks, and, as a consequence, surpassing the predictive power of any parametric approach and for any kind of stimulus. Besides, the possibility of working, for a particular tumor, with different boundary and initial conditions, permits the use of PGNNIV for defining virtual therapies and for drug design, thus making the first steps towards in silico personalised medicine

Data-Driven Computational Simulation in Bone Mechanics.

The data-driven approach was formally introduced in the field of computational mechanics just a few years ago, but it has gained increasing interest and application as disruptive technology in many other fields of physics and engineering. Although the fundamental bases of the method have been already settled, there are still many challenges to solve, which are often inherently linked to the problem at hand. In this paper, the data-driven methodology is applied to a particular problem in tissue biomechanics, a context where this approach is particularly suitable due to the difficulty in establishing accurate and general constitutive models, due to the intrinsic intra and inter-individual variability of the microstructure and associated mechanical properties of biological tissues. The problem addressed here corresponds to the characterization and mechanical simulation of a piece of cortical bone tissue. Cortical horse bone tissue was mechanically tested using a biaxial machine. The displacement field was obtained by means of digital image correlation and then transformed into strains by approximating the displacement derivatives in the bone virtual geometric image. These results, together with the approximated stress state, assumed as uniform in the small pieces tested, were used as input in the flowchart of the data-driven methodology to solve several numerical examples, which were compared with the corresponding classical model-based fitted solution. From these results, we conclude that the data-driven methodology is a useful tool to directly simulate problems of biomechanical interest without the imposition (model-free) of complex spatial and individually-varying constitutive laws. The presented data-driven approach recovers the natural spatial variation of the solution, resulting from the complex structure of bone tissue, i.e. heterogeneity, microstructural hierarchy and multifactorial architecture, making it possible to add the intrinsic stochasticity of biological tissues into the data set and into the numerical approach