Molecular Dynamics Simulations of PIP2 and PIP3 in Lipid Bilayers: Determination of Ring Orientation, and the Effects of Surface Roughness on a Poisson-Boltzmann Description

被引:56
|
作者
Li, Zheng [2 ]
Venable, Richard M. [1 ]
Rogers, Laura A. [2 ]
Murray, Diana [2 ]
Pastor, Richard W. [1 ]
机构
[1] NHLBI, Lab Computat Biol, NIH, Bethesda, MD 20892 USA
[2] Columbia Univ, Dept Pharmacol, Presbyterian Hosp, New York, NY USA
基金
美国国家卫生研究院;
关键词
ELECTROSTATIC SEQUESTRATION; LIQUID/LIQUID INTERFACES; COMPUTER-SIMULATION; INOSITOL LIPIDS; MEMBRANES; BINDING; MODEL; PEPTIDES; HEADGROUP; PROTEINS;
D O I
10.1016/j.bpj.2009.04.037
中图分类号
Q6 [生物物理学];
学科分类号
071011 ;
摘要
Molecular dynamics (MD) simulations of phosphatidylinositol (4,5)-bisphosphate (PIP2) and phosphatidylinositol (3,4,5) -trisphosphate (PIP3) in 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC) bilayers indicate that the inositol rings are tilted similar to 40 degrees with respect to the bilayer surface, as compared with 17 degrees for the P-N vector of POPC. Multiple minima were obtained for the ring twist (analogous to roll for an airplane). The phosphates at position 1 of PIP2 and PIP3 are within an Angstrom of the plane formed by the phosphates of POPC; lipids in the surrounding shell are depressed by 0.5-0.8 angstrom, but otherwise the phosphoinositides do not substantially perturb the bilayer. Finite size artifacts for ion distributions are apparent for systems of similar to 26 waters/lipid, but, based on simulations with a fourfold increase of the aqueous phase, the phosphoinositide positions and orientations do not show significant size effects. Electrostatic potentials evaluated from Poisson-Boltzmann (PB) calculations show a strong dependence of potential height and ring orientation, with the maxima on the -25 mV surfaces (17.1 +/- 0.1 angstrom for PIP2 and 19.4 +/- 0.3 angstrom for PIP3) occurring near the most populated orientations from MD. These surfaces are well above the background height of 10 angstrom estimated for negatively charged cell membranes, as would be expected for lipids involved in cellular signaling. PB calculations on microscopically flat bilayers yield similar maxima as the MD-based (microscopically rough) systems, but show less fine structure and do not clearly indicate the most probable regions. Electrostatic free energies of interaction with pentalysine are also similar for the rough and flat systems. These results support the utility of a rigid/flat bilayer model for PB-based studies of PIP2 and PIP3 as long as the orientations are judiciously chosen.
引用
收藏
页码:155 / 163
页数:9
相关论文
共 5 条
  • [1] The Effects of Ca2+ on the Dynamics of PIP2 Containing Lipid Bilayers
    Mc Cabe, Ian
    Forstner, Martin B.
    BIOPHYSICAL JOURNAL, 2013, 104 (02) : 85A - 85A
  • [2] The Effects of Ca2+ on the Dynamics of PIP2 Containing Lipid Bilayers
    Mc Cabe, Ian P.
    Forstner, Martin B.
    BIOPHYSICAL JOURNAL, 2015, 108 (02) : 242A - 242A
  • [3] Effects of cholesterol and PIP2 on interactions between glycophorin A and Band 3 in lipid bilayers
    Qin, Xiaoxue
    Tieleman, D. Peter
    Liang, Qing
    BIOPHYSICAL JOURNAL, 2022, 121 (11) : 2069 - 2077
  • [4] Novel PIP2 lipid binding to focal adhesion kinase probed through coarse-grained molecular dynamics simulations
    Feng, Jun
    Mertz, Blake
    Bennett, C. Eugene
    ABSTRACTS OF PAPERS OF THE AMERICAN CHEMICAL SOCIETY, 2013, 246
  • [5] Interactions of the EGFR juxtamembrane domain with PIP2-containing lipid bilayers: Insights from multiscale molecular dynamics simulations
    Abd Halim, Khairul Bariyyah
    Koldso, Heidi
    Sansom, Mark S. P.
    BIOCHIMICA ET BIOPHYSICA ACTA-GENERAL SUBJECTS, 2015, 1850 (05): : 1017 - 1025