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
Overview
Affiliations
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 approximately 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 Angström of the plane formed by the phosphates of POPC; lipids in the surrounding shell are depressed by 0.5-0.8 A, but otherwise the phosphoinositides do not substantially perturb the bilayer. Finite size artifacts for ion distributions are apparent for systems of approximately 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 A for PIP2 and 19.4 +/- 0.3 A for PIP3) occurring near the most populated orientations from MD. These surfaces are well above the background height of 10 A 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.
Nanoscopy reveals integrin clustering reliant on kindlin-3 but not talin-1.
Wu Y, Cao Z, Liu W, Cahoon J, Wang K, Wang P Cell Commun Signal. 2025; 23(1):12.
PMID: 39773732 PMC: 11707915. DOI: 10.1186/s12964-024-02024-8.
Atomistic Simulations and Analysis of Peripheral Membrane Proteins with Model Lipid Bilayers.
Campbell O, Allsopp R, Klauda J, Monje V Methods Mol Biol. 2024; 2888:281-303.
PMID: 39699738 DOI: 10.1007/978-1-0716-4318-1_19.
Daear W, Mundle R, Sule K, Prenner E BBA Adv. 2023; 1:100021.
PMID: 37082006 PMC: 10074965. DOI: 10.1016/j.bbadva.2021.100021.
Modeling the membrane binding mechanism of a lipid transport protein Osh4 to single membranes.
Karmakar S, Klauda J Biophys J. 2022; 121(8):1560-1575.
PMID: 35247338 PMC: 9072576. DOI: 10.1016/j.bpj.2022.03.001.
Simulations of Kindlin-2 PIP binding domains reveal protonation-dependent membrane binding modes.
Palmere R, Case D, Nieuwkoop A Biophys J. 2021; 120(24):5504-5512.
PMID: 34813727 PMC: 8715242. DOI: 10.1016/j.bpj.2021.11.021.