Application of finite-size scaling to the Pink model for lipid bilayers

Main transition in the Pink membrane model: finite-size scaling and the influence of surface roughness

We consider the main transition in single-component membranes using computer simulations of the Pink model [D. A. Pink et al., Biochemistry 19, 349 (1980)]. We first show that the accepted parameters of the Pink model yield a main transition temperature that is systematically below experimental values. This resolves an issue that was first pointed out by Corvera and co-workers [Phys. Rev. E 47, 696 (1993)]. In order to yield the correct transition temperature, the strength of the van der Waals coupling in the Pink model must be increased; by using finite-size scaling, a set of optimal values is proposed. We also provide finite-size scaling evidence that the Pink model belongs to the universality class of the two-dimensional Ising model. This finding holds irrespective of the number of conformational states. Finally, we address the main transition in the presence of quenched disorder, which may arise in situations where the membrane is deposited on a rough support. In this case, we observe a stable multidomain structure of gel and fluid domains, and the absence of a sharp transition in the thermodynamic limit.

First order melting transitions of highly ordered dipalmitoyl phosphatidylcholine gel phase membranes in molecular dynamics simulations with atomistic detail

Molecular dynamics simulations with atomistic detail of the gel phase and melting transitions of dipalmitoyl phosphatidylcholine bilayers in water reveal the dependency of many thermodynamic and structural parameters on the initial system ordering. We quantitatively compare different methods to create a gel phase system and we observe that a very high ordering of the gel phase starting system is necessary to observe behavior which reproduces experimental data. We performed heating scans with speeds down to 0.5 K/ns and could observe sharp first order phase transitions. Also, we investigated the transition enthalpy as the natural intrinsic parameter of first order phase transitions, and obtained a quantitative match with experimental values. Furthermore, we performed systematic investigations of the statistical distribution and heating rate dependency of the microscopic phase transition temperature.

Histogram method to obtain heat capacities in lipid monolayers, curved bilayers, and membranes containing peptides

Lipid monolayer chain melting transitions were simulated using a two-state Doniach model, and experimental melting profiles of lipid vesicles were analyzed. We sampled the information of a Monte Carlo simulation into a single broad histogram containing complete information about the distribution of states. The information of the monolayer histogram was first used to calculate the melting behavior of a bilayer constructed from two uncoupled monolayers. We then fitted calorimetric heat profiles of various preparations of dipalmitoyl phosphatidylcholine vesicles. This analysis was extended to lipid bilayers. A fixed mean bilayer curvature was shown to result in a broadening of bilayer melting profiles. We furthermore used the histogram method to obtain the chain melting behavior of simple lipid-peptide mixtures.

Simulation of the gel-fluid transition in a membrane composed of lipids with two connected acyl chains: application of a dimer-move step

Phospholipids have been treated as dimers on a hexagonal lattice, and a move has been introduced that allows the dimers to move and change their orientation on the lattice. Simulations have been performed in which phospholipid chains have been treated as being either independent or infinitely coupled thermodynamically with regard to their conformational state. Both types of simulation have reproduced well experimental heat-capacity curves of dipalmitoyl phosphatidylcholine small unilamellar vesicles. Apart from a different gel-fluid interaction parameter and a different number of unlike nearest-neighbor contacts, most of the averages and thermodynamic quantities were essentially the same in the two types of simulation. These results indicate that the transition is not first order and validate those of previous Monte Carlo simulations that have neglected the dimeric nature of phospholipids in the sense that they show that for the thermotropic transition the approximation of phospholipids as monomers is valid.

A Monte Carlo simulation study of protein-induced heat capacity changes and lipid-induced protein clustering

Monte Carlo simulations were used to describe the interaction of peripheral and integral proteins with lipids in terms of heat capacity profiles and protein distribution. The simulations were based on a two-state model for the lipid, representing the lipid state as being either gel or fluid. The interaction between neighboring lipids has been taken into account through an unlike nearest neighbor free energy term delta omega, which is a measure of the cooperativity of the lipid transition. Lipid/protein interaction was considered using the experimental observation that the transition midpoints of lipid membranes are shifted upon protein binding, a thermodynamic consequence of different binding constants of protein with fluid or gel lipids. The difference of the binding free energies was used as an additional parameter to describe lipid-protein interaction. The heat capacity profiles of lipid/protein complexes could be well described for both peripheral and integral proteins. Binding of proteins results in a shift and an asymmetric broadening of the melting profile. The model results in a coexistence of gel and fluid lipid domains in the proximity of the thermotropic transition. As a consequence, bound peripheral proteins aggregate in the temperature range of the lipid transition. Integral proteins induce calorimetric melting curves that are qualitatively different from that of peripheral proteins and aggregate in either gel or liquid crystalline lipid phase. The results presented here are in good agreement with calorimetric experiments on lipid-protein complexes and have implementations for the functional control of proteins.

The computer as a laboratory for the physical chemistry of membranes

A mini-review is given of some recent advances in the use of computer-simulation approaches to the study of physico-chemical properties of lipid bilayers and biological membranes. The simulations are based on microscopic molecular interaction models as well as random-surface models of fluid membranes. Particular emphasis is put on those properties that are controlled by the many-particle character of the lamellar membrane, i.e. correlations and fluctuations in density, composition and large-scale conformational structure. It is discussed how dynamic membrane heterogeneity arises and how it is affected by various molecular species interacting with membranes, such as cholesterol, drugs, insecticides, as well as polypeptides and integral membrane proteins. The influence of bending rigidity and osmotic-pressure gradients on large-scale membrane conformation and topology is described.