Temperature-Sensitive Properties and Drug Release Processes of Chemically Cross-Linked Poly(N-isopropylacrylamide) Hydrogel: A Molecular Dynamics Simulation

This study utilized a dynamic cross-linking algorithm to formulate a chemical cross-linked hydrogel model of poly(N-isopropylacrylamide) (PNIPAM) with N, N’-methylenebisacrylamide (BIS). Molecular dynamics (MD) simulations were conducted to investigate the temperature sensitivity and ibuprofen release mechanism of this hydrogel under varying cross-linking degrees and water contents. The low critical solution temperature (LCST) of the hydrogel was determined based on changes in solvent-accessible surface area (SASA) and hydrogen bond count. The LCST was found to be between 300 and 310 K. As the temperature increased, both SASA and hydrogen bond counts generally exhibited a gradual decrease. However, near the LCST, polymer chain collapse temporarily exposed the hydrophilic groups of the PNIPAM, forming hydrophilic regions that increased the contact area with water. This led to a transient increase in SASA (8% higher than that before 300 K) and hydrogen bond counts (6.25% higher than that at 290 K). Concurrently, Young’s modulus of the PNIPAM hydrogel was found to decrease with increasing water content (from 3.11 GPa to 2.59 GPa, representing a 16.7% decrease when water content increased from 0% to 50% for 80% cross-linking degree) and increase with rising cross-linking density (from 2.02 GPa to 2.94 GPa, representing a 45.5% increase when the cross-linking degree increased from 0% to 80% for 20% water content). These findings indicate that enhancing cross-linking density is an effective strategy for improving the hydrogel’s mechanical properties. A PNIPAM–ibuprofen delivery model was constructed and molecular dynamics (MD) simulations were conducted, revealing temperature dependence release behavior. Below the LCST, the PNIPAM hydrogel remains in a highly swollen state (PNIPAM single-chain radius of gyration, Rg = 0.64 nm at 290 K), with ibuprofen molecules adsorbed within the PNIPAM polymer chain network. Conversely, above the LCST, PNIPAM undergoes phase separation (Rg decreases to 0.56 nm at 320 K, representing a 12.5% decrease), resulting in volume contraction (cavity volume reduced by 35%) and disruption of the hydrogen bond network. This process results in the release of ibuprofen molecules, accompanied by an increase in their diffusion coefficient from 1.3817 × 10−9 (280 K) to 4.2847 × 10−9 m2/s (320 K). Concurrently, the interaction energy with PNIPAM experiences a decline, from −126.72 kcal/mol to −108.69 kcal/mol. The findings of this study provide insights into the optimization of the structural stability of ibuprofen delivery carriers.