Molecular Dynamics Simulation of Amorphous Hydroxypropylmethylcellulose and Its Mixtures With Felodipine and Water

Understanding drug-polymer molecular interactions, their miscibility, supersaturation potential, and the effects of water uptake may be invaluable for selecting amorphous polymer dispersions that can maximize the oral bioavailability of poorly water-soluble drugs. Molecular dynamics simulations were performed using a model for hydroxypropylmethylcellulose (HPMC) resembling the substitution patterns found experimentally. HPMC at low and high water contents (0.9%-23.0% wt/wt) and mixtures with a hydrophobic drug, felodipine (FEL), were constructed. T_g values and densities after ∼30 ns aging at 298 K were close to published results. Except for hydrogen bonds (HBs) between the 5-O- and a 3-OH group in a neighboring repeat unit, HPMC oxygen atoms have a low HB probability (p < 0.1) perhaps due to shielding by surrounding substituents. Water molecules tend to be isolated at low water content while clusters were prevalent at ≥10.7% water. The Flory-Huggins FEL-HPMC interaction parameter (-0.20 ± 0.07) predicts complete miscibility at all HPMC compositions, in agreement with experiments. However, HBs between the FEL-N-H and HPMC favoring miscibility are disrupted with increasing water. Apparent diffusion coefficients versus water content were generated for water and FEL and a theory for the non-Einsteinian nature of water diffusion is proposed.