Numerical Simulation of Swelling and Drug Release Processes for Weak Polyelectrolyte Hydrogels†

doi:10.7503/cjcu20170541

Abstract

Abstract:

Based on the three-phase model, a mathematical model of swelling and drug release for weak polyelectrolyte hydrogels was established by combination of Langmuir isothermal adsorption equation and drug release equation; the finite element method was adopted to investigate the influence of the material parameters on swelling and drug release in the case of drugs were distributed within the weak polyelectrolyte hydrogels homogeneously; the mechanism of drugs released from weak polyelectrolyte hydrogels in different moments was analyzed as well; then the impacts of drug distribution conditions within the hydrogels which consist of homogeneous distribution and inhomogeneous distribution on swelling and drug release were compared comprehensively, and it has been found that drug release process is consistent with zero order release kinetics when drugs presented the sine distribution in the hydrogels; finally, the simulated results were comparatively compared with Beebe’s experimental results to validate the accuracy of this simulation.

Key words: Weak polyelectrolyte hydrogel, Swelling, Drug release, Mathematical model, Finite element simulation

TrendMD:

YE Hui, LIU Yabo, JIA Yuxi. Numerical Simulation of Swelling and Drug Release Processes for Weak Polyelectrolyte Hydrogels^†[J]. Chem. J. Chinese Universities, 2018, 39(4): 817.

Figures/Tables 10

Table 1 Material parameters and values

Reference object	Parameter	Value
Free ions	Diffusion coefficient, D_i/(m²·s^-1)	1×1 $0 - 9 27$
Drugs	Initial concentration, c_d0/(kg·m^-3)	1.264
Solution	Viscosity, η/(Pa·s)	1×10^-3
	Fluid pressure, p₀	0
	Electric potential, ψ₀	0
Polymer network	Initial volume fraction of water, ϕ_w0	0.8
	Initial concentration of fixed charge, c_f0/(mmol·L^-1)	2.12

Table 1 Material parameters and values

Reference object	Parameter	Value
Free ions	Diffusion coefficient, D_i/(m²·s^-1)	1×1 $0 - 9 27$
Drugs	Initial concentration, c_d0/(kg·m^-3)	1.264
Solution	Viscosity, η/(Pa·s)	1×10^-3
	Fluid pressure, p₀	0
	Electric potential, ψ₀	0
Polymer network	Initial volume fraction of water, ϕ_w0	0.8
	Initial concentration of fixed charge, c_f0/(mmol·L^-1)	2.12

Fig.1 Change of swelling degree with time for different ƒsp ƒsp: a. 1×1013; b. 3×1013; c. 5×1013; d. 7×1013.

Fig.2 Change of drug release amount with time for different ƒspƒsp: a. 1×1013; b. 3×1013; c. 5×1013; d. 7×1013.

Fig.3 Change of swelling degree with time for different polymer density Polymer density/(kg·m-3): a. 1200; b. 1500;c. 1600.

Fig.4 Change of drug release amount with time for different polymer densityPolymer density/(kg·m-3): a. 1200; b. 1500;c. 1600.

Table 2 Exponential relation of drug release ratio to time

Diffusive type	n	Control step of drug release
Case Ⅰ diffusion(Fickian diffusion)	0.45	Drug diffusion
Case Ⅱ diffusion(Relaxation diffusion)	0.89	Polymer relaxation
Case Ⅲ diffusion(Anomalous diffusion)	0.45-0.89	Drug diffusion and polymer relaxation

Fig.5 Logarithmic diagram of drug release ratio to time for different data fitting models

Fig.6 Radial distribution of drug concentration for c0=3.5sin(r×π/4)(A) and c0=3.5[1-sin(r×π/4)](B)

Fig.7 Change of drug release amount with time for 3.5(a), 3.5[1-sin(r×π/4)](b) and 3.5sin(r×π/4)(c)

Fig.8 Change of swelling degree with pH valuesa. Simulated data; b. experimental data.

References 40

[1]	Longo G. S., Cruz M. O. D. L., Szleifer I., Macromolecules,2011, 44(1), 147-158
[2]	Adimi M., Attar H., Barati A., Seifkordi A., Bakhtiari M., Int. J. Biosci., 2014, 5(11), 129-137
[3]	Wilde L., Bock M., Glöckl G., Garbacz G., Weitschieset W., Int. J. Pharm., 2014, 461(1), 296-300
[4]	Longo G. S., Cruz M. O. D. L., Szleifer I., Soft Matter, 2016, 12(40), 8359-8366
[5]	Helminger M., Wu B., Kollmann T., Benke D., Schwahn D., Pipich V., Faivre D., Zahn D., Cölfen H., Adv. Funct. Mater., 2014, 24(21), 3187-3196
[6]	Chen S. C., Wu Y. C., Mi F. L., Lin Y. H., Yu L. C., Sung H. W., J. Control. Release, 2004, 96(2), 285-300
[7]	Singh N. K., Lee D. S., J. Control. Release, 2014, 193, 214-227
[8]	Cevik O., Gidon D., Kizilel S., Acta Biomater., 2015, 11, 151-161
[9]	Beebe D. J., Moore J. S., Bauer J. M., Yu Q., Liu R. H., Devadoss C., Jo B. H., Nature,2000, 404(6778), 588-590
[10]	Kurnia J. C., Birgersson E., Mujumdar A. S., Quah L. C., Adv. Mater. Res., 2009, 74, 33-36
[11]	Che Y., Zinchenko A., Murata S., J. Colloid Interf. Sci., 2015, 445, 364-370
[12]	Deligkaris K., Tadele T. S., Olthuis W., Berg A. V. D., Sensor. Actuat. B: Chem., 2010, 147(2), 765-774
[13]	Marshall A. J., Blyth J., Davidson C. A. B., Lowe C. R., Chem A., Anal. Chem., 2003, 75(17), 4423-4431
[14]	Kwon H. J., Yasuda K., Gong J. P., Ohmiya Y., Macromol. Res., 2014, 22(3), 227-235
[15]	Gupta P., Vermani K., Garg S., Drug Discov. Today, 2002, 7(10), 569-579
[16]	Zhang Z., Chen L., Zhao C. W., Jing X. B., Polymer,2011, 52(3), 676-682
[17]	Tanaka T., Hocker L. O., Benedek G. B., J. Chem. Phys., 1973, 59(9), 5151-5159
[18]	Tanaka T., Fillmore D. J., J. Chem. Phys., 1979, 70(3), 1214-1218
[19]	Li Y., Tanaka T., J. Chem. Phys., 1990, 92(2), 1365-1371
[20]	Wang C. J., Li Y., Hu Z. B., Macromolecules,1997, 30(16), 4727-4732
[21]	Wallmersperger T., Kröplin B., Gülch R. W., Mech. Mater., 2004, 36(5), 411-420
[22]	Attaran A., Brummund J., Wallmersperger T., Smart Mater. Struct., 2015, 24(3), 035021
[23]	Li H., Ng T. Y., Yew Y. K., Lam K. Y., Biomacromolecules,2005, 6(1), 109-120
[24]	Li H., Lai F., Luo R., Langmuir,2009, 25(22), 13142-13150
[25]	Feng L. G., Jia Y. X., Li X., An L. J., J. Mech. Behav. Biomed. Mater., 2011, 4(7), 1328-1335
[26]	Feng L. G., Jia Y. X., Chen X. L., Li X., An L. J., J. Chem. Phys., 2010, 133(11), 114904
[27]	Lee P. I., J. Control. Release, 1985, 2, 277-288
[28]	Xu Y. H., Jia Y. X., Wang Z., Wang Z. J., J. Pharm. Sci., 2013, 102(5), 1532-1543
[29]	Xu Y. H., Zhang Y. Y., Jia Y. X., An L. J., J. Control. Release, 2013, 1(172), e88
[30]	Liu Y. B., Xu Y. H., Zhao Y. R., Jia Y. X., Chem. Res. Chinese Universities, 2016, 32(2), 302-310
[31]	Liu Y. B., Xu Y. H., Zhu P., Jia Y. X., J. Control. Release, 2015, 213, e144
[32]	De S. K., Aluru N. R., Johnson B., Crone W. C., J. Microelectromech. S., 2002, 11(5), 544-555
[33]	Redlich O., Peterson D. L., J. Phys. Chem., 1959, 63(6), 1024-1024
[34]	Gillespie D., Nonner W., Eisenberg R. S., J. Phys. Condens. Mat., 2002, 14(46), 12129
[35]	Ninawe P. R., Parulekar S. J., Ind. Eng. Chem. Res., 2012, 51, 1741
[36]	Korsmeyer R. W., Meerwall E. V., Peppas N. A., J. Polym. Sci. B: Polym. Phys., 1986, 24(2), 409-434
[37]	Ye H., Liu Y. B., Gao L. L., Du T. X., Jia Y. X.,Polym. Int., 2017, doi:10.1002/pi.5429
[38]	Nanao H., Liu X., Mori S., Wear,2002, 257(7), 665-670
[39]	Argon A. S., Cohen R. E., Patel A. C., Comput. Theor. Polym. Sci., 1999, 9(3/4), 339-352
[40]	Namazi H., Kanani A., Carbohyd. Polym., 2009, 76(1), 46-50