Modeling the kinetics of acylation of insulin using a recursive method for solving the systems of coupled differential equations.

This paper describes a theoretical method for solving systems of coupled differential equations that describe the kinetics of complicated reaction networks in which a molecule having multiple reaction sites reacts irreversibly with multiple equivalents of a ligand (reagent). The members of the network differ in the number of equivalents of reagent that have reacted, and in the patterns of sites of reaction. A recursive algorithm generates series, asymptotic, and average solutions describing this kinetic scheme. This method was validated by successfully simulating the experimental data for the kinetics of acylation of insulin.