Computer simulation of the main gel-fluid phase transition of lipid bilayers
Research output: Contribution to journal › Journal article › Research › peer-review
Standard
Computer simulation of the main gel-fluid phase transition of lipid bilayers. / Mouritsen, O. G.; Boothroyd, A.; Harris, R.; Jan, N.; Lookman, T.; MacDonald, L.; Pink, D. A.; Zuckermann, M. J.
In: The Journal of Chemical Physics, Vol. 79, No. 4, 1983, p. 2027-2041.Research output: Contribution to journal › Journal article › Research › peer-review
Harvard
APA
Vancouver
Author
Bibtex
}
RIS
TY - JOUR
T1 - Computer simulation of the main gel-fluid phase transition of lipid bilayers
AU - Mouritsen, O. G.
AU - Boothroyd, A.
AU - Harris, R.
AU - Jan, N.
AU - Lookman, T.
AU - MacDonald, L.
AU - Pink, D. A.
AU - Zuckermann, M. J.
PY - 1983
Y1 - 1983
N2 - Monte Carlo techniques have been applied to a study of two related quasi-two-dimensional microscopic interaction models which describe the phase behavior of phospholipid bilayers. The two models are Ising-like lattice models in which (a) the acyl chains of the phospholipids interact via anisotropic van der Waals forces and (b) the rotational isomerism of the chains is accounted for by two and ten selected conformational states, respectively. Monte Carlo experiments are performed on both models so as to determine whether the static thermodynamic properties of lipid bilayers are most accurately represented by a simple two state gel-fluid concept or whether a more complicated melting process involving intermediate states takes place. To this purpose, the temperature dependence of several static thermodynamic properties has been calculated for both models. This includes the chain cross-sectional area, the internal and free energies, the coherence length, the lateral compressibility, and the specific heat. Particular care has been devoted to the transition region, since no analytical results are available in this region for either model. The comparison between the Monte Carlo results for the two models demonstrates that, whereas the two-state model has a first-order transition with jumplike behavior in the transition region the ten-state model exhibits a first-order transition associated with a closed hysteresis loop. Next, the Monte Carlo results for cross-sectional areas per lipid chain, coherence lengths and lateral compressibilities are discussed in the context of experimental results for dipalmitoyl phosphatidylcholine (DPPC). A detailed comparison is made with the results of molecular field calculations throughout the paper. Finally, a Monte Carlo analysis of bilayers composed of both DPPC and cholesterol shows that a two-state model does not adequately describe the thermodynamic behavior of lipid-cholesterol mixtures implying that intermediate states have to be introduced to account for the experimental data.
AB - Monte Carlo techniques have been applied to a study of two related quasi-two-dimensional microscopic interaction models which describe the phase behavior of phospholipid bilayers. The two models are Ising-like lattice models in which (a) the acyl chains of the phospholipids interact via anisotropic van der Waals forces and (b) the rotational isomerism of the chains is accounted for by two and ten selected conformational states, respectively. Monte Carlo experiments are performed on both models so as to determine whether the static thermodynamic properties of lipid bilayers are most accurately represented by a simple two state gel-fluid concept or whether a more complicated melting process involving intermediate states takes place. To this purpose, the temperature dependence of several static thermodynamic properties has been calculated for both models. This includes the chain cross-sectional area, the internal and free energies, the coherence length, the lateral compressibility, and the specific heat. Particular care has been devoted to the transition region, since no analytical results are available in this region for either model. The comparison between the Monte Carlo results for the two models demonstrates that, whereas the two-state model has a first-order transition with jumplike behavior in the transition region the ten-state model exhibits a first-order transition associated with a closed hysteresis loop. Next, the Monte Carlo results for cross-sectional areas per lipid chain, coherence lengths and lateral compressibilities are discussed in the context of experimental results for dipalmitoyl phosphatidylcholine (DPPC). A detailed comparison is made with the results of molecular field calculations throughout the paper. Finally, a Monte Carlo analysis of bilayers composed of both DPPC and cholesterol shows that a two-state model does not adequately describe the thermodynamic behavior of lipid-cholesterol mixtures implying that intermediate states have to be introduced to account for the experimental data.
U2 - 10.1063/1.445987
DO - 10.1063/1.445987
M3 - Journal article
AN - SCOPUS:0001319712
VL - 79
SP - 2027
EP - 2041
JO - The Journal of Chemical Physics
JF - The Journal of Chemical Physics
SN - 0021-9606
IS - 4
ER -
ID: 238392364