Network development in biological gels: Role in lymphatic vessel development

In this paper, we present a model that explains the prepatterning of lymphatic vessel morphology in collagen gels. This model is derived using the theory of two phase rubber material due to Flory and coworkers and it consists of two coupled fourth order partial differential equations describing the evolution of the collagen volume fraction, and the evolution of the proton concentration in a collagen implant; as described in experiments of Boardman and Swartz (Circ. Res. 92, 801-808, 2003). Using linear stability analysis, we find that above a critical level of proton concentration, spatial patterns form due to small perturbations in the initially uniform steady state. Using a long wavelength reduction, we can reduce the two coupled partial differential equations to one fourth order equation that is very similar to the Cahn-Hilliard equation; however, it has more complex nonlinearities and degeneracies. We present the results of numerical simulations and discuss the biological implications of our model.