Using molecular dynamics to study liquid phase behavior: simulations of the ternary sodium laurate/sodium oleate/water system.

The prediction of surfactant phase behavior has applications in a wide range of areas. An accurate modeling of liquid phase behavior can aid our understanding of colloidal process or be used to design phases that respond in a defined way to their environment. In this work, we use molecular dynamics to model the phase behavior of the ternary sodium laurate/sodium oleate/water system and compare the simulation results to experimental data. Simulations were performed with the GROMOS 53A6 united-atom force field and cover the entire ternary phase diagram, producing micellar, hexagonal, and lamellar phases. The aggregate simulation time for the 33 simulations performed during this study is 4.4 μs. We find that the simulations were able to model the experimentally observed liquid phase behavior accurately, showing that the carboxylate and lipid parameters of the 53A6 force field give very good quality results for the in silico prediction of liquid system phase behavior.

[1]  B. Smit,et al.  Molecular simulation of the DMPC-cholesterol phase diagram. , 2010, The journal of physical chemistry. B.

[2]  D. Tieleman,et al.  Molecular Dynamics Simulation of a Polyunsaturated Lipid Bilayer Susceptible to Lipid Peroxidation , 2004 .

[3]  Dusanka Janezic,et al.  Liquid-ordered phase formation in cholesterol/sphingomyelin bilayers: all-atom molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[4]  Markus Christen,et al.  The GROMOS software for biomolecular simulation: GROMOS05 , 2005, J. Comput. Chem..

[5]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[6]  F. Stanley,et al.  Rapid determination of surfactant critical micelle concentrations using pressure-driven flow with capillary electrophoresis instrumentation. , 2009, Journal of chromatography. A.

[7]  P. Hünenberger,et al.  Simulating the Transition between Gel and Liquid-Crystal Phases of Lipid Bilayers: Dependence of the Transition Temperature on the Hydration Level. , 2010, Journal of Chemical Theory and Computation.

[8]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[9]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[10]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[11]  Dallas B. Warren,et al.  Using polymeric precipitation inhibitors to improve the absorption of poorly water-soluble drugs: A mechanistic basis for utility , 2010, Journal of drug targeting.

[12]  Roland Faller,et al.  Simulation of domain formation in DLPC-DSPC mixed bilayers. , 2004, Langmuir : the ACS journal of surfaces and colloids.

[13]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[14]  Lukas D. Schuler,et al.  Molecular dynamics simulation of n-dodecyl phosphate aggregate structures , 2001, European Biophysics Journal.

[15]  I. Vattulainen,et al.  Influence of cis double-bond parametrization on lipid membrane properties: how seemingly insignificant details in force-field change even qualitative trends. , 2008, The Journal of chemical physics.

[16]  P. Coppock,et al.  Determination of phase transition temperatures for atomistic models of lipids from temperature-dependent stripe domain growth kinetics. , 2010, The journal of physical chemistry. B.

[17]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[18]  Gordon J. T. Tiddy,et al.  Phase behaviour of polyoxyethylene surfactants with water. Mesophase structures and partial miscibility (cloud points) , 1983 .

[19]  Haiyan Liu,et al.  Refining the description of peptide backbone conformations improves protein simulations using the GROMOS 53A6 force field , 2009, J. Comput. Chem..

[20]  C. Pouton,et al.  Enhancing intestinal drug solubilisation using lipid-based delivery systems. , 2008, Advanced drug delivery reviews.

[21]  Simon Hanna,et al.  Atomistic simulation of a model liquid crystal. , 2006, The Journal of chemical physics.

[22]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[23]  C. Macosko,et al.  Rheology of highly concentrated anionic surfactants , 2006 .

[24]  David K. Chalmers,et al.  Molecular dynamics simulations of spontaneous bile salt aggregation , 2006 .

[25]  Colin W Pouton,et al.  Formulation of poorly water-soluble drugs for oral administration: physicochemical and physiological issues and the lipid formulation classification system. , 2006, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[26]  V. Cherezov,et al.  Rational design of lipid for membrane protein crystallization. , 2004, Journal of structural biology.

[27]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[28]  Dallas B. Warren,et al.  Structure and dynamics of glyceride lipid formulations, with propylene glycol and water. , 2009, Molecular pharmaceutics.

[29]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[30]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[31]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[32]  Carl-Johan Högberg,et al.  A molecular dynamics investigation of the influence of hydration and temperature on structural and dynamical properties of a dimyristoylphosphatidylcholine bilayer. , 2006, The journal of physical chemistry. B.

[33]  G. Edwards,et al.  Evaluation of the impact of surfactant digestion on the bioavailability of danazol after oral administration of lipidic self-emulsifying formulations to dogs. , 2008, Journal of pharmaceutical sciences.

[34]  Claudio Zannoni,et al.  Towards in silico liquid crystals. Realistic transition temperatures and physical properties for n-cyanobiphenyls via molecular dynamics simulations. , 2009, Chemphyschem : a European journal of chemical physics and physical chemistry.

[35]  P. Coppock,et al.  Atomistic simulations of mixed-lipid bilayers in gel and fluid phases. , 2009, Langmuir : the ACS journal of surfaces and colloids.

[36]  J. Mccammon,et al.  Molecular Dynamics Simulations of a Polyalanine Octapeptide under Ewald Boundary Conditions: Influence of Artificial Periodicity on Peptide Conformation , 2000 .

[37]  Christopher J H Porter,et al.  Formulation of lipid-based delivery systems for oral administration: materials, methods and strategies. , 2008, Advanced drug delivery reviews.

[38]  Zhi‐Wu Yu,et al.  Structural characterization on the gel to liquid-crystal phase transition of fully hydrated DSPC and DSPE bilayers. , 2009, The journal of physical chemistry. B.

[39]  B. Smit,et al.  Phase behavior of model lipid bilayers. , 2005, The journal of physical chemistry. B.

[40]  Dmitry S Yakovlev,et al.  Molecular dynamics simulations of mixed cationic/anionic wormlike micelles. , 2007, Langmuir : the ACS journal of surfaces and colloids.

[41]  Siewert J Marrink,et al.  Lipids on the move: simulations of membrane pores, domains, stalks and curves. , 2009, Biochimica et biophysica acta.

[42]  B. Ladanyi,et al.  Molecular dynamics simulation of aerosol-OT reverse micelles. , 2009, The journal of physical chemistry. B.

[43]  D Peter Tieleman,et al.  Molecular dynamics simulation of spontaneous membrane fusion during a cubic-hexagonal phase transition. , 2002, Biophysical journal.

[44]  C. Pouton,et al.  Design of lipid-based formulations for oral administration of poorly water-soluble drugs: precipitation of drug after dispersion of formulations in aqueous solution. , 2009, Journal of pharmaceutical sciences.

[45]  A. Sum,et al.  Molecular studies of the gel to liquid-crystalline phase transition for fully hydrated DPPC and DPPE bilayers. , 2007, Biochimica et biophysica acta.

[46]  J. Kindt,et al.  Molecular dynamics simulations of glycocholate-oleic acid mixed micelle assembly. , 2010, Langmuir : the ACS journal of surfaces and colloids.

[47]  G. Edwards,et al.  Increasing the Proportional Content of Surfactant (Cremophor EL) Relative to Lipid in Self-emulsifying Lipid-based Formulations of Danazol Reduces Oral Bioavailability in Beagle Dogs , 2007, Pharmaceutical Research.

[48]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[49]  Mikko Karttunen,et al.  Structural properties of ionic detergent aggregates: a large-scale molecular dynamics study of sodium dodecyl sulfate. , 2007, The journal of physical chemistry. B.

[50]  Daniel E. Rosenfeld,et al.  Dynamics of water confined within reverse micelles. , 2006, The journal of physical chemistry. B.

[51]  P. V. Coveney,et al.  Large scale molecular dynamics simulation of self-assembly processes in short and long chain cationic surfactants , 1999 .