Mathematical modeling of drug release from bioerodible microparticles: effect of gamma-irradiation.

Bioerodible polymers used in controlled drug delivery systems, such as poly(lactic-co-glycolic acid) (PLGA) undergo radiolytic degradation during gamma-irradiation. In spite of the considerable practical importance, yet only little knowledge is available on the consequences of this sterilization method on the resulting drug release patterns in a quantitative way. The major objectives of the present study were: (i) to monitor the effects of different gamma-irradiation doses on the physicochemical properties of drug-free and drug-loaded, PLGA-based microparticles; (ii) to analyze the obtained experimental results using adequate mathematical models; (iii) to get further insight into the occurring physical and chemical phenomena; and (iv) to relate the applied gamma-irradiation dose in a quantitative way to the resulting drug release rate. 5-Fluorouracil-loaded, PLGA-based microparticles were prepared with an oil-in-water solvent extraction method and exposed to gamma-irradiation doses ranging from 0 to 33 kGy. Size exclusion chromatography, differential scanning calorimetry, scanning electron microscopy, particle size analysis, determination of the actual drug loading and in vitro drug release kinetics were used to study the effects of the gamma-irradiation dose on the physicochemical properties of the microparticles. Two mathematical models-a simplified and a more comprehensive one-were used to analyze the experimental results. The simplified model considers drug diffusion based on Fick's second law for spherical geometry and a Higuchi-like pseudo-steady-state approach. The complex model combines Monte Carlo simulations (describing polymer erosion) with partial differential equations quantifying drug diffusion with time-, position- and direction-dependent diffusivities. Interestingly, exponential relationships between the gamma-irradiation dose and the initial drug diffusivity within the microparticles could be established. Based on this knowledge both models were used to predict the resulting drug release kinetics as a function of the gamma-irradiation dose. Importantly, the theoretical predictions were confirmed by experimental results.