A theoretical study of the distribution of insulin-like growth factor in human articular cartilage

We present a mathematical simulation which integrates the mechanisms that are currently believed to govern the concentration of the growth factor, IGF1, in cartilage. Articular cartilage is treated as a two-layer continuum: a thin surface layer, exposed to synovial fluid, with a higher cell density, and a deeper layer with impermeable bony endplate. A system of differential equations accounts for diffusion of IGF1 from synovial fluid into, and throughout, the cartilage; IGF1 synthesis, its reactions with soluble binding protein, with cell receptors, and with immobile binding sites on the extracellular matrix. We have collected all available physiologic data relevant to the solution of these equations and used it to compute numerical solutions that yield time dependent profiles for free and complex IGF1 throughout the depth of normal cartilage. Equations for osteoarthritic cartilage were formulated as well. Numerical results indicate a time-scale of several days for IGF1 profiles to settle down after a disturbance. The number of cell receptors for IGF1 appears to be more important than their rate of internalization. There is a lower bound to the number of cell receptors and of immobile binding sites. Parameters that await experimental determination are identified.