Phosphatidylcholine reverse micelles on the wrong track in molecular dynamics simulations of phospholipids in an organic solvent

Ion-Specific Effects on Ion and Polyelectrolyte Solvation

Ion-specific effects on aqueous solvation of monovalent counter ions, Na+ ${^+ }$ , K+ ${^+ }$ , Cl- ${^- }$ , and Br- ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- ${^- }$ and Br- ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

Dioctyl sodium sulfosuccinate surfactant self-assembly dependency of solvent hydrophilicity: a modelling study

The self-assembly of dioctyl sodium sulfosuccinate (AOT) model surfactant in solvent environments of differing polarity is examined by means of dissipative particle dynamics (DPD) bead model parametrized against Hildebrand solubility parameters from atomistic molecular dynamics (MD) simulations. The model predicts that in hydrophobic solvents (e.g. dodecane) the surfactant forms small (Nagg ∼ 8) reverse micellar aggregates, while in a solvent corresponding to water lamellar assembly takes place, in good agreement with literature structural parameters. Interestingly, solvents of intermediate polarity lead to formation of large, internally structured aggregates. In these, the surfactant headgroups cluster within the aggregate, surrounded by a continuous phase formed by the hydrocarbon tails. We show that the partitioning of the headgroups between the aggregate surface layer and the inner clustered phase depends primarily on solvent polarity, and can be controlled by the solvent, but also system composition. Finally, we compare the DPD assembly response to simplified effective interaction potentials derived at dilute concentration limit for the interactions. The comparison reveals that the simplified effective potential descriptions provide good level of insight on the assembly morphologies, despite drastic, isotropic interactions simplification involved.

Spray Dried Rugose Lipid Particle Platform for Respiratory Drug Delivery

PURPOSE

To develop a new lipid-based particle formulation platform for respiratory drug delivery applications. To find processing conditions for high surface rugosity and manufacturability. To assess the applicability of the new formulation method to different lipids.

METHODS

A new spray drying method with a simplified aqueous suspension feedstock preparation process was developed for the manufacture of rugose lipid particles of 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC). A study covering a wide range of feedstock temperatures and outlet temperatures was conducted to optimize the processing conditions. Aerosol performance was characterized in vitro and in silico to assess the feasibility of their use in respiratory drug delivery applications. The applicability of the new spray drying method to longer-chain phospholipids with adjusted spray drying temperatures was also evaluated.

RESULTS

Highly rugose DSPC lipid particles were produced via spray drying with good manufacturability. A feedstock temperature close to, and an outlet temperature lower than, the main phase transition were identified as critical in producing particles with highly rugose surface features. High emitted dose and total lung dose showed promising aerosol performance of the produced particles for use as a drug loading platform for respiratory drug delivery. Two types of longer-chain lipid particles with higher main phase transition temperatures, 1,2-diarachidoyl-sn-glycero-3-phosphocholine (DAPC) and 1,2-dibehenoyl-sn-glycero-3-phosphocholine (22:0 PC), yielded similar rugose morphologies when spray dried at correspondingly higher processing temperatures.

CONCLUSIONS

Rugose lipid particles produced via spray drying from an aqueous suspension feedstock are promising as a formulation platform for respiratory drug delivery applications. The new technique can potentially produce rugose particles using various other lipids.

Physisorption of bio oil nitrogen compounds onto montmorillonite

We assess computationally the adsorption of a series of nitrogen containing heterocycles and fatty acid amides from bio-oil on a model clay surface, Na-montmorillonite. The adsorption energies and conformations predicted by atomistic detail molecular dynamics (MD) simulations are compared against density functional theory (DFT) based molecular electrostatic potentials (MEP) and Hirshfeld, AIM, Merz-Singh-Kollman, and ChelpG charges. MD predicts systematically adsorption via cation bridging with adsorption strength of the heterocycles following purine > pyridine > imidazole > pyrrole > indole > quinoline. The fatty acid amides adsorption strength follows the steric availability and bulkiness of the head group. A comparison against the DFT calculations shows that MEP predicts adsorption geometries and the MD simulations reproduce the conformations for single adsorption site species. However, the DFT derived charge distibutions show that MD force-fields with non-polarizable fixed partial charge representations parametrized for aqueous environments cannot be used in apolar solvent environments without careful accuracy considerations. The overall trends in adsorption energies are reproduced by the Charmm GenFF employed in the MD simulations but the adsorption energies are systematically overestimated in this apolar solvent environment. The work has significance both for revealing nitrogen compound adsorption trends in technologically relevant bio oil environments but also as a methodological assessment revealing the limits of state of the art biomolecular force-fields and simulation protocols in apolar bioenvironments.

pH-Induced Changes in Polypeptide Conformation: Force-Field Comparison with Experimental Validation

Microsecond-long all-atom molecular dynamics (MD) simulations, circular dichroism, laser Doppler velocimetry, and dynamic light-scattering techniques have been used to investigate pH-induced changes in the secondary structure, charge, and conformation of poly l-lysine (PLL) and poly l-glutamic acid (PGA). The employed combination of the experimental methods reveals for both PLL and PGA a narrow pH range at which they are charged enough to form stable colloidal suspensions, maintaining their α-helix content above 60%; an elevated charge state of the peptides required for colloidal stability promotes the peptide solvation as a random coil. To obtain a more microscopic view on the conformations and to verify the modeling performance, peptide secondary structure and conformations rising in MD simulations are also examined using three different force fields, i.e., OPLS-AA, CHARMM27, and AMBER99SB*-ILDNP. Ramachandran plots reveal that in the examined setup the α-helix content is systematically overestimated in CHARMM27, while OPLS-AA overestimates the β-sheet fraction at lower ionization degrees. At high ionization degrees, the OPLS-AA force-field-predicted secondary structure fractions match the experimentally measured distribution most closely. However, the pH-induced changes in PLL and PGA secondary structure are reasonably captured only by the AMBER99SB*-ILDNP force field, with the exception of the fully charged PGA in which the α-helix content is overestimated. The comparison to simulations results shows that the examined force fields involve significant deviations in their predictions for charged homopolypeptides. The detailed mapping of secondary structure dependency on pH for the polypeptides, especially finding the stable colloidal α-helical regime for both examined peptides, has significant potential for practical applications of the charged homopolypeptides. The findings raise attention especially to the pH fine tuning as an underappreciated control factor in surface modification and self-assembly.

Adsorption of impurities in vegetable oil: A molecular modelling study

Here, the adsorption of impurity species from triglyceride solvent representing a model vegetable oil is studied using atomistic molecular dynamics simulations. We compare the adsorption of water, glycerol, oleic acid, monoolein, and two types of phospholipids on model silica adsorbents differing in their OH-group density, i.e. hydrogen bonding ability, quartz and cristobalite. We find that the species containing charged groups, phospholipids DOPC and DOPE, adsorb significantly stronger than the nonionic impurities. Secondary contribution to adsorption arises from hydrogen bonding capability of the impurity species, the silica surface, and also the triglyceride solvent: in general, more hydrogen bonding sites in impurity species leads to enhanced adsorption but hydrogen bonding with solvent competes for the available sites. Interestingly, adsorption is weaker on cristobalite even though it has a higher hydrogen bonding site density than quartz. This is because the hydrogen bonds can saturate each other on the adsorbent. The finding demonstrates that optimal adsorption response is obtained with intermediate adsorbent hydrogen bonding site densities. Additionally, we find that monoolein and oleic acid show a concentration driven adsorption response and reverse micelle like aggregate formation in bulk triglyceride solvent even in the absence of water. The findings offer insight into adsorption phenomena at inorganic adsorbent - apolar solvent interfaces and provide guidelines for enhanced design of adsorbent materials for example for vegetable oil purification.

Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles

The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.

Hybrid Atomistic and Coarse-Grained Model for Surfactants in Apolar Solvents

Here, we develop and verify the performance of a hybrid molecular modeling approach that combines coarse-grained apolar solvent and atomistic solute or polar solvent description, for example, for description of reverse micellar systems. The coarse-grained solvent model is directly applicable to organic solvents encompassing alkane, alkene, and fatty acid ester functional groups and connects directly to both standard united-atom GROMOS 53A6 and all-atom CHARMM27 force fields, as well as the atomistic detail water models compatible with these force fields. The different levels of description are coupled via explicit, unscaled electrostatics, and scaled mixing rules for dispersive interactions. The hybrid model is in near-quantitative agreement with fully atomistic simulations when combined with the CHARMM27 model but underestimates modestly surfactant aggregation when using GROMOS 53A6 united-atom description. The use of truncated electrostatics affords up to a 9-fold increase in computational speed without significant loss of accuracy. However, long-range electrostatic calculations and load imbalance at high core counts can significantly degrade the performance. We demonstrate the usability of the hybrid model by assessing the reverse micelle formation of a homologous series of nonionic glycerolipids via large-scale self-assembly simulations. The presented model is demonstrated here for accurate description of surfactant systems in apolar solvents, with and without also polar solvent (water) in the system. The formulation can be expected to describe well also other solute species or interfaces with an apolar solvent in an apolar environment.

Phospholipid-Based Reverse Micelle Structures in Vegetable Oil Modified by Water Content, Free Fatty Acid, and Temperature

Colloidal assemblies of phospholipids in oil are known to be highly sensitive to changes in system composition and temperature. Despite the fundamental biological and high industrial relevance of these aggregates, the mechanisms behind the structural changes, especially in real oils, are not well understood. In this work, small-angle X-ray scattering (SAXS) was combined with molecular dynamics simulations to characterize the effects of oleic acid, water, and temperature on self-assembled structures formed by lecithin in rapeseed oil. SAXS showed that adding water to the mixtures caused the precipitation of liquid-crystalline phases with lamellar or hexagonal geometry. The combination of SAXS and molecular dynamics simulations revealed that stable spherical reverse micelles in oil had a core radius of about 2 nm and consisted of approximately 60 phospholipids centered around a core containing water and sugars. The presence of oleic acid improved the stability of reverse micelles against precipitation due to the increase in the water concentration in oil by allowing the reverse micelle cores to expand and accommodate more water. The shape and size of the reverse micelles changed at high temperatures, and irreversible elongation was observed, especially in the presence of oleic acid. The findings show the interdependency of the structure of the reverse micellar aggregates on system composition, in particular, oleic acid and water, as well as temperature. The revealed characteristics of the self-assembled structures have significance in understanding and tuning the properties of vegetable oil-based emulsions, food products, oil purification, and drug delivery systems.

Towards capturing cellular complexity: combining encapsulation and macromolecular crowding in a reverse micelle

This paper studies the orientational structure and dynamics of multi-protein systems under confinement and discusses the implications on biological cells.

Hydration and Temperature Response of Water Mobility in Poly(diallyldimethylammonium)-Poly(sodium 4-styrenesulfonate) Complexes

The combination of all-atom molecular dynamics simulations with differential scanning calorimetry (DSC) has been exploited to investigate the influence of temperature and hydration on the water distribution and mobility in poly(diallyldimethylammonium) (PDADMA) and poly(sodium 4-styrenesulfonate) (PSS) complexes. The findings show that the vast majority of the water molecules hydrating the polyelectrolyte complexes (PECs) with 18-30 wt % hydration are effectively immobilized due to the strong interactions between the PE charge groups and water. Temperature and hydration were found to decrease similarly the fraction of strongly bound water. Additionally, at low hydration or at low temperatures, water motions become dominantly local vibrations and rotations instead of translational motion; translation dominance is recovered in a similar fashion by increase of both temperature and hydration. DSC experiments corroborate the simulation findings by showing that nonfreezing, bound water dominates in hydrated PECs at comparable hydrations. Our results raise attention to water as an equal variable to temperature in the design and engineering of stimuli-responsive polyelectrolyte materials and provide mechanistic explanation for the similarity.

Evaluating the Effects of Geometry and Charge Flux in Force Field Modeling

We apply a model for analyzing the importance of conformational charge flux to 11 molecules with the R-(CH₂) _n-R structure (R = Cl, F, OH, SH, COOH, CONH₂, and NH₂ and n = 4-6). Atomic charges were obtained by fitting to results from density functional theory calculations using the HLY procedure, and their geometry dependence is decomposed into contributions from changes in bond lengths, bond angles, and torsional angles. The torsional degrees of freedom are the main contribution to the conformational dependence of atomic charges and molecular dipole moments, but indirect effects due to changes in bond distances and angles account for ∼15% of the variations. While the magnitude of charge flux and geometry effects have been found to be independent of the number of internal degrees of freedom, the nature of the R- group has a moderate influence. The indirect effects are comparable for all of the R-groups and are approximately one-half the magnitude of the corresponding effects in peptide models. However, the magnitudes are different, yet the relative importance of geometry and charge flux effects are completely similar to those of the peptide models, which suggests that modeling the charge flux effects for changes in bond lengths, bond angles, and torsional angles should be considered for developing improved force fields.

Effects of 1-hexanol on C12E10 micelles: a molecular simulations and light scattering study

The micelles of the non-ionic C₁₂E₁₀ surfactant and 1-hexanol as an aqueous solution additives are studied toward the purpose of understanding the role of alcohol additives in tuning the characteristics of alkyl-ethoxylate micellar systems. Our dynamic light scattering and cloud point experiments show that the addition of hexanol induces a response similar to an increase of temperature. We associate the change with increased attraction between the micelles at low to moderate hexanol loadings and a potential increase of the aggregate size at a high hexanol-to-surfactant ratio. Detailed molecular dynamic simulation characterization shows that hexanol solubilizes to a micelle palisade layer when the hexanol-to-C₁₂E₁₀ ratio is less than or equal to 0.5 while swollen micelles, in which a part of hexanol forms an oil core, are present when the ratio increases above approximately 1.5. The simulations indicate that the surface of the micelles is rough. Formation of reverse hexanol structures akin to those found in bulk octanol is observed in the oil core. Molecular simulations associate the increase in attraction between micelles observed via the experiments with decreased chain density in the headgroup region. This density decrease is caused by hexanol molecules solubilized between neighbouring surfactants. Altogether, these findings provide detailed physical characterization of the effect of an archetypal solution additive, hexanol, on an alkyl ethoxylate micelle system. These findings could bear a significance in designing micellar and emulsion based systems with desired solution characteristics or properties for e.g. drug delivery, catalysis, or platforms for green chemistry reactions.

Aggregation response of triglyceride hydrolysis products in cyclohexane and triolein

Aggregation mechanism and the existence of cmc depend on apolar solvent quality and surfactant head group polarity.

The development of phytosterol-lecithin mixed micelles and organogels

We demonstrate that by mixing the phytosterol-ester oryzanol with lecithin in an organic solvent, both components may be dispersed at much higher concentrations than they may be individually. Dynamic light scattering and molecular dynamics simulations show that the mechanism for this is the formation of r ∼ 4 nm mixed micelles. Infrared spectroscopy and simulations suggest that these micelles are formed due in part to hydrogen bonding of the phosphate of the lecithin head-group, and the phenol group of the oryzanol. Rheology shows that by mixing these materials at an equimolar ratio, highly viscous suspensions are created. Furthermore, by adding water to these samples, a solid-like gel may be formed which offers mechanical properties close to those desired for a margarine type spread, whilst still solubilizing the oryzanol.

Effect of temperature, water content and free fatty acid on reverse micelle formation of phospholipids in vegetable oil

The self-assembly of phospholipids in oil, specifically lecithin in rapeseed oil, was investigated by combining experimental and computational methods The influence of temperature, water, and free fatty acids on the onset of lecithin aggregation in the rapeseed oil was determined using the 7,7,8,8 -tetracyanoquinodimethane dye (TCNQ) solubilization method and the size and shape of the self-assembled lecithin structures were investigated by small-angle X-ray scattering and cryogenic transmission electron microscopy. In the absence of excess water in the system (0.03wt-% water in oil), stable cylindrical lecithin reverse micelles were observed above the critical micelle concentration (CMC). Comparing the aggregation response in room temperature and at 70°C revealed that CMC decreased with increasing temperature. Furthermore, already a modest amount of added water (0.3wt-% water in oil) was sufficient to induce the formation of lamellar lecithin structures, that phase separated from the oil. In low water content, oleic acid suppressed the formation of lecithin reverse micelles whereas in the presence of more water, the oleic acid stabilized the reverse micelles. Consequently, more water was needed to induce phase separation in the presence of oleic acid. Molecular dynamics simulations indicated that the stabilizing effect of oleic acid resulted from oleic acid enhancing phospholipid solubilization in the oil by forming a solvating shell around the phosphate head group. The findings showed that the response of the mixed surfactant system is a delicate interplay of the different components and variables. The significance of the observations is that multiple parameters need to be controlled for desired system response, for example towards vegetable oil purification or phospholipid based microemulsions.

Reverse Micelles As Antioxidant Carriers: An Experimental and Molecular Dynamics Study

Water-in-oil microemulsions with biocompatible components were formulated to be used as carriers of natural antioxidants, such as hydroxytyrosol (HT) and gallic acid (GA). The system was composed of a mixture of natural surfactants, lecithin and monoglycerides, medium chain triglycerides, and aqueous phase. A dual approach was undertaken to study the structure and dynamics of these complicated systems. First, experimental data were collected by using adequate techniques, such as dynamic light scattering (DLS) and electron paramagnetic resonance (EPR) spectroscopy. Following this, a coarse-grained molecular dynamics (CGMD) study based on the experimental composition using the MARTINI force field was conducted. The simulations revealed the spontaneous formation of reverse micelles (RMs) starting from completely random initial conformations, underlying their enhanced thermodynamic stability. The location of the bioactive molecules, as well as the structure of the RM, were in accordance with the experimental findings. Furthermore, GA molecules were found to be located inside the water core, in contrast to the HT ones, which seem to lie at the surfactant interfacial layer. The difference in the antioxidants' molecular location was only revealed in detail from the computational analysis and explains the RM's swelling observed by GA in DLS measurements.

The Size of AOT Reverse Micelles

Reverse micelles (RMs) made from water and sodium bis(2-ethylhexyl) sulfosuccinate (AOT) are commonly studied experimentally as models of aqueous microenvironments. They are small enough for individual RMs to also be studied by molecular dynamics (MD) simulation, which yields detailed insight into their structure and properties. Although RM size is determined by the water loading ratio (i.e., the molar ratio of water to AOT), experimental measurements of RM size are imprecise and inconsistent, which is problematic when seeking to understand the relationship between water loading ratio and RM size, and when designing models for study by MD simulation. Therefore, a systematic study of RM size was performed by MD simulation with the aims of determining the size of an RM for a given water loading ratio, and of reconciling the results with experimental measurements. Results for a water loading ratio of 7.5 indicate that the interaction energy between AOT anions and other system components is at a minimum when there are 62 AOT anions in each RM. The minimum is due to a combination of attractive and repulsive electrostatic interactions that vary with RM size and the dielectric effect of available water. Overall, the results agree with a detailed analysis of previously published experimental data over a wide range of water loading ratios, and help reconcile seemingly discrepant experimental results. In addition, water loss and gain from an RM is observed and the mechanism of water exchange is outlined. This kind of RM model, which faithfully reproduces experimental results, is essential for reliable insights into the properties of RM-encapsulated materials.

On the Structural and Dynamical Properties of DOPC Reverse Micelles

The structure and dynamics of phospholipid reverse micelles are studied by molecular dynamics. We report all-atom unconstrained simulations of 1,2-dioleoyl-sn-phosphatidylcholine (DOPC) reverse micelles in benzene of increasing sizes, with water-to-surfactant number ratios ranging from W₀ = 1 to 16. The aggregation number, i.e., the number of DOPC molecules per reverse micelle, is determined to fit experimental light-scattering measurements of the reverse micelle diameter. The simulated reverse micelles are found to be approximately spherical. Larger reverse micelles (W₀ > 4) exhibit a layered structure with a water core and the hydration structure of DOPC phosphate head groups is similar to that found in phospholipid membranes. In contrast, the structure of smaller reverse micelles (W₀ ≤ 4) cannot be described as a series of concentric layers successively containing water, surfactant head groups, and surfactant tails, and the head groups are only partly hydrated and frequently present in the core. The dynamics of water molecules within the phospholipid reverse micelles slow down as the reverse micelle size decreases, in agreement with prior studies on AOT and Igepal reverse micelles. However, the average water reorientation dynamics in DOPC reverse micelles is found to be much slower than in AOT and Igepal reverse micelles with the same W₀ ratio. This is explained by the smaller water pool and by the stronger interactions between water and the charged head groups, as confirmed by the red-shift of the computed infrared line shape with decreasing W₀.

Role of hydration in phosphatidylcholine reverse micelle structure and gelation in cyclohexane: a molecular dynamics study

Molecular simulations reveal lipid headgroup hydration is at key role in organogel transition in lecithin–water–cyclohexane systems.