Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2013 Jul 18;117(28):8322-9.
doi: 10.1021/jp401414y. Epub 2013 Jul 3.

Quantum and all-atom molecular dynamics simulations of protonation and divalent ion binding to phosphatidylinositol 4,5-bisphosphate (PIP2)

Affiliations

Quantum and all-atom molecular dynamics simulations of protonation and divalent ion binding to phosphatidylinositol 4,5-bisphosphate (PIP2)

David R Slochower et al. J Phys Chem B. .

Abstract

Molecular dynamics calculations have been used to determine the structure of phosphatidylinositol 4,5 bisphosphate (PIP2) at the quantum level and to quantify the propensity for PIP2 to bind two physiologically relevant divalent cations, Mg(2+) and Ca(2+). We performed a geometry optimization at the Hartree-Fock 6-31+G(d) level of theory in vacuum and with a polarized continuum dielectric to determine the conformation of the phospholipid headgroup in the presence of water and its partial charge distribution. The angle between the headgroup and the acyl chains is nearly perpendicular, suggesting that in the absence of other interactions the inositol ring would lie flat along the cytoplasmic surface of the plasma membrane. Next, we employed hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to investigate the protonation state of PIP2 and its interactions with magnesium or calcium. We test the hypothesis suggested by prior experiments that binding of magnesium to PIP2 is mediated by a water molecule that is absent when calcium binds. These results may explain the selective ability of calcium to induce the formation of PIP2 clusters and phase separation from other lipids.

PubMed Disclaimer

Figures

Figure 1
Figure 1
(A) The optimized geometry of PIP2 highlighting several intramolecular hydrogen bonds. (B) The atom-naming scheme used in the text.
Figure 2
Figure 2
(A) The results of a relaxed potential energy scan of the phosphodiester dihedral angle joining the headgroup of PIP2 to the glycerol moiety and the acyl chains. The dihedral angle is defined by the atoms labeled on the left of panel C. A value of 0 indicates that the plane defined by C1, O1, and P1 is aligned with the plane defined by O1, P1, and O11. (B) The angle between the ‘head’ and ‘tail’ of PIP2 measured after a geometry relaxation at the values of the phosphodiester dihedral angle reported in panel (A). The head-tail angle corresponds to the angle between the planes on the right of panel (C).
Figure 3
Figure 3
(A) The average surface area of PIP2 as calculated by squaring the maximum distance between oxygen atoms on two different phosphomonoester groups. The two histograms come from two separate QM/MM simulations with trajectories of 5 and 10 ps. (B) A comparison between the average surface area in simulations and Langmuir monolayer experiments of PIP2 without a divalent ion present, with calcium, or magnesium. The error bars are the standard deviation of the area. (C) The same data as in panel (B) but the values for calcium and magnesium have been normalized by the experimental or simulation reported average surface area in the absence of divalent ions.
Figure 4
Figure 4
(A) Free energy contours for calcium binding between the two phosphomonoester groups or to only the 4-phosphate group. (B) A scatter plot of the thermodynamic sampling path of the trajectory with renderings shown on the right for specific points along the path.
Figure 5
Figure 5
The free energy of spreading the phos phomonoester groups of PIP2 in the presence of calcium, magnesium, or potassium using classical US simulations. The insets show renderings from the simulations. Calcium (blue) binds tightly to the 4-phosphate group of PIP2. Magnesium (red) binds loosely to the 4-phosphate group and is coordinated by water molecules (TIP3 oxygens within 4 Å of the magnesium are shown as small blue spheres). Potassium does not appear to tightly bind PIP2 and is not shown. The shaded background represents the mean square difference between analyzing only the first half umbrella sampling trajectories and the full window.
Figure 6
Figure 6
The free energy of total deprotonation of the PIP2 phosphomonoester groups using QM/MM US simulations. The insets show the initial configuration, with proton H51 bound to the 5-phosphate, an intermediate structure with proton H51 bound to the 4-phosphate, and the final configuration with proton H51 (blue) forming a hydronium ion and water molecules hydrogen bonding (black dashed line) with the PIP2 phosphomonoester groups.
Figure 7
Figure 7
(A) The equilibrium distance to PIP2 obtained for the divalent ions, calcium or magnesium, in separate classical MD simulations and (B) the transient binding behavior of several potassium ions in a water sphere with a single PIP2 molecule. Renderings of snapshots from simulations are shown in the insets.
Figure 8
Figure 8
The free energy landscape for dissociating divalent ions from PIP2.
Figure 9
Figure 9
The free energy landscape for dissociating a proton from the 5-phosphate group of PIP2 in the presence of calcium or magnesium. The shaded background represents the mean square difference between analyzing only the first half umbrella sampling trajectories and the full window. In the insets, the proton to be dissociated is shown in black, magnesium in red, and calcium in blue.

Similar articles

Cited by

References

    1. Katso R, Okkenhaug K, Ahmadi K, White S, Timms J, Waterfield MD. Cellular Function of Phosphoinositide 3-Kinases: Implications for Development, Homeostasis, and Cancer. Annu. Rev. Cell Dev. Biol. 2001;17:615–675. - PubMed
    1. Yamaguchi H, Yoshida S, Muroi E, Kawamura M, Kouchi Z, Nakamura Y, Sakai R, Fukami K. Phosphatidylinositol 4,5-Bisphosphate and PIP5-Kinase Iα Are Required for Invadopodia Formation in Human Breast Cancer Cells. Cancer Sci. 2010;101:1632–1638. - PMC - PubMed
    1. Rodriguez S, Huynh-Do U. The Role of PTEN in Tumor Angiogenesis. J. Oncol. 2012;2012:141236. - PMC - PubMed
    1. Lassing I, Lindberg U. Polyphosphoinositide Synthesis in Platelets Stimulated with Low Concentrations of Thrombin Is Enhanced Before the Activation of Phospholipase C. FEBS Lett. 1990;262:231–233. - PubMed
    1. Levental I, Cēbers A, Janmey PA. Combined Electrostatics and Hydrogen Bonding Determine Intermolecular Interactions Between Polyphosphoinositides. J. Am. Chem. Soc. 2008;130:9025–9030. - PMC - PubMed

Publication types

LinkOut - more resources