Wavefront invasion for a chemotaxis model of Multiple Sclerosis

In this work we study wavefront propagation for a chemotaxis reaction-diffusion system describing the demyelination in Multiple Sclerosis. Through a weakly non linear analysis, we obtain the Ginzburg–Landau equation governing the evolution of the amplitude of the pattern. We validate the analytical findings through numerical simulations. We show the existence of traveling wavefronts connecting two different steady solutions of the equations. The proposed model reproduces the progression of the disease as a wave: for values of the chemotactic parameter below threshold, the wave leaves behind a homogeneous plaque of apoptotic oligodendrocytes. For values of the chemotactic coefficient above threshold, the model reproduces the formation of propagating concentric rings of demyelinated zones, typical of Baló’s sclerosis.