In silico modelling of drug–polymer interactions for pharmaceutical formulations

Selecting polymers for drug encapsulation in pharmaceutical formulations is usually made after extensive trial and error experiments. To speed up excipient choice procedures, we have explored coarse-grained computer simulations (dissipative particle dynamics (DPD) and coarse-grained molecular dynamics using the MARTINI force field) of polymer–drug interactions to study the encapsulation of prednisolone (logp = 1.6), paracetamol (logp = 0.3) and isoniazid (logp = −1.1) in poly(l -lactic acid) (PLA) controlled release microspheres, as well as the encapsulation of propofol (logp = 4.1) in bioavailability enhancing quaternary ammonium palmitoyl glycol chitosan (GCPQ) micelles. Simulations have been compared with experimental data. DPD simulations, in good correlation with experimental data, correctly revealed that hydrophobic drugs (prednisolone and paracetamol) could be encapsulated within PLA microspheres and predicted the experimentally observed paracetamol encapsulation levels (5–8% of the initial drug level) in 50 mg ml−1 PLA microspheres, but only when initial paracetamol levels exceeded 5 mg ml−1 . However, the mesoscale technique was unable to model the hydrophilic drug (isoniazid) encapsulation (4–9% of the initial drug level) which was observed in experiments. Molecular dynamics simulations using the MARTINI force field indicated that the self-assembly of GCPQ is rapid, with propofol residing at the interface between micellar hydrophobic and hydrophilic groups, and that there is a heterogeneous distribution of propofol within the GCPQ micelle population. GCPQ–propofol experiments also revealed a population of relatively empty and drug-filled GCPQ particles.