Equilibrium Distributions of Dipalmitoyl Phosphatidylcholine and Dilauroyl Phosphatidylcholine in a Mixed Lipid Bilayer: Atomistic Semigrand Canonical Ensemble Simulations

@article{dejoannis06a, abstract = {Abstract: Conventional molecular dynamics (MD) simulations are seriously limited by the slow rate of diffusive mixing in their ability to predict lateral distributions of different lipid types within mixed-lipid bilayers using atomistic models. A method to overcome this limitation, using configuration-bias Monte Carlo (MC) "mutation" moves to transform lipids from one type to another in dynamic equilibrium, is demonstrated in binary fluid-phase mixtures of lipids whose tails differ in length by four carbons. The hybrid MC-MD method operates within a semigrand canonical ensemble, so that an equilibrium composition of the mixture is determined by a constant difference in chemical potential () chosen for the components. Within several nanoseconds, bilayer structures initiated as pure dipalmitoyl phosphatidylcholine (DPPC) or pure dilauroyl phosphatidylcholine (DLPC) converge to a common composition and structure in independent simulations conducted at the same . Trends in bilayer thickness, area per lipid, density distributions across the bilayer, and order parameters have been investigated at three mixture compositions and compared with results from the pure bilayers at 323 K. The mixtures showed a moderate increase in DPPC acyl tail sites crossing the bilayer midplane relative to pure DPPC. Correlations between lateral positions of the two lipid types within or across the bilayer were found to be weak or absent. While the lateral distribution is consistent with nearly ideal mixing, the dependence of composition on indicates a positive excess free energy of mixing.}, address = {Department of Chemistry and Emerson Center for Scientific Computation, Emory University, 1515 Dickey Drive, Atlanta, Georgia, 30322}, author = {Dejoannis, J. and Jiang, Y. and Yin, F. and Kindt, J. T. }, citeulike-article-id = {1790168}, doi = {http://dx.doi.org/10.1021/jp065734y}, journal = {J. Phys. Chem. B}, keywords = {simulation}, month = {December}, number = {51}, pages = {25875--25882}, posted-at = {2007-10-19 19:49:15}, priority = {2}, title = {Equilibrium Distributions of Dipalmitoyl Phosphatidylcholine and Dilauroyl Phosphatidylcholine in a Mixed Lipid Bilayer: Atomistic Semigrand Canonical Ensemble Simulations}, url = {http://dx.doi.org/10.1021/jp065734y}, volume = {110}, year = {2006} }

TY - JOUR ID - dejoannis06a L3 - citeulike-article-id:1790168 TI - Equilibrium Distributions of Dipalmitoyl Phosphatidylcholine and Dilauroyl Phosphatidylcholine in a Mixed Lipid Bilayer: Atomistic Semigrand Canonical Ensemble Simulations JF - J. Phys. Chem. B VL - 110 IS - 51 SP - 25875 EP - 25882 AD - Department of Chemistry and Emerson Center for Scientific Computation, Emory University, 1515 Dickey Drive, Atlanta, Georgia, 30322 N2 - Abstract: Conventional molecular dynamics (MD) simulations are seriously limited by the slow rate of diffusive mixing in their ability to predict lateral distributions of different lipid types within mixed-lipid bilayers using atomistic models. A method to overcome this limitation, using configuration-bias Monte Carlo (MC) "mutation" moves to transform lipids from one type to another in dynamic equilibrium, is demonstrated in binary fluid-phase mixtures of lipids whose tails differ in length by four carbons. The hybrid MC-MD method operates within a semigrand canonical ensemble, so that an equilibrium composition of the mixture is determined by a constant difference in chemical potential () chosen for the components. Within several nanoseconds, bilayer structures initiated as pure dipalmitoyl phosphatidylcholine (DPPC) or pure dilauroyl phosphatidylcholine (DLPC) converge to a common composition and structure in independent simulations conducted at the same . Trends in bilayer thickness, area per lipid, density distributions across the bilayer, and order parameters have been investigated at three mixture compositions and compared with results from the pure bilayers at 323 K. The mixtures showed a moderate increase in DPPC acyl tail sites crossing the bilayer midplane relative to pure DPPC. Correlations between lateral positions of the two lipid types within or across the bilayer were found to be weak or absent. While the lateral distribution is consistent with nearly ideal mixing, the dependence of composition on indicates a positive excess free energy of mixing. KW - simulation AU - Dejoannis, J AU - Jiang, Y AU - Yin, F AU - Kindt, JT PY - 2006/12/28/ UR - http://dx.doi.org/10.1021/jp065734y ER -