In silico modeling of axonal reconnection within a discrete fiber tract after spinal cord injury

Following spinal cord injury (SCI), descending axons that carry motor commands from the brain to the spinal cord are injured or transected, producing chronic motor dysfunction and paralysis. Reconnection of these axons is a major prerequisite for restoration of function after SCI. Thus far, only modest gains in motor function have been achieved experimentally or in the clinic after SCI, identifying the practical limitations of current treatment approaches. In this paper, we use an ordinary differential equation (ODE) to simulate the relative and synergistic contributions of several experimentally-established biological factors related to inhibition or promotion of axonal repair and restoration of function after SCI. The factors were mathematically modeled by the ODE. The results of our simulation show that in a model system, many factors influenced the achievability of axonal reconnection. Certain factors more strongly affected axonal reconnection in isolation, and some factors interacted in a synergistic fashion to produce further improvements in axonal reconnection. Our data suggest that mathematical modeling may be useful in evaluating the complex interactions of discrete therapeutic factors not possible in experimental preparations, and highlight the benefit of a combinatorial therapeutic approach focused on promoting axonal sprouting, attraction of cut ends, and removal of growth inhibition for achieving axonal reconnection. Predictions of this simulation may be of utility in guiding future experiments aimed at restoring function after SCI.