Thermodynamics of glycophorin in phospholipid bilayer membranes

We have developed a model of glycophorin in a phospholipid bilayer membrane in order to study the thermodynamics of this system and to understand the detailed behavior of recent calorimetric data. We assume that the larger glycophorin polar group can be considered as either adopting a pancakelike conformation at the bilayer interface (D state) or be directed generally away from the interface (U state) [Ruppel, D., Kapitza, H.G., Galla, H.J., Sixl, F., & Sackmann, E. (1982) Biochim. Biophys. Acta 692, 1-17]. Lipid hydrocarbon chains are described either as excited (e state) with high energy and relatively many gauche conformers or as generally extended (g state) with low energy. We performed a Monte-Carlo simulation using the Glauber and Kawasaki procedures on a triangular lattice which represents the plane of half of the bilayer. Lattice sites can be occupied either by lipid hydrocarbon chains or by model glycophorin alpha-helical hydrophobic cores. The states D and U are represented by hexagons of different sizes in the plane of the lattice, and the hard core repulsion between two such polar groups is accounted for by forbidding hexagon-hexagon overlap. We have studied the effect of having the glycophorin polar group interact in various ways with the lipid bilayer. We find that the protein polar group in its D state interacts, either directly or indirectly, with the lipid bilayer so as to reduce the effective lateral pressure acting on the lipid hydrocarbon chains by about 3 dyn/cm. Polar groups in their U states do not reduce this lateral pressure.(ABSTRACT TRUNCATED AT 250 WORDS)