Multi-order fractional differential equations and their numerical solution

Relaxation processes and reaction kinetics of proteins deviate from exponential behavior because of their large amount of conformational substrates. The dynamics are governed by many time scales and, therefore, the decay of the relaxation function or reactant concentration is slower than exponential. Applying the idea of self-similar dynamics, we derive a fractal scaling model that results in an equation in which the time derivative is replaced by a differentiation (d/dt)beta of non-integer order beta. The fractional order differential equation is solved by a Mittag-Leffler function. It depends on two parameters, a fundamental time scale tau 0 and a fractional order beta that can be interpreted as a self-similarity dimension of the dynamics. Application of the fractal model to ligand rebinding and pressure release measurements of myoglobin is demonstrated, and the connection of the model to considerations of energy barrier height distributions is shown.