Dataset Information

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers.

ABSTRACT: Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4-5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations-except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.

SUBMITTER: Kim T

PROVIDER: S-EPMC3318127 | biostudies-literature | 2012 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers.

Kim Taehoon T Lee Kyu Il KI Morris Phillip P Pastor Richard W RW Andersen Olaf S OS Im Wonpil W

Biophysical journal 20120403 7

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain l ...[more]

PMID: 22500755

Similar Datasets

Project description:We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

Project description:Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.

Dataset Information

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers.

Publications

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets