Molecular dynamics study of peptide-bilayer adsorption

Two 6-ns simulations of the somatostatin analog sandostatin and a 1-palmityl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer are presented. In the first simulation, the peptide was placed in a region of bulk water density and allowed to spontaneously move toward and bind to the bilayer surface. An attractive force between the peptide and bilayer drove the binding process, which was opposed by a significant frictional force caused by the solvent (water). During the approach of the peptide toward the bilayer the area of the interacting surface between the species was inversely proportional to the distance between them, supporting the application of such a relationship in continuum calculations of peptide-bilayer binding free energies. In the second simulation, the N-terminus of the surface-bound peptide was deprotonated. Consistent with experiment, this strengthened interactions between the peptide and the bilayer. Details of both peptide-bilayer complexes, including the orientation, percent buried surface area, and orientation of the lipid headgroups are in good agreement with those obtained from experiment. The location of the different side chains in the bilayer is in direct correlation with an interfacial hydrophobicity scale developed using model peptides. The aromatic side chains of the Phe and Trp residues all lie flat with respect to the bilayer surface in both complexes. Changes in lipid and water ordering due to peptide binding suggest a possible domination of lipophobic over hydrophobic effects, as proposed by other workers. Where appropriate, peptide and lipid properties in the bound states are compared with separate simulations of sandostatin and the bilayer in water, respectively, so as to monitor the response of the system to the binding process.