Lateral diffusion of gramicidin S, M-13 coat protein and glycophorin in bilayers of saturated phospholipids: Mean field and Monte Carlo studies

Abstract

We have developed a general model that relates the lateral diffusion coefficient of one isolated large intrinsic molecule (mol. wt. ≳1000) in a phosphatidylcholine bilayer to the static lipid hydrocarbon chain order. We have studied how protein lateral diffusion can depend upon protein-lipid interactions but have not investigated possible non-specific contributions from gel-state lattice defects. The model has been used in Monte Carlo simulations or in mean-field approximations to study the lateral diffusion coefficients of Gramicidin S, the M-13 coat protein and glycophorin in dimyristoyl- and dipalmitoylphosphatidylcholine (DMPC and DPPC) bilayers as functions of temperature. Our calculated lateral diffusion coefficients for Gramicidin S and the M-13 coat protein are in good agreement with what has been observed and suggest that Gramicidin S is in a dimeric form in DMPC bilayers. In the case of glycophorin we find that the ‘ice breaker’ effect can be understood as a consequence of perturbation of the lipid polar region around the protein. In order to understand this effect is is necessary that the protein hydrophilic section perturb the polar regions of at least approx. 24 lipid molecules, in good agreement with the numbers of 29–30 measured using 31P-NMR. Because of lipid-lipid interactions this effect extends itself out to four or five lipid layers away from the protein so that the hydrocarbon chains of between approx. 74 and approx. 108 lipid molecules are more disordered in the gel phase, so contributing less to the transition enthalpy, in agreement with the numbers of 80–100 deduced from differential scanning calorimetry (DSC). An understanding of the abrupt change in the diffusion coefficient at a temperature below the main bilayer transition temperature requires an additional mechanism. We propose that this change may be a consequence of a ‘coupling-uncoupling’ transition involving the protein hydrophilic section and the lipid polar regions, which may be triggered by the lipid bilayer pretransition. Our calculation of the average number of gauche bonds per lipid chain as a function of temperature and distance away from an isolated polypeptide or integral protein shows the extent of statically disordered lipid around such molecules. The range of this disorder depends upon temperature, particularly near the main transition.