1-palmitoyl-2-oleoylphosphatidylcholine and 1-2-diarachidonoyl-glycero-3-phosphocholine

Coarse-grained molecular dynamics simulations are used to calculate the free energies of transfer of miltefosine, an alkylphosphocholine anticancer agent, from water to lipid bilayers to study its mechanism of interaction with biological membranes. We consider bilayers containing lipids with different degrees of unsaturation: dipalmitoylphosphatidylcholine (DPPC, saturated, containing 0%, 10%, and 30% cholesterol), dioleoylphosphatidylcholine (DOPC, diunsaturated), palmitoyloleoylphosphatidylcholine (POPC, monounsaturated), diarachidonoylphosphatidylcholine (DAPC, polyunsaturated), and dilinoleylphosphatidylcholine (DUPC, polyunsaturated). These free energies, calculated using umbrella sampling, were used to compute the partition coefficients (K) of miltefosine between water and the lipid bilayers. The K values for the bilayers relative to that of pure DPPC were found to be 5.3 (DOPC), 7.0 (POPC), 1.0 (DAPC), 2.2 (DUPC), 14.9 (10% cholesterol), and 76.2 (30% cholesterol). Additionally, we calculated the free energy of formation of miltefosine-cholesterol complexes by pulling the surfactant laterally in the DPPC + 30% cholesterol system. The free energy profile that we obtained provides further evidence that miltefosine tends to associate with cholesterol and has a propensity to partition into lipid rafts. We also quantified the kinetics of the transport of miltefosine through the various bilayers by computing permeance values. The highest permeance was observed in DUPC bilayers (2.28 × 10(-2) m/s) and the lowest permeance in the DPPC bilayer with 30% cholesterol (1.10 × 10(-7) m/s). Our simulation results show that miltefosine does indeed interact with lipid rafts, has a higher permeability in polyunsaturated, loosely organized bilayers, and has higher flip-flop rates in specific regions of cellular membranes.

Topics: 1,2-Dipalmitoylphosphatidylcholine; Antineoplastic Agents; Cholesterol; Kinetics; Lipid Bilayers; Membrane Microdomains; Molecular Dynamics Simulation; Phosphatidylcholines; Phosphorylcholine; Thermodynamics; Water

Hydrophobic mismatch which is defined as the difference between the lipid hydrophobic thickness and the peptide hydrophobic length is known to be responsible in altering the lipid/protein dynamics. Gramicidin A (gA), a 15 residue β helical peptide which is well recognized to form ion conducting channels in lipid bilayer, may change its structure and function in a hydrophobic mismatched condition. We have performed molecular dynamics simulations of gA dimer in phospholipid bilayers to investigate whether or not the conversion from channel to non-channel form of gA dimer would occur under extreme negative hydrophobic mismatch. By varying the length of lipid bilayers from DLPC (1, 2-Dilauroyl-sn-glycero-3-phosphocholine) to DAPC (1, 2-Diarachidoyl-sn-glycero-3-phosphocholine), a broad range of mismatch was considered from nearly matching to extremely negative. Our simulations revealed that though the ion-channel conformation is retained by gA under a lesser mismatched situation, in extremely negative mismatched situation, in addition to bilayer thinning, the conformation of gA is changed and converted to a non-channel one. Our results demonstrate that although the channel conformation of Gramicidin A is the most stable structure, it is possible for gA to change its conformation from channel to non-channel depending upon the local environment of host bilayers.

Topics: Dimyristoylphosphatidylcholine; Gramicidin; Hydrophobic and Hydrophilic Interactions; Ion Channels; Lipid Bilayers; Molecular Dynamics Simulation; Phosphatidylcholines; Protein Structure, Secondary; Thermodynamics