Folding and insertion thermodynamics of the transmembrane WALP peptide

• Published in: The Journal of chemical physics Vol. 143; no. 24; p. 243127
• Main Authors: Bereau, Tristan, Bennett, W. F. Drew, Pfaendtner, Jim, Deserno, Markus, Karttunen, Mikko
• Format: Journal Article
• Language: English
• Published: United States American Institute of Physics 28.12.2015
• Subjects: Amino acids; Anchors; Cell Membrane - chemistry; Chains; COMPARATIVE EVALUATIONS; Computer simulation; COMPUTERIZED SIMULATION; COOPERATIVES; Dependence; Folding; FREE ENERGY; Helices; INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; Insertion; INTERFACES; LAYERS; Lipid Bilayers - chemistry; LIPIDS; MEMBRANE PROTEINS; MEMBRANES; Molecular Dynamics Simulation; PEPTIDES; Peptides - chemistry; POTENTIALS; PROBES; Protein Folding; Protein Structure, Secondary; Proteins; RESIDUES; STABILITY; THERMODYNAMICS; WATER
• ISSN: 0021-9606; 1089-7690
• DOI: 10.1063/1.4935487

The anchor of most integral membrane proteins consists of one or several helices spanning the lipid bilayer. The WALP peptide, GWW(LA)n (L)WWA, is a common model helix to study the fundamentals of protein insertion and folding, as well as helix-helix association in the membrane. Its structural properties have been illuminated in a large number of experimental and simulation studies. In this combined coarse-grained and atomistic simulation study, we probe the thermodynamics of a single WALP peptide, focusing on both the insertion across the water-membrane interface, as well as folding in both water and a membrane. The potential of mean force characterizing the peptide’s insertion into the membrane shows qualitatively similar behavior across peptides and three force fields. However, the Martini force field exhibits a pronounced secondary minimum for an adsorbed interfacial state, which may even become the global minimum—in contrast to both atomistic simulations and the alternative PLUM force field. Even though the two coarse-grained models reproduce the free energy of insertion of individual amino acids side chains, they both underestimate its corresponding value for the full peptide (as compared with atomistic simulations), hinting at cooperative physics beyond the residue level. Folding of WALP in the two environments indicates the helix as the most stable structure, though with different relative stabilities and chain-length dependence.