Modeling the kinetics of acylation of insulin using a recursive method forsolving the systems of coupled differential equations

This paper describes a theoretical method for solving systems of coupled di fferential equations that describe the kinetics of complicated reaction net works in which a molecule having multiple reaction sites reacts irreversibl y with multiple equivalents of a ligand (reagent). The members of the netwo rk 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 met hod was validated by successfully simulating the experimental data for the kinetics of acylation of insulin.