Intersecting polymers in lipid bilayers: cliques, static order parameters and lateral diffusion

We have modelled a macrolipid polymer composed of lipid molecules (monomers) embedded in a lipid bilayer or monolayer and polymerized via their polar groups. Because of fluctuations perpendicular to the plane of the bilayer, the polar region occupied by the polymer chain possesses sufficient space so that the polymer might exhibit 'self-intersection' if its conformational state is projected onto the plane of the bilayer/monolayer. We represent the plane of the bilayer/monolayer by a triangular lattice. Each site can be occupied by a monomer or be empty (and thus occupied by one of the unpolymerizable lipids which make up the bilayer/monolayer). A macrolipid is represented by a sequence of N monomers connected by N-1 bonds. Bonds may be either short (connecting nearest neighbour monomers) or long (between second neighbour monomers), in accord with the average properties of the spacers between the polymerized lipids. We have carried out computer simulation of this system using the Carmesin-Kremer bond stretching algorithm. Although no two monomers can occupy the same site, bonds may cross each other. We analyzed the dependence of and approximately N2vc and + approximately N2 sigma c, where Nsc and Nmc are the number of bond-crossings in the same macrolipid ('self-crossing') or in two different macrolipids ('mutual-crossing'). For single macrolipids, we confirmed that vc = 3/4 and have found that sigma c approximately 0.52, which we consider supports that sigma c = 1/2. For the dense case with monomer concentration, c = 0.72, we found that vc = 1/2 and that sigma c approximately 0.52 supports that sigma c = 1/2. In the semi-dilute regime (c = 0.2) we found crossover behaviour, although sigma c = 1/2. The total number of bond crossings thus scale like N, independent of concentration. We studied the connectivity of the system by calculating the weight averaged cluster, or 'clique', size. Cliques are defined as being composed of all macrolipids which exhibit at least one crossing bond with one other member of the clique. We found that while the average clique contains about two macrolipids at low concentrations, the clique size approaches the maximum possible value at high concentrations if the macrolipids are sufficiently long. In the latter case a transition appears to occur as the macrolipid length increases. This transition occurs at length = 40 when c = 0.72. These observations should have experimental consequences for the viscoelastic properties of the system.(ABSTRACT TRUNCATED AT 400 WORDS)