A Hybrid Particle-Field Coarse-Grained Molecular Model for Pluronics Water Mixtures

Bacterial lipids drive compartmentalization on the nanoscale

The design of cellular functions in synthetic systems, inspired by the internal partitioning of living cells, is a constantly growing research field that is paving the way to a large number of new remarkable applications. Several hierarchies of internal compartments like polymersomes, liposomes, and membranes are used to control the transport, release, and chemistry of encapsulated species. However, the experimental characterization and the comprehension of glycolipid mesostructures are far from being fully addressed. Lipid A is indeed a glycolipid and the endotoxic part of Gram-negative bacterial lipopolysaccharide; it is the moiety that is recognized by the eukaryotic receptors giving rise to the modulation of innate immunity. Herein we propose, for the first time, a combined approach based on hybrid Particle-Field (hPF) Molecular Dynamics (MD) simulations and Small Angle X-Ray Scattering (SAXS) experiments to gain a molecular picture of the complex supramolecular structures of lipopolysaccharide (LPS) and lipid A at low hydration levels. The mutual support of data from simulations and experiments allowed the unprecedented discovery of the presence of a nano-compartmentalized phase composed of liposomes of variable size and shape which can be used in synthetic biological applications.

Molecular Insights on Morphology, Composition, and Stability of Mixed Micelles Formed by Ionic Surfactant and Nonionic Block Copolymer in Water Using Coarse-Grained Molecular Dynamics Simulations

The nanoscale association domains are the ultimate determinants of the macroscopic properties of complex fluids involving amphiphilic polymers and surfactants, and hence, it is foremost important to understand the role of polymer/surfactant concentration on these domains. We have used coarse-grained molecular dynamics simulations to investigate the effect of polymer/surfactant concentration on the morphology of poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO, i.e., pluronics or poloxamers) block copolymer, and ionic surfactants sodium dodecyl sulfate (SDS), mixed micelles in aqueous solution. The proclivity of the surfactant to form the mixed micelles is also probed using umbrella sampling simulations. In this study, we observed that the core of the pluronic + SDS formed mixed micelles consists of PPO, the alkyl tail of SDS, and some water molecules, whereas the PEO, water, and sulfate headgroups of SDS form a shell, consistent with experimental observations. The micelles are spherical at high-pluronic/low-SDS compositions, ellipsoidal at high-SDS/low-pluronic compositions, and wormlike-cylindrical at high-pluronic/high-SDS compositions. The transitions in micelle morphology are governed by the solvent accessible surface area of mixed aggregates, electrostatic repulsion between SDS-headgroups, and dehydration of PEO and PPO segments. The free energy barrier for SDS escape is much higher in mixed micelles than in pure SDS micelles, indicating a stronger tendency for SDS to form pluronic-SDS mixed micelles.

Insights into the Phase Diagram of Pluronic L64 Using Coarse-Grained Molecular Dynamics Simulations

We investigate the concentration-dependent phase diagram of pluronic L64 in aqueous media at 300 and 320 K using coarse-grained (CG) molecular dynamics (MD) simulations. The CG model is derived by adapting the Martini model for nonbonded interactions along with the Boltzmann inversion (BI) of bonded interactions from all-atom (AA) simulations. Our derived CG model successfully captures the experimentally observed micellar-, hexagonal-, lamellar-, and polymer-rich solution phase. The end-to-end distance reveals the conformational change from an open-chain structure in the micellar phase to a folded-chain structure in the lamellar phase, increasing the orientational order. An increase in temperature leads to expulsion of water molecules from the L64 moiety, suggesting an increase in L64 hydrophobicity. Thermodynamic analysis using the two-phase thermodynamics (2PT) method suggests the entropy of the system decreases with increasing L64 concentration and the decrease in free energy (F) with temperature is mainly driven by the entropic factor (-TS). Further, the increase in aggregation number at lower concentrations and self-assembly at very high concentrations is energetically driven, whereas the change from the micellar phase to the lamellar phase with increasing L64 concentration is entropically driven. Our model provides molecular insights into L64 phases which can be further explored to design functionality-based suprastructures for drug delivery and tissue engineering applications.

Computer modeling of polymer stars in variable solvent conditions: a comparison of MD simulations, self-consistent field (SCF) modeling and novel hybrid Monte Carlo SCF approach

Computer-aided modeling is a systematic approach to grasp the physics of macromolecules, but it remains essential to know when to trust the results and when not. For a polymer star, we consider three approaches: (i) Molecular Dynamics (MD) simulations and implementing a coarse-grained model, (ii) the self-consistent field approach based on a mean-field approximation and implementing the lattice model due to Scheutjens and Fleer (SF-SCF) and (iii) novel hybrid Monte Carlo self-consistent field (MC-SCF) method, which combines a coarse-grained model driven by a Monte Carlo method and a mean-field representation driven by SF-SCF. We compare the performance of these approaches under a wide range of solvent qualities. The MD approach is formally the most exact but suffers from reasonable convergence. The mean-field approach works similarly in all solvent qualities but is quantitatively least accurate. The MC-SCF hybrid allows us to combine the benefits of the simulation route and the effective performance of SCF. We consider the center-to-end distance R_ce, the radius of gyration R_g² of the star and the polymer density profiles φ(r) of polymer-segments in it. All three methods show a good qualitative agreement one to another. The MC-SCF method is in good agreement with the scaling predictions in the whole range of solvent quality values showing that it grasps the essential physics while remaining computationally in bounds.

Aggregation of Lipid A Variants: A Hybrid Particle-Field Model

Lipid A is one of the three components of bacterial lipopolysaccharides constituting the outer membrane of Gram-negative bacteria, and is recognized to have an important biological role in the inflammatory response of mammalians. Its biological activity is modulated by the number of acyl-chains that are present in the lipid and by the dielectric medium, i.e., the type of counter-ions, through electrostatic interactions. In this paper, we report on a coarse-grained model of chemical variants of Lipid A based on the hybrid particle-field/molecular dynamics approach (hPF-MD). In particular, we investigate the stability of Lipid A bilayers for two different hexa- and tetra-acylated structures. Comparing particle density profiles along bilayer cross-sections, we find good agreement between the hPF-MD model and reference all-atom simulation for both chemical variants of Lipid A. hPF-MD models of constituted bilayers composed by hexa-acylated Lipid A in water are stable within the simulation time. We further validate our model by verifying that the phase behavior of Lipid A/counterion/water mixtures is correctly reproduced. In particular, hPF-MD simulations predict the correct self-assembly of different lamellar and micellar phases from an initially random distribution of Lipid A molecules with counterions in water. Finally, it is possible to observe the spontaneous formation and stability of Lipid A vesicles by fusion of micellar aggregates.

Micellization of Pluronic P123 in Water/Ethanol/Turpentine Oil Mixed Solvents: Hybrid Particle-Field Molecular Dynamic Simulation

The hybrid particle-field molecular dynamics simulation method (MD-SCF) was applied to study the self-assembly of Pluronic PEO20-PPO70-PEO20 (P123) in water/ethanol/turpentine oil- mixed solvents. In particular, the micellization process of P123 at low concentration (less than 20%) in water/ethanol/turpentine oil-mixed solvents was investigated. The aggregation number, radius of gyration, and radial density profiles were calculated and compared with experimental data to characterize the structures of the micelles self-assembled from P123 in the mixed solvent. This study confirms that the larger-sized micelles are formed in the presence of ethanol, in addition to the turpentine oil-swollen micelles. Furthermore, the spherical micelles and vesicles were both observed in the self-assembly of P123 in the water/ethanol/turpentine oil-mixed solvent. The results of this work aid the understanding of the influence of ethanol and oil on P123 micellization, which will help with the design of effective copolymer-based formulations.

Biomembrane solubilization mechanism by Triton X-100: a computational study of the three stage model

The solubilization mechanism of lipid membranes in the presence of Triton X-100 (TX-100) is investigated at molecular resolution using hybrid particle field–self consistence field simulations.

Hybrid particle–field molecular dynamics simulation for polyelectrolyte systems

An effective hybrid computer simulation method combining molecular dynamics and self-consistent field theory is developed by including electrostatic interactions.

Mean-field coarse-grained model for poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) triblock copolymer systems

The microscopic modeling of surfactant systems is of the utmost importance in understanding the mechanisms related to the micellization process because it allows for prediction and comparison with experimental data of diverse equilibrium system properties. In this work, we present a coarse-grained model for Pluronics, a trademarked type of triblock copolymer, from simulations based on a single-chain mean-field theory (SCMF). This microscopic model is used to quantify the micellization process of these nonionic surfactants at 37 °C and has been shown to be able to quantitatively reproduce experimental data of the critical micelle concentration (CMC) along with other equilibrium properties. In particular, these results correctly capture the experimental behavior with respect to the lengths of the hydrophobic and hydrophilic moieties of the surfactants for low and medium hydrophobicities. However, for the more highly hydrophobic systems with low CMCs, a deviation is found which has been previously attributed to nonequilibrium effects in the experimental data (Garcı́a Daza, F. A.; Mackie, A. D. Low Critical Micelle Concentration Discrepancy between Theory and Experiment. J. Phys. Chem. Lett. 2014, 5, 2027-2032).

Generation of Well-Relaxed All-Atom Models of Large Molecular Weight Polymer Melts: A Hybrid Particle-Continuum Approach Based on Particle-Field Molecular Dynamics Simulations

A procedure based on Molecular Dynamics (MD) simulations employing soft potentials derived from self-consistent field (SCF) theory (named MD-SCF) able to generate well-relaxed all-atom structures of polymer melts is proposed. All-atom structures having structural correlations indistinguishable from ones obtained by long MD relaxations have been obtained for poly(methyl methacrylate) (PMMA) and poly(ethylene oxide) (PEO) melts. The proposed procedure leads to computational costs mainly related on system size rather than to the chain length. Several advantages of the proposed procedure over current coarse-graining/reverse mapping strategies are apparent. No parametrization is needed to generate relaxed structures of different polymers at different scales or resolutions. There is no need for special algorithms or back-mapping schemes to change the resolution of the models. This characteristic makes the procedure general and its extension to other polymer architectures straightforward. A similar procedure can be easily extended to the generation of all-atom structures of block copolymer melts and polymer nanocomposites.