Abstract
Model biological membranes consisting of peptide/lipid-bilayer complexes can nowadays be studied by classical molecular dynamics (MD) simulations at atomic detail. In most cases, the simulation starts with an assumed state of a peptide in a preformed bilayer, from which equilibrium configurations are difficult to obtain due to a relatively slow molecular diffusion. As an alternative, we propose an extension of reported work on the self-organization of unordered lipids into bilayers, consisting of including a peptide molecule in the initial random configuration to obtain a membrane-bound peptide simultaneous to the formation of the lipid bilayer. This strategy takes advantage of the fast reorganization of lipids, among themselves and around the peptide, in an aqueous environment. Model peptides of different hydrophobicity, CH3-CO-W2L18W2-NH2 (WL22) and CH3-CO-W2A18W2-NH2 (WA22), in dipalmitoyl-phosphatidylcholine (DPPC) are used as test cases. In the equilibrium states of the peptide/membrane complexes, achieved in time ranges of 50-100 ns, the two peptides behave as expected from experimental and theoretical studies. The strongly hydrophobic WL22 is inserted in a transmembrane configuration and the marginally apolar, alanine-based WA22 is found in two alternative states: transmembrane inserted or parallel to the membrane plane, embedded close to the bilayer interface, with similar stability. This shows that the spontaneous assembly of peptides and lipids is an unbiased and reliable strategy to produce and study models of equilibrated peptide/lipid complexes of unknown membrane-binding mode and topology.