Free energy simulations of amylin I26P mutation in a lipid bilayer

The amylin peptide in a dioleoylphosphatidylcholine (DOPC) bilayer is studied using united atom molecular dynamics (MD) simulations. Dynamics and transport properties of the peptide and the phospholipid bilayer are investigated. The lateral diffusion of DOPC is in the order of 10−8 cm2 s−1, which is in agreement with the experimental results. The order parameter and density profile for phospholipid molecules in the bilayer are calculated. The secondary structure of amylin peptide shows that the amino acids in two terminals are structureless and two α-helical segments in the peptide are connected through an unstructured link. This structure is similar to the experimental structure in the membrane-mimicking media. Free energy calculations of the Ile26 → Pro mutation in the amylin peptide are performed in the bilayer and in aqueous solution using molecular dynamics simulations and a thermodynamic cycle. It is shown that in the mutated peptide in aqueous solution, the α-helix structure changes to a 5-helix, whereas this configuration is preserved in the bilayer environment. It is interesting that the accessible surface area increases for hydrophobic residues in the bilayer and for hydrophilic residues in aqueous solution as the coupling parameter changes from 0 to 1. These results are significant to understanding the aggregation mechanism of human amylin monomers in membranes to the dimers, trimers, oligomers, and fibrils associated with the type 2 diabetes at the atomic level.