Coarse-Grained Simulations of Release of Drugs Housed in Flexible Nanogels: New Insights into Kinetic Parameters

Zwitterion-Modified Nanogel Responding to Temperature and Ionic Strength: A Dissipative Particle Dynamics Simulation

The self-assembly and stimuli-responsive properties of nanogel poly(n-isopropylacrylamide) (p(NIPAm)) and zwitterion-modified nanogel poly(n-isopropylacrylamide-co-sulfobetainemethacrylate) (p(NIPAm-co-SBMA)) were explored by dissipative particle dynamics simulations. Simulation results reveal that for both types of nanogel, it is beneficial to form spherical nanogels at polymer concentrations of 5-10%. When the chain length (L) elongates from 10 to 40, the sizes of the nanogels enlarge. As for the p(NIPAm) nanogel, it shows thermoresponsiveness; when it switches to the hydrophilic state, the nanogel swells, and vice versa. The zwitterion-modified nanogel p(NIPAm-co-SBMA) possesses thermoresponsiveness and ionic strength responsiveness concurrently. At 293 K, both hydrophilic p(NIPAm) and superhydrophilic polysulfobetaine methacrylate (pSBMA) could appear on the outer surface of the nanogel; however, at 318 K, superhydrophilic pSBMA is on the outer surface to cover the hydrophobic p(NIPAm) core. As the temperature rises, the nanogel shrinks and remains antifouling all through. The salt-responsive property can be reflected by the nanogel size; the volumes of the nanogels in saline systems are larger than those in salt-free systems as the ionic condition inhibits the shrinkage of the zwitterionic pSBMA. This work exhibits the temperature-responsive and salt-responsive behavior of zwitterion-modified-pNIPAm nanogels at the molecular level and provides guidance in antifouling nanogel design.

Interplay between Diffusion and Bond Cleavage Reaction for Determining Release in Polymer-Drug Conjugates

In conjugated polymeric drug delivery systems, both the covalent bond degradation rate and the diffusion of the freely moving drug particles affect the release profile of the formulation. Using Monte Carlo simulations in spherical matrices, the release kinetics resulting from the competition between the reaction and diffusion processes is discussed. For different values of the relative bond cleavage rate, varied over four orders of magnitude, the evolution of (i) the number of bonded drug molecules, (ii) the fraction of the freely moved detached drug within the polymer matrix, and (iii) the resulting fractional release of the drug is presented. The characteristic release time scale is found to increase by several orders of magnitude as the cleavage reaction rate constant decreases. The two extreme rate-limiting cases where either the diffusion or the reaction dominates the release are clearly distinguishable. The crossover between the diffusion-controlled and reaction-controlled regimes is also examined and a simple analytical formula is presented that can describe the full dependence of the release time on the bond cleavage rate constant. This simple relation is provided simply by the sum of the characteristic time for purely diffusional release and the bond cleavage decay time, which equals the inverse of the reaction rate constant.