Self-assembly mechanism of PEG-b-PCL and PEG-b-PBO-b-PCL amphiphilic copolymer micelles in aqueous solution from coarse grain modeling

We followed the self-assembly of high-molecular weight MePEG-b-PCL (poly(methyl ethylene glycol)-block-poly(epsilon-caprolactone)) diblock and MePEG-b-PBO-b-PCL (poly(methyl ethylene glycol)-block-poly(1,2-butylene oxide)-block-poly(epsilon-caprolactone)) into micelles using molecular dynamics simulation with a coarse grain (CG) force field based on quantum mechanics (CGq FF). The triblock polymer included a short poly(1,2-butylene oxide) (PBO) at the hydrophilic-hydrophobic interface of these systems. Keeping the hydrophilic length fixed (MePEG(45)), we considered 250 chains in which the hydrophobic length changed from PCL44 or PBO6-b-PCL43 to PCL62 or PBO9-b-PCL61. The polymers were solvated in explicit water for 2 mu s of simulations at 310.15 K. We found that the longer diblock system undergoes a morphological transition from an intermediate rod-like micelle to a prolate-sphere, while the micelle formed from the longer triblock system is a stable rod-like micelle. The two shorter diblock and triblock systems show similar self-assembly processes, both resulting in slightly prolate-spheres. The dynamics of the self-assembly is quantified in terms of chain radius of gyration, shape anisotropy, and hydration of the micelle cores. The final micelle structures are analyzed in terms of the local density components. We conclude that the CG model accurately describes the molecular mechanisms of self-assembly and the equilibrium micellar structures of hydrophilic and hydrophobic chains, including the quantity of solvent trapped inside the micellar core.