Molecular dynamics simulation of micelle formation in amphiphilic solution

Interrogating the relationship between the microstructure of amphiphilic poly(ethylene glycol-b-caprolactone) copolymers and their colloidal assemblies using non-interfering techniques

Understanding the microstructural parameters of amphiphilic copolymers that control the formation and structure of aggregated colloids (e.g., micelles) is essential for the rational design of hierarchically structured systems for applications in nanomedicine, personal care and food formulations. Although many analytical techniques have been employed to study such systems, in this investigation we adopted an integrated approach using non-interfering techniques - diffusion nuclear magnetic resonance (NMR) spectroscopy, dynamic light scattering (DLS) and synchrotron small-angle X-ray scattering (SAXS) - to probe the relationship between the microstructure of poly(ethylene glycol-b-caprolactone) (PEG-b-PCL) copolymers [e.g., block molecular weight (MW) and the mass fraction of PCL (f_PCL)] and the structure of their aggregates. Systematic trends in the self-assembly behaviour were determined using a large family of well-defined block copolymers with variable PEG and PCL block lengths (number-average molecular weights (M_n) between 2 and 10 and 0.5-15 kDa, respectively) and narrow dispersity (Ð < 1.12). For all of the copolymers, a clear transition in the aggregate structure was observed when the hydrophobic f_PCL was increased at a constant PEG block M_n, although the nature of this transition is also dependent on the PEG block M_n. Copolymers with low M_n PEG blocks (2 kDa) were observed to transition from unimers and loosely associated unimers to metastable aggregates and finally, to cylindrical micelles as the f_PCL was increased. In comparison, copolymers with PEG block M_n of between 5 and 10 kDa transitioned from heterogenous metastable aggregates to cylindrical micelles and finally, well-defined ellipsoidal micelles (of decreasing aspect ratios) as the f_PCL was increased. In all cases, the diffusion NMR spectroscopy, DLS and synchrotron SAXS results provided complementary information and the grounds for a phase diagram relating copolymer microstructure to aggregation behaviour and structure. Importantly, the absence of commonly depicted spherical micelles has implications for applications where properties may be governed by shape, such as, cellular uptake of nanomedicine formulations.

Self-assembly of peptide amphiphiles by vapor pressure osmometry and dissipative particle dynamics

Peptide amphiphiles are one of the most promising materials in the biomedical field, so much effort has been devoted to characterizing the mechanism of their self-assembly and thermosensitive gelation. In this work, vapor pressure osmometry measurements were carried out to parameterize the thermosensitivity of interactions between peptide amphiphiles in an aqueous solution. The osmometry measurement verified that the peptides became more hydrophobic as temperature increased, which was quantitatively described with the Flory-Huggins χ parameter. Thereafter, a coarse-grained molecular model was used to simulate peptide amphiphiles dissolved in an aqueous solution. The temperature sensitive coarse-grained parameter a HW, which is the repulsive force between the hydrophilic head of the peptide amphiphile and water was estimated from the aforementioned experimentally obtained χ. Furthermore, the effects of concentration and temperature on the self-assembly behavior of peptide amphiphiles were quantitatively studied by dissipative particle dynamics. The simulation results revealed that a HW plays an important role in self-assembly characteristics and in the resulting microstructure of the peptide amphiphiles, which coincides with previous experimental and computational findings. The methodology in quantitatively linking the coarse-grained parameter from experiment and theory provides a sensible foundation for bridging future simulation studies with experimental work on macromolecules.

Congruency between biophysical data from multiple platforms and molecular dynamics simulation of the double-super helix model of nascent high-density lipoprotein

The predicted structure and molecular trajectories from >80 ns molecular dynamics simulation of the solvated Double-Super Helix (DSH) model of nascent high-density lipoprotein (HDL) were determined and compared with experimental data on reconstituted nascent HDL obtained from multiple biophysical platforms, including small angle neutron scattering (SANS) with contrast variation, hydrogen-deuterium exchange tandem mass spectrometry (H/D-MS/MS), nuclear magnetic resonance spectroscopy (NMR), cross-linking tandem mass spectrometry (MS/MS), fluorescence resonance energy transfer (FRET), electron spin resonance spectroscopy (ESR), and electron microscopy. In general, biophysical constraints experimentally derived from the multiple platforms agree with the same quantities evaluated using the simulation trajectory. Notably, key structural features postulated for the recent DSH model of nascent HDL are retained during the simulation, including (1) the superhelical conformation of the antiparallel apolipoprotein A1 (apoA1) chains, (2) the lipid micellar-pseudolamellar organization, and (3) the solvent-exposed Solar Flare loops, proposed sites of interaction with LCAT (lecithin cholesteryl acyltransferase). Analysis of salt bridge persistence during simulation provides insights into structural features of apoA1 that forms the backbone of the lipoprotein. The combination of molecular dynamics simulation and experimental data from a broad range of biophysical platforms serves as a powerful approach to studying large macromolecular assemblies such as lipoproteins. This application to nascent HDL validates the DSH model proposed earlier and suggests new structural details of nascent HDL.